fig, ax = plt.subplots(1, 1, figsize=(10, 7))
ax.scatter(x, y, s=1)
for region in boundaries:
    region.plot_region(ax)