...
Code Block | ||||
---|---|---|---|---|
| ||||
# Define the Rijksdriehoek projection system (EPSG 28992) epsg28992 = SpatialReference() epsg28992.ImportFromEPSG(28992) rd2latlon = CoordinateTransformation(epsg28992, epsg4326) latlon2rd = CoordinateTransformation(epsg4326, epsg28992) # Check the transformation for the same point, without using towgs84 correction lonlatz = rd2latlon.TransformPoint(155000.0, 446000.0) print (lonlatz) # (5.387202946158022, 52.00237563479786, 43.6057764403522) |
...