先谢谢各位点进来我的求助帖子了,最近老师给的任务是处理CMORPH_V1.0三小时的数据,空间分辨率是0.25度,已经被校正过的。
我以前全然没有接触过matlab,在别人代码的帮助下通过matlab能够打开这个数据了,但是这个是三小时的数据,我不知道怎么能够把每个时段的数据都读出来,以及怎么能转换成栅格文件进行处理。。。还请各位老师、同学多多帮助。
下面附上我下载的数据,以及我用于打开文件的代码,这个代码只能保证把数据文件打开,但是不能够把三小时的数据栅格化提取出来,还请大家多多指导,谢谢大家了。。。更希望有老师或者同学能够提供完整的处理方法,实在不胜感激。。。
%%
clear;
close all;
clc;
%%
dir='I:\CMORPH_3HLY\1998';
file='CMORPH_19980101';
filename=[dir,'\',file];
fid=fopen(filename,'r','ieee-le');
data = fread( fid,'single' );
data = reshape( data, [1440,480,2,4 ] );
fclose( fid );
t1 = squeeze( data( :,:,1 ) );
t1( t1 == -999.0 ) = NaN;
wd = -59.875:0.25:-59.875 + ( 480-1 )*0.25; % lat
jd = 0.125:0.25:0.125 + ( 1440-1 )*0.25; % long
load coast
geoshow( lat,long,'color','k' );
hold on;
clf
axesm( 'miller','maplonlimit',[0 360],'maplatlimit',[-60 60],'frame','on',...
'grid','on','meridianlabel','on','parallellabel','on');
geoshow( lat,long,'color','k' );
pcolorm( wd,jd,t1');
colorbar;