docs(crs): 实证确定项目 CRS = EPSG:4547 (CGCS2000 3度带 GK CM 114E)

- 方法: 同一测线既有网格电极 lat/lon(EPSG:4326) 又有散点 projectX/Y;
  pyproj 把 lat/lon→EPSG:4547 得 dX=0.20m/dY=0.30m(亚米级), 反向~0.3m;
  Krassovsky(北京54)椭球差 44m 排除 → 椭球 GRS80/CGCS2000(与 WGS84 现差<1m)。
- proj4: +proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80
- 旧"EPSG:32649 解出 111°E"系因 32649 中央经线=111°E; 改 CM=114°E 即完美吻合。
- 影响: DEM/影像 与 dd_voxel 回归 不再被 CRS 阻塞(散点可经 4547→4326→GeoLocalFrame
  与帘面/地图配准)。影像 tfw 仍 EPSG:3857 须重投影; dem.tif CRS 待 GDAL 核对。
- 同步 STATUS §4/§5 + spec §5。

docs/superpowers/STATUS.md

@@ -45,8 +45,12 @@
 
 - 登录:基址 `http://tenant.geomative.cn/pop-api`;头 `geomativeauthorization: Geomative <token>`;`getImageCode` 验证码答案**明文** `data.code`;`verifyCodeCheck→login2` 靠**会话 cookie** 串联;login2 `checkCode` 空;token=`data.accessToken`。
 - **RSA 公钥已取**(Playwright route 注入 setPublicKey hook + 缓存绕过 + 真登录捕获),`resources/rsa_public_key.pem`(RSA-2048,PKCS#1 v1.5)。
-- **CRS ⚠️**:剖面 `projectX/Y` 真实 CRS **不是 EPSG:32649**(实测 111°E vs 数据 114°E/香港);影像 tfw 是 EPSG:3857。**做 DEM/底图配准前必须向客户确认项目 CRS**。
-- **坐标统一**:帘面/地图/数据详情已统一用 **lat/lon→GeoLocalFrame(等距圆柱)局部米**;但 **dd_voxel 输入(散点)只有 projX/projY**,与 lat/lon 无法配准(CRS 未确认)→ **体素暂搁置**(VoxelActor 代码保留,未接 UI)。
+- **CRS ✅ 已实证确定(2026-06-08)**:剖面 `projectX/Y` 的真实 CRS = **EPSG:4547「CGCS2000 / 3-degree Gauss-Kruger CM 114E」**(=`+proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80`)。
+  - **如何确定**:同一条测线既有网格电极 `lat/lon`(EPSG:4326)又有散点 `projectX/Y`;用 pyproj 把 lat/lon→EPSG:4547 得 `dX=+0.20m, dY=+0.30m`(亚米级),反向→经纬 `dlon/dlat≈3e-6°(~0.3m)`。Krassovsky(北京54)椭球差 44m 排除→**椭球=GRS80/CGCS2000(与 WGS84 现今差<1m,M1 可视化无影响)**。
+  - 旧结论"EPSG:32649 解出 111°E"是因 32649 中央经线=111°E;改用 CM=114°E 即完美吻合。
+  - **影像 `image.tfw` 仍是 EPSG:3857**(与项目系不同,DEM/底图须 PROJ 重投影);`dem.tif` 自身 CRS 待 GDAL 读时核对。
+  - **影响**:DEM/影像 与 **dd_voxel 回归 不再被 CRS 阻塞**——散点 projX/Y 可经 EPSG:4547→4326→`GeoLocalFrame` 与帘面/地图配准。
+- **坐标统一**:帘面/地图/数据详情已统一用 **lat/lon→GeoLocalFrame(等距圆柱)局部米**;dd_voxel 输入(散点 projX/projY)现可经 **EPSG:4547→lat/lon→GeoLocalFrame** 配准(CRS 已定);VoxelActor 代码就绪,待接入。
 - 网格 `v`/`z` 为 `[j=y][i=x]`,点序 i 最快;深度 `y` 越大越深→渲染时 **z 取负**;colorBar 网格色阶 alpha=0–255、LVL 色阶 alpha=0–1;等值线级用 **colorBar 真实分段值**(均匀分级会一片蓝)。
 
 ## 5. 正确的视图模型（已实现,**P4 续做务必照此,勿再走偏**）

docs/superpowers/specs/2026-06-07-geopro-desktop-m1-design.md

@@ -157,7 +157,7 @@ geopro/
 数据现实（已核验真实样本）：
 
 - 剖面/网格/异常：带 GIS 投影 `projectX`≈516868=**Easting**、`projectY`≈2494259=**Northing**；另带局部米 `xlist/ylist`（各数据集自原点起算）。
-  - ⚠️ **CRS 待确认(Phase 1 用 PROJ 实测纠正)**:`projectX/Y` 的真实 CRS **不是 EPSG:32649**。PROJ 实测 `(516868,2494259)` 在 EPSG:32649 下解出 lon≈**111.16°E**,而网格自带 lat/lon 是 **114.16°E**(docx 标明为**香港** Volia 数据,香港≈114°E)——真实 CRS 的中央经线在 ~114°E(疑为港式/自定义 TM)。**做底图/影像配准(M1.5)前必须向客户确认项目 CRS**。
+  - ✅ **CRS 已实证确定(2026-06-08)= EPSG:4547「CGCS2000 / 3-degree Gauss-Kruger CM 114E」**(`+proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80`)。方法:同测线既有 `lat/lon`(4326)又有 `projectX/Y`,pyproj 转换互校 `dX=0.20m/dY=0.30m`、反向 `~0.3m`;北京54(Krassovsky)差 44m 排除。旧"32649 解出 111°E"系因 32649 中央经线=111°E,改 CM=114°E 即吻合。**DEM/影像与 dd_voxel 配准不再被 CRS 阻塞**(散点 projX/Y 可经 4547→4326→GeoLocalFrame)。datum WGS84 vs CGCS2000 现差<1m,M1 无影响。
   - **对 M1 core 无影响**:`LocalFrame` 用相对米(减原点,CRS 无关),网格自带 lat/lon;`CrsTransform` 已实现并单测验证 PROJ 机制本身。
 - 影像 `image.tfw`：原点 (12708343, 2577685) = **EPSG:3857（Web 墨卡托）**，与剖面**不同投影**。
 - 网格另带 `elevation[100]` / `lat/lon`（经纬度，EPSG:4326）。