Guest User

Untitled

a guest
Jan 18th, 2018
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.46 KB | None | 0 0
  1. from osgeo import gdal,ogr
  2. import sys
  3. import os
  4. gdal.UseExceptions()
  5. os.chdir(path)
  6.  
  7. src_ds = gdal.Open("IND_ppp_2015_v2.tif")
  8.  
  9. if src_ds is None:
  10. print ("Unable to open Worldpop data")
  11. sys.exit(1)
  12.  
  13. try:
  14. srcband = src_ds.GetRasterBand(1)
  15. drv = ogr.GetDriverByName("ESRI Shapefile")
  16. dst_ds = drv.CreateDataSource(dst_layername + ".shp")
  17. dst_layer = dst_ds.CreateLayer(dst_layername,srs = None)
  18.  
  19. gdal.Polygonize(srcband, None,dst_layer, -1,[],callback = None)
Add Comment
Please, Sign In to add comment