- 积分
- 31801
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-7-5
- 最后登录
- 1970-1-1

|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 wlzhongouc 于 2013-11-30 15:18 编辑 以下无限制可见部分为管理员注
非常感谢猛牛的分享,也许你的积分还不够看到下面的具体内容,请理解作者文字和论坛大环境的来之不易,习惯于索取的人太多,如果你不是一个过客,你常来支持论坛,1000积分不是个很大的数字,如果你只是个过客,我们也只能对您说,我们的有些帖子不是为过客准备的,谢谢您的支持~论坛的发展需要大家的分享!
[ 批量读取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);
......(以下写你的程序啦 )
else
continue
end
end
end
end
@Aires
|
评分
-
查看全部评分
|