CN103083046A - 用于成像的方法和电子计算装置以及成像*** - Google Patents

用于成像的方法和电子计算装置以及成像*** Download PDF

Info

Publication number
CN103083046A
CN103083046A CN2012104439540A CN201210443954A CN103083046A CN 103083046 A CN103083046 A CN 103083046A CN 2012104439540 A CN2012104439540 A CN 2012104439540A CN 201210443954 A CN201210443954 A CN 201210443954A CN 103083046 A CN103083046 A CN 103083046A
Authority
CN
China
Prior art keywords
image data
data group
acoustic pressure
magnetic resonance
tissue
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.)
Pending
Application number
CN2012104439540A
Other languages
English (en)
Inventor
P.格罗斯
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.)
Siemens AG
Original Assignee
Siemens AG
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 Siemens AG filed Critical Siemens AG
Publication of CN103083046A publication Critical patent/CN103083046A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5261Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0035Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

在成像中通过应用超声(23)来产生声压(24),从而引起检查对象(4)的组织(22)的形变。依据至少一个超声图像数据组来确定第一数据,该第一数据位置分辨地表示作为对所述声压(24)的反应的组织(22)的形变。依据至少一个磁共振图像数据组来确定第二数据,该第二数据位置分辨地表示作为对所述声压(24)的反应的组织(22)的形变。依据所述第一数据和所述第二数据的比较来联合配准所述至少一个超声图像数据组和所述至少一个磁共振图像数据组。

Description

