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