gdal 重新投影的时候,注意不要包含极点坐标,会有神奇的bug

建议使用 gdal_translate 提前切除极点坐标.

gdal_translate -projwin 90 60 160 0 -of GTiff input.tif ouput.tif

推荐阅读

https://github.com/syncpoint/terrain-rgb