CN104968279B - 消除医学图像中由生理功能引起的运动影响 - Google Patents
消除医学图像中由生理功能引起的运动影响 Download PDFInfo
- Publication number
- CN104968279B CN104968279B CN201480005221.6A CN201480005221A CN104968279B CN 104968279 B CN104968279 B CN 104968279B CN 201480005221 A CN201480005221 A CN 201480005221A CN 104968279 B CN104968279 B CN 104968279B
- Authority
- CN
- China
- Prior art keywords
- image
- image data
- movement
- period
- sequence
- 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.)
- Active
Links
- 230000035790 physiological processes and functions Effects 0.000 title description 4
- 230000033001 locomotion Effects 0.000 claims abstract description 60
- 238000000034 method Methods 0.000 claims description 43
- 239000002872 contrast media Substances 0.000 claims description 39
- 230000003902 lesion Effects 0.000 claims description 30
- 230000000241 respiratory effect Effects 0.000 claims description 30
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 15
- 230000010412 perfusion Effects 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 8
- 210000002216 heart Anatomy 0.000 claims description 6
- 230000010410 reperfusion Effects 0.000 claims description 6
- 206010028980 Neoplasm Diseases 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 2
- 239000003814 drug Substances 0.000 claims 1
- 230000000116 mitigating effect Effects 0.000 claims 1
- 238000003384 imaging method Methods 0.000 abstract description 13
- 230000000694 effects Effects 0.000 abstract description 6
- 238000004445 quantitative analysis Methods 0.000 abstract 1
- 238000012546 transfer Methods 0.000 description 17
- 210000004185 liver Anatomy 0.000 description 11
- 239000000523 sample Substances 0.000 description 10
- 210000004204 blood vessel Anatomy 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 230000017531 blood circulation Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- 238000005406 washing Methods 0.000 description 6
- 210000001367 artery Anatomy 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000013139 quantization Methods 0.000 description 5
- 210000003484 anatomy Anatomy 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000002440 hepatic effect Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000033115 angiogenesis Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 3
- 210000000038 chest Anatomy 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 210000003240 portal vein Anatomy 0.000 description 3
- 241000415078 Anemone hepatica Species 0.000 description 2
- 206010027476 Metastases Diseases 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 230000007850 degeneration Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 208000019423 liver disease Diseases 0.000 description 2
- 230000009401 metastasis Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 230000017074 necrotic cell death Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 241000156948 Aphantopus hyperantus Species 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000004037 angiogenesis inhibitor Substances 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000009514 concussion Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000011393 cytotoxic chemotherapy Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 125000004051 hexyl group Chemical group [H]C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])C([H])([H])* 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 201000007270 liver cancer Diseases 0.000 description 1
- 208000014018 liver neoplasm Diseases 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 210000000496 pancreas Anatomy 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 201000001514 prostate carcinoma Diseases 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- 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/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
- A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5284—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving retrospective matching to a physiological signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- 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/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- 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/30096—Tumor; Lesion
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Quality & Reliability (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
描述了一种超声成像方法,其用于从图像的序列中的图像数据中消除诸如呼吸运动的生理运动的影响,所述图像的序列中的图像数据要用于量化分析,例如时间‑强度曲线计算。循环中的图像中的明亮反射物被识别并且所述明亮反射物的变化被用于识别运动的周期,所述明亮反射物反映运动影响,例如,图像与图像之间的亮度或位置的周期性变化。
Description
相关申请的交叉引用
本申请主张享有2013年1月17日提交的美国临时申请序列号61/753898的优先权,在此通过引用将其并入。
技术领域
本申请涉及一种超声诊断成像系统,并且具体而言,涉及使用超声诊断成像系统来评价肿瘤的治疗性处置的进展。
背景技术
国际专利公布WO 2006/090309(Bruce等人)描述了一种用于通过使用超声造影剂来检测肝脏中的病变的超声成像技术。一团造影剂被引入身体并且检测所述造影剂到达肝脏的时间。当一团造影剂行进通过身体的血管并且开始在身体中的特定器官或位置处出现时,图像中对比度的建立被称为造影剂的“洗入”。当造影剂的灌注在身体中的所述位置处到达平台并且然后其通过血流被从所述位置带走而下降,所述下降被称为造影剂的“洗出”。在前面提到的专利发表中,发明人利用了以下事实:到达肝脏的血流来自两个来源,即,肝动脉和门静脉。由于在血流的第一动脉相期间,血流将首先灌注HCC和转移的肝脏病变,发明人通过检测在血流的动脉相和随后的门相期间造影剂到达肝脏的时间来识别这样的病变。造影剂到肝脏的早期洗入的区可以是病变的症状。
一旦病变或转移通过该方法或其他方法而被识别出,则通常由医师开具处置方案。例如,治疗可以涉及高温/低温、细胞毒性化疗、或抗血管生成剂。治疗通常不是在单个疗程中执行的,而是在跨越数周或数月的时段的若干疗程中执行的。在每个治疗疗程中,医师通常期望评价治疗的进展以确定其对病人的有效性。例如,可以诊断性地对病变或转移进行成像以察看其是否正在收缩。但是经常是,处置的进展缓慢并且从前面的疗程以来仅在病变或转移中发生了小的变化。在这样的情况下,期望通过测量肿瘤的某些特性来定量地评价治疗的进展。一种这样的措施是肿瘤血管生成的退化(regression)。由于病变或转移随着其细胞的坏死而收缩,因而发展出的用于滋养病变的微血管将为病变提供更少的血液供应并且其自身可能会开始收缩。一种量化途径是评价该血管生成的退化、病变的微脉管的性能的下降。进行此的一种优选技术是利用造影剂,所述造影剂在对含有造影剂的血液的摄入期间明亮地照亮血流,以及在所述造影剂的洗出相期间亮度上的随后下降。对该再灌注的时间-强度曲线的测量使得能够对指向病变坏死的处置治疗的成功性进行精确评价。但是采集造影剂洗入和洗出的必要图像序列要求成像探头在造影剂的洗入和洗出时段期间(其可能持续长达100秒或更长)稳定地对图像平面中的病变进行成像。然而,肝脏邻近身体中的膈肌,引起患者的呼吸运动从而引起在成像过程期间解剖结构关于探头的图像采集平面而移动。该移动能够引起病变在过程期间移入和移出成像平面。因此,图像中病变及其造影剂的亮度将随着病变移动而变化,而不是仅由于造影剂洗入和洗出而变化。因此,期望在采集用于灌注评价的图像时消除这些不需要的运动的影响。
美国专利8,529,453(Averkiou等人)中描述了用于消除运动的影响的技术。该技术涉及采集目标解剖结构(在该情况下是肝脏)的图像,其也在图像中示出膈肌的部分。图像中的膈肌的尺寸和位置被用作界标。包括处于期望的尺寸和位置的膈肌的图像被用于灌注评价,而在其中膈肌改变了其尺寸和/或位置或其完全缺失的图像从随后的处理中被删除。这在效果上仅保留了病变在其中被示出在一致的成像平面的那些图像用于处理。但是从序列中识别和删除图像是手动过程,需要对每个图像以及它的膈肌的视图的单独的检查和评价。因此,期望以自动的过程消除运动的影响,所述过程既稳定又需要对图像的较低强度的用户分析。
发明内容
根据本发明的原理,描述了用于处理医学图像以使得由于特定生理功能(诸如由于呼吸或心脏运动的周期性移动)的移动能够被识别并且从后续处理图像中被选择性地移除的方法。这优选地通过仅处理属于运动周期的相同范围的阶段的那些图像帧来完成。通过仅处理在运动周期的相同阶段范围期间采集的那些图像帧,图像数据更加一致地受运动影响并且不同的运动的影响大部分被从数据中消除,允许对诸如灌注的解剖学性能的更加精确的量化评价。所述方法的临床应用包括肝脏、肾脏、胰腺、乳房、以及前列腺癌症评价。
在一种实现方式中,在运动的一个并且优选多个周期中采集图像的序列(循环)。例如,成人在休息时的正常呼吸率为每分钟8-16次呼吸,所以一分钟循环的采集将采集跨过呼吸运动的约8-16个循环的图像。采集的图像利用其采集时间被标记,使得这些时间能够与运动的阶段相关。在图像中识别强反射结构并且它们的运动遍及所述循环被跟踪,这通过它们所述循环中的位置和/或亮度的周期性变化来识别。例如,反射物可以开始在图像中的一个位置,跨图像的第一系列而移动到另一位置,然后跨图像的第二系列回到其原始位置,采集图像的两个系列的时间为运动的一个周期的时间。作为第二个范例,反射物可以以图像中的给定亮度开始,在图像的第一系列期间亮度增加,然后在图像的第二系列期间亮度下降到其原始强度。同样,采集图像的两个系列所需要的时间是运动周期的时间。所识别的运动周期的一个或多个时间可以与所研究的生理功能的典型周期时间进行比较以察看它们是否相关。例如,如果所识别的运动周期在每分钟5-20次呼吸的期望的范围之外,则周期测量结果将作为呼吸周期而被丢弃并且进行另一周期测量。
一旦识别出循环的运动周期,则在运动周期的期望阶段期间采集的那些图像被接受用于进一步处理,同时剩余的被丢弃。例如,可能期望仅接受在呼吸周期的阶段的25%范围期间采集的那些图像用于处理。在每个周期的呼吸阶段的该范围中采集的图像然后被接受用于处理以进行量化测量(诸如时间-强度曲线)和灌注测量。
附图说明
在附图中:
图1以方框图的形式图示了根据本发明的原理构造的超声诊断成像系统。
图2和图3图示了在针对灌注评价计算时间-强度曲线时图1中的QLab处理器的操作的细节。
图4a、图4b和图4c图示了随着图像平面由于呼吸运动而移动时肿瘤和血管的外观的变化。
图5以方框图的形式图示了根据本发明的原理如由图1中的超声系统执行的从图像特征量化中移除呼吸运动影响的第一方法。
图6以方框图的形式图示了根据本发明的原理如由图1中的超声系统执行的从图像特征量化中移除呼吸运动影响的第二方法。
图7a、图7b和图7c是根据本发明图示用于提取呼吸曲线的运动信息矩阵的形成的超声图像。
图8图示了根据本发明的使用呼吸曲线的呼吸门控。
图9图示了根据本发明的利用呼吸门控的和未利用呼吸门控处理的肝脏转移的时间-强度曲线。
具体实施方式
首先参考图1,以方框图的形式示出了根据本发明的原理构造的超声系统。该系统通过扫描正利用超声发射射束而被成像的身体的二维或三维区域而操作。随着每个射束被发射沿其导向路径通过身体,所述射束返回具有对应于所发射的频率分量的线性和非线性(基本频率和谐波频率)分量的回波信号。发射信号通过由射束遇到的造影剂微泡的非线性响应而被调制,从而生成具有谐波分量的回波信号。
图1中的超声系统利用在期望的射束方向上发射选定的调制特性的波或脉冲的发射器16,以从身体内的散射物返回谐波回波分量。所述发射器对若干控制参数做出响应,所述控制参数确定发射射束的特性,包括发射射束的频率分量、它们的相对强度或幅度、以及发射信号的相位或极性。所述发射器通过发射/接收开关14而被耦合到超声探头10的阵列换能器12的元件。所述阵列换能器能够是用于平面(二维)成像的一维阵列,或者是用于二维或体积(三维)成像的二维阵列。二维矩阵阵列也能够以xMatrix模式操作,在所述xMatrix模式中,由矩阵阵列探头来扫描单个平面(xPlane),所述单个平面能够在体积区域中被导向。
换能器阵列12接收来自身体的回波,所述回波包含处于换能器通带中的基本(线性)频率分量和谐波(非线性)频率分量。这些回波信号通过开关14被耦合到射束成形器18,射束成形器18合适地将来自不同换能器元件的回波信号延时,然后将它们组合以形成来自沿射束方向从浅到更深的深度的线性和谐波信号的序列。优选地,所述射束成形器是数字射束成形器,其在数字化的回波信号上操作以产生从图像的近场到远场深度的离散相干数字回波信号的序列。所述射束成形器可以是多线射束成形器,其响应于单个发射射束而产生沿多个空间上不同的接收扫描线的回波信号的两个或更多个序列,其对3D成像特别有用。射束形成的回波信号被耦合到集合存储器(ensemble memory)22。
在图1中的超声系统中,使用不同的调制技术,在每个射束方向上发射多个波或脉冲,造成针对图像域中每个被扫描点处的多个回波的接收。对应于共同的空间位置的回波在本文中被称为回波的集合,并且被存储在集合存储器22中,回波的集合可以从所述集合存储器22被检索并且一起被处理。集合的回波通过非线性信号分离器24以各种方式组合以产生期望的非线性或谐波信号。例如,具有不同的相位或极性调制的两个脉冲能够被发射到图像域中的每个点。当由所述两个脉冲造成的回波由超声系统接收到并且被相加地组合,则不同的调制使得回波的基本频率分量相互抵消并且谐波分量相互增强。这将回波信号的谐波分量分离出来。备选地,当两个回波相减时,基本频率分量被增强并且谐波分量抵消。这将基本频率分离出以构建标准B模式图像。该调制被称为“脉冲反转”,并且能够由如在美国专利5,706,819(Hwang等人)、5,951,478(Hwang等人)、以及5,577,505(BrockFisher等人)中所描述的相位、极性或幅度调制来完成。
分离的信号通过滤波器30被滤波以进一步移除不想要的频率分量,然后经受由检测器32进行的B模式或多普勒检测。检测到的信号被耦合到非线性信号组合器34以减少图像纹斑内容。信号然后在图像处理器36中被处理以形成二维、三维、谱学、参数化、或其他期望的图像,并且所述图像然后被显示在显示器38上。检测到的不需要纹斑降低处理的基本(线性)信号被直接耦合到图像处理器36以进行图像形成和显示。
在图1的系统中,超声图像数据也被耦合到QLab图像处理器40以产生时间-强度曲线和造影剂洗入和洗出特性。由QLab处理器产生的时间-强度曲线和特性被耦合回图像处理器,在所述图像处理器中,它们可以以数字或图形方式连同超声图像一起被显示在显示器38上。适合产生时间-强度曲线的标准QLab处理器可从马萨诸塞州安德沃的飞利浦医疗保健获得。
标准QLab处理器产生公知的时间-强度曲线,也被称为灌注曲线或再灌注曲线。参见美国专利5,833,613(Averkiou等人)、国际专利公布WO2005/099579(Rafter)、以及国际专利公布WO 2005/054898(Garge等人)。如这些公开所示出,在造影剂到达身体中的位置期间,监测造影剂在组织中的点(图像中的点)处的堆积。随着造影剂洗入组织,造影剂在一个点处的量由从每个点处的造影剂微泡返回的回波的强度来指示,并且被呈现在由低功率(低MI)发射采集的图像的序列中。能够针对组织中的每个点的造影剂强度的这种堆积和在造影剂的洗出期间它的随后下降来形成时间-强度曲线,所述组织中的每个点逐帧地返回回波的时间序列。对正被观看的整个组织的时间-强度曲线的量化呈现能够通过利用在图像中的每个点处表示时间-强度曲线的参数的颜色来对解剖图像中的每个像素进行着色来形成。例如,Garg等人的申请示出了心肌的参数化图像的形成,其中,图像中每个像素的颜色表示心肌中每个点处通过时间-强度曲线获得的峰值水平。备选地,时间-强度曲线的斜率代替峰值能够用于指示再灌注的速率。也参见美国专利6,692,438(Skyba等人)。
在本发明的实现方式中,随着造影剂到达身体中的转移的位置、堆积并且随后洗出,在图像的序列(循环)上采集造影剂灌注回波数据。回波的强度值将因此从不存在造影剂的基线水平开始、然后上升、达到平台、并且随着造影剂洗出而下降。曲线拟合算法然后将该数据变化拟合到如下定义的误差函数:
I(t)=A[erf{(t-t0)/T}+I0]
其中,I(t)是时间t处的线性强度,A是基线偏移上的最大强度,T是洗入时间参数,其与洗入时间线性地成比例(例如,从5%到95%),I0是基线偏移,并且t0是时间偏移。洗入时间优选地从经拟合的曲线而非有噪声的图像数据提取。优选地,造影剂回波数据在该处理之前不经历数据压缩,使得数据保持它的采集的线性关系。另一个方法是将整个时间-强度曲线(而非仅洗入部分)拟合到合适的数学模型作为对数正态分布,例如,定义为:
其中,μ和σ是正态分布的平均偏差和标准偏差,从所述正态分布获得对数变换。所述曲线能够通过改变μ而水平地缩放并且通过改变σ而关于偏斜变化。曲线下的面积是A,t0是时间偏移,并且C是基线强度偏移。指数正态拟合的曲线被用于提取洗入时间。
图1中的超声系统能够计算转移的灌注的量化度量,被称为洗入时间比(WITR)的参数。WITR如图2中的方框图所示出地被计算。根据造影剂洗入和洗出期间的转移或病变的图像的时间序列,针对转移的ROI转移(转移中的感兴趣区域)计算造影剂强度值,如方框72所指示。如以上所解释,能够通过针对在洗入期间采集的序列的每个图像而组合转移ROI的像素值来计算这些值。在方框74中,针对组织的正常实质的ROI实质计算强度值。这可以通过追踪肝脏转移的图像中的正常组织区域并且使用该第二追踪中的正常组织灌注像素值来完成。因此,这些值是正常组织的灌注值。在方框76中,时间-强度曲线被拟合到ROI转移的灌注值,并且在方框78中,时间-强度曲线被拟合到ROI实质的灌注值。该拟合并不总是必要的,但是其给出对可以从时间-强度曲线中提取的各种血液动力学参数的更好估计。尽管能够直接根据数据测量这些参数,但是数据中的噪声能够干涉测量的准确性,因此干涉对曲线拟合的偏好。针对每条曲线寻找洗入时间参数WIT,例如,通过使用以上所描述的误差函数或对数正态分布。这分别确定了针对转移和正常实质两者的洗入时间参数,WIT转移和WIT实质。如由方框70所示,然后根据两个洗入参数通过将WIT转移除以WIT实质来计算洗入时间比WITR。通过正常组织的洗入时间参数来对WIT转移进行归一化的效果是降低或消除过程中的变量(例如团剂尺寸、心脏输出、以及超声系统设置)的影响,所述变量从一个治疗疗程到另一个可以是不同的。因此,能够在患者正被处置的数周或数月的时段内针对每个治疗疗程来产生如由其血管生成所指示的转移的生长或收缩的可比较的量化度量。
图3中图示了对转移血管生成的另一量化度量,其减少或消除了团剂注入速率、患者的心脏输出、或机器设置的变化的影响。如在方框76和78中所示,针对转移和实质的每个ROI而拟合时间-强度曲线。在方框82和84中,每条时间-强度曲线的范围被归一化。一个方便的归一化尺度是零到一。在方框80中,差异曲线ΔT-I曲线被计算为两条归一化曲线T-I曲线转移和T-I曲线实质的差。时间-强度曲线的这一产生的进一步的细节可以在美国专利8,529,453(Averkiou等)中找到,在此通过引用将其并入。
如以上所提及,造影剂洗入和洗出的典型时段能够持续长达100秒或更多。这意味着采集图像数据的医师必须稳定地将病变的相同图像维持100秒,使得每个强度值属于病变的相同区域。例如,如果探头在采集期间移动,则病变能够移出图像平面并且采集的数据不能使用。即使探头靠紧患者的身体而保持稳定,但是病变由于患者的呼吸运动还是可以相对于探头视场而移动。这由图4a、图4b和图4c中的超声图像所图示。这些是由图像上的白色箭头所指示的肝脏病变的图像,其中,主血管(门静脉和肝动脉)在图像中可见,在病变下方并且在圆圈2上方。这些图像是在超声医师试图将探头紧靠在患者的腹部保持稳定时被采集的,使得相同的组织从一个图像到下一个图像连续地在图像中示出。但是图像受患者呼吸的影响并且呼吸运动令病变和期望的图像移入和移出探头的图像平面。这通过在这些图像中看到的肝动脉的变化而被最好的图示。在图4a中,能够看出,通过门静脉的图像切割面是通过血管的中央区域的,并且其在图像中显得是大的。但是在图4b中,血管的那一中央部分已经移出了成像平面,并且其显得细得多,因为通过血管的不同的切割面己移入成像平面。在图4c中,血管勉强可见,因为其己基本完全移出探头的图像平面。相应地,由白色箭头指示的病变的尺寸和位置也从一幅图像到另一幅图像变化。克服呼吸运动的不利效果的一种方式是针对呼吸周期对图像采集进行门控。呼吸信号能够通过已知的方式采集,诸如环绕患者的胸部的具有应力或压力传感器的弹性带。另一技术是在患者的胸部两侧的传感器之间发射小的信号并且测量患者的胸部阻抗变化。这些和其他技术能够产生呼吸周期的周期性信号并且能够被用于针对呼吸周期的相同阶段来对图像的采集进行门控。
根据本发明的原理,呼吸门控是通过图像后处理来执行的,如由图5中的方框图所示出,图5图示了本发明的一种实现方式。在步骤50中,在一个或多个运动周期中采集相继的图像的序列。该图像循环可以是在以上提及的造影剂洗入和洗出的100+秒期间采集的图像序列,并且可以包括数百幅或者甚至数千幅图像。显而易见,手动地处理这许多图像将是痛苦的并且可能不准确的任务。在每分钟十三次呼吸的名义呼吸下,这么长的图像循环将覆盖超过20个呼吸周期。在每秒一次搏动的名义心率下,这么长的图像循环将覆盖约100心跳周期。因此,该循环将被两者生理现象中的任一个的运动的很多周期所污染。
在步骤52中,识别周期性运动。对此的一种方法可以参考图4中的图像来理解。可以看出,在这些图像中,由于圆圈2上方的主肝血管被填充具有造影剂,其是图像中非常亮的反射物,实际上其是图像中最大的亮结构。当图像的切割平面通过动脉的中央时,看到血管处于最大和最亮,如图4a中所示,并且在其他图像中没这么大和亮。对这些图像中的在肝动脉亮度以下的像素值的进行阈值处理或抽样将从图像的中央提取肝动脉的像素,并且像素值被求和或积分以指示每幅图像中该结构的亮度。门静脉每次完全移动进入探头图像平面的切割平面时,门静脉的峰值亮度的再现指示呼吸运动移动图像平面的周期性,并且运动的周期的持续时间是采集两幅这样的图像之间的时间,如由它们的时间标签所示。随着这样识别出运动周期,在步骤54中识别在运动的期望阶段或阶段范围采集的图像。例如,可以决定使用每个呼吸周期中每幅峰值亮度的图像(图4a)和峰值亮度图像之前和之后的五幅图像。每个呼吸周期中的这十一幅图像然后被选择用于灌注评价处理,如步骤56中所示。在呼吸周期的较不期望的阶段期间采集的图像然后被从处理中忽略。结果是,灌注测量应该相对不受呼吸运动的影响所影响。
图6以方框图的形式图示了本发明的优选的实现方式。使用该附图中从步骤50采集的图像序列,图像序列中包含的运动信息被提取到在这里所称的“运动信息矩阵”中。图7a-图7d图示了形成运动信息矩阵(MIM)的一种方式。图7a图示了肝脏的超声图像,并且图7b图示的图像为包括图7a中的图像的循环中的所有图像帧的均值。循环中的每个帧,例如图7a中的帧,被从帧的均值(图7b)中减去,并且特定阈值以上的像素强度被赋予为1的值,并且所有其他像素强度被赋予为零的值。这产生每幅图像中的明亮结构的二值图像,例如图7c示出的一幅。得到的二值图像被加和到一起以形成如图7d中所示的MIM,能够在图像循环中容易地识别并且潜在地受运动影响的结构图像。
下一步骤62是从循环中选择触发帧。触发帧为要被包括在最终处理的图像中的帧。例如,其可以是这样一幅图像,即病变如通过图像中的高强度或大尺寸而被清晰地定义的图像。或者其可以为几乎与循环中先前和/或随后采集的图像一致的图像,并且因此表现出小的运动污染或者没有运动污染。触发帧可以由用户通过检查来指定,或者可以通过自动方法来选择,诸如通过将每对相继的帧相减并且使用使用产生最小的差异图像的帧中的帧。用于运动识别的候选结构然后被从触发帧中提取用于运动曲线形成,如步骤64所示。该提取可以包括如上所讨论的根据触发帧产生的二值图像。从触发图像提取的结构与MIM中的结构相比较,以选择显著地受运动影响的一个。例如,可以选择在大的强度范围或位置范围上变化的MIM中的结构。变化的频率可以被与生理现象的频率相比较,以确认它的变化对应于生理学功能的变化。如果其不对应,则忽略该结构并且选择MIM中的另一结构用于运动识别。为了辅助选择,能够根据MIM中的结构表现运动的影响的可能性对进行它们进行排序。被发现为用于运动提取的最感兴趣区域的结构然后被用作二值掩模,所述二值掩模被用于提取呼吸或其他身体运动的曲线。结构的掩模被应用到循环中的每幅图像中的对应位置,并且从所述图像帧中提取的所述掩模中的一系列平均强度值表示周期性运动,所述一系列平均强度值表现循环期间的强度变化。这些平均强度值被绘制为时间的函数(或者,等效地,序列中的帧号)以产出近似呼吸曲线。初始绘制可以被拟合到平均曲线以产生平滑地变化的呼吸曲线,如图8中的呼吸曲线100所示。可以使用己知的曲线拟合技术来产生平滑变化的曲线,例如,如拟合到近似曲线的波峰和波谷的分段立方厄米插值多项式。
图8图示了根据600个图像帧的循环中的结构的呼吸运动平均强度变化形成的呼吸曲线100。能够看到,曲线100指示在所述600个图像帧的采集期间发生的21个呼吸周期。图顶部的小圈110指示循环的一些离散帧。对应于呼吸周期中的点,圆圈被示出在呼吸曲线上。阈值102被应用到曲线100以描绘每个呼吸周期的阶段范围,在此期间,图像帧要被接受用于灌注处理。例如,假设呼吸曲线的下限峰为患者完全吸气的点,并且上限峰为患者完全呼气的点。临床医师可以决定使用均在完全呼气的时间附近采集的图像帧,例如,恰在患者完全呼气之前或恰在患者完全呼气之后。所述系统然后应用如附图中所示的阈值102,所述阈值在该范例中描绘每个呼吸周期的该阶段范围,看到其在该范例中包括每个周期的约30%。在该呼吸阶段范围期间采集的图像帧,由圆圈110表示的帧,然后被选择用于处理,如图6中的步骤66所示。步骤68中的经处理的图像帧因此是均在呼吸运动的相同阶段范围期间采集的并且应该大体空间地彼此相对应的一组帧,其中,在经处理的图像中的每幅中,图像中的结构和病变在尺寸和位置上呈现为大体相同。这将产出对诸如时间-强度曲线产生和灌注评价的生理度量和功能的准确的量化。用于这样的量化的数据可以使用触发帧而从接受的图像中选择。病变在触发帧中被识别为感兴趣区域(ROI),并且从被接受的图像帧中的每个中的ROI的位置提取数据。该数据然后如以上以及如在欧洲专利EP 2234544B1中所描述的被处理以产生时间-强度曲线,在此通过引用将其内容并入。
前述技术被应用于从若干肝癌患者临床采集的图像循环。图9图示了在这些情况中的一个中产生的时间-强度曲线,具有呼吸阶段门控和不具有呼吸阶段门控两者。图像中的感兴趣区域的强度数据被拟合到对数正态函数以产生再灌注的时间-强度曲线。拟合的优度使用R平方值和均方根误差值来量化。时间-强度曲线120根据不具有门控的所有图像帧的数据产生,并且被视为受到在图像的ROI处取样的强度值140的高度振荡特性的影响,所述强度值140由于由解剖结构的呼吸运动引起的在病变之外的对图像的取样而包含大的尖峰。当根据本发明应用门控时,仅使用在运动的较低阶段范围周围采集的图像帧110。强度数据的震荡效应在很大程度上被减小,这是因为由于在被接受用于处理的图像帧110上病变ROI的尺寸和位置更加均匀,采样在病变内更均匀地发生。
附录A中详记了本发明的实现方式中使用的数学公式,解释并且其在本发明的方法中的用途。
附录A
Claims (17)
1.一种用于减少从对象的医学图像的序列采集的图像数据中生理运动的影响的方法,包括:
在影响目标区域的周期性生理运动的发生期间采集(50)医学图像的序列,所述医学图像的序列含有对象的所述目标区域的图像数据;
使用(52)所述医学图像的图像数据来识别所述生理运动的周期;
识别(54)在识别出的运动周期的阶段的特定范围期间采集的图像;并且
处理(56)识别出的图像的图像数据以评价所述目标区域的属性;
其特征在于,
基于对如下图像数据的变化进行的跟踪来识别所述生理运动的周期:对应于在所述医学图像中识别的至少一个反射结构的图像数据。
2.如权利要求1所述的方法,其中,采集步骤还包括对受呼吸运动或心脏运动影响的医学图像的循环进行采集。
3.如权利要求1所述的方法,其中,使用图像数据还包括使用所述图像中的每幅中的相同的对应位置的图像数据来识别所述生理运动的周期。
4.如权利要求3所述的方法,其中,使用图像数据还包括使用强度图像数据。
5.如权利要求1所述的方法,其中,处理还包括处理来自识别出的图像中的每幅的所述目标区域的图像数据。
6.如权利要求1所述的方法,其中,采集还包括在通过将所述目标区域移入和移出图像平面来影响所述目标区域的生理运动的发生期间采集所述图像平面的二维医学图像的序列。
7.如权利要求1所述的方法,其中,采集还包括在通过在采集期间改变所述目标区域在图像中的强度或位置来影响所述目标区域的生理运动的发生期间采集图像平面的二维医学图像的序列。
8.如权利要求2所述的方法,其中,使用图像数据还包括从所述图像数据提取运动信息矩阵(60),所述运动信息矩阵包括具有至少一个反射结构的各结构的图像,所述结构能够在所述医学图像的所述循环中被易于识别并潜在地受到运动影响。
9.如权利要求1所述的方法,其中,使用图像数据还包括选择在所述周期中的一个的阶段的范围期间采集的触发帧(62)。
10.如权利要求8所述的方法,其中,使用图像数据还包括选择在所述周期中的一个的阶段的范围期间采集的触发帧(62)。
11.如权利要求10所述的方法,其中,使用图像数据还包括使用所述触发帧的图像数据和所述运动信息矩阵的数据以及所述触发帧的结构的运动的周期的频率来选择所述图像中的结构,所述结构的图像数据要被用于运动曲线识别。
12.如权利要求11所述的方法,其中,使用图像数据还包括使用图像的所述序列上的所选择的结构在亮度或位置上的变化来识别运动曲线。
13.如权利要求1所述的方法,其中,使用所述图像的图像数据来识别运动的周期还包括将运动的周期的频率与周期性生理现象的频率的期望范围进行比较。
14.如权利要求1所述的方法,其中,处理图像数据还包括处理来自识别出的图像中的所述目标区域的图像数据以产生灌注的时间-强度曲线。
15.如权利要求14所述的方法,其中,处理图像数据还包括处理来自识别出的图像中含有病变的图像区域的造影剂图像数据以产生造影剂再灌注的时间-强度曲线。
16.如权利要求14所述的方法,其中,处理图像数据还包括处理来自识别出的图像中含有肿瘤的图像区域的造影剂图像数据以产生造影剂再灌注的时间-强度曲线。
17.一种用于减轻医学图像的时间序列的图像数据中呼吸运动的影响的方法,所述医学图像的时间序列的图像数据要被用于产生时间-强度曲线,所述方法包括:
采集(50)含有目标区域的图像数据的医学图像的序列,针对所述图像数据要产生时间-强度曲线,所述序列是在影响所述目标区域的呼吸运动的发生期间被采集的;
使用(52)所述医学图像的图像数据来识别所述呼吸运动的周期;
选择(66)识别出的呼吸运动周期的阶段的范围;并且
处理(68)所述序列的图像的图像数据,所述序列的所述图像是在所述呼吸运动周期的所选择的阶段的范围的阶段中的一个阶段处被采集的,
其特征在于,
基于对如下图像数据的变化进行的跟踪来识别所述呼吸运动的周期:对应于在所述医学图像中识别的至少一个反射结构的图像数据。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361753898P | 2013-01-17 | 2013-01-17 | |
US61/753,898 | 2013-01-17 | ||
PCT/IB2014/058293 WO2014111860A2 (en) | 2013-01-17 | 2014-01-15 | Eliminating motion effects in medical images caused by physiological function |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104968279A CN104968279A (zh) | 2015-10-07 |
CN104968279B true CN104968279B (zh) | 2018-08-10 |
Family
ID=50114440
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480005221.6A Active CN104968279B (zh) | 2013-01-17 | 2014-01-15 | 消除医学图像中由生理功能引起的运动影响 |
Country Status (7)
Country | Link |
---|---|
US (1) | US9576357B2 (zh) |
EP (1) | EP2945544B1 (zh) |
JP (2) | JP6297593B2 (zh) |
CN (1) | CN104968279B (zh) |
BR (1) | BR112015016804A2 (zh) |
RU (1) | RU2015134348A (zh) |
WO (1) | WO2014111860A2 (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105940431B (zh) * | 2014-01-23 | 2019-10-01 | 皇家飞利浦有限公司 | 使用对比增强的超声成像对颈动脉斑块的评估 |
US9519966B2 (en) * | 2015-04-22 | 2016-12-13 | King Fahd University Of Petroleum And Minerals | Method, system and computer program product for breast density classification using parts-based local features |
KR102525616B1 (ko) * | 2015-10-08 | 2023-04-26 | 삼성메디슨 주식회사 | 조영제 초음파 진단 장치 및 방법 |
GB201519985D0 (en) | 2015-11-12 | 2015-12-30 | Respinor As | Ultrasonic method and apparatus for respiration monitoring |
CN107133549B (zh) | 2016-02-29 | 2020-11-24 | 上海联影医疗科技有限公司 | Ect运动门控信号获取方法及ect图像重建方法 |
WO2018086003A1 (en) * | 2016-11-09 | 2018-05-17 | Edan Instruments, Inc. | Systems and methods for temporal persistence of doppler spectrum |
US11116479B2 (en) * | 2017-01-04 | 2021-09-14 | Koninklijke Philips N.V. | Time-based parametric contrast enhanced ultrasound imaging system and method |
US11419584B2 (en) * | 2017-06-28 | 2022-08-23 | Tokitae, LLC | Devices, systems, and methods for diagnosis of pulmonary conditions through detection of b-lines in lung sonography |
CN110019855B (zh) * | 2017-12-04 | 2023-07-21 | 沪江教育科技(上海)股份有限公司 | 根据时间轴控制动效的方法、装置、存储介质及电子设备 |
JP7084193B2 (ja) | 2018-04-10 | 2022-06-14 | ザイオソフト株式会社 | 医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム |
US11497475B2 (en) * | 2020-01-31 | 2022-11-15 | Caption Health, Inc. | Ultrasound image acquisition optimization according to different respiration modes |
US11816832B2 (en) | 2020-11-18 | 2023-11-14 | Canon Medical Systems Corporation | Devices, systems, and methods for medical imaging |
US11712215B2 (en) | 2021-04-13 | 2023-08-01 | Canon Medical Systems Corporation | Devices, systems, and methods for motion-corrected medical imaging |
WO2023025837A1 (en) * | 2021-08-27 | 2023-03-02 | Koninklijke Philips N.V. | Method for use in analysing ultrasound image data of a subject |
CN114403924B (zh) * | 2022-01-19 | 2024-02-02 | 复旦大学附属中山医院 | 一种基于超声造影评估aip激素治疗疗效的方法 |
TW202342510A (zh) | 2022-02-18 | 2023-11-01 | 英商Rq生物科技有限公司 | 抗體 |
WO2024067948A1 (en) * | 2022-09-26 | 2024-04-04 | University College Dublin | Method and apparatus |
CN116306937B (zh) * | 2023-03-22 | 2023-11-10 | 中航信移动科技有限公司 | 一种基于时间序列离线数据的规则提取方法、介质及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1190573A (zh) * | 1997-02-13 | 1998-08-19 | 通用电器横河医疗系统株式会社 | 识别被观察器官状态时间相位的方法和基于该方法的超声波诊断装置 |
CN1814322A (zh) * | 2005-01-31 | 2006-08-09 | 重庆微海软件开发有限公司 | 利用呼吸信号控制图像同步的方法及其装置 |
CN1819798A (zh) * | 2004-05-31 | 2006-08-16 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及超声波图像处理方法 |
CN1846231A (zh) * | 2003-07-21 | 2006-10-11 | 派昂公司 | 用于在描绘运动器官的一系列图像内识别最佳图像的方法和系统 |
CN101147687A (zh) * | 2006-09-22 | 2008-03-26 | 株式会社东芝 | 超声波图像取得装置、诊断图像显示方法及医疗装置 |
EP1772103B1 (en) * | 2005-10-07 | 2009-05-13 | Medison Co., Ltd. | Ultrasound imaging system for extracting volume of an object from an ultrasound image and method for the same |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0770352B1 (en) * | 1995-10-10 | 2004-12-29 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging with contrast agents |
US5706819A (en) | 1995-10-10 | 1998-01-13 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging with harmonic contrast agents |
US5833613A (en) | 1996-09-27 | 1998-11-10 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging with contrast agents |
US5577505A (en) | 1996-02-06 | 1996-11-26 | Hewlett-Packard Company | Means for increasing sensitivity in non-linear ultrasound imaging systems |
JP2000210283A (ja) * | 1999-01-25 | 2000-08-02 | Shimadzu Corp | 医用動画像の運動周期求出装置 |
US6692438B2 (en) | 2001-12-18 | 2004-02-17 | Koninklijke Philips Electronics Nv | Ultrasonic imaging system and method for displaying tissue perfusion and other parameters varying with time |
JP2003210456A (ja) * | 2002-01-21 | 2003-07-29 | Toshiba Corp | 時系列画像の処理装置 |
US7731660B2 (en) * | 2003-07-25 | 2010-06-08 | Siemens Medical Solutions Usa, Inc. | Phase selection for cardiac contrast assessment |
EP1692543A1 (en) | 2003-12-03 | 2006-08-23 | Koninklijke Philips Electronics N.V. | Ultrasonic imaging system and method for simultaneous display of blood flow and perfusion parameters |
EP1713400B1 (en) * | 2004-01-16 | 2012-12-26 | Boston Scientific Scimed, Inc. | Apparatus for medical imaging |
US20050215904A1 (en) * | 2004-03-23 | 2005-09-29 | Siemens Medical Solutions Usa, Inc. | Ultrasound breathing waveform detection system and method |
WO2005099579A1 (en) | 2004-04-16 | 2005-10-27 | Koninklijke Philips Electronics, N.V. | Automated myocardial contrast echocardiography |
US20060074315A1 (en) * | 2004-10-04 | 2006-04-06 | Jianming Liang | Medical diagnostic ultrasound characterization of cardiac motion |
JP4685458B2 (ja) * | 2005-01-18 | 2011-05-18 | 株式会社東芝 | 超音波診断装置 |
WO2006090309A2 (en) | 2005-02-23 | 2006-08-31 | Koninklijke Philips Electronics, N.V. | Ultrasonic diagnostic imaging system and method for detecting lesions of the liver |
JP4206107B2 (ja) * | 2006-07-05 | 2009-01-07 | アロカ株式会社 | 超音波診断装置 |
WO2008107905A2 (en) * | 2007-03-08 | 2008-09-12 | Sync-Rx, Ltd. | Imaging and tools for use with moving organs |
JP2009082181A (ja) * | 2007-09-27 | 2009-04-23 | Aloka Co Ltd | 超音波診断装置 |
EP2234544B1 (en) | 2008-01-23 | 2011-09-07 | Michalakis Averkiou | Respiratory-gated therapy assessment with ultrasonic contrast agents |
US9289191B2 (en) * | 2011-10-12 | 2016-03-22 | Seno Medical Instruments, Inc. | System and method for acquiring optoacoustic data and producing parametric maps thereof |
US20130172730A1 (en) * | 2011-12-29 | 2013-07-04 | Amit Cohen | Motion-Compensated Image Fusion |
KR20150036650A (ko) * | 2012-07-16 | 2015-04-07 | 미라빌리스 메디카 인코포레이티드 | 초음파 유도 치료를 위한 인간 인터페이스 및 장치 |
-
2014
- 2014-01-15 CN CN201480005221.6A patent/CN104968279B/zh active Active
- 2014-01-15 EP EP14704901.9A patent/EP2945544B1/en active Active
- 2014-01-15 BR BR112015016804A patent/BR112015016804A2/pt not_active IP Right Cessation
- 2014-01-15 JP JP2015553207A patent/JP6297593B2/ja active Active
- 2014-01-15 US US14/760,736 patent/US9576357B2/en active Active
- 2014-01-15 WO PCT/IB2014/058293 patent/WO2014111860A2/en active Application Filing
- 2014-01-15 RU RU2015134348A patent/RU2015134348A/ru not_active Application Discontinuation
-
2018
- 2018-02-21 JP JP2018028589A patent/JP6542928B2/ja active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1190573A (zh) * | 1997-02-13 | 1998-08-19 | 通用电器横河医疗系统株式会社 | 识别被观察器官状态时间相位的方法和基于该方法的超声波诊断装置 |
CN1846231A (zh) * | 2003-07-21 | 2006-10-11 | 派昂公司 | 用于在描绘运动器官的一系列图像内识别最佳图像的方法和系统 |
CN1819798A (zh) * | 2004-05-31 | 2006-08-16 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及超声波图像处理方法 |
CN1814322A (zh) * | 2005-01-31 | 2006-08-09 | 重庆微海软件开发有限公司 | 利用呼吸信号控制图像同步的方法及其装置 |
EP1772103B1 (en) * | 2005-10-07 | 2009-05-13 | Medison Co., Ltd. | Ultrasound imaging system for extracting volume of an object from an ultrasound image and method for the same |
CN101147687A (zh) * | 2006-09-22 | 2008-03-26 | 株式会社东芝 | 超声波图像取得装置、诊断图像显示方法及医疗装置 |
Non-Patent Citations (4)
Title |
---|
A posteriori respiratory gating in contrast ultrasound for assessment of hepatic perfusion;G Renault et al.;《PHYSICS IN MEDICINE AND BIOLOGY,INSTITUTEOF PHYSICS PUBLISHING》;20051231;第4465–4480页 * |
An Original Methodology for Quantitative Assessment of Perfusion in Small Animal Studies Using Contrast-Enhanced Ultrasound;MULE S et al.;《2007 ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY》;20071130;摘要,正文第Ⅱ部分A、B小节,Fig.1-Fig.4 * |
Automatic phase determination for retrospectively gated cardiac CT;R.Manzke et al.;《MEDICAL PHYSICS》;20041231;第3345-3362页 * |
Image-Based Cardiac Phase Retrieval in Intravascular Ultrasound Sequences;AURA HERNANDEZ-SABATE et al.;《IEEE TRANSACTIONS ON ULTRASONICS FERROELECTRICS AND FREQUENCY CONTROL》;20110131;第60-72页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2014111860A3 (en) | 2014-09-18 |
JP6297593B2 (ja) | 2018-03-20 |
US9576357B2 (en) | 2017-02-21 |
JP6542928B2 (ja) | 2019-07-10 |
EP2945544A2 (en) | 2015-11-25 |
RU2015134348A (ru) | 2017-02-20 |
EP2945544B1 (en) | 2018-11-07 |
WO2014111860A2 (en) | 2014-07-24 |
JP2016503707A (ja) | 2016-02-08 |
US20150371379A1 (en) | 2015-12-24 |
CN104968279A (zh) | 2015-10-07 |
BR112015016804A2 (pt) | 2017-07-11 |
JP2018118061A (ja) | 2018-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104968279B (zh) | 消除医学图像中由生理功能引起的运动影响 | |
US11562463B2 (en) | Anatomically intelligent echochardiography for point-of-care | |
CN101951839B (zh) | 利用超声造影剂的呼吸门控治疗评估 | |
Golemati et al. | Using the Hough transform to segment ultrasound images of longitudinal and transverse sections of the carotid artery | |
CN103889337B (zh) | 超声波诊断装置以及超声波诊断装置控制方法 | |
JP5627706B2 (ja) | 超音波撮像装置、超音波撮像プログラム | |
JP2012508053A (ja) | 時間的パラメトリック・コントラスト画像を使った超音波病変識別 | |
CN105940431B (zh) | 使用对比增强的超声成像对颈动脉斑块的评估 | |
CN102573647A (zh) | 用于监测肝脏治疗的肝脏血流的对比增强超声评估 | |
JP7041147B6 (ja) | 造影剤流の肝灌流を特徴付けるシステム及び方法 | |
KR20140126815A (ko) | 호흡 주기 동안 체내 장기의 변화를 추적하는 방법, 장치 및 시스템. | |
US20100324420A1 (en) | Method and System for Imaging | |
JP2019198389A (ja) | 超音波診断装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム | |
EP2782070A1 (en) | Imaging method and device for the cardiovascular system | |
JP5885234B2 (ja) | 疾患判定装置および超音波画像形成装置における画像解析方法 | |
KR102250086B1 (ko) | 의료 영상 정합 방법, 이를 포함하는 장치 및 컴퓨터 기록 매체 | |
CN110167448B (zh) | 基于时间的参数对比增强超声成像系统和方法 | |
Pieragostini | Test of a new ultrasound system for monitoring the edges of the inferior vena cava | |
Fukuda et al. | Development of a Ubiquitous Echographic Device for Evaluating Human Body Composition | |
Abreu | Automatic Location of Atherosclerotic Plaques by Ultrasounds | |
Jia | Speckle tracking for cardiac strain imaging in ultrasound imaging and contrast enhancement in photoacoustic imaging | |
Raunig | A multivariate method for ultrasound tissue segmentation for biomarker analysis of tumor growth |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |