Matlab-错误使用 griddedInterpolant 网格向量未定义与给定值匹配的点网格

IDE： Matlab 2018a

BUG : 

错误使用 griddedInterpolant
网格向量未定义与给定值匹配的点网格。

出错 interp3 (line 146)
            F = griddedInterpolant(X, Y, Z, V,
            method,extrap);

出错 slice (line 134)
    vi = interp3(x,y,z,v,sx,sy,sz,method);

代码：

clear all;
load left_kidney.mat;
load p1.mat 
load p2.mat
load p3.mat
x1 = 0;y1 =0; z1 = 0;
x2 = 0;y2 =0; z2 = 0;
x3 = 0;y3 =0; z3 = 0;
x1=p1(1); y1=p1(2); z1 = p1(3);
x2=p2(1); y2=p2(2); z2 = p2(3);
x3=p3(1); y3=p3(2); z3 = p3(3);
ux=x2-x1;
uy=y2-y1;
uz=z2-z1;
vx=x3-x1; 
vy=y3-y1;
vz=z3-z1; 
u_cross_v = [uy*vz-vy*uz, uz*vx-vz*ux, ux*vy-vx*uy];
% a(x-x1)+b(y-y1)+ c(z-z1)=0；
[m,n,q] = size(left_kidney);
[xs,ys]=meshgrid(0:n,0:m);
zs = -1 * (u_cross_v(1)* (xs - x1) + u_cross_v(2)*(ys-y1) - u_cross_v(3)* z1) ./u_cross_v(3);
[X,Y,Z] = meshgrid(1:m,1:n,1:q);
figure;
slice(X,Y,Z,left_kidney, xs, ys, zs);
view(3);
shading interp;
colormap gray
hold on;

代码功能：在3D空间中任意选择三个点，以选定的三个点为基础构造出剖面，并可视化显示。

出错原因：

left_kidney 维度：117x93x102, 与X, Y, Z :93x117x102维度不一致导致。将[X, Y, Z] = meshgrid(1:m, 1:n, 1:q) 替换成 [X,Y,Z]=meshgrid(1:n, 1:m, 1:q)即可。

参考文献：

【1】已知三点求平面方程、平面法向量和点到平面的距离 https://blog.csdn.net/zhouschina/article/details/8784908

【2】Matlab 沿三维任意方向切割CT图的仿真计算 https://www.cnblogs.com/yangwenbo214/p/9836251.html