sinLcosL????\?\0???X?(N?H)cosB??L??(N?H)cosB?????B???sinBcosLsinBsinLcosB???Y???????\??\?\?(M?H)(M?H)(M?H)??????H??????Z?sinBsinLsinB???cosBcosL??N(1?e2)?HtgBcosL?N?H??(N?H)?Ne2sin2BsinL??M?H?2???NesinBcosBsinL0???N?2???esinBcosB?\?m?M??(N?H)?Ne2sin2B????N(1?e2)?HtgBsinL?1???x?N?H???(N?H)?Ne2sin2B???y?cosL0??M?H??z???02NesinBcosBcosL??00????22(2?esinB)N?sinBcosB?\???a?e2sinBcosB?\??Ma???1?f?N???f?M??(1?e2sin2B)(1?e2sin2B)sin2B???1?a?a?
?B,?L,?H 同一点位在两个坐标系下的纬度差、经度差、大地
高差，经纬度差单位为弧度，大地高差单位为米，
??180?3600/? 弧度秒
?a
椭球长半轴差，单位为米，
?f 扁率差，无量纲， ?X,?Y,?Z 平移参数，单位为米，
?x,?y,?z 旋转参数，单位为弧度，
m 尺度参数，无量纲。
（二）高斯正反算公式 1、高斯投影正算公式
l2?1l21l4?2242244x???NtcosB2?0.5?(5?t?9??4?)cosB2?(61?58t?t)cosB4?24720?????2－32－
l?1l21l4?222242224y?NcosB?1?(1?t??)cosB2?(5?18t?t?14??58?t)NcosB4???6120???子午线弧长X计算见附录。
2、高斯投影反算公式
?1y21y4?222224)?1?(5?3tf??f?9?ftf)()?(61?90tf?45tf)()?12N360Nf??f???y?1y21y4?2224222l?()?1?(1?2tf??f)()?(5?28tf?24tf?6?f?8?ftf)()?cosBfNf?6N120N?ff??BfBf?ftfyB?Bf?y(2MfNf?tf式中、分别为按值计算的相应量，的计算见附录。
（三）坐标转换精度评定和评估方法
对于1954年北京坐标系、1980西安坐标系与2000国家大地坐标系转换分区转换及数据库转换点位的平均精度应小于图上的0.1mm。具体：
对于1:5千坐标转换，1980西安坐标系与2000国家大地坐标系转换分区转换平均精度≤0.5m；1954年北京坐标系与2000国家大地坐标系转换分区转换平均精度≤1.0m；
1:5万基础地理信息数据库坐标转换精度≤5.0m； 1:1万基础地理信息数据库坐标转换精度≤1.0m； 1:5千基础地理信息数据库坐标转换精度≤0.5m。 依据计算坐标转换模型参数的重合点的残差中误差评估坐标转换精度。对于n个点，坐标转换精度估计公式如下：
①V（残差）=重合点转换坐标-重合点已知坐标 ②空间直角坐标X残差中误差
MX??[vv]Xn?1 －33－
③空间直角坐标Y残差中误差
MY??[vv]Yn?1 ④空间直角坐标Z 残差中误差点位中误差Mp222?MX?MY?MZMZ??[vv]Zn?1
[vv]xn?1⑤平面坐标x残差中误差⑥平面坐标y残差中误差⑦大地高H残差中误差
平面点位中误差为MpMx?? My??[vv]yn?1MH??[vv]Hn?1 22?Mx?My
－34－
附录：
1、常用量定义
a为椭球长半轴，1954年北京坐标系为6378245m
1980西安坐标系为6378140m
b为椭球短半轴
f为椭球扁率, 1954年北京坐标系为1/298.3
1980西安坐标系为1/298.257
f?a?ba
b?a1?e2
e─第一偏心率
2e?a?b2 e2a?2f?f2 e? ─第二偏心率 2e'?a?b2b ?2?e'2cos2B
t?tgB
V?1?e'2cos2B V2?1??2 W?1?e2sin2B
B为纬度，单位弧度
c?a2b
35－－
a(1?e2)c?M─子午圈曲率半径 M? W3V3N─卯酉圈曲率半径 ??ac? WV2、子午线弧长X
设有子午线上两点p1和p2，p1在赤道上p2的纬度为B，p1、p2间的子午线弧长X计算公式
X?a(1?e2)(A?arcB?B?sin2B?C?sin4B?D?sin6B
?E?sin8B?F?sin10B?G?sin12B)式中
3451756110258436591069369312A??1?e2?e4?e?e?e?e4642561638465536104857632154525622058727651029729712B??e?e?e?e?e?e832102440961310725242881541056220581039510148648512C??e?e?e?e?e256102416384655368388608356105810395105505512D??e?e?e?e
30724096262144104857631583465109909912E??e?e?e131072524288838860869310900912F??e?e13107205242880100112G??e83886083、底点纬度Bf迭代公式
B0?X?F(Bi)X ,B?B?i?1iF?(Bi)a(1?e2)A直到Bi?1?Bi小于某一个指定数值，即可停止迭代。 式中
－36－
F(B)?a(1?e2)[A?arcB?B?sin2B?C?sin4B?D?sin6B
?E?sin8B?F?sin10B?G?sin12B]F?(B)?a(1?e2)[A??2B?cos2B?4C?cos4B?6D?cos6B
????8Ecos8B?10Fcos10B?12Gcos12B]
－37－