用于成像的方法和电子计算装置以及成像***
技术领域
本发明涉及用于成像的方法和装置。本发明特别地涉及这样的方法和装置,通过其能够配准利用不同成像模态采集的图像。
背景技术
在医学技术中采用大量不同的成像方法。属于这些方法的例如有计算机断层造影、磁共振断层造影、超声以及正电子发射断层造影。这些成像模态位置分辨地提供关于检查对象的检查区域的不同物理特性的信息。由于不同的成像模态能够提供关于检查区域的不同特性的信息,期望利用不同的成像模态采集图像数据。
在利用多个成像模态对检查区域进行成像时,要求将不同获得的图像彼此配准。联合配准(Koregistrierung)要求以合适的方式对齐并且正确缩放利用不同成像模态所采集的图像。由于基于不同成像模态的不同物理过程,在利用不同成像模态采集的图像中会存在不同的失真。例如,利用磁共振成像装置所采集的图像会由于梯度场的非线性而具有失真,而在利用超声成像装置对同样区域所采集的图像中不存在该失真。
为了联合配准利用不同成像模态采集的图像,可以在不同图像中识别用户定义的解剖学标志。例如可以使用骨骼或导管作为标志。但这种标志在骨骼的情况下是与位置相连的或者在导管的情况下是仅能介入定位的并且取决于患者的解剖结构。用户定义的规定可以是耗费时间的。如果用户不能可靠地识别解剖学标志,则依据用户的经验而用户定义的规定会是不精确的或易于产生错误的。难点在于,确定足够多数量的配准点,以便能够在三维中进行非线性联合配准。在图像中能够良好识别的解剖学标志的数量对于这种目的来说太少。
US 2010/0087727A1和DE 10 2008 047 644 A1分别描述了一种方法和一种组合***,其中依据热信息对利用不同成像模态所采集的图像进行配准。这一点要求检查对象的可相应探测的发热。
还存在用于如下成像的方法和装置的需求,其中能够可靠地联合配准利用不同成像模态所采集的图像,以便能够彼此融合在不同图像数据中包含的信息。
发明内容
因此,本发明要解决的技术问题是,提供用于成像的一种方法、一种装置和一种***,其允许可靠地联合配准利用不同成像模态所采集的图像数据组。本发明要解决的技术问题是,实现借助多个配准点的配准,而不需要测量温度分布图。
根据实施例,给出一种用于成像的方法。在此,通过应用超声来产生声压,从而引起检查对象的组织形变。利用超声成像装置拍摄至少一个超声图像数据组。依据该至少一个超声图像数据组来确定第一数据,其中第一数据位置分辨地表示作为对声压的反应的组织形变。利用磁共振成像装置拍摄至少一个磁共振图像数据组。依据该至少一个磁共振图像数据组来确定第二数据,其中第二数据位置分辨地表示作为对声压的反应的组织形变。依据第一数据和第二数据的比较来联合配准该至少一个超声图像数据组和该至少一个磁共振图像数据组。
在该方法中,将超声成像和磁共振成像组合。超声和核自旋共振校验不同的物质特性。通过这种组合成像可以获得关于检查对象对超声的响应以及对核自旋动态的响应的信息。超声和核自旋共振仅具有小的或没有相互影响。这一点有助于组合成像。
在该方法中,所引起的形变或检查对象的组织的反应被作为用于配准的标记。由此可以由无需介入干预的成像***并且通过合适地产生声压来产生局部标志。由此实现了高的灵活性。可以这样设置标志的数量和密度,使得可以以期望的精度来进行非线性的联合配准。
可以这样产生声压,使得引起组织的局部形变。由此有助于在利用超声成像和利用磁共振成像所获得的数据中识别合适的标志。可以时间顺序地或时间并行地引起组织的局部形变。每次形变可以具有比利用成像模态成像的检查对象区域更小的,尤其是极其小的延伸。
这样确定第一数据,使得其表示组织的取决于位置的刚性(Steifigkeit)或粘性可以简单并可靠地确定这样的第一数据,例如通过使用超声检查弹性图(sonographische Elastographie)或剪切波成像(Scherwellenbildgebung)来确定第一数据。组织的刚性或粘性量化了作为对声压的反应的组织形变并且可以将其用于联合配准。
也可以这样确定第一数据,使得其表示在至少一个超声图像数据组的坐标系中的作为对声压的反应的组织取决于位置的偏移。为此,例如可以依据时间顺序地所采集的两个超声图像数据组来取决于位置地确定组织的偏移。这样的形变图有助于对利用不同成像模态所采集的图像的比较。
可以这样确定第二数据,使得其表示在至少一个磁共振图像数据组的坐标系中的组织的取决于位置的刚性或粘性。例如,可以通过组合用于产生声学辐射压力的超声与磁共振成像来采集组织的刚性或粘性。为此可以采用特别是扩散灵敏的磁共振成像序列,以便可以显示粘性的脉冲传输。组织的刚性或粘性量化了作为对声压的反应的组织的形变并且可以将其用于联合配准。
可以这样确定第二数据,使得其表示在至少一个磁共振图像数据组的坐标系中的作为对声压的反应的组织取决于位置的偏移。为此,例如可以依据时间顺序地所采集的两个磁共振图像数据组来取决于位置地确定组织的偏移。这样的形变图有助于对利用不同成像模态所采集的图像的比较。
依据第一数据和第二数据的比较可以确定多个,特别是多于三个配准点。联合配准可以是非线性的联合配准。联合配准可以包括至少一个磁共振图像数据组和/或至少一个超声图像数据组的非线性的变形或压缩。联合配准可以包括至少一个磁共振图像数据组和/或至少一个超声图像数据组的剪切。
可以依据配准点的所期望的数量和/或密度自动地通过应用超声来产生用于引起组织局部形变的声压。用于成像的***可以自动地选择配准点的数量和/或密度,并且通过应用超声相应地自动控制声压的产生。通过这种方式可以自动地确保对于非线性的联合配准存在足够数量和/或密度的配准点。
可以依据磁共振成像装置的特性通过应用超声来产生用于引起组织局部形变的声压。可以依据在磁共振成像装置的梯度场中的非线性通过应用超声来产生用于引起组织局部形变的声压。通过这种方式可以确保在其中至少一个磁共振图像数据组具有较强失真的区域中比在其中梯度场显示线性特性的区域中产生更多数量的用于配准的标志。
可以时间上错开地拍摄第一磁共振图像数据组和第二磁共振图像数据组。可以依据第一磁共振图像数据组和第二磁共振图像数据组的比较来确定第二数据。可以在产生声压之前拍摄第一磁共振图像数据组,并且可以在产生声压期间或者产生声压之后拍摄第二磁共振图像数据组。由此可以这样确定第二数据,使得其表示作为对声压的反应的组织取决于位置的偏移。
在该方法中,也可以在产生用于引起组织形变的声压期间仅拍摄一个磁共振图像数据组。可以通过扩散灵敏的磁共振成像序列来采集该磁共振图像数据组。通过这种方式可以在没有多次成像相同体素的情况下位置分辨地确定组织的刚性或粘性。其可以与第一数据相比较,以便进行联合配准。通过由多个磁共振图像数据组随后形成差由此可以避免多次扫描相同体素。
在该方法中,可以在产生声压之前拍摄第一超声图像数据组。可以在产生引起组织局部形变的声压之后或在产生声压期间拍摄第二超声图像数据组。可以依据第一超声图像数据组和第二超声图像数据组的比较来确定第一数据。由此可以这样确定第一数据,使得其表示作为对声压的反应的组织取决于位置的偏移以及由此产生的局部形变。
可以取决于为拍摄至少一个磁共振图像数据组所选择的检查对象的区域产生声压。例如,可以依据磁共振成像装置的运行参数来规定由超声转换器为产生声压而产生的超声的方向和/或焦点。可以这样选择超声的方向,使得其相应于如下方向,在该方向上磁共振成像序列是扩散灵敏的。在利用磁共振成像装置顺序读取检查对象的多个层的情况下,可以这样设置由超声转换器所以产生的声学辐射的焦点,使得其分别位于由磁共振成像装置读取的层中和/或距离利用层共振成像装置读取的层比阈值更小。
根据本发明的另一个实施例,给出一种计算机程序,其包括命令序列,该命令序列在执行中通过具有磁共振成像装置和超声成像装置的成像***的电子计算装置促使该成像***执行根据本发明的方面或实施例的方法。计算机程序例如可以加载到成像***的分析和控制计算机的存储器中。计算机程序可以作为源代码或者作为组合的命令序列呈现。通过计算机程序可以将成像***程序地构造为用于执行该方法。
根据本发明的另一个实施例,给出一种数据载体,在该数据载体上存储了包含命令序列的计算机程序,该计算机程序在执行中通过具有磁共振成像装置和超声成像装置的成像***的电子计算装置促使该成像***执行根据本发明的方面或实施例的方法。数据载体例如可以是CD-ROM、DVD、磁带、闪存或USB棒,在其上存储了作为电子可读的控制信息的计算机程序。如果从数据载体读取控制信息并且将其由成像***的计算装置执行,则可以执行根据成像***的不同方面或实施方式的方法。
根据本发明的另一个实施例,给出一种用于成像的电子计算装置。该电子计算装置包括用于接收至少一个超声图像数据组的第一接口和用于接收至少一个磁共振图像数据组的第二接口。该电子计算装置包括处理器单元,其被构造为,依据至少一个超声图像数据组这样确定第一数据,使得第一数据位置分辨地表示作为对声压的反应的检查对象的组织形变。该处理器单元被构造为,依据至少一个磁共振图像数据组这样确定第二数据,使得第二数据位置分辨地表示作为对声压的反应的检查对象的组织形变。该处理器单元被构造为,依据第一数据和第二数据的比较来执行至少一个超声图像数据组和至少一个磁共振图像数据组的联合配准。
根据本发明的另一个实施例,给出一种成像***,该成像***包括超声成像装置、磁共振成像装置和根据实施例的电子计算装置,该电子计算装置与超声成像装置和磁共振成像装置相耦合。
该电子计算装置可以被构造为,用于控制超声转换器,以便通过应用超声来产生声压。通过这种方式该电子计算装置可以控制动态地设置由于组织形变而产生的标志。
根据实施例的电子计算装置和成像***的实施方式以及由此达到的效果相应于根据实施例的实施方式和方法的效果。
成像***可以被构造为用于执行根据实施例的方法。
附图说明
本发明的上面描述的特性、特征和优点以及如其实现的那样的方式结合下面描述的结合附图更详细地解释的实施例可以更清楚和更明显地理解。
图1示出了成像***。
图2和图3示出了用于解释成像***工作原理的检查对象的检查区域。
图4示出了第一数据和第二数据,它们分别表示作为对声压的反应的组织形变。
图5示出了根据实施例的方法流程图。
图6示出了多个超声图像数据组和从中获得的第一数据。
图7示出了多个磁共振图像数据组和从中获得的第二数据。
图8示出了根据另一个实施例的方法的流程图。
图9示出了根据另一个实施例的方法的流程图。
图10示出了根据另一个实施例的方法的流程图。
具体实施方式
图1示出了根据实施例的成像***1。成像***1包括磁共振(MR)成像装置2和超声***5。MR成像装置2具有支撑装置3,例如检查卧榻,在该检查卧榻上安置检查对象4。超声***5可以具有布置在共同壳体中的超声(US)成像装置6和US治疗装置7。US治疗装置7被构造为用于产生具有振幅的超声脉冲,从而产生的声压导致检查对象4的软组织形变。可以这样选择振幅,使得患者的软组织具有一微米或几微米的偏移。其组织的偏移可以在一微米和二十微米之间。在运行中采用超声脉冲的情况下可以实现软组织的这种偏移,其中不超过常规声学能量。作为与US成像装置6一起布置在壳体中的US治疗装置7的替换或附加,超声***5可以具有独立的US转换器8,用来将超声耦合到检查对象中和/或采集超声。
MR成像装置2和超声***5还具有出于清楚起见而未示出的本身公知的元件。例如MR成像装置2具有至少一个高频线圈和梯度***。在一种实施方式中,超声***5的组件可以固定在机械臂上,该机械臂可以相对于MR成像装置2运动。通过这种机械臂可以自动地或用户定义地定位将超声耦合到检查对象4中的超声***5的组件。
成像***1具有电子计算装置10,该电子计算装置10与MR成像装置2和超声***5耦合。电子计算装置10满足分析功能。如还要详细描述的那样,电子计算装置10从利用US成像装置6采集的至少一个US图像数据组中来确定第一数据,该第一数据表示作为对利用超声***产生的声压的反应的组织的局部形变。第一数据可以包括在US图像数据组中的至少一个点的位置,其相对于其应用声压之前的位置具有最大偏移。第一数据也可以包括多个局部偏移或形变最大值的位置。电子计算装置10还被构造为从利用MR成像装置2采集的至少一个MR图像数据组中来确定第二数据,该第二数据表示作为对利用超声***产生的声压的反应的组织的局部形变。第二数据可以包括在MR图像数据组中的至少一个点的位置,其相对于其应用声压之前的位置具有最大偏移。第二数据也可以包括在MR图像数据组中的多个局部偏移或形变最大值的位置。电子计算装置10被构造为将这样识别的配准点彼此关联。根据通过这种方式识别的配准点来联合配准MR图像数据组和US图像数据组。
此外,电子计算装置10可以满足控制功能。可以设置单独的控制和分析计算机,其在组合中满足电子计算装置10的功能。电子计算装置10可以这样控制超声***5,使得产生导致检查对象的组织的局部形变的声压。可以这样产生声压,使得产生一个超声脉冲或多个短的超声脉冲来引起局部形变。超声脉冲可以具有少于500微秒的持续时间,例如从50至150微秒的间隔的持续时间。通过使用导致组织的局部偏移的脉冲可以保持小的能量输入和组织加热。超声脉冲可以在空间上聚焦。直至20微米范围的偏移,例如从1至20微米的偏移足以在US图像数据组中以及在MR图像数据组中被可靠地识别。利用由超声***使用的常规能量可以实现这样的偏移。可以这样产生为引起局部形变而产生的超声脉冲,使得在击中预定组织时的侧向延伸小于阈值。
电子计算装置10可以被构造为用于自动地确定配准点的数量和/或密度。电子计算装置10可以自动地确定数量和/或密度,其中配准点不必是均匀分布的。例如电子计算装置10可以根据利用MR成像装置2或利用US成像装置5采集的概览图自动地识别软组织。
通过声压可以在利用MR成像装置2和US成像装置5并行进行的图像拍摄中分别动态地产生其中引起了软组织形变的标志,该形变既由MR成像装置2也由US成像装置5来成像。
电子计算装置10可以控制超声***5依据MR成像装置2和/或超声***5的特性来产生组织的局部形变。例如MR成像装置2的梯度场可以具有非线性,其在成像范围的一些区域中更明显。在这种情况下电子计算装置10可以这样控制超声***5,使得作为用于联合配准的标志的组织的局部形变以足够的数量和/或密度设置在这些区域中,以便允许执行非线性的联合配准。例如,可以不进行MR图像数据的线性缩放而是进行非线性缩放,即非线性的压缩或变形。可以确定非线性缩放的参数,方法是通过局部形变设置多于三个配准点。替换地或附加地,为了联合配准可以剪切MR图像数据。可以确定剪切的参数,方法是通过局部形变设置多于三个配准点。
电子计算装置10具有第一接口16,经由该第一接口16处理器单元11从超声***5接收数据。如果电子计算装置10还执行控制功能,则经由第一接口16还可以向超声***5输出控制命令。通过控制命令可以控制经由超声来产生声压用于设置形变标志。电子计算装置10具有第二接口12,经由该第二接口12处理器单元11从MR成像装置2接收数据。如果电子计算装置10还执行控制功能,则经由第二接口12还可以向MR成像装置2输出控制命令。电子计算装置10被构造为,将产生用于设置形变标志的声压与通过MR成像装置2的成像时间上协调地进行。可以依据用于控制MR成像装置2的控制信号来控制超声***5。超声***5可以被这样控制,即通过为组织局部形变而产生的脉冲在检查对象的层中或部分体积中引起形变,该层或该部分体积通过MR成像装置2来成像。处理器单元11可以被程序技术地构造为,使得电子计算装置10执行在此描述的不同分析和控制功能。
图2和图3表示了采用声压来动态地产生作为用于联合配准的标志的形变。图2和图3分别示出了检查对象的检查区域21。具有微小刚性的组织22由于入射超声脉冲23而形变。通过粘性脉冲传输产生由入射到组织的声压24导致的力。超声脉冲23产生通过声学辐射产生的力。该通过声学辐射产生的力引起组织22的局部形变。利用MR成像以及US成像来采集通过声学辐射产生的力或由此引起的作为对其反应的组织形变。在MR成像数据组的坐标系中以及在US图像数据组的坐标系中识别形变组织的位置。由此提供配准点。该过程可以时间上并行地或时间上顺序地重复,其中在检查对象的不同的点处引起多个局部形变。
可以这样使用由声压24引起的形变,使得其由MR成像装置成像。如果MR成像关于特定的时间点空间选择性地成像层或部分体积25,则可以这样产生超声脉冲23,使得在层或部分体积25中引起形变。如果MR成像关于另一个时间点成像另一个层或另一个部分体积26,则可以这样产生超声脉冲23,使得在另一个层或另一个部分体积26中引起形变。
图4表示了US图像数据和MR图像数据的联合配准。从多组US图像数据中计算地获得第一形变图31。第一形变图31位置分辨地表示作为对局部起作用的力的反应的组织的位置偏移,该局部起作用的力通过超声产生。第一形变图31例如可以通过两个US图像数据组的比较获得,这两US图像数据组中的一个在由于声压引起的形变之前拍摄,而另一个在其期间或者在其之后拍摄。形变轮廓32位置分辨地示出了作为对声压的反应的组织偏移,如其在US图像数据组中成像的那样。最大偏移出现在图像点33处。
从多组MR图像数据中计算地获得第二形变图36。第二形变图36位置分辨地表示作为对局部起作用的力的反应的组织的位置偏移,该局部起作用的力通过超声产生。第二形变图36例如可以通过两个MR图像数据组的比较获得,这两个MR图像数据组中的一个在由于声压引起的形变之前拍摄,而另一个在其期间或者在其之后拍摄。形变轮廓37位置分辨地示出了作为对声压的反应的组织偏移,如其在MR图像数据组中成像的那样。最大偏移出现在图像点38处。
通过US图像数据中的图像点33和MR图像数据中的图像点38提供配准点。如果确定了形变轮廓,则可以根据形变轮廓进行联合配准。替换地或附加地,可以重复设置作为用于联合配准的标志的局部形变,以便获得其它配准点。
可以采用不同的措施来提高联合配准的鲁棒性。例如可以采用阈值比较。在从US图像数据组和MR图像数据组中来确定形变图的情况下可以分别考虑组织偏移超过阈值的区域。阈值可以取决于如下而自动地确定:置于声压下的组织是相对较软还是相对较硬。
图5示出了根据实施例的方法40的流程图。方法40可以由成像***1实施。在方法40中采集两个US图像数据组,以便在US图像数据组的坐标系中确定响应于所产生的声压的组织的偏移。在方法40中采集两个MR图像数据组,以便在MR图像数据组的坐标系中确定响应于所产生的声压的组织的偏移。
在41中执行第一US图像采集。在41中可以执行B图像超声检查法(B-Bild-Sonographie)。由此可以获得第一检查图像。在43中产生超声脉冲用来通过声压引起组织的局部形变。在44中进行第二US图像采集。第二US图像采集可以在43中产生声压期间或者在43中产生声压之后进行。第二US图像采集然后在超声脉冲结束之后的时间段内执行,该时间比由于超声脉冲而偏移的组织的弛豫时间要小。通过这种方式确保了组织没有弛豫回到其静止位置。44中的第二图像采集可以通过利用其进行超声检查跟踪的超声辐射执行。44中的第二US图像采集可以通过常规的超声检查技术来实施。
与US图像采集并行地还进行MR图像采集。在42中进行第一MR图像采集。第一MR图像采集可以是利用MR进行的概览图采集。在45中进行第二MR图像采集。第二MR图像采集可以在43中产生声压期间或者在43中产生声压之后进行。第二MR图像采集与第二US图像采集时间上重叠地进行。第二MR图像采集同样在超声脉冲结束之后的时间段内执行,该时间比由于超声脉冲而偏移的组织的弛豫时间要小。通过这种方式确保了组织没有弛豫回到其静止位置。
在46中依据在第一US图像采集中拍摄的第一US图像数据组和在第二US图像采集中拍摄的第二US图像数据组的比较来计算第一数据。第一数据位置分辨地给出在US图像数据组的坐标系中的组织的偏移。为此例如可以计算出灰阶图(Graustufenbild),该灰阶图位置分辨地描述了在US图像数据组中成像的组织的偏移。在47中依据在第一MR图像采集中拍摄的第一MR图像数据组和在第二MR图像采集中拍摄的第二MR图像数据组的比较来计算第二数据。第二数据位置分辨地给出在MR图像数据组的坐标系中的组织的偏移。为此例如可以计算出灰阶图,该灰阶图位置分辨地描述了在MR图像数据组中成像的组织偏移。在MR图像数据组和US图像数据组的坐标系中这样确定的形变图可以用于联合配准。例如,按照第一数据在US图像数据组中具有最大偏移的点和按照第二数据在MR图像数据组中具有最大偏移的点可以被用作配准点。如果例如利用第一数据和第二数据确定取决于位置的形变轮廓,则也可以将第一数据和第二数据中的多个点用作配准点。
在48中检查是否已经确定了足够数量的配准点。如果配准点还没有达到足够数量,则通过在其上施加声压可以动态地在另一个点处设置形变标志。可以重复在US图像数据和MR图像数据中的形变确定。在48中检查时可以确定是否已经达到预定数量的配准点,这根据用于联合配准的参数数量来选择。对于非线性的联合配准必须存在多于三对在利用US成像和MR成像所采集的图像数据组中的配准点。
在49中进行联合配准。动态地由***通过应用声压而产生的足够数量的局部组织形变在此作为用于联合配准的标志来使用。联合配准可以是非线性的联合配准。可以非线性地伸缩和/或剪切至少一个图像数据组,例如至少一个MR图像数据组。
图6示例性地示出了US图像数据组51、52及其在方法40中的处理。在通过超声脉冲引起组织的局部形变之前拍摄的第一US图像数据组51中包含组织的图像55。在将声压施加到组织来引起局部形变之后或期间采集的第二US图像数据组52中包含组织的图像56。局部应用的声压导致相对于原始图像55的偏移。可以计算地获得第一数据53,该第一数据量化了偏移。例如可以计算差图像。可以采用更复杂的技术来识别在US图像数据组的坐标系中的片段57,其中通过超声脉冲来偏移组织。
图7示例性地示出了MR图像数据组61、62及其在方法40中的处理。在通过超声脉冲引起组织的局部形变之前拍摄的第一MR图像数据组61中包含组织的图像65。在将声压施加到组织来引起局部形变之后或期间采集的第二MR图像数据组62中包含组织图像66。局部应用的声压导致相对于原始图像65的偏移。可以计算地获得第二数据63。该第二数据量化了偏移。例如可以计算差图像。可以采用更复杂的技术来识别在MR图像数据组的坐标系中的片段67,其中通过超声脉冲来偏移组织。
在利用不同成像模态采集的图像中的计算确定的片段57、67可以被用于联合配准。
对由于局部施加声压而偏移的组织的识别不仅可以通过采集多个图像数据组以及多个图像数据组的比较进行。还可以在唯一的成像中量化组织对声压的反应,其例如测量组织的刚性或粘性。在超声成像中这一点可以通过超声检查弹性图或剪切波成像来实现。在此,例如可以采用在DE 10 2009 033286 A1中描述的方法。也可以在超声图像中直接成像超声焦点,如其在US5,769,790中描述的那样。在MR成像中例如通过使用扩散灵敏的MR脉冲序列通过沿着扩散灵敏的方向入射超声来实现刚性或粘性的量化。例如在WO 2008/067905A1中描述了在此可以采用的方法。
在实施例中,可以利用至少一个成像模态来采集刚性图或粘性图。刚性图或粘性图位置分辨地量化在利用相应成像模态所采集的图像数据组中的刚性或粘性。可以将利用一个成像模态所采集的刚性图或粘性图与利用另一个成像模态所采集的刚性图或粘性图相比较。替换地或附加地,可以将利用一个成像模态所采集的刚性图或粘性图与从利用另一个成像模态所采集的两个图像数据组中计算出的形变图相比较。参考图8至图10描述了这些方法。
图8示出了根据另一种实施例的方法70的流程图。方法70可以由成像***1实施。
在方法70中,在72中采集US图像数据组。在此可以产生导致组织局部形变的声压。US图像数据组的采集可以通过超声检查弹性图或剪切成像进行。也可以采用在US 5,769,790中描述的方法来超声检查地成像超声脉冲的焦点。依据在72中所采集的US图像数据组可以在74中计算关于组织的刚性和粘性的位置分辨的信息。采用关于刚性和粘性的位置分辨的信息作为用于联合配准的第一数据。
在方法70中,在通过超声脉冲引起组织局部形变之前,在71中采集第一MR图像数据组。在73中采集第二MR图像数据组。可以在由于US图像采集72而产生声压期间或者在由于US图像采集72而产生声压之后进行第二MR图像采集。在由于声压而形变的组织弛豫回到其静止位置之前进行第二MR图像采集。在75中可以由第一MR图像数据组和第二MR图像数据组计算第二数据,该第二数据量化了在MR图像数据组的坐标系中的组织的偏移。这一点可以如参考图5至图7所描述的那样进行。
在76中依据在74中确定的第一数据与在75中确定的第二数据的比较进行联合配准。
方法70可以包含是否存在足够数量的配准点的检查。可以选择性地重复方法步骤,以便通过局部形变产生足够数量的配准点。
图9示出了根据另一种实施例的方法80的流程图。方法80可以由成像***1实施。
在方法80中,在81中采集第一US图像数据组。在82中产生引起组织局部形变的超声脉冲。在84中采集第二US图像数据组。在85中可以从第一US图像数据组和第二US图像数据组中计算第一数据,该第一数据量化了在US图像数据组的坐标系中的组织的偏移。采用该数据作为用于联合配准的第一数据。步骤81、82、84和85可以如参考图5所描述的那样执行。
与在82中产生超声脉冲并行地,在83中执行MR图像采集。在此可以是扩散灵敏的MR图像采集。在82中可以沿着如下方向入射超声脉冲,沿着该方向MR图像采集是扩散灵敏的。由于超声脉冲引起的粘性脉冲传输允许在MR图像采集中对组织的刚性或粘性的成像。在86中可以从MR图像数据中确定刚性或粘性作为在MR图像数据的坐标系中的位置的函数。采用该位置分辨的数据作为用于联合配准的第二数据。MR数据采集和分析可以如在WO 2008/067905A1中所描述的那样进行。
在87中依据在85中确定的第一数据与在86中确定的第二数据的比较进行联合配准。例如可以根据利用MR成像确定的刚性沿着传播方向的改变来确定超声脉冲在哪里击中具有最小刚性或粘性的组织。可以将该位置与在其中US图像数据组具有最大组织偏移的位置相比较。
方法80可以包括是否存在足够数量的配准点的检查。可以选择性地重复方法步骤,以便通过局部形变产生足够数量的配准点。
图10示出了根据另一种实施例的方法90的流程图。方法90可以由成像***1实施。
在方法90中,在91中进行US图像采集,通过其获得关于组织的刚性或粘性的位置分辨的信息。在此,采用导致组织局部形变的超声脉冲。可以通过超声检查弹性图或剪切波成像进行US图像数据组的采集。依据在91中采集的US图像数据组可以在93中计算关于组织的刚性和粘性的位置分辨的信息。采用关于刚性和粘性的位置分辨的信息作为用于联合配准的第一数据,该刚性和粘性量化了作为对辐射压力的反应的组织的形变。
与US图像采集并行地在92中进行MR图像采集。在此可以是扩散灵敏的MR图像采集。在US图像采集中采用的超声脉冲可以沿着如下方向入射,沿着该方向MR图像采集是扩散灵敏的。由于超声脉冲引起的粘性脉冲传输允许在MR图像采集中对组织的刚性或粘性的成像。在94中可以从MR图像数据中确定刚性或粘性作为在MR图像数据的坐标系中的位置的函数。采用该位置分辨的数据作为用于联合配准的第二数据。MR数据采集和分析可以如在WO 2008/067905A1中所描述的那样进行。
在95中依据在93中确定的第一数据与在94中确定的第二数据的比较进行联合配准。例如可以将利用91中US图像采集和94中MR图像采集确定的刚性或粘性的最大值作为配准点使用。
方法90可以包括对是否存在足够数量的配准点的检查。可以选择性地重复方法步骤,以便通过局部形变产生足够数量的配准点。
虽然通过优选的实施例详细解释和描述了本发明,但本发明不受到所公开的示例的限制并且可以由专业人员从中导出其它方案,而不脱离本发明的保护范围。

