置信度打点，点太密，如何处理？

liuzf406

各位高手，在置信度打点时，点太密，想用skip命令处理，但似乎这个命令不起作用，请教如何处理？

liuzf406

river 发表于 2014-10-27 15:24
试试这个命令，最好把你的gs也贴上来分享下吧，图挺好看
set shpopts

多谢river朋友，我是想算趋势，然后用线性回归方法来检验其显著性，显著的区域打点，gs文件如下，欢迎大家指点：
'run rgbset2'
'reinit'
'set grid on'
'set grads off'

'open year.ctl'
'sdfopen NA1901.2011.tmp.djf.nc'
'set grid on'
'set grads off'
'set poli on'


'set poli on'
'set mpdraw on'

'set vpage 0 11 0 8.5'

'set parea 1 10 1 8'
'set mpdset hires'
'set frame circle'
'set mproj nps'

'set lon -140 -60'
'set lat  15 75'
'set xlint 10'
'set ylint 10'

'set ylopts 1 4 0.20'
'set xlopts 1 4 0.20'

'set gxout shaded'
'set csmooth on'

'set rbcols 49 48 47 46 45 44 43 42 22 23 24 25 26 27 28 29'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'


'set t 1 110'
'ltrend djft.2 air slope rmse'
'set t 1'
'd slope*10'
'run cbarn'
'run gmt.gs 0.1 10 10'

'set x 1'
'set y 1'
'set z 1'
'set t 1 110'
'define year= aave(var.1,global)'
'set lon -140 -60'
'set lat  15 75'

'set z 1'
'set t 1'
'define rr=tcorr(year,djft.2, t=1, t=110)'
*'d rr'
********************************test the significance of correlation correlation coefficient （rr）
* degree of freedom=110-2=108 (from 1902 to 2011), 95% significance (both sides)=0.195
* degree of freedom=110-2=108 (from 1902 to 2011), 90% significance (both sides)=0.164

'define tt=rr*sqrt(108/(1-pow(rr,2)))'
'define ts=abs(tt)-0.195'

'set gxout shp'
'set shp -pt shppt'
'd skip(maskout(tt,ts),2,2)'
'set shpopts -1 2 0.001'
'draw shp shppt'
'drawbox -95 -75 30 40'
'drawbox -130 -100 45 65'

river

试试这个命令，最好把你的gs也贴上来分享下吧，图挺好看
set shpopts

set shpopts fillpoly <marktype> <marksize>

Sets current attributes for drawing shapefiles.

fillpoly   color number to use when drawing filled polygon elements of a shapefile. Default is -1, which draws polygons unfilled.
marktype   Mark type to use when drawing point elements of a shapefile. Please see the draw mark command for the list of options.
marksize   Mark size of use when drawing point elements of a shapefile.

Usage Notes

This command is available with GrADS version 2.0.a8 or later.

When drawing shapefiles that contain points, the default mark type is a closed circle (type 3) and the default size is 0.05.

When drawing shapefiles that contain polygons, the default behavior is to draw only the perimeter of each polygon element. Use the fillpoly option with set shpopts to draw filled polygons and set the fill color. The polygon perimeters will also be drawn when the fillpoly option is used. The color, style, and thickness of the polygon perimeters are controlled by the set line command.

Use the q shpopts command to see the current settings for drawing and writing shapefiles.

Please see the documentation page on shapefiles for more details.

Examples

set shpopts 15

liuzf406

river 发表于 2014-10-27 15:24
试试这个命令，最好把你的gs也贴上来分享下吧，图挺好看
set shpopts

刚试了下这个命令，不行啊！

river

liuzf406 发表于 2014-10-27 20:15
刚试了下这个命令，不行啊！

你说的密我不太理解，每个点都能看清楚啊。还是你说的是通过检验的区域太大？

liuzf406

river 发表于 2014-10-27 20:27
你说的密我不太理解，每个点都能看清楚啊。还是你说的是通过检验的区域太大？

您好，我意思是在slope显著区域，这些点分布太密了，掩盖了slope的底图色，可不要以象skip命令那样，让它隔一个显示一个点？

river

liuzf406 发表于 2014-10-27 20:32
您好，我意思是在slope显著区域，这些点分布太密了，掩盖了slope的底图色，可不要以象skip命令那样，让它 ...

我个人感觉一点儿都不密，正好。实在不行你把点的颜色稍微换一下。

river

liuzf406 发表于 2014-10-27 20:32
您好，我意思是在slope显著区域，这些点分布太密了，掩盖了slope的底图色，可不要以象skip命令那样，让它 ...

你去看有一下官方的例子吧。毕竟新功能还没有多少人有深入研究http://grads.iges.org/grads/gadoc/shapefiles.html

liuzf406

river 发表于 2014-10-27 20:40
我个人感觉一点儿都不密，正好。实在不行你把点的颜色稍微换一下。

也是，换个颜色看看，新功能不少，但怕去研究，以前觉gmf格式的图很好，但现在却没法输出了。

river

liuzf406 发表于 2014-10-27 20:54
也是，换个颜色看看，新功能不少，但怕去研究，以前觉gmf格式的图很好，但现在却没法输出了。

现在可以直接出PS之类的图啊，设置一下分辨率什么的，都挺好