如何使我的画图区域和成都边界贴合

蜻蜓队长

lqouc 发表于 2017-12-21 10:42
上gs和ctl

看这个形状估计你是站点插值格点之后用了maskout，不要用maskout，自己做个basemap文件就行了 ...

我好想没有找见，可以具体一些吗？或者提供些链接，我找到的是如何利用MI制作适用于Grads的maskout文件

蜻蜓队长

lqouc 发表于 2017-12-21 10:42
上gs和ctl

看这个形状估计你是站点插值格点之后用了maskout，不要用maskout，自己做个basemap文件就行了 ...

我好像没有找见，可以具体一些吗？或者提供些链接，我找到的是如何利用MI制作适用于Grads的maskout文件

river

都说了让你上脚本······

蜻蜓队长

river 发表于 2017-12-23 09:56
都说了让你上脚本······

'reinit'
'open E:\*****\grid.ctl'
'open E:\*****\wstemp.ctl'
'open E:\*****\chengduwind.ctl'
'open E:\*****\chengdumask.ctl'
'set font 5'
'set grid off'
'set grads off'
'set parea 1 10 2 8'
'set lon 102 105'
'set lat 30 32'
'set z 1'
'set t  1'
********temperature******************
'set gxout shaded'
'set csmooth on'
'define a=oacres(g,temp.2)'
*'define a1=maskout(a,g-0.05)'
*'define aa=smth9(a1)'
'd smth9(maskout(a,mask.4))'
********wind speed******************
'set gxout contour'
'set csmooth on'
'define b=oacres(g,ws.2)'
*'define b1=maskout(b,g-0.05)'
*'define bb=smth9(b1)'
'd smth9(maskout(b,mask.4))'
********wind vector******************
'set gxout vector'
'drawvect u.3 v.3 br 0.2 2 1 6 1 0.1 0 m/s'
********shijian******************
'q time'
time1=subwrd(result,3)
'draw title 'time1''
********.shpwenjian******************
'draw shp E:\*****\chengdugezhanming.shp'
'draw shp E:\*****\dijishi_2004.shp'
' cbarn 1 0 5.4 0.5'
'printim E:\*****\111.png white'
'reinit'
;

蜻蜓队长

river 发表于 2017-12-23 09:56
都说了让你上脚本······

还要ctl文件吗？

river

蜻蜓队长 发表于 2017-12-23 10:03
还要ctl文件吗？

贴子里就是用这个脚本画出来的？试试调整一下插值半径，类似于这样'define a=oacres(g,temp.2,50,25,10,5)'

蜻蜓队长

river 发表于 2017-12-23 10:31
贴子里就是用这个脚本画出来的？试试调整一下插值半径，类似于这样'define a=oacres(g,temp.2,50,25,10,5 ...

是的呢，我想知道这个差值函数的差值半径选取有没有什么讲究呢？在论坛上搜索了很多，都没有详细介绍oacres函数差值原理的帖子

river

蜻蜓队长 发表于 2017-12-23 10:36
是的呢，我想知道这个差值函数的差值半径选取有没有什么讲究呢？在论坛上搜索了很多，都没有详细介绍oacr ...

为啥要在论坛搜呢，这是grads的函数，去看官网介绍会详细有用得多

oacres()

oacres (gexpr, sexpr <,radii>)

A Cressman objective analysis is performed on the station data to arrive at a gridded result that represents the station data.

gexpr   a valid grid expression
sexpr   a valid station data expression
radii   optional radii of influence. Defaults are: 10,7,4,2,1
The Cressman Analysis scheme is described in a 1959 paper in Monthly Weather Review. Multiple passes are made through the grid with increasingly smaller radii of influence. At each pass, a new value is calculated for each grid point based on a correction factor that is determined by looking at each station within the radius of influence.

For each such station, an error is defined as the difference between the station value and a value arrived by interpolation from the grid to that station. The correction factor is based on a distance weighted formula applied to all such errors within the radius of influence. The correction factors are applied to each grid point before the next pass is made.

Any grid boxes that do not have stations within the third specified radius are set to the missing data value.

Usage Notes

The oacres function can be quite slow to execute, depending on grid and station data density.
The actual values of the gridded expression are ignored, but the grid itself is used as a template to perform the analysis. The scaling of the grid must be linear in lat-lon.
The Cressman Analysis scheme can be unstable if the grid density is substantially higher than the station data density (ie, far more grid points than station data points). In such cases, the analysis can produce extrema in the grid values that are not realistic. It is thus suggested that you examine the results of oacres and compare them to the station data to insure they meet your needs.
Objective analysis is a complex topic, and many schemes for doing it have been developed over the years. The oacres function is more of a quick-look function rather than a rigorous analysis scheme. If you have specific analysis requirements, consider doing your objective analysis outside of GrADS with a user-defined function.
See the related function oabin.
Examples

The simplest case:
oacres(ts,ts.2)
To specify your own radii of influence:
oacres(ts,ts.2,12,8,5,4,3,2,1)

蜻蜓队长

river 发表于 2017-12-25 22:17
为啥要在论坛搜呢，这是grads的函数，去看官网介绍会详细有用得多

oacres()

谢谢river，我把这招给忘记了

蜻蜓队长

lqouc 发表于 2017-12-21 10:42
上gs和ctl

看这个形状估计你是站点插值格点之后用了maskout，不要用maskout，自己做个basemap文件就行了 ...

想问下basemap文件中的out文件排列有什么要求吗？除了首尾要添加一些序列