我想将此地图中的图像：http://earthquake.usgs.gov/regional/nca/soiltype/map/从线性投影堆叠到传单地图中。源图块处于已知但非标准的缩放级别，并且传单贴图需要墨卡特mercator XYZ图块。原则上，我知道如何做到这一点 - 我有两个地图中的XY坐标更改为lat-lng坐标的功能，我只需要根据源地图中的像素映射目标地图的像素。
遗憾的是，这是非常重要的，因为源像素分布在数百个不同的图像文件中，我试图将它们分成数百个图像。是否有一个软件包使这更简单一点？如果没有用于处理此类数据的库，似乎确实应该......
答案 0 :(得分：1)
Postgis采用RT_ST_Transform方法，引擎盖使用GdalWarp。所以，你至少有这两个选项。如果你使用Postgis，你需要使用raster2pgsql实际注册/导入图像到Postgis，然后在每个上调用RT_ST_Transform然后再将它们转储出去 - 这可以在某种程度上使用plpgsql编写脚本（Postgres的脚本语言）。使用Postgis栅格需要学习曲线，如果您计划进行任何其他图像处理分析，这可能是值得的。如果你不想使用Postgis路线，你也可以写一个shell脚本（或类似的）来自动化gdalwarp。
答案 1 :(得分：1)
对于比gdalwarp（一个优秀的程序）更不正式的方法，你可以查看Leaflet插件Leaflet.imageTransform，它可以在浏览器中动态转换和成像。