gpt4 book ai didi

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

转载 作者:行者123 更新时间:2023-11-28 16:57:58 24 4
gpt4 key购买 nike

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

pyproj 中的Projtransform 函数似乎都适合这样的任务:

当我使用 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')

points

最佳答案

之间的转换

EPSG:4326

EPSG:28992

不仅需要投影输入坐标,在本例中使用 Oblique_Stereographic 投影。这可以通过 Proj 函数来完成...

...但也执行基准转换,WGS84 到阿默斯福特基准。这只能通过 transform 函数来完成。请注意,transform 执行 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]

...意味着需要执行数据转换...这不能用 Proj 函数完成,只能用 transform 完成。

结论:

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

关于python - 点 EPSG 4326 (WSG 84) 到 (EPSG 28992) 的 PROJ pyproj 转换,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/56528564/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com