Oblicz całkowity obszar wielokątów w pliku kształtu za pomocą GDAL?

Mam plik kształtu w projekcji 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",