利用fortran选取样本个例

大鹏

我的数据如（只是一部分来举个例子）下：
年        月  日 时分秒毫秒   维度            经度       电流幅值
2009 06 05 08 31 37 398 36.29283 120.35725   16.837
2009 06 05 09 45 15 442 35.97290 120.62507  -17.169
2009 06 05 09 49 15 084 35.91164 120.61160  -23.925
2009 06 05 09 50 24 234 35.90533 120.66304  -14.064
2009 06 05 09 51 45 562 35.91808 120.67847  -23.675
2009 06 05 09 51 46 012 35.98477 120.63826  -12.856
2009 06 05 09 52 45 014 35.91652 120.67787  -17.328
2009 06 05 09 53 34 797 36.01581 120.64723  -21.914
2009 06 05 09 56 24 814 35.88767 120.62961  -12.617
2009 06 05 10 00 57 545 35.83185 120.66916   -8.416
2009 06 05 10 12 46 961 35.89435 120.67093   15.645
2009 06 05 10 38 36 600 36.13832 120.02913   -9.030
2009 06 05 10 44 32 532 36.05127 120.19465   -8.982
2009 06 05 10 46 07 680 36.11611 120.22468  -30.509
2009 06 05 10 46 33 485 36.05568 120.17908  -17.203
2009 06 05 10 46 57 936 36.07071 120.17706  -20.544
2009 06 05 10 47 15 978 36.08936 120.25863  -19.719
2009 06 05 10 47 43 451 36.06186 120.16666  -30.630
2009 06 05 10 47 51 453 35.98332 120.04471  -12.008
2009 06 05 10 48 18 213 36.06453 120.17844  -24.195
2009 06 05 10 48 24 491 35.98157 120.04153  -13.267
2009 06 05 10 48 40 152 36.05930 120.15053  -16.545
2009 06 05 10 48 54 417 36.09039 120.22959  -16.996
2009 06 05 10 49 45 823 36.05396 120.16038  -12.895
2009 06 05 10 49 58 422 36.08897 120.08882  -18.794
2009 06 05 10 50 18 868 36.02868 120.12343  -16.988

现在想实现的功能是：这是一个闪电定位数据，要选取雷暴个例，条件是时间间隔小于15分钟的所有地闪都归为同一个雷暴，然后继续选取，能不能把各个样本分别整到单独的文档里显示出来 ，谢谢各位 帮帮忙，关键是怎么把每一个样本整出来

大鹏

大伙要给力啊

言深深

方法就是程序识别，找15分钟内的，但是你给的这一条件程序没有办法识别，存在bug；
原因在于如下符合你需要的方案：
   1，选取定点Z，然后观察+15分钟以内的数据量，可能得到m个；
   2，定点Z+1，再观察+15分钟以内的数据量，得到n个，其中k个重复与定点Z的选取结果，这是显然的，这时候两次选取结果总数目为m+n-k；
   3，按以上1、2方法Z+L，L遍列数据总量，这样可以得到极少的几次过程。

方案中，得到的结果一次过程可能大于15分钟，第一次和第二次之间的间隔可能跨度达到30分钟甚至更多，这时候如果仍然需要满足你的15分钟以内的发生数据量，需要添加条件！

一种可能的解决方案是：

对以上123选取出来的几次较大的过程，绘制概率分布图，选择概率发生最大的15分钟，作为一次过程。

done

不知道以上的分析过程你是否明了，我没有做过雷暴，仅仅是对于你的需要，做了一下简单的分析，希望对你有帮助！
希望你至少花 30分钟以上的时间 读懂上面的解决方案！

mofangbao

是你自己要给力

mofangbao

言深深 发表于 2012-12-6 12:30

                               
登录/注册后可看大图

方法就是程序识别，找15分钟内的，但是你给的这一条件程序没有办法识别，存在bug；
原因在于如下符合你需要 ...

深深说的我还没有去仔细看，不过我相信你应该好好看看，因为他至少花了30分钟来分析你的问题了。
先说一下，你这样提问是不合适的，至少你自己要先思考，给出你自己的解决方案，把问题描述的更清楚一些。
看了你的提问，我也有几个问题比较疑惑：
1、间隔小于15分钟时针对当前样本的第一个时间还是当前样本的最后一条记录的时间？
2、是不是需要设置一个经纬度的范围，该范围内的闪电才能作为同一次过程？或者你的记录里面一定是一个小区域的范围？
对于第一个问题，如果是针对该样本第一个时间，那程序就比较好办，直接读一条，认为是一个样本，然后继续一条条读取，与当前样本比较即可，超过15分钟的作为下一个样本。如果是第二种情况，就比较复杂，因为会存在交叉的情况，可能真的只有像深深说的那样进行概率上的选取了，比如和上一条的间隔是10分钟，和下一条的间隔是6分钟，那么上下两条记录的间隔就是16分钟，不能作为一次过程，但是该记录却即可以作为上一过程，又可以作为下一过程。所以一定要给一些额外的限定条件来处理这种情况。
对于第二个问题，如果这个文件中是一个小区域的话，应该就没什么问题，要不然北京和南京的15分钟内的记录是一定不能作为同一个过程的。
楼主还是自己先好好思考下吧！

言深深

mofangbao 发表于 2012-12-6 14:19

                               
登录/注册后可看大图

深深说的我还没有去仔细看，不过我相信你应该好好看看，因为他至少花了30分钟来分析你的问题了。
先说一 ...

就是这样的，主要是条件给的有问题，想法没错，但是需要给程序识别的空间。

关于你说的第2点空间应该还好，因为他的研究范围应该不会太大，估计楼主应该是地方性的，数据也是地方的测站数据。否则要是拿全球/全国的数据来做是必然需要在写程序的时候对空间予以考虑的。
祝你成功

大鹏

言深深 发表于 2012-12-6 14:23

                               
登录/注册后可看大图

就是这样的，主要是条件给的有问题，想法没错，但是需要给程序识别的空间。

关于你说的第 ...

首先要非常感谢您，然后我的条件是给的比较模糊，上面mofangbao分析的很对，我的研究范围是在一个探空站附近几十公里内的，再个我选取的15分钟间隔是对于样本的最后一条数据而言的；因为我对fortran不是很熟，还想再问问，程序识别可以把每一个样本单独存放一个文本文档吗 您可不可以给个类似的程序让我学习一下，谢谢

大鹏

mofangbao 发表于 2012-12-6 14:19

                               
登录/注册后可看大图

深深说的我还没有去仔细看，不过我相信你应该好好看看，因为他至少花了30分钟来分析你的问题了。
先说一 ...

谢谢mofangbao，您的分析条条切中要害，真心佩服

言深深

大鹏 发表于 2012-12-6 20:37

                               
登录/注册后可看大图

首先要非常感谢您，然后我的条件是给的比较模糊，上面mofangbao分析的很对，我的研究范围是在一个探空站附 ...

很可惜你没有看懂我说的！等你看懂的时候再问吧······