2. 钻孔点位坐标基础
在钻孔点位批量转换工作中,最先要弄清楚的并不是软件操作,而是“坐标到底属于哪一套系统”。同一组数字,如果坐标系判断错误,转换后可能整体偏移几十米甚至上百米,钻孔孔位就会落到错误的地段上。因此,本节从坐标的两种基本表达方式讲起,再落到我国常见的几套坐标系及其识别方法。
一、地理坐标系与投影坐标系:两种表达,一套地球
地理坐标系(Geographic Coordinate System,GCS)用经度、纬度(以及大地高)表示地面点的位置,单位是度(或度分秒)。它的本质是“椭球面上的角度坐标”:经度是过该点的子午面与本初子午面的夹角,纬度是该点法线与赤道面的夹角。地理坐标不涉及平面长度,任意两点间的距离需要按椭球面弧长计算,因此不适合直接在图上量距、放样。
投影坐标系(Projected Coordinate System,PCS)则是把椭球面按一定数学法则“展开”到平面上,用东坐标(Easting,常记 X 或 y)、北坐标(Northing,常记 Y 或 x)表示,单位是米。我国大比例尺测绘最常用的是高斯-克吕格投影(横轴墨卡托投影的一种),按经差分带,3° 带共 120 带、6° 带共 60 带,每带设一条中央子午线,投影后中央子午线为直线且长度比不变。
高斯投影有两个必须记住的约定:其一,为避免出现负值,北坐标(Y)在赤道处加 500 km 常数;其二,东坐标(X)前加带号,例如 3° 带第 38 带写作 38500000 一类的形式。判断一个平面坐标是否“加了带号”,看整数位数即可:6 位为不带带号,8 位通常为带号+6 位。这一条在批量转换中极为实用。
二、我国常用坐标系的区别
我国测绘历史上先后使用过北京54、西安80、CGCS2000 三套国家大地坐标系,它们最核心的差别在于所采用的参考椭球及定位定向参数不同,因而同一地面点在其中的经纬度和平面坐标都不相同。
| 坐标系 | 参考椭球 | 长半轴 a(m) | 扁率倒数 1/f | 特点与现状 |
|---|---|---|---|---|
| 北京54 | 克拉索夫斯基椭球 | 6378245 | 298.3 | 1954 年建立,参心坐标系,与我国大地水准面吻合较差,已基本停用 |
| 西安80 | IAG75 椭球 | 6378140 | 298.257 | 1980 年建立,参心坐标系,定位定向明显优于北京54,历史成果量大 |
| CGCS2000 | CGCS2000 椭球 | 6378137 | 298.257222101 | 2000 国家大地坐标系,地心坐标系,现行国家标准,新建项目应优先采用 |
需要强调的是,三套坐标系之间的差异不是简单的“平移一个常数”,而是椭球参数、原点定位、定向共同作用的结果。北京54 与西安80 之间、西安80 与 CGCS2000 之间,通常需要通过公共点求取七参数(或四参数)进行转换,且转换参数具有明显的区域性和精度限制。工程上常见的做法是:收集测区内至少 3 个(最好 4~6 个)同时具备两套坐标的控制点,用最小二乘求解转换参数,再用残差检验转换质量。
此外还有一类“地方坐标系”,例如某市独立坐标系、某矿区坐标系、工程独立坐标系。它们往往是在国家坐标系基础上做平移、旋转、缩放,或采用自定义中央子午线和抵偿高程面,目的是减小投影变形、便于施工放样。地方坐标系的坐标值常常“看起来很正常”,但和国家坐标对不上,这是钻孔点位转换中最容易出错的地方。
三、识别方法:从元数据、数值特征到实地校核
实际工作中,拿到一批钻孔点位坐标,应按“先查元数据、再看数值特征、最后用控制点验证”的顺序判断坐标系。
- 查元数据与成果说明:优先查看原始测量记录、成果表表头、GIS 数据属性、CAD 图纸说明,寻找“1980 西安坐标系”“CGCS2000”“3° 带第 XX 带”“中央子午线 XXX°”等字样。这是最可靠的依据。
- 看数值量级与位数:经纬度一般在 70°~140°E、15°~55°N 之间,带小数;平面坐标 X 通常为 6 位或 8 位整数,Y 为 7 位左右。若 X 为 8 位且前两位在 13~23(6° 带)或 25~45(3° 带)之间,说明带号已加。
- 反算中央子午线:由带号可反推中央子午线,6° 带 L0 = 6N − 3,3° 带 L0 = 3N。若坐标不带带号,则需结合点位经度判断所属带号。
- 用已知控制点比对:选取 1~2 个已知国家点,将待判坐标与其比较,若差值呈现系统性平移或旋转,说明坐标系或带号有误。
- 图形叠加检查:把钻孔点位叠加到带地理配准的底图或已有正确坐标的勘探线上,观察是否整体偏移、旋转或缩放,这是最直观的“肉眼验证”。
四、一个识别示例
假设某钻孔成果表中给出点位坐标:X = 38524678.35,Y = 4321987.62。按前述规则判断:X 为 8 位,前两位“38”可视为 3° 带带号,则中央子午线 L0 = 3 × 38 = 114°;去掉带号后东坐标为 524678.35,减去 500 km 得 24678.35 m,说明该点位于中央子午线以东约 24.7 km;Y 为 7 位,符合含假东移的北坐标形式。由此可初步判定这是 3° 带、中央子午线 114° 的投影坐标,但仍不能确定是西安80 还是 CGCS2000——两者在该区域的平面坐标可能只差几十厘米到数米,必须依靠成果说明或控制点比对确认。
若同一批数据中另有一些点 X 只有 6 位,如 524678.35,则说明该批数据未加带号,转换前必须统一补带号或统一按不带带号处理,否则批量转换后点位会整体错带。这类“同批数据带号不一致”的问题,在钻孔点位批量处理中非常常见,建议在转换前先做一次位数统计与异常点筛查。
综上,坐标系识别是钻孔点位批量转换的第一道关口。正确做法是:先确认地理坐标还是投影坐标,再确认属于哪套大地坐标系、哪个投影带、中央子午线是多少,最后用控制点或底图叠加验证。只有这四步都落实,后续的批量转换才有可靠的前提。