CN101283929A - 一种血管三维模型的重建方法 - Google Patents
一种血管三维模型的重建方法 Download PDFInfo
- Publication number
- CN101283929A CN101283929A CNA2008100550379A CN200810055037A CN101283929A CN 101283929 A CN101283929 A CN 101283929A CN A2008100550379 A CNA2008100550379 A CN A2008100550379A CN 200810055037 A CN200810055037 A CN 200810055037A CN 101283929 A CN101283929 A CN 101283929A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- image
- sigma
- path
- ivus
- 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.)
- Granted
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
一种血管三维模型的重建方法,属医学检测技术领域,用于解决血管重建中的精度问题。其技术方案是:它将由IVUS图像序列获取的血管横截面信息和由基于X射线造影图像的三维重建获得的超声导管空间几何信息结合起来,准确重建出血管的解剖结构。与分别单独利用两种图像获得的重建结果相比较,本发明能够更全面和准确地反映血管及可能存在的斑块的真实形态,从而为冠心病及其它血管病变的临床诊治提供更为可靠的依据。
Description
技术领域
本发明涉及一种基于X射线血管造影和血管内超声图像的血管三维重建方法,属医学检测技术领域。
背景技术
X射线冠状动脉造影(coronary angiography,CAG)检查是目前国际上公认的诊断冠心病的“金标准”。具体实施方法是:用特制的导管经皮穿刺从股动脉或上臂的肱动脉,***导管,在X射线投影图像的指导下推送至主动脉根部置于冠状动脉口,然后推注造影剂,进行电影摄影或连续摄片。重复数次不同角度的拍摄,就能发现由粥样硬化引起的狭窄性病变,根据静止或动态观察造影剂的充盈和消失情况来判断血管解剖学形态异常的部位、性质和程度。该技术不仅能诊断缺血性心脏病及冠状动脉畸形等疾病,而且对冠状动脉内溶栓、PTCA(经皮腔内冠状动脉成形术)等介入手术治疗具有重要意义。但是CAG只能反映血管腔被造影剂充填后的投影轮廓,因而不仅会存在盲区,而且不能提供血管壁的结构信息和病变程度。
血管内超声(intravascular ultrasound,IVUS)是近年来应用于临床诊断血管病变的一项新技术。IVUS和CAG成像是同时进行的:经皮穿刺从股动脉或肱动脉***导管,向目标血管内注入***后,在X射线透视图像的指导下将顶端带有高频超声探头的0.014英寸导引钢丝穿越病变部位,到达血管远端,将超声探头与超声成像仪连接去除伪影后,经马达控制匀回撤导引钢丝,并记录图像。该技术可清晰显示血管横断面,根据斑块声学特征进行组织学分型,发现CAG不能显示的血管病变,观察分叉处或血管重叠处的模糊病变等。但是由于采用高频探头,影响了探测深度,只能对某一段病变血管进行测量,不能进入严重狭窄的管腔,并且无法确定截面的轴向位置和空间方向。
与二维图像相比,医学图像的三维重建和可视化可以提供更为丰富的医疗信息。采用近似正交角度的CAG图像对进行冠状动脉的三维重建,不仅能为医生提供形象、直观的三维血管图像,而且可以对血管的有关参数(如直径、长度、曲率等)进行定量测量。但是,该重建结果仅描述了血管内腔的空间位置和形态,无法得到管腔截面的形态和硬化斑块的弥漫程度及组成。且重建时一般假设管腔横截面为椭圆,事实上冠状动脉发生狭窄时管腔形状复杂多样,狭窄多呈偏心型和不规则型,因此这种假设是不准确的,很难保证三维重建的精度。
与CAG相比,IVUS有其独特的图像方位,可以显示管腔的横断面、血管壁的形态结构、斑块的组织形态学特征等,因此IVUS图像的三维重建有利于评价沿血管长轴方向管腔的边界、截面轮廓以及斑块的几何形态等。目前常用的IVUS三维重建方法是把一系列的超声图像按采集顺序叠加起来形成一个三维直血管段。由于IVUS本身不能提供每一帧图像的空间几何信息,因此这种方法没有考虑在图像获取过程中导管的弯曲和扭曲,其结果也是不准确的。
总之,目前还没有一种更为精确的血管三维重建方法,能够克服IVUS和CAG各自显示血管形态的局限性,准确重建出血管的解剖结构和反映血管的真实弯曲和扭曲,得到病变的准确位置和形态。
发明内容
本发明的目的是克服现有技术的不足、提供一种既能描述血管内腔的空间位置和形态,又能得到管腔截面的形态以及病变的准确位置和形态的血管三维模型的重建方法。
本发明所称问题是以下述技术方案实现的:
一种血管三维模型的重建方法,它将由IVUS图像获取的血管横截面信息和由基于X射线造影图像的三维重建获得的超声导管空间几何信息结合起来,准确重建出血管的解剖结构,具体重建步骤如下:
a、同时采集感兴趣血管段的IVUS和CAG图像:
将机械式超声导管探头置于感兴趣血管段的远端,在匀速等距地回撤导引钢丝过程中,利用血管内超声成像仪以ECG门控的方式在相同的心脏相位处采集等距的IVUS图像序列。同时,利用C型臂单面X射线血管造影机在导管回撤路径的起点拍摄记录相同心脏状态的两幅近似垂直角度的CAG图像;
b、根据上述两幅近似垂直角度的CAG图像,三维重建出超声导管的回撤路径;
c、在管腔截面是椭圆的假设前提下,从CAG图像对中三维重建出管腔,为后续确定各帧IVUS图像的空间方位所用:
将3D导管路径向两个成像平面反投影,得到对应的2D路径,对于2D路径上的每个点,通过在垂直于路径的方向上,寻找灰度梯度的两个极大值,得到血管管腔投影的左右边缘,然后在假设管腔的横截面是椭圆的前提下,完成整个管腔的三维重建;
d、各帧IVUS图像中血管壁轮廓的提取:
采用snake模型与动态规划相结合的方法完成各帧超声图像中血管壁内外膜轮廓的提取。首先,在首帧IVUS图像中手动选择血管壁内膜和外膜轮廓上的几个点,以连接这些点所形成的多边形作为snake初始位置。然后通过snake变形获得血管壁内外膜的轮廓,分割出血管壁和可能存在的斑块。对于后续帧,则将前一帧的提取结果作为snake的初始位置,完成对IVUS序列中连续多帧图像的分割;
e、确定各帧IVUS图像的轴向位置:
按照IVUS图像的采集顺序和间距,沿重建出的3D导管回撤路径将IVUS图像顺序排列,确定出各帧图像的轴向位置;
f、确定各帧IVUS图像的空间方位:
在重建后的3D导管路径上建立各帧超声图像的Frenet-Serret标架,即局部坐标系,三个坐标轴分别为单位切矢量t、单位主法矢量n和单位副法矢量b,导管的位置位于IVUS图像的中心。用ρ表示血管壁轮廓的重心偏离导管位置的偏心向量,把从CAG图像对中重建出的管腔的椭圆截面轮廓投影到对应的超声图像上,用μ来表示椭圆轮廓中心偏离导管位置的偏心向量,ε为向量ρ的模,θ为ρ与μ的夹角,用统计优化方法使θ最小,确定超声图像序列的空间方位:
设定一个宽度为w的移动窗口,在该窗口中进行统计分析,对于N帧组成的超声图像序列,存在nw=N-(w-1)个移动窗口,在每个窗口位置m处,累计偏心距离∑εm、加权偏心夹角平均值θm以及偏心夹角的加权标准偏差σ(θm)可分别由下式计算:
利用这些数值,在每一个窗口位置处,计算可靠性权重因子:rm=∑εm/σ(εm),在偏心距离较大的位置给予较大的权重因子,同时限制σ(θm)较大的位置,通过下式计算出一个校正偏心角θcorr:
并将其应用到IVUS序列的所有图像中,从而获得各帧图像的空间方位;
g、利用基于NURBS曲面拟合的表面提取法完成血管表面的绘制。
上述血管三维模型的重建方法,所述超声导管回撤路径的三维重建方法是:首先建立CAG***在两个近似垂直角度的透视投影成像模型,并推导出成像***的几何变换矩阵,再对手动选取的导管路径上的采样点进行三维重建,并以连接各重建点所得折线作为初始位置,通过snake变形完成导管路径的三维重建。
上述血管三维模型的重建方法,为了消除心脏的周期性运动和呼吸的影响,获得对应于同一时刻的图像序列,应采用ECG(心电)门控的方法,在相同的心脏相位处采集IVUS图像。
上述血管三维模型的重建方法,所述冠状动脉造影图像序列的两个采集角度之间夹角的取值范围为60°至120°。
上述血管三维模型的重建方法,在超声导管回撤路径的三维重建过程中,手动选取的导管路径上的采样点包括起点、终点和3~6个中间点。
本发明将由两近似正交角度的单面造影图像得到的血管空间几何信息与由血管内超声图像获得的管腔横截面信息结合起来,充分利用两种成像手段的互补性,完成了血管的准确三维重建,该方法克服了IVUS和CAG各自显示冠状动脉形态的局限性,能准确重建出血管的解剖结构和反映血管的真实弯曲和扭曲,得到病变的准确位置和形态。与分别单独利用两种图像获得的重建结果相比较,本发明能够更全面和准确地反映血管的真实形态,从而为冠心病及其它血管病变的临床诊治提供更为可靠的依据。
附图说明
下面结合附图对本发明作进一步详述。
图1是本发明的流程图;
图2是本发明的CAG和IVUS图像采集示意图;
图3是本发明的造影***在两个角度的成像示意图;
图4是本发明的超声图像序列中相邻帧间相对方位的确定示意图;
图5是本发明的超声图像偏心距离和偏心夹角示意图。
图中各符号为: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导管路径的空间参数曲线。
文中所用符号:t、单位切矢量;n、单位主法矢量;b、单位副法矢量;ρ、超声图像中血管壁轮廓的重心偏离导管位置的偏心向量;μ、管腔椭圆轮廓中心偏离导管位置的偏心向量;ε、向量ρ的模;θ、ρ与μ的夹角;w、移动窗口宽度;θm、加权偏心夹角平均值;σ(θm)、偏心夹角的加权标准偏差;rm、可靠性权重因子;θcprr、校正偏心角。
具体实施方式
下面结合附图详细说明本发明的步骤:
(1)图像采集:
采集设备包括C型臂单面X射线血管造影机和血管内超声成像仪。
IVUS和CAG成像是同时进行的(附图2)。在X射线透视图像的指导下将顶端带有高频超声探头的导引钢丝穿越病变部位,到达血管远端,将超声探头与超声成像仪连接去除伪影后,经马达控制匀速等距的回撤导引钢丝,并记录图像。成像过程中记录造影角度和X射线源焦点至接收屏的距离。
对于IVUS图像的采集,由于心脏的周期性运动和呼吸的影响,很难获得对应于同一时刻的图像序列。因此本发明采用ECG(心电)门控的方法,在相同的心脏相位处采集图像。
对于CAG图像的采集,传统方法一般是在心导管回撤过程中,采用X射线血管造影持续监视导管尖端的运动。由于会增加X射线的剂量,延长医患暴露在射线下的时间,因而该方法的临床应用价值不高。本发明采用仅在导管回撤路径的起点拍摄一对造影图像的方法。由于本发明采用机械式超声导管探头,超声换能器位于一可弯曲的轴心头端,轴心在外鞘管内旋转,而鞘管是固定不动的。因此可保证在回撤导管的后续过程中,带有超声换能器的导管尖端的回撤轨迹不会偏离在回撤路径起点拍摄的造影图像中显示的导管影像。
对于单面造影***,两个角度的图像不是同时获取的,给空间信息的提取带来很大困难。本发明采用同步记录的ECG信号选取记录相同心脏状态的两个角度的图像,从而解决两个角度的单面造影图像对的匹配问题。
(2)造影图像中导管路径和管腔边缘的提取和三维重建:
本发明首先建立CAG***在两个近似垂直角度的透视投影成像模型(附图3)。之后,根据在造影过程中同步记录的距离和角度参数,得到成像***的几何变换矩阵。然后利用三维snake模型技术,snake的初始位置采用手动取点获得,即在导管的一个投影上手动选取若干采样点(一般选取回撤路径的起点、终点和3~6个中间点即可),然后根据外极约束得到这些点在另一投影上的对应点。由这几组对应点分别求出它们的三维坐标,用直线段连接这些3D点,所得折线作为3D snake的初始位置。之后,表示导管的snake曲线在内外力的共同左右下直接在空间中变形,完成导管的三维重建。其中内力保证曲线的连续和光滑,外力包括两部分,分别对应于左右投影,保证三维曲线在两个角度成像平面上的投影恰好位于对应的导管投影处。该方法避免了基于外极约束的两个角度间的逐点匹配,提高了重建精度和运算速度。
按照成像***的几何变换矩阵,将3D导管路径向两个成像平面反投影,得到对应的2D路径。对于2D路径上的每个点,通过在垂直于路径的方向上,寻找灰度梯度的两个极大值,完成对血管管腔左右边缘的提取。之后,在假设管腔的横截面是椭圆的前提下,完成整个管腔的三维重建,为确定各帧IVUS图像的空间方位所用。
(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图像的中心,分割出的目标轮廓的重心一般不与导管位置重合,如附图5所示,其中C点表示导管,Oc为椭圆轮廓的中心(即在假设血管横截面是椭圆时,基于CAG的三维重建中所对应的血管中心线位置),OI为从超声图像中提取出的血管截面轮廓的重心。采用偏心向量ρ表示轮廓的重心偏离导管位置的程度:ρ=OI-C。
由于血管中心线和导管路径不重合,在血管同一位置处的超声图像轮廓和基于造影图像重建出的椭圆轮廓方位不一致,把椭圆轮廓投影到对应的超声图像上。同样采用椭圆轮廓的偏心向量μ来表示血管中心线偏离导管位置的程度:μ=OC-C。
超声图像的匹配误差可用向量ρ的模ε和ρ与μ的夹角θ表示。本发明利用统计优化方法使椭圆轮廓和超声轮廓的离心向量间的夹角θ最小,确定超声图像序列的绝对方位。设定一个宽度为w的移动窗口,在该窗口中进行统计分析。对于N帧组成的超声图像序列,存在nw=N-(w-1)个移动窗口。在每个窗口位置m处,累计偏心距离∑εm、加权偏心夹角平均值θm以及偏心夹角的加权标准偏差σ(θm)可分别由下式计算:
利用这些数值,在每一个窗口位置处,计算可靠性权重因子:rm=∑εm/σ(εm)。在偏心距离较大的位置给予较大的权重因子,同时限制σ(θm)较大的位置。通过下式计算出一个校正偏心角θcorr:
并将其应用到IVUS序列的所有图像中,从而获得各帧超声图像的空间方位。
(5)融合结果的显示
在对IVUS图像序列完成边缘提取并确定各帧的空间位置后,利用基于NURBS曲面拟合的表面提取法完成血管表面的绘制。此后,利用虚拟现实造型语言来显示重建结果,不仅可显示重建后血管段的整体外观,而且可显示长轴纵剖面图像。
Claims (4)
1、一种血管三维模型的重建方法,其特征是,它将由IVUS图像获取的血管横截面信息和由基于X射线造影图像的三维重建获得的超声导管空间几何信息结合起来,准确重建出血管的解剖结构,具体重建步骤如下:
a、同时采集感兴趣血管段的IVUS和CAG图像:
将机械式超声导管探头置于感兴趣血管段的远端,在匀速等距地回撤导引钢丝过程中,利用血管内超声成像仪以ECG门控的方式在相同的心脏相位处采集等距的IVUS图像序列,同时,利用C型臂单面X射线血管造影机在导管回撤路径的起点拍摄记录相同心脏状态的两个近似垂直角度的CAG图像;
b、根据上述两个近似垂直角度的CAG图像,三维重建出超声导管的回撤路径;
c、在管腔截面是椭圆的假设前提下,从CAG图像对中三维重建出管腔,为后续确定各帧IVUS图像的空间方位所用:
将3D导管路径向两个成像平面反投影,得到对应的2D路径,对于2D路径上的每个点,通过在垂直于路径的方向上,寻找灰度梯度的两个极大值,得到血管管腔投影的左右边缘,然后在假设管腔的横截面是椭圆的前提下,完成整个管腔的三维重建;
d、各帧IVUS图像中血管壁轮廓的提取:
采用snake模型与动态规划相结合的方法完成各帧超声图像中血管壁内外膜轮廓的提取,首先,在首帧IVUS图像中手动选择血管壁内膜和外膜轮廓上的几个点,以连接这些点所形成的多边形作为snake初始位置,然后通过snake变形获得血管壁内外膜的轮廓,分割出血管壁和可能存在的斑块,对于后续帧,则将前一帧的提取结果作为snake的初始位置,完成对IVUS序列中连续多帧图像的分割;
e、确定各帧IVUS图像的轴向位置:
按照IVUS图像的采集顺序和间距,沿重建出的3D导管回撤路径将IVUS图像顺序排列,确定出各帧图像的轴向位置;
f、确定各帧IVUS图像的空间方位:
在重建后的3D导管路径上建立各帧超声图像的Frenet-Serret标架,即局部坐标系,三个坐标轴分别为单位切矢量t、单位主法矢量n和单位副法矢量b,导管的位置位于IVUS图像的中心,用ρ表示血管壁轮廓的重心偏离导管位置的偏心向量,把从CAG图像对中重建出的管腔的椭圆截面轮廓投影到对应的超声图像上,用μ来表示椭圆轮廓中心偏离导管位置的偏心向量,ε为向量ρ的模,θ为ρ与μ的夹角,用统计优化方法使θ最小,确定超声图像序列的空间方位:
设定一个固定宽度w的移动窗口,在该窗口中进行统计分析,对于N帧组成的超声图像序列,存在nw=/N-(w-1)个移动窗口,在每个窗口位置m处,累计偏心距离∑εm、加权偏心夹角平均值θm以及偏心夹角的加权标准偏差σ(θm),可分别由下式计算:
利用这些数值,在每一个窗口位置处,计算可靠性权重因子:rm=∑εm/σ(εm),在偏心距离较大的位置给予较大的权重因子,同时限制σ(θm)较大的位置,通过下式计算出一个校正偏心角θcorr:
并将其应用到IVUS序列的所有图像中,从而获得各帧图像的绝对方位;
g、利用基于NURBS曲面拟合的表面提取法完成血管表面的绘制。
2、根据权利要求1所述的血管三维模型的重建方法,其特征是,所述超声导管回撤路径的三维重建方法是:首先建立CAG***在两个近似垂直角度的透视投影成像模型,并推导出成像***的几何变换矩阵,再对手动选取的导管路径上的采样点进行三维重建,并以连接各重建点所得折线作为初始位置,通过snake变形完成导管路径的三维重建。
3、根据权利要求1或2所述的血管三维模型的重建方法,其特征是,所述冠状动脉造影图像序列的两个采集角度之间夹角的取值范围为60°至120°。
4、根据权利要求2所述的血管三维模型的重建方法,其特征是,在超声导管回撤路径的三维重建过程中,手动选取的导管路径上的采样点包括起点、终点和3~6个中间点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100550379A CN101283929B (zh) | 2008-06-05 | 2008-06-05 | 一种血管三维模型的重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100550379A CN101283929B (zh) | 2008-06-05 | 2008-06-05 | 一种血管三维模型的重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101283929A true CN101283929A (zh) | 2008-10-15 |
CN101283929B CN101283929B (zh) | 2010-06-16 |
Family
ID=40056495
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100550379A Expired - Fee Related CN101283929B (zh) | 2008-06-05 | 2008-06-05 | 一种血管三维模型的重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101283929B (zh) |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101964118A (zh) * | 2010-09-30 | 2011-02-02 | 华北电力大学(保定) | 一种血管内超声图像序列的三维分割方法 |
CN102289814A (zh) * | 2011-08-30 | 2011-12-21 | 北京理工大学 | 一种心脏核磁共振图像分割方法 |
CN101625757B (zh) * | 2009-08-13 | 2012-02-22 | 华北电力大学(保定) | 一种抑制冠状动脉内超声图像序列中刚性运动伪影的方法 |
CN102393970A (zh) * | 2011-12-13 | 2012-03-28 | 北京航空航天大学 | 一种物体三维建模与渲染***及三维模型生成、渲染方法 |
CN102411795A (zh) * | 2011-08-16 | 2012-04-11 | 中国人民解放军第三军医大学第一附属医院 | 一种关于冠状动脉的多模图像的联合显示方法 |
CN101799864B (zh) * | 2010-01-15 | 2012-05-09 | 北京工业大学 | 一种基于血管内超声波图像的动脉斑块类型自动识别方法 |
CN102509292A (zh) * | 2011-11-03 | 2012-06-20 | 北京理工大学 | 一种心脏核磁共振图像的快速分割方法 |
CN102521873A (zh) * | 2011-11-22 | 2012-06-27 | 中国科学院深圳先进技术研究院 | 血管的建模方法 |
CN102846337A (zh) * | 2011-06-29 | 2013-01-02 | 清华大学 | 三维超声***及其目标点的定位方法和装置 |
CN102968783A (zh) * | 2012-10-15 | 2013-03-13 | 深圳市旭东数字医学影像技术有限公司 | 从腹部图像数据中自动分割骨骼的方法及*** |
CN101770650B (zh) * | 2009-01-07 | 2013-04-24 | 深圳迈瑞生物医疗电子股份有限公司 | 三维超声实时成像方法与装置以及成像*** |
CN103218797A (zh) * | 2012-01-19 | 2013-07-24 | 中国科学院上海生命科学研究院 | 血管图像处理与分析的方法和*** |
CN103462647A (zh) * | 2013-09-17 | 2013-12-25 | 浙江大学 | 一种集成血管内超声(ivus)影像和数字减影(dsa)影像的一体化在线实时处理仪 |
CN103593869A (zh) * | 2013-10-12 | 2014-02-19 | 沈阳东软医疗***有限公司 | 一种扫描设备及其图像显示方法 |
CN104644195A (zh) * | 2013-11-21 | 2015-05-27 | 西门子公司 | 用于执行旋转血管造影的检查对象的血管造影检查方法 |
CN104822323A (zh) * | 2012-12-03 | 2015-08-05 | 皇家飞利浦有限公司 | 超声和x射线模态的整合 |
CN105096270A (zh) * | 2015-08-07 | 2015-11-25 | 北京欣方悦医疗科技有限公司 | 一种冠脉三维重建中的钙化斑块去除方法 |
CN105095534A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 对血管中血流进行模拟的方法和*** |
CN105095615A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 对血管中血流数据进行处理的方法和*** |
CN105167849A (zh) * | 2015-05-19 | 2015-12-23 | 上海大学 | 一种基于蚁群算法的血管三维路径规划方法 |
CN105336004A (zh) * | 2014-06-27 | 2016-02-17 | 华为技术有限公司 | 一种曲面模型创建方法及装置 |
CN105640534A (zh) * | 2014-11-28 | 2016-06-08 | 韦伯斯特生物官能(以色列)有限公司 | 身体器官的差别标测 |
CN105938513A (zh) * | 2015-03-04 | 2016-09-14 | 三星电子株式会社 | 为计算机辅助诊断提供可靠性的装置和方法 |
CN106649913A (zh) * | 2016-09-08 | 2017-05-10 | 上海术理智能科技有限公司 | 一种基于多角度图像的3d打印设备及3d打印方法 |
CN107067398A (zh) * | 2017-04-28 | 2017-08-18 | 青岛海信医疗设备股份有限公司 | 用于三维医学模型中缺失血管的补全方法及装置 |
CN107392994A (zh) * | 2017-06-30 | 2017-11-24 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
CN107665737A (zh) * | 2017-01-23 | 2018-02-06 | 上海联影医疗科技有限公司 | 血管壁应力应变状态获取方法、计算机可读介质及*** |
CN108090901A (zh) * | 2017-12-28 | 2018-05-29 | 西安中科微光影像技术有限公司 | 一种基于心血管oct影像的生物支架对齐方法及装置 |
CN108430330A (zh) * | 2015-12-03 | 2018-08-21 | 西门子保健有限责任公司 | 产生血管***的体积图像的序列的断层造影设备和方法 |
CN108648231A (zh) * | 2018-05-14 | 2018-10-12 | 合肥融视信息科技有限公司 | 基于三维医学影像的管状结构长度测量***及方法 |
CN108961224A (zh) * | 2018-06-20 | 2018-12-07 | 厦门大学 | 一种基于隐式拉伸曲面的冠脉血管几何模型重建方法 |
CN110494094A (zh) * | 2017-02-06 | 2019-11-22 | 克里夫兰诊所基金会 | 表征解剖学结构的行为 |
CN110946651A (zh) * | 2013-08-13 | 2020-04-03 | 波士顿科学国际有限公司 | 解剖项的计算机可视化 |
CN111009032A (zh) * | 2019-12-04 | 2020-04-14 | 浙江理工大学 | 基于改进外极线约束匹配的血管三维重建方法 |
CN111507886A (zh) * | 2019-01-31 | 2020-08-07 | 许斐凯 | 利用超音波与深度学习技术的气管模型重建方法及*** |
CN112669449A (zh) * | 2020-12-31 | 2021-04-16 | 浙江理工大学 | 基于3d重建技术的cag和ivus精准联动分析方法及*** |
CN112711831A (zh) * | 2020-12-07 | 2021-04-27 | 上海联影医疗科技股份有限公司 | 血管仿真分析方法、设备、装置、计算机设备和存储介质 |
CN112914605A (zh) * | 2021-03-11 | 2021-06-08 | 中南大学湘雅二医院 | 一种血管3d构型定位方法及装置 |
CN113160254A (zh) * | 2020-02-26 | 2021-07-23 | 医慧(嘉兴)科技有限公司 | 基于三次b样条插值的血管内超声图像分割与测量方法 |
CN113192137A (zh) * | 2021-04-27 | 2021-07-30 | 北京阅影科技有限公司 | 光源矫正方法及装置 |
CN113499099A (zh) * | 2021-07-21 | 2021-10-15 | 上海市同仁医院 | 一种颈动脉超声自动扫查及斑块识别***及方法 |
CN114224484A (zh) * | 2021-12-13 | 2022-03-25 | 杭州脉流科技有限公司 | 颅内动脉瘤的表面区域定位方法、装置、计算机设备和存储介质 |
US20220167938A1 (en) * | 2015-09-16 | 2022-06-02 | Koninklijke Philips N.V. | Apparatus for vessel characterization |
CN116524134A (zh) * | 2023-07-04 | 2023-08-01 | 天津恒宇医疗科技有限公司 | 基于血管内图像的三维血管建模方法及ffr计算方法和*** |
WO2023151280A1 (zh) * | 2022-02-08 | 2023-08-17 | 天津恒宇医疗科技有限公司 | 双模冠脉血管图像三维融合的方法和融合*** |
CN117064446A (zh) * | 2023-10-13 | 2023-11-17 | 山东大学 | 基于血管内超声的血管动态三维重建*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE60328983D1 (de) * | 2002-06-04 | 2009-10-08 | Koninkl Philips Electronics Nv | Hybride dreidimensionale rekonstruktion der koronararterien mittels rotationsangiographie |
CN1871998A (zh) * | 2006-04-20 | 2006-12-06 | 北京思创贯宇科技开发有限公司 | 一种血管三维重建及虚拟支架置入的方法和*** |
-
2008
- 2008-06-05 CN CN2008100550379A patent/CN101283929B/zh not_active Expired - Fee Related
Cited By (70)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770650B (zh) * | 2009-01-07 | 2013-04-24 | 深圳迈瑞生物医疗电子股份有限公司 | 三维超声实时成像方法与装置以及成像*** |
CN101625757B (zh) * | 2009-08-13 | 2012-02-22 | 华北电力大学(保定) | 一种抑制冠状动脉内超声图像序列中刚性运动伪影的方法 |
CN101799864B (zh) * | 2010-01-15 | 2012-05-09 | 北京工业大学 | 一种基于血管内超声波图像的动脉斑块类型自动识别方法 |
CN101964118A (zh) * | 2010-09-30 | 2011-02-02 | 华北电力大学(保定) | 一种血管内超声图像序列的三维分割方法 |
CN101964118B (zh) * | 2010-09-30 | 2012-09-19 | 华北电力大学(保定) | 一种血管内超声图像序列的三维分割方法 |
CN102846337B (zh) * | 2011-06-29 | 2015-09-02 | 清华大学 | 三维超声***及其目标点的定位方法和装置 |
CN102846337A (zh) * | 2011-06-29 | 2013-01-02 | 清华大学 | 三维超声***及其目标点的定位方法和装置 |
CN102411795A (zh) * | 2011-08-16 | 2012-04-11 | 中国人民解放军第三军医大学第一附属医院 | 一种关于冠状动脉的多模图像的联合显示方法 |
CN102411795B (zh) * | 2011-08-16 | 2014-12-10 | 中国人民解放军第三军医大学第一附属医院 | 一种关于冠状动脉的多模图像的联合显示方法 |
CN102289814A (zh) * | 2011-08-30 | 2011-12-21 | 北京理工大学 | 一种心脏核磁共振图像分割方法 |
CN102289814B (zh) * | 2011-08-30 | 2012-12-26 | 北京理工大学 | 一种心脏核磁共振图像分割方法 |
CN102509292A (zh) * | 2011-11-03 | 2012-06-20 | 北京理工大学 | 一种心脏核磁共振图像的快速分割方法 |
CN102509292B (zh) * | 2011-11-03 | 2013-09-25 | 北京理工大学 | 一种心脏核磁共振图像的快速分割方法 |
CN102521873A (zh) * | 2011-11-22 | 2012-06-27 | 中国科学院深圳先进技术研究院 | 血管的建模方法 |
CN102393970A (zh) * | 2011-12-13 | 2012-03-28 | 北京航空航天大学 | 一种物体三维建模与渲染***及三维模型生成、渲染方法 |
CN103218797A (zh) * | 2012-01-19 | 2013-07-24 | 中国科学院上海生命科学研究院 | 血管图像处理与分析的方法和*** |
CN103218797B (zh) * | 2012-01-19 | 2016-01-27 | 中国科学院上海生命科学研究院 | 血管图像处理与分析的方法和*** |
CN102968783B (zh) * | 2012-10-15 | 2016-09-14 | 深圳市旭东数字医学影像技术有限公司 | 从腹部图像数据中自动分割骨骼的方法及*** |
CN102968783A (zh) * | 2012-10-15 | 2013-03-13 | 深圳市旭东数字医学影像技术有限公司 | 从腹部图像数据中自动分割骨骼的方法及*** |
US11213273B2 (en) | 2012-12-03 | 2022-01-04 | Koninklijke Philips N.V. | Integration of ultrasound and x-ray modalities |
CN104822323A (zh) * | 2012-12-03 | 2015-08-05 | 皇家飞利浦有限公司 | 超声和x射线模态的整合 |
CN104822323B (zh) * | 2012-12-03 | 2020-06-30 | 皇家飞利浦有限公司 | 超声和x射线模态的整合 |
CN110946651A (zh) * | 2013-08-13 | 2020-04-03 | 波士顿科学国际有限公司 | 解剖项的计算机可视化 |
CN103462647A (zh) * | 2013-09-17 | 2013-12-25 | 浙江大学 | 一种集成血管内超声(ivus)影像和数字减影(dsa)影像的一体化在线实时处理仪 |
CN103462647B (zh) * | 2013-09-17 | 2015-01-07 | 浙江大学 | 一种集成血管内超声(ivus)影像和数字减影(dsa)影像的一体化在线实时处理仪 |
CN103593869A (zh) * | 2013-10-12 | 2014-02-19 | 沈阳东软医疗***有限公司 | 一种扫描设备及其图像显示方法 |
CN103593869B (zh) * | 2013-10-12 | 2016-08-10 | 沈阳东软医疗***有限公司 | 一种扫描设备及其图像显示方法 |
CN104644195A (zh) * | 2013-11-21 | 2015-05-27 | 西门子公司 | 用于执行旋转血管造影的检查对象的血管造影检查方法 |
US9888895B2 (en) | 2013-11-21 | 2018-02-13 | Siemens Aktiengesellschaft | Angiographic examination method to implement a rotational angiography |
CN105095615A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 对血管中血流数据进行处理的方法和*** |
CN105095534A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 对血管中血流进行模拟的方法和*** |
CN105095615B (zh) * | 2014-04-23 | 2019-05-21 | 北京冠生云医疗技术有限公司 | 对血管中血流数据进行处理的方法和*** |
CN105336004A (zh) * | 2014-06-27 | 2016-02-17 | 华为技术有限公司 | 一种曲面模型创建方法及装置 |
CN105640534A (zh) * | 2014-11-28 | 2016-06-08 | 韦伯斯特生物官能(以色列)有限公司 | 身体器官的差别标测 |
CN105938513A (zh) * | 2015-03-04 | 2016-09-14 | 三星电子株式会社 | 为计算机辅助诊断提供可靠性的装置和方法 |
CN105167849B (zh) * | 2015-05-19 | 2017-07-25 | 上海大学 | 一种基于蚁群算法的血管三维路径规划方法 |
CN105167849A (zh) * | 2015-05-19 | 2015-12-23 | 上海大学 | 一种基于蚁群算法的血管三维路径规划方法 |
CN105096270B (zh) * | 2015-08-07 | 2018-04-06 | 北京欣方悦医疗科技有限公司 | 一种冠脉三维重建中的钙化斑块去除方法 |
CN105096270A (zh) * | 2015-08-07 | 2015-11-25 | 北京欣方悦医疗科技有限公司 | 一种冠脉三维重建中的钙化斑块去除方法 |
US20220167938A1 (en) * | 2015-09-16 | 2022-06-02 | Koninklijke Philips N.V. | Apparatus for vessel characterization |
CN108430330A (zh) * | 2015-12-03 | 2018-08-21 | 西门子保健有限责任公司 | 产生血管***的体积图像的序列的断层造影设备和方法 |
CN108430330B (zh) * | 2015-12-03 | 2021-10-01 | 西门子保健有限责任公司 | 产生血管***的体积图像的序列的断层造影设备和方法 |
CN106649913A (zh) * | 2016-09-08 | 2017-05-10 | 上海术理智能科技有限公司 | 一种基于多角度图像的3d打印设备及3d打印方法 |
CN107665737A (zh) * | 2017-01-23 | 2018-02-06 | 上海联影医疗科技有限公司 | 血管壁应力应变状态获取方法、计算机可读介质及*** |
CN110494094B (zh) * | 2017-02-06 | 2022-10-25 | 克里夫兰诊所基金会 | 表征解剖学结构的行为 |
CN110494094A (zh) * | 2017-02-06 | 2019-11-22 | 克里夫兰诊所基金会 | 表征解剖学结构的行为 |
CN107067398A (zh) * | 2017-04-28 | 2017-08-18 | 青岛海信医疗设备股份有限公司 | 用于三维医学模型中缺失血管的补全方法及装置 |
WO2019000479A1 (zh) * | 2017-06-30 | 2019-01-03 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
CN107392994A (zh) * | 2017-06-30 | 2017-11-24 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
CN107392994B (zh) * | 2017-06-30 | 2018-11-06 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
CN108090901B (zh) * | 2017-12-28 | 2021-11-19 | 中科微光医疗研究中心(西安)有限公司 | 一种基于心血管oct影像的生物支架对齐方法及装置 |
CN108090901A (zh) * | 2017-12-28 | 2018-05-29 | 西安中科微光影像技术有限公司 | 一种基于心血管oct影像的生物支架对齐方法及装置 |
CN108648231A (zh) * | 2018-05-14 | 2018-10-12 | 合肥融视信息科技有限公司 | 基于三维医学影像的管状结构长度测量***及方法 |
CN108961224A (zh) * | 2018-06-20 | 2018-12-07 | 厦门大学 | 一种基于隐式拉伸曲面的冠脉血管几何模型重建方法 |
CN111507886A (zh) * | 2019-01-31 | 2020-08-07 | 许斐凯 | 利用超音波与深度学习技术的气管模型重建方法及*** |
CN111009032A (zh) * | 2019-12-04 | 2020-04-14 | 浙江理工大学 | 基于改进外极线约束匹配的血管三维重建方法 |
CN111009032B (zh) * | 2019-12-04 | 2023-09-19 | 浙江理工大学 | 基于改进外极线约束匹配的血管三维重建方法 |
CN113160254A (zh) * | 2020-02-26 | 2021-07-23 | 医慧(嘉兴)科技有限公司 | 基于三次b样条插值的血管内超声图像分割与测量方法 |
CN112711831A (zh) * | 2020-12-07 | 2021-04-27 | 上海联影医疗科技股份有限公司 | 血管仿真分析方法、设备、装置、计算机设备和存储介质 |
CN112669449A (zh) * | 2020-12-31 | 2021-04-16 | 浙江理工大学 | 基于3d重建技术的cag和ivus精准联动分析方法及*** |
CN112914605A (zh) * | 2021-03-11 | 2021-06-08 | 中南大学湘雅二医院 | 一种血管3d构型定位方法及装置 |
CN113192137A (zh) * | 2021-04-27 | 2021-07-30 | 北京阅影科技有限公司 | 光源矫正方法及装置 |
CN113499099A (zh) * | 2021-07-21 | 2021-10-15 | 上海市同仁医院 | 一种颈动脉超声自动扫查及斑块识别***及方法 |
CN114224484A (zh) * | 2021-12-13 | 2022-03-25 | 杭州脉流科技有限公司 | 颅内动脉瘤的表面区域定位方法、装置、计算机设备和存储介质 |
CN114224484B (zh) * | 2021-12-13 | 2023-05-30 | 杭州脉流科技有限公司 | 颅内动脉瘤的表面区域定位方法、装置、计算机设备和存储介质 |
WO2023151280A1 (zh) * | 2022-02-08 | 2023-08-17 | 天津恒宇医疗科技有限公司 | 双模冠脉血管图像三维融合的方法和融合*** |
CN116524134A (zh) * | 2023-07-04 | 2023-08-01 | 天津恒宇医疗科技有限公司 | 基于血管内图像的三维血管建模方法及ffr计算方法和*** |
CN116524134B (zh) * | 2023-07-04 | 2023-08-29 | 天津恒宇医疗科技有限公司 | 基于血管内图像的三维血管建模方法及ffr计算方法和*** |
CN117064446A (zh) * | 2023-10-13 | 2023-11-17 | 山东大学 | 基于血管内超声的血管动态三维重建*** |
CN117064446B (zh) * | 2023-10-13 | 2024-01-12 | 山东大学 | 基于血管内超声的血管动态三维重建*** |
Also Published As
Publication number | Publication date |
---|---|
CN101283929B (zh) | 2010-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101283929B (zh) | 一种血管三维模型的重建方法 | |
CN101422352B (zh) | 一种交互式冠状动脉虚拟血管镜的实现方法 | |
CN107392994B (zh) | 冠脉血管的三维重建方法、装置、设备及存储介质 | |
US6148095A (en) | Apparatus and method for determining three-dimensional representations of tortuous vessels | |
US5771895A (en) | Catheter for obtaining three-dimensional reconstruction of a vascular lumen and wall | |
Roelandt et al. | Three-dimensional reconstruction of intracoronary ultrasound images. Rationale, approaches, problems, and directions. | |
JP5281740B2 (ja) | 管状臓器の3次元表現決定システム | |
CN112652052B (zh) | 一种基于血管分支配准的冠状动脉三维重建方法及*** | |
CN103077550B (zh) | 一种非门控icus图像序列中血管的四维重建方法 | |
US20100189337A1 (en) | Method for acquiring 3-dimensional images of coronary vessels, particularly of coronary veins | |
CN101953696B (zh) | 一种icus图像序列中血管的三维形态参数测量方法 | |
JP5536678B2 (ja) | 医用画像表示方法、医用画像診断装置、及び医用画像表示装置 | |
US20200352542A1 (en) | Guided-transcranial ultrasound imaging using neural networks and associated devices, systems, and methods | |
Di Mario et al. | Three dimensional reconstruction of cross sectional intracoronary ultrasound: clinical or research tool? | |
CN114145719B (zh) | 双模冠脉血管图像三维融合的方法和融合*** | |
CN103247071B (zh) | 一种构建三维血管模型方法及设备 | |
US12039685B2 (en) | Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device | |
CN113057677B (zh) | 超声和ct图像融合的心脏影像建模方法、***及设备 | |
Shekhar et al. | Fusion of intravascular ultrasound and biplane angiography for three-dimensional reconstruction of coronary arteries | |
Wahle et al. | Four-dimensional coronary morphology and computational hemodynamics | |
CN112669449A (zh) | 基于3d重建技术的cag和ivus精准联动分析方法及*** | |
Dijkstra et al. | Quantitative coronary ultrasound: State of the art | |
CN114664455A (zh) | 一种冠状动脉血流储备分数计算方法及装置 | |
Zheng et al. | Reconstruction of coronary vessels from intravascular ultrasound image sequences based on compensation of the in-plane motion | |
Chandran et al. | Coronary arteries: imaging, reconstruction, and fluid dynamic analysis |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100616 Termination date: 20140605 |