from shapely import Polygon
from sectionproperties.pre import Geometry, CompoundGeometry

sq1 = Polygon([(0, 0), (2, 0), (2, 2), (0, 2)])
sq2 = Polygon([(2, 0), (6, 0), (6, 4), (2, 4)])
geom_sq1 = Geometry(geom=sq1)
geom_sq2 = Geometry(geom=sq2)
geom = CompoundGeometry(geoms=[geom_sq1, geom_sq2])
geom.plot_geometry()