Guest User

Untitled

a guest
Sep 20th, 2018
103
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.47 KB | None | 0 0
  1. def transform_rd(point):
  2. ''' Please note this returns an OGRGeometry, not a point '''
  3. src_string = '+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 +no_defs no_defs <>'
  4. src_srs = SpatialReference(src_string)
  5. geom = OGRGeometry(point.wkt, src_srs)
  6. geom.transform(4326)
  7. return geom
Add Comment
Please, Sign In to add comment