COLMAP模型转换中纯缩放矩阵导致相机姿态异常的技术解析

问题背景

在使用COLMAP进行三维重建时，用户可能会遇到需要调整模型比例的情况。特别是当使用带有已知基线的立体相机系统时，重建结果的比例往往需要与真实世界尺度对齐。COLMAP提供了model_transformer工具来实现这一目的，但在使用纯缩放矩阵时却出现了意外的相机姿态变换问题。

现象描述

当用户尝试使用一个纯缩放矩阵（仅包含对角线缩放值，无旋转分量）通过model_transformer工具对稀疏重建模型进行缩放时，发现相机姿态发生了不合理的旋转和扭曲。具体表现为：

1. 原始重建模型中的相机姿态与点云对齐良好
2. 应用纯缩放变换后，相机姿态出现明显倾斜
3. 点云本身缩放正确，但与相机姿态不再匹配

技术分析

根本原因

经过深入分析，发现问题出在COLMAP内部对变换矩阵的解析方式上。model_transformer工具期望的输入文件格式实际上是7个参数的表示法（1个缩放因子+4个旋转四元数+3个平移向量），而非4×4的变换矩阵。

当用户提供4×4的纯缩放矩阵时，COLMAP的Sim3d::FromFile方法无法正确解析出旋转分量，导致生成的相似变换中的旋转四元数被错误地设置为零向量（0,0,0,0），而非预期的单位四元数（1,0,0,0）。

内部机制

COLMAP使用相似变换（Sim3）来表示模型的变换，包含三个分量：

1. 缩放因子（scale）
2. 旋转（表示为四元数）
3. 平移（3D向量）

当使用4×4矩阵作为输入时，系统尝试从中提取这些分量。对于纯缩放矩阵，正确的处理应该是：

• 缩放因子取对角线元素
• 旋转设置为单位四元数（无旋转）
• 平移设置为零向量

然而实际实现中，旋转分量的提取出现了错误，导致最终得到的变换包含了无效的旋转信息。

解决方案

正确使用方法

要正确使用model_transformer工具进行纯缩放变换，应该：

1. 创建一个文本文件，包含7个参数：

   scale qw qx qy qz tx ty tz
   

2. 对于纯缩放，设置为：

   缩放因子 1 0 0 0 0 0 0
   

临时修复方案

如果确实需要使用4×4矩阵作为输入，可以修改COLMAP源代码，在应用变换时强制将旋转分量设置为单位四元数：

Sim3d temp_transform = new_from_old_world;
temp_transform.rotation.setIdentity();
frame.SetRigFromWorld(
    TransformCameraWorld(temp_transform, frame.RigFromWorld()));

最佳实践建议

1. 明确输入格式：在使用COLMAP工具时，务必查阅文档确认输入参数的准确格式要求
2. 验证变换结果：应用任何变换后，都应该检查相机姿态与点云的相对关系是否合理
3. 比例校正策略：对于立体相机系统，建议在初始重建时就通过相机标定参数确保正确比例，而非事后缩放
4. 开发注意事项：当实现相似变换相关功能时，需要特别注意对纯缩放、纯旋转等特殊情况的处理

总结

COLMAP作为强大的三维重建工具，其model_transformer功能在正确使用时能够有效调整重建模型的比例和姿态。理解工具期望的输入格式和内部处理机制，可以帮助用户避免类似问题。对于开发者而言，这也提醒我们在处理特殊变换情况时需要格外小心，确保数学表示的准确性。

通过本文的分析，希望读者能够更好地理解COLMAP中模型变换的工作原理，并在实际应用中避免类似的陷阱，获得更准确的三维重建结果。