import geopandas
import matplotlib.pyplot as plt
from geodatasets import get_path
path_to_data = get_path("nybb")
gdf = geopandas.read_file(path_to_data)
ax = gdf.plot(column="Shape_Area", cmap="viridis", legend=True, figsize=(8, 8))
ax.set_axis_off()
plt.tight_layout()
plt.show()