图幅理论面积与图斑椭球面积计算公式及要求

1图幅理论面积与图斑椭球面积计算公式及要求一、图幅理论面积计算公式4pb2L1PAsin(BB)cosB3Bsin(BB)cos3BCsin5(BB)cos5B36060221m221m221m79Dsin(BB)cos7BEsin(BB)cos9B221m221m(1)式中：椭球长半轴单位：米）a椭球扁率，椭球短半轴（单位：米）。e2=2-22。=+()e2+(/)8e04+(35/)1e162+(630/)2eo=()e2+(/)8e04+(21/)1e162+(420/)2eo=()e+(1)1e6+(180/)2eo=(1)1e26+(45/2)e。，=(5/2)3e4图幅东西图丿郭的经差（单位：分）。-图幅南北图廓的纬差（单位：弧度），=m(B+B)/)二、椭球面上任意梯形面积计算公式135cos5BmS2b皿卜幻(B2叫C0SBm珈阪他叫C0S3BmCSin2(B27DSin2(B2叫)C0S7Bm9Esi巧(B2B)C0S9Bm2）其中：为常数，按下式计算e2=(a2-b2)/a2A=1+(3/6)e2+(30/80)e4+(35/112)e6+(630/2304)e8B=(1/6)e2+(15/80)e4+(21/112)e6+(420/2304)e8C=(3/80)e4+(7/112)e6+(180/2304)e81/112)e6+(45/2304)e85/2304)e8式中：一椭球长半轴单位：米，一椭球短半轴单位：米；图块经差单位：弧度；-一图块纬差单位：弧度=（+）。三、高斯投影反解变换（x,yB丄）模型yy500000带号oooooo（若坐标不带带号，则不需减去带号x;EK0xBfEcosE(KsinE叫sin3EK3sin5E%sin7E)BB121f2720吨245t4!LbsB61osB击“282224t42g2.osB+中央子午线经度值（弧度）ttgB.2ellcosBNC/VCa2/bV.12/fK,K,K,K,K为与椭球常数有关的量。01234公式说明：若坐标为没有带号前缀格式，则不需减去带号X0若坐标为有带号前缀格式，则需减去带号X。四、计算用到的常数、椭球参数在计算图幅理论面积与任意图斑椭球面积时，有关常数及保留的位数按给定数值计算常数：n=椭球常数：椭球长半轴a椭球扁率,椭球短半轴b=6356755.29椭球第一偏心率e2=6.69438499958795E-03椭球第二偏心率e!2极点子午圈曲率半径c相关常数:五、计算中的取位及要求 高斯投影反解变换后的保留到秒后位小数，四舍五入。 采用计算机计算时，所有变量数据类型均要定义为双精度。 面积计算结果以平方米为单位，保留一位小数，四舍五入。 各种比例尺标准分幅图经差、纬差见表1。 在用大地坐标生成标准分幅图框时，要求在每条边框线的整秒处插入加密点。表1各种比例尺标准分幅图经差、纬差表比例尺万万万万万万万千经差OO/570351”纬差OOO0052015六、任意图斑椭球面积计算方法任意封闭图斑椭球面积计算的原理：将任意封闭图斑高斯平面坐标利用高斯投影反解变换模型，将高斯平面坐标换算为相应椭球的大地坐标，再利用椭球面上任意梯形图块面积计算模型计算其椭球面积，从而得到任意封闭图斑的椭球面积。1、计算方法：任意封闭区域总是可以分割成有限个任意小的梯形图块，因此，任意封闭区域的面积Ps，式中Si为分割的任意小的梯形图块面积（）用公式（2）计算。ii求封闭区域（多边形如图1）ABCD的面积，其具体方法为：(1)对封闭区域(多边形)的界址点连续编号(顺时针或逆时针)ABCD,提取各界址点的高斯平面坐标A(X1,Y1),B(X2,Y2),C(X3,Y3),D(X4,Y4)；(2) 利用高斯投影反解变换模型公式(3),将高斯平面坐标换算为相应椭球的大地坐标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,用公式(2)计算出各小梯形图块AEiFiA1的面积Si,然后累加Si就得到梯形图块ABB1A1的面积，同理，依次计算出梯形图块BCC1B1、CDD1C1、DAA1D1的面积(注：用公式(2)计算面积时，B1、B2分别取沿界址点编号方向的前一个、后一个界址点的大地纬度，DL为沿界址点编号方向的前一个、后一个界址点的大地经度的平均值与L0的差)；(5) 多边形ABCD的面积就等于4个梯形图块(ABB1A1、BCC1B1、CDD1C1、DAA1D1)面积的代数和。图1椭球面上任意多边形计算面积则任意多边形ABCD的面积P为：P=ABCD=BCC1B1+CDD1C1+DAA1D1-ABB1A12、计算要求 利用图形坐标点将高斯坐标系下的几何图形反算投影到大地坐标系，进行投影变换。 任意指定一条经线L0,从选定多边形几何形状的起始点开始，沿顺时针方向依次计算相邻两点构成的线段，以及两点到指定经线的平行线构成的梯形面积。 计算过程中应顺同一方向依坐标点逐个计算相邻两点连线与任意经线构成的梯形面积，坐标点不得有遗漏。若多边形包含内多边形（洞），则该多边形面积为外多边形面积减去所有内多边形面积之和。 计算所有梯形面积的代数和即为该多边形的面积。七、算法伪代码描述为了确保编程使用的参数、算法一致，保证不同软件计算的椭球面积一致，我们用算法伪代码描述的方法对编程进行统一，在利用计算机编制椭球面积计算软件时，计算参数与计算顺序应严格按照以下代码执行。1、概述计算规则：两个绝对值很大的数或两个绝对值很小的数相乘时，不能用幂计算符，并且尽量不使用连乘，应使用下述的方法：A*B*A*B；数据类型：当使用.net环境时，用Decimal数据类型代替Double类型；2、参数说明双精度类型：圆周率值：中央经线0椭球常数椭球长半轴：a椭球短半轴：b椭球扁率：椭球第一偏心率椭球第二偏心率极点子午圈曲率半径：3、算法描述初始化参数参数初始化结束中央经线转换为弧度选定本初子午线为参考经线起始点倒数第二点由高斯坐标反解计算经纬度值将经纬度转换为弧度值7计算梯形面积梯形面积值经差纬差纬度和计算变量按/照以下计算顺序：短半径*经差*短半径4、高斯坐标反解算法8秒/保/留到小数点后6位，四舍五入9