GDAL Python 创建轮廓线

Question

我想在 Python 中从 SRTM 图像生成轮廓线。它似乎可以计算，但如果我想添加轮廓线，则不会显示任何内容，并且属性表也为空。请看一下我的代码：

    from osgeo import gdal, gdal_array
    from osgeo.gdalconst import *
    from numpy import *
    from osgeo import ogr

    #Read in SRTM data
    indataset1 = gdal.Open( src_filename_1, GA_ReadOnly)
    in1 = indataset1.GetRasterBand(1)
    
    #Generate layer to save Contourlines in
    ogr_ds = ogr.GetDriverByName("ESRI Shapefile").CreateDataSource(dst_filename)
    contour_shp = ogr_ds.CreateLayer('contour')
    
    field_defn = ogr.FieldDefn("ID", ogr.OFTInteger)
    contour_shp.CreateField(field_defn)
    field_defn = ogr.FieldDefn("elev", ogr.OFTReal)
    contour_shp.CreateField(field_defn)
    
    #Generate Contourlines
    gdal.ContourGenerate(in1, 100, 0, [], 0, 0, contour_shp, 0, 1)
    ogr_ds.Destroy()

字段 ID 和字段高程似乎是空的，但轮廓形状文件相当大~100MB。

知道可能出了什么问题吗？

更新：我明白了！我忘记关闭数据源：ogr_ds.Destroy()

Answer 1

不要使用该Destroy()方法，如 GDAL/OGR Python 陷阱中所述。

要保存并关闭数据集，请取消引用该变量，并可以选择将其删除。

我通常在最后使用它来保存/关闭 GDAL 或 OGR 数据集：

ogr_ds = None
del ogr_ds