province-basemap抠图失败

breeze_li

自己用陕西省shp文件制作的 grads map文件和shannx_out.txt（边界数据），使用province-basemap命令画图，结果抠图失败，图中地图边界也是有的，请问有人遇到类似情况吗？（在别的软件测试，out边界文件应该是正常的），补充说明：我用之前grads自带的out文件和map文件绘图是正常的，说明脚本是没问题的，后面客户要求shp文件使用CGCS2000坐标系，不清楚之前的是什么坐标系，重新用arcgis处理了shp文件，然后重新用meteoinfo制作的map文件和out文件，同样的代码绘的图，图一是自带的地图数据，图二是后边制作的（之前同样的方法制作地图文件也是成功的）。
gs文件如下：

1. 'reinit'
   
2. 'open D:/www/cgi-bin/610000.ctl'
   
3. 'open D:/www/cgi-bin/610000sta.ctl'
   
4. 'set parea 0.5 7.5 0.3 9.7'
   
5. 'set grads off'
   
6. 'set grid off'
   
7. 'set lon 105 112'
   
8. 'set lat 31 40'
   
9. 'set xlopts -20'
   
10. 'set ylopts -20'
    
11. 'define rgrid=smth9(oacres(grid,t.2,30,20,15,10,5,1))'
    
12. 'set gxout shaded'
    
13. 'D:/www/ana/color/TEM.gs'
    
14. 'd rgrid'
    
15. 'set gxout contour'
    
16. 'D:/www/ana/color/TEM.gs'
    
17. 'set ccols 0'
    
18. 'set ccolor 0'
    
19. 'set clopts 1'
    
20. 'set clab masked'
    
21. 'set gxout contour'
    
22. 'province-basemap shannx rgrid'
    
23. 'D:/www/ana/city-cntyName/610000.gs'
    
24. 'cbar_interp 1 1 1 7.4787 5'
    
25. 'printim D:/Rending/temp.png white'
    
26. ;

复制代码

附图：

river

这个有论坛里有很多，不需要自己制作shannx_out.txt文件了。而且你的脚本也没有上传，也有可能是用法有问题

breeze_li

river 发表于 2018-8-18 19:49
这个有论坛里有很多，不需要自己制作shannx_out.txt文件了。而且你的脚本也没有上传，也有可能是用法有问题

你好，我更新了一下内容，代码附上了，麻烦有空看看

river

你这个后面直接是等值线也没画出来，我觉得是后面的这个province-basemap命令没有发挥作用。你改一下写法试试
'run province-basemap.gs shannx rgrid'

breeze_li

river 发表于 2018-8-18 22:57
你这个后面直接是等值线也没画出来，我觉得是后面的这个province-basemap命令没有发挥作用。你改一下写法试 ...

这个等值线没出来是我画图的时候，把那一部分代码注释了

四叶草

这个效果，out文件就没有起作用！用以前的out文件也不行噶？还是一定要新做的？

breeze_li

四叶草 发表于 2018-8-19 18:22
这个效果，out文件就没有起作用！用以前的out文件也不行噶？还是一定要新做的？

以前的是好的，新的是用arcgis转换坐标系了，说是要CGCS2000坐标系

四叶草

breeze_li 发表于 2018-8-19 18:42
以前的是好的，新的是用arcgis转换坐标系了，说是要CGCS2000坐标系

新坐标系的底图都是正常的，你用以前的out文件试试嘛。如果正常，新做的out文件有问题的嫌疑大

breeze_li

四叶草 发表于 2018-8-19 18:55
新坐标系的底图都是正常的，你用以前的out文件试试嘛。如果正常，新做的out文件有问题的嫌疑大

起初也怀疑是out文件的问题，因为用新的map文件和旧的out文件绘图也正常，后边把out文件的经纬度写成geojson文件，在openlayers中画的边界线也看着没问题，所以不知道问题出在哪，至于坐标系的问题，grads之前的坐标系用的是WSG84吗？

四叶草

breeze_li 发表于 2018-8-19 19:10
起初也怀疑是out文件的问题，因为用新的map文件和旧的out文件绘图也正常，后边把out文件的经纬度写成geoj ...

哈哈，你这个太高深了，不懂。