Total face area of selected meshes
Calculates the total global face area of all mesh objects selected.
Uses bmesh from object to get modifier deforms, and transforms using the world matrix to take into account scale.
import bpy
import bmesh
context = bpy.context
dg = context.evaluated_depsgraph_get()
bm = bmesh.new()
areas = []
for o in context.selected_objects:
if o.type != 'MESH':
continue
bm.from_object(o, dg)
bm.transform(o.matrix_world)
area = sum(f.calc_area() for f in bm.faces)
print(f"{o.name} faces: {len(bm.faces)} area: {area}")
areas.append(area)
bm.clear()
print(f"Total Area: {sum(areas)}")