MATLAB的图像处理介绍解析.pptVIP

对得到的中轴线坐标进行拟合的方法，在当年的论文中很多，这里只介绍使用多项式进行拟合，matlab中命令是polyfit。 m从零递增时， 但当m达到某一个值M之后，对mM，有 即m继续增加时 与 的值相差不大。此时就是所要确定的最小二乘多项式次数。 若m过小则可能不足以表示函数的固有规律，而m过大时又失去了光滑性，且产生抖动现象，实践中往往采用下面的方法确定最小二乘多项式的次数m，令： 显著减少。 根据上述理论定出多项式阶数：f1为8阶多项式，f2为6阶多项式。 YZ平面上：y=f1(z) XZ平面上：x=f2(z) Xy平面上：由上面两式联立得到。 具体程序 绘图结果： 拟合中轴线在yoz面投影 拟合中轴线在xoz面投影 拟合中轴线 具体程序 拟合中轴线在xoy面投影 将拟合前和拟合后的图放在一起比较： 包络球沿拟合中轴线形成血管的空间效果图 由于前面模拟出来的血管图像都是原始切片图，现在考虑已经得到了血管拟合中轴线以及血管半径，那么可以得到一个较为形象的模拟血管，下面就模拟血管。 下面就使用matlab结合圆形算法和曲面管道画法得到图形： 下面是将中轴线也标注在血管中得到图形： 下面是将中轴线也标注在血管中并将血管半透明化得到图形： 灵敏度分析的思路与方法 （1）由于机器中的像素点是由单个像素构成所以是一种离散的状态，所以距离以及边缘的判别都存在误差，可以增加灵敏度分析，思路是：在取边缘和中轴线坐标时建立扰动机制，分别把像素点活动一个位置，有8种选择，再对结果进行分析。 （2）图像中的截面，圆面和距离都是由像素点计算和表示的，存在误差，也可以采用上述的灵敏度分析的方法进行判别。 1 2 3 4 5 6 7 8 （3）待选点理论上都是中轴线的坐标点，都可以选择，所以选择那个点作为坐标点也会引起误差，灵敏度分析的方法可以采用逐个实验，再分析结果。 其他算法的总结 题目所使用的工具多种多样，有matlab，mathmatic，fortain语言和c语言，也有用photoshop叠加100幅图的。 题目所使用的方法也很多： （1）有用最小二乘法拟合每张图切片的上下曲线，并计算图片与实物的实际比例的。 （2）求最大内切球方法：血管中轴线与截面的交点必落在截面边界上的一对点的连线的中点上，过这对点的截面边界的两条切线是相互平行的。也就是最大内切球的原理。 （3）除去“弯月”的两端外，图形的宽度保持不变，这一宽度就是直径，选择均匀的中段在左侧定下A，在另一边找B可以得到直径。 （4）先求粗略的半径最大值rmax和最小值rmin，再求出切片的边缘点的位置，根据rmin,求出切片内所含的标准圆的圆心O可能存在的位置区域，用二分法不断缩小标准圆半径r0可能的取值范围，从而求出r0。最后用求出切片内所含的标准圆的圆心的位置。 （5）将图的边界点分为上弧、下弧。对于上弧固定一点，求出该点到下弧点的距离，当取遍上弧中每一个点，得到的最短距离就是管道直径；处理多个待选点的方法：取一个待选点到其余待选点距离之和最小的点作为截面的中心坐标点。 （6）最大覆盖法：找到的最佳圆心位置和半径长度是由此所形成的圆面能最大面积的覆盖截面的图形，具体实现是对区域内部的点进行遍历找到能覆盖面积最大的内切圆。 （8）此外还有截面投影法；搜索点法；搜索点法里的步长减半法；边缘相交相切法，最大最小半径法；枚举法；外推法；滚球法；变换法；细化法等大家也可以拓展自己的思路想出更好的方法。 （7）投影法：将各截面的图像叠加在xoy平面上，形成血管再xy平面上的投影，其中心线就是血管中轴线在xy平面上的投影。类似的再xoz和yoz投影。并在投影中计算半径值。 在拟合的方法选取中，方法也比较多，多是采用多项式拟合，并计算了多项式的阶，阶不同结果也不同；也有的采用Bezier曲线进行拟合。 求解管道半径有多种方法，建模竞赛的优秀论文中以“平均法”和“抽样法”居多。平均法是求出每张图像内的最大内切圆半径，再取管道半径为它们的算术平均值。但此方法是一种搜索算法，圆心遍历区域内点，半径由小到大搜索，其计算量相当庞大。 　　抽样法就是利用滚动球半径是常数，取前几片横断面图像内的最大内切圆半径的平均值为管道半径的值。此方法计算较少，但是抽样样本数目选取的合理性较难确定，存在计算误差。 　　一个启发算法：先将100张图片叠合形成一张图片，由于切片垂直Z轴的，此图片是血管在XOY面上的投影图像。因此，它也是滚动球在XOY面上的投影――滚动圆，沿中轴线的投影线滚动形成的二维包络图，且滚动球的半径与滚动圆的半径相同，因此只需求出滚动圆半径。具体图像叠合的方法可以用matlab，也可以用Photoshop控件进行叠加。然后在此投影面上使用上下弧搜索的方法得到半径。　　 　　这样做