IDL编程学习之风场图

15195775117

风场图是气象上典型的向量图，IDL里使用VECTOR函数
图例：

其中，向量大小可以用箭头长度和颜色来表示。

风之精灵

支持原创，鼓励楼主继续给出好的帖子

15195775117

fig.vector_style=1设置风向为气象上的F形标志：

15195775117

iVector还可以画流场：
IDL> u=randomu(1,20,20)-0.5
IDL> v=randomu(2,20,20)-0.5
IDL> ivector,u,v,/streamlines,x_streamparticles=10,y_streamparticles=10,head_size=0.1,streamline_nsteps=200

wuyula

{:eb502:}{:eb502:}{:eb502:}

wuyula

{:5_213:}{:5_213:}{:5_213:}{:5_213:}

15195775117

                               
登录/注册后可看大图

IDL绘制F形风场图，是按m/s的单位绘制的，由于平时风并不大，所以F形风标常常是一道横杆
-------------------
F形风场案例：
数据见附件sav文件，读取之：
restore,'C:\Users\Administrator\Desktop\book.sav'
纬向风：
u=transpose(book[0,*])
经向风：
v=transpose(book[1,*])
风速标量：
spd=transpose(book[2,*])
风速方位角：
an=transpose(book[3,*])
经度：
lon=transpose(book[4,*])
纬度：
lat=transpose(book[5,*])
存放图的路径：
figfile='C:\Users\Administrator\Desktop\000.png'
画个风速分布图，作为风向图的衬底，这样风向图看起来明显一点：
fig=contour(spd,lon,lat,rgb_table=33,/fill,/buffer)
;风标设置为黑色，vector_style=1表示F形，vector_style=0表示箭头形
fig = VECTOR(u,v,lon,lat,head_size=1,axis_style=0,position=[0,0,1,1],$
  LENGTH_SCALE=0.5,sym_color='k',vector_style=1,/overplot,/buffer)
  ;AUTO_COLOR=1,RGB_TABLE=33
存在桌面：
  fig.save,figfile
  fig.close
将图片半透明化，这是我自己写的函数，其实可以略去：
；  transpng,figfile
  制作KMZ：
  arctic = IMAGE(figfile,MAP_PROJECTION='Stereographic',GRID_UNITS=2,$
    LIMIT=[min(lat),min(lon),max(lat),max(lon)],$
    IMAGE_LOCATION=[min(lon),min(lat)],/buffer,$
    IMAGE_DIMENSIONS=[max(lon)-min(lon),max(lat)-min(lat)],$
    /CURRENT,NAME='Arctic Research')
  arctic.SAVE,'C:\Users\Administrator\Desktop\contour_kml.kmz'
  arctic.close
删除桌面图片
  file_delete,figfile
  ;打开kmz：
  spawn,'C:\Users\Administrator\Desktop\contour_kml.kmz'

qxjy009

15195775117 发表于 2019-2-20 14:28
IDL绘制F形风场图，是按m/s的单位绘制的，由于平时风并不大，所以F形风标常常是一道横杆
-------------- ...

请教一下，您用的是idl哪个版本，或者是导入了什么函数库吗，我用的7.1不能识别fig和vector

qxjy009

我用您前面提供的脚本做的，编译就报错了
C:\Users\Administrator\Desktop\TIM图片20190808095016.png

15195775117

qxjy009 发表于 2019-8-8 09:39
请教一下，您用的是idl哪个版本，或者是导入了什么函数库吗，我用的7.1不能识别fig和vector

我用的是8.5，7.1我没用过，会不会7.1没有快速可视化？