清风的双线性插值方法，为什么得到的TXT全是-999.00的缺测值呢？

kaiwang44 · 发表于 2016-7-22 10:47:46

好高深，看不懂

lqouc · 发表于 2016-7-22 12:50:04

经纬度的顺序反一下
就这句read(27,*) stid,newPoint(i,1),newPoint(i,2)按我写的改了试试
调程序别光动嘴，动手试试，到现在为止，你也没有把你的站点数据格式贴出来，谁也没法确定是什么问题。只能说，清风的程序应该是不会出错的。

jolincai · 发表于 2016-7-22 18:39:40

lqouc 发表于 2016-7-22 12:50
经纬度的顺序反一下
就这句read(27,*) stid,newPoint(i,1),newPoint(i,2)按我写的改了试试
调程序别光动 ...

好的我试试，以下是我站点的站号和经纬度。
50353  126.65 51.72
50434  121.47 50.50
50527  119.75 49.22
50557  125.23 49.17
50564  127.35 49.43
50632  121.92 48.77
50658  125.88 48.05
50727  119.95 47.17
50745  123.92 47.38
50756  126.97 47.43
50788  131.98 47.23
50854  125.32 46.38
50915  116.97 45.52
50949  124.83 45.12
50953  126.62 45.68
50963  128.73 45.97
50968  127.97 45.22
50978  130.95 45.28
51076 88.08 47.73
51087 89.52 46.98
51156 85.72 46.78
51243 84.85 45.60
51334 82.90 44.62
51379 89.57 44.02
51431 81.33 43.95
51463 87.62 43.78
51573 89.20 42.93
51644 82.95 41.72
51709 75.98 39.47
51716 78.57 39.80
51765 87.70 40.63
51777 88.17 39.03
51811 77.27 38.43
51828 79.93 37.13
52203 93.52 42.82
52418 94.68 40.15
52436 97.03 40.27
52495  104.50 40.75
52533 98.48 39.77
52602 93.38 38.83
52681  103.08 38.63
52713 95.37 37.85
52754  100.13 37.33
52787  102.87 37.20
52818 94.90 36.42
52836 98.10 36.30
52866  101.77 36.62
52884  103.93 36.35
52955  100.75 35.58
53068  112.00 43.65
53192  114.95 44.02
53276  112.90 42.40
53336  108.52 41.57
53352  110.43 41.70
53391  114.00 41.90
53463  111.68 40.82
53487  113.33 40.10
53502  105.75 39.78
53529  107.98 39.10
53614  106.22 38.48
53646  109.70 38.23
53673  112.70 38.75
53698  114.43 38.07
53723  107.40 37.78
53772  112.55 37.78
53845  109.50 36.60
53863  111.93 37.05
53898  114.37 36.12
53915  106.67 35.55
53959  111.02 35.03
54012  117.60 44.58
54026  120.90 44.57
54027  119.40 43.98
54094  129.60 44.57
54096  131.15 44.38
54102  116.07 43.95
54115  118.07 43.60
54135  122.27 43.60
54157  124.33 43.18
54161  125.22 43.90
54208  116.47 42.18
54218  118.97 42.27
54236  122.53 42.42
54292  129.47 42.88
54324  120.45 41.55
54337  121.12 41.13
54342  123.43 41.77
54346  123.78 41.32
54374  126.92 41.72
54405  115.50 40.40
54423  117.83 40.97
54471  122.27 40.67
54497  124.33 40.05
54511  116.47 39.80
54527  117.17 39.10
54539  118.90 39.42
54618  116.27 38.03
54662  121.63 38.90
54725  117.53 37.50
54776  122.68 37.40
54823  116.98 36.68
54843  119.08 36.70
54909  115.58 35.07
54916  116.85 35.57
55591 91.13 29.67
56004 92.43 34.22
56021 95.48 34.55
56029 97.02 33.02
56033 98.22 34.92
56046 99.65 33.75
56080  102.90 35.00
56096  104.92 33.40
56146 99.98 31.63
56172  102.48 31.67
56182  103.57 32.65
56187  103.83 30.70
56257  100.27 30.00
56444 98.90 28.50
56462  101.55 28.98
56492  104.53 28.82
56571  102.30 27.88
56651  100.43 26.87
56671  102.25 26.68
56739 98.50 25.02
56768  101.53 25.02
56778  102.68 25.02
56951  100.22 23.95
56954 99.93 22.57
56964  101.40 22.67
56985  103.38 23.38
57014  105.88 34.55
57067  111.02 34.00
57083  113.65 34.72
57127  107.03 33.07
57237  108.03 32.07
57265  111.67 32.38
57290  114.05 32.98
57297  114.08 32.12
57411  106.08 30.80
57447  109.37 30.27
57461  111.08 30.70
57494  114.07 30.63
57516  106.47 29.58
57633  108.77 28.80
57662  111.68 29.05
57687  112.92 28.22
57707  105.23 27.30
57713  106.88 27.70
57745  109.68 27.45
57799  114.97 27.12
57816  106.72 26.58
57866  111.62 26.23
57902  105.18 25.43
57957  110.30 25.33
57993  114.95 25.85
58027  117.30 34.28
58040  119.13 34.83
58102  115.78 33.88
58221  117.37 32.95
58238  118.80 32.00
58251  120.30 32.85
58314  116.33 31.40
58321  117.28 31.85
58362  121.47 31.42
58424  117.05 30.53
58457  120.17 30.23
58477  122.12 30.03
58527  117.20 29.30
58606  115.92 28.60
58633  118.87 28.97
58715  116.65 27.58
58752  120.62 27.80
58834  118.17 26.65
58847  119.28 26.08
58921  117.35 25.97
59023  108.05 24.70
59082  113.58 24.80
59134  118.07 24.45
59211  106.53 23.92
59254  110.08 23.40
59265  111.30 23.48
59287  113.32 23.13
59293  114.68 23.73
59316  116.68 23.40
59417  106.75 22.37
59431  108.35 22.82
59501  115.37 22.78
59632  108.60 21.95
59663  111.97 21.87
59758  110.35 20.03
59838  108.62 19.10
59855  110.47 19.23

