你几乎把它弄好了.您需要定义坐标的3D网格.创建单个向量不是正确的方法.你当然可以在这里使用interp3.尝试做：

[X,Y,Z] = meshgrid(1:213, 1:100, 1:140);

Vq = interp3(M, X, Y, Z);

请注意,我已经交换了行(100)和列(213)限制,因为第一个参数水平前进,而第二个参数垂直前进.

此外,通过以这种方式使用interp3,我们假设X,Y和Z的限制在1：213,1：100和1：140之间.如果您提供超出这些限制的任何值,您将获得NaN.有几种方法可以避免这种情况：

>在末尾指定样条标志以允许样条外推

>如果要调整矩阵的大小(就像调整图像大小一样),那么目前没有可以通过这种方式调整3D矩阵大小的内置方法.你必须自己写这个.

如果要执行步骤2,可以执行以下操作.

首先,您需要确定每个维度的比例因子.基本上,这是每个维度的输出大小与原始输入大小的比率.

在此之后,您创建一个2D网格,其限制受输入矩阵的原始大小限制,但此网格的大小将是输出矩阵的大小.比例因子在这里很有用,因为这有效地为我们提供了网格中每个值应该插值的内容.我们将创建新的坐标,从1到每个维度的输出大小,以1 / scaleFactor为增量.例如,如果我们想要将矩阵的大小加倍,则这是2的因子.如果我们的X和Y坐标分别从1到3和1到3,原始网格看起来像这样：

X = Y =

1 2 3 1 1 1

1 2 3 2 2 2

1 2 3 3 3 3

要加倍,这只是：

X = Y =

1 1.5 2 2.5 3 1 1 1 1 1

1 1.5 2 2.5 3 1.5 1.5 1.5 1.5 1.5

1 1.5 2 2.5 3 2 2 2 2 2

1 1.5 2 2.5 3 2.5 2.5 2.5 2.5 2.5

1 1.5 2 2.5 3 3 3 3 3 3

请注意,这会创建一个5 x 5的输出网格.为了使其加倍为6 x 6,您可以执行任何操作,但为了简单起见,只需复制最后一行和最后一列,以便：

X = Y =

1 1.5 2 2.5 3 3 1 1 1 1 1 1

1 1.5 2 2.5 3 3 1.5 1.5 1.5 1.5 1.5 1.5

1 1.5 2 2.5 3 3 2 2 2 2 2 2

1 1.5 2 2.5 3 3 2.5 2.5 2.5 2.5 2.5 2.5

1 1.5 2 2.5 3 3 3 3 3 3 3 3

1 1.5 2 2.5 3 3 3 3 3 3 3 3

这定义了用于调整大小的2D列网格.现在,这是在3D中调整大小的问题.我们能做的是在切片之间进行插值.我们可以使用MATLAB中的permute轻松完成此操作,我将在稍后向您展示如何执行此操作.因此,基本算法是这样的：

>确定所需输出矩阵的输出大小

>确定每个维度中的比例因子

>按照上述步骤为每个维度创建插值访问值的2D网格

>对于矩阵中的每个2D切片,使用interp2使用上述2D网格将每个切片的大小调整为输出行和列.

>之后,使用interp1和permute来调整第三维的大小.

不用多说了,这是执行此操作的代码：

%// Specify output size of your matrix here

outputSize = [100 213 140];

%//Figure out size of original matrix

d = size(M);

%//Scaling coefficients

scaleCoeff = outputSize ./ d;

%//Indices of original slices in 3D

z = 1:d(3);

%//Output slice indices in 3D

zi=1:1/scaleCoeff(3):d(3);

%//Create gridded interpolated co-ordinates for 1 slice

[X,Y] = meshgrid(1:1/scaleCoeff(2):d(2), 1:1/scaleCoeff(1):d(1));

%//We simply duplicate the last rows and last columns of the grid if

%//by doing meshgrid, we don't get exactly the output size we want

%//This is due to round off when perform 1/scaleCoeff(2) or

%//1/scaleCoeff(1). We would be off by 1.

if size(X,1) ~= outputSize(1)

X(end+1,:) = X(end,:);

Y(end+1,:) = Y(end,:);

end

if size(X,2) ~= outputSize(2)

X(:,end+1) = X(:,end);

Y(:,end+1) = X(:,end);

end

%//For each slice...

M2D = zeros(outputSize(1), outputSize(2), d(3));

for ind = z

%//Interpolate each slice via interp2

M2D(:,:,ind) = interp2(M(:,:,ind), X, Y);

end

%//Now interpolate in 3D

MFinal = permute(interp1(z,permute(M2D,[3 1 2]),zi),[2 3 1]);

%//If the number of output slices don't match after we interpolate in 3D, we

%//just duplicate the last slice again

if size(MFinal,3) ~= outputSize(3)

MFinal(:,:,end+1) = MFinal(:,:,end);

end

MFinal将是您最终的插值/调整大小的3D矩阵.在3D中插值的关键方法是置换方法.这将做的是,对于z的每个值,我们将生成2D切片值.因此,如果我们在z = 1处有切片而在z = 2处有一个切片,如果我们想要找到切片z = 1.5处的值的2D网格,这将生成使用信息创建这些插值的2D切片.在z = 1和z = 2之间.我们首先调用permute来执行此操作,然后进行另一个permute调用以撤消我们的排列并获取原始尺寸.