Advertisement
Guest User

Untitled

a guest
Jul 16th, 2019
101
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.91 KB | None | 0 0
  1. src_ds = gdal.Open('input.tif')
  2. if src_ds is None:
  3. print ('Unable to open %s' % src_filename)
  4. sys.exit(1)
  5. try:
  6. srcband = src_ds.GetRasterBand(1)
  7. except RuntimeError:
  8. print ('Band 1 not found')
  9. sys.exit(1)
  10.  
  11. dst_layername = "output/path/filename"
  12. drv = ogr.GetDriverByName("ESRI Shapefile")
  13. if os.path.exists(dst_layername + ".shp"):
  14. drv.DeleteDataSource(dst_layername + ".shp")
  15.  
  16. # create spatial reference for the new vector based on the raster
  17. srs = osr.SpatialReference()
  18. srs.ImportFromWkt(src_ds.GetProjection())
  19. dst_layer = dst_ds.CreateLayer(dst_layername, srs = srs)
  20.  
  21. # create polygon layer
  22. dst_ds = drv.CreateDataSource(dst_layername + ".shp")
  23. gdal.Polygonize(srcband, srcband, dst_layer, -1, [], callback=None)
  24.  
  25. ### fix geometries
  26. processing.run("native:fixgeometries",
  27. {'INPUT':"output/path/filename.shp",
  28. 'OUTPUT':"output_fixed.shp"})
  29. # > produces a file, but shows nothing
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement