geodata.cutout.calc_shp_area# geodata.cutout.calc_shp_area(shp, shp_projection='+proj=latlon')# calculate area in km^2 of the shapes for each shp object