MATLAB地图上绘制矩形框

haobang008

由于研究的需要，需要用矩形框给出研究的范围，写了一个函数，和大家分享以下：
效果图如下，后面是调用的代码和函数；

矩形框示意图

附件是basemap函数以及使用到的数据mat文件；

basemap.m (2.9 KB, 下载次数: 58)

2016-5-26 14:12 上传

点击文件名下载附件
源代码

basemap.mat (367.01 KB, 下载次数: 50)

2016-5-26 14:12 上传

点击文件名下载附件
数据文件

1. 
   
2. 
   
3. % 利用basemao 绘制地图模板
   
4. % 单张图片m_map效果会更好,如果利用subplot函数,最好使用basemap 进行绘制进行考虑;
   
5. 
   
6. % contour 画等值线图;
   
7. % hold on;
   
8. % 需要注意的是basemap的语法,分别是起点和数据的范围;
   
9. 
   
10. % 该范围并没有实际用处,只是示意而已;
    
11. lon_range = [110 290];
    
12. lat_range = [-54 65];
    
13. % 1代表糊住陆地 0 代表画岸线
    
14. basemap(-54,119,110,180,1,[0.7 0.7 0.7]);
    
15. hold on;
    
16. basemap(-54,119,110,180,0,[0.0 0.0 0.0]);
    
17. 
    
18. set(gca,'XTick',120:30:270)
    
19. set(gca,'XTickLabel',{'120°E','150°E','180°','150°W','120°W','90°W'})
    
20. set(gca,'YTick',-50:20:50)
    
21. set(gca,'YTickLabel',{'60°S','45°S','30°S','15°S','0°','15°N','30°N','45°N','60°N'});
    
22. 
    
23. % 给出纬向风指数位置;
    
24. latm_c=[-5 5];
    
25. lonm_c=[155 180];
    
26. line_color_c='b';
    
27. latm_w=[-5 5];
    
28. lonm_w=[120 155];
    
29. line_color_w='r';
    
30. latm_e=[-5 5];
    
31. lonm_e=[180 260];
    
32. line_color_e='c';
    
33. 
    
34. line_width=1.5;
    
35. 
    
36. research_range(latm_w,lonm_w,line_color_w,line_width);
    
37. research_range(latm_e,lonm_e,line_color_e,line_width);
    
38. research_range(latm_c,lonm_c,line_color_c,line_width);
    
39. 
    

复制代码

下面是调用的函数

1. function ok=research_range(latm,lonm,line_color,line_width)
   
2. 
   
3. hold on;
   
4. %
   
5. x(1)=lonm(1);x(2)=lonm(1);
   
6. y(1)=latm(1);y(2)=latm(2);
   
7. plot(x,y,line_color,'LineWidth',line_width);
   
8. 
   
9. x(1)=lonm(2);x(2)=lonm(2);
   
10. y(1)=latm(1);y(2)=latm(2);
    
11. plot(x,y,line_color,'LineWidth',line_width);
    
12. 
    
13. x(1)=lonm(1);x(2)=lonm(2);
    
14. y(1)=latm(1);y(2)=latm(1);
    
15. plot(x,y,line_color,'LineWidth',line_width);
    
16. 
    
17. x(1)=lonm(1);x(2)=lonm(2);
    
18. y(1)=latm(2);y(2)=latm(2);
    
19. plot(x,y,line_color,'LineWidth',line_width);
    
20. 
    
21. end

复制代码

浅笑无痕

{:5_213:}感谢！！

haobang008

浅笑无痕 发表于 2016-5-26 16:06
感谢！！

有个小小的瑕疵，希望注意哦

wyx123

学习中，谢谢分享

sid_nancy

刚刚发帖求助如何画矩形框，就看到了你的帖子~~很有用~~谢谢

jsdfzb

学习中，值得借鉴，好材料

aas314

学习啦 谢谢楼主

haobang008

sid_nancy 发表于 2016-5-27 09:58
刚刚发帖求助如何画矩形框，就看到了你的帖子~~很有用~~谢谢

非常基础的内容哈，相互交流进步，有什么好想法可以挂上来

sid_nancy

haobang008 发表于 2016-5-29 15:26
非常基础的内容哈，相互交流进步，有什么好想法可以挂上来

我在已经画好的图上想给研究区域加矩形框，所以直接调用了你的research_range函数，但是出来的图里面并没有加上矩形框。。我是刚刚开始学习MATLAB，不太懂到底是哪里出错，能否帮我看看我的脚本到底是哪里出了问题~~~感激不尽
    figure
    m_proj('Miller','lon',[90 140],'lat',[15 60])
    hold on;
    m_pcolor(X,Y,vec3(:,:,1)')
    shading interp
    colormap(othercolor('BuDRd_18'));
    hold on;
    m_coast('line','color','k');
     m_grid('linest','none','tickdir','out','linewi',1,...
            'xtick',[90:5:140],'ytick',[15:5:60]);
    shading interp;
    colorbar('southoutside');
    k=caxis;
    kk=max(abs(k));
    caxis([-kk kk]);
    hold on;
    latm=[35 43];
    lonm=[110 120];
    line_color='r';
    line_width=1.5;
    research_range(latm,lonm,line_color,line_width);
print -djpeg -r300 test.jpeg

haobang008

sid_nancy 发表于 2016-5-29 16:52
我在已经画好的图上想给研究区域加矩形框，所以直接调用了你的research_range函数，但是出来的图里面并没 ...

并没有调用m_map的指令，而是采用了basemap的方法
所以你在使用m_map绘图的时候，是不能用basemap的方法加上矩形框的
但是解决办法也是有的，你把函数research_range里面的plot全部替换成m_plot试试即可。
下面是效果图：