matlab存超过3g的文件,基于matlab的GIMMS3G-NDVI数据处理-从nc文件解析到月和年最大合成tif...

GIMMS的NDVI数据常常用来进行长时间序列的植被情况分析，最新的NDVI3g数据集，时间跨度为1982-2015年，半月合成，空间分辨率大致为8km(0.08333度)，下载地址https://nex.nasa.gov/nex由于下载下来的是nc文件，部分用户可能不熟悉该数据格式，无法很好的使用该数据集，本文提供一套基于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\月最低气温数据集_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数据，而不用引入样例数据了

更多需求，请查看个人介绍