library(sf)
## Linking to GEOS 3.8.1, GDAL 3.1.3, PROJ 7.1.1
# Sample point
pnt = st_point(c(34.810741, 31.897418))
pnt = st_sfc(pnt, crs = 4326)
# Transform 1 - Python
string_python = "+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"
# Transform 2 - Harel
string_harel = "+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.0024,-17.1032,-17.8444,-0.33009,-1.85269,1.66969,5.4248 +units=m +no_defs"
# Compare
coords_python = pnt %>% st_transform(string_python) %>% st_coordinates
coords_harel = pnt %>% st_transform(string_harel) %>% st_coordinates
as.character(coords_python - coords_harel)
## [1] "-0.00902634940575808" "0.012025561532937"