Linear interpolation between two surfaces

[GeoHamster]

Hi, I'm trying to create an interpolated volume from two surfaces. Let me explain exactly what I'm doing.

I am trying to obtain a rough estimate for the temperature in a certain geographical area and at depth. I have the temperature of the rocks at the surface in the form T1=T(x,y,z) where z is the elevation above sea level at the surface and T is the temp (varying between 5°C and 20°C). I also have the depth to the Curie temp, obtained from spectral analysis of aeromagnetic data, as another surface, much deeper, in the form Tc=T(x,y,z), where Tc is a constant temperature (600°C) and z is the depth to this temperature (varying between -20km and -30km above sea level, but never folding over itself).

I want to take every T1(x0,y0,z1) location and interpolate vertically to the corresponding Tc(x0,y0,zc) point assuming a linear geothermal gradient ([Tc-T1]/[zc-z1]).

With this, I should be able to obtain a rough 3D temperature model of the area.

I am a beginner with MATLAB so any workflow using it is truly appreciated!

Thanks!

 

[amirkhanas80]

Hi GeoHamster
I am studying geophysics. I encountered a problem the same as yours 'Linear interpolation between two surfaces'. I need to generate two interfaces as 'isosurfaces', then apply vertical linear interpolation between them. I think it is exactly the same problem as you described in your post. As your post is for one year ago, and also no one answered it, did you find any solution to it?
I deeply appreciate your help and anyone else who can give me some clue. My need is emergency.
Thank you very much.