企业官网建设流程全解析

简介：面向ArcGIS用户的WGS84经纬度坐标转CGCS2000投影坐标系操作说明文档，适合测绘、GIS数据处理及国土相关技术人员参考。文档完整梳理了从奥维地图导出shp文件，到ArcGIS中进行投影与变换，再到通过ITRF2000过渡并最终转换为CGCS2000坐标系的全流程，涵盖WGS84、ITRF2000、CGCS2000三种坐标系的椭球参数与转换逻辑，以及带号选择、3度带/6度带中央经线计算、坐标字段添加与十进制单位切换等关键细节。包体仅含1个doc文件，大小3.06MB，为图文步骤说明，可直接对照ArcToolbox操作。该文档已有646人学习下载，适合需要处理不同坐标系数据转换、规避常见误差的GIS使用者，是一份实用且可复用的操作笔记。对于将奥维兴趣点、测量成果等外部数据转入国家2000坐标系的应用场景，按此文档操作可减少转换错误；步骤按菜单层级逐步展开，也便于零基础读者对照完成。

1. 经纬度坐标系转CGCS2000，差的不是度数，是基准

许多人以为经纬度是全球通用的，拿到GPS坐标直接当成CGCS2000用。实际上普通GPS接收机默认输出WGS84经纬度，而CGCS2000是一个独立的地心坐标参考框架，二者椭球参数、参考历元和框架实现都有细微差异，造成的平面偏移在多数地区有0.5～1.5米。这在小比例尺地图上不明显，但施工放样、地籍测量、CAD与GIS拼图时足以让线位错位。这篇文章只解决一件事：如何把经纬度坐标系的数据正确转成CGCS2000坐标系，包括原理、七参数计算、复用工具和最终验证。适合测绘内业、GIS开发者、CAD转GIS的数据处理人员。

2. 经纬度转CGCS2000前，先把椭球、基准和投影分开算

2.1 经纬度本质是大地坐标，必须先锁定参考椭球

经纬度(B,L,H)本身只是表达方式，不是基准。同一坐标数值，放在不同椭球上，对应的空间位置不同。CGCS2000使用CGCS2000椭球，长半轴a=6378137m，扁率f=1/298.257222101，而WGS84椭球扁率为1/298.257223563，长半轴相同，扁率差很小但存在。转换第一步不是套公式，而是确认源和目标各处于哪个椭球/框架。一般从在线地图或手持GPS获取的经纬度可视为WGS84；国土、规划数据中的经纬度很可能已经是CGCS2000。

要把经纬度转成空间直角坐标，常见做法是先用椭球参数计算卯酉圈曲率半径N，再按下面关系式分解：

X = (N + H) * cos(B) * cos(L) Y = (N + H) * cos(B) * sin(L) Z = (N * (1 - e^2) + H) * sin(B)

其中e^2 = 2f - f^2，H是大地高（不是海拔）。这段公式可以直接用Python实现：

def geodetic_to_cartesian(lat_deg, lon_deg, h_m, a=6378137.0, rf=298.257222101): from math import radians, sin, cos, sqrt f = 1.0 / rf e2 = 2 * f - f * f b = radians(lat_deg) l = radians(lon_deg) n = a / sqrt(1 - e2 * sin(b) ** 2) x = (n + h_m) * cos(b) * cos(l) y = (n + h_m) * cos(b) * sin(l) z = (n * (1 - e2) + h_m) * sin(b) return x, y, z

这里把纬度、经度、大地高传进去，默认使用CGCS2000椭球。函数返回的是以地心为原点的XYZ直角坐标，单位米。rf是扁率倒数，e2是第一偏心率平方，这两个是椭球转换里最容易抄错的地方。如果你手里的是海拔高程，还要在转换前加大地水准面差距，否则转出的XYZ会整体偏高或偏低，影响后续七参数。

2.2 CGCS2000和WGS84只差一个基准转换

椭球参数只决定了数学形状，不决定地球框架。CGCS2000采用ITRF97参考框架，参考历元2000.0；WGS84在G1762版本后与ITRF对齐到厘米级。实际工作中，很多地方直接忽略转换，导致局部有规律性偏移。如果项目精度要求优于1米，就必须做基准转换。国内常用的是布尔沙七参数模型。

七参数包括平移参数DX、DY、DZ（米），旋转参数RX、RY、RZ（角秒或弧度），以及尺度因子m（ppm）。转换关系可以写成矩阵形式，小角度情况下常见做法是展开为：

X2 = DX + (1 + m) * X + RZ * Y - RY * Z Y2 = DY + (1 + m) * Y - RZ * X + RX * Z Z2 = DZ + (1 + m) * Z + RY * X - RX * Y

也就是说，把源坐标X、Y、Z做一次线性变换，得到目标框架下的坐标。注意旋转参数单位要换算成弧度，而且公式的旋转符号约定在不同地区可能相反，需要小范围试算确认。

七参数的来源要非常谨慎。示例参数如下，仅为写法示意，不能用于生产：

