1. 你滤波后会得到一个*.tts后缀的文件。第一步在TIMESAT中用TSF_fit2img将.tts转换为二进制文件。
2. 点击后需要输入一些参数，官方文档如下。
3. 你会得到一个或多个无后缀的二进制文件，最后用下面这段代码转制为.tif文件。
%% timesat二进制结果转换为tif
clear
input_name='I:\modis_data\hubei.tif'; %设置参考影像
[rows,cols]=size(input_name);
[~,R]=geotiffread(input_name); % 读取参考影像的投影信息
info=geotiffinfo(input_name);
sss = importdata('list.txt'); % 读取待转换文件清单
n = length(sss);
for i = 1:n
output_name = char(sss(i)); % 待转换的文件
fid = fopen(output_name,'r');
a1 = fread(fid,'*int16'); % 以int16读取数据
a2 = reshape(a1,cols,rows); % 转置
clear a1;
a2 = a2'; % 旋转
geotiffwrite(output_name,a2,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
end
2. 点击后需要输入一些参数，官方文档如下。
3. 你会得到一个或多个无后缀的二进制文件，最后用下面这段代码转制为.tif文件。
%% timesat二进制结果转换为tif
clear
input_name='I:\modis_data\hubei.tif'; %设置参考影像
[rows,cols]=size(input_name);
[~,R]=geotiffread(input_name); % 读取参考影像的投影信息
info=geotiffinfo(input_name);
sss = importdata('list.txt'); % 读取待转换文件清单
n = length(sss);
for i = 1:n
output_name = char(sss(i)); % 待转换的文件
fid = fopen(output_name,'r');
a1 = fread(fid,'*int16'); % 以int16读取数据
a2 = reshape(a1,cols,rows); % 转置
clear a1;
a2 = a2'; % 旋转
geotiffwrite(output_name,a2,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
end
