CN102292744B - 心脏图像处理和分析 - Google Patents
心脏图像处理和分析 Download PDFInfo
- Publication number
- CN102292744B CN102292744B CN201080005093.7A CN201080005093A CN102292744B CN 102292744 B CN102292744 B CN 102292744B CN 201080005093 A CN201080005093 A CN 201080005093A CN 102292744 B CN102292744 B CN 102292744B
- Authority
- CN
- China
- Prior art keywords
- curved surface
- image
- motion
- cardiac
- heart
- 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
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims description 23
- 230000008569 process Effects 0.000 title description 8
- 238000004458 analytical method Methods 0.000 title description 2
- 230000000007 visual effect Effects 0.000 claims abstract description 47
- 238000012800 visualization Methods 0.000 claims abstract description 28
- 238000005070 sampling Methods 0.000 claims abstract description 27
- 210000004165 myocardium Anatomy 0.000 claims abstract description 24
- 230000000295 complement effect Effects 0.000 claims abstract description 9
- 210000005240 left ventricle Anatomy 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims description 7
- 230000002123 temporal effect Effects 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 239000000203 mixture Substances 0.000 description 5
- 210000005003 heart tissue Anatomy 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 238000009472 formulation Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 230000002969 morbid Effects 0.000 description 3
- 238000002600 positron emission tomography Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000009206 nuclear medicine Methods 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 238000007794 visualization technique Methods 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 208000006083 Hypokinesia Diseases 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000013184 cardiac magnetic resonance imaging Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000003483 hypokinetic effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Clinical applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Cardiology (AREA)
- Optics & Photonics (AREA)
- Dentistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
一种用于可视化由心脏图像表示的心肌的系统,包括再采样装置和可视化装置。所述再采样装置在多个曲面上的采样点处再采样强度水平,每个曲面包围心腔的至少一部分以及所述多个曲面中的零个或多个,并被所述多个曲面中的其余曲面包围,所述多个曲面一起覆盖心脏图像中的中空区域,所述中空区域包括一组至少一个心腔的外腔壁。所述可视化装置被布置成用于使用从所述再采样装置获得的再采样强度水平使所述多个曲面中的至少一个的至少一部分可视化。所述组至少一个心腔可以仅是左心室。其也可以是整个心脏。
Description
技术领域
本发明涉及对心脏图像的可视化(visualizing)。本发明还涉及对心脏图像序列的运动补偿。本发明还涉及应变成像。
背景技术
心脏组织的运动可以通过使用一个或多个不同的2D和3D成像模式(特别是示踪或未示踪的心脏MRI,和/或不同模式的超声成像)观察心脏组织的医学图像来评估。根据这样的动态成像,能够提取组织速度、应变和/或应变率。可以使用参数图像来使这样的应变和应变率可视化,其中,与位置相关联的强度表示参数值,诸如应变或应变率。这样的参数成像技术能够提供有价值的诊断信息,其中,所述诊断信息涉及心肌特定部分的运动功能减退、运动不能(akinesy)、运动障碍(diskinesy)或者心肌异步。然而,这种参数成像的有效性可能取决于图像对比度以及取决于估计的运动场(field)的精确性。
Ledesma-Carbayo、M.J.Santos、A.Kybic、J.Mahia-Casado、P.Garcia-Fernandez、M.A.Malpica、N.Perez-David、E.Desco,M.于2004年9月发表在ComputersinCardiology第313-316页上的“Myocardialstrainanalysisofechocardiographicsequencesusingnonrigidregistration”(下文称为Ledesma-Carbayo等人)一文描述了评估所提取的运动场的精确性的方法。使用所提取运动场,以这样的方式生成运动补偿的图像序列,即,使得该序列中的所有图像类似于所选择的参考帧(frame)。通常,在心脏舒张末期的帧中选择参考帧。
用于可视化心脏图像的技术是已知的,例如使用直接体绘制或多平面重构(reformat)。然而,这些可视化未提供详细检查心肌的有效方法。
发明内容
提供一种使心脏图像表示的心肌可视化的经改进的方法将较为有利。为了更好地解决这一问题,在本发明的第一方面中提供了一种系统,其包括:
-再采样装置(means),其用于在多个曲面上的采样点处再采样强度水平,每个曲面包围所述多个曲面中的零个或多个并且被所述多个曲面中的其余曲面包围,曲面包围心腔的至少一部分,所述多个曲面一起覆盖心脏图像中的中空区域,所述中空区域包括一组至少一个心腔的外腔壁;以及
-可视化装置,其用于使用从再采样装置获得的再采样的强度水平使所述多个曲面中的至少一个的至少一部分可视化。
由于曲面包围心腔的一部分,并且曲面彼此包围并一起覆盖包括腔壁的区域,这些表面或多或少地遵循如心脏图像中所示的包含心肌的心脏壁的形状。这些表面中的至少一些与心肌交叉,并且局部地近似平行于心肌。对这样的表面的可视化允许检查心肌层。所述表面例如等距地分布于中空区域上。所述表面可以形成适应于心肌形状的同心外罩(shell)。该中空区域可以是三维区域或体积。中空区域可以在其中具有孔。例如,中空区域包括所述组至少一个心腔的外腔壁,其中,至少一个心腔的至少一部分从所述区域去除以形成孔。例如,可能从所述区域去除多于一半的心腔以形成所述区域的孔。还可以去除整个心腔;在那种情况下,所述区域的内边界与腔壁的内边界吻合。曲面可以分布在中空区域中。例如,曲面可以采取围绕孔的同心外罩的形式。
例如,所述组至少一个心腔可以仅是左心室。在这种情况下,所述外壁是左心室的壁。或者,所述组至少一个心腔可以是一组四个心腔室,外壁是心脏的外壁。在后者的情况中,连同曲面一起覆盖从心外膜延伸至心内膜的区域。
可以提供适配装置,用于使曲面的形状适应于腔壁的形状,使得所述曲面的形状相似于腔壁。因此,较大部分的表面可以与心肌吻合。当使这样的表面可视化时(例如以展平的方式),可以使相对大部分的心肌、或者甚至整层心肌可视化。
可以提供展平装置,用于展开和展平多个曲面中的曲面,以获得展平的表面,其中,曲面上的采样点被映射到展平的表面上的点。展平装置将曲面变形(或弯曲)成由展平的表面所指示的基本上为平面的表面。可以通过使至少部分展平的表面可视化来实现对多个曲面中的至少一个的至少一部分的可视化。这种展平的表面可视化提供了对心脏周围的心脏图像的快速概览。
可以提供收集装置,用于组合多个展平的表面,以形成包括一堆(stack)展平的表面的体积图像。可以通过堆叠连续展平的表面而将展平的表面组合成体积图像。所述堆中的连续展平的表面对应于连续的较大曲面。
可以提供体积可视化装置,用于使体积图像可视化。原则上,任何体积可视化模式都可以用于使体积图像可视化,包括例如切片视图和最大或最小强度投影(MIP)
体积可视化装置可以包括正交查看器(orthoviewer)。该正交查看器是本领域已知的观察模式,其可以允许用户浏览体积图像的切片,一次使体积图像的三个正交定向的切片可视化。
可视化装置可以被布置成用于使多个曲面中的曲面可视化,该曲面被可视化为与多个曲面中的其余曲面隔离。该曲面可以例如以其原始曲面、或者以如上所述的展平之后而被可视化。这种曲面的隔离视图可以允许精密地检查心肌的薄外罩。
根据本发明一方面,所述系统可以被布置成用于操作覆盖心跳的至少一部分的心脏图像的时间序列,并且所述系统还可以包括:
-运动补偿装置,其用于通过变换心脏图像的时间序列中的心脏图像来匹配参考心脏图像从而补偿心脏运动,以获得经运动补偿的图像序列;
-再采样装置,其被布置成用于使用同样的多个曲面执行对多幅经运动补偿的图像的再采样。
可以例如在统计分析中使用由此生成的再采样点的数据集,以识别任何残余运动。可以使用上述可视化模式而动态地可视化所述数据集,继而显示时间序列的顺序(subsequent)再采样图像。这种可视化允许识别残余运动。所述可视化因而使得用户能够评估运动补偿是否适当。运动补偿装置本身在本领域中是已知的。
运动补偿装置可以被布置成用于采用刚性或仿射变换,从而将心脏图像变换成经运动补偿的图像。这种刚性或仿射变换本身在本领域中是已知的。通过将变换限制为刚性或仿射变换,可以补偿大部分“健康”类型的心脏运动。病态或无规律的心脏运动不能由这种变换补偿。例如,异步性通常不能通过使用刚性或仿射变换而补偿。因此,这些病态或无规律的心脏运动保留在经运动补偿的图像中作为残余运动。与上述再采样和可视化模式组合,用户能够相对有效地识别残余运动,并且因而能够相对有效地识别病态或无规律的心脏运动或异步性。
可以提供顺序可视化装置,用于按照心脏图像序列的时序顺序动态可视化与所述心脏图像序列对应的至少一些再采样的点。所述可视化允许识别残余运动。
可以提供指示装置,用于使用户能够指示包括残余运动的区域;而运动补偿装置可以被布置成用于进一步运动补偿用户所指示的区域。对如上所述的心脏图像序列的可视化允许容易地识别残余运动。通过使得用户能够向系统指明这种残余运动,该系统能够施加专用于该区域的其他运动补偿。
所述指示装置可以被布置成用于进一步使用户能够指示残余运动的方向,所述其他运动补偿装置被布置成用于根据用户指示的残余运动的方向补偿用户所指示的区域。通过使用该方向,进一步改善了运动补偿。
所述系统可以被整合到医学工作站或成像采集设备中。
一种可视化由心脏图像表示的心肌的方法,所述心脏图像将强度水平与体积中的三维点相关联,所述方法可以包括:
-在多个曲面上的采样点处再采样强度水平,每个曲面包围多个曲面中的零个或多个并且被所述多个曲面中的其余曲面包围,所述多个曲面一起覆盖心脏图像中的中空区域,所述中空区域包括一组至少一个心腔的外腔壁;以及
-使用从再采样装置获得的再采样的强度水平使所述多个曲面中的至少一个的至少一部分可视化。
可以提供一种计算机程序产品,其包括用于令处理器系统执行所述方法的指令。
本领域技术人员应当认识到,本发明的两个或多个上述实施例、实施方式和/或方面可以用任何有用的方式组合。
本领域技术人员基于本说明可以实现图像采集设备、工作站、系统和/或与系统的所述修改和变型对应的计算机程序产品的修改和变型。
本领域技术人员应当认识到,所述方法可以被应用于多维图像数据,例如应用到2维(2-D)、3维(3-D)或4维(4-D)图像,所述图像可以由各种采集器械采集,诸如,但不限于,标准X射线成像,计算机断层摄影(CT)、磁共振成像(MRI)、超声(US)、正电子发射断层摄影(PET)、单光子发射计算机断层摄影(SPECT)和核医学(NM)。
本发明由独立权利要求限定。在从属权利要求中限定了有利实施例。
附图说明
将参考附图进一步阐明和描述本发明的这些和其他方面,在附图中
图1图示了用于使心脏图像表示的心肌可视化的系统;
图2图示了用于使心脏图像表示的心肌可视化的方法;
图3图示了确认和改进运动场的过程(process);
图4图示了曲面;
图5图示了正交查看器的三个象限;以及
图6图示了适于心脏形状的经调整的曲面。
具体实施方式
例如,可以从成像模式(多普勒超声(US)或专用的MRI模式)中或者通过使用例如Ledesma-Carbayo等人所述的刚性或非刚性配准技术直接提取心脏组织运动场。然而,需要解决若干困难,以便校正和评估估计的运动场的有效性。
所采集的图像对比度在需要评估运动的区域可能较低。当使用3D时间序列(4D可视化)时,检查经运动补偿的图像序列可能是单调而冗长的任务。
图1图示了用于可视化由心脏图像表示的心肌的系统,所述心脏图像将强度水平与体积中的三维点相关联。这种心脏图像被称为体积心脏图像。正如本领域中已知的,所述系统可以包括具有中央处理单元、存储器和存储装置的计算机系统。本文所述的一些或所有装置和模块可以实施为软件。
可以提供再采样装置10,用于在多个曲面上的采样点处再采样强度水平。每个曲面包围心腔的至少一部分以及多个曲面中的零个或多个,并且被所述多个曲面的其余曲面所包围。所述多个曲面一起覆盖心脏图像中的中空区域。所述中空区域可以包括一组至少一个心腔的外腔壁。例如,所述区域的形状为绕一个或多个心腔的外壳的形状,但所述区域可以部分地与所述一个或多个心腔重叠。
可以提供可视化装置11,用于使用从再采样装置获得的再采样的强度水平使所述多个曲面中的至少一个的至少一部分可视化。可视化装置11可以包括展平装置12,其用于展开和展平所述多个曲面中的曲面,以获得展平的表面,其中,将曲面上的采样点映射到展平的表面上的点。所述系统可以包括收集装置13,其用于组合多个展平的表面,以形成包括一堆展平的表面的体积图像。这将在下文中阐述。可视化装置11可以包括体积可视化装置14,其用于例如使用正交查看器可视化体积图像。
可视化装置14可以被布置成用于可视化多个曲面中的曲面,所述曲面被可视化为与所述多个曲面中的其余曲面隔离。例如,这发生在正交查看器中。
所述系统可以被布置成用于处理覆盖心跳的至少一部分的心脏图像的时间序列。这种图像序列可以经由例如网络连接的输入接收,或者可以存储在局部存储器中。所述系统还可以包括运动补偿装置16,其用于通过变换心脏图像的时间序列中的心脏图像来匹配参考心脏图像,从而获得经运动补偿的图像序列,以补偿心脏运动。再采样装置10可以被布置成用于使用相同的多个曲面在多幅经运动补偿的图像中执行再采样。因此,在相同的点处采样顺序的经运动补偿的图像。这样,例如,可以创建关于所述曲面中的一个的强度的影像。利用收集装置13,还可能收集顺序经运动补偿的图像的再采样数据,以形成体积图像序列。这一序列可以使用诸如正交查看器的体积观察器进行动态观察。
运动补偿装置16可以被布置成用于采用刚性或仿射变换,从而将心脏图像变换成经运动补偿的图像。可视化利用所述刚性或仿射变换形成的经运动补偿的图像允许观察者检查病态或异步性。
可视化装置11可以包括顺序可视化装置17,其用于按照心脏图像序列的时序顺序动态地可视化与所述心脏图像序列对应的至少一些再采样的点。
可以提供指示装置18,其用于使得用户能够指示包括残余运动的区域。运动补偿装置16可以被布置成用于对用户所指示的区域进行进一步的运动补偿。
指示装置18可以被布置成用于进一步使得用户能够指示残余运动的方向,运动补偿装置16被布置成用于根据用户所指示的残余运动的方向而补偿用户所指示的区域。
例如,所述系统可以包括在医学工作站中。所述系统还可以包括在医学图像采集设备中。
图2图示了可视化由心脏图像表示的心肌的方法,所述心脏图像将强度水平与体积中的三维点相关联。所述方法在步骤51中包括在多个曲面上的采样点处再采样强度水平,每个曲面包围心腔的至少一部分以及多个曲面中的零个或多个,并且被所述多个曲面中的其余曲面包围,所述多个曲面一起覆盖心脏图像中的中空区域,所述中空区域包括一组至少一个心腔的外腔壁。所述方法还包括在步骤52中使用从再采样装置获得的再采样强度水平使所述多个曲面中的至少一个的至少一部分可视化。
可以使用于精确运动估计和参数成像的可视化适于心脏形状。尤其是在以下情况下,即,如果可以针对参考图像获得心脏的分割。
可以应用可视化技术,以使得用户能够评估估计的组织运动场的质量或可靠性。此外,可以自动高亮显示其估计值可疑的区域。可以采用交互程序在感兴趣的小型区域或体积(分别是ROI或VOI)中校正组织运动场,并局部地改善ROI/VOI中的运动场。使用对心脏的分割,可以使得可视化几何结构适于心脏形态。然而,这并非是限制。例如,可以使用心脏的模型调整可视化几何形状。
图3图示了核实和改善运动场的处理。在步骤1中,可以执行和保存对心脏图像序列的采集。例如,采取3D超声采集或者心脏CT图像采集。
在步骤2a中,可以估计运动场。在步骤2b中,可以使用估计的非刚性运动生成经运动补偿的图像序列。可以使用能够评估和校正特定ROI/VOI中的运动场的任何非刚性图像配准技术。在S.Makram-Ebeid和O.Somphone于2007年10月14日至17日在里约热内卢(RiodeJaneiro)的ProceedingsoftheICCV2007Conference上发表的“Non-RigidImageRegistrationusingaHierarchicalPartitionofUnityFiniteElementMethod”(下文称为:Makram-Ebeid等)一文中,公开了一种合适的配准算法。可以用下列方式执行补偿,即,使得序列中的所有图像相似于所选择的图像(通常是心脏舒张末期图像),换言之,已经移动的心脏的部分“弯曲”至参考图像中的相应位置。
在步骤3中,可以使用经运动补偿的图像的动态可视化,以使得用户能够交互地检查解剖相关的部分,从而查看心脏组织是否显现为不移动。附加地或备选地,可以通过计算局部空间和/或时间图像变量而检测外观上的运动。可以在与每个图像元素的特定相邻的b体素内评估空间变量,其中b例如可以是3至5。同样地,可以在当前、前一或后一时间帧或多个时间帧中计算时间变量。残余外观运动的区域可以是图像强度的时间变量大于其空间变量的那些区域。那些区域可以在动态可视化中自动高亮显示。可以有利地使用在参考帧中对心脏的分割,以便限定适合于心脏解剖结构的坐标系。使用这一坐标系可以重构和可视化图像数据。
根据用户的请求,可以从估计的运动场提取运动、应变和/或应变率参数,并且可以使用参数成像程序加以可视化。还可以定量估计和显示应变和/或运动同步性或异步性。用于定量估计应变和/或运动的技术在本领域中是已知的。
如步骤3a中所示,如果图像质量足够好并且在经补偿的图像序列中基本上没有外观运动,那么可以维持所得到的参数成像中的置信度。在已经维持置信度的情况下,用户能够决定可视化用于诊断目的的参数图像(步骤4)。否则,过程可以继续进行至步骤5a。
可能由对比度不高的图像区域引起步骤3中高亮显示的残余运动。在包括用户交互作用的处理中,可以移除这种残余运动。在步骤5a中,用户可以指示需要校正运动评估的区域,例如通过限定覆盖该区域的空间ROI/VOI。运动矢量需要更新的方向还可以由用户指示(横向和/或纵向和/或环向)。ROI//VOI允许限定模糊窗口函数,在所述模糊窗口函数内,可以再计算运动场校正(在步骤5b)。使用Makram-Ebeid等人所公开的技术可以再计算所示区域中的运动场校正。在步骤6中,随后将局部校正应用于对运动场的先前的估计,以获得经改进的整个运动场。
所述过程可以迭代地进行步骤2b、3、3a、5a、5b和6。如果确定图像质量不合适或者如果上述校正程序未提供可接受的结果,那么可以执行新的采集(步骤1)。
图4图示了坐标系,其允许遵循心脏或心腔室(诸如左心室)的解剖形状。在附图中,相同的部分通常已经标注为相同的附图标记。图4A示出了曲面401的横截面视图。然而,该曲面仅是范例,而非限制。曲面401是闭合的表面。然而,这并非是限制。例如,可以使用曲面401的一部分。此外,附图图示了轴404,其从点405延伸至点406、近似平行于连接心尖和心底的线。参数值r(在图4中的若干位置处图示的)表示与轴404的距离。曲面401在其各处与轴404具有相同的距离r。此外,参数θ表示沿箭头408指示的方向围绕轴404的曲面401上的位置。例如,这可以用毫米表达为在逆时针方向上沿着该曲线404与在该曲线401上的预定参考点的距离。还可以对该参数θ使用其他单位。图4B图示了曲面401的另一横截面视图,在图4A中图示了图4B的横截面视图的取向。参数表示如箭头409所指示的围绕轴的角位置(例如以度为单位,从预定参考角位置开始)。可以选择轴404的长度、位置和取向以及半径r,使得曲面401或多或少地对应于曲面401,并且尤其使得由此选定给定的轴404,则能够选择最小半径r,其具有与之相关联的第一曲面401,并且可以选择最大半径r2,其具有与之相关联的第二曲面407,使得心室的心肌(或者在需要时是整个心脏的外壁)完全包含在两个曲面401和407之间。
可以使用以上引入的参数θ、和r,以规律(regular)的间隔,例如在规律的栅格点或者矩形的栅格点,再采样所有或部分心脏图像。这是“展开”表面的方式。可以对时间图像序列中的所有图像使用相同的坐标系。特别是如果时间图像序列包括经运动补偿的图像,可以实现第一曲面401和第二曲面407之间的相对窄的带,同时捕获每幅图像中的心肌。再采样可以形成新的矩形体积图像,其中,参数θ、和r形成正规的矩形三维坐标空间。该矩形体积可以例如在正交查看器中进行观察。可以在该正交查看器中显示经运动补偿、再采样图像的时间序列。用户可以在所述正交查看器中指示感兴趣区域(ROI),应当在所述感兴趣区域改进运动估计。使用所述坐标系,可以在原始(经运动补偿的)心脏图像中识别这一感兴趣区域,这使得能够在该感兴趣区域中改进运动补偿。
可以针对心脏的形状,具体而言,针对心壁或心肌的形状调整曲面401。这一调整可以基于对心脏的分割。获得这种适于心脏形状的表面的技术是适应性形状模型。适应性形状模型在本领域中是已知的。
图5示出了正交查看器。显示了三个象限。象限501示出了沿着r和θ轴的再采样体积的切片。可以改变以观察不同的切片。象限502示出了沿着r和轴的再采样的体积的切片。可以改变θ以观察不同的切片。象限503示出了沿着θ和轴的再采样的体积的切片。可以改变r以观察不同的切片。用户可以选择他想要观察的切片。例如,可以使用上右侧象限用于不同的可视化模式。正交查看器是可视化展平的表面的一种方式。
如上所示,表面可以适于心壁的形状。在这种情况下,在图5的正交查看器的象限中显示的图像的区域遵循表面所覆盖的区域;具体而言,所显示的r的值的范围可以根据参数θ和而改变。所述范围中的最小值对应于最接近心脏中心(或者心腔的中心,因为这种情况是可能的)的表面,而所述范围中的最大值可以对应于离该中心最远的表面。
图6示出了已经针对心脏形状进行调整的曲面601的交叉。dθ指示沿着θ轴的一阶(step)的意思。在附图中还已经指出了半径r的一些实例。
可以在例如移动工作站中的车载(on-cart)和/或非车载(off-cart)超声回波描记系统中使用上述技术。相同的方法可以用于可视化或采集心脏成像工作站,例如MRI心脏成像工作站。
应当认识到,本发明还扩展到计算机程序,具体扩展到在载体上或载体内的计算机程序,其适于将本发明付诸实践。所述程序可以是源代码、目标代码媒介源和目标代码的形式,诸如部分编译的形式,或者是适于在根据本发明的实施中使用的任意其他形式。还应当认识到,这样的程序可以具有许多不同的架构设计。例如,实施根据本发明的方法或系统功能的代码可以细分成一个或多个子进程。在这些子进程之间分配功能的许多不同方式对于本领域技术人员而言显而易见。所述子进程可以用可执行文档的形式存储在一起以形成自包含(self-contained)的程序。这样的可执行文档可以包括计算机可执行指令,例如处理器指令和/或解释指令(例如Java解释指令)。备选地,子进程中的一个、多个或全部存储在至少一个外部库文档中,并例如在运行时与主程序静态地或动态地链接。主程序包含对子进程中的至少一个的至少一次调用。同样地,子进程可以包括对彼此的功能调用。一个涉及计算机程序产品的实施例,包括与前述方法中的至少一个的处理步骤中的每个对应的计算机可执行指令。这些指令可以被细分为子进程和/或存储在可以静态或动态链接的一个或多个文档中。涉及计算机程序产品的另一实施例包括与前述系统和/或产品中的至少一个的模块中的每个对应的计算机可执行指令。这些指令可以被细分成子进程和/或存储在可以静态和/或动态链接的一个或多个文档中。
计算机程序的载体可以是能够承载程序的任何实体或装置。例如。所述载体可以包括诸如ROM存储介质,例如CDROM或半导体ROM;或磁记录介质,例如软盘或硬盘。此外,所述载体可以是可传输的载体,诸如电或光信号,其可以经由电缆或光缆或者通过射频或其他手段传递。当所述程序被嵌入到这种信号中时,可以通过线缆或其他装置或模块指定载体。备选地,所述载体可以是其中嵌入有所述程序集成电路,所述集成电路适于执行、或用于执行相关的方法。
应当认识到,上述实施例用于说明而非是限制本发明,并且本领域技术人员能够在不背离本发明权利要求范围的情况下设计许多备选实施例。在权利要求中,置于圆括号之间的任何参考标记不应当解释为对权力要求构成限制。“包括”一词以及其结合的使用并不排除权利要求中所述的那些元件或步骤之外元件或步骤。不定冠词“一”或“一个”并不排除多个此类元件的存在。本发明可以通过包括若干分立元件的硬件方式实施,以及通过适当编程的计算机的方式实施。在装置权利要求中列举了若干模块,这些模块中的一些可以通过硬件中的一个或相同内容实现。在相互不同的从属权利要求中所应用的特定措施并不指示不能有利的使用这些措施的组合。
Claims (14)
1.一种用于可视化由心脏图像表示的心肌的系统,所述心脏图像将强度水平与体积中的三维点相关联,所述系统包括:
-再采样装置(10),其用于在多个曲面上的采样点处再采样所述强度水平,每个曲面包围所述多个曲面中的零个或多个并且被所述多个曲面中的其余曲面包围,所述曲面包围心腔的至少一部分,所述多个曲面一起覆盖所述心脏图像中的中空区域,所述中空区域包括至少一个心腔的外腔壁;以及
-可视化装置(11),其用于使用从所述再采样装置获得的再采样的强度水平,使所述多个曲面中的至少一个的至少一部分可视化。
2.根据权利要求1所述的系统,其中,所述至少一个心腔是左心室,并且所述外腔壁是所述左心室的外壁。
3.根据权利要求1所述的系统,还包括适配装置,其用于针对所述外腔壁的形状调整所述曲面的形状,令所述曲面具有与所述外腔壁相似的形状。
4.根据权利要求1所述的系统,其中,所述可视化装置(11)包括展平装置(12),其用于展开和展平所述多个曲面中的曲面,以获得展平的表面,并且其中,将所述曲面上的所述采样点映射到所述展平的表面上的点。
5.根据权利要求4所述的系统,还包括收集装置(13),其用于组合多个展平的表面,以形成包括一堆所述展平的表面的体积图像。
6.根据权利要求5所述的系统,其中,所述可视化装置(11)包括体积可视化装置(14),其用于可视化所述体积图像。
7.根据权利要求6所述的系统,其中,所述体积可视化装置(14)包括正交查看器。
8.根据权利要求1所述的系统,其被布置成用于操作覆盖心跳的至少一部分的心脏图像的时间序列,所述系统还包括:
-运动补偿装置(16),其用于通过变换所述心脏图像的时间序列中的心脏图像来匹配参考心脏图像,以补偿心脏运动,从而获得经运动补偿的图像的序列;
-所述再采样装置(10),其被布置成用于使用相同的多个曲面执行对多幅所述经运动补偿的图像的所述再采样。
9.根据权利要求8所述的系统,其中,所述运动补偿装置(16)被布置成用于采用刚性或仿射变换将心脏图像变换成经运动补偿的图像。
10.根据权利要求8所述的系统,其中,所述可视化装置(11)包括顺序可视化装置(17),其用于按照心脏图像的序列的时序顺序动态可视化与所述心脏图像的序列对应的至少一些再采样的点。
11.根据权利要求10所述的系统,还包括:
-指示装置(18),其用于使得用户能够指示包括残余运动的区域;并且
-所述运动补偿装置(16)被布置成用于进一步运动补偿由所述用户指示的区域。
12.根据权利要求11所述的系统,其中,所述指示装置(18)被布置成用于进一步使得所述用户能够指示所述残余运动的方向,所述运动补偿装置(16)被布置成用于根据所述用户指示的所述残余运动的方向补偿由所述用户指示的所述区域。
13.一种包括根据权利要求1所述的系统的医学工作站。
14.一种可视化由心脏图像表示的心肌的方法,所述心脏图像将强度水平与体积中的三维点相关联,所述方法包括:
-在多个曲面上的采样点处再采样(51)所述强度水平,每个曲面包围所述多个曲面中的零个或多个并且被所述多个曲面中的其余曲面包围,所述曲面包围心腔的至少一部分,所述多个曲面一起覆盖所述心脏图像中的中空区域,所述中空区域包括至少一个心腔的外腔壁;以及
-使用从所述再采样获得的再采样强度水平,使所述多个曲面中的至少一个的至少一部分可视化(52)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP09305062.3 | 2009-01-23 | ||
EP09305062 | 2009-01-23 | ||
PCT/IB2010/050210 WO2010084446A1 (en) | 2009-01-23 | 2010-01-18 | Cardiac image processing and analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102292744A CN102292744A (zh) | 2011-12-21 |
CN102292744B true CN102292744B (zh) | 2015-11-25 |
Family
ID=42115677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080005093.7A Expired - Fee Related CN102292744B (zh) | 2009-01-23 | 2010-01-18 | 心脏图像处理和分析 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8848989B2 (zh) |
EP (1) | EP2389662A1 (zh) |
JP (1) | JP5996870B2 (zh) |
CN (1) | CN102292744B (zh) |
WO (1) | WO2010084446A1 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012132307A1 (en) | 2011-03-25 | 2012-10-04 | Semiconductor Energy Laboratory Co., Ltd. | Lithium-ion secondary battery |
EP2903510B1 (en) | 2012-10-02 | 2019-12-04 | NewStim, Inc. | Ventricular pacing in cardiac-related applications |
JP6401459B2 (ja) * | 2014-02-14 | 2018-10-10 | キヤノン株式会社 | 画像処理装置、画像処理方法 |
WO2015124388A1 (en) * | 2014-02-19 | 2015-08-27 | Koninklijke Philips N.V. | Motion adaptive visualization in medical 4d imaging |
KR101738364B1 (ko) | 2015-12-02 | 2017-05-22 | 숭실대학교산학협력단 | 4D(3D+t) CTA 영상 정보를 이용한 4D(3D+t) 심혈관 모델 보간 방법, 이를 수행하기 위한 기록 매체 및 장치 |
US10321878B2 (en) | 2016-12-22 | 2019-06-18 | Biosense Webster (Israel) Ltd. | Pulmonary vein display in two dimensions |
CN110383345B (zh) | 2017-01-12 | 2023-08-15 | 纳维斯国际有限公司 | 用于内腔导航的扁平化视图 |
US11282191B2 (en) | 2017-01-12 | 2022-03-22 | Navix International Limited | Flattened view for intra-lumenal navigation |
JP7462672B2 (ja) * | 2019-04-02 | 2024-04-05 | コーニンクレッカ フィリップス エヌ ヴェ | 超音波撮像におけるセグメンテーション及びビューガイダンス並びに関連するデバイス、システム及び方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5797843A (en) * | 1992-11-03 | 1998-08-25 | Eastman Kodak Comapny | Enhancement of organ wall motion discrimination via use of superimposed organ images |
US6757414B1 (en) * | 2000-08-31 | 2004-06-29 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for segmentation of a left ventricular epicardium |
CN1969298A (zh) * | 2004-04-23 | 2007-05-23 | 美国西门子医疗解决公司 | 用于多对象体数据可视化的方法和系统 |
CN101110124A (zh) * | 2006-02-21 | 2008-01-23 | 美国西门子医疗解决公司 | 利用虚拟切割的背景中体积可视化的系统及方法 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3045642B2 (ja) | 1994-01-25 | 2000-05-29 | アロカ株式会社 | 超音波診断装置 |
US5711305A (en) | 1995-02-17 | 1998-01-27 | Ep Technologies, Inc. | Systems and methods for acquiring endocardially or epicardially paced electrocardiograms |
CA2368237A1 (en) * | 1999-04-22 | 2000-11-02 | The Johns Hopkins University | Cardiac motion tracking using cine harmonic phase (harp) magnetic resonance imaging |
AU1459501A (en) | 1999-11-05 | 2001-06-06 | Zynergy Cardiovascular (Zvc) Inc. | Cardiac mapping systems |
JP4589490B2 (ja) | 2000-07-12 | 2010-12-01 | 株式会社東芝 | 超音波診断装置及び記録媒体 |
DE10161381A1 (de) * | 2001-12-14 | 2003-06-18 | Philips Intellectual Property | Patientendatenverarbeitungssystem und -verfahren |
US7520857B2 (en) * | 2002-06-07 | 2009-04-21 | Verathon Inc. | 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume |
JP4347319B2 (ja) * | 2002-06-25 | 2009-10-21 | アロカ株式会社 | 超音波診断装置 |
US7347821B2 (en) * | 2003-06-26 | 2008-03-25 | Koninklijke Philips Electronics N.V. | Adaptive processing of contrast enhanced ultrasonic diagnostic images |
WO2005030057A1 (en) | 2003-09-29 | 2005-04-07 | Koninklijke Philips Electronics, N.V. | Ultrasonic cardiac volume quantification |
US20050096543A1 (en) * | 2003-11-03 | 2005-05-05 | Jackson John I. | Motion tracking for medical imaging |
US7231076B2 (en) | 2004-06-30 | 2007-06-12 | Accuray, Inc. | ROI selection in image registration |
JP4901273B2 (ja) * | 2006-04-06 | 2012-03-21 | 株式会社東芝 | 超音波診断装置およびその画像処理プログラム |
US20080043024A1 (en) * | 2006-06-26 | 2008-02-21 | Siemens Corporate Research, Inc. | Method for reconstructing an object subject to a cone beam using a graphic processor unit (gpu) |
JP5536974B2 (ja) * | 2006-11-08 | 2014-07-02 | 株式会社東芝 | X線診断装置及び画像処理装置 |
US8923952B2 (en) * | 2006-12-11 | 2014-12-30 | Mayo Foundation For Medical Education And Research | System and method for quantitative molecular breast imaging |
RU2009128709A (ru) * | 2006-12-26 | 2011-02-10 | Конинклейке Филипс Электроникс Н.В. (Nl) | Медицинская система формирования изображений |
-
2010
- 2010-01-18 CN CN201080005093.7A patent/CN102292744B/zh not_active Expired - Fee Related
- 2010-01-18 JP JP2011547015A patent/JP5996870B2/ja not_active Expired - Fee Related
- 2010-01-18 EP EP10701935A patent/EP2389662A1/en not_active Ceased
- 2010-01-18 WO PCT/IB2010/050210 patent/WO2010084446A1/en active Application Filing
- 2010-01-18 US US13/145,920 patent/US8848989B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5797843A (en) * | 1992-11-03 | 1998-08-25 | Eastman Kodak Comapny | Enhancement of organ wall motion discrimination via use of superimposed organ images |
US6757414B1 (en) * | 2000-08-31 | 2004-06-29 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for segmentation of a left ventricular epicardium |
CN1969298A (zh) * | 2004-04-23 | 2007-05-23 | 美国西门子医疗解决公司 | 用于多对象体数据可视化的方法和系统 |
CN101110124A (zh) * | 2006-02-21 | 2008-01-23 | 美国西门子医疗解决公司 | 利用虚拟切割的背景中体积可视化的系统及方法 |
Non-Patent Citations (1)
Title |
---|
Non-grid image registration: theory and practice;W R CRUM et al.;《The British journal of radiology》;20140101;第77卷(第2期);S140-S153 * |
Also Published As
Publication number | Publication date |
---|---|
WO2010084446A1 (en) | 2010-07-29 |
US20110274326A1 (en) | 2011-11-10 |
JP2012515598A (ja) | 2012-07-12 |
US8848989B2 (en) | 2014-09-30 |
EP2389662A1 (en) | 2011-11-30 |
CN102292744A (zh) | 2011-12-21 |
JP5996870B2 (ja) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102292744B (zh) | 心脏图像处理和分析 | |
CN102473300B (zh) | 多模态乳房成像 | |
US11672505B2 (en) | Correcting probe induced deformation in an ultrasound fusing imaging system | |
US8837789B2 (en) | Systems, methods, apparatuses, and computer program products for computer aided lung nodule detection in chest tomosynthesis images | |
US9741131B2 (en) | Anatomy aware articulated registration for image segmentation | |
US10803354B2 (en) | Cross-modality image synthesis | |
US20070167784A1 (en) | Real-time Elastic Registration to Determine Temporal Evolution of Internal Tissues for Image-Guided Interventions | |
AU2006282500A1 (en) | Image processing method, image processing program, and image processing device | |
US20010007593A1 (en) | Method and unit for displaying images | |
JP2012515598A5 (zh) | ||
CN117813055A (zh) | 用于从快速spect扫描和ct图像合成spect图像的多模态和多尺度特征聚合 | |
Lappas et al. | Automatic contouring of normal tissues with deep learning for preclinical radiation studies | |
Mezheritsky et al. | Population-based 3D respiratory motion modelling from convolutional autoencoders for 2D ultrasound-guided radiotherapy | |
WO2011041473A1 (en) | Medical image analysis system for anatomical images subject to deformation and related methods | |
EP2483866A1 (en) | Medical image analysis system using n-way belief propagation for anatomical images subject to deformation and related methods | |
US9839404B2 (en) | Image data Z-axis coverage extension for tissue dose estimation | |
EP2483864A1 (en) | Medical image analysis system for displaying anatomical images subject to deformation and related methods | |
EP4364665A1 (en) | Disease label creation device, method, and program, learning device, and disease detection model | |
EP4128145B1 (en) | Combining angiographic information with fluoroscopic images | |
Pansiot et al. | CBCT of a moving sample from X-rays and multiple videos | |
Kan | Organ segmentation of pediatric computed tomography (ct) with generative adversarial networks | |
CN118985004A (zh) | 根据展开扫描图像的自动肋骨骨折检测 | |
Mazaheri et al. | Research Article Hybrid Pixel-Based Method for Cardiac Ultrasound Fusion Based on Integration of PCA and DWT |
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: 20151125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |