模拟2---A题--血管的三维重建论文.docVIP

模拟2---A题--血管的三维重建论文

A三维血管的重建

摘要

本文根据三维重建理论的知识，运用计算机就血管的二维平行切片重建成血管的三维形态结构问题给予研究。

问题一：计算管道的中轴线与半径，运用MATLAB将附件A中的图像全部传化成0-1的二维矩阵，再通过MATLAB编程，得到每一张图片的边缘轮廓线，找到其边缘轮廓线的最大内切圆，找出圆心（即半径固定的球的球心），将每一张图片得到的圆心拟合起来，就是管道的中轴线(见图4)；将每张图片的最大内切圆的半径相比取其较大者，所得即为管道的半径R，得R=29.2476

问题二：绘制中轴线在XY、YZ、ZX平面的投影图，由第一问得到函数，运用MATLAB绘出图像（见正文）（）

关键词：三维重建理论，二维矩阵，四邻域

一问题的重述

断面可用于了解生物组织、器官等的形态。例如，将样本染色后切成厚约1um的切片，在显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十、成百的平行切片，可依次逐片观察。根据拍照并采样得到的平行切片数字图象，运用计算机可重建组织、器官等准确的三维形态。

假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。例如圆柱就是这样一种管道，其中轴线为直线，由半径固定的球滚动包络形成。

现有某管道的相继100张平行切片图像（模拟2\附件A题bmp.zip），记录了管道与切片的交。图像文件名依次为0.bmp、1.bmp、…、99.bmp，格式均为BMP，宽、高均为512个象素（pixel）。为简化起见，假设：管道中轴线与每张切片有且只有一个交点；球半径固定；切片间距以及图像象素的尺寸均为1。

取坐标系的Z轴垂直于切片，第1张切片为平面Z=0，第100张切片为平面Z=99。Z=z切片图像中象素的坐标依它们在文件中出现的前后次序为

（-256，-256，z），（-256，-255，z），…（-256，255，z），

（-255，-256，z），（-255，-255，z），…（-255，255，z），

……

（255，-256，z），（255，-255，z），…（255，255，z）。

试计算管道的中轴线与半径，给出具体的算法，并绘制中轴线在XY、YZ、ZX平面的投影图。

二模型的假设

1.管道中轴线与每张切片有且只有一个交点；

2.球半径固定；

3.切片间距以及图像象素的尺寸均为1。

ddd=[ddd;j,k];

end

end

end

c{i}=ddd;

end

ddd=[];

save(bianjieshuju.mat,a,b,c,d););

步骤三：中轴线的生成

取出第一张图片并在区域d1中任意取一点dij1找出它与区域边界的最小距离h1，然后再在dij附近找点（如：d(i+1)j1,d(i-1)j1...d(i+1)(j+1)1）确定其与边界的最小距离h’n直到h’n比其附近所有点对应的最小距离都大时，此时的点d’ij1即为要找的焦点；然后取出第二张图片找出d’ij2到边界的最短距离，然后再在d’ij2附近找点直到h”n比其附近所有点对应的最小距离都大时，并取出焦点d”ij2；再取出第三张图片在点2(d”ij2)-d’ij1附近找点直到h(3)n比其附近所有点对应的最小距离都大时，并取出焦点d”’ij3；以此类推即可得到所有的焦点。图片中点到边界最短距离的最大值即为图片内接圆的半径，全部图片的半径的均值可看做是球的半径。

将每张图片的焦点用matlb画出散点图，然后拟合为函数，此函数即为管道的中轴线函数。

血管的中轴线：

0

0

100

200

300

400

500

250

300

350

400

450

0

20

40

60

80

100

x

y

z

图4

具体编程如下（见附录1）

绘制中轴线在XY、YZ、ZX平面的投影图

在上面我们已经找出了每个焦点的坐标,当我们将每个焦点z取0时焦点将回落在xy平面上，在平面直角坐标系中将每一个改进后的点在标注出来就可以拟合出中轴线在xy平面上的投影，同理当我们依次将x,y取0后，就可以拟合出中轴线在yz平面xz平面上的投影图。

其拟合方程如下：

所得图像：

六模型的改进与评价

优点：该模型结合多张切片的信息来确定血管中轴线，可以较精确的从血管切片数字图像中得到血管中轴线和管半径，有较高稳定性和精度。

缺点：本模型的求解过程需要借助大量编程，且未知量多，数字处理量大，在编程过程中程序的调试需要花费大量精力，且易出错。

参考文献

【1】/view/4b0ae621bcd126fff7050b8d.html

【2】