利用插值函数求三维图像的二分点,三分点与其到定曲面等距离的曲面.docVIP

二元函数的分片段线性插值 一 简介 由一元函数如y=f(x)的分段线性插值，我们知道，插值的关键就是要找到基函数，对于二元函数也一样，我们需要找到基函数，但此时的基函数，应该是x，y的函数，而不再是简单的x的函数，因此我们需要找到另外的方法来将三维的曲面分为小片段，从而构造基函数。 二 基函数和插值函数的构造 1. 三维曲面的基函数 由曲面的分片段构造的基函数方法有分片段线性插值基函数，分片段样条插值基函数两大类。其中分片段线性插值又分为三角形，矩形等等，分片段样条插值又分为B样条，三次样条等等。这里根据给出来的数据的特点，我们选择了分片段线性插值的三角形插值来计算相邻节点的二分点，三分点的值。 线性三角插值 如上图所示的三角形，已知三角形的三个顶点P1， P2，P3处的函数值分别为Z1=f(x1,y1)，Z2=f(x2,y2)，Z3=f(x3,y3)。我们需要求出线性函数P(x,y)=ax+by+c，使其满足P(x1,y1)=Z1，P(x2,y2)=Z2，P(x3,y3)=Z3。由此得线性方程组问题： 显然该方程组有唯一解的充要条件是系数矩阵行列式不等于0，该行列式恰好是三角形面积两倍。即只要当P1，P2，P3不共线时，三角形区域上的二元线性插值函数是存在的也是唯一的。 三角形的面积的两倍我们可以用行列式来表示，，那么可以利用行列式的性质来构造基函数如下： 显然，三个基函数分别满足插值条件： 三角形区域上线性插值函数的图形是空间上的一个三角形，此空间三角形的顶点恰好是(x1,y1,Z1)，(x2,y2,Z2），(x3,y3,Z3)，由此我们就构造出来了三角形插值的基函数。 相邻节点的二分点，三分点的x，y坐标 由于所给的数据比较有特点，它的数据结构上的相邻也是空间结构上的相邻，因此我们在数据结构中求解相邻节点的二分点，三分点。 所给的数据x,y,z结构均为矩阵，因此如下面的图形所示，我们用矩形表示x，y的矩阵，矩形的顶点为P1 P2 P3 P4,我们计算四个顶点的二分点，即中点的坐标值。如果我们连接P2 P3，会有两种形状的三角形，且M5位于三角形斜边的中点处。 二分点 三分点 依照同样的道理，我们可以做出三分点。M4,M8为M1，M11的三分点，M5,M9为M2，M12的三分点，那么在计算时，我们需要首先计算矩形四条边上的三分点，然后才能计算矩形中的三分点。连接P2 P3我们依然可以得到两种形状的三角形，此时对角线上有两个元素M5 M8，而且分成的两个三角形内部各有一个点,如第一种三角形内部点为M4，第二种三角形内部点为M9。当然上面给出的矩形相当于一个二阶矩阵，类似可得[m,n]阶矩阵的情形。 从上面我们可以看出，插值完之后矩阵的阶数要变大。因此需要一个大矩阵来存储插值节点的坐标。要注意在大矩阵中相应点处的值还是原来小矩阵中的元素的值。 4. 二分点，三分点坐标的matlab代码实现 function xx1=inter2xy(x) %计算相邻点的二分点的坐标 [m,n]=size(x); xx=zeros(2*m-1,2*n-1); %扩大矩阵来存储节点的坐标 for i=1:m %将原矩阵在矩形顶点的x,y放在扩大后的矩阵相应处 for j=1:n xx(2*i-1,2*j-1)=x(i,j); end end for i=1:m %计算原始行相邻点中点的x,y坐标 for j=1:n-1 xx(2*i-1,2*j)=1/2*(xx(2*i-1,2*j-1)+xx(2*i-1,2*j+1)); end end for i=1:m-1 %计算所有列的相邻点中点的x,y坐标 for j=1:2*n-1 xx(2*i,j)=1/2*(xx(2*i-1,j)+xx(2*i+1,j)); end end xx1=xx; %得到二分点插值后的大矩阵 end function xx1=inter3xy(x) %计算相邻点的三分点的坐标 [m,n]=size(x); xx=zeros(3*m-2,3*n-2); %扩大矩阵来存储节点的坐标 for i=1:m %原始位置 for j=1:n xx(3*i-2,3*j-2)=x(i,j); end end for i=1:m %原始行位置 for j=1:n-1 xx(3