Advertisement
Guest User

Python code

a guest
Mar 14th, 2023
120
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.61 KB | None | 0 0
  1. from osgeo import ogr, osr
  2.  
  3. driver = ogr.GetDriverByName('GeoJSON')
  4. dataset = driver.Open('geojson.json')
  5. layer = dataset.GetLayer()
  6. feature = layer.GetNextFeature()
  7. geom = feature.GetGeometryRef()
  8.  
  9. source = osr.SpatialReference()
  10. source.ImportFromProj4("+proj=stere +lat_0=90 +lat_ts=60 +lon_0=264 +x_0=0 +y_0=0 +R=6371229 +units=m +no_defs")
  11. target = osr.SpatialReference()
  12. target.ImportFromEPSG(4326)
  13. transform = osr.CreateCoordinateTransformation(source, target)
  14. geom.Transform(transform)
  15.  
  16. result_json = geom.ExportToJson()
  17.  
  18. with open("result_python.json", "w") as outfile:
  19. outfile.write(result_json)
  20.  
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement