Create spatial index for efficient rectangle querying to avoid 10000 limit

[kinverarity1]

Obviously the simplest strategy would be dividing up 1 degree squares by four.

But a much faster way would be to generate a most efficient way i.e. Some kind of equal density method limited to rectangles.

This would be very useful, because you could then use the drill date in the State wide cache layer to decide for an efficient strategy of updating the different bulk data download functions. Or the observation well network fields in the unit number search json