利用插值函数求三维图像的二分点及到定曲面等距离的曲面.docxVIP

二元函数的分片段线性插值 一函数构造 1?线性三角插值 如上图所示的三角形，已知三角形的三个顶点P1,P2, P3 如上图所示的三角形，已知三角形的三个顶点P1, P2, P3 处的函数值分别为Zl=f(xl,yl), Z2=f(x2,y2), Z3二f(x3,y3)。我们需要求出线性函 数 P(x,y)二ax+by+c,使其满足 P(xl,yl)=Zl, P(x2,y2)=Z2, P(x3,y3)二Z3。 三角形的面积的两倍我们可以用行列式来表示， 西X 1 D= x2 y2 1,我们假设P(x,y) = 1,(x,y)Zj +12(x,y)Z2 +13(x,y)Z3, 么可以 兀3力1 构造基 函数如 丿2(兀,刃=万 y 3 兀1 x2 X 显然，三个基函数分别满足插值条件: 2?相邻节点的二分点，三分点的x, y坐标 所给的数据x,y,z结构均为矩阵，如下面的图形所示，我们用矩形表示x, y的矩阵，矩形的顶点为Pl P2 P3 P4,我们计算四个顶点的二分点。如果我们连 接P2 P3,会得到两个三角形。 PlMlP2 Pl Ml P2 M4 M5 M2 P3 M3 P4 二分点 二分点坐标的mat lab代码实现 function xxl=inter2xy(x) %计算相邻点的：分点的坐标 [m,n]=size(x); xx=zeros(2*m-l,2*n-l); %扩人矩阵来存储节点的坐标 fori=l:m %将原矩阵在矩形顶点的x,y放在扩大后的矩阵相应处 for j=l:n xx(2*i-l ,2*j-l)=x(i,j); end end for i=l:m %计算原始行相邻点屮点的x,y坐标 for j=l:n-l xx(2*i-l ,2*j)= l/2*(xx(2*i-l ,2*j-l )+xx(2*i? 1,2*j+1)); end end for i=l:m-l %计算所有列的相邻点屮点的x,y坐标 for j=l:2*n-l xx(2*i,j)= l/2*(xx(2*i-l ,j)+xx(2*i+1 ,j)); end end xxl=xx; %得到二分点插值后的人矩阵 end 二分点插值函数的matlab代码实现 然后进行插值计算出二分点处的Z的值。由上面所说的 P(x, y) = I】(x, y)Z, +l2(x, y)Z2 +13 (x, y)Z3 即为插值函数。 Matlab代码实现如下： function zzz=Linear(x,y,z,xx,yy) % 计算二分点插值后的Z的值 [m,n]=size(xx); zzz=zeros(m,n); %用来存储Z值的矩阵 DO二ones(3,3); %三阶矩阵，用于基函数计算 [m,n]=size(x); for j=l:n-l for i=l:m-l D0(l,l)=x(i,j); % DO即上面所说的第一种三角形 D0(l,2)二 y(i,j); D0(2,l)二x(i,j+l); D0(2,2)=y(i,j+l); D0(3,l)=x(i+l,j); D0(3,2)=y(i+l,j); D1=DO; Dl(l,l)=xx(2*i-l,2*j); %D1 用于计算 ll(x,y) Dl(l,2)=yy(2*i-l,2*j); D2=D0; D2(2,1 )=xx(2*i? 1,2*j); % D2用于计算 12(x,y) D2(2,2)=yy(2*i-l,2*j); D3=D0; D3(3,l)=xx(2*i-l,2*j); % D3用于计算 13(x,y) D3(3,2)=yy(2*i-l,2*j); ll=det(Dl)/det(DO); 12=det(D2)/det(D0); 13=det(D3)/det(D0); zzz(2*i-l ,2*j)=ll *z(i,j)+12*z(i,j+1 )+13*z(i+ l,j); %M1 的 Z值 Dl(l,l)=xx(2*i,2*j); Dl(l,2)=yy(2*i,2*j); D2(2,l)=xx(2*i,2*j); D2(2,2)=yy(2*i,2*j); D3(3,l)=xx(2*i,2*j); D3(3,2)=yy(2*i,2*j); ll=det(Dl)/det(DO); 12=det(D2)/det(D0); 13=det(D3)/det(DO); zzz(2*i,2*j)=Il *z(i,j)+12*z(i,j+1 )+13*z(i+1 ,j); % M5 的Z值 Dl(l,l)=xx(2*i,2*j-1); Dl(l,2)=yy(2*i,2*j-1); D2(2,l)=xx(2*i,2*j?l); D2(2,2)=yy(2*i,2*j-1); D3(3,l)=xx(2*i,2*