CN100595606C - 用于探测目标运动的超声波方法和设备 - Google Patents
用于探测目标运动的超声波方法和设备 Download PDFInfo
- Publication number
- CN100595606C CN100595606C CN200480044281.5A CN200480044281A CN100595606C CN 100595606 C CN100595606 C CN 100595606C CN 200480044281 A CN200480044281 A CN 200480044281A CN 100595606 C CN100595606 C CN 100595606C
- Authority
- CN
- China
- Prior art keywords
- volume
- position information
- layer position
- time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 19
- 230000000747 cardiac effect Effects 0.000 claims abstract description 17
- 210000004027 cell Anatomy 0.000 claims description 23
- 239000013598 vector Substances 0.000 claims description 21
- 238000005516 engineering process Methods 0.000 claims description 12
- 230000005484 gravity Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 210000004165 myocardium Anatomy 0.000 claims description 5
- 241001465754 Metazoa Species 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims description 4
- 230000000241 respiratory effect Effects 0.000 claims description 4
- 238000002847 impedance measurement Methods 0.000 claims description 3
- 210000001519 tissue Anatomy 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 238000000844 transformation Methods 0.000 claims description 2
- 210000004204 blood vessel Anatomy 0.000 abstract description 2
- 238000009877 rendering Methods 0.000 abstract description 2
- 238000000638 solvent extraction Methods 0.000 abstract 1
- 239000010410 layer Substances 0.000 description 47
- 230000008602 contraction Effects 0.000 description 27
- 239000012634 fragment Substances 0.000 description 13
- 238000003384 imaging method Methods 0.000 description 5
- 239000000523 sample Substances 0.000 description 5
- 210000005240 left ventricle Anatomy 0.000 description 4
- 238000002604 ultrasonography Methods 0.000 description 4
- 230000004069 differentiation Effects 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 210000002837 heart atrium Anatomy 0.000 description 2
- 230000010247 heart contraction Effects 0.000 description 2
- 230000035479 physiological effects, processes and functions Effects 0.000 description 2
- 231100000241 scar Toxicity 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 238000007514 turning Methods 0.000 description 2
- 206010019280 Heart failures Diseases 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 238000013155 cardiography Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000003238 esophagus Anatomy 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000001308 heart ventricle Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000004118 muscle contraction Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 230000000661 pacemaking effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 210000005241 right ventricle Anatomy 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8993—Three dimensional imaging systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
- G01S7/52087—Details related to the ultrasound signal acquisition, e.g. scan sequences using synchronization techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
-
- 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
-
- 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/504—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 blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
- A61B8/543—Control of the diagnostic device involving acquisition triggered by a physiological signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
- G06T2219/20—Indexing scheme for editing of 3D models
- G06T2219/2012—Colour editing, changing, or manipulating; Use of colour codes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Graphics (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Hardware Design (AREA)
- Architecture (AREA)
- Physiology (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及用于检测目标运动的方法和装置,优选用于在心动周期中检测血管运动的超声方法和装置,包括步骤:使用表面重构和/或体积描绘技术,通过二维或三维超声图像扫描装置(4)扫描目标(1),以离散采集时间(t1-t6)采集数据集,将三维数据集分成多个体积单元(21),其中各个体积单元(21)具有单元表面(20),计算各个体积单元(21)在所述采集时间(t1-t6)上的体积变化,以获得在时间上各个相应单元表面(20)的层位信息(13-15),以及以二维或三维方式显示所述三维数据集及其相对时间的层位信息(13-15)。
Description
发明领域
本发明涉及用于探测目标运动的超声波方法和设备,优选用于心动周期中肌肉组织的运动,更优选涉及用于探测诸如心室和心房的收缩图之类的肌肉组织的一部分的运动的超声方法和设备。更具体地,本发明涉及一种用于视频显现人或动物体的心脏的收缩图的方法和设备,进而呈现例如心室的表面的图形化表示,以探测它是否具有同步或异步收缩。
背景技术
目前,用于超声三维成像的方法和设备是已知的,其中通过二维或三维超声成像扫描装置扫描目标以获得目标的横截面图像,其中各个横截面图像的离散采集时间连同图像信息一起被记录。此后,各个所述横截面图像被数字化且连同其相应位置和采集时间一起被记录。还有,本领域众所周知的是,使用表面重构和/或体积描绘(volume rendering)技术将所述横截面图像转变为二维或三维数据集,进而显示随时间的三维数据集,从而以三维方式显示被扫描目标的图形化表示,该图形化表示通过显示这些离散采集时间点的三维表示可以是活动的。
EP 0961 135A1公开了一种快速生成被扫描目标的线帧(wire-frame)体积模型以及如何通过在线帧体积模型上插值而显示目标表面的方法和设备。从而,使三维视频显示例如病人心脏及其在虚拟真实表面包层中随时间的运动成为可能。
然而,因为在心动周期中发生的相当复杂的心脏运动和心脏结构的动态改变,心脏的超声扫描(超声波心动描记术)遇到特殊的挑战和问题。因这些问题或其它问题,心动超声扫描在历史上已经被限制至二维或三维成像,其中不可能以一种可接受的方式在整个心动周期中观察或显现心肌收缩。对于心脏病专家和电生理专家,极其重要的是理解心脏肌肉尤其是心室和心房的收缩图。有规律的同步收缩产生心动功能的理想效果。心室间或心室内传导扰动可以导致心室的异步收缩,所以降低心脏的功能和效率。评估异步收缩的目前方法是有限的且劳动强度大的。至今,通过其它成像形态例如CT或MRI(磁共振成像)来检测心室内表面的收缩传播。然而,因对心脏心室壁的运动的采集和分析均是耗时且复杂的,因此这些技术并不普及。
目前,MRI和CT技术能够提供静态参数分布,从而通过使用彩色编码(功能成像)二维显现心室的收缩图。缺点在于可以由人眼区分的有限的颜色数量。此外,关于局部差异,对颜色图的解释通常不是直观的。例如,没有对于“黄色”大于或小于“绿色”的当然解释。工作区(workaround)是具有有限数量颜色(例如,兰和红)和平滑过渡的颜色图的使用。然而,尽管解释变得更直观,但是毗邻色度间的区分变得更困难。因此,在某一时刻,仅仅对有限的数值间隔进行描绘,而不具有对心脏表面的收缩传播的整体观察。
发明目标
因此,本发明的一个目的是公开一种以快速可靠方式探测目标运动的超声方法和设备,尤其是描述一种允许快速和可再现地评价诸如心肌的血管的收缩的方法和设备,尤其是视频显示人心脏表面上的收缩波随时间的传播,尤其在心室的内表面。
同时,本发明的目的是公开一种方法和设备,其非侵入和快速地获得上述信息以减少检测时间,从而尽可能地准确。
发明内容
本发明表述如下:
所要求的用于探测目标运动的超声方法包括使用图像扫描装置扫描目标的步骤,优选二维或三维超声图像扫描装置,以离散采集时间采集目标的至少二维横截面图像,数字化各个横截面图像,记录各个这些横截面图像及其相应的位置和采集时间,使用表面重构和/或体积描绘技术将这些图像转变为三维数据集,以及相对时间显示三维数据集,其中三维数据集被分成多个体积单元,各个体积单元具有单元表面,以及其中计算各个体积单元在所述采集时间上的体积改变以获得各个相应单元表面在相对间的层位(level)信息,以及然后显示所述层位信息。这种过程还应用于三维数据集被直接生成的采集方法。
优选地,二维超声图像扫描装置被用于采集目标的二维横截面图像。然而,还可能的是使用三维超声图像扫描装置采集作为所述目标的体积的三维横截面图像。因此,术语“横截面图像”涵盖两种可能性。尽管本发明使用已知技术获得三维数据集以显示被检测的目标的轮廓,例如心动周期中的心脏,现在令人意外的是,使用超声设备显示心动周期中心室内表面的运动是可能的。为了获得心脏表面的子集或片段的心脏收缩的同步性的个别信息,本发明将所述三维数据集分成多个体积单元并计算体积变化,以获得目标的分散的小表面单元的层位信息。通过这种方法,使人心室的收缩正面以三维或四维方式显现是可行的。
优选实施例的描述
优选地,各个单元表面的层位信息为一系列颜色或一系列灰度值,以容易显现各个体积单元的表面所处的“层位或位置”。优选通过使用被检测目标内的重心及连接所述单元表面的所有角和重心,将三维数据集分成多个体积单元。因此,整个目标被分成所述多个体积单元,各个单元为被检测目标的三维片段(segment)。
同时,有益的是,使用离散层位边界划分层位信息,以及将优选两个或三个的有限数量的颜色分配给所述离散层位信息。因此,可能的是以简单方式显现各个体积单元的表面的两个或三个位置,从而指示表面是否是处于“收缩”状态或“扩张”状态。通过二维或三维显示对时间的所述层位信息,可能显现目标表面的运动,如同在心动周期中收缩波在心脏表面上的传播一样。
更优选地,通过使用本身已知的二维变换方法,单元表面连同相应层位信息被二维显示在二维图上。
通过二维地显现三维目标(象地球的二维世界图一样),可能的是,通过片段的二维表面的颜色表示显现目标表面的各个片段的运动。由于颜色间的边界将跨越二维图像,因此可能的是,遵循例如心室收缩的传播波。
依据本发明的优选实施例,所述层位信息的传播随时间显示在所述二维图上,其中还可以显示另外的二维传播矢量。使用所述二维传播矢量的长度和方向,可能的是在所述图上显现各体积单元的体积变换的方向、数量和速度以及它们各自的主要位置和它们的数量。如果多个毗邻体积单元同时移入同一方向,但比目标中其它位置的体积单元早或晚点,则传播矢量将具有从第一组体积单元至第二组体积单元的方向,其中传播矢量的长度为波传播速度和/或体积变化差的量值。这里,各个体积单元的体积变化优选通过应用降噪算法、边缘增强算法和/或空间伪影减小算法来计算。
优选依据考虑通过心电图门控的得到的心脏周期变化和/或通过阻抗测量得到的呼吸周期变化的算法,选择用于心脏目标的所述离散采集时间。这些在本技术领域中是公知的。
依据本发明的设备包括用于扫描目标以在离散采集时间处采集所述目标横截面图像的二维或三维超声扫描装置以及用于数字化各个所述横截面图像的数字化装置、记录各个所述横截面图像及其相应位置和采集时间的记录装置、使用表面重构和/或体积描绘技术以将所述横截面图像转变为三维数据集的转换装置以及相对时间显示所述三维数据集的显示装置,还有划分装置、计算装置,该划分装置将三维数据集分成多个体积单元,各个体积单元具有单元表面,该计算装置计算各个体积单元在所述采集时间上的体积变化以获得各个相应单元表面对时间的层位信息,其中显示装置显示所述层位信息。
计算装置优选能够定义一个重心,其是所有体积单元的一个角,其它角全部被包括在目标的表面中。因此,保证各个体积单元的体积变化导致各个单元表面的层位变化。所述划分装置划分所述层位信息,优选使用离散层位边界和将优选两个或三个的有限数量的颜色分配给所述层位信息。同时使用已知的二维转换方法,转换装置将所述单元表面及其层位信息转换为对时间的二维图,以及所述显示装置优选二维地显示所述层位信息。
该设备另外包括心电图门控装置和呼吸触发装置,用于动态扫描心脏目标。同时,所述设备具有能够显示本发明的二维图和心电图影像***的二维图的显示装置,从而能够比较心脏表面中各个片段的电势值的传播和收缩值的传播,以对医生提供用于心脏诊断的可靠信息。
结合附图,将解释本发明的一些优选实施例,附图如下:
图1a显示内诊镜超声波心动描记法采集,
图1b显示心脏的连续横截面图像采集的原理,
图2显示三维换能器探头的例子,
图3示意性显示三维目标的片段表面,
图3a显示被分为多个体积单元的图3中的一个片段,
图3b显示图3a中的一个体积单元,
图4显示图3b中的体积单元随时间的体积变化,
图5显示图3中目标的一个局部片段,具有在6个不同时间点处具有相等收缩水平的线,
图6显示图3/图5中目标的二维表示,以及
图7显示在6个不同采集时间处的图3的目标的三维和二维表示。
图1a显示类似内部器官或组织的目标1(例如心脏),其由来自如换能器之类的扫描装置4的超声束8扫描,该换能器被结合在位于临近目标1的静脉或动脉(例如临近心室的食道)中的探头内,其中扫描装置4通过内窥镜通道2内的内窥镜连接部3与处理装置(未显示)相连接。跨食道的(transesophageal)探头4在标准内窥镜远端容纳可旋转阵列的换能器。扫描平面可以连续旋转180°,从纵向成像位置开始经由内窥镜手柄上的控制把手。心脏横截面包围锥状体积,其锥点始于换能器(同样比较图1b)。然后,通过控制逻辑激励步进电动机,控制逻辑使用考虑通过ECG(心电图)门控5得到的心脏周期变化和通过阻抗测量得到的呼吸周期变化的算法控制在给定平面中的图像采集。对于心脏的各个位置(例如,心脏收缩或心脏舒张),通过获得一系列横截面6来扫描心脏,从而获得横截面图像集7,该横截面图像集7属于心脏的特定位置(其类似于表示在各个采集时间t1、t2、……t6的心脏特定位置的频闪观测图片)。
图1b显示对连续横截面6的采集原理的放大视图。通过内窥镜连接3被链接至计算机的换能器4扫描由左心室11和右心室10构成的心脏1,通过超声束8采集所有采集平面9,从而以方向S扫描超声束。
图2显示三维超声图像扫描装置的例子,例如能够采集3-D(三维)图像元素的矩阵23的超声探头。因此,可以采集三维横截面图像,也就是体积。其它图像采集技术也是可行的,例如超声导管探头。
图3示意显示目标1表面的三维表示,该目标1被分为带有片段边界12的片段11。三维体积16可被进一步分成多个体积单元21,各个体积单元具有单元表面20。优选通过连接表面上各个拐角(图3a中,各个体积单元具有表面上的四个拐角)与被定位在三维体积16内的重心18,实施这种划分。如图3b中所示,由于重心18处于固定位置,各个体积单元21具有通过各种层位扩展或收缩的单元表面20。只要改变三维体积16的体积,各个体积单元21的体积改变ΔV将导致单元表面20的层位改变。
如图3b中所示,通过可预先确定的层位边界10划分三个层位13、14和15。只要时间变化ΔT足以将单元表面升高至层位边界10上方,则依据本发明的优选实施例,该表面将接收例如另外颜色的另外离散层位信息。因此,可能的是,例如将兰颜色分配给第一层位13,将灰颜色分配给第二层位15,以及将红颜色分配给第三层位14。
只要单元表面20处于三种上述层位13、14或15之一内,则表面将分配特定层位信息例如,分配兰颜色为第一层位13。然后,如图3a中所示的体积单元21的各个单元表面20将分配兰色。
图4显示三种体积单元20a、20b和20c(图3b中的20b)随时间的体积变化ΔV。如果检测心脏的心动周期,则t1为在心脏舒张结束后不久的时间,tp为心脏收缩结束处的时间,以及t6为舒张末期之前不久的时间。在tp时,体积单元20b达到其最低体积,以使单元表面20b达到其最低层位。各个时间tp不是所有体积单元20a、20b或20c达到其最小体积时的时间,而仅仅是一个单个体积单元20b达到其最小体积的时间。在t1和t6之间的其它最低体积时间tp2、tp3处,其它体积单元21达到其最低“位置”。
如可以从图4中看出,所有体积单元20a、20b和20c接收不同层位信息,也就是第一体积单元20a接收信息14,因为最低体积的时间tp2处于区域14中,(图3b中的)第二体积单元20b接收信息15,因为最低体积的时间tp1处于区域15中,以及第三体积单元20c接收信息13,因为最低体积的时间tp3处于区域13中。在随时间“移动”层位信息时,各个体积单元将随时间接收不同层位信息,例如颜色,因此,可能的是使收缩波象动画一样显示。静态地,各个体积单元具有特定“时间印记”ΔT,也就是达到最低体积时间tp所需要的时间。
图5显示在三维空间中二维表面的动态显示。带有层位边界10的体积16被部分显示在图5中,层位边界10显示各种层位13、14和15的边界。第一层位13指示已经达到其峰值收缩的壁片段。从t1至t2,颜色界限变换,以使用于层位13的颜色被赋予(address)层位14,以及对于t3,被赋予层15等。这样,特定颜色指示瞬时处于其峰值收缩阶段的区域,以及随后,颜色在时间上的传播显示心室的机械肌肉收缩面。心室的收缩面随时间传播,由于心室内表面的各个片段20在不同时间在不同方向中移动。
在将如图5中所示的体积16的三维表示转变为如图6中所示的二维图时,一种二维转变方法被使用。优选地,三维体积的中心轴22被用作二维图的中心,所述二维图展开三维体积的表面,如由片段11和片段边界12示意所示。
因此,图6显示图3或图5的三维体积的这些片段11,单元表面20的第一子集(仅仅局部被显示在图6中)指示第一层位13,单元表面20的第二子集指示第三层位14,以及单元表面20的其余部分指示第三子集,也就是第二层位15,例如为灰色。图6中的体积表示显示离散采集时间t1,……,t6,目标的特定范围/区域处于第三层位14,以及另外的范围或区域处于第一层位13(例如,兰色),从而指示第一区域具有处于第一层位13的状态,同时第二区域具有处于第三层位14的另外状态(例如红色)。
图6还显示源自中心轴22的传播矢量19,并且指向收缩传播方向。由于第一区域处于第一层位13中,传播矢量19指向该区域的收缩源的方向。
图7显示六个不同采集时间t1、……、t6,其中在三维体积16的表面的二维表示上显示依据图3中的三维体积16。第一层位13、第二层位15和第三层位14通过二维区域在二维图上以及因此在三维表面16上显示。
图7显示人心室的左心室的圆形二维收缩面图,从而显示处于第一层位13的区域(收缩)和处于第三层位14的区域(扩张)。通过选择时间t1……t6,颜色的动态范围在上、下界限处被混叠。因此,二维显示仅仅显示两种颜色,也就是由第二层位灰或白所区分的兰色和红色。通过心动循环使时间窗口进一步向前,可以观察到以前进方式达到峰值收缩tp的区域。然后,这可被显示为数字动画,将机械峰值收缩显示为在整个心脏左心室中的连续波。心脏表面的各个静态位置的彩色编码实现所有单元表面20的彩色影像和心脏收缩波传播的连续影像。因此,可以根据需要尽可能精细地调节单帧内的分辨率,以通过增加采集时间的数量及因此的帧数、以及通过增加体积单元的数量来解决局部差异。
图7的例子显示6个静态帧的序列。然而,本发明能够将这些帧显示为动态剪辑的部分,从而被作为动画观看。特定颜色随时间在二维图上传播显示收缩图案的暂时演变,从而允许对心脏周期的同步或异步及其传播容易地评价。这种显示可以容易地与电生理电势随时间变化的影像相比较,从而比较心脏肌肉的激励和收缩。如图6中所示,传播矢量19还可被用于这种动画表示。除了暂时评价之外可以考虑收缩空间图案,以产生这种传播矢量。二维矢量表示已经达到峰值收缩的左心室肌层数量和收缩波的平均方向。这可以被实现为叠加在收缩面传播上的移动箭头。矢量原点优选基于参考图像的中心22中。
可以增加另外的矢量,用于以同样方式设计的辅助收缩源。例如,这是必需的,如疤痕区域的明显二级失调(电信号在它们围绕疤痕时将延迟,以及可能导致几乎同时开始收缩的两个分离区域)产生的心力衰竭之后或在二室起搏(biventricular pacing)中所显现的,一些病人可能显示两种或多种收缩波。
传播矢量19的可替换的矢量表示还可通过以下方式实现:定义其原点处于最早峰值收缩区域的矢量,即该区域首先达到第一层位13,矢量的长度和方向跟随通过左心室的峰值收缩波。在多个波以各种方向从原点区域传播的情况中,可以以同样方式显示多传播矢量19(多收缩矢量影像-MCVM)。
上述方法还可被用于使用CT或MR技术所采集的图像,只要数据能够使用表面重构和/或体积描绘技术被转换为三维数据集,以及能够被显示为时间上的三维数据集。
Claims (15)
1.一种用于检测人或动物体的心脏的移动的方法,包括步骤:
-通过图像扫描装置(4)扫描心脏(1),以在离散采集时间(t1-t6)处采集所述心脏(1)的至少二维横截面图像(7),
-数字化各个所述横截面图像(7),
-记录各个所述横截面图像(7)及其相应位置和采集时间(t1-t6),
-使用表面采集和/或体积描绘技术,将所述横截面图像(7)转变为三维数据集,以及
-相对时间显示所述三维数据集,
特征在于,
-将三维数据集分成多个体积单元(21),其中各个体积单元(21)具有单元表面(20),
-计算各个体积单元(21)在所述采集时间(t1-t6)上的体积变化,以获得各个相应单元表面(20)对时间的层位信息(13-15),以及
-显示所述层位信息(13-15)。
2.如权利要求1所要求的方法,特征在于,
所述层位信息(13-15)为一系列颜色或一系列灰度值。
3.如权利要求1或2所要求的方法,特征在于,
通过使用所述心脏(1)内的重心(18),并将所述单元表面(20)的所有角与所述重心(18)连接,将所述三维数据集分成多个体积单元(21)。
4.如权利要求1或2所要求的方法,特征在于,
使用离散层位边界(10)将所述层位信息(13-15)划分成离散层位信息以及将有限数量的颜色分配给所述离散层位信息。
5.如权利要求1或2所要求的方法,特征在于,
相对时间三维地显示所述层位信息(13-15)。
6.如权利要求1或2所要求的方法,特征在于,
通过使用二维变换方法,二维地在图上显示所述单元表面(20)及其相应层位信息(13-15)。
7.如权利要求6所要求的方法,特征在于,
在所述图上显示所述层位信息(13-15)相对时间的传播。
8.如权利要求6所要求的方法,特征在于,
附加地显示二维传播矢量。
9.如权利要求8所要求的方法,特征在于,
相关于各体积单元(21)的体积改变的方向、数量及速度及其在所述图上的位置和数量,显示所述二维传播矢量(19)的长度和方向。
10.如权利要求1或2所要求的方法,特征在于,
通过降噪算法、边缘增强算法和/或空间伪影减小算法的应用,计算各个体积单元(21)的体积改变。
11.如权利要求1或2所要求的方法,特征在于,
依据考虑通过心电图门控获得的心脏周期变化和通过阻抗测量获得的呼吸周期变化的算法,选择用于运动中的心脏肌肉组织(1)的所述离散采集时间(t1-t6)。
12.用于检测人或动物体的心脏移动的设备,包括用于扫描人或动物体的心脏(1)以在离散采集时间(t1-t6)处采集所述心脏(1)的至少二维横截面图像(7)的超声扫描装置(4)、用于数字化各个所述横截面图像(7)的数字化装置、记录各个所述横截面图像(7)及其相应位置和采集时间(t1-t6)的记录装置、使用表面重构和/或体积描绘技术以将所述横截面图像(7)转变为三维数据集的转换装置以及相对时间显示所述三维数据集的显示装置,
特征在于,
划分装置,其将三维数据集分成多个体积单元(21),其中各个体积单元(21)具有单元表面(20),
计算装置,其计算各个体积单元(21)在所述采集时间(t1-t6)上的体积变化,以获得各个相应单元表面(20)相对时间的层位信息(13-15),以及
所述显示装置显示所述层位信息(13-15)。
13.如权利要求12所要求的设备,特征在于,
所述划分装置使用离散层位边界(10)将所述层位信息(13-15)划分成离散层位信息以及将有限数量的颜色分配给所述离散层位信息,
转换装置通过使用二维转换方法在图上二维地随时间转换所述单元表面(20)及其各个层位信息(13-15),以及
所述显示装置二维地显示所述层位信息(13-15)。
14.如权利要求13所要求的设备,特征在于,
所述显示装置显示二维传播矢量(19),从而相关于各体积单元(21)的体积改变的方向、数量及速度及其在所述图上的位置和数量显示所述二维传播矢量(19)的长度和方向。
15.如权利要求12或13所要求的设备,特征在于,
所述设备进一步包括心电图门控装置和呼吸触发装置,用于动态扫描心脏肌肉组织(1)。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/EP2004/013387 WO2006056221A1 (en) | 2004-11-25 | 2004-11-25 | Ultrasonic method and apparatus for detecting movements of objects |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101052893A CN101052893A (zh) | 2007-10-10 |
CN100595606C true CN100595606C (zh) | 2010-03-24 |
Family
ID=34959454
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200480044281.5A Expired - Fee Related CN100595606C (zh) | 2004-11-25 | 2004-11-25 | 用于探测目标运动的超声波方法和设备 |
Country Status (7)
Country | Link |
---|---|
US (1) | US8096948B2 (zh) |
EP (1) | EP1680690B1 (zh) |
CN (1) | CN100595606C (zh) |
AT (1) | ATE481652T1 (zh) |
DE (1) | DE602004029174D1 (zh) |
ES (1) | ES2350045T3 (zh) |
WO (1) | WO2006056221A1 (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1933703A2 (en) | 2005-09-16 | 2008-06-25 | The Ohio State University | Method and apparatus for detecting interventricular dyssynchrony |
WO2007035688A2 (en) | 2005-09-16 | 2007-03-29 | The Ohio State University | Method and apparatus for detecting intraventricular dyssynchrony |
WO2008073491A1 (en) * | 2006-12-11 | 2008-06-19 | University Of Florida Research Foundation, Inc. | System and method for analyzing progress of labor and preterm labor |
JP5238201B2 (ja) * | 2007-08-10 | 2013-07-17 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
JP5454844B2 (ja) * | 2008-08-13 | 2014-03-26 | 株式会社東芝 | 超音波診断装置、超音波画像表示装置及び超音波画像表示プログラム |
US8745536B1 (en) | 2008-11-25 | 2014-06-03 | Perceptive Pixel Inc. | Volumetric data exploration using multi-point input controls |
DE102011081987B4 (de) * | 2011-09-01 | 2014-05-28 | Tomtec Imaging Systems Gmbh | Verfahren zur Erzeugung eines Modells einer Oberfläche einer Hohlraumwand |
KR101386102B1 (ko) * | 2012-03-09 | 2014-04-16 | 삼성메디슨 주식회사 | 초음파 영상 제공 방법 및 그를 위한 초음파 장치 |
CN103077550B (zh) * | 2012-12-10 | 2016-04-20 | 华北电力大学(保定) | 一种非门控icus图像序列中血管的四维重建方法 |
EP2939217B1 (en) * | 2012-12-27 | 2022-06-22 | Koninklijke Philips N.V. | Computer-aided identification of a tissue of interest |
JP6535097B2 (ja) * | 2015-09-29 | 2019-06-26 | 古野電気株式会社 | 超音波組織検出装置、超音波組織検出方法、および、超音波組織検出プログラム |
IT201800009482A1 (it) | 2018-10-16 | 2020-04-16 | Italia Tech Alliance Srl | Macchina di avvolgimento con dispositivi per il calcolo del coefficiente di poisson e metodo |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5622174A (en) * | 1992-10-02 | 1997-04-22 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus and image displaying system |
NO943696D0 (no) | 1994-10-04 | 1994-10-04 | Vingmed Sound As | Fremgangsmåte ved ultralydavbildning |
JPH11221210A (ja) * | 1998-02-09 | 1999-08-17 | Aloka Co Ltd | 超音波診断装置 |
ATE228252T1 (de) | 1998-03-30 | 2002-12-15 | Tomtec Imaging Syst Gmbh | Verfahren und vorrichtung zur bild-aufnahme mit ultraschall |
JP3872424B2 (ja) | 2002-12-20 | 2007-01-24 | アロカ株式会社 | 超音波診断装置 |
US8064983B2 (en) * | 2003-11-26 | 2011-11-22 | General Electric Company | Method and system for prospective gating using multiple inputs |
-
2004
- 2004-11-25 AT AT04821950T patent/ATE481652T1/de active
- 2004-11-25 ES ES04821950T patent/ES2350045T3/es active Active
- 2004-11-25 WO PCT/EP2004/013387 patent/WO2006056221A1/en active Application Filing
- 2004-11-25 US US10/558,602 patent/US8096948B2/en not_active Expired - Fee Related
- 2004-11-25 EP EP04821950A patent/EP1680690B1/en not_active Not-in-force
- 2004-11-25 CN CN200480044281.5A patent/CN100595606C/zh not_active Expired - Fee Related
- 2004-11-25 DE DE602004029174T patent/DE602004029174D1/de active Active
Also Published As
Publication number | Publication date |
---|---|
ES2350045T3 (es) | 2011-01-17 |
DE602004029174D1 (de) | 2010-10-28 |
WO2006056221A1 (en) | 2006-06-01 |
EP1680690A1 (en) | 2006-07-19 |
US20080132788A1 (en) | 2008-06-05 |
ATE481652T1 (de) | 2010-10-15 |
EP1680690B1 (en) | 2010-09-15 |
US8096948B2 (en) | 2012-01-17 |
CN101052893A (zh) | 2007-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3187148B2 (ja) | 超音波診断装置 | |
CN100525714C (zh) | 图像处理设备和超声诊断设备 | |
US11317896B2 (en) | Ultrasound diagnosis apparatus and image processing apparatus | |
JP4734395B2 (ja) | 超音波診断装置 | |
US5916168A (en) | Three dimensional M-mode ultrasonic diagnostic imaging system | |
CN106456128B (zh) | 医学图像处理设备和方法 | |
JP4223775B2 (ja) | 超音波診断装置 | |
EP0920639B1 (en) | System and method for the dynamic display of three-dimensional image data | |
CN100595606C (zh) | 用于探测目标运动的超声波方法和设备 | |
US8540636B2 (en) | Ultrasonic diagnostic apparatus and medical image processing apparatus | |
Ofili et al. | Three-dimensional and four-dimensional echocardiography | |
JP2004195082A (ja) | 超音波診断装置 | |
Roelandt et al. | Precordial three‐dimensional echocardiography with a rotational imaging probe: Methods and initial clinical experience | |
CN106030657B (zh) | 医学4d成像中的运动自适应可视化 | |
Pini et al. | Echocardiographic three-dimensional visualization of the heart | |
JP2001128982A (ja) | 超音波画像診断装置および画像処理装置 | |
JP7019594B2 (ja) | 撮像システムおよび方法 | |
Bruining et al. | Three-dimensional echocardiography paves the way toward virtual reality | |
Tamura et al. | Three-dimensional reconstruction of echocardiograms based on orthogonal sections | |
Brekke et al. | Volume stitching in three-dimensional echocardiography: distortion analysis and extension to real time | |
Nadkarni et al. | In vitro simulation and quantification of temporal jitter artifacts in ECG-gated dynamic three-dimensional echocardiography | |
JP2004195028A (ja) | 超音波診断装置 | |
CN116263948A (zh) | 用于图像融合的***和方法 | |
Chen et al. | Three-dimensional left ventricular wall motion analysis using multiplane transesophageal echocardiography | |
Spicer et al. | Three-dimensional echocardiography: research toy or clinical tool? |
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: 20100324 Termination date: 20211125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |