Skip to content

Instantly share code, notes, and snippets.

@hoogenm

hoogenm/RD_to_wgs.py

Last active October 21, 2018 20:14
Show Gist options
  • Save hoogenm/4f679f66790823d4a5ec2d22ce877615 to your computer and use it in GitHub Desktop.
Save hoogenm/4f679f66790823d4a5ec2d22ce877615 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
RD_to_wgs.py
import pyproj
def RD_to_wgs(x, y):
rd_projection = pyproj.Proj("+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812 +units=m +no_defs")
wgs84_projection = pyproj.Proj(proj='latlong',datum='WGS84')
lat, lon, _ = pyproj.transform(rd_projection, wgs84_projection, x, y, 0.0)
return [lat, lon]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment