python - 如何将 epsg:6991 以色列网格转换为 epsg:4326 wgs 84 python

Question

我得到了带有 ["x","y"] 列的数据框。

“x”和“y”是 epsg:6991 以色列网格的坐标

例如我喜欢转换x=[212513], y=[617119]坐标，

至： epsg:4326 wgs 84

x=[35°7'2.734"] , y=[31°38'56.771"]

如何在两种类型的坐标之间进行转换？

Answer 1

import pyproj
from pyproj import CRS
from pyproj import Transformer

crs = CRS.from_epsg(6991)
crs.to_epsg()
crs = CRS.from_proj4("+proj=tmerc +lat_0=31.7343936111111 +lon_0=35.2045169444445 +k=1.0000067 +x_0=219529.584 +y_0=626907.39 +ellps=GRS80 +towgs84=-24.002400,-17.103200,-17.844400,-0.33007,-1.852690,1.669690,5.424800 +units=m +no_defs")
transformer = Transformer.from_crs("EPSG:6991", "EPSG:4326")
transformer.transform(220080, 634451)