DX = 1.20 m DY = -0.80 m DZ = -1.50 m RX = -0.02" RY = -0.01" RZ = 0.03" m = 0.5 ppm

提示：生产项目的七参数必须向当地自然资源或测绘主管部门索取，或使用GNSS连续运行参考站提供的本地参数，不要从网上随便抄一组。参数对不上时，结果会偏离数百米。

2.3 七参数怎么套进经纬度里

手工搞定的套路是：源经纬度转XYZ -> 应用七参数得到目标XYZ -> 目标XYZ反算目标经纬度。如果目标要投影平面坐标，再对目标经纬度做高斯-克吕格投影。整个过程只有三步，看起来不复杂，但每一步都依赖椭球参数和投影带号，其中一项选错后面全错。

反算XYZ到经纬度需要迭代，因为大地纬度出现在N的计算式里。常见做法是用初始纬度迭代几次，收敛到毫米级就停。

2.4 转到平面时要区分3度带和6度带

当目标CGCS2000坐标是平面坐标（如CAD图纸里的x=3456789.12, y=20612345.67），实际是高斯-克吕格投影平面坐标。Y坐标前两位数“20”是6度带带号，实际横坐标是612345.67（东偏移500km后）。3度带或6度带的中央经线计算公式：

• 6度带：L0 = 6 * n - 3，例如带号20，中央经线L0 = 117°E
• 3度带：L0 = 3 * n，例如带号39，中央经线L0 = 117°E

下表是东经114°附近几个带号的中央经线对照：

选择带号时按目标点经度算，不要照抄图纸里任意坐标的带号，尤其跨带图幅要额外说明。很多CAD转GIS的“6位坐标”问题，本质上就是Y坐标省略了带号，后面第4章会单独讲。

3. 用 pyproj 和手动七参数把经纬度转成 CGCS2000

3.1 首选 pyproj：两行代码先解决 WGS84 到 CGCS2000 的经纬度

pyproj是GDAL生态里最顺手的坐标转换库。如果源数据是WGS84经纬度，目标只需要CGCS2000经纬度，直接定义Transformer即可。CGCS2000的地理坐标EPSG代码是4490，WGS84地理坐标是4326。

from pyproj import Transformer t = Transformer.from_crs(4326, 4490, always_xy=True) # 输入经度、纬度，输出CGCS2000经纬度 lon, lat = t.transform(117.123456, 39.654321) print(lon, lat)

always_xy=True表示输入输出都是经度在前、纬度在后，避免和传统纬度在前习惯混在一起。转换时pyproj会按两个坐标系对应的椭球和基准做换算，默认没有使用七参数时，它做的是“忽略基准面差异”的转换，极端情况下会产生米级偏差。如果你的数据源本身是CGCS2000下的经纬度，直接跳过这步即可。

3.2 立即转平面坐标：高斯投影的 PROJ 字符串

很多项目最终要的是CGCS2000平面坐标。比如把WGS84经纬度转成6度带20带的高斯平面坐标，可以自定义一个PROJ坐标系，使用+proj=tmerc，中央经线设置为117°E。这样既能同时完成基准转换和投影：

from pyproj import CRS, Transformer source = CRS.from_epsg(4326) target = CRS.from_proj("+proj=tmerc +lat_0=0 +lon_0=117 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs") t = Transformer.from_crs(source, target, always_xy=True) x, y = t.transform(117.123456, 39.654321) print(x, y)

这里的x是高斯平面北坐标，y是横坐标，注意+x_0=500000表示横坐标加500公里偏移，输出Y在500000左右。如果你需要带带号，在Y前拼接带号20，例如结果y=206123456.78。CGCS2000椭球比GRS80扁率略大，但在毫米级精度要求不高的场景下，用+ellps=GRS80通常不影响工程应用；若需要严格匹配，把+ellps=GRS80替换为+a=6378137 +rf=298.257222101。

3.3 手动七参数转换：不依赖 PROJ 也能算

在某些离线环境里没有pyproj，或者甲方只给了一组七参数，我一般会直接手写转换函数。完整过程是：源BLH转XYZ，XYZ过七参数，再反算目标BLH。下面的函数实现了三步：

import math def blh_to_cgcs2000(lat_deg, lon_deg, h_m, dx, dy, dz, rx, ry, rz, m_ppm, rf_src, rf_dst): # 1. 源椭球BLH -> XYZ x, y, z = geodetic_to_cartesian(lat_deg, lon_deg, h_m, rf=rf_src) # 2. 七参数转换 rx = math.radians(rx / 3600.0) # 角秒转弧度 ry = math.radians(ry / 3600.0) rz = math.radians(rz / 3600.0) s = 1 + m_ppm * 1e-6 x2 = dx + s * (x + rz * y - ry * z) y2 = dy + s * (-rz * x + y + rx * z) z2 = dz + s * (ry * x - rx * y + z) # 3. XYZ -> BLH（迭代求纬度） lon = math.atan2(y2, x2) p = math.sqrt(x2 * x2 + y2 * y2) lat = math.atan2(z2, p * (1 - 2.0 / rf_dst)) for _ in range(5): n = 6378137.0 / math.sqrt(1 - (1 - 1.0 / rf_dst) ** 2 * math.sin(lat) ** 2) h = p / math.cos(lat) - n lat = math.atan2(z2, p * (1 - (1 - 1.0 / rf_dst) * n / (n + h))) lon_deg = math.degrees(lon) lat_deg = math.degrees(lat) return lat_deg, lon_deg, h

