- 积分
- 50
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
新手求教!处理GPM卫星降水数据,格式为HDF5,用m_map出图时遇到如下错误,找了半个月的错误,没找到,针对错误变量视窗也是对的,
烦请MATLAB大神解答!
[url=]下附部分输出部分代码[/url]
%output file name
file0='output\';
if exist(file0,'dir')==0
mkdir(file0); %create results directory
end
file1=FILE_NAME(1:8);
file2='.tiff';
file3='_';
file=[file0,file1,file3,year,mons,days,file3,DATAFIELD_NAME,file2];
m_proj('Lambert Conformal Conic','lon',lonlim,'lat',latlim) %select projection
m_pcolor(ChLon,ChLat,Chdata); %plot data
shading flat;
hold on
%China map
a=shaperead('bou2_4l.shp'); %read shp file
bou1_4lx=[a(:).X]; %longitude
bou1_4ly= [a(:).Y]; %latitude
m_plot(bou1_4lx,bou1_4ly,'-k') %plot map
for m=1:ns
m_text(slon(m),slat(m),'\bullet','color','red','fontsize',14);
m_text(slon(m)+0.3,slat(m),sname(m),'color','red','fontsize',8);
m_text(slon(m)+0.3,slat(m)-0.2,['(',num2str(round(myear_sdata(m),2)),')'],'color','red','fontsize',8);
end
%m_coast; %plot coast line
m_grid; %add lat/lon
hold off
%show color bar
colormap(jet(256));
h=colorbar();
set(get(h, 'title'), 'string', unit, 'FontSize', 12, 'FontWeight','bold');
title({[ave_date(1,:),'-',ave_date(nf,:),' '];DATAFIELD_NAME},'FontSize',16,'FontWeight','bold');
print('-dtiff','-r600',file);
clf;
出图也是出了,但是明显不对,烦请解答!
|
-
|