地质工作中平面坐标的辨识与转换.docVIP

地质工作中平面坐标的辨识与转换 坐标是地质工作中必不可少的因素。实际工作中有时候业主提供的资料不足以让我们的工作顺利展开，比如：某条公路压覆矿调查项目，业主一般会提供一个公路走向的CAD图件，图件里的坐标格式为（58XXXX,334XXXX），而国土资源厅矿业权管理系统查询坐标格式为Y坐标带有带号的完整坐标，此时就要求我们根据实际工作经验正确辨识并将已有资料转化为需要的坐标格式。笔者根据实际工作经历及同事在工作中咨询的一些问题结合网上收集整合资料跟大家分享下地质工作中坐标问题的解决方法。 1、简单概念 ⑴目前我们使用的坐标为大地坐标经过高斯—克吕格投影处理后得到的平面直角坐标。 ⑵平面直角坐标分带是将地球视为一个椭球面，按3°或者6°分割并展开投影到圆柱上。湖北3°带Y坐标带号为36—39，6°带Y坐标带号为19—20。 ⑶北京1954和西安1980为2个不同的椭球标准，湖北地区2个坐标系数值相差10到20米，而业主通过设计院提供的坐标一般为北京1954坐标，矿产资源管理系统使用的是西安1980坐标，所以我们查矿业权要先把北京坐标系转化为西安坐标系，不过差距并不是很大。 2、坐标辨识 图1 坐标分带图 ⑴3°分带和6°分带都是以东经3°为起始带的中央子午线，然后往后面顺延（见图1）。 ⑵在一个6°带区域内，以中央子午线为纵坐标轴（X轴），赤道为横坐标轴（Y轴），这样中央子午线左边的横坐标就为负的。为了解决这个问题，把北半球所有的纵坐标轴（X轴）都左移500公里（3°分带以及1.5°分带坐标都使用这个设定），也就是数值+500000，所以中央子午线附近的横坐标Y值通常为XX500000±（XX为带号）。 图2 X轴左移500公里以解决横坐标负值 ⑶纵坐标X都是以赤道为横坐标轴投影计算的，也就是横坐标轴固定不变，所以3°分带和6°分带以及1.5°分带X坐标是没有变化的。 ⑷平面直角坐标的格式为（X1X2X3X4X5X6X7, Y1Y2Y3Y4Y5Y6Y7Y8）, Y1Y2为带号，Y3Y4Y5Y6Y7Y8为投影面上坐标点到中央子午线（X轴）的矢量距离（米）+500000，X1X2X3X4X5X6X7为投影面上坐标点到赤道的距离。 ⑸3°带带号Y1Y2=中央子午线经度/3，6°带带号Y1Y2=（中央子午线经度+3）/6。 图2 湖北省各地区经度线及中央子午线 ⑹湖北省区域内6°分带的中央子午线有111°一条，3°分带的中央子午线有111°和114°， 对应的带号为19，区域为108°—114°，也就是图2中114°经线左边6°带坐标其Y坐标前2位必定为19， 114°经线右边6°坐标其Y坐标前2位为20。 ⑺3°分带的中央子午线有111°和114°，对应的带号分别为37和38，区域分别为109.5°—112.5°和112.5°—115.5°，即图2中粉红色虚线划定的区域，3°坐标中109.5°经线左边其Y坐标前2位为36，109.5°—112.5°区间Y前2位为37，依次往后112.5°—115.5°区间Y前2位为38，115.5°右边Y前2位为39。其对应的带号及区间见表1。 表1 湖北省中央子午线及其区域分带表 坐标带 省内中央子午线 省内区域 带号 Y前2位 6°带 111° 114°左边 19 19 114°右边 20 20 3°带 109.5°左边 36 36 111° 109.5°—112.5° 37 37 114° 112.5°—115.5° 38 38 115.5°右边 39 39 ⑻1.5°带坐标主要是为了提高精度及解决3°带坐标跨带问题，除了3°带坐标原有中央子午线外，将3°带分带线额外作为中央子午线，这样解决了跨带线附近敏感区的精度问题，在地方城市测量中用的较多。如果在3°分带线附近（如112.5°）Y坐标数值在500000左右，就是1.5°分带坐标，我们只需设置中央子午线（112.5°）在MAPGIS里面转换为3°或者6°坐标即可。 3、坐标转换 MAPGIS和coordinate转化软件都能实现坐标间的转换。 ⑴地理坐标（经纬度）与平面坐标的转换 在实际工作中，由于地理坐标（度分秒）误差太大，所以一般用于面积大、平面形状简单的探矿权区域。可根据投影原理计算出转换公式，再写个C语言小程序实现，转换公式为一个椭圆柱面（投影面）与通过原点的射线（经纬度）之交点，比较复杂。推荐软件转换。 ⑵换带转换 与6°带中央子午线相同的3°带区域坐标，除了Y坐标带号外，其他完全相同，带号按照表1由中央子午线计算。其他区的6°带与3°带坐标转换可由坐标转换软件coordinate以及MAPGIS实现，转换中非常要注意的是设置好中央子午线，带号不在转换之列，需单独计算。 图3 coordinate坐标换