arcgis 椭面积计算公式.docVIP

图幅理论面积与图斑椭球面积计算公式及要求 一、 图幅理论面积计算公式 （1） 式中： a—椭球长半轴 单位：米 ，α—椭球扁率，b—椭球短半轴 单位：米 。 е2﹦ a2﹣b2 /a2。 A﹦1﹢（3/6）е2﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8。 B﹦ （1/6）е2﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8。 C﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8。 D﹦ （1/112）е6﹢ （45/2304）е8。 E﹦ （5/2304）е8。 ΔL—图幅东西图廓的经差（单位：分）。 B2﹣B1 —图幅南北图廓的纬差（单位：弧度），Bm﹦（B1﹢B2）/2。 二、椭球面上任意梯形面积计算公式 （2） 其中：A,B,C,D,E 为常数，按下式计算: е2﹦ a2﹣b2 /a2 A﹦1﹢（3/6）е2﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8 B﹦ （1/6）е2﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8 C﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8 D﹦ （1/112）е6﹢（45/2304）е8 E﹦ （5/2304）е8 式中：a—椭球长半轴 单位：米 ，b—椭球短半轴 单位：米 ； ΔL—图块经差 单位：弧度 ； B2﹣B1 —图块纬差 单位：弧度 Bm﹦（B1﹢B2）/2。 三、高斯投影反解变换（）模型 （若坐标不带带号，则不需减去带号×1000000；） +中央子午线经度值（弧度） （3） 式中： 公式说明：若坐标为没有带号前缀格式，则不需减去带号×1000000；若坐标为有带号前缀格式，则需减去带号×1000000。 四、计算用到的常数、椭球参数 在计算图幅理论面积与任意图斑椭球面积时，有关常数及保留的位数按给定数值计算。 常数： π﹦3.1座机电话号码58979 206264.座机电话号码 80椭球常数： 座机电话号码 1/ 298.257 座机电话号码.29 6.座机电话号码958795E-03 6.座机电话号码947292E-03 座机电话号码.座机电话号码 相关常数: k0 1.座机电话号码472752E-07 k1 5.0座机电话号码91393E-03 k2 2.座机电话号码966158E-05 k3 2.座机电话号码981336E-07 k4 2.座机电话号码461273E-09 五、计算中的取位及要求 ①高斯投影反解变换后的B,L保留到后位，四舍五入② 采用计算机计算时，所有变量数据类型均要定义为双精度。 ③ 面积计算结果以平方米为单位，保留一位小数，四舍五入 六、任意图斑椭球面积计算任意封闭图斑椭球面积计算的原理：将任意封闭图斑高斯平面坐标利用高斯投影反解变换模型，将高斯平面坐标换算为相应椭球的大地坐标，再利用椭球面上任梯形图块面积计算模型计算其椭球面积，从而得到任意封闭图斑的椭球面积。 任封闭区域总是可以分割成有限个任意小的梯形图块，因此，任封闭区域的面积，式中Si为分割的任意小的梯形图块面积（i 1,2,…n）用公式 计算。 求封闭区域（多边形如图1）ABCD的面积 ，其具体方法为： （1）对封闭区域（多边形）的界址点连续编号（顺时针或逆时针）ABCD，提取各界点的高斯平面坐标A X1,Y1 ，B X2,Y2 ，C X3,Y3 ，D X4,Y4 ； （2）利用高斯投影反解变换模型公式（），将高斯平面坐标换算为相应椭球的大地坐标A B1,L1 ，B B2,L2 ，C B3,L3 ，D B4,L4 ； （3）任意给定一经线L0（如L0＝60°），这样多边形ABCD的各边AB、BC、CD、DA与L0就围成了4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）； （4）由于在椭球面上同一经差随着纬度升高，梯形图块的面积逐渐减小，而同一纬差上经差梯形图块的面积相等，所以，将梯形图块ABB1A1按纬差分割成许多个小梯形图块AEiFiA1，用公式（）计算出各小梯形图块AEiFiA1的面积Si，然后累加Si就得到梯形图块ABB1A1的面积，同理，依次计算出梯形图块BCC1B1、CDD1C1、DAA1D1的面积（注：用公式（）计算面积时，B1、B2分别取沿界址点编号方向的前一个、后一个界址点的大地纬度，ΔL为沿界址点编号方向的前一个、后一个界址点的大地经度的平均值与L0的差）； （5）多边形ABCD的面积就等于4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）面积的代数和。 图1 椭球面上任意多边形计算面积 则任意多边形ABCD的面积P为： P ABCD BCC1B1