wrf模式高度场等值线的平滑问题

半颗麦芽糖

王老师您好，今天我使用脚本画了一下我的模式高度场的图，但是等值线挺乱的，所以考虑使用平滑，但是在论坛上没有找到相关的语句，只找到了wContour，所以想问问王老师改用什么语句，附上脚本和图，请王老师指导，谢谢！
f = addfile('E:/gongzuo/3.8.1/test38101.ctl')
t = datetime.datetime(2014,10,26,12)
ph = f['PH'][[t],[500],:,:]
phb = f['PHB'][[t],[500],:,:]
height = (ph+phb)/98.1
#Plot
axesm()
lchina = shaperead('F:/11/MeteoInfo/map/bou2_4p.shp')
ltibet = shaperead('F:/11/MeteoInfo/map/DBATP/DBATP_Polygon.shp')
geoshow(lchina, edgecolor='k')
geoshow(ltibet, edgecolor='r',size=2)
levs = arange(544,588,2)
layer = contourm(height,levs, proj=f.proj)
colorbar(layer)
xlim(70,110)
ylim(25,45)

超级大坏蛋

{:eb502:}{:eb502:}

半颗麦芽糖

超级大坏蛋 发表于 2017-1-3 08:31

伟哥给点指导意见

超级大坏蛋

半颗麦芽糖 发表于 2017-1-3 09:15
伟哥给点指导意见

我这种菜鸟对你这高深程序看不懂，何来指导之说！

MeteoInfo

有空了我加上平滑的函数吧，不过平滑会改变原始值。目前你可以对数组进行稀疏化处理（比如隔4个点取一个值），然后再绘制等值线。

例如：
height = height[::4,::4]

半颗麦芽糖

MeteoInfo 发表于 2017-1-3 10:34
有空了我加上平滑的函数吧，不过平滑会改变原始值。目前你可以对数组进行稀疏化处理（比如隔4个点取一个值 ...

按照王老师所说的加上了height=height[::4,::4]进行稀疏化处理，效果好了很多，期待加上平滑函数的效果，附上修改后的图。

听海落雪

格点太密，就是平滑处理了效果也不太好。就像用 2.5×2.5 的再分析资料画出的高度场 比 1.0×1.0 再分析资料画出的高度场平滑。

半颗麦芽糖

听海落雪 发表于 2017-1-5 21:50
格点太密，就是平滑处理了效果也不太好。就像用 2.5×2.5 的再分析资料画出的高度场 比 1.0×1.0 再分析资 ...

恩恩，只能说尽可能让图好看一点，因为模式结果的格点很密，这也是没有办法的

FrankieLJY

想请问一下，是全国所有地区的等高线值WRF都有吗？都是在geo_em.d01这个数据里吗？

chsolomon

MeteoInfo 发表于 2017-1-3 10:34
有空了我加上平滑的函数吧，不过平滑会改变原始值。目前你可以对数组进行稀疏化处理（比如隔4个点取一个值 ...

老师，我用的wrf的grib2的数据，也是想稀疏化处理，附上我的GS
'reinit'
'open d:\wrf1km\wrfprs_d01.2017080212.18.ctl'
'set mpdset cnworld'
'set lon 105 140'
'set lat 24 50'
'set grid off'
'set csmooth on'
'set cterp on'
'PRESsfc=PRESsfc[::3,::3]'
'set cint 50'
'd PRESsfc/100'
'print'
'printim d:\wrf1km\xs\dm\0212.gif white'
*'disable print'
但是不成功，是不是数据的维数有问题。