Interpolation of layers to obtain a layer in the orthogonal plan

Question

Ricardo Duarte 2021-9-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1457739-interpolation-of-layers-to-obtain-a-layer-in-the-orthogonal-plan

评论： Ricardo Duarte 2021-9-23

Dear all,

I have three variables S1, S2 and S3 with the size (1000x1000) all that represent the sound levels in three different layers at three different depths. Now, I want to use those layers to calculate the sound levels in an orthogonal plan (kind of perpendicular cut of those layers), by using the interpolation of those three layers (see the figure below).

What and how can I do that?

Thank you all in advance.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Adam Danz 2021-9-21

In your description, it states that you have three layers at three different depths which makes the data set 3D.

If that's not the case, could you provide an explicit example?

Ricardo Duarte 2021-9-22

Yes, I have three slices (see the file attached).

What I want is to do a cut in the perpendicular axes (cutting plan) and extract the resulting new layer.

Thank you for your help.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Adam Danz 2021-9-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1457739-interpolation-of-layers-to-obtain-a-layer-in-the-orthogonal-plan#answer_792899

编辑：Adam Danz 2021-9-22

在 MATLAB Online 中打开

Load data

data1 = load('S1.mat');
data2 = load('S2.mat');
data3 = load('S3.mat');
depths = [1 5 10]; 

Concatenate layers to 3D array

layers = cat(3, data1.Layer1_1m_, data2.Layer2_5m_, data3.Layer3_10m_);

Select a slice along the x axis. The select the a YZ plane at x=5

xIdx = 5; % x-slice

Interpolate

yzSlice = squeeze(layers(xIdx, :, :));  % yz plane at x=xIdx
newZ = min(depths) : max(depths);       % interpolated z values
sliceInterp = interp1(depths, yzSlice', newZ)';  % interplated plane

Plot results. Circles and lines (-o) are the raw data. Dots are the interpolation.

figure()

plot(depths, yzSlice', 's-')

hold on

set(gca, 'ColorOrderIndex', 1) % reset color order to match

plot(newZ, sliceInterp', '.')