import momepy
import geopandas as gpd
import osmnx as ox

#Get building geometries from OSM
BN = ox.geometries_from_bbox(40.550055, 40.436808, -88.917847, -89.032860,
tags={'building':True})

#The data cannot be in a geographic coordinate system for momepy to run correctly
BN_projected = ox.projection.project_gdf(BN)

#Explode muti-part geometries into multiple single geometries.
BN_exploded = BN_projected.explode()

#Strip out polygons from exploded gdf
BN_exploded_poly = BN_exploded.loc[BN_exploded.geometry.type=='Polygon',:]

#Check to make sure all non-polygon geometries have been removed
BN_exploded_poly.geometry

#Reset index - this is a necessary step due to a bug in geopandas
BN_exploded_poly = BN_exploded_poly.reset_index()

print(type(BN_exploded_poly))

#Run building data preprocessing
BN_buildings = momepy.preprocess(BN_exploded_poly, size=30, compactness=0.2, islands=True)