Claims (19)

1.一种用于成像的方法,包括:
-通过应用超声(23)来产生声压(24),从而引起检查对象(4)的组织(22)的形变,
-利用超声成像装置(6)拍摄至少一个超声图像数据组(51,52),
-依据所述至少一个超声图像数据组(51,52)来确定第一数据(31;53),其中,所述第一数据(31;53)位置分辨地表示作为对所述声压(24)的反应的组织(22)的形变(32;57),
-利用磁共振成像装置(2)拍摄至少一个磁共振图像数据组(61,62),
-依据所述至少一个磁共振图像数据组(61,62)来确定第二数据(36;63),其中,所述第二数据(36;63)位置分辨地表示作为对所述声压(24)的反应的组织(22)的形变(37;67),并且
-依据所述第一数据(31;53)和所述第二数据(36;63)的比较来联合配准所述至少一个超声图像数据组(51,52)和所述至少一个磁共振图像数据组(61,62)。
2.根据权利要求1所述的方法,其中,这样产生所述声压(24),使得引起所述组织(22)的局部形变。
3.根据上述权利要求中任一项所述的方法,其中,这样确定第一数据,使得其表示所述组织(22)的取决于位置的刚性或粘性。
4.根据权利要求1或权利要求2所述的方法,其中,这样确定所述第一数据(31;53),使得其表示作为对所述声压(24)的反应的所述组织(22)取决于位置的偏移。
5.根据上述权利要求中任一项所述的方法,其中,这样确定所述第二数据,使得其表示所述组织(22)的取决于位置的刚性或粘性。
6.根据权利要求1至4中任一项所述的方法,其中,这样确定所述第二数据(36;63),使得其表示作为对所述声压(24)的反应的所述组织(22)取决于位置的偏移。
7.根据上述权利要求中任一项所述的方法,
其中,依据所述第一数据(31;53)和所述第二数据(36;63)的比较确定多于三个配准点(33,38);并且
其中,联合配准包括:
-依据所述第一数据(31;53)和所述第二数据(36;63)的比较实施非线性的联合配准。
8.根据上述权利要求中任一项所述的方法,其中,所述至少一个磁共振图像数据组(61,62)的拍摄包括:
-时间上错开地拍摄第一磁共振图像数据组(61)和第二磁共振图像数据组(62)。
9.根据权利要求8所述的方法,其中,依据所述第一磁共振图像数据组(61)和所述第二磁共振图像数据组(62)的比较来确定第二数据(36;63)。
10.根据权利要求8或权利要求9所述的方法,
其中,在产生声压(24)之前拍摄所述第一磁共振图像数据组(61),并且
其中,在产生声压(24)期间或者产生声压(24)之后拍摄所述第二磁共振图像数据组(62)。
11.根据权利要求1至7中任一项所述的方法,
其中,所述至少一个磁共振图像数据组的拍摄包括:
-在产生声压(24)期间仅拍摄一个磁共振图像数据组;
其中,依据该磁共振图像数据组确定第二数据。
12.根据上述权利要求中任一项所述的方法,其中,所述至少一个超声图像数据组(51,52)的拍摄包括:
-在产生声压(24)之前拍摄第一超声图像数据组(51),并且
-在产生声压(24)之后拍摄第二超声图像数据组(52)。
13.根据权利要求1至11中任一项所述的方法,其中,所述至少一个超声图像数据组(51,52)的拍摄包括:
-在产生声压(24)之前拍摄第一超声图像数据组(51),并且
-在产生声压(24)期间拍摄第二超声图像数据组(52)。
14.根据权利要求12或权利要求13所述的方法,其中,依据所述第一超声图像数据组(51)和所述第二超声图像数据组(52)的比较来确定第一数据(31;53)。
15.根据上述权利要求中任一项所述的方法,其中,取决于为拍摄所述至少一个磁共振图像数据组(61,62)所选择的检查对象(4)的区域(25,26)执行声压(24)的产生。
16.一种用于成像的电子计算装置,包括:
用于接收至少一个超声图像数据组(51,52)的第一接口(16);
用于接收至少一个磁共振图像数据组(61,62)的第二接口(12);
处理器单元(11),其被构造为,
-依据所述至少一个超声图像数据组(51,52)这样确定第一数据(31;53),使得所述第一数据(31;53)位置分辨地表示作为对声压(24)的反应的检查对象(4)的组织(22)的形变(32,57),
-依据所述至少一个磁共振图像数据组(61,62)这样确定第二数据(36;63),使得所述第二数据(36;63)位置分辨地表示作为对声压(24)的反应的检查对象(4)的组织(22)的形变(37;67),并且
-依据所述第一数据(31;53)和所述第二数据(36;63)的比较来执行所述至少一个超声图像数据组(51,52)和所述至少一个磁共振图像数据组(61,62)的联合配准。
17.一种成像***,包括:
超声成像装置(6);
磁共振成像装置(2);和
根据权利要求16所述的电子计算装置,该电子计算装置与所述超声成像装置(6)和磁共振成像装置(2)相耦合。
18.根据权利要求17所述的成像***,其中,所述电子计算装置(10)被构造为,用于控制超声转换器(8)来产生声压(24)。
19.根据权利要求17或权利要求18所述的成像***,所述成像***被构造为用于执行根据权利要求2至15中任一项所述的方法。
CN2012104439540A 2011-11-08 2012-11-08 用于成像的方法和电子计算装置以及成像*** Pending CN103083046A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102011085894.6A DE102011085894B4 (de) 2011-11-08 2011-11-08 Verfahren und elektronische Recheneinrichtung zur Bildgebung und Bildgebungssystem
DE102011085894.6 2011-11-08

