Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from osgeo import gdal,ogr
- import sys
- import os
- gdal.UseExceptions()
- os.chdir(path)
- src_ds = gdal.Open("IND_ppp_2015_v2.tif")
- if src_ds is None:
- print ("Unable to open Worldpop data")
- sys.exit(1)
- try:
- srcband = src_ds.GetRasterBand(1)
- drv = ogr.GetDriverByName("ESRI Shapefile")
- dst_ds = drv.CreateDataSource(dst_layername + ".shp")
- dst_layer = dst_ds.CreateLayer(dst_layername,srs = None)
- gdal.Polygonize(srcband, None,dst_layer, -1,[],callback = None)
Add Comment
Please, Sign In to add comment