improve zoom region handling.

going from O(Nzoom) to O(Nzoom**(1/3)) by sweeping the region
array for each individual axis.

The sweeping can be further accelerated with a augumented BST for
each axis, especially for the inner most axis.