MATLAB coding/loop question

8 次查看（过去 30 天）

Benjamin Davidson 2019-4-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/453682-matlab-coding-loop-question

评论： Benjamin Davidson 2019-4-8

Hi there, I have a coding quesion! I'm writing a scrip for neuroimaging analysis. Some subjects have 1 functional scan, some have 2, it's about 50-50 and there's no way to predict. Is there any easy way to allow my code to be flexible for this? My code is below

for i = 1:2 
    for j = 1:2 
        filename = sprintf('/media/nir/SharedDrive/Ben/BATCHtester1/Sub014%d_Ses1/Sub014%d_Ses1_Scan_0%d_BOLD%d.nii', i, i, j + 1, j)
        BATCH.Setup.functionals{i}{j}= filename;
        clear filename;
    end
    end

Thank you!

8 个评论
显示 6更早的评论隐藏 6更早的评论

Guillaume 2019-4-1

在 MATLAB Online 中打开

Benjamin Davidson's comment mistakenly posted as an answer moved (and reformated) here:

Thanks for the reply! yes this code is just intended to build a path. Basically it's part of a larger code that goes through a large database of patient MRI scans, building a path to each patient's scans. In the code I have there, j = 1:2 indicates that many of the patients have 2 MRI's in their indirectory, but some of them actually have only 1.

Here's the entire code, and thanks in advance for any help!

clear BATCH;
BATCH.filename= '/media/nir/SharedDrive/Ben/BATCHtester1/conn_2.mat';
BATCH.Setup.isnew=1;
BATCH.Setup.nsubjects=2;
for i = 1:2 
    for j = 1:2 
        filename = sprintf('/media/nir/SharedDrive/Ben/BATCHtester1/Sub014%d_Ses1/Sub014%d_Ses1_Scan_0%d_BOLD%d.nii', i, i, j + 1, j)
        BATCH.Setup.functionals{i}{j}= filename;
        clear filename;
    end
end
for i = 1:2
    filename = sprintf('/media/nir/SharedDrive/Ben/BATCHtester1/Sub014%d_Ses1/Sub014%d_Ses1_Scan_01_ANAT1.nii', i, i)
    BATCH.Setup.structurals{i}= filename;
    clear filename;
end
BATCH.Setup.preprocessing.steps={'default_mni'};
BATCH.Setup.preprocessing.fwhm=8;
BATCH.Setup.preprocessing.voxelsize_func=2;
BATCH.Setup.preprocessing.sliceorder=[1:2:47,2:2:47];
BATCH.Setup.RT=3.0;
BATCH.Setup.analyses=[1,2];
BATCH.Setup.voxelmask=1;
BATCH.Setup.voxelresolution=1;
BATCH.Setup.outputfiles=[0,1,0];
BATCH.Setup.roi.names={'LeftVentralStriatum','RightVentralStriatum', 'LeftDorsalPutamen', 'RightDorsalPutamen', 'LeftMedialDorsalThalamus', 'RightMedialDorsalThalamus', 'LeftVentralPutamen', 'RightVentralPutamen', 'LeftDorsalStriatum', 'RightDorsalStriatum'};
conn_batch(BATCH);

Stephen23 2019-4-6

编辑：Stephen23 2019-4-6

在 MATLAB Online 中打开

Simplify the code by defining the root path before the loop:

N = 1400;
R = '/media/nir/SharedDrive/Ben/BATCHtester1';
for ii = 1:N
    D = sprintf('Sub014%d_Ses1', ii);
    S = dir(fullfile(R,D,'*.nii'));
    for jj = 1:numel(S)
        F = fullfile(R,D,S(jj).name);
        ... your code
    end
end