- 积分
- 15400
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-9-8
- 最后登录
- 1970-1-1
|
发表于 2017-3-6 11:51:30
|
显示全部楼层
MATLAB可以处理,这里发一个读月降水数据的代码:
fid=fopen('gpcp_v2.2_psg.2014','r','b','US-ASCII');
fseek(fid,576,'bof');
for i = 1:12
gpcp(:,:,i) = fread(fid,[144 72],'real*4');
end
lat = 88.75:-2.5:-88.75;
lon = 1.25:2.5:358.75;
[LON,LAT]=meshgrid(lon,lat);
m_proj('Equidistant cylindrical','lon',[0,360],'lat',[-90,90]);
m_coast('Color',[.65 .65 .65]);
hold on
m_contourf(LON,LAT,gpcp(:,:,4)','linestyle','none');
colorbar
m_grid('linestyle','none','tickdir','out'); |
|