计算3D结构的分形维数

MATLABcalcFD：计算3D结构的分形维数(神经科学、图像处理、分形维数)

function [fdsubjects] = calcFD（subjectssubjectpathoptions）
% Calculate the fractal dimensionality of a 3D structure.
% Designed to work with intermediate files from FreeSurfer analysis pipeline
%   （ribbon.mgz aparc.a2009s+aseg.mgz and others）.
% Also can use other mgz volume as input （e.g. see ‘benchmark folder‘）.
%
% See ‘wrapper_sample.m‘ for an example of how to use the calcFD toolbox.
%
% REQUIRED INPUTS:
% subjects      = list of subjects names in a cell array
%                 alternatively accepts {‘.‘} to run on all subjects in folder
%
% subjectpath   = FreeSurfer ‘SUBJECTDIR‘ where standard directory structure is
%
% options       = specify details of running the analysis
%
% options.alg   = ‘dilate‘ | ‘boxcount‘
%
% options.countFilled = 0 | 1
%                       0   == Surface-only （FDs）
%                       1   == Filled volume （FDf）
%
% options.aparc = ‘Ribbon‘ | ‘Dest_aparc‘ | ‘Dest_select‘ | ‘DKT‘ | ‘Economo‘ | ‘none‘
%                 ‘Ribbon‘      == Cortical Ribbon （unparcellated）
%                 ‘Dest_aparc‘  == Parcellated cortical regions （Destrieux）
%                                   ** requires options.input.
%                 ‘Dest_select‘ == Any region in the aparc.a2009s+aseg.mgz volume
%                                   ** requires options.input.
%                 ‘DKT‘         == Parcellated cortical regions （DKT）.
%                                   ** requires aparc.DKTaltas+aseg.mgz （FS 6） or 
%                                      aparc.DKTaltas40+aseg.mgz （FS 5.3） to exist.
%                                   The volume can be generated （FS 5.3） using:
%                                   mri_aparc2aseg --s [SUBJECTID] --annot aparc.DKTatlas40
%                                   See Madan & Kensinger （2017 Brain Informatics） for further details.
%                 ‘Economo‘     == Parcellated cortical regions （von Economo-Koskinas）.
%                                   ** requires economo+aseg.mgz to exist.
%                                   The volume can be generated using:
%                                   mris_ca_label mris_anatomical_stats 
%                                   See Scholtens et al. （2018 NeuroImage） and 
%                                   Madan & Kensinger （2018 Eur J Neurosci） for further details.
%                 ‘none‘        == Binarized volume to be manually entered 
%                                   （e.g. benchmark volumes）.
%
% options.input = filename string required for ‘Dest_aparc‘ and ‘Dest_select
%               if options.aparc == ‘Dest_aparc‘
%                   This should be a file with the name ‘mask_*.txt‘
%                       where * is the value in options.input.
%                   File should have either 74 or 148 rows only 1 column.
%                       If only 74 values labels are assigned bilaterally.
%                   Value in each row is the label to assign to that parcellated region
%                       based on the Destri

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-06-02 20:16  cMadan-calcFD-7160b02\
     文件        1113  2018-06-02 20:16  cMadan-calcFD-7160b02\CITATION
     文件       35141  2018-06-02 20:16  cMadan-calcFD-7160b02\LICENSE
     文件        1374  2018-06-02 20:16  cMadan-calcFD-7160b02\PUBLICATIONS.md
     文件        6090  2018-06-02 20:16  cMadan-calcFD-7160b02\README.md
     文件          25  2018-06-02 20:16  cMadan-calcFD-7160b02\_config.yml
     目录           0  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\
     文件      137359  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_Sarmadillo.mgz
     文件      190870  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_Sbunny.mgz
     文件      226490  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_cube.mgz
     文件       93622  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_fibercup.mgz
     文件       76780  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_fibercup_disc.mgz
     文件      279022  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_menger01.mgz
     文件      317358  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_menger02.mgz
     文件      389660  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_menger04.mgz
     文件      139838  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_mug.mgz
     文件      209647  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_sphere.mgz
     文件      138054  2018-06-02 20:16  cMadan-calcFD-7160b02\benchmark\example_teapot.mgz
     目录           0  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\
     文件       12298  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\calcFD.m
     文件        1057  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\calcFD_boxcount.m
     文件         510  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\calcFD_dilate.m
     文件         695  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\calcFD_hollowVol.m
     文件         990  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\calcFD_save.m
     文件         768  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\calcFD_volCrop.m
     文件        7498  2018-06-02 20:16  cMadan-calcFD-7160b02\calcFD\load_mgh.m
     目录           0  2018-06-02 20:16  cMadan-calcFD-7160b02\examples\
     文件       14086  2018-06-02 20:16  cMadan-calcFD-7160b02\examples\calcFD_mask.xlsx
     文件      190347  2018-06-02 20:16  cMadan-calcFD-7160b02\examples\example_Sbunny.nii.gz
     文件         707  2018-06-02 20:16  cMadan-calcFD-7160b02\examples\example_nii.m
     文件          37  2018-06-02 20:16  cMadan-calcFD-7160b02\examples\lobes_legend.txt
............此处省略7个文件信息