Calculer la superficie totale des polygones dans le fichier de formes à l'aide de GDAL?

J'ai un fichier de formes dans la projection British National Grid: 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",