import ogr, os # Input data pointCoord = -124.4577,48.0135 fieldName = 'test' fieldType = ogr.OFTString fieldValue = 'test' outSHPfn = 'test.shp' # Create the output shapefile shpDriver = ogr.GetDriverByName("ESRI Shapefile") if os.path.exists(outSHPfn): shpDriver.DeleteDataSource(outSHPfn) outDataSource = shpDriver.CreateDataSource(outSHPfn) outLayer = outDataSource.CreateLayer(outSHPfn, geom_type=ogr.wkbPoint ) #create point geometry point = ogr.Geometry(ogr.wkbPoint) point.AddPoint(pointCoord[0],pointCoord[1]) # create a field idField = ogr.FieldDefn(fieldName, fieldType) outLayer.CreateField(idField) # Create the feature and set values featureDefn = outLayer.GetLayerDefn() outFeature = ogr.Feature(featureDefn) outFeature.SetGeometry(point) outFeature.SetField(fieldName, fieldValue) outLayer.CreateFeature(outFeature)
这个方法中创建一个新的形状文件,添加一个指向它,并添加一个值的属性栏吧。