GDAL boş poligon oluşturmak python poligonize?

Python Polygonize işlevini kullanarak sorun yaşıyorum. Bunun için yemek kitabı örneğini burada bulabilirsiniz .

Kodumun ilgili kısmı:

sourceRaster = gdal.Open('myraster.tif')
band = sourceRaster.GetRasterBand(1)
bandArray = band.ReadAsArray()
outShapefile = "polygonized"
driver = ogr.GetDriverByName("ESRI Shapefile")
if os.path.exists(outShapefile+".shp"):
    driver.DeleteDataSource(outShapefile+".shp")
outDatasource = driver.CreateDataSource(outShapefile+ ".shp")
outLayer = outDatasource.CreateLayer("polygonized", srs=None)
gdal.Polygonize( band, None, outLayer, -1, [], callback=None )
outDatasource.Destroy()
sourceRaster = None

Grubun ilgili bilgilere sahip olduğunu biliyorum, burada bir snippet bandArray:

array([[ 4.,  4.,  3.,  3.,  3.,  2.,  2.,  2.,  2.,  3.,  3.,  3.,  3.,
         3.,  3.,  3.,  3.,  4.,  4.,  4.,  4.,  4.,  4.,  4.,  4.,  4.,
         4.,  4.,  4.,  4.],

Nitelik tablosunu QGIS'de açtığımda boş:

Düzenle:

Raster -> Dönüştürme -> Poligonize aracını kullanarak dönüştürme QGIS'de iyi çalışır

Poligonize edilecek olan raster ekran görüntüsü:

Ve QGIS aracından elde edilen dönüşümün ekran görüntüsü:

Windows 7, GDAL sürüm 1.10.0-3'te Enthought dağıtımını kullanıyorum

Sorun, GDAL ve yemek kitabı örneğini kullanarak python'da bir raster poligonize edemem, QGIS GUI'de sorun olmadan bu aynı raster poligonize edebilirim

Sorun şu ki, tarama bandını saklamak için bir alan oluşturmuyordum. Gdal_polygonize.py dosyasını inceledikten sonra, bunun burada bulunan işlevi kullanan gdal.Polygonize çağrılırken otomatik olarak yapılmadığını fark ettim .

Alan oluşturmak ve alana bant yazmak için gereken ek adım:

newField = ogr.FieldDefn('MYFLD', ogr.OFTInteger)
outLayer.CreateField(newField)

Daha sonra grubu bu alana 0 diziniyle yazabiliriz:

gdal.Polygonize(band, None, outLayer, 0, [], callback=None )