파이썬에서 빈 다각형을 만드는 GDAL?

파이썬에서 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에서 속성 테이블을 열면 비어 있습니다.

편집하다:

래스터-> 변환-> 다각형 도구를 사용하여 QGIS에서 변환이 제대로 작동합니다.

다각형 화 될 래스터의 스크린 샷 :

그리고 QGIS 도구의 결과 변환 스크린 샷 :

Windows 7에서 Enthought 배포판을 사용하고 있습니다. GDAL 버전 1.10.0-3

문제는 GDAL과 요리 책 예제를 사용하여 파이썬에서 래스터를 다각형 화 할 수 없다는 것입니다 .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 )