如何使用gdalwarp提取尺寸一致的城市缩略图？

命令构建逻辑 #

给定参数：巴黎中心坐标48.8566°N, 2.3522°E，范围100km×100km，GSD=1km，对应输出100×100像素。需先将距离转换为地理坐标边界，再通过gdalwarp提取目标区域。

最终命令 #

gdalwarp -te $(echo "2.3522 - 50/111.32" | bc) $(echo "48.8566 - 50/111.13" | bc) $(echo "2.3522 + 50/111.32" | bc) $(echo "48.8566 + 50/111.13" | bc) -tr 0.00898315 0.00898315 -ts 100 100 -r bilinear input_europe.tif paris_thumbnail.tif

参数详解 #

• -te：定义输出影像的地理范围（xmin ymin xmax ymax），通过计算得到巴黎中心向东西南北各延伸50km的经纬度边界（WGS84坐标系下，1°经度约为111.32km，1°纬度约为111.13km）
• -tr：指定输出影像的分辨率（GSD），1km对应约0.00898315°
• -ts：强制输出影像尺寸为100×100像素，确保与GSD匹配
• -r bilinear：采用双线性重采样算法，平衡缩略图的清晰度和平滑度
• input_europe.tif：输入的欧洲地理参考TIFF文件路径
• paris_thumbnail.tif：输出的巴黎缩略图文件路径

额外说明 #

• 若原影像采用非WGS84坐标系（如UTM），需先将巴黎中心坐标转换为对应坐标系的坐标，再计算边界偏移，同时可添加-t_srs参数指定输出坐标系
• 如需更精准的坐标计算，可使用gdaltransform工具完成距离到地理坐标的转换

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