grads绘制槽线的简单方法

zhi726

由于在工作用经常需要在grads中绘制槽线，笔者想到一个简单的办法：利用贝塞尔曲线来绘制槽线，使用贝塞尔曲线最简单的二次方公式，由人工确定三点位置，再根据公式绘制槽线，槽线的起点和终点都在槽线上，但中间点一般不在槽线上，而是起到调整槽线曲率的作用，gs如下：
'reinit'
'open C:\dcz\20160627\fnl\fnldata\fnl.0627.ctl'
'enable print C:\dcz\20160627\fnl\hgt500.wind.gmf'
'set rgb 22 153 96 53'

'set lon 90 125'
'set lat 20 50'
'set t 1'
'set grads off'
'set grid off'
'set gxout contour'
'set lev 500'
'set cint 2'
'set ccolor 4'
'set cthick 6'
'set clab forced'
'set clopts 4 4 0.12'
'set clskip 2'
'set csmooth on'
'define hh=smth9(HGTprs/10)'
'd hh'
'set clevs 588'
'set ccolor 4'
'set cthick 12'
'd hh'
*******draw Trough line**********
*******need to locate three dots of the trough position, the first and the last is on the trough****
'q w2xy 108 33'
ln1=subwrd(result,3)
lt1=subwrd(result,6)
'q w2xy 106 30'
ln2=subwrd(result,3)
lt2=subwrd(result,6)
'q w2xy 101 28.5'
ln3=subwrd(result,3)
lt3=subwrd(result,6)
'set line 22 1 12'
rc = bezier(ln1,lt1,ln2,lt2,ln3,lt3)
'print'
'disable print'
'c'
'reinit'

function bezier(x1,y1,x2,y2,x3,y3)
   bt=0
   bt1=0
   bt2=0
   dt=0.1
   while(bt<=1-dt)
      bt1=bt
      bt2=bt+dt
      xs=(1-bt1)*(1-bt1)*x1 + 2*bt1*(1-bt1)*x2 + bt1*bt1*x3
      ys=(1-bt1)*(1-bt1)*y1 + 2*bt1*(1-bt1)*y2 + bt1*bt1*y3
      xe=(1-bt2)*(1-bt2)*x1 + 2*bt2*(1-bt2)*x2 + bt2*bt2*x3
      ye=(1-bt2)*(1-bt2)*y1 + 2*bt2*(1-bt2)*y2 + bt2*bt2*y3
      'draw line 'xs' 'ys' 'xe' 'ye''
      say 't= 'bt2''
      bt=bt+dt
   endwhile
return
*******************************

feiji158

厉害~~能全自动就更好了~

nunu18

看槽线的位置和风场有出入啊，出图以后直接在画图板画上可能最准确

rabin_xu

feiji158 发表于 2018-1-23 12:59
厉害~~能全自动就更好了~

能自动的是micaps吧

weinihou

确实和风场有出入诶

我这种非技术党都会优先考虑，出底图后PPT之类的手动加线

mayifan18

谢谢，值得学习，谢谢

小其其格

值得鼓励和借鉴，我现在后期都是再用photoshop绘制槽线切变线

zhi726

nunu18 发表于 2018-1-23 13:55
看槽线的位置和风场有出入啊，出图以后直接在画图板画上可能最准确

其实图片只是试验，因为槽线的起点和终点是人工确定的，所以可以修改它的终点就ok了

nunu18

zhi726 发表于 2018-1-24 00:13
其实图片只是试验，因为槽线的起点和终点是人工确定的，所以可以修改它的终点就ok了

哈哈~~~~厉害

親親ˇ/亲亲

感谢分享！