关于用grads2ascii将grads计算结果输出为txt的问题

liuzf406

刚看了兰溪的贴子, 如何用grads2ascii 将grads计算结果输出为txt (http://bbs.06climate.com/forum.php?mod=viewthread&tid=4764&extra=&page=1), 有个问题不太明白. 如果我要输出的数据是不同格点的, 里面是否包括经纬度信息啊?
刚也看了清风的贴子:(http://bbs.06climate.com/forum.php?mod=viewthread&tid=7357, 二、提取格点数据插值为站点数据的结果). 你那个结果很好, 但我不需要插值, 只想输出经度, 纬度, 变量值.  运行了一下grads2ascii, 我的输出结果如下面, 如何将经纬度和变量值对应起来呢? 谢谢!

-38.3378	-43.376	26.7141	28.2916	52.2796	-9.99E+08	-9.99E+08	-83.5929
-62.3533	-9.99E+08	-9.99E+08	-9.99E+08	-91.5844	9.29299	-33.5595	-20.0744
-63.1749	-73.5768	-27.563	-73.8197	-31.4969	-100.741	-71.9224	-66.7743
-9.99E+08	-51.7005	-9.99E+08	-44.1601	-48.2441	16.6732	13.3509	9.31587
18.9312	-16.6866	-17.147	-9.99E+08	-95.6972	-104.77	-9.99E+08	-45.7616
-78.1954	-23.9763	-13.8892	-1.64904	-90.35	-27.9558	-15.9525	-29.0194
115.983	-60.0411	-33.0559	-9.99E+08	-9.99E+08	-9.99E+08	-87.4122	-51.5873
-9.26044	-70.7823	-24.6502	78.073	13.7077	20.5778	-31.8791	5.68892
27.8533	-9.99E+08	-9.99E+08	-9.99E+08	-9.30971	-47.9302	-47.442	39.9804
-57.1412	-80.4545	-20.3351	-47.3948	-76.7288	-59.478	-17.1421	-52.7949
-49.0865	-9.99E+08	-40.8377	-37.8093	12.971	23.676	-19.6591	30.1537
15.1241	26.8828	10.0179	45.9156	16.4469	9.60464	-15.8429	-2.00067
-9.99E+08	24.6976	-81.7517	3.88164	-9.99E+08	-49.9524	-63.8028	-19.1774
-54.4036	-9.48039	-10.9256	-22.4313	-6.52395	26.7909	-32.5655	-38.7375
-9.99E+08	2.42511	5.73159	-74.2469	31.2473	22.5957	33.9881	27.1418
14.8415	-7.70484	-19.9919	-8.46245	-10.7335	-110.487	-50.8192	-22.3071
-71.9879	-61.247	-69.8177	-4.66317	-26.4542	-10.5157	-32.1769	-17.8654
-32.9285	-25.8344	-51.9634	-45.7308	-26.4005	-10.909	-35.0275	-13.8601
-52.3206	26.7932	25.6668	17.9571	18.6333	-5.28089	-32.5338	-83.7465
-60.6701	-139.24	-91.1389	-9.99E+08	-9.99E+08	24.6441	17.7584	-13.4419
-82.9134	-63.9637	-17.2323	-41.1879	18.8448	-27.2295	-55.4815	-41.6375

liao001

多谢。。。。。。。。。。。。。。。。。。。。。。。。。。

mofangbao

liuzf406 发表于 2012-11-13 22:13

                               
登录/注册后可看大图

搞定了. 你那个输出结果里面包括经纬度信息, 但我的里面没有.
分别输出经度, 纬度 和 变量值即可, 如下: ...

合适的话，把gs贴出来让大家参考下吧，后面的人也许有同样的疑问~谢谢了~

mofangbao

不需要插值就别插值啊，你这个是格点数据了，我那个帖子里面不就有站点输出的例子么？只是因为我拿了同一份数据做例子，来回插值了一下，你要看本质啊！

liuzf406

mofangbao 发表于 2012-11-13 21:06

                               
登录/注册后可看大图

不需要插值就别插值啊，你这个是格点数据了，我那个帖子里面不就有站点输出的例子么？只是因为我拿了同一份 ...

是啊, 我就是按你那个站点输出的例子来操作的, 我的结果却不同啊, 看不出经度和纬度来啊.

随缘

你可以自己设置下输出格式的，最后两个参数，一个是列数，最后一个是间隔

liuzf406

随缘 发表于 2012-11-13 21:20

                               
登录/注册后可看大图

你可以自己设置下输出格式的，最后两个参数，一个是列数，最后一个是间隔

谢谢啊, 出来了向你汇报!

mofangbao

liuzf406 发表于 2012-11-13 21:18

                               
登录/注册后可看大图

是啊, 我就是按你那个站点输出的例子来操作的, 我的结果却不同啊, 看不出经度和纬度来啊.

我那个是插值过的，你直接输出原始变量就行了，怎么还没看懂。。。

liuzf406

mofangbao 发表于 2012-11-13 21:50

                               
登录/注册后可看大图

我那个是插值过的，你直接输出原始变量就行了，怎么还没看懂。。。

搞定了. 你那个输出结果里面包括经纬度信息, 但我的里面没有.
分别输出经度, 纬度 和 变量值即可, 如下:
15.00         70.00         -38.34
15.00         72.50         -43.38
15.00         75.00         26.71
15.00         77.50         28.29
15.00         80.00         52.28
15.00         82.50         -999000000.00
15.00         85.00         -999000000.00
15.00         87.50         -83.59
15.00         90.00         -62.35
15.00         92.50         -999000000.00
15.00         95.00         -999000000.00
15.00         97.50         -999000000.00
15.00         100.00         -91.58
15.00         102.50         9.29
15.00         105.00         -33.56
15.00         107.50         -20.07
15.00         110.00         -63.17
15.00         112.50         -73.58
15.00         115.00         -27.56
15.00         117.50         -73.82
15.00         120.00         -31.50
15.00         122.50         -100.74
15.00         125.00         -71.92
15.00         127.50         -66.77
............

liuzf406

mofangbao 发表于 2012-11-13 22:25

                               
登录/注册后可看大图

合适的话，把gs贴出来让大家参考下吧，后面的人也许有同样的疑问~谢谢了~

我就是分别输出了lat, lon 和 变量值, 各一列, 命令如下:
'grads2ascii lat c:\zhongfang\CNiso\cnpmiso\lat.txt %g 1'
'grads2ascii lon c:\zhongfang\CNiso\cnpmiso\lon.txt %g 1'
'grads2ascii wv_isotope c:\zhongfang\CNiso\cnpmiso\wv_isotope.txt %g 1'
就这样的的, 应该是正确的, 因为我用另外一种方法验证了一下. 这个方法是这里的老板写的一个小代码txtout, 结果是一样的!

谢谢你啊,清风大侠! 真的是在论坛上学了不少东西.

另外, 我将所有提取的格点数据(没插值的)输入arcgis, 用克里格进行了插值,做了个平滑, 其结果比grads里面插值看上去要好了很多.

river

liuzf406 发表于 2012-11-13 22:52

                               
登录/注册后可看大图

我就是分别输出了lat, lon 和 变量值, 各一列, 命令如下:
'grads2ascii lat c:\zhongfang\CNiso\cnpmiso ...

不能把这几个变量输出在一个txt里吗？