把.txt转化成.grd文件中确定格点数的问题

葡萄儿

用surfer画图时能识别的数据类型是.grd。我现在是.txt的文件。里面有站点的经纬度，但是以两位小数的形式存放的。我现在编个程序要把.txt转为，grd。里面的格点数怎么设置啊？

数据文件内容依次是：站名，经纬度，12个特征向量场，有27个站。

我看网上人家写的程序要定义：经纬度的格点数，层次，数据容量——我这里即是27*12。

因此可以要写成这样：parameter(nx=？,ny=？,nz=1,nt=360)

人家的nx和ny是这样算的：例如：经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度 2*2 ，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，根据数据的层次确定几层。
还有格点精度要怎么确定呢？

river

你的是站点资料，人家算的那个是格点资料。你首先要搞清楚surfer能识别怎样格式的站点资料，然后用fortran写成那样的就行了

nuister

楼上正解。楼主这是要画什么图啊，你要是想按照“人家”的方法计算就得把你的站点资料改成格点资料。

葡萄儿

river 发表于 2013-4-25 16:38

                               
登录/注册后可看大图

你的是站点资料，人家算的那个是格点资料。你首先要搞清楚surfer能识别怎样格式的站点资料，然后用fortran写 ...

人家的也是由经纬度算出来的格点数啊，从哪里判断已经是格点资料了呢？

葡萄儿

nuister 发表于 2013-4-25 18:45

                               
登录/注册后可看大图

楼上正解。楼主这是要画什么图啊，你要是想按照“人家”的方法计算就得把你的站点资料改成格点资料。

所以应该怎么换呢

nuister

葡萄儿 发表于 2013-4-25 22:08

                               
登录/注册后可看大图

所以应该怎么换呢

你有马红云那本grads软件基础教程么，第八章有详细的教程，一定要按照那个步骤来，先处理成二进制文件，再建立映射文件，然后生成格点文件，书上有很详细的说明。

nuister

不知道楼主要画什么样的图，个人觉得格点数的确定是和你站点所在范围和你要区的精度有关啊，你要用1*1的，格点数自然多，2.5*2.5的格点数就多。

葡萄儿

nuister 发表于 2013-4-27 09:07

                               
登录/注册后可看大图

不知道楼主要画什么样的图，个人觉得格点数的确定是和你站点所在范围和你要区的精度有关啊，你要用1*1的，格 ...

这个问题我懂了，谢谢啦

chengxf

葡萄儿 发表于 2013-4-27 09:39

                               
登录/注册后可看大图

这个问题我懂了，谢谢啦

你是懂了可我还没懂。Surfer中一个站点.txt文件转化成.grd文件中确定格点数的问题，还需要fortran程序来计算？你需要生成的grd文件的密度，也就是格点数在对话框中是可以设定的，很简单事情搞复杂了。对不起，也许是我没搞懂你的问题实质，说得不对请谅解。

葡萄儿

chengxf 发表于 2013-4-28 12:44

                               
登录/注册后可看大图

你是懂了可我还没懂。Surfer中一个站点.txt文件转化成.grd文件中确定格点数的问题，还需要fortran程序来计 ...

恩恩，你说的对。就是简单问题复杂化了。我之前没用过surfer，这一点看手册的时候没搞明白，其实就是软件自己插值。O(∩_∩)O谢谢你的热心，以后可能还要多多请教了~~~