Mapping toolbox绘制地理栅格数据

转自：http://blog.csdn.net/rumswell/article/details/16930823

Mapping toolbox绘制地理栅格数据:

[plain] view plain copy

1. %% worldmap1  
2. clc;clear;clf;  
3. load geoid; load coast  
4. figure(1); axesm robinson;framem on;  
5. set(gca,'Visible','off');  
6. %geoid(180*360)是地理高度(m),分辨率是1*1degree  
7. %geoidlegend也就是geoidrefvec=[1,90,0],是地理参考向量  
8. geoshow(geoid,geoidlegend,'DisplayType','texturemap')  
9. %设置一下色标  
10. cb=colorbar('southoutside');  
11. p=get(cb,'Position');%[x起点,y起点,宽度,高度]  
12. set(cb,'Position',[p(1)+0.2,p(2),p(3)-0.4,p(4)-0.01]);  
13. %绘制海岸线  
14. geoshow(lat,long,'color','k')  

使用matlab来绘制带地图背景的地理数据时很方便的.下面换一个投影:

[plain] view plain copy

1. %% worldmap1  
2. clc;clear;  
3. %海岸线数据  
4. coast = load('coast');  
5. figure('Color','w')  
6. axesm('mercator','MapLatLimit',[-90 90],'MapLonLimit',[-60 300])  
7. axis off; framem on; gridm on; mlabel on; plabel on;  
8. setm(gca,'MLabelLocation',60)  
9. setm(gca,'PLabelLocation',[-75:30:75])%设置纬度标签值  
10. geoshow(coast.lat,coast.long,'DisplayType','polygon')