IDL采用axis给IMAGE绘图绘制经纬度
发布时间: 2013-11-26
所属分类: IDL修炼秘籍
在之前的文章里,我提到可以用IMAGE函数自带的MAP_PROJECTION属性来使IDL的图像显示经纬度。如图:
但是之后老师告诉我说,这是错误的,倒不是说MAP_PROJECTION的使用是错误的,而是这幅图像本来就是投影好的。所以考虑另一种方法。代码如下:
spdimg = image(wdspd,rgb_table=2,title=diplayname) grdx= axis('x',location=[0,0],axis_range=[0,360],coord_transform=[0,0.25],minor=0,major=19) grdy= axis('y',location=[0,0],axis_range=[-90,90],coord_transform=[-90,0.25],minor=0,major=7)
主要的实现是下面的两行代码,效果如下:
这样的效果就是正确的了。
相关阅读
声明
1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。
2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。
手机阅读
公众号关注
知识星球
手机阅读
最新GIS干货
私享圈子
有个情况咨询一下,这种方法,只能在坐标轴显示数字,如果是小范围的一个区域,坐标得用到分和秒的话,要怎么做?
你可以参考这篇文章:http://malagis.com/idl-on-axis-function.html,通过设置axis的major设置一级刻度线的数目,通过minor设置二级刻度
十分感谢!
不必客气~