Bagaimana cara mengubah geometri Shapefile ke WKB menggunakan OGR?

Saya mendapat geometri dari fitur shapefile, dan saya ingin menyimpan geometri itu di postgis (dalam format WKB seperti itu terjadi ketika kita mengimpor shapefile menggunakan perintah shp2pgsql dan psql). Bagaimana saya mengonversi itu

Untuk mendapatkan Geometri saya telah menggunakan perpustakaan OGR OSGeo misalnya:

feat = layer.GetFeature(0)
geometry = feat.GetGeometryRef()

dan saya punya

<osgeo.ogr.Geometry; proxy of <Swig Object of type 'OGRGeometryShadow *' at 0x0096A2D8> >

jadi bagaimana akan mengubahnya menjadi geometri WKB? Saya menggunakan Python untuk ini.

Anda hampir sampai. Anda hanya perlu memanggil ExportToWkbfungsinya.

import ogr
# Get the driver
driver = ogr.GetDriverByName('ESRI Shapefile')
# Open a shapefile
shapefileName = "D:/temp/myshapefile.shp"
dataset = driver.Open(shapefileName, 0)

layer = dataset.GetLayer()
for index in xrange(layer.GetFeatureCount()):
    feature = layer.GetFeature(index)
    wkb = feature.GetGeometryRef().ExportToWkb()