def select_area_of_interest(blocks, poly):
    """Select all blocks that are within the boundaries given. Utility function."""
    return blocks[blocks.geometry.map(lambda block: poly.contains(block))]