Làm cách nào để thêm thuộc tính Tính năng tùy chỉnh vào Shapefile bằng Python?

Tôi đang tìm cách lấy một Shapefile hiện có có Bộ tính năng gồm 200 quốc gia. Mỗi tính năng quốc gia có một thuộc tính của "NAME." Mục tiêu của tôi là tạo ra một tập lệnh Python có thêm một thuộc tính bổ sung tùy ý (bây giờ), nói, "DÂN SỐ".

Tất nhiên tôi đã cài đặt các mô-đun OSGeo và GeoDjango. Tôi xa như:

 from osgeo import ogr

    infile = ogr.Open('sample.shp', 1) #'sample.shp' is a pre-existing ESRI shapefile described above
    inlyr = ogr.GetLayerByIndex(0)

Tôi có thiếu chức năng OGR cho phép tôi chèn các trường thuộc tính Feature vào Shapefile hiện tại không?

Tôi tin rằng mẫu Đa giác Assemble có những gì bạn đang tìm kiếm:

# Open the datasource to operate on.

ds = ogr.Open( infile, update = 0 )

poly_layer = ds.GetLayerByName( 'Polygon' )

#############################################################################
#   Create output file for the composed polygons.

nad83 = osr.SpatialReference()
nad83.SetFromUserInput('NAD83')

shp_driver = ogr.GetDriverByName( 'ESRI Shapefile' )
shp_driver.DeleteDataSource( outfile )

shp_ds = shp_driver.CreateDataSource( outfile )

shp_layer = shp_ds.CreateLayer( 'out', geom_type = ogr.wkbPolygon,
                                srs = nad83 )

src_defn = poly_layer.GetLayerDefn()
poly_field_count = src_defn.GetFieldCount()

for fld_index in range(poly_field_count):
    src_fd = src_defn.GetFieldDefn( fld_index )

    fd = ogr.FieldDefn( src_fd.GetName(), src_fd.GetType() )
    fd.SetWidth( src_fd.GetWidth() )
    fd.SetPrecision( src_fd.GetPrecision() )
    shp_layer.CreateField( fd )

Có thể thêm một trường vào một shapefile hiện có bằng Python OGR ..

from osgeo import ogr
driver = ogr.GetDriverByName('ESRI Shapefile')
dataSource = driver.Open(“c:/test/Test2.shp”, 1) #1 is read/write

#define floating point field named DistFld and 16-character string field named Name:
fldDef = ogr.FieldDefn('DistFld', ogr.OFTReal)
fldDef2 = ogr.FieldDefn('Name', ogr.OFTString)
fldDef2.SetWidth(16) #16 char string width

#get layer and add the 2 fields:
layer = dataSource.GetLayer()
layer.CreateField(fldDef)
layer.CreateField(fldDef2)