rasterio项目中的WarpedVRT与merge功能兼容性问题解析

在rasterio 1.4.0版本中，用户发现了一个重要的功能兼容性问题：rasterio.merge.merge函数无法正确处理rasterio.vrt.WarpedVRT对象。这个问题影响了TorchGeo等依赖rasterio进行地理空间数据处理的项目。

问题背景

地理空间数据处理中，经常需要处理以下几种情况：

1. 合并相邻的栅格文件
2. 将不同坐标参考系统(CRS)的数据重新投影到统一CRS
3. 对超出原始数据范围的区域进行采样

在rasterio 1.3.x版本中，这些功能可以协同工作。用户可以通过WarpedVRT实现内存中的重投影，使用merge函数合并数据，并通过设置boundless=True实现扩展边界读取。

问题表现

升级到rasterio 1.4.0后，当尝试合并经过WarpedVRT处理的数据时，系统会抛出ValueError: WarpedVRT does not permit boundless reads错误。这个错误甚至在没有显式请求扩展边界读取的情况下也会出现。

技术分析

问题的根源在于rasterio 1.4.0中merge函数的内部实现发生了变化。新版本在合并操作中默认使用了扩展边界读取，而WarpedVRT从6年前开始就禁止了扩展边界读取，因为这种操作存在大量潜在错误。

这种变化实际上是一个回归问题(regression)，因为在1.3.x版本中，这种使用方式是可行的。开发团队承认他们没有预见到用户会以这种方式使用WarpedVRT和merge的组合。

解决方案

rasterio团队在1.4.3版本中通过修改merge函数的实现解决了这个问题。新版本不再使用扩展边界读取，因此可以与WarpedVRT正常配合工作。

对于用户而言，如果必须使用1.4.0-1.4.2版本，可以考虑以下替代方案：

1. 放弃内存中的合并操作，改为先单独处理每个文件，最后合并处理结果
2. 使用GDAL的warp工具进行重投影和合并
3. 评估是否真的需要扩展边界读取功能

最佳实践建议

对于需要同时使用重投影、合并和扩展边界读取的场景，建议：

1. 升级到rasterio 1.4.3或更高版本
2. 如果必须使用旧版本，可以尝试修改代码，显式设置boundless=False
3. 考虑将数据处理流程拆分为多个步骤，降低各步骤间的耦合度

这个问题提醒我们，在升级地理空间数据处理库时，需要特别注意功能兼容性变化，特别是那些没有明确文档说明的行为变化。