I'm trying to use those parameters and I receive a difference of ~0.5m
ALON and BSHM APN Points:
coords = [
(dms2rad([31,42,28.5694]), dms2rad([34,36,23.92196]), 55.671),
(dms2rad([32,46,44.34472]), dms2rad([35,1,22.74061]), 225.046),
]
Transform:
def lla_to_ig05(lon, lat, alt):
import pyproj
lla = pyproj.Proj(proj='latlong', ellps='WGS84', datum='WGS84')
ig05 = pyproj.Proj("+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.33077,-1.852690,1.669690,5.424800 +units=m +no_defs")
x, y, h = pyproj.transform(lla, ig05, lon, lat, alt, radians=True)
return x, y, h
for lat, lon, alt in coords:
print(pd.DataFrame(lla_to_ig05(lon, lat, alt))) # pd for pandas
print("")
Output:
0
0 162790.443392
1 624087.953905
2 55.671000
0
0 202456.903152
1 742717.095893
2 225.046000
Real numbers:
0
0 162789.564
1 624088.355
2 55.671000
0
0 202456.453
1 742718.779
2 225.046000
Any hint will be very helpful!