Bagaimana cara menambahkan atribut Fitur khusus ke Shapefile menggunakan Python?

Saya mencari cara untuk mengambil Shapefile yang ada yang memiliki serangkaian Fitur dari 200 negara. Setiap Fitur negara memiliki atribut "NAME." Tujuan saya adalah membuat skrip Python yang menambahkan atribut tambahan sembarang (untuk saat ini), katakan, "POPULASI".

Tentu saja saya sudah menginstal modul OSGeo dan GeoDjango. Saya sejauh:

 from osgeo import ogr

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

Apakah saya kehilangan fungsi OGR yang akan memungkinkan saya untuk memasukkan bidang atribut Fitur ke dalam Shapefile yang ada?

Saya yakin sampel Assemble TIGER Polygons memiliki apa yang Anda cari:

# 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 )

Apakah mungkin untuk menambahkan bidang ke shapefile yang ada menggunakan 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)