python求正多边形面积,不规则多边形面积公式 python

　　我用爪哇岛重写了矩阵实验室的" areaint "函数，得到了完全相同的结果。

　　" areaint "计算"每单位表面",所以我把答案乘以地球的表面积(5.10072e14平方米)。私有双精度区（数组列表lats，数组列表lons)

　　{

　　double sum=0；

　　double prev colat=0；

　　double prevaz=0；

　　双colat0=0

　　双az0=0；

　　for(int I=0；我

　　{

　　double colat=2 * math。atan2(数学。sqrt(数学。pow(数学。罪(拉特。得到(I)yxdhy/180/2)，2)数学。cos(拉特。get(I)yxdhy/180)* math。pow(数学。罪(lons。get(I)yxdhy/180/2)、2))、math。sqrt(1-数学。pow(数学。罪(拉特。get(I)yxdhy/180/2)，2)-数学

　　double az=0；

　　if (lats.get(i)=90)

　　{

　　az=0；

　　}

　　else if (lats.get(i)=-90)

　　{

　　az=数学；

　　}

　　其他

　　{

　　az=数学。atan2(数学。cos(拉特。get(I)yxdhy/180)* math。罪(lons。get(I)yxdhy/180)，数学。罪(拉特。get(I)yxdhy/180))%(2 yxdhy)；

　　}

　　如果(i==0)

　　{

　　colat0=colat

　　az0=az

　　}

　　如果(i0我

　　{

　　sum=sum(1-数学。cos(prev colat(colat-prev colat)/2))yxdhy *(数学。ABS(az-prevaz)/数学.PI)-2 *数学。ceil((数学。ABS(az-prevaz)/数学.PI)-1)/2))*数学。signum(az-prevaz)；

　　}

　　prevcolat=colat

　　}

　　sum=sum(1-数学。cos(prev colat(colat 0-prev colat)/2))*(az0-prevaz)；

　　返回5.10072 e14 *数学。最小(数学。ABS(sum)/4/Math .PI，1-Math.abs(sum)/4/Math .PI)；

　　}

郑重声明：本文由网友发布，不代表盛行IT的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。