麻辣GIS微信平台

更多 GIS 干货

微信关注不错过

IDL遥感应用入门(9):改进利用QuikScat数据绘制风速矢量图

在之前的文章IDL遥感应用入门(5):利用QuikScat数据绘制风速矢量图中,我绘制了QuikScat的风速矢量图,但是有一个缺陷就是没有实现没有实现图例和和说明信息,为此我特意请教了IDL的@杨铭伦大神,整理下了这篇文章。

QuikScat数据说明

QuikScat提供3个方面数据信息:表面风速,表面风向和降雨信息。

其中数据在0到250属于有效信息,大于250的值均属于无效信息。

计算风速的时候需要乘以0.2,风向需要乘以1.5

更多详细信息可以查看这里

代码实现

;Author:杨铭伦
;2014-01-13

pro yangminglun
filename=file_search('D:\3','*.gz')
device,decomposed=1
!p.color=0
!p.background='FFFFFF'xl
;读取某一目录下的所有数据
for data=0,0 do begin
  a=bytarr(1440,720,3)
  openr,lun,filename[data],/compress,/get_lun
  readu,lun,a
  free_lun,lun
  ws=congrid(a(*,*,0),180,90)*0.2
  dir=congrid(a(*,*,1),180,90)*1.5/180*!pi

  ;构建vector函数所需参数
  u=ws*sin(dir)
  v=ws*cos(dir)
  x=findgen(1440)*0.25
  y=findgen(720)*0.25-90
  x=congrid(x,180)
  y=congrid(y,90)

  ;构建图例
  u(50,70)=30
  v(50,70)=0

  ;出图显示
  window,data,xsize=1440,ysize=720
  name=strarr(7)
  name[data]=string(20060700+data+1)
  map_set,0,180,0,title=name(data),/continents,/grid,/label
  velovect,u,v,x,y,missing=250*0.2,/overplot

  xyouts,103,49,'10m/s',color='FF0000'xl
 endfor
end  

效果图

一些说明

文中file_search('D:\3','*.gz') 是定义的数据路径,这里共享下我用的实验数据 。猛击这里下载

声明:本文所有的代码来自@杨铭伦大神,非本人所写,在此表示强烈的感谢之情。

麻辣GIS-Sailor
作者:
GIS爱好者,学GIS,更爱玩GIS。
微博关注
手机阅读
赞赏支持
手机阅读
微信捐助麻辣GIS
微信打赏
支付宝捐助麻辣GIS
支付宝打赏

声明

1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。

2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。

3.如果本文对您有所帮助,请您对我们团队进行 打赏捐助让我们在传播3S的路上可以走得更远,不胜感激。PS:如果你是学生党,请优先把经费用于购买学习资料 以及 与小哥哥/小姐姐约会上:-)

如您有疑问,可在文末留言,或到麻辣GIS QQ群中提问。

QQ群(一):291616564 加入QQ群

QQ群(二):166408035 加入QQ群

QQ群(三):627853279 加入QQ群

QQ群(四):436386604 加入QQ群

QQ群(五):606176554 加入QQ群

QQ群(六):946178380 加入QQ群

QQ群(七):861271808 加入QQ群

新疆互助群:910717627 加入QQ群

微信公众号:malagis,扫描右边二维码直接关注。

发表评论