grads屏蔽不掉地图外的填色，求指导

guotingsaijin

大家好，我用grads画填色图，先用meteoinfo出的grads的maskout的底图分辨率0.05的格点画出来的图却填不满，由于meteoinfo出不了分辨率0.01以下的底图，我就用fortran来出分辨率0.005的格点背景场文件，同样用的一样的gs，画出来的却屏蔽不掉地图外的填图，有大神知道原因吗？用的一样的gs
第一张图是格点0.005背景场分辨率（fortran出的底图）的图，屏蔽不了底图外的。第二张图是格点0.05背景场分辨率的（meteoinfo出的底图）图，但是填不满地图，就想换分辨率更小的，却又屏蔽不了。
这是我的gs文件
'reinit'
'open H:\grads\BENGBU\grid.ctl'  
'open H:\grads\BENGBU\first.ctl'
'set mpdset BENGBU'
'set lon 116.5 118.5';'set lat 32.5 34'
'define a=oacres(mask.1,day.2,30)'
'define a1=maskout(a,mask.1-1.0)'
'define aa=smth9(a1)'
'set grads off'
'set grid off'
'set xlopts 1 6 0.18';'set ylopts 1 6 0.18'
'set xlint 1';'set ylint 1'
'set rgb 20 255 244 197'
'set rgb 21 255 227 101'
'set rgb 22 255 198 42'
'set rgb 23 255 167 42'
'set rgb 24 255 144 42'
'set rgb 25 255 100 42'
'set rgb 26 255 74 0'
'set rgb 27 227 32 23'
'set rgb 28 194 23 15'
'set rbcols 20 21 22 23 24 25 26 27 28 '
'set clevs 20 22 24 26 28 30 32 34'
'set gxout shaded'
'd aa'
'cbar 1.0 1'
'printim H:\grads\BENGBU\grid.png white'
;

river

0.05的分辨率应该足够了，你修改插值半径试试，往小了调调
'define a=oacres(mask.1,day.2,30)'

guotingsaijin

river 发表于 2025-1-13 12:51
0.05的分辨率应该足够了，你修改插值半径试试，往小了调调
'define a=oacres(mask.1,day.2,30)'

我调了还是不行

guotingsaijin

guotingsaijin 发表于 2025-1-13 14:55
我调了还是不行

'reinit'
'open H:\grads\BENGBU\bb0.05.ctl'  
'open H:\grads\BENGBU\first.ctl'
'set mpdset BENGBU'
'set lon 116.5 118.5';'set lat 32.5 34'
'define a=oacres(mask.1,day.2,10)'
'define a1=maskout(a,mask.1-1.0)'
'define aa=smth9(a1)'
'set grads off'
'set grid off'
'set xlopts 1 6 0.18';'set ylopts 1 6 0.18'
'set xlint 1';'set ylint 1'
'set rgb 20 255 244 197'
'set rgb 21 255 227 101'
'set rgb 22 255 198 42'
'set rgb 23 255 167 42'
'set rgb 24 255 144 42'
'set rgb 25 255 100 42'
'set rgb 26 255 74 0'
'set rgb 27 227 32 23'
'set rgb 28 194 23 15'
'set rbcols 20 21 22 23 24 25 26 27 28 '
'set clevs 20 22 24 26 28 30 32 34'
'set gxout shaded'
'd aa'
'cbar 1.0 1'
'printim H:\grads\BENGBU\bb.png white'
;[img]H:\GRADS\Bengbu\bb.pngH:\GRADS\Bengbu\bb.png[/img]

四叶草

Fortran出的背景场屏蔽底图外信息建议用province-basemap脚本，meteoinfo出的底图锯齿终究会存在，强迫症不建议用。

jiliguluwulala

maskout时候数据和mask的分辨率应该是要一样的，不知道你的a插值完后是不是0.005的分辨率

jiliguluwulala

另外我看meteoinfo可以设置分辨率为0.005，是输出的文件会出问题吗

guotingsaijin

jiliguluwulala 发表于 2025-1-14 11:00
另外我看meteoinfo可以设置分辨率为0.005，是输出的文件会出问题吗

是的，输出的文件名有点奇怪

river

guotingsaijin 发表于 2025-1-13 14:55
我调了还是不行

那你还是在论坛搜索province-basemap脚本吧，那个可以做到严丝合缝的