上一篇推送介绍了如何将下载的nc格式的NDVI数据转换成tif格式，这样结果得到的是旬NDVI【GIMMS NDVI数据处理-nc格式转tif】。今天分享旬NDVI的月最大值合成和月NDVI的年最大值合成。

1旬NDVI最大值合成月NDVI

[a,R]=geotiffread('E:\data\GIMMS NDVI3g V1.0 1982-2015\tif_data\ndvi\1982_1_ndvi.tif');info=geotiffinfo('E:\data\GIMMS NDVI3g V1.0 1982-2015\tif_data\ndvi\1982_1_ndvi.tif');for year=1982:1982  for month=1:12        datasum=zeros(size(a,1)*size(a,2),2);         data1=importdata(strcat('E:\data\GIMMS NDVI3g V1.0 1982-2015\tif_data\ndvi\',int2str(year),'_',int2str(month),'.5_ndvi.tif'));        data2=importdata(strcat('E:\data\GIMMS NDVI3g V1.0 1982-2015\tif_data\ndvi\',int2str(year),'_',int2str(month),'_ndvi.tif'));        data3=reshape(data1,size(a,1)*size(a,2),1);        data4=reshape(data2,size(a,1)*size(a,2),1);        datasum(:,1)=data3;datasum(:,2)=data4;        maxndvi=zeros(1,size(a,1)*size(a,2))+NaN;        for n=1:length(datasum)            data=datasum(n,:);            datamax=max(data);            if datamax>-1                maxndvi(1,n)=datamax;            end        end        maxndvi=reshape(maxndvi,size(a,1),size(a,2));        filename=strcat('E:\data\GIMMS NDVI3g V1.0 1982-2015\MVC_month\',int2str(year),'_',int2str(month),'mvc.tif');        geotiffwrite(filename,maxndvi,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);    endend

关键是将每个像元当月上下旬的NDVI值读取成数组，然后用max函数获取两者最大值，作为月最大NDVI。2月NDVI最大值合成年NDVI

[a,R]=geotiffread('E:\data\GIMMS NDVI3g V1.0 1982-2015\MVC\1982_1ndvi.tif');info=geotiffinfo('E:\data\GIMMS NDVI3g V1.0 1982-2015\MVC\1982_1ndvi.tif');for year=1982:2015    datasum=zeros(size(a,1)*size(a,2),12);         for month=1:12            data=importdata(strcat('E:\data\GIMMS NDVI3g V1.0 1982-2015\MVC\',int2str(year),'_',int2str(month),'ndvi.tif'));            data=reshape(data,size(a,1)*size(a,2),1);            datasum(:,month)=data;        end        maxndvi=zeros(1,size(a,1)*size(a,2))+NaN;        for n=1:length(datasum)            data=datasum(n,:);            datamax=max(data);            if datamax>-1                maxndvi(1,n)=datamax;            end        end        maxndvi=reshape(maxndvi,size(a,1),size(a,2));        filename=strcat('E:\data\GIMMS NDVI3g V1.0 1982-2015\MVC年\',int2str(year),'_mvc.tif');        geotiffwrite(filename,maxndvi,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);end

关键是将每个像元12个月的NDVI值读取成数组，然后用max函数获取其中最大值，作为年最大NDVI。

说说线性规划

VOSviewer文献综述

泰勒图的MATLAB实现

深度学习--卷积神经网络基础

分享一套数据---GIMMS NDVI

ArcPy批量定义投影和批量投影转换

机器人局部规划算法--DWA算法原理

ArcGIS时间滑块实现车辆轨迹动态展示

GPS数据处理---在野外采样寻点中的应用

Python散点图矩阵--以空气质量数据为例

▼更多精彩推荐，敬请关注我们▼