ndvi matlab,基于Matlab的NDVI最大合成

采用最大合成法将月数据ndvi合成年数据NDVI,本文以spot的月ndvi为例进行,利用matlab平台,代码如下所示:

%author yinlichang3064@163.com

[a,R]=geotiffread('F:\项目\data\spot_月ndvi\ndvi200001.tif');%先投影信息

info=geotiffinfo('F:\项目\data\spot_月ndvi\ndvi200001.tif');

for i=2000:2015

datasum=zeros(size(a,1)*size(a,2),12); %size(a,1)和size(a,2)分别是进行合成的图像的行列号

for j=1:12

if j<10

data=importdata(strcat('F:\项目\data\spot_月ndvi\ndvi',int2str(i),'0',int2str(j),'.tif'));

else

data=importdata(strcat('F:\项目\data\spot_月ndvi\ndvi',int2str(i),int2str(j),'.tif'));

end

data=reshape(data,size(a,1)*size(a,2),1);

datasum(:,j)=data;

end

maxndvi=zeros(size(a,1),size(a,2))+NaN;

for n=1:length(datasum)

data=datasum(n,:);

datamax=max(data);

if datamax>-1

maxndvi(n)=datamax;

end

end

filenameet=strcat('F:\项目\data\2000-2015年NDVI\spot_ndvi_',int2str(i),'.tif');

geotiffwrite(filenameet,maxndvi,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);

end

通过上述方式就能够批量输出tif格式的2000-2015年的NDVI最大合成数据了

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值