求FORTRAN程序——在地图上，求一条等值线内的面积

yuxiao19910

           最初的是一个全球的处理过的数据a，想求：
           在90E-90W,0-90N 范围内，a>=3200的面积，也就是a=3200这条等值线内部围成的总面积。
           数据是2.5*2.5的，应该可以把它当成求若干个方格的面积。
           但是，纬度加权肯定是必须的。

           求给力啊！！！~~~~(>_<)~~~~

Aires

我也想要 呵呵

yuxiao19910

Aires 发表于 2013-3-31 13:40

                               
登录/注册后可看大图

我也想要 呵呵

呵呵呵，哪位大神赶紧出现啊·······

言深深

程序的编写不是特别的难，但是有一个问题就是设计上面是不是有缺陷，比如存在某一个小区域，只有有限的1个或者2个点超过3200，那么是否需要呗识别出来：

整体程序可以用下的识别方案：
定义一个纬向数组，存放每一个纬度上面的个点数目，然后和全球的格点作对比，并上纬向面积权重就ok了
假定数据构成x(m,n),m=73,n=144分别是纬向和经向格点数目，主要s是动态的···

do i=1,m
s=0
do j=1,n
if(x(i,j).ge.3200) then
s=s+1
endif
enddo
y(i)=s
enddo

下面进入面积计算
s=0
do i=2,72
s=s+y(i)/n*pi*r/36*2*pi*r*cos(90-i*2.5)
enddo
极点的就算，面积只考虑一般
s=s+y(1)/n*pi*r/36*2*pi*r*cos(90-1*2.5)*0.5
s=s+y(73)/n*pi*r/36*2*pi*r*cos(90-73*2.5)*0.5

print*,"计算面积是：",s
end

刚刚拿笔算了一下，差不多这样设计应是ok的。思想就是计算每一个纬圈5°范围内的面积，计算的时候有点儿误差，误差的产生在于用的是对应纬圈长度乘以5°的长度······个人感觉是可以忽略的，当然楼主可以改良一下。

yuxiao19910

言深深 发表于 2013-3-31 14:59

                               
登录/注册后可看大图

程序的编写不是特别的难，但是有一个问题就是设计上面是不是有缺陷，比如存在某一个小区域，只有有限的1个或 ...

恩恩，差不多了然了，但是有两个没看懂：
s=s+y(i)/n*pi*r/36*2*pi*r*cos(90-i*2.5)
pi*r/36 这一个是什么意思啊？
90-i*2.5 这个不用换成弧度么？

yuxiao19910

言深深 发表于 2013-3-31 14:59

                               
登录/注册后可看大图

程序的编写不是特别的难，但是有一个问题就是设计上面是不是有缺陷，比如存在某一个小区域，只有有限的1个或 ...

哦~~~ 第一个 我懂了，是2*pi*r/72 的意思，小方块的高，差不多2.5度是吧，最后有一个 *0.5 是考虑到方块向极地的收缩么？

言深深

yuxiao19910 发表于 2013-3-31 15:36

                               
登录/注册后可看大图

哦~~~ 第一个 我懂了，是2*pi*r/72 的意思，小方块的高，差不多2.5度是吧，最后有一个 *0.5 是考虑到方块 ...

是的，你的理解是正确的哈！然后，那边确实应该换算成弧度的···因为直接在回复里面敲的，所以想的漏了···惭愧···本来真的担心你看不到，现在看来你的基础还是很不错的哈···nice···

yuxiao19910

言深深 发表于 2013-4-1 10:23

                               
登录/注册后可看大图

是的，你的理解是正确的哈！然后，那边确实应该换算成弧度的···因为直接在回复里面敲的，所以想的漏了 ...

我聪明嘛
  我弄出来了，thanks~~~