I want to polygonize a raster by using GDAL in python. I found this webpage: Polygonize a Raster Band. But it creates a polygon feature layer. I want every band with the same value to create one feature.
For instance, even if band values are far away from each other but the same value, they should be one feature (maybe multipolygon)
EDIT Here is my polygonize code:
from osgeo import gdal, ogr import sys # this allows GDAL to throw Python Exceptions gdal.UseExceptions() # # get raster datasource # src_ds = gdal.Open( "INPUT.tif" ) if src_ds is None: print 'Unable to open %s' % src_filename sys.exit(1) try: srcband = src_ds.GetRasterBand(3) except RuntimeError, e: # for example, try GetRasterBand(10) print 'Band ( %i ) not found' % band_num print e sys.exit(1) # # create output datasource # dst_layername = "POLYGONIZED_STUFF" 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 )
If there are 500 different band values I should get a layer with 500 features.