GrADS

rookie

作为一个新人，有些关于GrADS问题向各位请教下！
ctl文件：
PDEF  90  90  LCC  36.000  105.000  .......................
XDEF  270  LINEAR  88.55157  0.12162162
YDEF  200  LINEAR  23.49021  0.12162162


1、ctl文件中PDEF变量中有当前网格的行列以及中心经纬度坐标，XDEF、YDEF中按照我自己的摸索和理解是将90*90的网格扩展为270*200，XDEF和YDEF组成的大网格包含了90*90的小网格，那么对于XDEF*YDEF网格中有些值是依据90*90网格中的值经过插值得到，那么是如何插值的，或者源码有没有关于插值方法的文件。
听一个使用GrADS的介绍说插值使用的双线性插值，是否和以下的插值一样
http://baike.baidu.com/link?url=GDgL2TKdCrhdqNncQjSS8Ck1hpPta6kCIkwgcrDjqkKSpNfcOw0V_407qqjpo8oW_mMJF_QQVFz0TyCd4hGkw_

river

给个例子，楼主看一下
pdef 130 120 lcc 31.0 110. 50.0 63.0 30. 60. 110. 30000. 30000.
原始数据水平有130×120格点，lcc:采用Lambert投影；模式的参考点在北纬
31度，东经110度；该点对应的网格点坐标为水平第50个格点，y方向第63格点；
Lambert投影对应的两个标准纬度为北纬30度和60度，标准经度为东经110度；
网格距为（x/y）30000M×30000M。
以下定义标准的经纬度网格，网格的范围间距由用户自己定义，可以与原
网格无关。只要能保证原数据能落在网格范围内即可。

rookie

river 发表于 2015-9-7 20:03
给个例子，楼主看一下
pdef 130 120 lcc 31.0 110. 50.0 63.0 30. 60. 110. 30000. 30000.
原始数据水平 ...

感谢river对pdef的详细说明，我后面查阅的GrADS的说明文档，对我的问题有所理解。90*90插值为270*200网格，按文档上说明，针对270*200网格中的点，先映射到90*90网格上，然后找到周围4个点，做双线性插值。
具体正在尝试。

river

rookie 发表于 2015-9-8 15:06
感谢river对pdef的详细说明，我后面查阅的GrADS的说明文档，对我的问题有所理解。90*90插值为270*200网格 ...

不客气，这个确实是GrADS的一个难点，一般用的少的人，都不太懂

rookie

river 发表于 2015-9-8 16:53
不客气，这个确实是GrADS的一个难点，一般用的少的人，都不太懂

river，您好，还有个问题向问下您。原始网格为80*80大小 ，插值为200*200，由于80*80为lcc投影，插值后再图片上有数值的区域为扇形，也就是原来80*80网格的范围。
参考GrADS说明文档，插值采用双线性插值：
Bi-linear interpolation is used within this grid box to get down to the         position 31.24,67.88. The interpolation is linear within the i,j grid.


那么，如果我想获取 （106.23,38.44）点的值，GrADS是先在80*80网格中找到离（106.23,38.44）最近的周围四个点，然后用这四个点采用双线性插值获取（106.23,38.44）点的值。
不知道我的理解对不对？
还有插值时位置用网格点的行列号还是经纬度值？


PDEF  80  80  LCC  36.000  105.000  40.500  40.500  46.00000  26.00000  105.00000  30000.000  30000.000
XDEF  200  LINEAR  88.55157  0.12162162
YDEF  200  LINEAR  23.49021  0.12162162

river

rookie 发表于 2015-9-10 17:59
river，您好，还有个问题向问下您。原始网格为80*80大小 ，插值为200*200，由于80*80为lcc投影，插值后再 ...

太高级了，我觉得你的理解是对的

str@nger