My first post in the thread I linked to above adds a bounding box check to Jerry's original query. x_delta_max and y_delta_max are half the respective side lengths of the bounding box, so if your distance unit is feet and you want a 500 foot bounding box you would use 250 for each of x_delta_max and y_delta_max.
