CN101422352B - 一种交互式冠状动脉虚拟血管镜的实现方法 - Google Patents
一种交互式冠状动脉虚拟血管镜的实现方法 Download PDFInfo
- Publication number
- CN101422352B CN101422352B CN2008100800020A CN200810080002A CN101422352B CN 101422352 B CN101422352 B CN 101422352B CN 2008100800020 A CN2008100800020 A CN 2008100800020A CN 200810080002 A CN200810080002 A CN 200810080002A CN 101422352 B CN101422352 B CN 101422352B
- Authority
- CN
- China
- Prior art keywords
- intravascular ultrasound
- ivus image
- path
- frame
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
一种交互式冠状动脉虚拟血管镜的实现方法,属医学成像技术领域。目的是解决冠状动脉可视化诊治问题。技术方案是:它是通过将由近似正交的X射线冠状动脉造影图像获得的管腔三维几何形态信息与由血管内超声获得的管腔横截面数据相融合,得到血管的三维模型,然后运用虚拟现实建模语言交互地描述血管模型,实现内镜漫游模式的冠状动脉可视化。本发明实现了对三维血管模型的交互式访问和显示,为冠状动脉粥样硬化病变的发展、冠心病的可视化诊治、对介入治疗效果评价等的研究,以及医务人员的培训提供了一个理想的平台。
Description
技术领域
本发明涉及一种基于多成像方法融合的交互式冠状动脉虚拟血管镜的实现方法,属医学成像技术领域。
背景技术
虚拟内窥镜技术(Virtual Endoscope,VE)是利用医学影像作为原始数据,综合利用数字图像处理、计算机图形学、科学计算可视化、虚拟现实等技术,重建三维图像,形成虚拟人体组织;然后把视点置入重建出的器官空腔内,借助导航或漫游技术以及伪彩技术来逼真地模拟腔道内镜检查。
目前临床普遍采用的冠状动脉介入影像手段是X射线冠状动脉造影(CAG,Coronary ArteryAngiography)和血管内超声(IVUS,Intravascular Ultrasound),二者是同时进行的,血管造影和血管内超声分别同步显示导管超声探头在管腔内的部位和相应血管壁的结构形态。CAG和IVUS具有优势和不足互补的特点:CAG反映血管腔被造影剂充填后的投影轮廓,能诊断缺血性心脏病及冠状动脉畸形等疾病,而且对冠状动脉内溶栓、PTCA(经皮腔内冠状动脉成形术)等介入手术治疗具有重要意义,但不能提供血管壁的结构信息和病变程度;IVUS可清晰显示血管横断面,根据斑块声学特征进行组织学分型,发现CAG不能显示的血管病变,观察分叉处或血管重叠处的模糊病变等。但是由于采用高频超声探头,影响了探测深度,只能对某一段病变血管进行测量,不能进入严重狭窄的管腔,并且无法确定截面的轴向位置和空间方向。
此外介入成像检测还包括冠状动脉血管镜,它是利用光纤技术的一种微小内窥镜成像技术。但该技术在临床上并未得到广泛接受,原因包括:只能提供管腔表面的形态学资料,不能观察到管壁内的病变深部结构,也不能进行狭窄程度和血流的定量分析;不能用于显像主动脉-冠状动脉开口处的病变和前降支及回旋支近端的病变;从侧孔进入的血流会使视野模糊;导管缺乏可操纵性,限制了显像范围;检查过程中需要暂时堵塞血流,会导致心肌缺血的发生等。
目前,无创性的心血管影像检查主要包括CTA(CT Angiography)和MRCA(Magnetic ResonanceCoronary Angiography)。但是心脏CT检查的主要局限性在于容易产生伪像,影响图像质量。对于MRCA,由于冠脉血管本身较细、扭曲和结构较复杂,且有心脏搏动和呼吸的影响,冠脉周围脂肪组织和心肌组织等信号可影响其显像的结果。同时MR检查过程中虽然没有放射线,相对安全,但有噪音的影响,一些金属植入物(如人工金属瓣、心脏起搏器等)的安全性也受到关注。总之,由于成像原理所造成的不足和技术上的局限,使得到目前为止CTA和MRCA一般仅可作为对心脏综合评价的一种选择方法,或作为冠心病导管造影检查的筛查措施,减少不必要的创伤性检查,以及对心脏手术或介入治疗效果的无创性随访研究,在对冠心病的临床诊治上并不能完全取代介入性的影像检查方法。
综上所述,介入性的CAG和IVUS仍然是临床诊治冠心病的主要影像方法,而且二者具有优势与不足互补的特点。目前还没有一种基于CAG和IVUS图像融合的交互式冠状动脉虚拟血管镜***,能够实现内镜漫游模式的冠状动脉可视化。
发明内容
本发明的目的是克服现有技术的不足、提供一种能够实现内镜漫游模式的冠状动脉可视化的交互式冠状动脉虚拟血管镜的实现方法。
本发明所称问题是以下述技术方案实现的:
一种交互式冠状动脉虚拟血管镜的实现方法,它是通过将由近似正交的X射线冠状动脉造影图像获得的管腔三维几何形态信息与由血管内超声获得的管腔横截面数据相融合,得到血管的三维模型,然后运用虚拟现实建模语言(VRML)交互地描述血管模型,实现内镜漫游模式的冠状动脉可视化,具体步骤如下:
A、同时采集感兴趣血管段的血管内超声IVUS和X射线冠状动脉造影CAG图像:
将机械式超声导管探头置于感兴趣血管段的远端,在匀速等距回撤导引钢丝的过程中,利用血管内超声成像仪以心电ECG门控的方式采集等距的IVUS图像序列,即以心电信号的R波作为触发,仅在每个心动周期内相同的心脏相位处采集图像,可解决冠状动脉IVUS图像序列中的运动伪影问题。同时,利用C型臂单面X射线血管造影机在导管回撤路径的起点拍摄记录相同心脏状态的两幅近似垂直角度的CAG图像;
B、利用采集的IVUS和CAG图像建立血管的三维模型;
C、运用虚拟现实建模语言实现冠状动脉血管重建结果的内镜漫游模式的可视化。
上述交互式冠状动脉虚拟血管镜的实现方法,所述利用采集的IVUS和CAG图像建立血管的三维模型的具体步骤如下:
a、根据两个近似垂直角度的CAG图像,三维重建出超声导管的回撤路径;
b、从CAG图像中三维重建出血管管腔:
将重建出的3D导管路径向左右两个CAG成像平面反投影,得到对应的2D路径,对于2D路径上的每个点,通过在垂直于路径的方向上,寻找灰度梯度的两个极大值,得到血管管腔左右边缘,然后在假设管腔的横截面是椭圆的前提下,完成整个血管管腔的三维重建,该结果用于后续确定各帧超声图像的空间方向;
c、血管内超声图像序列中血管壁轮廓的提取:
在首帧图像中手动选择血管壁内、外膜轮廓上的几个特征点,以连接这些点所形成的多边形作为初始位置,通过snake变形获得血管壁内、外膜的轮廓,分割出血管壁和可能存在的斑块,对于后续帧,则将前一帧的提取结果作为snake的初始位置,完成对连续多帧IVUS图像的分割;
d、确定各帧血管内超声图像的轴向位置:
按照IVUS图像的采集顺序和间距,沿重建出的3D导管回撤路径将各帧IVUS图像顺序排列,确定出各帧图像的轴向位置;
e、确定各帧血管内超声图像的空间方位:
在重建后的3D导管路径上建立各帧超声图像的局部坐标系,即Frenet-Serret标架,三个坐标轴分别为单位切矢量t、单位主法矢量n和单位副法矢量b,导管的位置位于IVUS图像的中心;
将各帧超声图像绕导管旋转至其正确的方向以确定血管内超声图像的空间方位:
用ρ表示血管壁轮廓的重心偏离导管位置的离心向量,把从CAG图像重建出的血管管腔的椭圆轮廓投影到对应的超声图像上,用μ来表示血管管腔椭圆轮廓中心线偏离导管位置的离心向量,ε为向量ρ的模,θ为ρ与μ的夹角,用统计优化方法确定超声图像序列的空间方位,目的是使θ最小:
设定一个固定宽度w的移动窗口,在该窗口中进行统计分析,对于N帧组成的超声图像序列,存在nw=N-(w-1)个移动窗口,在每个窗口位置m处,累计偏心距离∑εm、加权偏心夹角平均值θm以及偏心夹角的加权标准偏差σ(θm)可分别由下式计算:
利用这些数值,在每一个窗口位置处,计算可靠性权重因子:rm=∑εm/σ(εm),在偏心距离较大的位置给予较大的权重因子,同时限制σ(θm)较大的位置,通过下式计算出一个校正偏心角θcorr:
并将其应用到序列的所有图像中,从而获得各帧图像的空间方位;
f、利用基于NURBS曲面拟合的表面提取法完成血管管腔内外表面的绘制。
上述交互式冠状动脉虚拟血管镜的实现方法,所述运用虚拟现实建模语言实现冠状动脉血管重建结果的内镜漫游模式的可视化的具体步骤如下:
①漫游路径的计算:
根据超声导管回撤路径的三维重建结果,以沿导管的第i帧IVUS图像采集点的坐标Pi为当前视点位置,Pi+1为下一个视点位置,以z轴负半轴的方向矢量 为初始视点方向,则在VRML中,当前视点的位置矢量为 旋转轴为当前视点方向为 的单位矢量 旋转轴为 旋转角为
②将IVUS像素数据***到虚拟场景中,并采用半透明的显示方式显示;
③开发交互式的用户图形接口。
上述交互式冠状动脉虚拟血管镜的实现方法,所述超声导管回撤路径的三维重建方法是:首先建立CAG***在两个近似垂直角度的透视投影成像模型,再根据在造影过程中同步记录的距离和角度参数,得到成像***的几何变换矩阵,然后采用三维snake模型技术,snake直接在空间中变形,完成导管回撤路径的三维重建;
上述交互式冠状动脉虚拟血管镜的实现方法,所述冠状动脉造影图像的两个采集角度之间夹角的取值范围为60°至120°,且仅在超声导管回撤路径的起点拍摄一对造影图像。
本发明将由两近似正交角度的单面造影图像得到的血管空间几何信息与由血管内超声图像获得的管腔横截面信息结合起来,充分利用两种成像手段的互补性,完成了血管的准确三维重建,并运用虚拟现实建模语言实现内镜漫游模式的冠状动脉可视化。本发明实现了对三维血管模型的交互式访问和显示,为冠状动脉粥样硬化病变的发展、冠心病的可视化诊治、对介入治疗效果评价等的研究,以及医务人员的培训提供了一个理想的平台。
附图说明
下面结合附图对本发明作进一步详述。
图1是根据本发明方法的三维重建血管的流程图;
图2是根据本发明方法的CAG和IVUS图像采集示意图;
图3是根据本发明方法的造影***在两个角度的成像示意图;
图4是根据本发明方法的各帧超声图像相对方位的确定示意图;
图5是根据本发明方法的超声图像偏心距离和偏心夹角示意图;
图6是根据本发明方法的漫游视点位置的确定示意图;
图7是根据本发明方法的漫游视点方向的确定示意图。
图中各符号为:Image A、Image B、成像平面;s1、s2、两次造影过程中X射线源焦点的位置;s1x1y1z1、以s1为原点的空间坐标系;s2x2y2z2、以s2为原点的空间坐标系;U1V1O1、成像平面A上的直角坐标系;U2V2O2、成像平面B上的直角坐标系;D1、s1到成像平面A的垂直距离;D2、s2到成像平面B的垂直距离;P、空间血管上的点;p1、P点在成像平面A上的投影;p2、P点在成像平面B上的投影;u1、p1在坐标系U1V1O1内的横坐标;v1、p1在坐标系U1V1O1内的纵坐标;u2、p2在坐标系U2V2O2内的横坐标;v2、p2在坐标系U2V2O2内的纵坐标;c(s)、表示3D导管路径的空间参数曲线;C、超声图像中导管的位置,它也是超声图像的中心;OC、椭圆截面轮廓的中心(即在假设血管横截面是椭圆时,基于CAG的三维重建中所对应的血管中心线位置);OI、从超声图像中提取出的管腔截面轮廓的重心;ρ、ρ=OI-C是超声轮廓的重心偏离导管的离心向量;μ、μ=OC-C是血管中心线偏离导管的离心向量;θ、ρ与μ的夹角;Pi、Pi+1、视点;位置矢量;旋转轴;φi、旋转角;z轴负半轴的方向矢量;待求的视点方向;y轴的单位矢量,也即VRML中默认的向上的方向;εi、在x-y平面内的旋转角。
文中所用符号:t、单位切矢量;n、单位主法矢量;b、单位副法矢量;ε、向量ρ的模;w、移动窗口宽度;θm加权偏心夹角平均值;σ(θm)、偏心夹角的加权标准偏差;rm、可靠性权重因子;θcorr、校正偏心角。
具体实施方式
下面结合附图详细说明本发明的步骤:
(1)图像采集:
采集设备包括C型臂单面X射线血管造影机和血管内超声成像仪。
参看图2,IVUS和CAG成像是同时进行的。常规经右股动脉或上臂的肱动脉穿刺,行选择性冠脉造影。在X射线透视图像的指导下***高频超声探头导管,至血管远端。将超声探头与超声成像仪连接去除伪影后,经马达控制匀速等距地回撤导管。当探头导管以1800转/分作360°旋转时连续获得30帧/秒的实时血管切面图像。采用临床常用的、让病人在导管回撤过程中屏住呼吸的方法,减小呼吸运动的影响。采用ECG门控的方式采集超声图像,从而减小心脏运动的影响。
仅在导管回撤路径的起点,采用ECG门控的方式,在相应的心脏相位处拍摄一对近似垂直角度的造影图像。由于采用机械式超声导管探头,超声换能器位于一可弯曲的轴心头端,轴心在外鞘管内旋转,而鞘管是固定不动的,因此可保证回撤路径的稳定。成像过程中记录造影角度和X射线源焦点至接收屏的距离。
(2)造影图像中导管路径和管腔边缘的提取和三维重建:
本发明首先建立CAG***在两个近似垂直角度的透视投影成像模型(附图3)。之后,根据在造影过程中同步记录的距离和角度参数,得到成像***的几何变换矩阵。然后利用三维snake模型技术,完成导管路径的三维重建。
snake模型又称活动轮廓模型(active contour model),是由Kass等在1987年提出的一种变形模型技术(Kass M,Witkin A,Terzopoulos T.Snakes:active contour models.International Journal of Computer Vision,1987,1(4):321-331),近年来在图像处理领域中应用十分广泛,完成图像分割、匹配和运动跟踪。
具体实现方法为:snake的初始位置采用手动取点获得,即在导管的一个投影上手动选取若干采样点(一般选取回撤路径的起点、终点和3~4个中间点即可),然后根据外极约束得到这些点在另一投影上的对应点。由这几组对应点分别求出它们的三维坐标,用直线段连接这些3D点,所得折线作为3D snake的初始位置。
snake模型的能量函数为:
其中c(s)=(x(s),y(s),z(s)),s∈[0,1]是表示导管的三次B样条曲线。式(1)中内部能量Eint的表达式为:
Eint(c(s))=(α|c′(s)|2+β|c"(s)|2)/2 (2)
其中c′(s)和c″(s)分别为c(s)的一阶和二阶导数。内部能量保证曲线的连续和光滑。
外部能量函数Eext是保证snake收敛的外部力,包括两部分,分别对应于左右投影,保证三维曲线在两个角度成像平面上的投影恰好位于对应的导管投影处:
其中IL(u1,v1)和IL(u1,v1)分别是左投影点的灰度和灰度梯度值;IR(u2,v2)和IR(u2,v2)分别是右投影点的灰度和灰度梯度。由于造影图像中,血管的灰度值比背景小,所以权重系数γ取正值。根据透视投影成像的几何关系和外极线约束关系,可推导出u1、v1、u2和v2都是空间点三维坐标c=(x1,y1,z1)的函数:
[u1v1]T=FL(c),[u2v2]T=FR(c) (4)
之后,通过使式(1)的能量函数最小化,snake曲线的最终位置就确定了导管的三维轴线。该方法避免了基于外极约束的两个角度间的逐点匹配,提高了重建精度和运算速度。
按照成像***的几何变换矩阵,将重建出的3D导管路径向左右两个CAG成像平面反投影,得到对应的2D路径。对于2D路径上的每个点,通过在垂直于路径的方向上,寻找灰度梯度的两个极大值,完成对血管管腔左右边缘的提取。之后,在假设管腔的横截面是椭圆的前提下,完成整个管腔的三维重建,该结果为后续确定各帧超声图像的空间方向所用。
(3)血管内超声图像序列中血管壁轮廓的提取:
本发明采用结合动态规划的snake模型完成对IVUS图像序列中血管壁内外膜轮廓的提取。操作者只需在首帧中手动选择目标轮廓上的几个特征点,连接这些点所形成的多边形作为snake的初始位置。对于后续帧,将前一帧的提取结果作为下一帧snake的初始位置,完成对连续多帧IVUS图像的分割,可大大节省计算时间。
(4)IVUS与CAG数据的融合:
这里主要需解决两个问题:确定各IVUS帧的3D轴向位置和空间方位。
(4.1)超声图像三维轴向位置的确定:
在采集超声图像的过程中,采用马达驱动的方式,匀速等距的从远端向近端连续拉出导管。调节拉出导管的速度,即可根据需要调节切面间距。采用CAG图像重建出导管的轴线之后,根据已知的切面间距依轴向将各帧IVUS图像顺序排列,即可确定出各帧图像的轴向位置。
(4.2)各帧超声图像空间方位的确定
本发明利用一种非迭代的统计最优化方法来计算各帧超声图像的空间方位。首先在重建后的3D导管路径上建立各帧超声图像的局部坐标系,即Frenet-Serret标架,三个坐标轴分别为单位切矢量t、单位主法矢量n和单位副法矢量b(附图4),坐标原点是IVUS图像中导管的位置。在完成导管路径的三维重建后,可得到其3D曲线方程c(s),根据微分几何的知识,t、n和b可根据曲线方程计算如下:
其中“×”表示向量的叉乘;“·”表示向量的点乘;c′(s)和c″(s)分别为c(s)的一阶和二阶导数。
导管的位置位于IVUS图像的中心,分割出的目标轮廓的重心一般不与导管位置重合,如附图5所示,其中C点表示导管,OC为椭圆轮廓的中心(即在假设血管横截面是椭圆时,基于CAG的三维重建中所对应的血管中心线位置),01为从超声图像中提取出的血管截面轮廓的重心。采用离心向量ρ表示轮廓的重心偏离导管位置的程度:ρ=OI-C。
由于血管中心线和导管路径不重合,在血管同一位置处的超声图像轮廓和基于造影图像重建出的椭圆轮廓方位不一致,把椭圆轮廓投影到对应的超声图像上。同样采用椭圆轮廓的离心向量μ来表示血管中心线偏离导管位置的程度:μ=OC-C。
超声图像的匹配误差可用向量ρ的模ε和ρ与μ的夹角θ表示。本发明利用统计优化方法确定超声图像序列的绝对方位,目的是使椭圆轮廓和超声轮廓的离心向量间的夹角θ最小。设定一个固定宽度w的移动窗口,在该窗口中进行统计分析。对于N帧组成的超声图像序列,存在nw=N-(w-1)个移动窗口。在每个窗口位置m处,累计偏心距离∑εm、加权偏心夹角平均值θm以及偏心夹角的加权标准偏差σ(θm)可分别由下式计算:
利用这些数值,在每一个窗口位置处,计算可靠性权重因子:rm=∑εm/σ(εm)。在偏心距离较大的位置给予较大的权重因子,同时限制σ(θm)较大的位置。通过下式计算出一个校正偏心角θcorr:
并将其应用到序列的所有图像中,从而获得各帧图像的空间方位。
(5)血管腔内外表面的拟合
在对IVUS图像序列完成边缘提取并确定各帧的空间位置后,本发明采用NURBS(非均匀有理B样条)曲面拟合沿三维回撤路径正确排列的各横截面上的采样点,得到连续的三维血管表面。
(6)交互式冠状动脉虚拟内窥镜***
利用虚拟现实造型语言来显示内镜漫游模式的重建结果,不仅可显示重建后血管段的整体外观,而且可显示长轴纵剖面图像。
包括漫游路径的计算、重建出的管腔内外表面的绘制、虚拟场景中IVUS图像数据的显示以及交互式的用户图形接口的开发。
(6.1)漫游路径的计算
漫游路径是在目标血管腔内的一系列视点组成的序列。对于每个视点,都需要确定其位置和方向,其中方向采用观察者所在局部坐标系绕任意轴的旋转表示。
虚拟观察者在目标血管内的漫游既可沿IVUS导管的回撤路径进行,也可沿管腔中心线进行。故 可根据两者的三维重建结果直接得到:对于前者而言,就是沿导管的第i帧IVUS图像采集点的坐标,也即该帧图像的中心坐标;对于后者,就是从第i帧IVUS图像中分割出的管腔轮廓的重心坐标。由于导管的刚性和连续性都比计算出的管腔中心线要好,因此本发明采用第一种方法确定从而得到更为连续光滑的动画效果。
视点方向的确定:视点方向的初始值设定为沿z轴的负半轴,如附图6所示,其中 表示z轴负半轴的方向矢量,也即初始视点方向。由于漫游是沿导管的回撤路径进行的,因此根据上述视点位置的确定方法,由当前视点Pi可知下一个视点Pi+1的位置矢量 从而得到向量 其单位矢量即为待求的视点方向 与的叉积即为旋转轴:
(6.2)VRML中IVUS像素数据的显示
本发明对虚拟场景中***的IVUS像素数据采用半透明的显示方式,即一帧IVUS图像中各像素的透明度值不是同一个常数,而是取决于像素在图像中的位置和其灰度值。超声图像中除了血管壁和斑块以外,其它结构在虚拟内镜场景中都应该是不可见的,此时利用前面对超声图像的二维分割结果,将表示管腔和外膜以外回声信号的像素设置为全透明,允许漫游路径无阻挡地穿越这些区域。而对于管壁和斑块这些感兴趣区域,其透明度值取决于像素的灰度值,例如,亮回声信号表示可能存在的斑块,因此其不透明度值应设置为较高的数值;图像中的暗区可能表示其它血管的管腔或者没有产生回声的其它结构,其不透明度值应设为较低的数值。
(6.3)交互式的用户图形接口
本发明在VRML环境中设计开发出一个简明清晰、方便灵活的用户控制面板,使其能够完成以下功能:①用户可随时开启和关闭控制面板,并且开启时,以尽量不遮挡目标场景为原则。②虚拟观察者沿漫游路径前进时,在某个视点处,用户可在不同的显示模式之间进行切换,例如:按照正确的方向和位置显示在该点获取的IVUS图像;或者仅显示该点处的管腔表面(可同时开启或关闭半透明的IVUS帧);或者显示完成了彩色编码的管腔表面,其中彩色编码表示量化测量结果等。同时,用户可以随时进入或退出虚拟内镜的观察模式,显示重建后血管段的整体外观,或者长轴纵切面图像。③可调整漫游的速度和方向,虚拟观察者可在管腔内的任意位置停留。
Claims (3)
1.一种交互式冠状动脉虚拟血管镜的实现方法,其特征是,它通过将由近似正交的X射线冠状动脉造影CAG图像获得的管腔三维几何形态信息与由血管内超声获得的管腔横截面数据相融合,得到血管的三维模型,然后运用虚拟现实建模语言VRML交互地描述血管模型,实现内镜漫游模式的冠状动脉可视化,具体步骤如下:
A、利用同时采集的血管内超声IVUS图像和X射线冠状动脉造影CAG图像建立血管的三维模型;
所述血管内超声IVUS图像是在将机械式超声导管探头置于感兴趣血管段的远端,匀速等距回撤导引钢丝的过程中,利用血管内超声成像仪以ECG门控的方式在相同的心脏相位处采集到的等距的血管内超声IVUS图像序列,所述X射线冠状动脉造影CAG图像是利用C型臂单面X射线血管造影机在导管回撤路径的起点拍摄记录相同心脏状态的两幅近似垂直角度的X射线冠状动脉造影CAG图像;
所述建立血管三维模型的具体步骤如下:
a、根据两个成近似垂直角度的X射线冠状动脉造影CAG图像,重建出三维超声导管回撤路径;
b、从X射线冠状动脉造影CAG图像中重建出三维血管管腔:
将重建出的三维超声导管回撤路径向左右两个X射线冠状动脉造影CAG图像的成像平面反投影,得到对应的两维超声导管回撤路径,对于两维超声导管回撤路径上的每个点,通过在垂直于两维超声导管回撤路径的方向上,寻找灰度梯度的两个极大值,得到血管管腔左右边缘,然后在假设管腔的横截面是椭圆的前提下,完成整个三维血管管腔的重建,所述三维血管管腔用于后续确定各帧血管内超声IVUS图像的空间方向;
c、血管内超声IVUS图像序列中血管壁轮廓的提取:
在首帧血管内超声IVUS图像中手动选择血管壁内、外膜轮廓上的几个特征点,以连接这些特征点所形成的多边形作为初始位置,通过snake变形获得血管壁内、外膜的轮廓,分割出血管壁和可能存在的斑块,对于后续帧,则将前一帧的提取结果作为snake的初始位置,完成对连续多帧血管内超声IVUS图像的分割;
d、确定各帧血管内超声IVUS图像的轴向位置:
按照血管内超声IVUS图像的采集顺序和间距,沿重建出的三维超声导管回撤路径将各帧血管内超声IVUS图像顺序排列,确定出各帧血管内超声IVUS图像的轴向位置;
e、确定各帧血管内超声IVUS图像的空间方位:
在重建后的三维超声导管回撤路径上建立各帧血管内超声IVUS图像的局部坐标系,即Frenet-Serret标架,三个坐标轴分别为单位切矢量t、单位主法矢量n和单位副法矢量b,导管的位置位于血管内超声IVUS图像的中心;
将各帧血管内超声IVUS图像绕导管旋转至其正确的方向以确定血管内超声IVUS图像的空间方位:
用ρ表示血管壁轮廓的重心偏离导管位置的离心向量,把从X射线冠状动脉造影CAG图像重建出的血管管腔的椭圆轮廓投影到对应的血管内超声IVUS图像上,用μ来表示血管管腔椭圆轮廓中心线偏离导管位置的离心向量,ε为向量ρ的模,θ为ρ与μ的夹角,用统计优化方法确定血管内超声IVUS图像序列的空间方位,目的是使θ最小:
设定一个固定宽度w的移动窗口,在该窗口中进行统计分析,对于N帧组成的血管内超声IVUS图像序列,存在nw=N-(w-1)个移动窗口,在每个窗口位置m处,累计偏心距离∑εm、加权偏心夹角平均值以及偏心夹角的加权标准偏差σ(θm)分别由下式计算:
利用这些数值,在每一个窗口位置处,计算可靠性权重因子:rm=∑εm/σ(εm),在偏心距离较大的位置给予较大的权重因子,同时限制σ(θm)较大的位置,通过下式计算出一个校正偏心角
并将其应用到序列的所有血管内超声IVUS图像中,从而获得各帧血管内超声IVUS图像的空间方位;
f、利用基于NURBS曲面拟合的表面提取法完成血管表面的绘制;
B、运用虚拟现实建模语言VRML实现冠状动脉血管重建结果的内镜漫游模式的可视化:
具体步骤如下:
①漫游路径的计算:
根据三维超声导管回撤路径的重建结果,以沿导管的第i帧血管内超声IVUS图像采集点的坐标Pi为当前视点位置,Pi+1为下一个视点位置,以z轴负半轴的方向矢量为初始视点方向,则在VRML中,当前视点的位置矢量为旋转轴为当前视点方向为的单位矢量旋转轴为旋转角为
②将血管内超声IVUS图像像素数据***到虚拟场景中,并采用半透明的显示方式显示;
③开发交互式的用户图形接口。
2.根据权利要求1所述交互式冠状动脉虚拟血管镜的实现方法,其特征是,所述三维超声导管回撤路径的重建方法是:首先建立X射线冠状动脉造影CAG图像***在两个近似垂直角度的透视投影成像模型,再根据在造影过程中同步记录的距离和角度参数,得到成像***的几何变换矩阵,然后采用三维snake模型技术,snake直接在空间中变形,完成三维超声导管回撤路径的重建。
3.根据权利要求2所述交互式冠状动脉虚拟血管镜的实现方法,其特征是,所述X射线冠状动脉造影CAG图像的两个采集角度之间夹角的取值范围为60°至120°,且仅在三维超声导管回撤路径的起点拍摄一对造影图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100800020A CN101422352B (zh) | 2008-12-10 | 2008-12-10 | 一种交互式冠状动脉虚拟血管镜的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100800020A CN101422352B (zh) | 2008-12-10 | 2008-12-10 | 一种交互式冠状动脉虚拟血管镜的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101422352A CN101422352A (zh) | 2009-05-06 |
CN101422352B true CN101422352B (zh) | 2011-07-13 |
Family
ID=40613397
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100800020A Expired - Fee Related CN101422352B (zh) | 2008-12-10 | 2008-12-10 | 一种交互式冠状动脉虚拟血管镜的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101422352B (zh) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102289814B (zh) * | 2011-08-30 | 2012-12-26 | 北京理工大学 | 一种心脏核磁共振图像分割方法 |
CN102509267B (zh) * | 2011-11-08 | 2013-11-06 | 华北电力大学(保定) | 一种血管内超声图像序列的回顾性脱机门控方法 |
CN103247071B (zh) * | 2013-03-29 | 2015-11-11 | 哈尔滨工业大学深圳研究生院 | 一种构建三维血管模型方法及设备 |
KR20160107528A (ko) * | 2015-03-04 | 2016-09-19 | 삼성전자주식회사 | 컴퓨터 보조 진단을 위한 신뢰도 제공 장치 및 방법 |
CN106909770B (zh) * | 2015-12-21 | 2020-11-03 | 佳能株式会社 | 医疗图像处理装置、其控制方法和存储介质 |
CN114903591A (zh) * | 2016-03-21 | 2022-08-16 | 华盛顿大学 | 3d医学图像的虚拟现实或增强现实可视化 |
US11468111B2 (en) | 2016-06-01 | 2022-10-11 | Microsoft Technology Licensing, Llc | Online perspective search for 3D components |
CN107689072A (zh) * | 2016-06-12 | 2018-02-13 | 中慧医学成像有限公司 | 一种三维图像成像方法和*** |
CN108428210B (zh) * | 2017-02-15 | 2021-11-16 | 浙江京新术派医疗科技有限公司 | 血管图像重建方法和重建装置 |
CN107330236B (zh) * | 2017-05-11 | 2021-06-11 | 青岛大学附属医院 | 具有改良漫游效果的虚拟内窥镜*** |
CN107468334B (zh) * | 2017-08-01 | 2019-07-16 | 强联智创(北京)科技有限公司 | 一种三维微导管塑形辅助设计***及设计方法 |
DE102017214447B4 (de) * | 2017-08-18 | 2021-05-12 | Siemens Healthcare Gmbh | Planare Visualisierung von anatomischen Strukturen |
KR101910871B1 (ko) * | 2017-08-31 | 2018-12-31 | 주식회사 화승알앤에이 | 3차원 좌표를 이용한 튜브의 레이아웃을 해석하는 방법 및 기록매체 |
CN107945169B (zh) * | 2017-12-01 | 2022-02-15 | 中国人民解放军第三军医大学 | 一种冠状动脉影像分析方法 |
DE102018108643A1 (de) * | 2018-04-11 | 2019-11-14 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Positionsbestimmungsvorrichtung zum Bestimmen einer Position eines Gegenstands innerhalb einer tubulären Struktur |
CN109087352B (zh) * | 2018-08-16 | 2021-07-13 | 数坤(北京)网络科技股份有限公司 | 一种心脏冠脉优势型自动判别方法 |
CN109584148A (zh) * | 2018-11-27 | 2019-04-05 | 重庆爱奇艺智能科技有限公司 | 一种在vr设备中处理二维界面的方法和装置 |
CN109800814B (zh) * | 2019-01-25 | 2022-08-09 | 西南科技大学 | 叶片曲线测量定位的不变特征量提取方法 |
CN112057167B (zh) * | 2019-05-22 | 2022-05-24 | 深圳市德力凯医疗设备股份有限公司 | 一种血管手术超声导航方法以及超声导航设备 |
CN110458853B (zh) * | 2019-08-01 | 2021-02-26 | 北京灵医灵科技有限公司 | 一种医学影像中脚踝韧带分离方法和分离*** |
CN110522418A (zh) * | 2019-08-30 | 2019-12-03 | 深圳市中科微光医疗器械技术有限公司 | 一种心血管多模态融合分析方法、***、装置及存储介质 |
CN111009032B (zh) * | 2019-12-04 | 2023-09-19 | 浙江理工大学 | 基于改进外极线约束匹配的血管三维重建方法 |
CN112151180B (zh) * | 2019-12-05 | 2024-03-08 | 苏州润迈德医疗科技有限公司 | 具有狭窄病变的血管数学模型的合成方法和装置 |
CN111553979B (zh) * | 2020-05-26 | 2023-12-26 | 广州雪利昂生物科技有限公司 | 基于医学影像三维重建的手术辅助***及方法 |
CN112529976B (zh) * | 2020-11-26 | 2024-06-07 | 上海商汤智能科技有限公司 | 目标显示方法及装置、电子设备和存储介质 |
CN112652052B (zh) * | 2020-12-15 | 2022-07-22 | 山东大学 | 一种基于血管分支配准的冠状动脉三维重建方法及*** |
CN113409275B (zh) * | 2021-06-22 | 2022-07-01 | 青岛海信医疗设备股份有限公司 | 基于超声图像确定胎儿颈后透明层厚度的方法及相关装置 |
CN114145719B (zh) * | 2022-02-08 | 2022-04-26 | 天津恒宇医疗科技有限公司 | 双模冠脉血管图像三维融合的方法和融合*** |
-
2008
- 2008-12-10 CN CN2008100800020A patent/CN101422352B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101422352A (zh) | 2009-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101422352B (zh) | 一种交互式冠状动脉虚拟血管镜的实现方法 | |
CN101283929B (zh) | 一种血管三维模型的重建方法 | |
US6148095A (en) | Apparatus and method for determining three-dimensional representations of tortuous vessels | |
Roelandt et al. | Three-dimensional reconstruction of intracoronary ultrasound images. Rationale, approaches, problems, and directions. | |
JP5474342B2 (ja) | 3−d画像および表面マッピングによる解剖学的モデル化 | |
Pandian et al. | Dynamic three‐dimensional echocardiography: Methods and clinical potential | |
EP1088515B1 (en) | Vascular reconstruction | |
US6251072B1 (en) | Semi-automated segmentation method for 3-dimensional ultrasound | |
US20100189337A1 (en) | Method for acquiring 3-dimensional images of coronary vessels, particularly of coronary veins | |
EP0884975A1 (en) | Catheter for obtaining three-dimensional reconstruction of a vascular lumen and wall, and method of use | |
JP2006516440A (ja) | 三次元物体の再構築法 | |
Zheng et al. | Sequential reconstruction of vessel skeletons from X-ray coronary angiographic sequences | |
Di Mario et al. | Three dimensional reconstruction of cross sectional intracoronary ultrasound: clinical or research tool? | |
CN114145719B (zh) | 双模冠脉血管图像三维融合的方法和融合*** | |
Weng et al. | Three-dimensional surface reconstruction using optical flow for medical imaging | |
CN112669449A (zh) | 基于3d重建技术的cag和ivus精准联动分析方法及*** | |
Wahle et al. | Four-dimensional coronary morphology and computational hemodynamics | |
Weichert et al. | Virtual 3D IVUS vessel model for intravascular brachytherapy planning. I. 3D segmentation, reconstruction, and visualization of coronary artery architecture and orientation | |
Lengyel et al. | Three-dimensional reconstruction and volume rendering of intravascular ultrasound slices imaged on a curved arterial path | |
Zheng et al. | Reconstruction of coronary vessels from intravascular ultrasound image sequences based on compensation of the in-plane motion | |
Rotger et al. | Internal and external coronary vessel images registration | |
Garcia | Three-dimensional imaging for coronary interventions | |
Chandran et al. | Coronary arteries: imaging, reconstruction, and fluid dynamic analysis | |
Zheng | Three-dimensional reconstruction of vessels from intravascular ultrasound sequence and X-ray angiograms | |
US20240070855A1 (en) | Method and System of Calculating the Cross-section Area and Included Angle of Three-dimensional Blood Vessel Branch |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110713 Termination date: 20141210 |
|
EXPY | Termination of patent right or utility model |