Read and Write Spatialite with Python

read_and_write_spatialite.py

from osgeo import ogr, osr

def read_splt(db):
    drvr = ogr.GetDriverByName('SQLite')
    ds = drvr.Open(db, 0)
    layer=ds.GetLayer('pts')
    for feature in layer:
        print(feature)
        break

    ds.Destroy()

def write_splt(db):
    drvr = ogr.GetDriverByName('SQLite')
    ds = drvr.CreateDataSource(db, options=['SPATIALITE=yes'])
    ds = drvr.Open(db, 1)
    
    srs = osr.SpatialReference()
    srs.ImportFromEPSG(27700)
    
    lr = ds.CreateLayer("pts", srs, geom_type=ogr.wkbPoint)
    # Add an ID field
    lr.CreateField(ogr.FieldDefn("id", ogr.OFTInteger))

    # Create the feature and set values
    featureDefn = lr.GetLayerDefn()
    feature = ogr.Feature(featureDefn)
    feature.SetGeometry(ogr.CreateGeometryFromWkt('POINT (200000 300000)'))
    feature.SetField("id", 1)
    lr.CreateFeature(feature)
    del lr
    ds.Destroy()
    
if __name__ == "__main__":
    db = '/temp/mydb.sqlite'
    write_splt(db)
    read_splt(db)