I'm interested in automating spatial interpolation of groundwater level measurements. I'm wondering if anyone knows a method for automatically taking into account things like dry wells, where I can't get a depth to water measurement, but where I do know that the depth to water is greater than the depth of the screened interval of the well.
The interpolation would be constrained at that location so that the water table would be below the bottom of the well. In addition to saying things like "the value at $(x_1,y_1) = z_1$" I could also specify "the value at $(x_2, y_2) < z_2$" or "at $(x_3, y_3) > z_3$".
I'd like to know if anyone knows any interpolation algorithms that can do this, and if there are any tools out there to implement them.