matlab实现tca（temporal cluster analysis）算法

function addict_otca

clear
fm_path =‘D:\testdata\task240\s\‘    % 设定图像数据目录          

time_num = 240;              % 设定采样时间点数，即图像数目 11-260
%points_num = 53*63*46;%64*64*64;%16;       % 设定一时刻采样像素点数
points_num = 64*64;
slice_num = 16;

%---------------- 计算每个像素点的-----------------
baseline = zeros（points_num slice_num）;
for i=1:time_num          % 每次读取一副图像
    %---------------- 构造图像文件名-----------------
    zerofill = ‘00‘;
    if （i>=1） && （i<80）
        zerofill = ‘0‘;
    elseif （i>=80）
        zerofill = ‘‘;        
    end
    
    %filename1 = [fm_path‘20040424_K21_S012_B.MR.HEAD_FMRI.3.1.2004.4.24.13.28.32.703000.63088175x‘int2str（i）‘.img‘]
    filename1 = [fm_path‘sns‘int2str（205+i）‘a‘zerofillint2str（i+20）‘.img‘];
    %--------------------------------------------------------
    
    %---------------- 读入一副图像  -----------------
    fid1=fopen（filename1‘rb‘）;
    image_buffer1=fread（fid1[points_numslice_num]‘uint16‘）;    
    
    %image_buf