CN101594824A - 基于预确定的最优视图映射图对ra确定最优旋转轨迹 - Google Patents

基于预确定的最优视图映射图对ra确定最优旋转轨迹 Download PDF

Info

Publication number
CN101594824A
CN101594824A CNA2007800242615A CN200780024261A CN101594824A CN 101594824 A CN101594824 A CN 101594824A CN A2007800242615 A CNA2007800242615 A CN A2007800242615A CN 200780024261 A CN200780024261 A CN 200780024261A CN 101594824 A CN101594824 A CN 101594824A
Authority
CN
China
Prior art keywords
type arm
ray
optimal
interest
freedom
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
Application number
CNA2007800242615A
Other languages
English (en)
Other versions
CN101594824B (zh
Inventor
B·莫瓦萨吉
O·温克
S-Y·陈
J·A·加西亚
J·D·卡罗尔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
University of Colorado
Original Assignee
Koninklijke Philips Electronics NV
University of Colorado
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics NV, University of Colorado filed Critical Koninklijke Philips Electronics NV
Publication of CN101594824A publication Critical patent/CN101594824A/zh
Application granted granted Critical
Publication of CN101594824B publication Critical patent/CN101594824B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/44Constructional features of apparatus for radiation diagnosis
    • A61B6/4429Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4435Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
    • A61B6/4441Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/46Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/467Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50Clinical applications
    • A61B6/504Clinical applications involving diagnosis of blood vessels, e.g. by angiography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/027Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents

Abstract

描述了一种为C型臂X射线系统确定三维旋转X射线冠状动脉血管造影的最优轨迹(25)的方法。该C型臂X射线系统具有至少两个自由度。它们由以左/右冠状动脉斜角表示的C型臂(11)的旋转运动和以后侧角/前侧角表示的C型臂(11)的滚转运动所定义。该方法依次按顺序下述步骤。首先,执行感兴趣区域内体血管中心线的三维表示的生成。其次,生成至少一个最优视图映射图。最后,在最优视图映射图内确定X射线系统的最优轨迹(25),其中最优轨迹(25)至少由C型臂在其两个自由度内的运动所确定,该轨迹使得图像投影具有最小的投影缩短和/或重叠,同时使对X射线的暴露最小化。

Description

基于预确定的最优视图映射图对RA确定最优旋转轨迹
本发明涉及旋转X射线血管造影(RA)领域,并且更具体地涉及基于最优视图映射图(optimal view map,OVM)来确定RA的最优旋转轨迹。
在介入神经放射学中,神经放射学家或血管造影师在任何时候知道导管以毫米精度位于病人体内的位置是重要的。该信息从数字减影血管造影(DSA)图像推导出,由于他/她具有解剖知识,他/她心理上将该减影图像与术前三维图像(例如,磁共振(MR)图像)联系起来。
现今,3D X射线旋转血管造影(3D-RA)重建的体积常规地由旋转血管造影(RA)序列生成。尽管这样的体积已被证实为给医师带来实际的辅助性帮助,DSA仍然是介入神经放射学的基础。结果,具有3D-RA体积的DSA图像的配准似乎是非常有前途的特征。
当前的导管室,即冠状动脉疾病或其他血管疾病治疗中心,基于二维(2D)投影图像来执行诸如定性的狭窄确定、球囊扩张、支架治疗等的介入过程。近年来,旋转血管造影(RA)已被引入,其中C型臂X射线系统围绕病人旋转,同时从填充造影剂的冠状动脉采集投影图像。如J.T.Maddux、O.Wink、J.C.Messenger、B.M.Groves、R.Liao、J.Strzelczyk、S.Y.Chen、J.D.Carroll在已出版的2004年的Catheterization andCardiovascular Interventions中的“A Randomized Study of the Safety andClinical Utility of Rotational Angiography versus Standard Angiography in theDiagnosis of Coronary Artery Disease”中所描述的,这些数据集可以用于诊断。
如B.Movassaghi、V.Rasche、M.Grass、M.Viergever、W.Niessen在2004年的IEEE Trans.Med.Imag第12卷第23期第1517-1531页的“Aquantitative analysis of 3D coronary modeling from two or more projectionimages”中所描述的,该数据集还可以用于3D冠状动脉建模。
此外,如V.Rasche、A.Buecker、M.Grass、R.Suurmond、R.Koppe.、H.Kuehl在2003年的RSNA第83届科学年会第C19-382页的“ECG-gated3D Rotational CoronaryAngiography”中所述的,所采集的数据集还用于3D冠状动脉重建程序。上述出版物的主题被视为本申请的整体部分并且应以引用的方式将其包括。
基于医师的经验来选择当前临床应用的旋转采集方案并且该方案不受到任何科学背景的影响。典型地,医师将X射线系统定位在特定坐标以采集投影图像,将X射线系统定位在下一特定坐标以采集下一投影图像,等。因此,用于采集所述投影图像的位置确定由于人体的差异而不能被最优化。为了获得冠状动脉树的最精确的三维模型,应该使用C型臂血管造影系统的病人特定位置。
冠状动脉树的三维(3D)特性由于投影几何结构而导致任何投影中各种片段的投影缩短(foreshortening)。因此,各种旋转采集方案包括具有或多或少血管投影缩短和血管重叠的大量投影图像。如果感兴趣对象没有被定位为与X射线探测器的投影平面平行,而是如图1所示地与该平面成某一角度,则发生投影缩短。
医师基于经验从实验值中选择坐标以减少感兴趣区域(ROI)的各个冠状动脉树的投影缩短和重叠。实际上,得到多于理论上必需的图像以便对最终的冠状动脉树重建选择最佳图像投影。另一方面,将所采集的投影数量保持得较低以便将病人对X射线的暴露保持为最低限度对病人是有利的。另一方面,医师需要视图和图像投影以便重建冠状动脉树的最可能的三维模型以供其诊断和/或治疗。
S.James Chen和John D.Carrol在2000年4月第19卷第4期的IEEEtransaction on medical imaging中的“3D Reconstruction of Coronary ArterialTree to Optimize Angiographic Visualization”中提出一种基于对特定心脏相位从所确定的冠状动脉3D中心线得出的计算机生成的2D中心线模型,对每个任意投影角确定血管投影缩短和血管重叠的定量值的方法。该文献描述了由于血管重叠和投影缩短,多个投影对于用血管造影充分评估冠状动脉树是必要的。只有当成功解决这些可视化问题时,才可能最优地执行基于导管的介入。传统方法提供多个选择的视图,其中基于二维(2D)投影主观地使重叠和投影缩短最小化。根据常规血管造影研究利用单平面成像系统以任意取向采集的一对图像被选择用于三维(3D)重建。在选择感兴趣动脉片段(例如,单个冠状动脉狭窄或分叉病变)之后,计算一组使片段投影缩短最小化的C型臂成角。然后将多个具有最小化片段投影缩短的计算机生成的投影图像用来选择视图,该视图具有相对于感兴趣片段或区域(ROI)最小重叠的血管。然后利用最优视图或最优视图映射图来引导随后的血管造影采集和介入。
然后增强该方法以便生成合并有冠状动脉树的4D(3-D加时间)特性的完整最优视图映射图。医师可能使用如此生成的最优视图映射图(OVM)来选择具有最小的血管投影缩短和重叠的静态视图以用于介入过程。典型地,所生成的最优视图映射图是彩色的。明亮区域典型地指示具有最小的投影缩短和/或重叠的区域。在图2的示例中,白色区域示出了具有0-10%投影缩短或重叠的区域。因此,血管造影系统的这些位置表示相对于病人身体的最优位置从而对冠状动脉树或其部分的最佳三维重建生成图像投影。
可以基于Chen SYJ、Carroll JD、Messenger JC发表在2002年的IEEETrans.Med.Imag.第21期第724-740页的“Quantitative analysis ofreconstructed 3-D coronary arterial tree and intra-coronary devices”中描述的方法来确定血管重叠。具体血管片段Ck的重叠值被定义为相对所有其他动脉Ci的重叠传播(propagation of overlap):
Figure A20078002426100081
其中Ci和Ck分别表示第i动脉的动脉管腔和在第k动脉处所选择的感兴趣片段。∏α,β(C)表示算子,其基于C型臂角度(α,β)在将对象C投影到图像平面之后对像素数量进行记数。
然而,因为基于医师的经验来选择当前临床应用的旋转采集方案,所以确定最优的一组采集图像投影的问题仍存在。
因此需要一种确定C型臂X射线系统最优旋转运行的方法,以便对感兴趣区域采集具有最小血管投影缩短和最小血管重叠的投影图像。
根据本发明的第一方面,提供了一种为C型臂X射线系统针对血管状结构确定旋转X射线血管造影的最优轨迹的方法,该C型臂X射线系统具有至少两个自由度,所述自由度由以左/右冠状动脉斜角α表示的C型臂的螺旋桨式运动,和以后侧角/前侧角β表示的C型臂的滚转运动所定义。该方法包括以下步骤:
(a)生成感兴趣区域内的体血管中心线的三维表示;
(b)生成在x和y方向由α和β的最大值限制的最优视图映射图和/或选择由临床医生确定的优选视角;并且
(c)在最优视图映射图的限制内为X射线系统C型臂计算最优轨迹,其中最优轨迹至少由C型臂在其两个自由度内的运动所确定,该轨迹使得投影图像具有最小投影缩短和重叠,同时使感兴趣区域对X射线的暴露最小化。
所述方法的优点在于,在最小的血管投影缩短和血管重叠方面,自动确定C型臂X射线系统的旋转血管造影的最优轨迹,从而确保二维血管造影具有最优诊断值和/或最优三维重建图像质量。感兴趣区域典型地是病人的特定血管片段或完整冠状动脉树。任何其他体血管也可以是感兴趣区域。那些区域包括冠状动脉、头部内的冠状叶(coronary vane)或血管或结构。应该明确提及的是,病人体内所有血管均不应被排除在使用所要求保护的方法之外。该方法可以应用到活体的任何部分。
当前使用的轨迹或C型臂位置基于医师的经验,不能保证最优图像投影,并且因此不能向医师递交所需质量的图像。这限制了图像投影或重建的冠状动脉树或其部分对于医师的价值。
本发明的方法使用现代X射线血管造影系统,该系统可以同时以两个自由度移动血管造影系统的C型臂。这是必需的以便在C型臂力学的限制内实现X射线系统的任何可能轨迹。传统上,3D-RA系统(三维X射线血管造影系统)允许每次仅围绕一个角度的旋转。C型臂允许螺旋桨式运动,其中临床相关角度范围从120°LAO(左动脉倾斜)到120°RAO(右动脉倾斜)。此外,C型臂的滚转运动也是可能的。临床相关角度的范围从60°CRAN(前侧角)到60°CAUD(后侧角)。这相当于围绕C型臂的中立位相对于其旋转枢轴点或中心旋转+/-60°。
所要求保护的方法的步骤(b)允许两个可替代的选项:生成最优视图映射图,或临床医生具有他/她自己的能从其构造图像的一组最优视角。然后可以计算最优轨迹以便访问这些以视角表示的点。此外,许多临床医生可以对最优地满足他们需求的轨迹达成共识。
根据本发明的另一方面,提供了一种X射线C型臂系统,其具有运行根据用于对旋转X射线冠状动脉血管造影确定最优轨迹的方法所确定的轨迹的能力。
根据本发明的另一方面,提供了一种计算机系统,其用于根据用于对旋转X射线冠状动脉血管造影确定最优轨迹并控制C型臂X射线血管造影系统的运动的方法来计算轨迹。
本发明的另一方面是提供一种计算机程序产品,其具有用于根据用于对旋转X射线冠状动脉血管造影确定最优轨迹并控制C型臂X射线血管造影系统的运动的方法来计算轨迹的指令。
本发明的另一方面是提供一种计算机可读介质,其存储用于计算机系统并能够根据用于对旋转X射线冠状动脉血管造影确定最优轨迹并控制C型臂X射线血管造影系统的运动的方法来计算轨迹的一组指令。
上述方面的优点在于,使得RA系统全自动地运行用于提供具有最小可能的血管投影缩短和/或血管重叠的最优图像投影的轨迹。
在本发明的一个实施例中,确定最优轨迹的步骤包括使下列等式最小化的步骤:
F ( κ , λ , α , β ) = Σ α = - 120 α = 120 Σ β = - 60 β = 60 κf ( α , β ) + λO ( α , β ) - - - ( 2 )
其中:
κ是加权参数;
λ是加权参数;
α是左冠状动脉斜角或右冠状动脉斜角的角度值;
β是后侧角或前侧角的角度值;
f(α,β)是与血管投影缩短相关的函数;以及
O(α,β)是与血管重叠相关的函数。
κ和λ是典型地从0到1变化的参数,并且根据本发明的另一实施例,它们的数学和是1。
在一个实施例中,血管造影系统的C型臂安装到L型臂上,产生组合的L型臂/C型臂血管造影系统并能够实现一个额外的自由度,从而在相对于病人定位造影系统上给予血管造影师更大的灵活性。
在另一实施例中,生成感兴趣区域中的血管中心线的3D表示的步骤包括基于两个或更多个采集的投影图像的建模方法。如上所述,这些投影图像可以通过B.Movassaghi、V.Rasche、M.Grass、M.Viergever、W.Niessen在2004年的IEEE Trans.Med.Imag第12卷第23期第1517-1531页的“Aquantitative analysis of 3D coronary modeling from two or more projectionimages”中所描述的方法来采集。该文献的主题也应以引用的方式将其并入。
在本发明的另一实施例中,生成感兴趣区域的中心线的三维表示的步骤包括基于预采集的不同模态的图像数据集的建模方法。因此该实施例使用来自先前检查的预采集的数据集。两类数据集应被提及为示例:计算机断层摄影(CT)扫描和磁共振(MR)数据集。这些数据对于有心脏病的病人通常是可获得的并可被重复使用。
在又一实施例中,生成感兴趣区域内血管中心线的三维表示的步骤包括以基于非最优的任意采集轨迹的预采集的三维旋转血管造影图像为基础的建模方法。如上所述,可以基于医师的经验来采集这些图像。这些图像中的一些对于有规划地治疗病人来说可能足够好,但它们表示无法保证以下方面的点数据:即它们表示用于重建血管的三维模型的最佳可用数据。
还有另一方式来执行生成感兴趣区域内血管中心线的三维表示的步骤。在该实施例中,该步骤包括基于表示平均的人体冠状动脉的体模的建模方法。熟知的体模是包含从大约二十四个病人获得的经验的Dudge模型。该模型是对不同测试个人之间的变化进行平均的抽象模型。然而,该模型也不给出完善的体模。
必须注意,已经参考不同主题描述了本发明的各实施例。具体地,已经参考装置类型权利要求描述了一些实施例,而参考方法类型权利要求描述了其他实施例。然而,本领域的技术人员可从上面和下面的描述中获悉,除非特别指出,除属于一种主题类型的特征的任意组合外,涉及不同主题的特征之间的任意组合,具体是方法类型权利要求的特征和设备类型权利要求的特征之间的任意组合也被认为是由本申请所公开的。
基于用于为组合C型臂X射线或组合L型臂/C型臂系统确定旋转X射线冠状动脉血管造影的最优轨迹的方法的上面已给出的和下面的解释,技术人员能将该方法的所述步骤转换成执行该方法的计算机程序。
根据下文中要描述的实施例的示例,本发明的上面限定的方面和另外的方面将变得明显并参考实施例的示例进行解释。但本发明不限于这些示例。
附图的图示说明是示意性的。应该注意的是,在不同的附图中,类似的或同一元件具有相同的附图标记。附图示出了:
图1a是血管0%投影缩短的示意图;
图1b是血管30%投影缩短的示意图;
图1c是示出了三维视图中血管0%和30%投影缩短的示意图;
图2是特定血管的最优视图映射图的示例,其中不同灰度阴影示出了不同的投影缩短和/或重叠值;
图3a是使右动脉斜角(RAO)的含义可视化的组合L型臂/C型臂血管造影系统;
图3b是在中立中间位置的具有所有可能运动角度的组合L型臂/C型臂血管造影系统;
图3c是使左动脉斜角(LAO)的含义可视化的组合L型臂/C型臂血管造影系统;
图4a是使前侧角(CRAN)的含义可视化的组合L型臂/C型臂血管造影系统;
图4b是在中立中间位置的具有所有可能运动角度的组合L型臂/C型臂血管造影系统;
图4c是使后侧角(CRAD)的含义可视化的组合L型臂/C型臂血管造影系统;
图5a是示出了L型臂向左旋转运动90°的组合L型臂/C型臂血管造影系统;
图5b是在中立中间位置的具有所有可能运动角度的组合L型臂/C型臂血管造影系统;
图5c是使L型臂向右旋转运动90°可视化的组合L型臂/C型臂血管造影系统;
图6a是建模的LCA(左冠状动脉)的次优视图;
图6b是根据接近完善的视角的建模的LCA的视图;
图6c是根据图6a和图6b的建模的LCA的位置的最优视图映射图;
图7a-f是C型臂血管造影系统的旋转式滚转采集轨迹的图示,结合有经过OVM的轨迹的图示;
图8a-f是C型臂血管造影系统的倾斜的螺旋桨式旋转运行的图示,结合有经过OVM的轨迹的图示;
图9a-f是C型臂血管造影系统的滚转运动的图示,结合有在螺旋桨式位置固定且L型臂向右定位90°的情况下经过OVM的轨迹的图示;
图10a-f是具有3D-RA系统的C型臂的同步滚转和螺旋桨式运动的双运动采集的图示,结合有经过OVM的轨迹的图示;并且
图11是经过最优视图映射图的C型臂系统的确定的最优轨迹的图示;
图12示出了根据本发明的示例性实施例的简化的计算机系统,该计算机系统可以用作根据本发明的示例性实施例的X射线系统的控制单元。
附图标记列表:
1  血管位置
3  血管位置
5  投影图像平面
7  投影对象
9  病人
11 C型臂
13 L型臂
15 支架
17 L型臂顶点
19 几乎没有投影缩短和重叠的OVM区域
21 白色十字
23 黑色十字
25 C型臂的最优轨迹
27 重叠和投影缩短的灰度标度
29 中央处理单元
31 输入/输出设备
33 主存储器
35 大容量存储器
37 控制单元
39 计算机系统
对于某些治疗和诊断步骤,医师需要针对血管即冠状动脉血管系统的最佳模型。关键问题是如图1a、1b、1c示出的血管投影缩短。如典型地由X射线血管造影系统生成的,三维血管投影在二维投影图像上产生不同的投影。在图1c中的两个位置1和3示出了血管。如果要投影的血管部分与图像平面5平行,则投影达到其最大值。在血管或血管部分的纵轴与图像平面之间的角度大于0的情况下,在图1c的示例中投影缩短了30%的系数。图1a和图1b示出了二维示意图中的效应。取决于对象7和图像平面5之间的不同距离和角度,对象7的投影通过投影缩短效应被缩短。图1a的示例给出了具有0%投影缩短的最佳图像投影。
图3a、图3b、图3c和图4a、图4b、图4c解释了典型C型臂血管造影系统的运动。病人平躺在支架15上。病人身体的纵轴形成用于取向的天然轴线。如果如图3b示出的,血管造影系统处在中立位,则C型臂被安装使其枢轴点在病人纵轴的线上或靠近病人纵轴的线。C型臂11可以围绕该枢轴点以螺旋浆式运动方式旋转。如图3a所示,从如图3b所示的中立位测量的角度典型地被称为RAO(右前斜角),其意指从病人的脚看去的逆时针方向旋转;而如图3c所示,从如图3b所示的中立位测量的角度典型地被称为LAO(左前斜角),其意指与RAO角度旋转相比在相反方向的旋转。此外,示出了L型臂13。L型臂13承载血管造影系统的C型臂11并且处在中立中间位置。
如本领域所公知的,血管造影系统的焦点即X射线发射机和X射线接收机阵列或图像平面定位于C型臂11的相应末端并随C型臂11运动。在血管造影系统的工作周期期间,感兴趣区域(ROI)即病人的心脏典型地处在恒定的位置。
图4a、图4b、图4c示出了C型臂11的另一可能运动。其是滚转运动。该运动用角度CRAN(后侧角)和CAUD(前侧角)来测量。它们典型地从中间位置起达到60°。例如如图3b或图4b所示,假定C型臂11从中立位开始,C型臂11末端的运动处于由病人纵轴和重力矢量所定义的平面。
最后,现代血管造影系统还可以使得L型臂13围绕其顶端17上的枢轴点旋转。该旋转轴被限定在处于检查室天花板上的安装装置和L型臂13之间。L型臂13的典型运动是如图5a所示地向左运动90°,或如图5c所示地向右运动90°。图5b再次示出了整个系统的中立位。
所述的血管造影系统表示具有三个自由度的系统:该三个自由度中的两个由C型臂11的可能运动所定义,而另一个由L型臂13的旋转所定义。然而,这里所讨论的方法不限于具有三个自由度的血管造影系统,其还适用于传统血管造影系统的C型臂的两个自由度的运动。在这种情况下,L型臂13被认为是固定的。利用这样的设置,还可以利用血管造影系统(如基底可移动血管造影系统)的实施例来执行本发明的方法,该可移动血管造影系统没有被安装在房间的天花板上。
图2示出了血管造影系统的典型结果。其是最优视图映射图(OVM)。OVM典型地是彩色的,其中每个像素表示在ROI中0-100%之间的重叠或投影缩短值。OVM中的每个像素由两个角度值表示:一对RAO/LAO,其定义图2的图表中的x方向,以及一对CRAN和CAUD,其定义图表中的y轴。图2是OVM的原始彩色图的黑白转换。根据图2,在OVM中具有浅色区域和深色区域。简单地讲,浅色区域表示具有较低血管投影缩短和/或重叠的区域,但更深色的区域表示具有较高值的投影缩短和/或重叠的区域。
这意指针对给定的角度对,例如RAO=LAO=0°且CRAN=30°,对由血管造影系统产生的二维投影图像建立的三维模型的良好图像重建得到的结果与针对例如RAO=100°且CRAN=30°的RAO/LAO和CRAN/CAUD组合产生的浅色区域(其在图2中由数字19表示)相比会差得多。
这对于医师意味着如果他想得到C型臂血管造影系统的最优位置以便产生为诊断或治疗目的的血管(例如冠状动脉管)的最佳三维重建,他要将C型臂11定位在OVM中的白色区域内的区域。
图6a和图6b示出了血管树的三维模型的示例。图6a是较差的示例,因为医师不能真正想象血管树实际上看起来会怎样。图6b是图像扫描生成的好得多的三维模型。在图6b中三维模型由于非常小的减小后的投影缩短和重叠值,其能够好得多地解释实际冠状动脉树。
在该情况下,假定医师想得到重建血管系统的三维模型的最佳可能图像投影,图6c中的白色十字21可以表示C型臂11的较差位置的RAO/LAO和CRAN/CAUD组合。产生如图6b中所示的血管树的三维模型与图6c中的黑色十字23相关。其表示血管造影系统采集图像的较好的位置以便得到血管树的最佳可解释的三维模型,因为其具有小得多的投影缩短和重叠。
现在,参考图7a-图7e讨论C型臂11的运动和OVM的位置。由于L型臂13向左运动90°,C型臂11的滚转运动不表示其在LAO/RAO方向的旋转,因为这是针对L型臂13处于中立中间位置的情形。这是因为必须根据ROI的病人视图来解释坐标系。由于L型臂13的取向,C型臂11末端的X射线发射机/探测器在从病人看去的CRAN/CAUD方向上运动。因此通过使L型臂13运动90°,LAO/RAO和CRAN/CAUD方向可相互交换。在图7f的OVM上,这产生沿x轴的运动,如由水平白色虚线所指示的。图7c表示具有LAO/RAO=0°和CRAN=CAUD=0°的OVM的中间点。图7a和图7e中的位置将标识图7f中白线的端点。
图8f中示出了类似的结果,即OVM上的轨迹。这里L型臂13再次处在中立中间位置。C型臂11作出滚转运动并且将其保持在固定于CAUD=30°处,同时C型臂11围绕ROI执行螺旋桨式运动。这产生图8f中OVM上的白色虚线。应该注意,示出的位置大致是从RAO=90°到LAO=90°,同时图8f的图表中示出了从RAO=120°到LAO=120°的线。但这没有改变OVM上直线的基本原理。还应该注意,C型臂11的轨迹在明亮区域所指示的CAUD大约30°且LAO大约10°至20°处的几乎没有投影缩短和重叠的区域上运动。然而,通过C型臂11的运动在OVM上也经过深色区域,因此通过使等式(2)最小化不能执行最优轨迹。
图9a-图9f示出了血管造影系统的C型臂在两个自由度上的仿真组合运动。在图9a中,LAO大约90°且CRAN大约10°。C型臂11的运动原则上产生图9f中的白色虚线。应该注意,图9a-图9e并不精确地与轨迹相关,该轨迹产生图9e中的白线,其以RAO=120°开始,在CAUD大约30°时经过y轴并且以LAO=120°且CRAN=CAUD=0°在图9e的右侧结束。然而,重要的是要实现C型臂造影系统能运行这样的轨迹。还应该注意,在图9f中白色虚线指示的血管造影系统的轨迹运动经过OVM的浅色区域和深色区域。
这意指从最小的重叠和投影缩短的视角出发所选择的轨迹不是最优的,因为轨迹没有经过图9f的OVM上的所有明亮区域。
在图10f中这是不同的。这里根据图10a-图10e,白色虚线是C型臂11的双重运动结果。C型臂血管造影系统的角运动以如图10f所示的经过OVM中多于一个白色区域的方式被操纵。
图11最终示出了C型臂11的另一最优轨迹25。所得到的轨迹或多或少地仅在基础OVM的白色区域运动。其基本上从OVM左上侧的象限1移动经过定义为OVM的左下象限的象限4,到定义为右下象限的象限3。轨迹既不接触也不接近由OVM中深色区域象征的高度重叠和/或投影缩短的区域。那些区域中的一个示出为在所示出的OVM的象限2(定义为右上象限)的深色斑。图表左侧的标度27解释了重叠和/或投影缩短的程度。可以看出,轨迹仅位于具有低于10%的重叠和投影缩短的区域。
如上所述,所得到的轨迹基于使等式(2)最小化的本发明的理念。当然该计算必须考虑实际血管造影系统的物理上的限制因素。
图12示出了根据本发明的示例性实施例的计算机系统39,以执行根据本发明的示例性实施例的方法。根据本发明的示例性实施例的软件可使计算机系统执行本发明的示例性实施例的方法步骤,该计算机系统包括用于所有处理目的的中央处理单元29、输入输出设备31、主存储器33、大容量存储器35和控制单元37。该输入/输出设备31典型地包括用于输入命令到该计算机系统的键盘和诸如计算机屏幕或任何其他显示器的可视化设备。主存储器33与CPU一起工作并支持该CPU同时存储用于CPU或数据值的可执行命令。这些组件的设置等价于本领域熟知的冯·诺依曼计算机。大容量存储器37可以存储如从C型臂X射线系统所得到或接收的图像数据的海量数据或任何其他数据、命令和程序代码。并不是经典冯·诺依曼计算机的部分的任选控制单元可以控制C型臂X射线系统的引擎和其他设备或从各个探测器接收数据。由此可以控制和监控X射线血管造影系统的C型臂的运动。计算机系统39还可以与支持X射线血管造影师的必需的其他电子设备具有通信链路。这些通信链路可以连接到输入/输出设备31或控制单元37。连接到其他电子系统的其他连接类型是任选的并且在本领域是熟知的。
应该注意,术语“包括”不排除其他元件或步骤,并且“一”或“一”不排除多个。结合不同实施例描述的元件还可以进行组合。还应该注意,权利要求中的附图标记不应该被解释为限制权利要求的范围。
为了概括本发明的上述各实施例,可以陈述为:中心理念是从已知的最优视图映射图开始,确定C型臂血管造影系统或3D-RA系统的轨迹。经过OVM的C型臂11的轨迹以仅经过具有最小的投影缩短和重叠的区域的方式被操纵。这最终产生用于冠状动脉树或其他血管或其部分的三维重建的最可能的图像投影或二维图像投影以用于治疗。

Claims (13)

1、一种用于为C型臂X射线系统针对血管状结构确定旋转X射线血管造影的最优轨迹(25)的方法,
其中,所述C型臂X射线系统的C型臂具有至少两个自由度,所述自由度由以右冠状动脉斜角或左冠状动脉斜角α表示的所述C型臂(11)的螺旋桨式旋转运动,和以后侧角或前侧角β表示的所述C型臂(11)的滚转运动所定义,所述方法包括下列步骤:
(a)生成感兴趣区域内的体血管中心线的三维表示;
(b)生成最优视图映射图和/或选择由临床医生确定的优选视角;并且
(c)计算最优轨迹(25),其中,最优轨迹(25)至少由所述C型臂(11)在其至少两个自由度内的运动所定义,所述轨迹使得图像投影具有最小的投影缩短和/或重叠,同时使感兴趣区域对X射线的暴露最小化。
2、如权利要求1所述的方法,其中,计算最优轨迹(25)的所述步骤包括使下列等式最小化的步骤:
F ( κ , λ , α , β ) = Σ α = - 120 α = 120 Σ β = - 60 β = 60 κf ( α , β ) + λO ( α , β )
其中:
κ是加权系数;
λ是加权系数;
α是所述左冠状动脉斜角或所述右冠状动脉倾斜角的角度值;
β是所述后侧角或所述前侧角的角度值;
f(α,β)是与血管投影缩短相关的函数;并且
O(α,β)是与血管重叠相关的函数。
3、如权利要求1所述的方法,其中,所述加权系数κ和λ满足下列等式:
κ+λ=1。
4、如权利要求1所述的方法,其中,生成感兴趣区域的中心线的三维表示的所述步骤包括基于两个或更多个采集的投影图像的建模方法。
5、如权利要求1所述的方法,其中,生成感兴趣区域的中心线的三维表示的所述步骤包括基于预采集的不同模态的图像数据集的建模方法。
6、如权利要求5所述的方法,其中,所述不同的模态从包括计算机断层图和磁共振数据的组中选出。
7、如权利要求1所述的方法,其中,生成感兴趣区域的中心线的3D表示的所述步骤包括以基于非最优采集轨迹预采集的三维旋转血管造影图像为基础的建模方法。
8、如权利要求1所述的方法,其中,生成感兴趣区域的中心线的3D表示的所述步骤包括基于表示平均的人体冠状动脉的体模的建模方法。
9、如权利要求8所述的方法,其中,所述体模从包括Dudge模型的一组体模中选出。
10、一种X射线C型臂系统,具有:
C型臂,其具有至少两个自由度,所述自由度由以左冠状动脉斜角或右冠状动脉斜角α表示的所述C型臂(11)的螺旋桨式旋转运动,和以后侧角或前侧角β表示的所述C型臂(11)的滚转运动所定义;
至少一个马达,用于所述C型臂在其至少两个自由度内的运动;以及
控制单元,用于控制所述至少一个马达,从而使所述C型臂的所述运动表示最优轨迹(25);
其中,控制单元适用于基于下列步骤来确定所述最优轨迹(25):
(a)生成感兴趣区域内的体血管中心线的三维表示;
(b)生成最优视图映射图;以及
(c)计算所述C型臂运动的所述最优轨迹(25),所述轨迹使得图像投影具有最小投影缩短和/或重叠,同时使感兴趣区域对X射线的暴露最小化。
11、一种为C型臂(13)X射线系统确定三维旋转X射线冠状动脉血管造影的最优轨迹(25)的计算机系统,其中,所述C型臂X射线系统的C型臂具有至少两个自由度,所述自由度由以左冠状动脉斜角或右冠状动脉斜角α表示的所述C型臂(11)的螺旋桨式旋转运动,和以后侧角或前侧角β表示的所述C型臂(11)的滚转运动所定义,其中,所述计算机系统适用于通过至少执行下列步骤来确定所述C型臂(13)的所述最优轨迹(25):
(a)生成感兴趣区域内的体血管中心线的三维表示;
(b)生成最优视图映射图;以及
(c)计算所述C型臂(11)的所述最优轨迹(25),所述轨迹使得图像投影具有最小投影缩短和/或重叠,同时使感兴趣区域对X射线的暴露最小化。
12、一种在计算机系统中使用的计算机程序产品,其用于为C型臂(11)X射线系统确定三维旋转X射线冠状动脉血管造影的最优轨迹(25),其中,所述C型臂X射线系统的C型臂具有至少两个自由度,所述自由度由以左冠状动脉斜角或右冠状动脉斜角α表示的所述C型臂(11)的螺旋桨式旋转运动,和以后侧角或前侧角β表示的所述C型臂(11)的滚转运动所定义;其中,所述计算机程序产品适用于使所述计算机系统执行下列步骤以确定所述C型臂(11)的所述最优轨迹(25):
(a)生成感兴趣区域内的体血管中心线的三维表示;
(b)生成最优视图映射图;以及
(c)计算所述C型臂(13)的所述最优轨迹(25),所述轨迹使得图像投影具有最小投影缩短和/或重叠,同时使感兴趣区域对X射线的暴露最小化。
13、一种用于存储如权利要求12所述的计算机程序产品的计算机可读介质。
CN2007800242615A 2006-06-28 2007-06-18 基于预确定的最优视图映射图对ra确定最优旋转轨迹 Expired - Fee Related CN101594824B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP06116183 2006-06-28
EP06116183.2 2006-06-28
PCT/IB2007/052324 WO2008001260A2 (en) 2006-06-28 2007-06-18 Optimal rotational trajectory determination for ra based on pre-determined optimal view map

Publications (2)

Publication Number Publication Date
CN101594824A true CN101594824A (zh) 2009-12-02
CN101594824B CN101594824B (zh) 2012-01-11

Family

ID=38846051

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007800242615A Expired - Fee Related CN101594824B (zh) 2006-06-28 2007-06-18 基于预确定的最优视图映射图对ra确定最优旋转轨迹

Country Status (6)

Country Link
US (1) US8463014B2 (zh)
EP (1) EP2046202B1 (zh)
JP (1) JP5415943B2 (zh)
CN (1) CN101594824B (zh)
PL (1) PL2046202T3 (zh)
WO (1) WO2008001260A2 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103126692A (zh) * 2011-11-29 2013-06-05 株式会社东芝 X射线摄影装置以及医用图像处理装置
CN103299346A (zh) * 2010-12-15 2013-09-11 皇家飞利浦电子股份有限公司 用于从3d或4d数据集生成和显示2d投影的系统和方法
CN103356212A (zh) * 2012-04-09 2013-10-23 株式会社东芝 X 射线诊断装置
CN104080405A (zh) * 2011-11-11 2014-10-01 皇家飞利浦有限公司 扩大的角度门控窗口c 型臂图像采集
CN105877767A (zh) * 2016-03-31 2016-08-24 北京思创贯宇科技开发有限公司 一种冠脉图像造影方法和装置
CN106157320A (zh) * 2016-07-29 2016-11-23 上海联影医疗科技有限公司 一种图像血管分割方法及装置
CN107106106A (zh) * 2015-01-14 2017-08-29 皇家飞利浦有限公司 针对具有减小的角范围的旋转c型臂计算机断层摄影的自适应分割
CN107913078A (zh) * 2016-10-11 2018-04-17 西门子保健有限责任公司 用于确定灌注数据组的方法

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008001260A2 (en) * 2006-06-28 2008-01-03 Koninklijke Philips Electronics N. V. Optimal rotational trajectory determination for ra based on pre-determined optimal view map
US20080275467A1 (en) * 2007-05-02 2008-11-06 Siemens Corporate Research, Inc. Intraoperative guidance for endovascular interventions via three-dimensional path planning, x-ray fluoroscopy, and image overlay
JP5491700B2 (ja) 2008-02-14 2014-05-14 株式会社東芝 データ処理装置及びx線装置
EP2317924A1 (en) * 2008-08-11 2011-05-11 Koninklijke Philips Electronics N.V. Determining foreshortening optimal view maps taking into account the shape of the vessel cross-section in cardiovascular x-ray systems
EP2408375B1 (en) 2009-03-20 2017-12-06 Orthoscan Incorporated Moveable imaging apparatus
WO2011038236A2 (en) 2009-09-25 2011-03-31 Ortho Kinematics, Inc. Systems and devices for an integrated imaging system with real-time feedback loops and methods therefor
WO2012082799A1 (en) 2010-12-13 2012-06-21 Orthoscan, Inc. Mobile fluoroscopic imaging system
DE102011006484A1 (de) * 2011-03-31 2012-10-04 Siemens Aktiengesellschaft Angiographiesystem zur angiographischen Untersuchung eines Untersuchungsobjekts und angiographisches Untersuchungsverfahren
EP2570079B1 (en) * 2011-09-13 2017-06-14 Pie Medical Imaging BV Method and apparatus for determining optimal 3D reconstruction of an object
US9105200B2 (en) 2011-10-04 2015-08-11 Quantant Technology, Inc. Semi-automated or fully automated, network and/or web-based, 3D and/or 4D imaging of anatomy for training, rehearsing and/or conducting medical procedures, using multiple standard X-ray and/or other imaging projections, without a need for special hardware and/or systems and/or pre-processing/analysis of a captured image data
US10734116B2 (en) 2011-10-04 2020-08-04 Quantant Technology, Inc. Remote cloud based medical image sharing and rendering semi-automated or fully automated network and/or web-based, 3D and/or 4D imaging of anatomy for training, rehearsing and/or conducting medical procedures, using multiple standard X-ray and/or other imaging projections, without a need for special hardware and/or systems and/or pre-processing/analysis of a captured image data
JP6662943B2 (ja) * 2012-04-09 2020-03-11 キヤノンメディカルシステムズ株式会社 X線診断装置
US9977104B2 (en) 2012-06-04 2018-05-22 Koninklijke Philips N.V. Magnetic resonance imaging along energy-delivering device axis
JP6222801B2 (ja) * 2012-10-26 2017-11-01 東芝メディカルシステムズ株式会社 医用画像処理装置、x線診断装置及び医用画像処理プログラム
JP2014168571A (ja) * 2013-03-04 2014-09-18 Toshiba Corp X線診断装置
US9135272B2 (en) 2013-10-04 2015-09-15 International Business Machines Corporation Finding similar coronary angiograms in a database
WO2015121301A1 (en) * 2014-02-13 2015-08-20 Siemens Aktiengesellschaft Medical imaging optimization
JP5776807B2 (ja) * 2014-03-10 2015-09-09 株式会社島津製作所 X線撮影装置
EP3061970A1 (de) * 2015-02-27 2016-08-31 Arno Hofmann Membranpumpe, insbesondere zur Verwendung im Abgasstrang eines Verbrennungsmotors, und Verbrennungsmotor mit Membranpumpe
JP6625347B2 (ja) 2015-06-01 2019-12-25 キヤノンメディカルシステムズ株式会社 X線血管撮影装置
US9895127B2 (en) * 2015-08-31 2018-02-20 General Electric Company Systems and methods of image acquisition for surgical instrument reconstruction
WO2018001099A1 (zh) 2016-06-30 2018-01-04 上海联影医疗科技有限公司 一种血管提取方法与系统
US10667869B2 (en) 2017-05-17 2020-06-02 General Electric Company Guidance system for needle procedures
EP3459463A1 (en) 2017-09-26 2019-03-27 Koninklijke Philips N.V. Device and method for determining a volume of projection of a dual-axis computed tomography system
CN109949265B (zh) * 2017-12-20 2023-10-03 东芝医疗系统株式会社 医用图像处理装置及医用图像处理方法
US10922897B2 (en) * 2018-01-15 2021-02-16 Canon Medical Systems Corporation Medical information processing apparatus, X-ray diagnostic system, and medical information processing method
EP3705046A1 (en) * 2019-03-07 2020-09-09 Koninklijke Philips N.V. Apparatus for determining a control protocol for controlling a c-arm system
EP3834732A1 (en) * 2019-12-10 2021-06-16 Koninklijke Philips N.V. Field of view matching for mobile 3d imaging
US20210282731A1 (en) * 2020-03-10 2021-09-16 GE Precision Healthcare LLC Systems and methods for registration of angiographic projections with computed tomographic data
US20230190218A1 (en) * 2021-12-20 2023-06-22 GE Precision Healthcare LLC Methods and system for positioning a c-arm

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5381791A (en) * 1992-03-10 1995-01-17 Siemens Medical Systems, Inc. Automatic indentification of anatomical features of interest from data acquired in nuclear medicine studies and automatic positioning of scintillation cameras to carry out such studies at optimal positions
US6047080A (en) * 1996-06-19 2000-04-04 Arch Development Corporation Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images
FR2809849B1 (fr) * 2000-06-05 2002-11-01 Ge Medical Systems Global Tech Cy Llc Procede et dispositif d'acquisition d'images
US7180976B2 (en) * 2002-06-04 2007-02-20 Koninklijke Philips Electronics N.V. Rotational angiography based hybrid 3-D reconstruction of coronary arterial structure
US7020510B2 (en) 2002-07-25 2006-03-28 Koninklijke Philips Electronics, N.V. Optimal view map V.0.01
JP2004329729A (ja) * 2003-05-12 2004-11-25 Toshiba Corp X線診断装置及びx線撮影方法
JP2007502140A (ja) * 2003-08-12 2007-02-08 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 心臓の画像を製造する装置及び方法
JP5296379B2 (ja) 2004-08-13 2013-09-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 冠動脈造影法のための選択的取得
WO2008001260A2 (en) * 2006-06-28 2008-01-03 Koninklijke Philips Electronics N. V. Optimal rotational trajectory determination for ra based on pre-determined optimal view map

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103299346A (zh) * 2010-12-15 2013-09-11 皇家飞利浦电子股份有限公司 用于从3d或4d数据集生成和显示2d投影的系统和方法
CN104080405A (zh) * 2011-11-11 2014-10-01 皇家飞利浦有限公司 扩大的角度门控窗口c 型臂图像采集
CN103126692A (zh) * 2011-11-29 2013-06-05 株式会社东芝 X射线摄影装置以及医用图像处理装置
CN103356212A (zh) * 2012-04-09 2013-10-23 株式会社东芝 X 射线诊断装置
US9259202B2 (en) 2012-04-09 2016-02-16 Kabushiki Kaisha Toshiba X-ray diagnostic apparatus
CN107106106A (zh) * 2015-01-14 2017-08-29 皇家飞利浦有限公司 针对具有减小的角范围的旋转c型臂计算机断层摄影的自适应分割
CN105877767A (zh) * 2016-03-31 2016-08-24 北京思创贯宇科技开发有限公司 一种冠脉图像造影方法和装置
CN106157320A (zh) * 2016-07-29 2016-11-23 上海联影医疗科技有限公司 一种图像血管分割方法及装置
CN106157320B (zh) * 2016-07-29 2019-02-01 上海联影医疗科技有限公司 一种图像血管分割方法及装置
CN107913078A (zh) * 2016-10-11 2018-04-17 西门子保健有限责任公司 用于确定灌注数据组的方法

Also Published As

Publication number Publication date
WO2008001260A3 (en) 2008-07-10
US8463014B2 (en) 2013-06-11
WO2008001260A2 (en) 2008-01-03
US20100014740A1 (en) 2010-01-21
JP5415943B2 (ja) 2014-02-12
CN101594824B (zh) 2012-01-11
EP2046202B1 (en) 2014-06-11
PL2046202T3 (pl) 2015-04-30
EP2046202A2 (en) 2009-04-15
JP2009542283A (ja) 2009-12-03

Similar Documents

Publication Publication Date Title
CN101594824B (zh) 基于预确定的最优视图映射图对ra确定最优旋转轨迹
US20180204330A1 (en) Tomography apparatus and method for reconstructing tomography image thereof
CN102598059B (zh) 利用偏移探测器几何结构生成计算断层摄影图像的方法和设备
JP5491700B2 (ja) データ処理装置及びx線装置
US7180976B2 (en) Rotational angiography based hybrid 3-D reconstruction of coronary arterial structure
JP6175073B2 (ja) 最適なデバイスナビゲーションのための血管系ビューのリアルタイム表示
US8625862B2 (en) Image display method and image display apparatus
US9576391B2 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
CN101313334B (zh) 高对比度对象的运动补偿ct重建
US10229516B2 (en) Method and apparatus to improve a 3D + time reconstruction
KR101725891B1 (ko) 단층 영상 장치 및 그에 따른 단층 영상 복원 방법
CN101874740B (zh) 图像显示装置以及x射线诊断装置
EP3107457B1 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
CN107427274A (zh) 断层扫描设备及其用于重构断层扫描图像的方法
CN101902967B (zh) 心脏ct中非自主呼吸运动的校正
CN102232842B (zh) 用于心脏计算机断层造影拍摄的改进的多片段图像重建
Wink et al. Intra-procedural coronary intervention planning using hybrid 3-dimensional reconstruction techniques1
CN106651985A (zh) Ct图像的重建方法和装置
CN101868809A (zh) 具有计算机辅助检测和/或诊断的图像处理
US20230230243A1 (en) Methods, devices, and systems for dynamic fluoroscopy of c-shaped arm devices
Wu et al. Validation study of 3-D reconstruction of coronary arteries

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120111

Termination date: 20180618