[最新]血管三维重建.doc

血管的三维重建 摘要 本题要求我们通过100张血管的平行切片图像，运用matlab进行数字图像处理,计算血管的中轴线和半径，绘制中轴线在XY、YZ、ZX平面的投影图，以达到将血管的三位形态表现出来的目的。 血管可视为一种特殊的管道，而管道的表面可视为由球心沿着某一曲线(中轴线)的球滚动包络而成。血管中轴线与每张切片有且只有一个交点，且等径。首先，我们确定每个切片中截面的轮廓线所包含区域内最大内切圆的圆心即为切片和血管中轴线的交点，半径即为血管半径；然后，我们用imread函数读取100张血管的切面图,位图文件就转换成了二值矩阵（像素点取值只限于0，1）,用函数edge可以提取出边界矩阵，用函数bwmorph可以提取出骨架矩阵。我们抽取其中部分平行切片图，matlab编程计算，得出20组截面的最大内切圆半径和坐标（表1）的数据，以这些数据为基础，运用高次多项式拟合，绘制出中轴线在XY、YZ、ZX平面的投影图，对投影图进行3D旋转，即可以看到血管中轴线的空间大致形态。 关键词：MATLAB数字图像处理 最大内切圆 三维重建 matlab 一 问题重述 断面可用于了解生物组织、器官等的形态。例如，将样本染色后切成厚约如m的切片，在显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十、成百的平行切片，可依次逐片观察。根据拍照并采样得到的平行切片数字图象，运用计算机可重建组织、器官等准确的三维形态。 假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线(称为中轴线)的球滚动包络而成。例如圆柱就是这样一种管道，其中轴线为直线，由半径固定的球滚动包络形成。 现有某管道的相继100张平行切片图象，记录了管道与切片的交。图象文件名依次为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平面的投影图。 二 符号说明 截面内某点的X轴坐标 截面内某点的y轴坐标 切片轮廓线上某点的X轴坐标 切片轮廓线上某点的y轴坐标 坐标为的内点到轮廓线的距离 第张切片图的最大内切圆半径 三 模型假设 1、假设某些血管可视为一类特殊的管道,该管道的表面是由球心沿着某一曲线(称为中轴线)的球滚动包络而成； 2、管道中轴线与每张切片有且只有一个交点； 3、球半径固定; 4、切片间距以及图象象素的尺寸均为1； 5、血管无严重扭曲。 四 问题分析 首先，要计算血管的中轴线和半径，要从100张平行切片图象入手。因为血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线(称为中轴线)的球滚动包络而成，而球的任意截面都是圆，圆的圆心是一个非常重要的特征，可以从圆心找到中轴线，平行切片中血管的截面粗细不均，经查阅资料（见参考文献）发现：每个切片中截面的轮廓线所包含区域内最大内切圆的圆心即为切片和血管中轴线的交点。因此解决问题的关键就是找到每张平行切片图象的最大内切圆的圆心坐标和半径。 我们抽取其中部分数据，用matlab编程进行计算，得到20组截面的最大内切圆半径和坐标（表1）；然后，运用多项式拟合，即可绘制出中轴线在XY、YZ、ZX平面的投影图，对投影图进行3D旋转可以看到血管中轴线的空间大致形态。 五 模型建立与求解 1、血管三维重建： 象素，是图像在计算机显示中的度量单位。BMP格式，即位图文件，整幅图可视为一个数字矩阵。根据数据矩阵和图像象素颜色匹配关系，MATLAB中图像可分为：索引图像、灰度图像、二值图像和RGB图像。该题中像素点取值只有0和1，因此为二值图像。在matlab中利用函数imread可将图像读成一个矩阵的形式。编写matlab程序（见附录一），将每张平行切片图象机器轮廓线转化为0-1矩阵，0为黑色，1为白色。 （1）将存放在E:\data1的图0-99读入matlab中： for b=0:99 m(:,:,b+1)=imread([E:\data1\,int2str(