请教数据筛选的Fortran程序

双木林 · 发表于 2012-6-12 09:42:54

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

用Anuspline插值出来的数据，想用Fortran写个程序提取中国范围内的插值结果，请教各位高手，程序怎么写？

mofangbao · 发表于 2012-6-12 10:59:44

你先自己想把，写的过程中有什么具体问题再讨论，思路就是先获得中国的经纬度范围，然后循环你的每一个点，判断你的点是否在这个经纬度范围之内，这个经纬度范围要是矩形才行，如果是不规则的（例如国界），而且是用来画图，那就没有必要提取了，否则你的插值算是白做了，因为插值完就是矩形，你再给提出一个不规则的，又成了站点数据，何必呢。
不知道你提取是为了什么目的，有必要说一下，如果是为了屏蔽画图的外部区域，那直接用其他方法屏蔽就好了。如果只是为了插值然后，看一下中国范围内的所有格点上的值，那么就是判断点是否在多边形内，这个算法在不少文献中会有，你也可以谷歌一下相关的算法。或者使用MeteoInfo来提取多边形内的点集，方法自己在MeteoInfo版块找一下

mofangbao · 发表于 2012-6-12 10:15:40

如果你是筛选不规则的点的话，筛选出来又是站点数据了，如果只是筛选一个矩形区域，那很好办，论坛的言深深发过相关的帖子，你搜索一下

topmad · 发表于 2012-6-12 10:23:10

格点插值还是站点呢？站点的话可以用站好进行筛选，格点的话，就像楼上所言，是个矩形框

双木林 · 发表于 2012-6-12 10:45:31

topmad 发表于 2012-6-12 10:23

登录/注册后可看大图

格点插值还是站点呢？站点的话可以用站好进行筛选，格点的话，就像楼上所言，是个矩形框

插值后的grid文件，是个矩形，但是中国范围内的格点是不规则的。

双木林 · 发表于 2012-6-12 13:30:10

mofangbao 发表于 2012-6-12 10:59

登录/注册后可看大图

你先自己想把，写的过程中有什么具体问题再讨论，思路就是先获得中国的经纬度范围，然后循环你的每一个点， ...

写了个循环判断每个点是否在中国范围的程序，费时，所以想是不是有省时间的方法。我是要提取中国范围的0.1度格点数据作为模型的输入。在ArcGIS里可以实现，但是一次做太多会出问题，所以想写个程序一次把几十年的数据一次处理了。

mofangbao · 发表于 2012-6-12 13:56:25

再省时的方法也要判断在不在中国范围内的吧，不判断怎么可能知道是舍弃还是保留呢，程序运行耗时是很正常的，如果超过了你能接受的程度可能是你程序设计上存在问题。

言深深 · 发表于 2012-6-13 10:16:02

没有必要单单筛选不规则区域，大气现象是没有国界的

		自动登录	找回密码
密码			立即注册

[求助] 请教数据筛选的Fortran程序

登录后查看更多精彩内容~

相关帖子