python-点EPSG 4326(WSG 84)到(EPSG 28992)的PROJ pyproj转换

给定一个(lon,lat)点(5.068913,52.067567)我想使用pyproj从EPSG 4326转换为EPSG 28992.

pyproj中的Proj和transform函数似乎都适合这样的任务：

> https://pyproj4.github.io/pyproj/dev/api/proj.html?highlight=proj#pyproj-proj
> https://pyproj4.github.io/pyproj/dev/api/transformer.html?highlight=transform#pyproj-transform

当我使用Proj函数时,与使用transform会得到不同的结果,为什么呢？

例如

from shapely.geometry import Point
from pyproj import Proj, transform

from matplotlib import pyplot as plt

x1, y1 = 5.068913, 52.067567

in_proj = Proj(init='epsg:4326') 
out_proj = Proj(init='epsg:28992')

point1 = Point(out_proj(x1, y1))
point2 = Point(transform(in_proj, out_proj, x1 ,y1))

print(point1 == point2)

fig, ax = plt.subplots()
x, y = point1.xy
ax.plot(x, y, 'ro')
x, y = point2.xy
ax.plot(x, y, 'ro')

解决方法:

之间的转换

EPSG:4326

和

EPSG:28992

在这种情况下,不仅需要使用Oblique_Stereographic投影来投影输入坐标.这可以通过Proj函数来完成…

…但还执行了基准转换,WGS84到Amersfoort基准.这只能通过转换功能来实现.请注意,变换会执行Proj所做的所有事情以及基准变换.不需要两步转换.

举例说明,我们为EPSG：28992提供了以下WKT,

PROJCS["Amersfoort / RD New",
    GEOGCS["Amersfoort",
        DATUM["Amersfoort",
            SPHEROID["Bessel 1841",6377397.155,299.1528128,
                AUTHORITY["EPSG","7004"]],
            TOWGS84[565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725],
            AUTHORITY["EPSG","6289"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4289"]],
    PROJECTION["Oblique_Stereographic"],
    PARAMETER["latitude_of_origin",52.15616055555555],
    PARAMETER["central_meridian",5.38763888888889],
    PARAMETER["scale_factor",0.9999079],
    PARAMETER["false_easting",155000],
    PARAMETER["false_northing",463000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    AUTHORITY["EPSG","28992"]]

字符串

TOWGS84[565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725]

…意味着需要执行原点转换…这是无法通过Projfunction完成的,而只能通过transform完成.

结论：

假设您正在使用的pyproj安装配置了正确的TOWGS84字符串(并非总是如此),那么转换结果应该被认为是正确的.

标签：coordinate-systems,gis,pyproj,python
来源： https://codeday.me/bug/20191210/2104334.html