使用gdal2tiles.py和Leaflet，如何预处理旋转地理参考影像的瓦片？

我来帮你捋捋这个旋转影像瓦片预处理的问题，核心矛盾就是Leaflet里的地图瓦片是对齐地图坐标系的矩形，而你的原始影像瓦片是带旋转的，直接用load事件拿地图瓦片没法对应到原始影像空间对吧？给你几个可行的方向：

你可以在使用gdal2tiles.py切割瓦片的同时，给每个地图瓦片额外生成一份元数据文件（比如同名的.json），里面记录这个地图瓦片对应的原始影像坐标系下的覆盖范围，以及它对应的原始影像瓦片索引。

• 操作方式：要么修改gdal2tiles.py的代码，在生成每个瓦片时顺便输出元数据；要么切完瓦片后用GDAL工具批量处理，给每个瓦片计算并写入对应的影像空间信息。
• 后续在Leaflet加载瓦片前，先请求对应的元数据文件，拿到信息后就能精准定位到原始影像的对应区域，再做预处理就没问题了。

比起找现成插件，自定义L.GridLayer反而更灵活，能完全适配你的旋转场景：

• 重写createTile方法：当Leaflet需要创建某个地图瓦片时，先把这个瓦片的四个角坐标，通过原始影像的地理变换参数（比如GDAL的GeoTransform）转换成原始影像的像素坐标。
• 根据像素坐标确定需要哪些原始影像瓦片来覆盖这个地图瓦片区域（毕竟旋转后一个地图瓦片可能对应多个原始影像瓦片）。
• 加载这些原始影像瓦片，先做你需要的预处理，再裁剪拼接成符合地图坐标系的矩形瓦片，最后返回给Leaflet渲染。

你说仅靠逆变换不够，确实是因为旋转后地图瓦片的矩形区域，在原始影像空间是个平行四边形，单个地图瓦片可能覆盖多个原始影像瓦片。这时候可以：

• 用影像的旋转矩阵，把地图瓦片的四个顶点都转换成原始影像的像素坐标。
• 根据这些像素坐标，计算出覆盖这个平行四边形的所有原始影像瓦片（比如用像素坐标除以256取整，再扩展边界）。
• 把这些原始瓦片加载后，先预处理，再拼接成地图需要的瓦片形状。

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