[分享] 椭球面积计算公式

admin

1. double ParaE1 = m采用de11;//WGS2000采用E1;//0.00673949677548;  //WGS-2000椭球偏心率e的平方
2.         double E1 = ParaE1;
3.         // ParamA—ParamE
4.         double ParamA = 1.0 + 0.5 * E1 + 0.375 * E1 * E1 + 0.3125 * E1 * E1 * E1 + 0.2734375 * E1 * E1 * E1 * E1;
5.         double ParamB = 0.1666666666666666666666666667 * E1 + 0.1875 * E1 * E1 + 0.1875 * E1 * E1 * E1 + 0.1822916666666666666666666667 * E1 * E1 * E1 * E1;
6.         double ParamC = 0.0375 * E1 * E1 + 0.0625 * E1 * E1 * E1 + 0.078125 * E1 * E1 * E1 * E1;
7.         double ParamD = 0.0089285714285714285714285714 * E1 * E1 * E1 + 0.01953125 * E1 * E1 * E1 * E1;
8.         double ParamE = 0.0021701388888888888888888889 * E1 * E1 * E1 * E1;
11.         // 两点围成梯形面积
12.         double area=0;
13.         double lDiference = 0.0;    // 经差
14.         double bDiference = 0.0;    // 纬差
15.         double bSum = 0.0;          // 纬度和
17.         // 纬差||纬度和||经差
18.         bDiference = (B1 - B)/ 2.0;
19.         bSum = (B1 + B) / 2.0;
20.         lDiference = (L1 + L) / 2.0;
23.         double ItemValue[5];
24.         double bRadius=m采用db;//WGS2000采用eb;
25.         // 按照以下计算顺序：短半径*经差*短半径
26.         double num1 = 2.0 * bRadius * lDiference * bRadius;
27.         double num2 = cos(bSum);
28.         double num3 = sin(bDiference);
29.         double d1,d2,d3,d4,d5;
30.         d1=ItemValue[0] = num1 * ParamA * num2 * num3;
31.         d2=ItemValue[1] = num1 * ParamB * sin(3.0 * bDiference) * cos(3.0 * bSum);
32.         d3=ItemValue[2] = num1 * ParamC * sin(5.0 * bDiference) * cos(5.0 * bSum);
33.         d4=ItemValue[3] = num1 * ParamD * sin(7.0 * bDiference) * cos(7.0 * bSum);
34.         d5=ItemValue[4] = num1 * ParamE * sin(9.0 * bDiference) * cos(9.0 * bSum);
35.         area = ItemValue[0] - ItemValue[1] + ItemValue[2] - ItemValue[3] + ItemValue[4];

复制代码