- 积分
- 31301
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-7-5
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2016-5-6 11:03:18
|
显示全部楼层
看不到吗?设置修改了啊
内容如下:
批量读取nc格式的数据(针对不同日期的数据)示例:
for year=2003:2012
for month=1:12
mon=num2str(month+100);
strmon=mon(2:3);
% To Judge How many days in a month
if (month==1)||(month==3)||(month==5)||(month==7)||(month==8)||(month==10)||(month==12)
day=31;
elseif (month==4)||(month==6)||(month==9)||(month==11)
day=30;
elseif (mod(year,400)==0||mod(year,4)==0 && mod(year,100)~=0) && (month==2)
day=29;
else
day=28;
end
for id=1:day
mday=num2str(id+100);
strday=mday(2:3);
name=[F:\**,num2str(year),num2str(strmon),num2str(strday),'.nc']; %按时间变化的不同数据
out=netcdf.open(name,'nc_nowrite');
temp=netcdf.getvar(out,要导入的变量);
lat=netcdf.getvar(out,纬度);
lon=netcdf.getvar(out,经度);
......(以下自己设计吧,如作图什么的)
end
end
end
-----------------------------------------------------------------------------------------------------------------------------------------------------
批量读取hdf格式的数据示例:
for year=2003:2012
for month=1:12
for day=1:31
filename=[**,[num2str(i),num2str(month,'%.2d'),num2str(day,'%.2d')],'.hdf'];
if exist(filename,'file')
fileinfo=hdfinfo(filename);
data=hdfread(fileinfo.SDS);
......(以下写你的程序啦{:soso_e100:})
else
continue
end
end
end
end
|
|