Publications (1)

Publication Number Publication Date
CN103083046A true CN103083046A (zh) 2013-05-08

Family

ID=48128840

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012104439540A Pending CN103083046A (zh) 2011-11-08 2012-11-08 用于成像的方法和电子计算装置以及成像***

Country Status (3)

Country Link
US (1) US8897855B2 (zh)
CN (1) CN103083046A (zh)
DE (1) DE102011085894B4 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105338905A (zh) * 2013-06-26 2016-02-17 皇家飞利浦有限公司 用于多模态组织分类的方法和***
CN111493931A (zh) * 2019-08-01 2020-08-07 深圳迈瑞生物医疗电子股份有限公司 一种超声成像方法及设备、计算机可读存储介质

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3013244B1 (en) * 2013-06-26 2019-01-16 Koninklijke Philips N.V. System and method for mapping ultrasound shear wave elastography measurements
KR102406937B1 (ko) * 2014-11-07 2022-06-10 삼성메디슨 주식회사 초음파 영상장치 및 그 제어 방법
US10842445B2 (en) * 2018-11-08 2020-11-24 General Electric Company System and method for unsupervised deep learning for deformable image registration

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0659387A2 (en) * 1993-12-24 1995-06-28 Olympus Optical Co., Ltd. Ultrasonic diagnosis and therapy system in which focusing point of therapeutic ultrasonic wave is locked at predetermined position within observation ultrasonic scanning range
CN1342057A (zh) * 1998-10-26 2002-03-27 株式会社日立制作所 超声波治疗装置
US20070016031A1 (en) * 2000-11-28 2007-01-18 Allez Physionix Limited Systems and methods for making noninvasive assessments of cardiac tissue and parameters
US20070167705A1 (en) * 2005-08-04 2007-07-19 Chiang Alice M Integrated ultrasound imaging system
CN101053531A (zh) * 2007-05-17 2007-10-17 上海交通大学 基于多模式增敏成像融合的早期肿瘤定位跟踪方法
US20080039723A1 (en) * 2006-05-18 2008-02-14 Suri Jasjit S System and method for 3-d biopsy
CN101152646A (zh) * 2006-09-27 2008-04-02 香港理工大学 柔性超声换能器阵列及其应用装置
WO2009023432A2 (en) * 2007-08-09 2009-02-19 Juvent Inc. Method and apparatus for ultrasound therapy
US20110044521A1 (en) * 2009-06-26 2011-02-24 Regents Of The University Of Minnesota Tracking and representation of multi-dimensional organs

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5769790A (en) 1996-10-25 1998-06-23 General Electric Company Focused ultrasound surgery system guided by ultrasound imaging
US7257244B2 (en) * 2003-02-24 2007-08-14 Vanderbilt University Elastography imaging modalities for characterizing properties of tissue
US8233681B2 (en) * 2004-09-24 2012-07-31 The University Of North Carolina At Chapel Hill Methods, systems, and computer program products for hierarchical registration between a blood vessel and tissue surface model for a subject and a blood vessel and tissue surface image for the subject
DE102006058264A1 (de) 2006-12-08 2008-06-12 Rheinische Friedrich-Wilhelms-Universität Bonn MR-Bildkontrast mittels Ultraschall
US20080242979A1 (en) * 2007-03-30 2008-10-02 Rayette Ann Fisher Combined X-ray detector and ultrasound imager
US10028722B2 (en) * 2007-09-25 2018-07-24 Hospital For Special Surgery Methods and apparatus for assisting cartilage diagnostic and therapeutic procedures
US8187187B2 (en) 2008-07-16 2012-05-29 Siemens Medical Solutions Usa, Inc. Shear wave imaging
DE102008047644B4 (de) 2008-09-17 2015-09-10 Siemens Aktiengesellschaft Verfahren zur Registrierung zweier Bildgebungsmodalitäten
US9652696B2 (en) * 2010-07-09 2017-05-16 Temple University-Of The Commonwealth System Of Higher Education Apparatus and method for surface and subsurface tactile sensation imaging
JP5373749B2 (ja) * 2010-11-30 2013-12-18 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 医用画像表示装置及びその制御プログラム
JP5846411B2 (ja) * 2011-06-16 2016-01-20 学校法人上智学院 イメージング方法及び変位計測方法及び装置、並びに、超音波画像診断装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0659387A2 (en) * 1993-12-24 1995-06-28 Olympus Optical Co., Ltd. Ultrasonic diagnosis and therapy system in which focusing point of therapeutic ultrasonic wave is locked at predetermined position within observation ultrasonic scanning range
CN1342057A (zh) * 1998-10-26 2002-03-27 株式会社日立制作所 超声波治疗装置
US20070016031A1 (en) * 2000-11-28 2007-01-18 Allez Physionix Limited Systems and methods for making noninvasive assessments of cardiac tissue and parameters
US20070167705A1 (en) * 2005-08-04 2007-07-19 Chiang Alice M Integrated ultrasound imaging system
US20080039723A1 (en) * 2006-05-18 2008-02-14 Suri Jasjit S System and method for 3-d biopsy
CN101152646A (zh) * 2006-09-27 2008-04-02 香港理工大学 柔性超声换能器阵列及其应用装置
CN101053531A (zh) * 2007-05-17 2007-10-17 上海交通大学 基于多模式增敏成像融合的早期肿瘤定位跟踪方法
WO2009023432A2 (en) * 2007-08-09 2009-02-19 Juvent Inc. Method and apparatus for ultrasound therapy
US20110044521A1 (en) * 2009-06-26 2011-02-24 Regents Of The University Of Minnesota Tracking and representation of multi-dimensional organs

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105338905A (zh) * 2013-06-26 2016-02-17 皇家飞利浦有限公司 用于多模态组织分类的方法和***
US11660065B2 (en) 2013-06-26 2023-05-30 Koninklijke Philips N.V. Method and system for multi-modal tissue classification
CN111493931A (zh) * 2019-08-01 2020-08-07 深圳迈瑞生物医疗电子股份有限公司 一种超声成像方法及设备、计算机可读存储介质

