很有趣 --还是置信度打点的问题

liuzf406

请大家看图，我的置信度区域可以用线绘出来，为何却无法打点输出，问题出在哪里？
显示错误为： Error in gashpwrt: mproj latlon required for gxout shapefile.
难道置信度打点不适用于mproj nps？
以下是显示图及打点的代码，请过来人指点：
'page q1'
'd anzpt'
'cbarn'
**打点**
'define tp=abs(tanzpt)-1.330'
'set gxout shp'
'set shpopts 15'
'set shp -pt shppt'
'set shpopts -1 2 0.001'
'd maskout(tanzpt,tp)'
'draw shp shppt'

*绘线*
'set gxout contour'
'set clab off'
'set  clevs 0'
'set cthick 5'
'set ccolor 15'
'd abs(tanzpt)-1.330'

'page q2'
'd anznt'
'cbarn'

**打点**
'define tn=abs(tanznt)-1.328'
'set gxout shp'
'set shpopts 15'
'set shp -pt shppt'
'set shpopts -1 2 0.001'
'd maskout(tanznt,tn)'
'draw shp shppt'

*绘线*
'set gxout contour'
'set clab off'
'set  clevs 0'
'set cthick 5'
'set ccolor 15'
'd abs(tanznt)-1.328'

liuzf406

river 发表于 2014-10-29 14:54
哦，这个gs我也有，但是从来没用过。楼主可以发帖分享一下经验啊

'drawmark hgt 1 1 0.05 0'
Usage: drawmark <var> <mark> <color> <size> [<magnitude> [<append> [<text>]]]]
      <var>: variable name.
      <mark>: mark type.
      <color>: mark color.
      <size>: mark size.
      <magnitude>: magnitude of <var> corresponding to <size>. If =0 then all marks will have the same size of <size>. Default=0.
      <append>: 0 (default) or 1. Set to 1 if appending to an existing plot. (Run "legend.gs" only once after all data are plotted.)
      <text>: Text to be shown in the legend (see "legend.gs"). Text beginning with a minus sign or containing spaces must be double quoted.

运行时需要qdims.gsf和parsestr.gsf 这两个辅助程序，均可在binguan 网页上下载。

liuzf406

牛人快来啊，我一直在等, 急！

river

看错误提示的意思是需要在输出成shp之前设置经纬度之类的

liuzf406

river 发表于 2014-10-29 09:21
看错误提示的意思是需要在输出成shp之前设置经纬度之类的

是投影的问题，好象只能为latlon投影才可以！

liuzf406

river 发表于 2014-10-29 09:21
看错误提示的意思是需要在输出成shp之前设置经纬度之类的

是投影的问题，好象只能为latlon投影才可以！

river

liuzf406 发表于 2014-10-29 09:47
是投影的问题，好象只能为latlon投影才可以！

那么你这个就不能用打点来表示置信区了

liuzf406

river 发表于 2014-10-29 12:18
那么你这个就不能用打点来表示置信区了

搞定了，用版上一位朋友发的drawmark来完成的。不错！

river

liuzf406 发表于 2014-10-29 14:06
搞定了，用版上一位朋友发的drawmark来完成的。不错！

哦，这个gs我也有，但是从来没用过。楼主可以发帖分享一下经验啊

river

liuzf406 发表于 2014-10-29 15:02
'drawmark hgt 1 1 0.05 0'
Usage: drawmark     [ [ []]]]
      : variable name.

你真实在，我也知道gs里面有用发。我的意思是你把你的整个画图的gs和出来图的效果什么的分享一下。毕竟论坛关于打点这一类的分享比较少