2. 矿山测量坐标系与高程基准
无人机影像测量得到的成果,最终要落到“坐标”和“高程”这两个数值体系上。影像本身只是像素,只有赋予正确的平面坐标系与高程基准,空三加密、DOM、DSM、点云和方量计算才具有工程意义。矿山测量常见的麻烦往往不在飞行或建模,而在基准选错、转换参数用错,导致成果整体偏移或高程系统性偏差。本节围绕常用坐标系、高程基准和坐标转换三件事展开。
一、平面坐标系:国家2000与地方坐标
我国现行国家大地坐标系为2000国家大地坐标系(CGCS2000),它是一种地心坐标系,椭球参数为CGCS2000椭球,与WGS84在厘米级上接近但不等同。无人机RTK/PPK获取的原始定位通常为CGCS2000或WGS84经纬度,经投影后得到平面坐标。矿山成果若需上报自然资源部门、与国土“一张图”套合,应统一采用CGCS2000下的高斯-克吕格投影,并明确中央子午线与投影带号(3°带或6°带)。
地方坐标则是矿山为满足本地生产而建立的独立平面坐标系,常见做法是选取矿区中央子午线做高斯投影,或直接采用矿区假定坐标,并对长度进行投影变形改正,使图上距离与实地距离尽量一致。地方坐标的优点是与采区、台阶、运输道路等生产图件贴合紧密;缺点是脱离国家体系,跨部门交换时必须转换。实际作业中要牢记:同一批成果不能混用两套平面坐标,DOM、DSM、点云、控制点必须同基准。
二、高程基准:1985国家高程与矿山独立高程
高程方面,国家采用1985国家高程基准,以青岛验潮站资料确定的黄海平均海水面为起算面。矿山生产图件、台阶标高、储量计算常使用矿山独立高程系统,即以矿区某固定点假定高程为起算,或与国家高程保持一个固定差值。两者之间通常只差一个常数,但这个常数必须通过联测确定并记录在案。
需要区分的是大地高与正常高:无人机RTK直接给出的是大地高(相对椭球面),而工程使用的高程是正常高(相对似大地水准面),二者相差高程异常ζ。若测区高程异常未改正,高程可能出现米级系统偏差。因此,用RTK/PPK数据做高程时,应使用测区似大地水准面模型或至少用已知水准点做拟合改正。
三、坐标转换怎么做
平面坐标转换常用四参数(平移ΔX、ΔY,旋转ε,尺度K)或七参数(增加Z平移与三个旋转)。四参数适用于小范围、两平面坐标系间的转换;七参数用于不同椭球基准间的三维转换。转换参数必须由公共点求解:在测区内选取至少3个(四参数)或3个以上(七参数)同时具有两套坐标的控制点,用最小二乘求参,再用未参与解算的检核点验证残差。残差应控制在厘米级,若明显偏大,说明公共点分布不均或存在粗差。
高程转换相对简单:若已知矿山独立高程与国家高程的固定差值ΔH,则 H矿山 = H国家 + ΔH。但要注意,若独立高程是假定起算且未做水准联测,差值可能随距离变化,此时应在矿区布设若干水准点分段拟合,而不是全矿用一个常数。
四、一个教学示例
假设某露天矿需将无人机成果从CGCS2000转换到矿山地方坐标。选取矿区周边4个控制点,其中3点求四参数、1点检核。示例参数(合理假设,非实测):ΔX = −1250.342 m,ΔY = 860.117 m,旋转 ε = 0.00231 rad,尺度 K = 1.000012。用检核点计算转换后坐标与原地方坐标之差,若平面残差约0.02 m,可认为参数可用;若某点残差达0.15 m,应排查该点是否被碰动或坐标录入错误。高程方面,联测得ΔH = −1.236 m,则国家高程1200.000 m对应矿山高程1198.764 m。
总之,矿山无人机测量的基准问题,核心是“平面选对带、高程联测准、转换有检核”。把坐标系、高程基准和转换参数在技术设计书中写清楚,并在成果元数据中注明,才能保证后续方量计算、变形监测和成果汇交的可靠性。