三维血管的重建
基于256CT影像数据的冠脉血管三维实体模型重建
基于256CT影像数据的冠脉血管三维实体模型重建摘要：目的：讨论应用基于256CT影像数据的冠脉血管三维实体模型重建技术对于冠脉血管疾病的诊断与治疗的作用。
方法：选取本院心内科2017年1月至2017年10月收治的心脏冠脉血管疾病患者45例作为研究对象，选取同期入院的无心血管疾病患者45例作为健康组，使用星云工作平台对256CT的图像结果进行三维重建，评价重建后的诊断符合率。
结果：根据三维图像判定的结果，对于45名病人判定的正确率为100%，而对于健康人，其中3人被错误的判定为冠心病患者，正确率为93.3%，两组正确率的差异具有统计学意义。
结论：应用基于256iCT的心脏冠脉血管三维重建技术，对于冠脉血管状态的重现性好，并且更加直观，根据三维重建结果对于疾病判定的灵敏度与特异度均较高，提示此技术的优越性，对于其进一步的应用仍然值得探索。
关键词：影像数据；三维重建；冠脉血管心脏冠脉血管为心脏的营养血管，其血管分布繁杂并且在个体之间走行多变，而临床对于冠脉血管病变的诊治一直为技术性非常强的工作，并且对于冠脉血管畸形的患者，其病变往往十分细小，诊断更加困难，而数字三维成像技术为心脏冠脉血管疾病的影像诊治带来了新的思路，本研究讨论了应用基于256CT影像数据的冠脉血管三维实体模型重建技术对于冠脉血管疾病的诊断与治疗的作用，为技术的发展以及冠脉血管疾病的诊治提供了思路。
1.材料与方法1.1研究对象的选取选取本院心内科2017年1月至2017年10月收治的心脏冠脉血管疾病患者45例作为研究对象，选取同期入院的无心血管疾病患者45例作为健康组。
纳入标注为（1）既往无严重的心瓣膜病，无恶性心律失常者；（2）病人病情能够耐受CT诊断，对于碘海醇无过敏者；（3）自愿参与本次研究，并且签署了知情同意书的患者。
1．2 方法所有患者均由护理人员辅助使用飞利浦256Ict进行心脏扫描，扫描前嘱患者保持心情平静，并且由护理人员教授患者在扫描之中进行屏气，待患者心率稳定于55到70次每分时进行扫描，静脉输入碘海醇作为造影剂。
血管的三维重建课件
m(:,:,b+1)=imread([int2str(b),'.bmp']); end
血管的三维重建
-10-
建模方法思想
需考虑的细节：
2）何谓边界点？
四邻域的概念 找边界点坐标的算法
血管的三维重建
-2-
Z=0
Z=49
Z=98
Z=1
Z=50
血管的三维重建
Z=99
-3-
假设
1）血管的表面是由半径固定、球心沿着某一曲 线（称为中轴线）的球滚动而形成的包络面。 2）中轴线上任两点处的法截面圆不相交。 3）管道中轴线与每张切片平面有且只有一个交 点。
血管的三维重建
-4-
图象的矩阵表示
后面的切片。
Z=57
Z=60
血管的三维重建
-22-
误差大的原因和改进途径
误差大的原因：
1）图象误差 实际图象边界上的点是连续的，在转换成
bmp图象时，象素表示的图象边界是离散的，成 锯齿状，与实际图象有误差（舍入误差）。 2）同一张切片上的最大内切圆不唯一
解决办法：
1）方法一：取平均 求出同一张切片上的所有最大内切圆的圆心，然 后求平均值。
问题重述
断面可用于了解生物组织、器官等的形态。例如，将 样本染色后切成厚约1m m的切片，在显微镜下观察该横断 面的组织形态结构。如果用切片机连续不断地将样本切成 数十、成百的平行切片， 可依次逐片观察。根据拍照并 采样得到的平行切片数字图象，运用计算机可重建组织、 器官等准确的三维形态。
假设某些血管可视为一类特殊的管道，该管道的表面是 由球心沿着某一曲线（称为中轴线）的球滚动包络而成。 例如圆柱就是这样一种管道，其中轴线为直线，由半径固 定的球滚动包络形成。
2006三维血管重构_汤志高
6.2 模型分析
目标分析： 根据 6.1.1 的证明可知。 切片最大内切圆等价于管道的横截面， 最大内切圆的半径等 于管道的半径。取坐标系的 Z 轴垂直于切片，即切片平行与 XY 平面。只要在切片平面 内求解最大内切圆的半径即可。 约束分析： (1) 在切片内部寻求一个圆，使它成为最大内切圆，圆心一定在切片边界组成的封闭 区域 Di （ i 1 N ，N 表示切片的数量）内 ( xi , yi ) Di (2) 最大内切圆的半径要小于圆心到边界上任意一点的距离。 设第 i 个切片的最大内切 圆的圆心坐标为 ( xi , yi ) ，第 i 张切片边界的第 j 个点的坐标为 ( X ij , Yij ) ，已经由模型准备 6.1.2 中获得。所以，最大内切圆半径应该小于圆心到这些离散点的任意一点的距离
血管三维重建
作者：汤志高 等 hihihi@ ［摘要］
本文通过对样本断面的数字化处理，运用计算机重建血管三维形态，得出 血管轴线方程和管体半径。 对断面特征图像文件数值化：视位图图像黑色象素点为 0，白色为 1，在 MATLAB 中建立图像矩阵 A。对 A 进行转化，建立图像轮廓矩阵 X ij , Yij ，通过
X ij , Yij 绘制出血管三维重建图形。
在求解管道半径时，结合三维图形与断面图形特征，发现管道半径 Rg 即 为断面图形最大内切圆半径 r ，轴线与断面的交点即为最大内切圆圆心，并在 6.1.1 中给出了严格证明。就此，将问题转化为求解断面的最大内切圆的半径。 基于转化的目标，首先加入内圆约束（内切圆半径不大于圆心到断面边界 上任意一点的距离） ； 然后采用约束算法 “四向扫描法” （见 6.4.1） 使约束加细， 限制圆心坐标 ( xi , yi ) D 。主要原理是对图象矩阵 A 再次转化，将每个断面的 边界都分割成四部分，分离出各断面图形的上、下、左、右轮廓，建立四向轮 廓矩阵，对圆心坐标从四个方向约束，当圆心坐标同时满足四个约束时，圆心 一定位于断面内部。最后，以 ( xi , yi ) 为决策变量，内切圆半径 r 最大为目标建 立非线性规划模型，采用 LINGO 软件编程求解，最终管道半径为 29.7069； 。 ( xi , yi ) 坐标集（见附录 10.1） 在求解中轴线的过程中，首先建立 zi 向量（即各内切圆圆心 ( xi , yi ) 的 z 坐 标） ，运用 100 张断面与中轴线的交点在 xoz 和 yoz 平面的投影，建立 xoz 面的 投影坐标集 xi , zi ，采用高次非线性多项式拟合得到平面曲线方程 x f z ， 同理对 yoz 面的投影坐标集 yi , zi 采用同样方法得到平面曲线方程 y f z ， 以 z 为参数， 由经过上述两条曲线的空间曲面方程的联立， 建立了中轴线方程。 最后根据空间中轴线的方程，可以分别画出中轴线在 xoy 面， xoz 面， yoz 面的轴线投影图和中轴线的空间图象。
【 数学建模竞赛】血管的三维重建模型g精品
血管的三维重建模型摘要：本文对血管三维重建中，中轴线及球的半径确定问题进行了讨论。
首先，根据问题及图象处理提取有效数据，给出两种可行算法，利用上述数据建立了最大最小方法和二次规划方法。
搜索中心点，并给出全局和局部搜索，得到各切片中心点坐标(见表1)，并通过插值方式得到中轴线图象及其各投影。
最后对模型给出检验方式。
一 、问题的重述假设某些血管可视为一类特殊的管道,该管道的表面是由球心沿着某一曲线 (称为中轴线)的球（命名为包络球）滚动包络而成。
现有某管道的相继100张平行切片图象,记录了管道与切片的交。
假设：管道中轴线与每张图片有且只有一个交点；球半径固定；切片间距以及图象象素的尺寸均为1.取坐标的Z 轴垂直于切片，第1张切片为平面0=Z ，第100张切片为平面99=Z . 计算管道的中轴线与半径，给出具体的算法，并绘制中轴线在XY 、YZ 、ZX 平面的投影图。
二、模型假设与符号说明1、 基本假设：（1） 该管道的表面为一定长半径的球沿一固定的曲线运动所得曲面族包络的光滑表面。
（2） 该管道的中轴线连续而且光滑。
（3） 该管道的中轴线与每个切面有且只有一个交点。
（4） 图象象素的尺寸为1. （5） 切片的间距尺寸为1.2、 符号说明：L 中轴线R 包络球的半径()z y x O i ,, 中轴线与第i 个切片的交点（定为此切片的中心）i S 第i 个切片切得的图形 i D 第i 个切片的图象数据矩阵三、问题分析及建模准备 问题分析：通常血管的表面可认为是连续且光滑的曲面，断面可用于了解其形态等特性。
本问题给出的是一些离散的切面，要求重建出原图中轴线和求出包络球半径。
因为每一个切面与中轴线L 有且只有一个交点i O ，如果找出所有i O ,就可以用插值或拟合的方式作出L 的近似图象，其在坐标平面上的投影就很容易画出。
问题的关健转变为求每个平面上的i O . 建模准备：1、 图象的读取由于切片图象中只有黑、白两种颜色的象素,而且所给的BMP 格式图象文 件是512×512象素的.因此,把图象读取为一个512×512的数字矩阵；用数字1表示黑色的象素，用数字0表示白色的象素。
数学建模血管的三维重建问题
A题血管的三维重建问题摘要：本论文讨论基于切片的血管三维重建问题。
其背景是：采取存储二维切片信息，使用时再利用切片信息重建原物体三维形态的方法，可以有效地保存和利用三维信息。
此技术在实际中有很大的用途，在医学和其他领域有广泛的应用。
如要将人体全部三维信息，包含内部错综复杂的结构，完整地存储在计算机中，以现在的技术也是有一定难度的，但若改用存储人体切片信息，使用时重建再现的方法，则是利用现有技术可以解决的。
本论文基于题中对血管形态的假设，建立管道中轴线参数方程，并综合考虑实际情况中由于切片厚度及数字图像离散化带来的偏差，通过在每张切片图像中搜索其中阴影区域所能包含的最大圆面，确定管半径为R=29，在此基础上，将每张切片图像中阴影区域所能包含的半径大于等于R的圆面圆心作为中轴线与各切片交点（即中心点）的候选点集合。
本模型使用了三种改进算法对该候选点集进行筛选以确定实际交点。
最终迭代算法简述如下：1．对每个切片，建立中心点的候选点集，并取点集的中位点为中心点初值2．利用得到的中心点建立中轴线方程3．利用中轴线方程推导导数信息，根据导数信息比例选取中心点的候选点集的某点作为中心点的新值4．重复步骤2、3，直至结果达到较稳定状态为止5．输出中心点及中轴线方程在模型建立中，对选取侯选点集、求中位点、利用导数信息进行比例选取均给出完整的算法，并且对半径确定、候选点选取、采用导数作为比例选取依据等问题给出详尽的证明。
考虑到实际血管的中轴线应充分光滑，计算最终中轴线参数表达式时采取了六阶多项式拟合。
最后用还原的血管形态模拟切片过程可以得到一系列数字图像，与原切片图像进行比较，可以检验模型的合理性及精度。
该模型最终计算结果如下。
血管中轴线示意图从模型结果中看出，中心点分布均匀稳定，模拟检验的切片数字图像与原切片的数字图像吻合较好，模型结果精度及稳定性符合要求。
本模型算法简明，理论严密，比例选取算法使结果中心点尽可能收敛于真实中心点，迭代算法保证了结果的精度和稳定性，符合题目要求。
【 数学建模竞赛】血管的三维重建模型g
血管的三维重建模型摘要：本文对血管三维重建中，中轴线及球的半径确定问题进行了讨论。
首先，根据问题及图象处理提取有效数据，给出两种可行算法，利用上述数据建立了最大最小方法和二次规划方法。
搜索中心点，并给出全局和局部搜索，得到各切片中心点坐标(见表1)，并通过插值方式得到中轴线图象及其各投影。
最后对模型给出检验方式。
一 、问题的重述假设某些血管可视为一类特殊的管道,该管道的表面是由球心沿着某一曲线 (称为中轴线)的球（命名为包络球）滚动包络而成。
现有某管道的相继100张平行切片图象,记录了管道与切片的交。
假设：管道中轴线与每张图片有且只有一个交点；球半径固定；切片间距以及图象象素的尺寸均为1.取坐标的Z 轴垂直于切片，第1张切片为平面0=Z ，第100张切片为平面99=Z . 计算管道的中轴线与半径，给出具体的算法，并绘制中轴线在XY 、YZ 、ZX 平面的投影图。
二、模型假设与符号说明1、 基本假设：（1） 该管道的表面为一定长半径的球沿一固定的曲线运动所得曲面族包络的光滑表面。
（2） 该管道的中轴线连续而且光滑。
（3） 该管道的中轴线与每个切面有且只有一个交点。
（4） 图象象素的尺寸为1. （5） 切片的间距尺寸为1.2、 符号说明：L 中轴线R 包络球的半径()z y x O i ,, 中轴线与第i 个切片的交点（定为此切片的中心）i S 第i 个切片切得的图形 i D 第i 个切片的图象数据矩阵三、问题分析及建模准备 问题分析：通常血管的表面可认为是连续且光滑的曲面，断面可用于了解其形态等特性。
本问题给出的是一些离散的切面，要求重建出原图中轴线和求出包络球半径。
因为每一个切面与中轴线L 有且只有一个交点i O ，如果找出所有i O ,就可以用插值或拟合的方式作出L 的近似图象，其在坐标平面上的投影就很容易画出。
问题的关健转变为求每个平面上的i O . 建模准备：1、 图象的读取由于切片图象中只有黑、白两种颜色的象素,而且所给的BMP 格式图象文 件是512×512象素的.因此,把图象读取为一个512×512的数字矩阵；用数字1表示黑色的象素，用数字0表示白色的象素。
血管的三维重建数学建模
血管的三维重建数学建模
首先，血管的三维重建通常是通过医学影像学来实现的。
医学
影像学包括CT、MRI等技术，这些技术可以提供血管的断层扫描图像。
在这些图像的基础上，可以利用图像处理的方法，如边缘检测、分割等技术，来提取血管的形状和结构信息。
其次，几何建模是血管三维重建的关键环节。
在图像处理的基
础上，需要进行几何建模，将提取到的血管形状转化为数学模型。
这涉及到曲面重建、体素网格生成等技术，以及对血管内部结构的
建模。
另外，数学算法在血管三维重建中也起着重要作用。
例如，曲
面重建可以利用曲面拟合算法，体素网格生成可以利用体细胞自动
机等算法。
此外，对血管的分支、扭曲等特征的识别和建模也需要
借助数学算法来实现。
除此之外，血管的三维重建数学建模还涉及到计算机图形学、
计算几何学等领域的知识。
这些知识和技术的综合运用，可以实现
对血管形状、结构和特征的全面建模和重建。
总的来说，血管的三维重建数学建模是一个复杂而多样化的过程，涉及到多个学科和领域的知识。
通过综合运用图像处理、几何建模、数学算法等技术，可以实现对血管的全面、准确的三维重建和建模。
下腰椎前血管解剖形态的三维重建分析
下腰椎前血管解剖形态的三维重建分析目的：分析下腰椎前方血管走行特点，为下腰椎前路椎间融合术及微创手术提供解剖学依据。
方法：回顾性分析60例（男女各30例）腹部三维CT血管造影（3D-CTA ）患者的影像，观察下腰椎前方血管在椎体前的走行特点，测量与腰椎前方血管位置相关的解剖学参数。
结果：三维CT重建图像发现腹主动脉分叉点变异较多，髂总静脉汇合点相对恒定。
腹主动脉分叉角男性（56.4±8.1）°，女性（61.9±8.7）°，髂总静脉汇合角男性（71.1±10.2）°，女性（77.7±11.31）°，腹主动脉分叉点到L5/S1间盘上缘的距离男性（3.39 ±1.02）cm ，女性（3.42±1.34）cm；髂总静脉汇合点到L5/S1间盘上缘的距离男性（1.71±0.76）cm，女性（1.90±1.15）cm；L5/S1椎间隙手术窗大小男性（3.07±0.85）cm ，女性（3.16±1.11）cm。
结论：下腰椎前方血管解剖位置具有多变性，以腹主动脉分叉点较为多见，术前行腰椎3D-CTA检查能明确血管解剖，可为微创手术提供影像学依据。
标签：下腰椎前血管；解剖形态；三维重建随着医疗技术的发展，脊柱前路手术的开展越来越普遍，其中通过腰椎前路行椎间融合术治疗下腰椎不稳、椎间盘源性腰痛等脊柱疾患应用广泛。
脊柱手术的微创化如腹部小切口和腹腔镜技术也成为趋势。
许多学者研究认为腰椎前路和小切口微创手术的主要困难和风险与椎前血管相关[1-4]。
术前确定椎前血管解剖结构有助于减少手术风险。
本研究旨在通过分析60例行腹部三维CT血管造影（3D-CTA）检查患者的三维影像，观察下腰椎前方血管位置及毗邻结构，以期为腰椎前路微创手术提供术前的解剖学影像评估依据。
1 资料与方法1.1 一般资料回顾性分析2008年1月-2009年7月在本院接受多层螺旋CT 检查，行腹部血管造影的60例患者的资料。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
血管的三维重建摘要对于血管的三维重建，本文研究了血管这一类特殊管道的中轴线及其半径的算法，绘制中轴线在XY 、YZ 、ZX 平面的投影图这些问题，问题分为三部分。
针对第一部分，先将100张切片图片在MATLAB 中导出生成0-1矩阵数据，在计算100张切片的最大内切圆半径及对应圆心坐标，为减小误差求100张切片最大内切圆的平均半径41666.29 d 。
中轴线的曲线方程可在MATLAB 中拟合得到。
针对第二部分，得到中轴线曲线方程在MATLAB 中绘制出中轴线方程的空间曲线，之后将其投影在XY 、YZ 、ZX 平面上。
针对第三部分，对100张切片进行叠加重合，得到血管的三维立体图，再通过MATLAB 对血管的三维立体图进行优化完成血管的三维重建。
关键词：MATLAB 软件管道半径中轴线曲线方程一、问题重述1.1基本情况断面可用于了解生物组织、器官等的形态。
如果用切片机连续不断地将样本切成数十、成百的平行切片，可依次逐片观察。
根据拍照并采样得到的平行切片数字图象，运用计算机可重建组织、器官等准确的三维形态。
1.2相关信息假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。
现有某管道的相继100张平行切片图象，记录了管道与切片的交。
图象文件名依次为0.bmp、1.bmp、…、99.bmp，格式均为BMP，宽、高均为512个象素（pixel）。
取坐标系的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）。
1.3提出的问题问题一：计算出管道的中轴线与半径，给出具体的算法。
问题二：绘制中轴线在XY、YZ、ZX平面的投影图。
问题三：绘制血管的三维重建立体图。
二、问题分析血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。
那么求管道半径问题可转化为求每张切片最大内切圆半径问题，中轴线曲线方程的模型建立可用问题一中求出的数据在MATLAB中拟合得到，题中所要求绘制的图像可通过MATLAB绘制得到。
针对问题一，因为血管可视为一种特殊的管道，管道的表面是由球心沿某一中轴线的球滚动包络，因此过球心的截圆半径就是血管半径，即为切片中包含的最大内切圆半径，再取100个最大内切圆半径平均值以缩小误差。
首先将100张图片通过MATLAB导入将其转化为512×512的二维0-1矩阵，并运用MATLAB 中的edge函数求出边界。
在第k张切片中，求出任意内点到边界上所有点的最小距离，在这些最小距离中取最大值，最大值即为最大内切圆半径，对应的内点即为每个切片的最大内切圆的圆心。
对于中轴线的求解，上述问题分析可分别求出100张切片最大内切圆的圆心坐标及半径，为简化方程设t为参数变量，通过MATLAB拟合工具箱进行曲线拟合可以的到关于参数t的中轴线曲线方程。
针对问题二，基于问题一的基础上，通过MATLAB中拟合工具画出中轴线在的三维图像[1]，将其投影在XY、YZ、ZX平面上。
针对问题三，在MATLAB中用拟合工具将100张切片图像进行叠加重合，绘制出血管的三维图像，再用拟合工具进行优化还原得到血管的拟合还原图。
三、模型的假设结合本题实际，为了确保模型求解的准确性和合理性，我们排除了一些因素的干扰，提出以下几点假设：1.管道中轴线与每张切片有且只有一个交点。
2.球半径固定；切片间距以及图象象素的尺寸均为1。
3.假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。
4.所有数据均是准确的，根据像素能够近似的画出图形。
5.对切片的拍照过程中不存在误差。
6.中轴线上任意两点处的法截面圆不相交。
四、符号说明符号符号说明d任意内点到边界点的距离d100张切片最大内切圆半径平均值dd最大内切圆半径d第k张切片kX中轴线在X轴的参数方程)(tY中轴线在Y轴的参数方程)(tZ中轴线在Z轴的参数方程T中轴线参数方程五、模型的建立与求解通过建立数学模型求解出每张切片的最大内切圆半径，并得出最大内切圆的圆心坐标，借此拟合出血管中轴线的曲线方程，同时可绘制出中轴线在各坐标面上的投影，从而实现血管的三维重建。
5.1问题一5.1.1血管半径的模型建立题中将血管通过切片机切成100个平行切片，并给出这100张平行切片的二维平面图，这些平面图是由经过球心的球截面所截而得，其球截面是所有截圆中半径最大的圆，半径最大圆的半径即为血管半径，再求出这100个最大内切圆半径的平均值以减小误差。
求切片的最大内切圆半径，要先找出切片内任意一点到切片边界距离最小值，在这些最小值中找到最大值，由此建出血管半径的模型。
第k 张切片图像第i 个内点))2,(),1,((i i nd 到第j 个边界点))2,(),1,((j j bjd 的距离：)1())2,()2,(())1,()1,((22j bjd i nd j bjd i nd d -+-=利用MATLAB 软件得到图像的最大内切圆半径：(2))max(=]1,1[[],2),,min(=dd zb d d dd 为了减小误差，对100张切片的最大内切圆半径求平均值：)3(1001100∑=kK d d 5.1.2血管半径的求解（1）首先将切片图像导入到MATLAB 中。
（2）提取切片图像的轮廓，并对图像进行处理，将图像像素点转化为512×512的矩阵[2]。
（3）将矩阵数值导入求切片图像内切圆的程序中。
（4）对切片图像内切圆进行筛选，得到最大值。
（5）将最大内切圆的半径和圆心的坐标记录下来。
表一：部分切片的最大内切圆的半径和圆心的坐标X Y Z 半径-1601129.06888-1600228.28427-1602329.01724-1602429.06888-1602529.06888-1602629.06888-1601729.00000-1604829.01724-1601929.00000-16011028.86174-16071128.86174-16081228.86174-16091328.86174-160101429.01724-160121529.01724-160131629.01724-160141729.01724-160161829.01724-160171929.01724-160182029.01724-160192129.01724-160202229.01724-160212329.01724-160222429.01724（注：表一为部分最大内切圆半径及坐标，具体表格详见附录）分析：将100张切片图像的最大内切圆半径进行求平均值，得到血管管道半径为41666.29=d 5.1.3中轴线模型建立在血管管道模型建立中我们可以知道，每张切片的最大内切圆的圆心位于血管的中轴线上。
以题中给的坐标系为模型,根据空间中轴线是由散点组成且z 轴值是逐层单调递增,为简化方程的计算,取为参数变量t ，在MATLAB 中通过拟合得到中轴线的模型[3]：(4))99,210=为参数，()()( ，,,t t t Z t Y Y t X X ⎪⎩⎪⎨⎧===5.1.4管道中轴线的求解（1）根据切片图像最大内切圆圆心的x 轴和y 轴坐标的数据，使用MATLAB 进行拟合得到关于参数t 的方程。
（2）将拟合曲线函数与Z=t 联立，得到中轴线曲线方程。
X 轴坐标随参数t 变化的曲线方程：)5(9.641422.000121.010669.710494.610313.4)(23648510+++⨯-⨯-⨯=---t t t t t t X Y 轴坐标随参数t 变化的曲线方程为：)6(66.18682.54082.001331.010*******.010062.6)(234657-+-+⨯-⨯=--t t t t t t Y 将x 、y 对应坐标与切片数绘制三维散点图，并将x 轴坐标和y 轴坐标的变化曲线X,Y 与切片数Z 联立，得到三维曲线，即中轴线，其方程为：(7))99，，2，0,1=为参数，(66.18682.54082.001331.010*******.010062.69.641422.000121.010669.710494.610313.4234657)(23648510)( t t t Z t t t t t Y t t t t t X T t t ⎪⎩⎪⎨⎧=-+-+⨯-⨯=+++⨯-⨯-⨯==-----通过MATLAB做出中轴线的三维图像：图1中轴线的三维图像5.2问题二在MATLAB中用plot函数可画出中轴线的立体图及在YZ、YX、ZX平面投影拟合曲线。
图2中轴线在X-Y面上的投影图3中轴线在Y-Z面上的投影图4中轴线在X-Z面上的投影5.3问题三在MATLAB中通过plot3函数可将100张平行切片进行累计叠加重合得到血管管道的三维图像，如图5所示。
图5切片叠加而成的血管三维图形根据问题一及问题二中所求得的管道半径和中轴线曲线方程，可绘制出拟合后的三维血管图形（图5），在次通过拟合工具拟合得到血管拟合还原图（图6）。
图6血管拟合还原图六、模型的优点与缺点6.1模型的优点（1）应用MATLAB软件对图片进行了处理，得出了大量数据。
（2）采用平均法对数据进行了科学精确地处理，保证了数据计算的精准度。
（3）采用巧妙的算法，求解结果可以达到很高的精度，能够很好的重建包络，在医学方面有很大作用。
（4）对原始图片进行操作能较为准确的得到血管的三维图，模型结果直观形象，具有非常高的实用和推广价值。
6.2模型的缺点（1）题目自身的像素所决定的误差不可避免。
（2）模型针对性强，不能解决较复杂的三维重建问题，并有一定的局限性。
参考文献[1]邓巍,丁为民,张浩.MATLAB在图像处理和分析中的应用[J].农机化研究,2006(06):194-198.[2]金献珍,吴艳.MATLAB实现数字图像锐化处理[J].商场现代化,2008(36): 12-13.[3]胡庆婉.使用MATLAB曲线拟合工具箱做曲线拟合[J].电脑知识与技术, 2010,06(21):5822-5823.附录附录1：导入图片转化出0-1矩阵im=strcat(num2str('0.bmp'));tp=imread(im);bj=edge(tp);bjd=[];for i=1:512for j=1:512if bj(i,j)==1bjd=[bjd;i,j];%边界点的坐标endendendnd=[];for i=1:512for j=1:512if tp(i,j)==0nd=[nd;i,j];%内点的坐标endendendm=length(nd);n=length(bjd);d=zeros(m,n);for i=1:mfor j=1:nd(i,j)=sqrt((nd(i,1)-bjd(j,1))^2+(nd(i,2)-bjd(j,2))^2);%任意内点到边界点的距离endenddd=min(d,[],2);[d1,zb1]=max(dd);%最大内切圆的圆心附录2100张切片的最大内切圆圆心坐标及半径X Y Z半径-1601129.06888-1600228.28427-1602329.01724-1602429.06888-1602529.06888-1602629.06888-1601729.00000-1604829.01724-1601929.00000-16011028.86174-16071128.86174-16081228.86174-16091328.86174-160101429.01724-160121529.01724-160131629.01724-160141729.01724-160161829.01724-160171929.01724-160182029.01724-160192129.01724-160202229.01724-160212329.01724-160222429.01724-160212529.06888-160212629.06888-160212729.06888-159302829.15476 -159302929.27456 -159293029.27456 -158353129.42788 -157403229.61419 -157403329.61419 -157403429.61419 -156443529.61419 -153553629.73214 -153553729.73214 -153553829.73214 -152583929.73214 -152584029.61419 -150634129.54657 -149664229.54657 -148684329.52965 -148684429.52965 -143784529.52965 -137884629.41088 -137884729.41088 -1161154829.69848 -1151164929.69848 -1151165029.69848 -1141175129.69848 -1141175229.69848 -1131185329.69848 -1121195429.69848 -1111205529.68164 -1111205629.20616 -631515729.41088 -751455829.52965 -811425929.52965 -511566029.54657 -511566129.54657-311626229.61419 -311626329.61419 -311626429.61419 -351616529.61419 -351616629.61419 -261636729.42788 -351616829.41088 -261636929.27456 461637029.42788 461637129.61419 461637229.61419 461637329.61419 651587429.61419 681577529.73214 651587629.73214 811527729.54657 811527829.52965 811527929.52965 1351188029.41088 1361178129.69848 1361178229.69848 1371168329.69848 1381158429.69848 1381158529.69848 1391148629.69848 1391148729.69848 1391148829.69848 1401138929.69848 1401139029.68164 172679129.52965 172679229.52965 172679329.52965 172679429.52965 182439529.73214187249629.61419187249729.61419187249829.61419187249929.614191881810029.42788附录3:绘制中轴线在yz,xy,xz平面上的投影t=linspace(0,99)X=9.665*10^(-9)*t.^6+(-2.245*10^(-6))*t.^5+0.0001684*t.^4+(-0.005139) *t.^3+0.1188*t.^2+t.*(-0.9457)+258.7;Y=(-1.191*10^(-8))*t.^6+3.724*(10)^(-6)*t.^5+(-0.0004097)*t.^4+0.0188 1*t.^3+(-0.381)*t.^2+t.*2.944+411.2;Z=t;plot(Z,Y,'-R');ylabel('Y轴');zlabel('Z轴');title('中轴线在yoz平面上的投影');t=linspace(0,99)X=9.665*10^(-9)*t.^6+(-2.245*10^(-6))*t.^5+0.0001684*t.^4+(-0.005139) *t.^3+0.1188*t.^2+t.*(-0.9457)+258.7;Y=(-1.191*10^(-8))*t.^6+3.724*(10)^(-6)*t.^5+(-0.0004097)*t.^4+0.0188 1*t.^3+(-0.381)*t.^2+t.*2.944+411.2;Z=t;plot(X,Y,'-R');ylabel('Y轴');xlabel('X轴');title('中轴线在xoy平面上的投影');t=linspace(0,99)X=9.665*10^(-9)*t.^6+(-2.245*10^(-6))*t.^5+0.0001684*t.^4+(-0.005139) *t.^3+0.1188*t.^2+t.*(-0.9457)+258.7;Y=(-1.191*10^(-8))*t.^6+3.724*(10)^(-6)*t.^5+(-0.0004097)*t.^4+0.0188 1*t.^3+(-0.381)*t.^2+t.*2.944+411.2;Z=t;plot(X,Z,'-R');zlabel('Z轴');xlabel('X轴');title('中轴线在zox平面上的投影');附录4：绘制中轴线的三维图像t=linspace(0,99)X=9.665*10^(-9)*t.^6+(-2.245*10^(-6))*t.^5+0.0001684*t.^4+(-0.005139) *t.^3+0.1188*t.^2+t.*(-0.9457)+258.7;Y=(-1.191*10^(-8))*t.^6+3.724*(10)^(-6)*t.^5+(-0.0004097)*t.^4+0.0188 1*t.^3+(-0.381)*t.^2+t.*2.944+411.2;Z=t;plot3(X,Y,Z);xlabel('X轴')ylabel('Y轴');zlabel('Z轴');title('中轴线的三维图像');scatter(x,y,z,'.');附录5：切片累加形成的血管三维图for k=0:99photo=imread(strcat(int2str(k),'.bmp'));bianjie(:,:,k+1)=edge(photo,'sobel');endfor k=0:99for i=1:512for j=1:512if bianjie(i,j,k+1)==1plot3(i,j,k+1,'r-');hold onendendendendgrid ontitle('切片累加形成的血管三维图')rotate3dhold off附录6：导入数据jg=ones(100,4);for k=0:99data=imread(strcat(int2str(k),'.bmp'));%导入图片和数据for i=1:512for j=1:512data(i,j)=1-data(i,j);endendlk=edge(data,'sobel');gj=bwmorph(data,'skel','inf');[x0,y0]=find(lk);[a0,b0]=find(gj);m=length(a0);n=length(x0);d=zeros(m,n);za=ones(m,1);for i=1:mfor j=1:nd(i,j)=sqrt((a0(i)-x0(j))^2+(b0(i)-y0(j))^2);endza(i,1)=min(d(i,:));end[zd,zdxh]=max(za);x=a0(zdxh)-256;y=b0(zdxh)-256;z=k+1;jg(k+1,1)=x;jg(k+1,2)=y;jg(k+1,3)=z;jg(k+1,4)=zd;end附录7：绘制血管拟合还原图cenn=zeros(100,3);cenn(:,1)=jg(:,1);cenn(:,2)=jg(:,2);cenn(:,3)=jg(:,3);t=linspace(0,pi,50);p=linspace(0,2*pi,50);[theta,phi]=meshgrid(t,p);for i=1:100;x=29.4159*sin(theta).*cos(phi)+cenn(i,1);y=29.4159*sin(theta).*sin(phi)+cenn(i,2);z=29.4159*cos(theta)+cenn(i,3);hold on;surf(x,y,z)axis([-200200-50200-50150])endXlabel('x轴');Ylabel('y轴');Zlabel('z轴'); title('血管拟合还原图');hold offshading interp。