Also Published As

Publication number Publication date
DE102011085894B4 (de) 2014-03-27
DE102011085894A1 (de) 2013-05-08
US20130116542A1 (en) 2013-05-09
US8897855B2 (en) 2014-11-25

Similar Documents

Publication Publication Date Title
CN110248606B (zh) 气穴定位
CN102641137B (zh) 使用幅度-相位调制超声波的粘弹性测量
JP6430498B2 (ja) 超音波剪断波エラストグラフィ測定のマッピングのためのシステムおよび方法
AU2012326218B2 (en) Transmissive imaging and related apparatus and methods
CN103908256B (zh) 用于支持多参数图像的获取的设备和方法
CN103083046A (zh) 用于成像的方法和电子计算装置以及成像***
CN108135565A (zh) 用于配准使用各种成像模态获得的图像并验证图像配准的***和方法
CN102686278B (zh) 治疗设备
CN103356189A (zh) 磁共振和超声参数化图像融合
CN101902970B (zh) 弹性成像中的应变增益补偿的方法和***
Solovchuk et al. Temperature elevation by HIFU in ex vivo porcine muscle: MRI measurement and simulation study
CA2902550A1 (en) Transmissive imaging and related apparatus and methods
KR20090119860A (ko) 수차 삽입 부재를 통해 파동의 초점 집중을 최적화하는 방법
CN103826534A (zh) 避免共存的***的mri干扰
CN106999162A (zh) 超声波诊断装置以及弹性评价方法
CN103239258A (zh) 采用超声波的同轴切变波表征
US9508154B2 (en) Medical imaging apparatus and method of operating the same
Ruiter et al. USCT data challenge
Mougenot et al. High intensity focused ultrasound with large aperture transducers: a MRI based focal point correction for tissue heterogeneity
US20160133003A1 (en) Method and apparatus for image registration
CN107106128A (zh) 用于分割解剖目标的超声成像装置和方法
KR101027599B1 (ko) 고 프레임율 탄성 초음파 영상을 제공하는 초음파 진단 시스템 및 방법
CN106102576A (zh) 动态医疗图像中运动的检测
CN110418609A (zh) 一种超声弹性测量装置及弹性对比测量方法
Yu et al. Feasibility study of 3D time‐reversal reconstruction of proton‐induced acoustic signals for dose verification in the head and the liver: A simulation study

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130508

WD01 Invention patent application deemed withdrawn after publication