Versions Compared

Old Version 5

changes.mady.by.user Kees den Heijer

Saved on

compared with

New Version 6

changes.mady.by.user Kees den Heijer

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Code Block
languagepy
titleConvert coordinates using towgs84 parameter
from osgeo.osr import SpatialReference, CoordinateTransformation

# Define the Rijksdriehoek projection system (EPSG 28992)
epsg28992 = SpatialReference()
epsg28992.ImportFromEPSG(28992)

# correct the towgs84
epsg28992.SetTOWGS84(565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812)

# Define the wgs84 system (EPSG 4326)
epsg4326 = SpatialReference()
epsg4326.ImportFromEPSG(4326)

rd2latlon = CoordinateTransformation(epsg28992, epsg4326)
latlon2rd = CoordinateTransformation(epsg4326, epsg28992)

# Check the transformation for a point close to the centre of the projected grid
lonlatz = rd2latlon.TransformPoint(155000.0, 446000.0)
print (lonlatz) # (5.387203018813555, 52.002375635973344, 43.614926571026444)

...