GIMMS的NDVI數(shù)據(jù)常常用來進行長時間序列的植被情況分析，最新的NDVI3g數(shù)據(jù)集冰抢，時間跨度為1982-2015年，半月合成，空間分辨率大致為8km（0.08333度）忆畅，下載地址https://nex.nasa.gov/nex由于下載下來的是nc文件，部分用戶可能不熟悉該數(shù)據(jù)格式尸执，無法很好的使用該數(shù)據(jù)集家凯，本文提供一套基于matlab的處理方式缓醋。

% @author yinlichang3064@163.com
lon=ncread('F:\Global\NDVI3g\ndvi3g_geo_v1_2015_0712.nc4','lon');
lat=ncread('F:\Global\NDVI3g\ndvi3g_geo_v1_2015_0712.nc4','lat');
%tmn=ncread('C:\Users\lenovo\Documents\Tencent Files\1079192939\FileRecv\月最低氣溫數(shù)據(jù)集_2000_2002\tmn_2000_2002.nc','tmn');
lonlim_left=-180;
lonlim_right=180;
latlim_top=90;
latlim_bottom=-90;
GeoRef = georasterref('Rastersize',[2160,4320],'Latlim',[latlim_bottom,latlim_top],'Lonlim',[lonlim_left,lonlim_right]);
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('F:\Global\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=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,GeoRef);
    end
    m=m-12;
    i=i+1;
end

通過上述代碼可直接輸出每個月的TIF數(shù)據(jù)，而不用引入樣例數(shù)據(jù)了

更多需求绊诲，請查看個人介紹