jolincai · 发表于 2016-7-27 12:29:29

lqouc 发表于 2016-7-22 12:50
经纬度的顺序反一下
就这句read(27,*) stid,newPoint(i,1),newPoint(i,2)按我写的改了试试
调程序别光动 ...

大神前几天回家了一直耽误着没上网，嗯程序确实把经纬度的赋值弄反了，我按您的意思把读取经纬度的顺序换了一下，生成的TXT除了后十几个站点有正常的值外，前一百多个站点的值依然还是缺测-999.0的值。。。这是什么原因呢。。这是后几十个站点的值，前面全是-999.0
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  0.1081238
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  0.1081238
  -999.0000
  -999.0000
  -999.0000
  -999.0000
  -999.0000
2.443799
  -999.0000
  -999.0000
3.523591
  -999.0000
4.387475
4.532664
  -999.0000
  -999.0000
1.305725
  -999.0000
3.327590
3.058903
5.630439
4.686519
4.586181
  0.7537445
8.422973
4.030970
4.813225
5.201714
5.677449
5.550224
10.46564
4.615253
4.617191
10.94399
7.024707
8.326277
16.18757
16.66358
17.59061

lqouc · 发表于 2016-7-27 12:50:04

jolincai 发表于 2016-7-27 12:29
大神前几天回家了一直耽误着没上网，嗯程序确实把经纬度的赋值弄反了，我按您的意思把读取经纬度的顺序 ...

网格数据都是有数值的么？能插值的那些站点结果大概对不对？
提问尽量先把有用的信息都贴出来，那么长的程序，很难直接看出来哪里有问题。

zhunuo1990 · 发表于 2016-7-27 14:33:54

赞！赞！赞！

jolincai · 发表于 2016-7-28 10:18:58

lqouc 发表于 2016-7-27 12:50
网格数据都是有数值的么？能插值的那些站点结果大概对不对？
提问尽量先把有用的信息都贴出来，那么长的 ...

嗯网格点都是有数据的，能插值出来到站点的这些点数据是：
  3.327590
3.058903
5.630439
4.686519
4.586181
  0.7537445
8.422973
4.030970
4.813225
5.201714
5.677449
5.550224
10.46564
4.615253
4.617191
10.94399
7.024707
8.326277
16.18757
16.66358
17.59061
对应的站点数据是：
2.71
3.60
5.90
3.04
6.19
4.07
7.63
7.63
7.33
5.50
7.51
6.12
8.66
8.00
6.80
9.57
8.00
9.01
13.15
13.97
13.55
我是做资料对比，所以要相互插值。这个插值结果数值上有出入，不知道算不算结果正确。

jolincai · 发表于 2016-7-28 10:20:05

jolincai 发表于 2016-7-28 10:18
嗯网格点都是有数据的，能插值出来到站点的这些点数据是：
3.327590
3.058903

上面是插值结果，下面是对应站点的原始数据。

lqouc · 发表于 2016-7-28 10:53:18

jolincai 发表于 2016-7-28 10:20
上面是插值结果，下面是对应站点的原始数据。

......我问的是你差值的结果和网格数据附近位置上的值是否相符。
传一份网格数据和站位数据上来吧。

jolincai · 发表于 2016-7-28 22:06:01

一个是192个站点的站号经纬度一个是欧洲中心2.5*2.5 从70°E到140 ， 15°N到55的29*17的格点数据。

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

[求助] 清风的双线性插值方法，为什么得到的TXT全是-999.00的缺测值呢？