使用gdal_translate转换红外ECW至TIF时遇Read/Write错误求助

问题背景 #

• 处理大型ECW文件裁剪转TIFF时，RGB文件执行以下命令可正常完成转换：

gdal_translate -of GTiff -projwin 101000.0 444000.0 101750.0 443250.0 -co COMPRESS=DEFLATE -co TILED=YES -co predictor=2 2023_PLRR_tm_hr_RGB_o.ecw test.ecw

• 红外文件执行相同命令时出现无明确指向的错误：

gdal_translate -of GTiff -projwin 101000.0 444000.0 101750.0 443250.0 -co COMPRESS=DEFLATE -co TILED=YES -co predictor=2 2022_PLRR_tm_hr_CIR_o.ecw test.ecw

错误输出：

Input file size is 420000, 360000
0ERROR 1: Could not perform Read/Write on file

• 两类文件系统权限一致，仅PDAL可正常读取红外文件，尝试指定CRS后使用gdalwarp仍报错。

排查与解决步骤 #

1. 检查ECW关联文件完整性
   ECW格式可能依赖.eww索引辅助文件，需确保2022_PLRR_tm_hr_CIR_o.eww与对应.ecw文件同名且在同一目录，且文件未损坏。若RGB文件无对应.eww，可尝试删除红外文件的.eww，让GDAL重新生成索引文件。
   同时用gdalinfo验证红外文件有效性：

   gdalinfo 2022_PLRR_tm_hr_CIR_o.ecw
   

   若gdalinfo报错，说明文件本身损坏或GDAL驱动不兼容。

2. 验证GDAL ECW驱动兼容性
   ECW为专有格式，GDAL对不同版本ECW的支持存在差异：

   • 检查GDAL是否启用ECW驱动：

     gdalinfo --formats | grep ECW
     

     若输出为空，需重新编译GDAL并加入ECW驱动支持（需获取ECW SDK授权）。
   • 尝试更新GDAL至最新稳定版，或匹配对应版本的LibECW库。

3. 修复缺失CRS的问题
   -projwin参数依赖正确的坐标系统，先为红外文件赋予CRS再执行转换：

   # 替换EPSG:XXXX为实际坐标系编码
   gdal_edit.py -a_srs EPSG:XXXX 2022_PLRR_tm_hr_CIR_o.ecw
   # 重新执行转换
   gdal_translate -of GTiff -projwin 101000.0 444000.0 101750.0 443250.0 -co COMPRESS=DEFLATE -co TILED=YES -co predictor=2 2022_PLRR_tm_hr_CIR_o.ecw test_cir.tif
   

   若仍报错，先跳过裁剪直接转换小范围测试，排除区域坐标问题：

   gdal_translate -of GTiff -srcwin 0 0 1000 1000 -co COMPRESS=DEFLATE -co TILED=YES -co predictor=2 2022_PLRR_tm_hr_CIR_o.ecw test_small.tif
   

4. 排查系统资源与权限细节

   • 确认输出目录的写入权限，以及磁盘剩余空间是否充足（大型文件转换需临时空间）。
   • 针对大文件添加BIGTIFF参数，避免超过4GB的TIFF文件报错：

     gdal_translate -of GTiff -projwin 101000.0 444000.0 101750.0 443250.0 -co COMPRESS=DEFLATE -co TILED=YES -co predictor=2 -co BIGTIFF=YES 2022_PLRR_tm_hr_CIR_o.ecw test_cir.tif
     

5. 通过PDAL中转转换
   利用PDAL可读取红外文件的特性，先转为TIFF再用GDAL裁剪：

   # 替换EPSG:XXXX为实际坐标系编码
   pdal translate 2022_PLRR_tm_hr_CIR_o.ecw temp_cir.tif --writers.gdal.crs=EPSG:XXXX
   # 执行裁剪
   gdal_translate -of GTiff -projwin 101000.0 444000.0 101750.0 443250.0 -co COMPRESS=DEFLATE -co TILED=YES -co predictor=2 temp_cir.tif test_cir.tif
   

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