GDAL在Python中多边形创建空白多边形？

我在python中使用Polygonize函数遇到麻烦。可以在此处找到此食谱的示例。

我的代码的相关部分是：

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

我知道乐队有相关信息，这是一个片段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.],

当我在QGIS中打开属性表时，它是空的：

编辑：

使用Raster-> Conversion-> Polygonize工具，该转换在QGIS中工作正常

栅格化的屏幕截图：

以及来自QGIS工具的最终转换的屏幕截图：

我正在Windows 7 GDAL 1.10.0-3版上使用Enthought发行版

问题是我无法使用GDAL和Cookbook示例在python中对栅格进行多边形处理，我可以在QGIS GUI中对同一栅格进行多边形处理而没有问题

问题是我没有创建一个字段来存储栅格波段。深入研究gdal_polygonize.py文件后，我意识到调用gdal.Polygonize时并不会自动完成此操作，而是使用此处找到的功能。

这是创建字段并向该字段写入乐队所需的额外步骤：

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

然后，我们可以将频段写入索引为0的该字段：

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