Рассчитать общую площадь полигонов в шейп-файле, используя GDAL?

У меня есть шейп-файл в британской национальной проекции сетки: 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",...