登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
老师说我做的剖面是图中各线条在经度方向上的投影,并不代表图8中a,b,,c,d,e,f所对应的剖面。求问大神们如何修改,TRMM的HDF数据 Latitude = hdfread('G:\论文数据\数据库\天兔\2A25.20131007.90525.7.HDF', '/Swath/Latitude', 'Index', {[1 1],[1 1],[9246 49]});Longitude = hdfread('G:\论文数据\数据库\天兔\2A25.20131007.90525.7.HDF', '/Swath/Longitude', 'Index', {[1 1],[1 1],[9246 49]});rain = hdfread('G:\论文数据\数据库\天兔\2A25.20131007.90525.7.HDF', '/Swath/rain', 'Index', {[1 1 1],[1 1 1],[9246 49 80]});lat=double(Latitude);lon=double(Longitude);R=rain(5716,:,:);J3=Longitude(5716,:);A2=[20:-0.25:0.25];R=double(R);C=reshape(R,49,80);B=C';[Bsize1,Bsize2]=size(B);for i=1:Bsize1 for j=1:Bsize2 if B(i,j)<=0 dBB(i,j)=0; else dBB(i,j)=(B(i,j))/100; end endendpcolor(J3,A2,dBB);shading interp;colorbar;xlabel('longitude (degree)');%xlabel('latitude (degree)');ylabel('high-km (degree)');title('geo-RT-剖面')
|