基于matlab的GIMMS3G-NDVI数据处理-从nc文件解

GIMMS的NDVI数据常常用来进行长时间序列的植被情况分析，最新的NDVI3g数据集，时间跨度为1982-2015年，半月合成，空间分辨率大致为8km（0.08333度），下载地址https://nex.nasa.gov/nex由于下载下来的是nc文件，部分用户可能不熟悉该数据格式，无法很好的使用该数据集，本文提供一套基于matlab的处理方式。

[~,R]=geotiffread('J:\全球性数据集\NDVI3g\example.tif');
info=geotiffinfo('J:\全球性数据集\NDVI3g\example.tif');
i=1;
k=[1:12];
m=1;
year=[1982  1982    1983    1983    1984    1984    1985    1985    1986    1986    1987    1987    1988    1988    1989    1989    1990    1990    1991    1991    1992    1992    1993    1993    1994    1994    1995    1995    1996    1996    1997    1997    1998    1998    1999    1999    2000    2000    2001    2001    2002    2002    2003    2003    2004    2004    2005    2005    2006    2006    2007    2007    2008    2008    2009    2009    2010    2010    2011    2011    2012    2012    2013    2013    2014    2014    2015    2015
];
format='monthNDVI3g.tif';
montha=[1:6];
monthb=[7:12];
e=dir(fullfile('J:\全球性数据集\NDVI3g','*nc4'));
while i<=length(e)
    a=ncread(e(i).name,'ndvi');
    while m<=length(k)
    j= m+1;
    h=j/2;
    b1=max(a(:,:,k(j)),a(:,:,k(m)));
    b=flipud(rot90(b1));
    m=m+2;
    if mod(i,2)==0
        g=monthb;
    else
        g=montha;
    end
    name1=strcat(int2str(year(i)),int2str(g(h)),format);
    geotiffwrite(name1,b,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
    end
    m=m-12;
    i=i+1;
end

本文中样例数据example.tif的制作方式可查看https://www.jianshu.com/p/1cae03e4452e不过需要对立面的一些参数进行修改，或给直接作者发邮件，看到后直接发送样例tif数据