对matlab里自带的类似m_map函数的疑问，内附该函数的调用方式，Matlab7.0以上测试可用

容与

卫星实习课虽然学校很坑爹（向来如此），老师很无趣，课程很无聊，但学到了我觉得是新东西的几个函数（这门课对我唯一的用途，妈蛋，1学分80块钱呢），可以不用调用m_map文件夹的函数来绘图，具体是：
axesm( 'giso','Grid','on','MapLonLimit',[-120 -50],'MapLatLimit',[15 60]);
plotm(lat,long,'k'); % plot sea-land boundry
setm(gca,'parallellabel','on','meridianlabel','on','plabellocation',5,'mlabellocation',5,'labelformat','signed'); setm(gca,'mlabelparallel',-90);
surfm(plat,plon,bt);
效果如下：

但不知道这一类的是否有其他的函数（已测试有contourm），求大神指教 satellite-practise.rar (533.12 KB, 下载次数: 11)

2015-7-10 11:29 上传

点击文件名下载附件
程序与文件

1. <p><span style="line-height: 1.5;">% read brightness temperature value [K]</span></p><p>fid=fopen('bt.txt','r');</p><p>bt=fscanf(fid,'%f',[98 351]);  % [98 351] 数据的大小</p><p>fclose(fid);</p><p>
   
2. </p><p>% read latitude value </p><p>fid=fopen('latitude.txt','r');</p><p>plat=fscanf(fid,'%f',[98 351]);</p><p>fclose(fid);</p><p>
   
3. </p><p>% read longitude value </p><p>fid=fopen('longitude.txt','r');</p><p>plon=fscanf(fid,'%f',[98 351]);</p><p>fclose(fid);</p><p>
   
4. </p><p>%</p><p>
   
5. </p><p>load coast;   %装入全球的矢量海岸线 COAST.MAT (9589*9589)经纬度值</p><p>
   
6. </p><p>%%亮度的选择性输出</p><p> %bt(find(bt>240 ))=NaN ;</p><p>%bt(find(bt<280))=NaN;</p><p>%bt(find(bt>240 & bt<280))=NaN ;</p><p> </p><p>%----------plot using auto colorbar--------------</p><p>
   
7. </p><p>figure;axesm( 'giso','Grid','on','MapLonLimit',[-120 -50],'MapLatLimit',[15 60]);  %地图的范围</p><p>plotm(lat,long,'k'); % plot sea-land boundry   %plotm 在前面map范围内的画图</p><p> setm(gca,'parallellabel','on','meridianlabel','on','plabellocation',5,'mlabellocation',5,'labelformat','signed');</p><p> setm(gca,'mlabelparallel',-90);</p><p> surfm(plat,plon,bt);</p><p> set(gca,'CLim',[200 300]); %色标的范围 </p><p> colorbar;tightmap;</p><p>%---------plot using defined colorbar---</p><p>
   
8. </p><p>levels=200:12:300;  % every 10 K one color; you can modify this interval</p><p>[colormap_levels]=colormap_user_define(levels);</p><p>
   
9. </p><p>figure;axesm( 'giso','Grid','on','MapLonLimit',[-120 -50],'MapLatLimit',[15 60]);</p><p>plotm(lat,long,'k'); % plot sea-land boundry</p><p> setm(gca,'parallellabel','on','meridianlabel','on','plabellocation',5,'mlabellocation',5,'labelformat','signed');</p><p> setm(gca,'mlabelparallel',-90);</p><p> surfm(plat,plon,bt);</p><p> set(gca,'CLim',[200 300]);</p><p>colormap(colormap_levels);</p><p>caxis([levels(1) levels(end)]);</p><p>ch=colorbar;</p><p>set(ch,'ytick',levels);</p><p>set(ch,'yticklabel',levels); </p><p>tightmap; </p><p>
   
10. </p><p>%-----------using RGB color---------</p><p>color_levels =[</p><p>     0         0    0.6  </p><p>      %1    1            1</p><p>     % 0.5    0.5         0.5</p><p>      0         0    1.0000       </p><p>      0    1.0000    1.0000 </p><p>      0    0.6        0</p><p>      0    1.0000         0</p><p>      1.0000    1.0000         0</p><p>      1.0000    0.6000         0 </p><p>      0.8000    0.4000         0</p><p>      1.0000         0         0</p><p>      0.6000         0         0</p><p>      0.6000         0    0.6000</p><p>      1.0000         0    1.0000];</p><p>  </p><p>  levels=200:10:300; %  every 10 K one color; you can modify this interval</p><p> % levels=200:20:300;</p><p>  figure;axesm( 'giso','Grid','on','MapLonLimit',[-120 -50],'MapLatLimit',[15 60]);</p><p>  plotm(lat,long,'k'); % plot sea-land boundry</p><p> setm(gca,'parallellabel','on','meridianlabel','on','plabellocation',5,'mlabellocation',5,'labelformat','signed');</p><p> setm(gca,'mlabelparallel',-90);</p><p> surfm(plat,plon,bt);</p><p> set(gca,'CLim',[200 300]);</p><p> pleng=length(levels);</p><p> colormap(color_levels(1:pleng-1,:));</p><p> caxis([levels(1) levels(end)]);</p><p> ch=colorbar; %set(ch,'Location','SouthOutside'); </p><p> %set(ch,'Location','SouthOutside');</p><p> set(ch,'YTick',levels,'FontSize',14,'fontname','Times New Roman');</p><p> set(ch,'YTickLabel',levels,'FontSize',14,'fontname','Times New Roman');</p><p> tightmap;</p><p> </p><p> break</p>

复制代码

沙颖凯

搜 MATLAB mapping toolbox

但m_map要好用的多而且免费。

smart_abell

这个相当好呢，有数据有例子，更容易看懂。感谢分享

容与

沙颖凯 发表于 2015-7-11 06:48
搜 MATLAB mapping toolbox

但m_map要好用的多而且免费。

找到了，谢啦！

aas314

值得学习 赞大赞

aas314

好东西，正在学习中