如何将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)