Rewrite function subdividing a region to look for empty spaces

[rdpickard]

The find_non_empty_regions_bounded function looks really janky.

There is probably a better way to do it in numpy.