用EXCEL处理高斯坐标计算问题.docVIP

图 1 EXCEL 工作簿中分别建立正算和反算两个工作表Fig.1 Separate establishment of positive and reverse 图 1 EXCEL 工作簿中分别建立正算和反算两个工作表 Fig.1 Separate establishment of positive and reverse calculation work table in EXCEL workbook 第 22 卷第 1 期 2012 年 3 月 安徽地质 Vol.22 No.1 Geology of Anhui March 2012 文章编号:1005 - 6157(2012)01 - 70 - 03 用 EXCEL 处理高斯坐标计算问题 吕建国 （安徽省煤田地质局第三勘探队，安徽 宿州 234000） 摘 要：本文介绍了使用 EXCEL 处理高斯坐标的正算、反算问题，着重阐述了计算公式的编制，具有使用方便、实 用性较强的特点。 关键词：EXCEL；高斯坐标；地理坐标；正算；反算 中国分类号：TP317.3 文献标识码：A 地理坐标的单位一般是由度分秒组成的，而 EX? CEL 中角度单位是弧度制，所以首先必须将以度分秒 为单位的坐标转换为以弧度为单位的坐标。在 EX? CEL 工作簿中分别建立正算和反算两个工作表，如下 图 1 所示，在 B2、C2 两个单元格中分别输入纬度和 经度坐标，特别注意，输入纬度和经度的坐标是以度 为整数位，小数点后两位为十分位和分位，十分位和 分位是零的用 0 补齐，分位后为秒位。为了便于说 明，对建立的每个单元格公式按表 1 进行说明： 2 反算 高斯投影反算是由高斯平面坐标 （X，Y） 求地 理坐标 （L，B）。 建立反算工作表后，在 B2、C2 单元格中分别输 入 X、Y 值，象正算一样，按照高斯投影反算公式逐 个建立单元格公式，如表 2 所示： 0 引言 在地质勘探、地形图测绘和公路、铁路等工程 施工中，经常涉及地理坐标与高斯坐标的换算，即 高斯投影坐标计算，包括高斯投影正算和高斯投影 反算。无论是正算还是反算，为了确保精度要求， 常常会用到高阶导数 （一般算至六阶） 按照公式计 算，或者使用 《高斯、克吕格投影计算表》 进行计 算，不管是利用公式计算还是查表计算，计算都比 较复杂，而且精度较低，容易出错，计算速度慢。 随 着 计 算 机 的 广 泛 应 用 ， 利 用 EXCEL 电 子 表 格 计 算，不仅速度快而且精度高。 1 正算 高斯投影正算是由地理坐标 （L，B） 求高斯平 面坐标 （X，Y）。 收稿日期：2010-08-16 (H2,6)+0.004*POWER(H2,8)表 2 高斯反算单元格公式及说明Table 2 Notes on Gauss reverse calculation cell formula*PI()/180)*COS((E2/3600)*PI()/180)*206264.8063分秒的形式。分秒的形式。充，正算时对 D 至 N 列单元格填充，反算时对 D 至 R列单元格填充，计算时只需在 B 列、C 列输入已知参数，在正算工作表中，从 O 列、P 列中分别得到相应的X、Y 值，在反算工作表中，可以在 Q 列、R 列分别得到相应的地理坐标 (H2,6)+0.004*POWER(H2,8) 表 2 高斯反算单元格公式及说明 Table 2 Notes on Gauss reverse calculation cell formula *PI()/180)*COS((E2/3600)*PI()/180)*206264.8063 分秒的形式。 分秒的形式。 充，正算时对 D 至 N 列单元格填充，反算时对 D 至 R 列单元格填充，计算时只需在 B 列、C 列输入已知参 数，在正算工作表中，从 O 列、P 列中分别得到相应的 X、Y 值，在反算工作表中，可以在 Q 列、R 列分别得 到相应的地理坐标 B、L 的值。对单元格填充可以很方 便地达到批量计算的目的，大大提高了工作效率。 （3） 采用 EXCEL 公式计算，不仅计算方便，出 错率很低，而且精度较高，正算时其误差不超过± 0.5m，反算时地理坐标的误差达到±0.001″ 。 3 几点说明 （1） 3°投影时，可以根据中央子午线计算公式 对正算工作表中 F2 单元格和反算工作表中 R2 单元 格作如下修改即可利用： F2=INT((E2+1.5)/3) R2=D2*3 + INT(P2/3600) + INT((P2/3600-INT(P2/ 3600))*60)/100+(P2/60-INT(P2/60))*60/10000 （2） 无论是正