如何将Colmap生成的位姿/坐标转换为GPS坐标？技术求助

Colmap输出的位姿属于局部三维坐标系，和GPS的WGS84大地坐标系没有直接映射关系，必须通过「坐标对齐」建立转换规则，以下是两种可行方案：

前提判断 #

先确认你具备以下条件之一：

• 图像带有EXIF GPS信息（拍摄时记录了经纬度、海拔）
• 能提供至少3个已知精确GPS坐标的场景控制点（比如墙角、固定路牌，需在Colmap模型中找到对应3D坐标）

方案1：利用图像EXIF GPS（最简便） #

如果图像自带GPS，按以下步骤操作：

1. 提取EXIF GPS并转十进制
   从所有图像的EXIF中导出GPS数据，整理成图像文件名, 纬度, 经度, 海拔表格。注意把度分秒格式转成十进制（例：北纬30°12'30" = 30 + 12/60 + 30/3600 = 30.2083°）。
2. 提取Colmap相机中心坐标
   Colmap导出的images.txt中，奇数行是图像元数据，格式为：ID, Qw, Qx, Qy, Qz, tx, ty, tz, 相机ID, 文件名。相机中心坐标计算方式：
   • 把四元数Qw,Qx,Qy,Qz转换为3x3旋转矩阵R
   • 提取平移向量t = [tx, ty, tz]^T
   • 相机中心C = -R^T * t（旋转矩阵转置后乘以-t）
3. 坐标系转换与拟合
   • 先把GPS经纬度转成UTM平面坐标（米级，方便做刚体变换），可用Python的pyproj库实现
   • 用最小二乘法拟合旋转矩阵R和平移向量T，满足：GPS_UTM_i = R * Colmap_Center_i + T
4. 批量转换
   用拟合出的R和T把所有Colmap相机中心转成UTM坐标，再转回WGS84经纬度，得到对应GPS坐标。

方案2：使用控制点对齐（无EXIF时） #

如果图像没有GPS，需手动标记控制点：

1. 采集控制点GPS
   在场景中选3个以上特征明显的固定点（如墙角、井盖），用高精度GPS设备（或手机开高精度模式）记录它们的经纬度、海拔。
2. 在Colmap中标记控制点
   打开Colmap GUI加载模型，通过Model Control的Add Point工具，在至少2张图像中标记同一个控制点，完成后可在points3D.txt中提取该点的Colmap局部3D坐标。
3. 拟合变换关系
   收集至少3组Colmap局部坐标 <-> GPS坐标对应对，重复方案1的「转UTM→拟合刚体变换」步骤。
4. 批量转换
   用拟合的变换规则，将所有Colmap相机中心或3D点转换为GPS坐标。

实用代码片段 #

经纬度与UTM互转（Python） #

from pyproj import Transformer

# WGS84经纬度转UTM（以中国50N带为例，EPSG代码32650）
transformer_to_utm = Transformer.from_crs("EPSG:4326", "EPSG:32650")
utm_x, utm_y = transformer_to_utm.transform(latitude, longitude)

# UTM转WGS84经纬度
transformer_to_wgs = Transformer.from_crs("EPSG:32650", "EPSG:4326")
latitude, longitude = transformer_to_wgs.transform(utm_x, utm_y)

注意事项 #

• Colmap局部坐标系的尺度是任意的，若没有已知尺度的控制点（如两个点的实际距离），转换后的GPS坐标会存在尺度误差，需用至少一组已知距离的控制点校准。
• 仅2个控制点只能做平移+旋转，无法校准尺度；3个及以上才能完成完整的刚体变换（含尺度）。
• GPS采集尽量在开阔无遮挡区域，误差控制在5米内，否则转换结果精度会大幅下降。

内容的提问来源于stack exchange，提问作者Techlover