I have some data on 3D models of soil's moisture content, the model is in the finite-element grid. I have converted the program outputs to a table in four columns x, y, z, Theta where Theta is the moisture content. I want to take cross sections from the 3D domain in some specific positions (section ABCD in the figure). I want to calculate the value of Theta in a 5*5 grid in the cross-section, but the points around each node of the grid are not coplanar with the unknown point. I did this before for 2D domains but the 3D domains seem more complicated for me.
I wonder if someone could help me with a method to calculate the moisture of point P0 depending on some neighboring points P1, P2, ...
I found an answer to my question. The answer is to used Gridded interpolation like the Kriging method or any similar nD interpolation especially the trilinear interpolation.
I have posted an implementation in
pythonon this page at StackOverflow.Hope everyone benefits from it.
Best regards