请把这个输入格式改成合适的clc clear [aa,R]=geotiffread('G:\ERA5hr\fc\modis\20010306.tif');%先投影信息 info=geotiffinfo('G:\ERA5hr\fc\modis\20010306.tif'); [m,n]=size(aa); begin_year=2011;%用户修改 end_year=2021;%用户修改 long=end_year-begin_year+1; ndvisum=zeros(m*n,long); k=1; for year=begin_year:end_year ndvi=importdata(['G:\ERA5hr\fc\modis\',int2str(year),'.tif']) ; %根据自己名称修改，本数据名称为'anveage_ndvi.tif' %注意数据的有效范围 ndvi(ndvi<-10000)=NaN; %有效范围是-1到1 ndvisum(:,k)=reshape(ndvi,m*n,1); k=k+1; end

好的，这是一个MATLAB代码段，但输入格式需要稍作修改，以下是修改后的代码段：

clc
clear

% 读取地理信息TIFF文件
[aa, R] = geotiffread('G:\ERA5hr\fc\modis\20010306.tif');
info = geotiffinfo('G:\ERA5hr\fc\modis\20010306.tif');

[m, n] = size(aa);
begin_year = 2011; % 用户修改
end_year = 2021; % 用户修改
long = end_year - begin_year + 1;
ndvisum = zeros(m * n, long);
k = 1;

for year = begin_year:end_year
    % 导入NDVI数据
    ndvi = importdata(['G:\ERA5hr\fc\modis\', int2str(year), '.tif']);
    % 注意数据的有效范围
    ndvi(ndvi < -10000) = NaN; % 有效范围是-1到1
    ndvisum(:, k) = reshape(ndvi, m * n, 1);
    k = k + 1;
end

希望这个修改能满足你的需求。