grads设置风矢量箭头颜色问题，set ccolor无法使用

happymgj

由于雷达反演风速的产品有一些瑕疵，想通过限制风速的大小来掩盖瑕疵产生的坏点（就是很多大值），如图一，查到可以用mag来实现，但使用这个命令之后，矢量箭头的颜色没有办法通过set ccolor进行控制，如图二，这该怎么解决呢？附上程序：
'reinit'
'open F:/MCS725/Radar/yi/dBz.ctl'
'open F:/MCS725/Radar/yi/U.ctl'
'open F:/MCS725/Radar/yi/V.ctl'

i = 1
while(i<=60)

'set xlopts 1 8 0.16'
'set ylopts 1 8 0.16'
'set lon 116 119.2'
'set lat 38.2 40.5'
'set t 'i''
'set z 6'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set map 15 1 0.5'

'set gxout shaded'
'set csmooth on'
'set rgb 21 1 160 246'
'set rgb 22 0 236 236'
'set rgb 23 0 216 0'
'set rgb 24 1 144 0'
'set rgb 25 255 255 0'
'set rgb 26 231 192 0'
'set rgb 27 255 144 0'
'set rgb 28 255 0 0'
'set rgb 29 214 0 0'
'set rgb 30 192 0 0'
'set rgb 31 255 0 240'
'set rgb 32 150 0 180'
'set rgb 33 173 144 240'
'set clevs 15 20 25 30 35 40 45 50 55 60 65 70'
'set ccols 0 22 23 24 25 26 27 28 29 30 31 32 33'
'd dbz.1'
'cbarn.gs'

'set gxout vector'
'set arrscl 0.4 12'
'set cthick 6'
'set arrowhead 0.04'
'set cmax 22'
'set cmin 0'
'set ccolor 1'
'd skip(U.2,5,5);V.3;mag(U.2,V.3)'

'q time'
x=subwrd(result,3)
'draw title 'x'UTC   3km'
'printim F:/MCS725/Radar/yi/test/'%i%'.png white'
'c'
i=i+1
endwhile
'disable print'
'reinit'

river

'd skip(U.2,5,5);V.3;mag(U.2,V.3)'
这个写法本身就是利用最后一个表达式给矢量箭头上色，来显示给矢量箭头的大小。确实没办法控制颜色了。
但是你这个可以利用另一个函数实现，那就是maskout，这个也是通过设置一个阀值，达到出图的时候不画大于阀值的数据的目的。你去看一下具体用法就会了。

happymgj

river 发表于 2017-8-21 21:28
'd skip(U.2,5,5);V.3;mag(U.2,V.3)'
这个写法本身就是利用最后一个表达式给矢量箭头上色，来显示给矢量箭 ...

嗯嗯，多谢提示，确实可以实现

luxj

请问一下画近三十年一月份平均风场时可以画出图但无法导出PNG图片是怎么回事呢？
以下是gs文件：
'clear'
'sdfopen D:\uwnd.mon.mean.nc'
'sdfopen D:\vwnd.mon.mean.nc'
'set lat -90 90'
'set lon -180 180'
'set lev 850'
'set t 25 384'
'set vpage 0 11 0 8.5'
'set parea 1 10.5 1 8'
'set gxout vector'
'aveu=ave(uwnd.1,t=25,t=373,12)'
'avev=ave(vwnd.2,t=25,t=373,12)'
'display skip(aveu,3,3);skip(avev,3,3)'
'draw title wind in 850 hpa'
'draw xlab lon'
'draw ylab lat'
'printim D:\11.png white'

兰北

可以通过set line 来进行设置

river

luxj 发表于 2018-4-21 00:03
请问一下画近三十年一月份平均风场时可以画出图但无法导出PNG图片是怎么回事呢？
以下是gs文件：
'clear' ...

这个问题应该已经解决了······

药草

set arrscl 0.4 12 这个表示的含义以及单位是多少啊

Linxinrui

药草 发表于 2022-8-15 10:16
set arrscl 0.4 12 这个表示的含义以及单位是多少啊

set arrscl size <magnitude>  - 设置矢量的长度和大小。选项size的值设定矢量的长度，其单位为英寸（虚页英寸），取值范围为0.5到1.0。选项magnitude的值设定矢量的大小，在其设定后矢量的长度将会随之适当的调整。当magnitude缺省时所有矢量同长

chenhongfengye

请问解决了吗？我用了maskout之后缺测出来的是灰色是什么原因？