TIMESAT 无格式文件迭代转tif

%% timesat二进制结果转换为tif
clear
input_name='2014_01_01.tif'; %设置参考影像
[~,R]=geotiffread(input_name); % 读取参考影像的投影信息
rows = R.RasterSize(1);
cols = R.RasterSize(2);
info=geotiffinfo(input_name);

maindir = '.\DL_result';
subdir  = dir( maindir );

for i = 1 : length( subdir )
    if( isequal( subdir( i ).name, '.' )||...
        isequal( subdir( i ).name, '..'))               % 如果不是目录则跳过
        continue;
    end
    dirpath = fullfile( maindir, subdir( i ).name);
    fid = fopen(dirpath,'r');
    a1 = fread(fid,'float'); % 以int16读取数据
    a2 = reshape(a1,cols,rows); % 转置
    clear a1;
    a2 = a2'; % 旋转
    geotiffwrite( subdir( i ).name,a2,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
end