如何利用province-basemap绘制多行政区域（京津冀地区）

李文帅

在绘制京津冀地区变量时，利用'province-basemap heb oacres(g.2,pm25.1)'只能绘制河北地区而掩盖了北京和天津地区，如何解决呢？脚本及出图如下：
'reinit'
'open D:\Mine\Work\Program\2022.02.08_JJJ\grads_plot\monthly_201501.ctl'
'open D:\Mine\Work\Program\2022.02.08_JJJ\grads_plot\grid_0.5.ctl'

'set grid off'
'set grads off'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set map 1 1 2'
'set mpdset hires'

'set lon 113 120'
'set lat 35 43'
'set ylabs 35N|37N|39N|41N|43N'
'set xlabs 113E|115E|117E|119E|121E'
'set csmooth on'
'set gxout shaded'
'D:\Research_Room\Data\ThirdGroup\gs\colorbars\PM10_28colors.gs'
'set clevs 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210'
'set t 1'
'province-basemap heb oacres(g.2,pm25.1)'
*'cnbasemap oacres(g.2,pm25.1)'
*    'basemap O 15 0 M'
'set font 1'
'run cbar_interp 2 1 1 9 4.3'
*    run cbarn_interp jian'ge size vert/hori xmid ymid

'draw title 2015.01'
'printim D:\Mine\Work\2022.02.08_JJJ\pictures\distribution\2015.01_0.5.png white'   
;

river

两个方法：
1 修改beij heb tianj 三个文件形成一个京津冀地区的经纬度范围，然后使用province-basemap来出图
2 利用meteoinfo软件或者其他软件制作京津冀地区dat 、ctl文件形成掩膜，然后用maskout命令来出图。可参考@传说中的谁  的帖子 http://bbs.06climate.com/forum.php?mod=viewthread&tid=893

李文帅

river 发表于 2022-2-16 13:08
两个方法：
1 修改beij heb tianj 三个文件形成一个京津冀地区的经纬度范围，然后使用province-basemap来 ...

感谢