Gesamtfläche von Polygonen im Shapefile mit GDAL berechnen?

Ich habe ein Shapefile in der British National Grid-Projektion: Geometry: 3D Polygon Feature Count: 5378 Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952) Layer SRS WKT: PROJCS["British_National_Grid", GEOGCS["GCS_airy", DATUM["OSGB_1936",