大地测量相关公式及参数docWord文件下载.docx
《大地测量相关公式及参数docWord文件下载.docx》由会员分享，可在线阅读，更多相关《大地测量相关公式及参数docWord文件下载.docx（12页珍藏版）》请在冰点文库上搜索。
△L—图幅东西图廓的经差（单位：
弧度）。
（B2-Bi）—图幅南北图廓的纬差（单位：
弧度）,Bm=（BmB2）/2o
二、椭球面上任意梯形面积计算公式
013
S=2/rALAsin—（B2—B|）cos^m-Bsin-（^2-
22
cos3B|]]+Csin：
（〃2一用）cos5Bm
79
-Dsin—（Z?
2-B|）cos7+Esin—（Z?
2-）cos9Z?
m
（2）
（1）
其中：
A,B,C,D,E为常数,按下式计算:
C2=（a2-b2）/a2
A=1+
（3/6）e2+（30/80）
e4+（35/112）e6+
（630/2304）
e8
（1/6）e2+（15/80）
e4+（21/112）e6+
（420/2304）
（3/80）
e4+（7/112）e6+
（180/2304）
（1/112）e6+
（45/2304）
（5/2304）
c8
米），b一椭球短半轴（单位：
米）；
△L—图块经差（单位：
弧度）；
（B2-BJ-图块纬差（单位：
弧度）
Bm=（Bi+B2）/2o
三、高斯投影反解变换（模型
矿=y-500000-带号x1000000（若坐标不带带号，则不需减去带号X1000000；
）
E=K*
Bf=E+cosE（K{sinE-K2sin3E+K3sin5E-K4sin7E）
+中央子午线经度值（孤度）（3）
t=tgBf7]2=e,2cos2BfN=C/VC=a2/hV+fKo，Ki，K>
K3，K4为与椭球常数有关的量。
公式说明：
若坐标为没有带号前缀格式，则不需减去带号X1000000；
若坐标为有带号前缀格式，则需减去带号XIOOOOOOo
四、计算用到的常数、椭球参数
在计算图幅理论面积与任意图斑椭球面积时，有关常数及保留的位数按给定数值计算。
常数：
JI=3.14159265358979
p=206264.8062471
80椭球常数：
轴a=6378140椭球扁率°
=1/298.257
椭球短半轴》=6356755.29
椭球第二偏心率=6.73950181947292E-03
极点子午圈川1率半径c=6399596.65198801
相关常数：
ko=1.57048687472752E-07
ki=5.05250559291393E-03
k2=2.98473350966158E-05
k3=2.41627215981336E-07
k4=2.22241909461273E-09
五、计算中的取位及要求
%1高斯投影反解变换后的B,L以秒为单位，保留到小数点后6位，四舍五入。
%1采用计算机计算时，所右变量数据类型均要定义为双精度。
%1面积计算结果以平方米为单位，保留一位小数，四舍五入。
%1各种比例尺标准分幅图经差、纬差见表1。
%1在用大地坐标生成标准分幅图框时，要求在每条边框线的整秒处插入加密点。
表1各种比例尺标准分幅图经差、纬差表
比例尺
1:
100万
50万
25万
10万
1：
5万
2.5万
1万
5千
经差
6°
3°
1。
30，
30'
15，
1!
30〃
3,45〃
52.5”
纬差
4°
2°
1°
20，
5，
2Z30”
15"
六、任意图斑椭球面积计算方法
任意封闭图斑椭球面积计算的原理:
将任意封闭图斑高斯平而坐标利用高斯投影反解变换模型，将高斯平而坐标换算为相应椭球的大地坐标，再利用椭球面上任意梯形图块面积计算模型计算其椭球血积，从而得到任意封闭图班的椭球血积。
1、计算方法：
任意封闭区域总是可以分割成有限个任意小的梯形图块，因此，任意封闭区域的而积
P=式中Si为分割的任意小的梯形图块面积（i=l,2,...n）用公式
（2）计算。
/=1
求封闭区域（多边形如图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）任意给定一经线LO（如L0=60。
），这样多边形ABCD的各边AB、BC、CD、DA与L（）就围成了4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）；
（4）由于在椭球面上同一经差随着纬度升高，梯形图块的面积逐渐减小，而同一纬差上经差梯形图块的面积相等，所以，将梯形图块ABB1A1按纬差分割成许多个小梯形图块AEiFiAl,用公式
（2）计算出各小梯形图块AEiFiAl的面积Si,然后累加Si就得到梯形图块ABB1A1的面积,同理,依次计算出梯形图块BCC1B1、CDDICkDAA1D1的面积（注:
用公式
（2）计算面积时，Bl、B2分别取沿界址点编号方向的前一个、后一个界址点的大地纬度，AL为沿界址点编号方向的前一个、后一个界址点的大地经度的平均值与L0白笼）；
（5）多边f?
ABCD的面积就等于4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）而积的代数和。
图1椭球面上任意多边形计算面积
则任意多边形ABCD的面积P为：
P=ABCD=BCC|B|+CDD|C）+DAA】D|.ABB）A）
2、计算要求
%1利用图形坐标点将高斯坐标系下的儿何图形反算投影到大地坐标系，进行投影变
换。
%1任意指定一条经线L0,从选定多边形儿何形状的起始点开始，沿顺时针方向依次计算相邻两点构成的线段，以及两点到指定经线的平行线构成的梯形而积。
将该梯形沿纬度变化方向（Y轴）进行切割，至少需切割为2个部分。
%1计算过程中应顺同一方向依坐标点逐个计算相邻两点连线与任意经线构成的梯形面积，坐标点不得有遗漏。
若多边形包含内多边形（洞），则该多边形面积为外多边形面积减去所有内多边形面积之和。
%1计算所有梯形面积的代数和即为该多边形的面积。
七、算法伪代码描述
为了确保编程使用的参数、算法一致，保证不同软件计算的椭球面积一致，我们用算法伪代码描述的方法对编程进行统一，在利用计算机编制椭球面积计算软件时，计算参数与计算顺序应严格按照以下代码执行。
1、参数说明
双精度类型：
圆周率值:
PI=3.14159265358979
中央经线:
CenterL
RHO=206264.8062471
A:
ParamA
B:
ParamB
C:
ParamC
D:
ParamD
E:
ParamE
ConstZEROAsDouble=0.000000000001
80椭球常数
椭球长半轴：
aRadius=6378140
椭球短半轴：
bRadius=6356755.29
椭球扁率：
ParaAF=1/298.257
椭球第一偏心率：
ParaEl=6.69438499958795E-03椭球第二偏心率:
ParaE2二6.73950181947292E-03极点子午圈曲率半径：
ParaC=6399596.65198801
k0:
ParakO
kl:
Parakl
=1.57048687472752E-07
=5.05250559291393E-03
k2:
Parak2=2.98473350966158E-05
k3:
Parak3=2.41627215981336E-07
k4:
Parak4=2.22241909461273E-09
2、算法描述
初始化参数
Doublee;
Doublea;
e=ParaE2;
ParaC=aRadius/（1-ParaAF）;
ParamA=1+（3/6）*o+（30/80）*Power（o,2）+（35/112）*Power（e,3）+（630/2304）*Power（e,4）;
ParamB=（1/6）*e+（15/80）*Power（e,2）+（21/112）*Power（e,3）+（420/2304）*Power（c,4）;
ParamC=（3/80）*Power（e,2）+（7/112）*Power（e,3）+（180/2304）*Power（e,4）;
ParamD=（1/112）*Power（e,3）+（45/2304）*Power（e,4）;
ParamE=（5/2304）*Power（e,4）;
参数初始化结束
中央经线转换为弧度
ContcrL=TransDegrecToArc（CcnterL）
选定本初子午线为参考经线
StandardLat=0
For起始点To倒数第二点
由高斯坐标反解计算经纬度值
ComputeXYGeo（PntColl.Point（i）.y,PntColl.Point（i）.x,B,L,CenterL）
ComputeXYGeo（PntColl.Point（i+1）.y,PntColl.Point（i+1）.x,Bl,LI,CenterL）
将经纬度转换为弧度值
B=B/RII0
L=L/RHO
Bl=Bl/RHO
LI=LI/RHO
计算梯形面积
DoubleAreaVal;
//梯形面积值
DoubleIDiference;
〃经差
DoublebDiference;
〃纬差
DoublebSum;
//纬度和
DoubleItcmValue（5）；
//计算变量
bDiference=（Bl-BO）；
bSum=（Bl+BO）/2；
IDiference=（LI+L）/2;
ParamA*Sin（bDiference/2）*Cos（bSum）:
=ParamI3*Sin（3*bDiference/2）*Cos（3*bSum）;
AreaVal=2*bRadius*IDiference*bRadius*（TtemValue（0）-TtemValuc（l）+ItemValue
（2）-ItcmValuc（3）+ItemValue（4））;
areaSum二areaSum+AreaVal;
Next
EndSub
3、高斯坐标反解算法
PublicSubComputeXYGco（xAsDouble,yAsDouble,BAsDouble,LAsDouble,centerAsDouble）
DimylAsDoubleDimbfAsDoubleyl=y-500000
DimeAsDouble
e=ParakO*x
DimseAsDouble
se=Sin（e）
bf=e+Cos（e）*（Parakl*se-Parak2*Power（se,3）+Parak3*Power（se,5）-Parak4*Power（sc,7））
DimvAsDouble
DimtAsDouble
DimNAsDouble
DimnlAsDouble
DimvtAsDouble
DimynAsDouble
Dimt2AsDouble
DimgAsDouble
t=Tan（bf）
nl=ParaEl*Power（Cos（bf）,2）
v=Sqr（1+nl）
N=ParaC/v
yn=yl/N
vt=Power（v,2）*t
t2=Power（t,2）
B=bf-vt*Power（yn,2）/2+（5+3*t2+nl-9*nl*t2）*vt*Power（yn,4）/24
-（61+90*t2+45*Power（t2,2））*vt*Power（yn,6）/720
B=TransArcToDegree（B）
DimcbfAsDouble
cbf=1/Cos（bf）
L=cbf*yn-（1+2*t2+nl）*cbf*Power（yn,3）/6+（5+28*t2+24*Power（t2,
2）+6*nl+8*nl*12）*cbf*Power（yn,5）/120+center
L=TransArcToDegree（L）
弧度转换为度
PublicFunctionTransArcToDegree（arcAsDouble）AsDouble
DimdegreeAsDouble
DimminAsDouble
DimseeAsDouble
DimretAsDouble
DimtmpAsDouble
ret=arc♦180/PT
degree=FormatValue（ret,100,100）tmp=（ret-degree）*60
min二FormtitValue（tmp,100,100）sec=（tmp-min）*60
//秒保留到小数点后6位，四舍五入
sec=Format（sec,〃####.000000〃）1EormatValue（sec,10000000,100）
TransArcToDegree=degree*3600+min*60+sec
EndFunction
PrivateFunctionFormatValuc（inputValAsDouble,procsionAsLong,scalcNumAsLong）AsDouble
FormtitValue二（Int（inputVcil*precsion）-Int（inputVal*precsion）ModsealeNum）/precsion
升级会员 