GDAL中DIMAP格式投影信息处理问题解析

问题背景

在GDAL 3.4.3之后的版本中，DIMAP格式驱动在处理传感器产品时出现了一个重要变化。DIMAP是一种用于存储卫星遥感数据的格式，特别适用于Pleiades等卫星的传感器产品。这些产品原本处于焦平面坐标系（Sensor坐标系），本不应包含投影参考信息。

问题表现

在GDAL 3.4.3版本中，处理DIMAP格式的传感器产品时不会添加投影参考信息，这与产品的实际特性相符。然而，从GDAL 3.7.1版本开始（包括3.10版本），驱动会自动为这些产品添加投影参考信息，这与产品的物理特性不符。

技术分析

DIMAP格式文件中包含多种地理信息元素：

1. KMZ预览文件：用于帮助用户定位产品
2. RPC模型：用于几何定位
3. 产品足迹：在WGS84坐标系中计算得出

关键问题在于RASTER_GEOMETRY元素的处理。在DIMAP文件中，该元素明确指定了"SENSOR"值，表明数据处于传感器坐标系而非投影坐标系。正确的处理方式应尊重这一设置，不强制添加投影参考。

解决方案

GDAL开发团队已经修复了这个问题。修复方案主要涉及：

1. 正确解析DIMAP文件中的RASTER_GEOMETRY元素
2. 当该元素值为"SENSOR"时，不自动添加投影参考
3. 确保JP2文件中的GML描述与传感器坐标系一致

对用户的影响

这一修复确保了：

1. 传感器产品能够保持其原始坐标系特性
2. 应用软件（如Orfeo Toolbox）能够正确处理这些数据
3. 避免了因错误添加投影参考导致的后续处理问题

最佳实践

用户在处理DIMAP格式的传感器产品时应注意：

1. 确认GDAL版本是否包含此修复
2. 检查产品的RASTER_GEOMETRY设置
3. 对于需要地理参考的应用，应使用产品提供的RPC模型而非假设的投影

这一修复体现了GDAL对遥感数据特性的精确把握，确保了数据处理的准确性。