grads绘制市县地图问题

breeze_li

看了几位前辈的帖子，http://bbs.06climate.com/forum.php?mod=viewthread&tid=21512（meteoinfo制作市县地图），http://bbs.06climate.com/forum.php?mod=viewthread&tid=2263&extra=&page=1（grads中shp地图使用），底图文件使用的是从帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=2754&extra=page%3D1上下载的shp文件，在动手做时出现了点问题，请各位帮忙看看问题出在哪？
这是我的数据：

数据描述文件：
DSET   E:/project/temp/sta.grd
DTYPE  station   
STNMAP E:/project/temp/sta.map
UNDEF  -9.99E33  
TITLE  Temperature Data Sample  
TDEF   1 linear 01may2011 1mo
VARS 1  
temper    0  99  Temper Data
ENDVARS

meteoinfo制作的西安市地图ctl文件：
DSET ^xian.dat
TITLE Mask data
UNDEF -9999.0
XDEF 31 LINEAR 107.0 0.1
YDEF 21 LINEAR 33.0 0.1
ZDEF 1 LINEAR 1.0 1.0
TDEF 1 linear 01may2011 1mo
VARS 1
  mask 0 99 background mask data (99)
ENDVARS

gs文件：
'reinit'  
'open e:/project/temp/sta.ctl'
'open e:/project/temp/xian.ctl'
'set mpdraw off'
'set lon 107 110'
'set lat 33 35'
'set gxout shaded'
'd maskout(oacres(mask.2,temper),mask.2(t=1))'
'draw shp e:/project/temp/xian.shp'
'set gxout contour'
'printim e:/project/temp/xian.png white'
;
不适用set gxout shaded语句与使用set gxout shaded语句的截图如下：


这个gxout shaded怎么回事，为什么填色附近有空白，而且覆盖了等直线？

独奏古筝到天明

set gxout shaded是画填色图的意思，必须要先shaded再contour，最好再看看手册哈。而且color bar 也没有画上

breeze_li

独奏古筝到天明 发表于 2016-9-6 16:01
set gxout shaded是画填色图的意思，必须要先shaded再contour，最好再看看手册哈。而且color bar 也没有画 ...

我知道这个是画填色图，我这个shaded是在contour之前的啊，color bar画不画不影响出图吧？

四叶草

没看出你要表达什么问题？你是想表达图锯齿太“丑”？？？还是另外的问题？从你的信息看，没什么问题。

breeze_li

四叶草 发表于 2016-9-6 16:56
没看出你要表达什么问题？你是想表达图锯齿太“丑”？？？还是另外的问题？从你的信息看，没什么问题。

这个shaded没有贴合地图边界，而且把等值线盖住了，想绘制之前用basemap做出来的那种效果，如下图：

四叶草

breeze_li 发表于 2016-9-6 17:06
这个shaded没有贴合地图边界，而且把等值线盖住了，想绘制之前用basemap做出来的那种效果，如下图：

你的做法是对的，没有遇到问题，仅仅是做的图和你预想的不一致。我大概懂你的意思了。
这个问题，家园里有很多帖子，说的也很多。你要贴合地图边界，就是选择province-basemap最给力（个人认为），边界有锯齿，属于正常现象（原因就是站点插值格点，格点就是方形的小区域，有锯齿很正常，这个和分辨率有关。）。
如果你用province-basemap，切记不要用meteoinfo制作的格点背景场，不然province-basemap的作用就体现不出来了。格点背景场你就自己编个小程序生成就行，有什么问题我们再交流。

四叶草

breeze_li 发表于 2016-9-6 17:06
这个shaded没有贴合地图边界，而且把等值线盖住了，想绘制之前用basemap做出来的那种效果，如下图：

http://bbs.06climate.com/forum.php?mod=viewthread&tid=37504
估计你是想要绘出这个帖子里的图，边界光滑无锯齿。

breeze_li

四叶草 发表于 2016-9-6 17:19
你的做法是对的，没有遇到问题，仅仅是做的图和你预想的不一致。我大概懂你的意思了。
这个问题，家园里 ...

之前使用province-basemap绘制的省地图，现在需要绘制市以及县的地图，不知道basemap该怎么处理边界

四叶草

breeze_li 发表于 2016-9-6 17:25
之前使用province-basemap绘制的省地图，现在需要绘制市以及县的地图，不知道basemap该怎么处理边界

可以绘制省地图，也可以绘制市地图，更可以绘制乡、村……只要你想要的都可以绘制。
只要你有边界数据，全部问题都能解决。
如果楼主有兴趣，直接看懂basemap脚本，那就什么都理解了。
边界数据来源：县级及以上，网上都能找到，乡村等的，那就自己用鼠标慢慢手工绘制，运用软件surfer和ARCGIS都行，我用的是ARCGIS。

breeze_li

四叶草 发表于 2016-9-6 17:31
可以绘制省地图，也可以绘制市地图，更可以绘制乡、村……只要你想要的都可以绘制。
只要你有边界数据， ...

像cnbasemap与province-basemap需要out文件，gs文件，市县需要哪种数据，该怎样用？看脚本手册，没有例子的话，对我这个小白来说理解起来有点吃力