画站点分布图

勇敢的悲伤

我要画江淮地区的站点分布图。结果出来的图只有地图，没有那些圆点，是GS文件哪儿出错了吗，高手们帮忙看看啊~~

GS:
'reinit'
'open f:\huatu\grid.ctl'
'open f:\huatu\r.ctl'
'set lon 110 124'
'set lat 27 35'
'set t 1'
'set mpdset hires'
'define a=oacres(g,r.2)'
'set grid off'
'set gxout contour'
'set cthick 10'
'set ccolor 1'
'set cmax -100'
'd a'
'q w2xy 34.29 110.05'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.48 111.12'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.03 111.02'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.50 112.26'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.38 112.28'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.47 111.36'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.43 113.39'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.01 113.51'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.46 114.23'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.18 111.30'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.02 112.35'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.53 113.03'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.47 114.31'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.00 110.25'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.51 110.49'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.03 110.44'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.23 111.40'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.09 112.45'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.00 114.01'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.08 114.03'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.04 110.24'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.10 112.34'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.37 113.49'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.11 115.01'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.47 110.14'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.11 110.41'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.42 111.18'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.52 111.05'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.20 112.11'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.40 113.10'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.37 114.08'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.24 110.10'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.35 111.22'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.22 112.24'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.59 113.55'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.23 113.05'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.02 114.35'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.28 110.24'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.03 111.41'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.23 111.13'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.51 112.22'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.12 113.05'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.43 113.34'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.13 112.55'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   27.52 113.10'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.27 115.40'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.25 116.20'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.17 117.09'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   34.50 119.07'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.52 115.46'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.38 116.59'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.59 118.31'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.38 119.01'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.36 119.02'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   33.46 120.15'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.55 115.49'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.10 115.40'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.33 116.47'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.57 117.23'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.18 118.18'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.00 118.48'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.48 119.27'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.52 120.19'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.01 120.51'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   32.04 121.36'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.45 116.30'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.24 116.19'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.04 116.57'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.52 117.14'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.37 117.52'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.20 118.22'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.09 118.35'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.46 119.56'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.26 119.29'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.04 120.26'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.24 121.29'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   31.10 121.26'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.44 115.40'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.15 115.03'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.32 117.03'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.37 118.59'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.08 118.09'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.21 119.25'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.14 120.10'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.37 121.05'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.16 121.13'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.44 122.27'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.44 122.49'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   30.02 122.06'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.44 116.00'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.35 115.59'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.00 116.41'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.18 117.12'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.43 118.17'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.37 119.01'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.07 119.39'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.36 120.49'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.52 121.34'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   29.12 121.57'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.36 115.55'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.04 115.33'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.18 117.13'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.58 118.52'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.41 118.15'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.27 119.55'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.05 119.08'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.49 120.55'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.00 120.40'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.51 121.08'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.37 121.25'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.27 121.54'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   28.05 121.16'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   27.57 117.50'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   27.55 118.32'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy   27.50 121.09'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'printim f:\huatu\r.png white'
;

river

论坛里很多讲这个的帖子，楼主搜一下看一下绝对会做了

勇敢的悲伤

river 发表于 2012-12-8 22:51

                               
登录/注册后可看大图

论坛里很多讲这个的帖子，楼主搜一下看一下绝对会做了

哈哈，弄好了。原来是我的经纬度放反了，所以导致黑点点出不来~~

mofangbao

你说你贴这么长的gs谁有心思给你从头看到尾呢？

Eudora99

应该是'q w2xy lon lat '吧，你好像是把经纬度的位置整反了

TrialBlazer

不需要打开格点文件，也不需要插值。直接打开站点文件就行了

风儿飘飘

如何将符合条件的点提取出来呢，请教

望晴檐

谢谢楼主的知识普及了~~

shirly

勇敢的悲伤 发表于 2012-12-9 09:08
哈哈，弄好了。原来是我的经纬度放反了，所以导致黑点点出不来~~

楼主你好，我现在在做和你一样的东西，出了问题，我看你gs里有下面一段，这是为了进行Cressman插值才写的是不是，我现在只需要把经纬度的地点在图中画出来看看分布，但老出错。我的gs除了没有下面这段程序之外，其它的都一样，只是经纬度最后保留的是4位数小数点。命令之间都是间隔一个空格这样的。如果知道是什么问题的话，还麻烦你指点一下~
'open f:\huatu\grid.ctl'
'open f:\huatu\r.ctl'
'set lon 110 124'
'set lat 27 35'
'set t 1'
'set mpdset hires'
'define a=oacres(g,r.2)'
'set grid off'
'set gxout contour'
'set cthick 10'
'set ccolor 1'
'set cmax -100'
'd a'

勇敢的悲伤

shirly 发表于 2014-3-27 17:11
楼主你好，我现在在做和你一样的东西，出了问题，我看你gs里有下面一段，这是为了进行Cressman插值才写的 ...

你现在写的这个只是插值到网格点上了吧。。如果你也是要把经纬度的地方放在图中，应该下面接着写我的那段。