python - 使用 tilestache/modestmaps 将 wgs84 中的 lon/lat 重新投影到目标投影？

Question

tilestache 的文档有点稀疏。有没有办法将 wgs84 中的点（经度/纬度）转换为目标投影？

在哪里：

target_srid = 3095
target_proj4_str = '+proj=utm +zone=54 +ellps=bessel +towgs84=-146.414,507.337,680.507,0,0,0,0 +units=m +no_defs'

我注意到 ModestMaps 中有一个转换对象：

ModestMaps.Geo.Transform()

但我不清楚如何使用它。

Answer 1

你最好的选择是“Goodies”下的 Proj4Projection 类：

http://tilestache.org/doc/TileStache.Goodies.Proj4Projection.html

我目前只在 TileStache 的核心中包含 WGS84 和 Spherical Mercator，但上面提供的模块应该对您有所帮助。请注意，除了投影之外，您还需要计算一个转换，这在 Per Liedman 的出色指南中进行了解释：http ://blog.kartena.se/local-projections-in-a-world-of-spherical-mercator /

Answer 2

Proj4Projection 有一个依赖项 (pyproj)。我不确定您具体要做什么，但是如果您试图使数据与 TileStache 的预期兼容，我只会ST_Transform在实际 SQL 查询中使用将 PostGIS 项目设置为 4326。这不是最干净的方法来做到这一点，但它至少会防止数据冗余。

另见：http: //vis4.net/blog/posts/no-more-mercator-tiles/