求助，利用fortran的数据选取问题

大鹏

34.718         118.83297 -5.156
34.91903 118.62345 -13.028
34.89554 118.63566 16.263
34.86106 118.40833 -4.479
34.90326 118.59366 26.529
34.88237 118.59694 6.928
34.72429 118.62522 -4.469
34.7278   118.63825 -11.215
35.17694 118.05029 -3.813
34.61176 118.61868 -5.979
34.79099 118.45542 -2.865
34.69324 118.93502 -4.036
34.50609 118.80571 -6.385
34.64913 118.64637 -2.889
34.15522 119.05854 -11.394
34.514         118.69115 -7.257
33.70742 120.38807 -19.24
34.45707 118.88397 -12.443
34.46544 118.8114 -14.403
34.18099 118.91912 -7.991
31.5869         116.75586 -5.374
31.57114 117.90841 -34.245
31.65319 117.26082 21.861
35.41562 117.58068 -6.23
37.39619 120.68105 -14.844
35.83979 117.7482 -5.092
37.42891 121.07681 -8.816
37.29512 121.17327 22.912
37.41264 120.75953 -7.938
37.34418 120.80621 -15.01
37.3038         120.79034 -9.26
36.99829 120.81897 -7.3
40.26443 119.96423 -30.914
36.14707 116.56629 -19.611
37.25105 122.02326 -12.51
37.08149 120.94093 19.794
37.11503 120.98639 21.823
36.99805 120.93431 18.176
37.24081 121.00047 19.334
37.26819 121.01183 19.826
37.24779 121.09838 25.224
37.22193 121.52854 -21.998
37.01967 120.91542 20.888
第一列是纬度，第二例是经度第三列是强度数据，这是09年.07.01当天山东全省的闪电分布的数据，因为山东有十七个地市，如何利用经纬度来选取每一个地市当天的数据，这是第一个需要求助的程序（如果直接利用各个市的经纬度作为帅选条件的话，这好像就成了一个长方形的网格了，与实际的曲线边界有些不同）；再个如果想得到以一个经纬度点为圆心，5或10km为半径的圆形区域内的频数，应该怎样编程，真心的求助啊 小论文急用 谢谢大家啊

大鹏

{:e311_b:}

topmad

第一个就是循环 然后判断筛选吧  不适合明白你的意思  试一下

第二个可以根据经纬度计算每个站点到其他所有站点的距离，然后筛选出符合条件的  输出就可以了

言深深

顶个板凳
注：x(i,j,1)是经度x(i,j,2)是纬度,x(i,j,3)，ab分别是定点经纬度,r是半径
do i=1,m
do j=1,n
if((x(i,j,1)-a)**2+(x(i,j,2)-b)**2.le.r) then
x(i,j,3)就是在半径r圆心ab内部，可以自己运算
endif
enddo
enddo

大鹏

topmad 发表于 2012-10-16 17:45

                               
登录/注册后可看大图

第一个就是循环 然后判断筛选吧  不适合明白你的意思  试一下

第二个可以根据经纬度计算每个站点到其他所 ...

第二个有点不懂 ，能不能给个例子啊

大鹏

topmad 发表于 2012-10-16 17:45

                               
登录/注册后可看大图

第一个就是循环 然后判断筛选吧  不适合明白你的意思  试一下

第二个可以根据经纬度计算每个站点到其他所 ...

请问您知道到哪里去获取各地市的经纬度的范围吗 在百度上没有这个范围

大鹏

言深深 发表于 2012-10-16 18:40

                               
登录/注册后可看大图

顶个板凳
注：x(i,j,1)是经度x(i,j,2)是纬度,x(i,j,3)，ab分别是定点经纬度,r是半径
do i=1,m

if((x(i,j,1)-a)**2+(x(i,j,2)-b)**2.le.r)在这个判定条件中，前面是经纬度的平方和 后面r是距离两个量可以比较吗是不是需要把距离r换算成经纬度的间隔呢

topmad

大鹏 发表于 2012-10-16 19:34

                               
登录/注册后可看大图

if((x(i,j,1)-a)**2+(x(i,j,2)-b)**2.le.r)在这个判定条件中，前面是经纬度的平方和 后面r是距离两个量可 ...

看看地球科学概论之类的书，看看球面距离求算，要跟地球半径结合起来的   
你前面算到的是平方和，后面的是距离  不能比较，要开方

言深深

topmad 发表于 2012-10-16 20:01

                               
登录/注册后可看大图

看看地球科学概论之类的书，看看球面距离求算，要跟地球半径结合起来的   
你前面算到的是平方和，后面的 ...

好吧，这边严格来说是不对的

言深深

大鹏 发表于 2012-10-16 19:34

                               
登录/注册后可看大图

if((x(i,j,1)-a)**2+(x(i,j,2)-b)**2.le.r)在这个判定条件中，前面是经纬度的平方和 后面r是距离两个量可 ...

只是做一个粗略的意思，需要进行转换的，有公式的