3. 坐标系统与投影基础
在把地质钻孔CSV转成Shapefile之前，必须先回答一个看似简单却极易出错的问题：钻孔的X、Y坐标到底属于哪一套坐标系统？如果这一步判断错误，后续的投影校验、图层叠加、图幅拼接都会建立在错误的基础上。本节围绕三个核心概念展开：地理坐标系与投影坐标系的区别、常见大地基准的识别，以及投影参数（中央经线、比例因子、假东假北）的解读。
一、地理坐标系与投影坐标系：球面位置与平面位置
地理坐标系（Geographic Coordinate System，GCS）用经度、纬度（以及可选的高程）描述地球表面点的位置，单位是度（或度分秒）。它的本质是“角度坐标”，依附于一个旋转椭球体和一个大地基准。经纬度本身不能直接用于量算平面距离，因为1度经度在赤道约111 km，而在纬度60°处只有约55.5 km。
投影坐标系（Projected Coordinate System，PCS）则是把椭球面按某种数学规则展开到平面上，得到以米（或英尺）为单位的X、Y坐标。它由“地理坐标系 + 投影方法 + 投影参数”三部分共同定义。地质工作中常见的投影是高斯-克吕格投影（横轴墨卡托的一种），我国国家基本比例尺地形图普遍采用3°带或6°分带的高斯-克吕格平面直角坐标。
对钻孔数据而言，最直接的判断依据是数值量级：经纬度坐标的X（经度）通常在70～140之间、Y（纬度）在15～55之间，且带小数；而投影坐标的X、Y通常是6位或7位以上的大数，例如东坐标约在数百千米量级、北坐标约在数百万米量级。若CSV中X=38421567.32、Y=2954321.87，基本可判定为投影坐标，且东坐标前缀“38”往往对应3°带的带号。
二、常见大地基准：CGCS2000、WGS84、西安80
大地基准定义了椭球体的形状、大小及其与地球的相对定位。同一经纬度在不同基准下对应的实际地面点可能相差几十米到上百米，因此基准识别是投影校验的第一道关口。
- CGCS2000（2000国家大地坐标系）：我国现行法定国家大地坐标系，采用ITRF2000框架、2000.0历元，椭球参数为长半轴6378137 m、扁率1/298.257222101。新采集的地质钻孔数据、国家基本比例尺地形图应优先使用CGCS2000。
- WGS84：GPS/GNSS直接输出的全球坐标系，椭球长半轴同为6378137 m，扁率1/298.257223563，与CGCS2000在厘米级上接近，但在严格制图与成果汇交中不能混用。
- 西安80：以IAG75椭球为基础的参心坐标系，长半轴6378140 m、扁率1/298.257。大量历史地质资料、老钻孔台账仍为西安80坐标，转换到CGCS2000需借助七参数或经国家转换参数。
实际判断时，可结合数据来源与数值特征：若CSV来自近年GNSS采集，多为CGCS2000或WGS84；若来自上世纪地质队台账，且坐标数值与现行图幅存在系统性偏移，应怀疑西安80或北京54。切勿仅凭“看起来差不多”就认定基准一致。
三、投影参数解读：中央经线、比例因子、假东假北
以高斯-克吕格3°带为例，带号n与中央经线的关系为L₀ = 3n。若东坐标前缀为38，则n=38，中央经线为114°E，覆盖范围约112.5°E～115.5°E。6°带的中央经线为L₀ = 6n − 3，例如n=19时L₀=111°E。中央经线决定了投影的对称轴，离中央经线越远，长度变形越大，因此钻孔数据应落在其所属带的有效范围内。
比例因子k₀表示投影面上长度与椭球面长度的比值。高斯-克吕格投影在中央经线上k₀=1，向两侧逐渐增大；UTM投影则取k₀=0.9996，使中央经线附近略有负变形、整体变形更均匀。若把UTM参数误用于高斯成果，长度量算会出现约0.04%的系统偏差，在千米级剖面上可达数十厘米。
假东（False Easting）与假北（False Northing）是为了避免出现负坐标而人为加的常数。我国高斯-克吕格投影通常给东坐标加500000 m，使中央经线以西的负值变为正值；北半球假北取0，南半球常取10000000 m。因此，一条东坐标“38421567.32”可拆解为：带号38 + 假东500000 + 实际偏移−78432.68 m，说明该点位于中央经线以西约78 km处。
四、落到钻孔CSV上的判断流程
拿到一份钻孔CSV，建议按以下顺序核对：先看X、Y数值量级，区分经纬度与平面坐标；若为平面坐标，再看东坐标前缀是否与工作区经度对应的带号一致；然后核对数据来源所声明的大地基准；最后在GIS中加载已知控制点或地形图进行套合检查。若套合后整体偏移稳定在数十米，多为基准差异；若偏移随位置变化，则可能是带号或中央经线选错。只有把这三层信息——坐标系类型、大地基准、投影参数——都确认清楚，后续的Shapefile生成与投影校验才有可靠前提。