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

图幅理论面积与图斑椭球面积计算公式及要求 一、 图幅理论面积计算公式 （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.14159265358979 206264.8062471 80椭球常数： = 6378140 = 1/ 298.257 = 6356755.29 = 6.69438499958795E-03 = 6.73950181947292E-03 = 6399596相关常数: k0 = 1.57048687472752E-07 k1 = 5.05250559291393E-03 k2 = 2.98473350966158E-05 k3 = 2.41627215981336E-07 k4 = 2.22241909461273E-09 五、计算中的取位及要求 = 1 \* GB3 ① 高斯投影反解变换后的B,L以秒为单位，保留到小数点后6位，四舍五入。 = 2 \* GB3 ② 采用计算机计算时，所有变量数据类型均要定义为双精度。 = 3 \* GB3 ③ 面积计算结果以平方米为单位，保留一位小数，四舍五入。 ④ 各种比例尺标准分幅图经差、纬差见表1。 ⑤ 在用大地坐标生成标准分幅图框时，要求在每条边框线的整秒处插入加密点。 表1 各种比例尺标准分幅图经差、纬差表 比例尺 1:100万 1:50万 1:25万 1:10万 1:5万 1:2.5万 1:1万 1:5千 经差 6o 3o 1o30′ 30′ 15′ 7′30″ 3′45″ 1′52.5″ 纬差 4o 2o 1o 20′ 10′ 5′ 2′30″ 1′15″ 六、任意图斑椭球面积计算方法 任意封闭图斑椭球面积计算的原理：将任意封闭图斑高斯平面坐标利用高斯投影反解变换模型，将高斯平面坐标换算为相应椭球的大地坐标，再利用椭球面上任意梯形图块面积计算模型计算其椭球面积，从而得到任意封闭图斑的椭球面积。 1、计算方法： 任意封闭区域总是可以分割成有限个任意小的梯形图块，因此，任意封闭区域的面积，式中Si为分割的任意小的梯形图块面积（i=1,2,…n）用公式(2)计算。 求封闭区域（多边形如图1）ABCD的面积