【请教】怎么将站点资料缺测补全啊？

孤鹜

       有这么一份西南地区各旬土壤湿度资料：一共97个站，时间从07年到11年。问题在于资料有很多缺测的，例如A站第一个数据是08年5月第1旬，B站是从07年11月第1旬开始记录的，记录也有部分中断。那我怎么将缺测资料补全，并分别做各站旬平均呢？谢谢各位~~~~

言深深

哎呀，刚刚看了你的数据，其中有几列完全都是-9999，这样的数据肯定是舍弃的。
那些好的数据应该就是你说的含有缺测的列了，却测数值用-9999标记，有这么几种个人觉得合理的方案你可以选择性的使用：
1，用气候态平均代替，就是这一段较长时间的数据平均代替这个缺测值；
2，用缺测月份平均，比如3月份有一个数据缺测了，那么用统计的这个月份中的数据进行平均算这个缺测的；
3，用该天平均，历史上除去缺测这一天以外每一年对应月份这一天数据求和代替缺测值；
4，多点滑动，缺测点前后数据求平均代替这个数值。
个人推崇2,4两种方法

叩容

因为土壤湿度的变化是缓慢的，个人觉得楼上第4种方法更好。

孤鹜

言深深 发表于 2012-4-18 19:47

                               
登录/注册后可看大图

哎呀，刚刚看了你的数据，其中有几列完全都是-9999，这样的数据肯定是舍弃的。
那些好的数据应该就是你说的 ...

那几列全是-9999的确是舍弃不进行分析的。可是有些是直接整旬资料缺失，在资料里直接没有（连-9999也不显示）。那么我怎么进行统计呢？

mofangbao

看了一下你的数据，给你一个我的思路：

1、新建一个文件用于存放挑选出来的结果；
2、确定你第一个时间应该是2008年5月第一旬（假设）
3、循环读取文件，每次一行，分别提取出数据中的年、月、旬，如果这个时间晚于第二步中确定的时间，把这个时间之前的数据写入新的文件，数据值当然都是缺测，再把这行记录写在最后，如果刚好是同一时间，就直接写出，然后把第一个时间赋值为本行记录的下一个时次，回到循环开始继续读取下一条记录，再和第一个时间（已变为上一条记录）比较，如果是晚于这个时间，那么把这两个时间段内的数据写出，数据值为缺测，如果刚好是同一时间，就直接写出，如此循环，直到数据读取完毕，退出循环，关闭文件。
4、对新生成的这个文件，采用深深介绍的其中某种方法，对缺测值进行处理，因为现在你的数据是完整的了，接下也就方便了。

思路大概这样，实现就要靠你自己了

孤鹜

mofangbao 发表于 2012-4-19 12:20

                               
登录/注册后可看大图

看了一下你的数据，给你一个我的思路：

1、新建一个文件用于存放挑选出来的结果；

哦，我明白了~~~这应该能写出来，主要是没找到思路。谢谢~~

言深深

孤鹜 发表于 2012-4-19 11:22

                               
登录/注册后可看大图

那几列全是-9999的确是舍弃不进行分析的。可是有些是直接整旬资料缺失，在资料里直接没有（连-9999也不显 ...

额，程序判断一下，如果确实的话可以用历史旬去补

孤鹜

mofangbao 发表于 2012-4-19 12:20

                               
登录/注册后可看大图

看了一下你的数据，给你一个我的思路：

1、新建一个文件用于存放挑选出来的结果；

我写了两个小程序来分步输出，可能是读取、判断语句写得复杂了，程序运行了十几分钟才结束。不过，总算完成了TT

孤鹜

言深深 发表于 2012-4-19 16:11

                               
登录/注册后可看大图

额，程序判断一下，如果确实的话可以用历史旬去补

嗯，那个已经做好了。下面要进行平均了。谢谢~

言深深

孤鹜 发表于 2012-4-19 16:26

                               
登录/注册后可看大图

嗯，那个已经做好了。下面要进行平均了。谢谢~

nice， comeon