GIS坐标系与面积计算类型对照
一句话总结 ：面积计算结果的正确性，取决于你用什么坐标系 + 用什么数学模型。
下表一次性理清所有组合关系。
一、核心概念速查
| 面积类型 | 英文 | 计算基准 | 要求数据在 | 是否受投影变形影响 | 是否考虑地球曲率 |
|---|---|---|---|---|---|
| 投影面积 | Projected / Planar Area | 平面直角坐标 (X,Y) | 投影坐标系(PCS) | ✅ 是 | ❌ 否 |
| 椭球面积 | Ellipsoidal Area | 参考椭球面 | 地理坐标系(GCS) | ❌ 否 | ✅ 是 |
| 测地面积 | Geodesic Area | 椭球面(测地线) | 地理坐标系(GCS) | ❌ 否 | ✅ 是（最严格） |
| °² 面积 | Degree-Squared Area | 经纬度当平面 | 地理坐标系(GCS) | ❌（但无物理意义） | ❌ 否 |
注意 ：在 ArcGIS / QGIS 等主流软件中，椭球面积 ≈ 测地面积，两者视为同一结果。
二、坐标系 ⇔ 可计算的面积类型对照表
2.1 地理坐标系（GCS）------ 经纬度
| 坐标系名称 | EPSG | 椭球 | ✅ 可算的面积 | ❌ 不能算的面积 | 备注 |
|---|---|---|---|---|---|
| CGCS2000 | 4490 | CGCS2000 / GRS80 | ✅ 椭球面积 ✅ 测地面积 | ❌ 投影面积（需先投影） | 中国法定基准，推荐用于国土统计 |
| WGS84 | 4326 | WGS84 | ✅ 椭球面积 ✅ 测地面积 | ❌ 投影面积 | GPS / 全球通用 |
| Beijing 1954 | 4214 | Krasovsky | ✅ 椭球面积（旧椭球） | ❌ 投影面积 | 已淘汰，不建议新数据使用 |
| Xian 1980 | 4610 | IAG75 | ✅ 椭球面积（旧椭球） | ❌ 投影面积 | 已淘汰，不建议新数据使用 |
⚠️ 在 GCS 下直接算平面面积（°²）是错误的，仅供程序内部提示，不能用于正式成果。
2.2 投影坐标系（PCS）------ 平面米单位
| 投影类型 | 典型 EPSG | 基于的 GCS | ✅ 可算的面积 | ❌ 不能算的面积 | 面积变形程度 | 推荐用途 |
|---|---|---|---|---|---|---|
| CGCS2000 高斯--克吕格 3°带 | 4523~4554（如 4547） | EPSG:4490 | ✅ 投影面积 | ❌ 椭球面积（需反投影） | 极小（<0.001%） | ⭐国土/工程首选 |
| CGCS2000 高斯--克吕格 6°带 | 4491~4500 | EPSG:4490 | ✅ 投影面积 | ❌ 椭球面积 | 小（<0.01%） | 中小比例尺制图 |
| WGS84 / UTM | 326xx / 327xx | EPSG:4326 | ✅ 投影面积 | ❌ 椭球面积 | 小（中央经线0.04%） | 国际项目 |
| Web Mercator | 3857 | WGS84 | ✅ 投影面积（数值偏大） | ❌ 用作统计面积 | 极大（高纬严重） | ❌ 仅限Web显示 |
✅ 投影面积 ≈ 椭球面积，前提是投影本身面积变形足够小（如 GK 3°带）。
三、面积计算合法性判定矩阵
| 数据当前坐标系 | 你想算的面积类型 | 是否需要额外操作 | 结果是否合法 |
|---|---|---|---|
| GCS（4490/4326） | 椭球/测地面积 | ❌ 不需要 | ✅ 合法，推荐 |
| GCS（4490/4326） | 投影面积 | ✅ 需要先投影到 PCS | ✅ 合法 |
| GCS（4490/4326） | °² 面积 | ❌ 不需要 | ❌ 非法，无物理意义 |
| PCS（GK/UTM） | 投影面积 | ❌ 不需要 | ✅ 合法 |
| PCS（GK/UTM） | 椭球/测地面积 | ✅ 需要反投影回 GCS | ✅ 合法 |
| PCS（3857） | 投影面积 | ❌ 不需要 | ⚠️ 数值偏大，不推荐统计 |
四、实战最佳实践（以 CGCS2000 为例）
✅ 推荐方案 A ------ 官方统计（椭球面积）
数据存储格式：EPSG:4490（CGCS2000 地理坐标系）
面积计算方法：椭球面积 / 测地面积
软件操作：
ArcGIS → Calculate Geometry → Area (Ellipsoidal, CGCS2000)
QGIS → $area（图层设为 GCS + 勾选 Ellipsoidal）
Python → geographiclib / pyproj.Geod.geometry_area_perimeter
结果：✅ 最接近真实地表面积，不受投影影响✅ 推荐方案 B ------ 工程制图（投影面积）
数据存储格式：EPSG:45xx（CGCS2000 高斯--克吕格 3°带）
面积计算方法：投影面积（平面鞋带公式）
软件操作：
ArcGIS → Calculate Geometry → Area (Planar)
QGIS → $area（图层在 PCS 下）
Python → shapely.Polygon.area
结果：✅ 数值接近椭球面积，可用于工程红线/地块❌ 禁止做法
❌ 数据在 EPSG:4490 下直接算平面面积（°²）
❌ 用 EPSG:3857 的投影面积做统计
❌ 混用不同基准（如 WGS84 ↔ CGCS2000）不做变换五、快速选择流程图
你的数据是经纬度还是米？
│
├── 经纬度（GCS） ──→ 直接算椭球/测地面积 ✅（推荐）
│ │
│ └── 也可投影到 GK 3°带 → 算投影面积（作为近似）
│
└── 米（PCS） ──→ 直接算投影面积 ✅
│
├── 如果是 GK/UTM → 数值可靠
└── 如果是 Web Mercator → ❌ 不可靠，需反投影六、国内常用 EPSG 编码速查
| 用途 | EPSG | 名称 | 说明 |
|---|---|---|---|
| 存储/分析（推荐） | 4490 | CGCS2000 地理坐标系 | 算椭球面积用这个 |
| GPS/全球交换 | 4326 | WGS84 地理坐标系 | 与 4490 几乎一致 |
| 投影（华北 3°带） | 4547 | CGCS2000 / 3° GK zone 35 | 中央经线 105°E |
| 投影（华东 3°带） | 4549 | CGCS2000 / 3° GK zone 39 | 中央经线 117°E |
| 投影（华南 3°带） | 4548 | CGCS2000 / 3° GK zone 37 | 中央经线 111°E |
| Web 显示 | 3857 | Web Mercator | 仅用于地图瓦片 |
核心记住三点：
- 要算"真实面积" → 用 GCS + 椭球面积
- 要算"工程面积" → 用 GK 3°带 + 投影面积
- 永远别用 °² 和 Web Mercator 做统计