- 积分
- 3400
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-4-21
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
程序如下。看了很多帖子 没有弄明白地图白化该怎么弄,因为例子大都是contour,照猫画虎之后发现也没报错但没有效果,图上还是没有白化。
黑体是m_pcolor和maskmap(我放进m_map工具箱里了,名字重命名了,m_maskmap,这里名字没问题)。求解~~
%%画图
%lat=hdfread(['E:\Modis\' 'MOD' num2str(2008) '(' num2str(i) ')' '.hdf'],'YDim');
%lon=hdfread(['E:\Modis\' 'MOD' num2str(2008) '(' num2str(i) ')' '.hdf'],'XDim');
lat2 = [15.5:1:55.5];
lon2 = [72.5:1:135.5];
lon1 = [72.5 135.5];
lat1 = [15.5 55.5];
guojie=shaperead('d:\calipsodata\southchinasea\bou2_4l.shp');%读取国界shp文件的内容
guojie1=shaperead('d:\calipsodata\southchinasea\bou1_4l.shp');%读取国界shp文件的内容
shengjie = shaperead('d:\calipsodata\southchinasea\bou2_4p.shp');
bou2_4lx=[guojie(:).X];%提取经度信息
bou2_4ly=[guojie(:).Y];%提取纬度信息
cn_provincex=[shengjie(:).X];
cn_provincey=[shengjie(:).Y];
m_proj('miller','lon',[72.5,136.5],'lat',[15.5,55.5]);
m_pcolor(lon2,lat2,aaod2);
shading interp
colormap jet;
colorbar
caxis([0,0.05])
hold on;
%m_maskmap('d:\calipsodata\southchinasea\bou2_4p.shp', true, 'lon', lon1, 'lat', lat1,'linewidth',0.5);
m_grid('linestyle','none','box','on','tickdir','out','fontsize',10,'fontname', 'Times New Roman','Xtick',[80 90 100 110 120 130], 'Ytick', [20 25 30 35 40 45 50 55] );
hold on;
m_plot(cn_provincex,cn_provincey,'k','LineWidth',0.5);
%m_plot(bou2_4lx,bou2_4ly,'k','LineWidth',0.5);
hold on
axes('Position',[0.114 0.175 0.09 0.2]);%创建坐标系时返回它的句柄
%set(gcf,'PaperPositionMode','auto');
m_proj('miller','lon',[106,122],'lat',[2,26]) %设置南海区域
m_pcolor(lon2,lat2,aaod2);
shading interp
m_plot(bou2_4lx,bou2_4ly,'k','LineWidth',0.5);
%m_plot(cn_provincex,cn_provincey,'k');%画省界
m_grid('linest','none','xticklabels',[],'yticklabels',[]);
|
|