七参数中的旋转单位在这里按角秒处理，角秒转弧度要除以3600再乘π/180，rx参数传入时是角秒值。第三步反算纬度用了迭代近似，5次循环已足够收敛到毫米级。这个函数把源椭球扁率rf_src和目标椭球扁率rf_dst都显式传进去，避免把大地高当海拔。

3.4 转换结果对拍表：到底差多少米

下面用一组示意数据演示转换前后差异。输入WGS84经纬度，CGCS2000近似等于其目标值；平面坐标只是示例说明量级。

表格里的差值是示意性的，实际差异由七参数决定。判断结果是否合理，关键是看Y坐标是否落在带号对应的中央经线附近，以及同一坐标在不同方法下回算残差是否小于0.01米。

4. 不写代码的路径：ogr2ogr、GeoHey 在线转换与 CAD 的 6 位坐标

4.1 用 ogr2ogr 批量转矢量文件

如果你手头是Shapefile或GeoPackage，不想写Python，用GDAL自带的ogr2ogr最直接。把矢量数据从WGS84经纬度转成CGCS2000地理坐标，命令：

ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:4490 output.shp input.shp

-s_srs指定源坐标系，-t_srs指定目标坐标系。要转到高斯平面坐标，需要写更完整的参数：

ogr2ogr -s_srs EPSG:4326 -t_srs "+proj=tmerc +lat_0=0 +lon_0=117 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs" output.shp input.shp

批量转换时建议先转一份小数据验证字段和坐标，再用循环处理整个目录。ogr2ogr默认会重写OGR字段，不会改变属性表结构。

4.2 GeoHey 在线坐标转换做单点抽查

需要快速核对单个经纬度时，在线坐标转换是很多人的习惯做法。GeoHey在线坐标转换工具里，通常要先选源坐标系和目标坐标系。这里有个容易忽略的细节：工具里如果没有“CGCS2000地理位置”选项，可以选“CGCS2000 / 3-degree Gauss-Kruger CM 117E”这类投影坐标系，结果会直接是平面坐标。单点抽查只用来排查方向性错误，不建议大批量生产。在线工具大多默认不公开具体七参数，所以要确认结果是否满足你项目的平面精度要求。

4.3 CAD到GIS：6位坐标其实是投影坐标没带带号

CAD图纸经常出现“6位坐标转换”的疑问：CAD里标出的点比如X=3456789.12，Y=612345.67，导入GIS后跑到别的位置。常见原因是这个Y是高斯平面坐标去掉了带号。比如6度带20带的实际横坐标应为20612345.67，Y=612345.67只是EASTING的十位到个位部分。对应处理方法是先补带号再转经纬度，或者补带号后直接定义成CGCS2000平面坐标。

我一般先看图纸说明：如果坐标范围只有6位，很可能是省略了中央经线前的带号；如果Y值大于500000则横坐标正常。补带号后在ArcGIS中定义坐标时，选择CGCS2000 6度带对应的投影坐标系，例如20带。如果图纸是3度带，就把带号换成39之类的3度带带号。这样处理后CAD到GIS的6位坐标问题就解决了一半，剩下需要确认长度单位是米还是毫米。

5. 三个必查参数和一组回算验证转换结果

5.1 必查参数：源椭球、目标框架和转换参数

转换结果异常时，先检查三个参数：源坐标的椭球和基准、目标坐标系有没有带投影带号、七参数是否匹配区域。按顺序排查，多数问题会浮出来。源坐标如果是RTK导出的经纬度，通常已落入CGCS2000；如果是手机GPS，基本是WGS84。目标坐标如果要求带带号的Y值，就不要选不带带号的定义，否则绘图软件会把它当普通坐标。

5.2 用已知控制点做残差回算

最可靠的验证方法是回算已知点。选一个当地已知的CGCS2000控制点（平面坐标或经纬度），把它和待转换点用同一套流程处理。具体做法是：先把你的参数和代码转出一个值，再用反函数转回原坐标系，比较初始值和回算值的差。也可以取两个已知点，一个做参数拟合，一个做独立验证。编写一个小脚本计算残差：

for pt in known_points: lat2, lon2 = transform(pt.lat, pt.lon) dx_m = (lon2 - pt.lon_cgcs) * 111000 * math.cos(math.radians(pt.lat)) dy_m = (lat2 - pt.lat_cgcs) * 110946 print(pt.name, dx_m, dy_m)

这里用近似公式把经度纬度差换算成平面米数，适合快速检查。正常残差应在厘米到分米级；如果出现米级以上偏差，回到5.1重新核对参数。验证完成后，保留脚本和参数表，作为交付文档的一部分。

本文还有配套的精品资源，点击获取