画站点风场图

指甲钳

1. 
   

复制代码

正在学NCL，画了一个比较简单的站点数据的风场，用的是CSV格式，第一列是经度，第二列是维度，第三列是风向，第四列是风速，用";"做间隔的。希望对初学者有用吧，也希望高手们指点一下，大家可以多点分享脚本，共同进步。
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

begin

wks = gsn_open_wks("png","wind_1")     ; Open a workstation
;-----------------------------设置色标-------------------------

gsn_define_colormap(wks,"default") ;


;-----------------读取数据-----------------------------
fpath = "E:\sql_data\wind1.csv"

lines = asciiread(fpath,-1,"string")

delim = ";"
lon  = tofloat(str_get_field(lines,1,delim))
  lat  = tofloat(str_get_field(lines,2,delim))
  d = tofloat(str_get_field(lines,3,delim))
  s = tofloat(str_get_field(lines,4,delim))
;data@_FillValue = -9999




;-------------------------------------参数设置-----------------------------------------


mpres              = True

mpres@gsnMaximize   = True
mpres@gsnDraw       = False
mpres@gsnFrame      = False
  mpres@mpMinLonF    = 112.8
  mpres@mpMaxLonF    =  113.3
  mpres@mpMinLatF    =  22.75
  mpres@mpMaxLatF    =  23.35




  map = gsn_csm_map(wks,mpres)



;-------shp----------------
filename="gis/nhreverse/nhall.shp"

resp                  = True
resp@gsFillColor = (/"White"/)
polygon = gsn_add_shapefile_polygons(wks,map,filename,resp)


lnfilename="gis/nanhai/nanhai.shp"
lnres                  = True
lnres@gsEdgesOn = True
lnres@gsEdgeColor      = "white"
lnres@gsLineThicknessF = 1.0

polyline = gsn_add_shapefile_polylines(wks,map,lnfilename,lnres)

;-----------------------主要部分------------------------------------------------

draw(map)

uv=wind_component(s,d,0);将风场分解为U;v

u = uv(0,:)*1.94386 ;转单位
v = uv(1,:)*1.94386

wmsetp("col", 2)   
wmsetp("wbs", .03)
wmbarb(wks, lon, lat,u, v)

frame(wks)

end

贝頥嚫

指甲钳 发表于 2013-9-9 14:47

                               
登录/注册后可看大图

我的数据是内部数据啊。。网上当然找不到

哦哦 我说找找那种风向杆的数据资料 完全没有

bigbigice

用了wmbarb函数，还可以实用wm家族的其他函数，如wmvector画风场。

/猪头哼哼

-.-怎么只有第一行

指甲钳

/猪头哼哼 发表于 2013-9-6 11:21

                               
登录/注册后可看大图

-.-怎么只有第一行

不好意思。那个代码的不会弄，现在直接放在帖子里了

kongfeng0824

感谢分享！                              

woshixiancai

感谢分享

贝頥嚫

楼主 请问你的数据是在哪里下的什么样子的啊 找不到。。

指甲钳

贝頥嚫 发表于 2013-9-9 14:36

                               
登录/注册后可看大图

楼主 请问你的数据是在哪里下的什么样子的啊 找不到。。

我的数据是内部数据啊。。网上当然找不到

singapple

这样多多分享脚本确实很好，像我这样的菜鸟就能从中学习很多，多谢楼主~~