CN108135566B - 一种增强超声造影图像的方法、系统及超声造影成像设备 - Google Patents
一种增强超声造影图像的方法、系统及超声造影成像设备 Download PDFInfo
- Publication number
- CN108135566B CN108135566B CN201680039006.7A CN201680039006A CN108135566B CN 108135566 B CN108135566 B CN 108135566B CN 201680039006 A CN201680039006 A CN 201680039006A CN 108135566 B CN108135566 B CN 108135566B
- Authority
- CN
- China
- Prior art keywords
- contrast
- information
- region
- tissue
- nonlinear
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 230000002708 enhancing effect Effects 0.000 title claims abstract description 20
- 239000002872 contrast media Substances 0.000 claims abstract description 97
- 238000002604 ultrasonography Methods 0.000 claims description 45
- 230000001629 suppression Effects 0.000 claims description 20
- 238000010586 diagram Methods 0.000 claims description 18
- 230000011218 segmentation Effects 0.000 claims description 12
- 230000004927 fusion Effects 0.000 claims description 10
- 239000000523 sample Substances 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000002592 echocardiography Methods 0.000 claims description 6
- 230000009022 nonlinear effect Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 4
- 230000005764 inhibitory process Effects 0.000 abstract 2
- 231100000732 tissue residue Toxicity 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 6
- 230000001965 increasing effect Effects 0.000 description 4
- 230000002401 inhibitory effect Effects 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 210000003743 erythrocyte Anatomy 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000010412 perfusion Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 239000002961 echo contrast media Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000011551 log transformation method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Images
Classifications
-
- G06T5/90—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- 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/10132—Ultrasound image
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
Abstract
一种增强超声造影图像的方法、系统和超声造影成像设备，根据从超声回波信号中获得的非线性的造影信息和线性的组织信息，计算非线性参量，得到非线性参量分布图，分别采用两极阈值对非线性参量分布图进行分割，每次划分都分出造影剂区域和组织残留区域，对造影剂区域进行增强处理，对组织残留区域抑制处理，从而可实现对绝对属于造影剂区域进行大程度的增强，对绝对属于组织残留的区域进行大程度的抑制，对介于造影剂和组织残留的区域增强和抑制的程度稍微弱些，最后再将两次分类和权重调整后的结果进行融合处理，使对造影剂的增强、对造影剂的抑制，对介于造影剂和组织残留的保守处理，都能达到比较理想的造影图像增强效果。
Description
技术领域
本发明涉及超声造影成像设备，具体涉及增强超声造影图像的方法及系统。
背景技术
超声造影成像与常规超声成像相比，能增强超声回波信号强度，显示常规超声无法看见的细小血管，随着临床对超声检查的要求越来越高，超声造影成像越来越多地被临床应用，且用来检查器官也越来越广，其在临床诊断中的地位越来越高。
超声造影成像是注入造影剂后再超声成像，是一种增强超声成像，造影剂是一种能增强超声回波信号强度的增强剂。超声造影剂的种类多样，应用较多的是一种大小和红细胞接近的空心球体样的微泡，由于大小和红细胞接近所以微泡能到达全身所有血管。由于是空心球样所以微泡的声阻抗差大，极大的增加反射回波信号强度，使微小血管或毛细血管能显像。大多情况下造影剂微泡通过静脉注入，随血流分布到全身，造影微泡与红细胞行为相似，检测或者跟踪微泡，能反映出正常和异常组织的血流灌注情况，超声造影成像技术也可以称之为检测微泡信号技术。
衡量造影图像的好坏，一般会从以下几个指标考虑，SNR(signal noise ratio)、CTR(contrast to tissue ratio)、时间分辨力、空间分辨力。SNR也可以理解为灵敏度或敏感性，就是对微泡的检测程度；CTR是造影图像特有的特征，是造影图像中造影剂信号和组织信号强度之比。理论上造影图像只是对微泡显影，没有注入造影剂时造影图像是没有内容的，但是实际情况下由于系统的相消性限制，即使在没有注入造影剂时造影图像也有少量显影，该少量显影称为组织残留，这些组织残留在造影剂微泡灌注后也存在，与微泡信号混在一起，影响用户对造影剂微泡的识别判断。CTR越大，组织残留越小，微泡强度和组织残留强度的差距越大，纯造影剂信号的动态范围更大，造影图像的层次更丰富，微泡的灌注和效果过程显示的更清楚，图像的对比分辨力也更好。时间分辨力和空间分辨力和常规超声图像的要求一样，分辨力越高越好。
随着超声造影的使用越来越广，临床对造影图像的要求也越来越高，不再是能看见微泡显影，而是希望能显示更多细节。这就要求提高造影图像质量以能够显示更多的细节。
在超声造影图像中，提高造影图像质量的关键是提高SNR和CTR。SNR是由探头灵敏度和系统平台决定的，SNR高，意味着能检测出更微小的微泡，而微泡越小信号幅度越小，有些小微泡的信号幅度和组织残留幅度相当甚至更小，因此，即使SNR很高，但如果CTR不够高，这些小微泡反映出的信号将淹没在组织残留信息中，在造影图像中也观察不到。因此提高造影图像质量最关键的是要尽量提高CTR。
发明内容
本发明要解决的技术问题是提供另一种对造影信息进行增强的方案，从而提高造影图像的CTR。
根据第一方面，一种实施例中提供一种增强超声造影图像的方法,包括：
获取信息，根据造影成像过程中的超声回波信号获得非线性的造影信息和线性的组织信息；
计算非线性参量，比较所述造影信息和组织信息得到一帧造影信息的非线性参量分布图，所述非线性参量是衡量超声波在媒质中传播时产生的非线性效应大小的参量；
第一阈值分割，采用第一阈值对非线性参量分布图进行分割，小于第一阈值的区域为组织残留区域，大于第一阈值的区域为造影剂区域，所述第一阈值使得非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧中对应的点属于组织残留信息；
形成第一图像，对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第一图像数据；
第二阈值分割，采用第二阈值对非线性参量分布图进行分割，小于第二阈值的区域为组织残留区域，大于第二阈值的区域为造影剂区域，所述第二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧中对应的点属于造影剂信息，所述第二阈值大于第一阈值；
形成第二图像，对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第二图像数据；
图像融合，将第一图像数据和第二图像数据进行融合，形成增强后的造影图像数据。
根据第二方面，一种实施例中提供一种超声造影图像增强系统,包括：
信息获取单元，用于根据造影成像过程中的超声回波信号获得非线性的造影信息和线性的组织信息；
非线性参量计算单元，用于比较所述造影信息和组织信息得到一帧造影信息的非线性参量分布图，所述非线性参量是衡量超声波在媒质中传播时产生的非线性效应大小的参量；
第一阈值分割单元，用于采用第一阈值对非线性参量分布图进行分割，小于第一阈值的区域为组织残留区域，大于第一阈值的区域为造影剂区域，所述第一阈值使得非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧中对应的点属于组织残留信息；
第一图像形成单元，用于对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第一图像数据；
第二阈值分割单元，用于采用第二阈值对非线性参量分布图进行分割，小于第二阈值的区域为组织残留区域，大于第二阈值的区域为造影剂区域，所述第二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧中对应的点属于造影剂信息，所述第二阈值大于第一阈值；
第二图像形成单元，用于对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第二图像数据；
图像融合单元，用于将第一图像数据和第二图像数据进行融合，形成增强后的造影图像数据。
根据第三方面，一种实施例中提供一种超声造影成像设备，包括：
探头，用于向组织发射超声波并接收组织反射回的超声回波；
前信号处理模块，对超声回波信号进行处理，从超声回波信号中提取反映造影剂微泡信息的非线性成分和反映组织解剖特征的线性成分；
第一信号处理模块，用于对线性成分进行处理；
组织图像生成模块，用于将第一信号处理模块的输出处理成组织图像数据；
第二信号处理模块，用于对非线性成分进行处理；
造影图像生成模块，用于将第二信号处理模块的输出处理成造影图像数据；
用于对造影信息进行增强的超声造影图像增强系统，超声造影图像增强系统如上述所述，以得到增强后的造影图像数据；
显示模块，用于对组织图像数据和增强后的造影图像数据进行可视化显示。
本发明实施例中，对造影信息帧的非线性参量分布图采用两极阈值分割，每次划分都分出造影剂区域和组织残留区域，对造影剂区域进行增强处理，对组织残留区域抑制处理，从而可实现对绝对属于造影剂区域进行大程度的增强，对绝对属于组织残留的区域进行大程度的抑制，对介于造影剂和组织残留的区域增强和抑制的程度稍微弱些，最后再将两次分类和权重调整后的结果进行融合处理，使对造影剂的增强、对造影剂的抑制，对介于造影剂和组织残留的保守处理，都能达到比较理想的造影图像增强效果。
附图说明
图1为超声造影成像设备的第一种结构示意图；
图2为超声造影成像设备的第二种结构示意图；
图3为超声造影成像设备的第三种结构示意图；
图4为一种实施例中超声造影图像增强系统的结构示意图；
图5为一种实施例中超声造影图像增强系统的流程图；
图6为一种实施例中采用直方图方式表示的非线性参量分布图；
图7为另一种实施例中采用二维方式表示的非线性参量分布图。
具体实施方式
下面通过具体实施方式结合附图对本发明作进一步详细说明。
如图1所示，超声造影成像设备包括探头10、前信号处理模块13、第一信号处理模块16、第二信号处理模块17、组织图像生成模块18、造影图像生成模块19、超声造影图像增强系统20和显示模块100。
探头10用于向人体组织11发射超声波并接收组织11反射回的超声回波，超声回波可存在信号接收模块12中。
前信号处理模块13用于对超声回波信号进行处理，提取有用的信息。前信号处理模块13主要包含放大器、波束合成器和滤波器，提取更多有用的、不同成分的信号。在超声造影成像设备中，前信号处理模块13需要提取两种不同的信号，一种是线性成分14，反映组织解剖特征；一种是非线性成分15，反映造影剂微泡信息。
第一信号处理模块16用于对线性成分进行处理，第二信号处理模块17用于对非线性成分进行处理。第一信号处理模块16和第二信号处理模块17的处理类似，主要包括解调、滤波、降采样、求包络、求Log、动态范围变换等，只是参数的选择不同。
组织图像生成模块18用于将第一信号处理模块16的输出处理成组织图像数据；造影图像生成模块19用于将第二信号处理模块17的输出处理成造影图像数据。
超声造影图像增强系统20用于对造影信息进行增强，以得到增强后的造影图像数据；显示模块100用于对组织图像数据和增强后的造影图像数据进行可视化显示。大多数造影图像旁会同时显示一幅组织图像，用于定位病灶位置。为了得到微泡的非线性特征，造影成像通常会在同一位置发射多次相位或者幅度不同的脉冲波，组织和造影图像都是用这多次发射波形中的部分进行成像，所以组织图像和造影图像显示的是同一时刻同一位置的信息，组织图像显示的是回波的线性特性，造影图像显示的是回波中的非线性特性。
在采用造影增强技术对造影图像进行增强时，需要通过阈值将图像划分成组织残留区域和造影剂区域，然后对划归为造影剂区域的点进行增强，但在对造影增强的继续研究中，发明人认识到组织残留区域和造影剂区域并没有明确的边界，任何阈值都无法正确的区分组织残留区域和造影剂区域。因此，在本发明实施例中，首先得到组织的线性信息和造影剂微泡的非线性信息，再比较非线性信息和线性信息得到非线性参量值，接着根据非线性参量值对造影图像进行两次分类，每次都分成组织残留和造影剂两类，对两类按照不同的权重增强或抑制，根据分类的不同，权重系数也不同；最后将两次分类和权重调整后的结果进行融合处理。
在具体实施例中，超声造影图像增强系统20可以对组织图像生成模块18和造影图像生成模块19输出的图像域的造影信息进行增强处理，如图1所示，超声造影图像增强系统18的输入端分别连接组织图像生成模块18和造影图像生成模块19，输出端连接到显示模块100，根据组织图像生成模块18和造影图像生成模块19输出的组织图像数据和造影图像数据对造影图像数据进行增强处理，并将增强结果输出到显示模块100。超声造影图像增强系统20还可以对前信号处理模块13输出的射频域的造影信息进行增强处理。如图2所示，超声造影图像增强系统20连接在前信号处理模块13和第二信号处理模块17之间，接收前信号处理模块13输出的非线性成分15和线性成分14，根据非线性成分15和线性成分14对非线性成分进行增强处理，并将增强结果输出到第二信号处理模块17。超声造影图像增强系统18还可以对第二信号处理模块中任一环节的造影信息进行增强处理。如图3所示，超声造影图像增强系统20的输入端分别连接第一信号处理模块16和第二信号处理模块17的任一环节，例如解调后信号，求包络后，Log变换后等，超声造影图像增强系统20的输出端连接造影图像生成模块19。
超声造影图像增强系统20的结构图如图4所示，包括信息获取单元201、非线性参量计算单元202、第一阈值分割单元203、第一图像形成单元204、第二阈值分割单元205、第二图像形成单元206和图像融合单元207。信息获取单元201用于根据造影成像过程中的超声回波信号获得非线性的造影信息和线性的组织信息，即造影信息和线性的组织信息是从同一超声回波信号中提取出的两个信息，非线性信息为造影信息，线性信息为组织残留信息；非线性参量计算单元202用于比较同一时刻的造影信息和组织信息得到一帧造影信息的非线性参量分布图，该非线性参量分布图称为该造影信息帧的非线性参量分布图，所述非线性参量是衡量超声波在媒质中传播时产生的非线性效应大小的参量；本实施例中，非线性参量分布图根据同一时刻的一帧非线性造影信息和一帧线性组织信息的比值或差值计算得出。第一阈值分割单元203用于采用第一阈值对非线性参量分布图进行分割，小于第一阈值的区域为组织残留区域，大于第一阈值的区域为造影剂区域，所述第一阈值使得非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧中对应的点属于组织残留信息；第一图像形成单元204用于对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域信息进行增强处理，并基于处理后的结果形成第一图像数据；第二阈值分割单元205用于采用第二阈值对非线性参量分布图进行分割，小于第二阈值的区域为组织残留区域，大于第二阈值的区域为造影剂区域，所述第二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧中对应的点属于造影剂信息，所述第二阈值大于第一阈值；第二图像形成单元206用于对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第二图像数据；图像融合单元207用于将第一图像数据和第二图像数据进行融合，形成增强后的造影图像数据。
在另一实施例中，超声造影图像增强系统20还包括噪声识别单元208，噪声识别单元208用于在非线性参量计算单元计算非线性参量之前判断造影信息中的噪声，将一帧造影信息分为噪声信号和除去噪声的非噪声造影信号；非线性参量计算单元202在计算非线性参量时采用非噪声造影信号计算非线性参量。
超声造影图像增强系统的具体处理流程如图5所示，包括以下步骤：
一、去噪
去噪的目的是还原信号自身大小，避免TGC(信号增益补偿)及噪声对后面处理带来影响。一般超声系统中用ATGC(analog time gain compensat ion,前端模拟信号增益补偿)、DTGC(digital time gain compensat ion,后端数字域增益补偿)用来补偿信号传播过程中深度方向的衰减来补偿信号沿传播方向的信号衰减，如果没有ATGC和DTGC的作用，探头空打时的噪声是近似为白噪声，加上了ATGC和DTGC的影响后，噪声均值曲线是随深度递增的曲线。为了得到信号真实的大小，减去噪声均值曲线不仅仅是去噪，也消除了ATGC和DTGC对信号的影响。
步骤21到24是实现去噪。去噪首先要获得噪声信号，本实施例中获得噪声信号的方法是暂停发射，在超声造影图像增强系统20前的环节执行步骤21和22，缓存一帧非线性造影噪声信号和一帧线性组织噪声信号，即控制探头不发射，然后根据回波得到一帧非线性造影噪声信号和一帧线性组织噪声信号。处理的具体步骤是计算一帧噪声信号每行的均值得到噪声纵向均值曲线，再经过低通滤波得到光滑的噪声均值曲线。
在步骤23中，将超声造影图像增强系统20输入的造影信息减去非线性造影噪声信号的噪声均值曲线，得到去噪后的造影信号。
在步骤24中，将超声造影图像增强系统20输入的组织信息减去线性组织噪声信号的噪声均值曲线，得到去噪后的组织信号。
二、检测噪声
步骤25，检测去噪后造影图像中的噪声，将一帧造影图像划分成噪声区域26和非噪声造影剂区域27。噪声检测可用阈值分割的方法，小于阈值的就是噪声，阈值是根据系统特性经验确定的。
三、计算非线性参量
本实施例中，非线性参量是衡量超声波在媒质中传播时产生的非线性效应大小的重要参量。
步骤28中，比较同一时刻的一帧造影信息和一帧组织信息，得到该造影信息帧的非线性参量分布图。理论上，非线性参量为媒质物态方程泰勒展开式中多次项系数与线性系数的比值或差值，例如，媒质物态方程泰勒展开式中，线性系数为A，二次项系数为B，三次项系数为C，非线性参量可以是B/A、C/A或(B+C)/A，及B、C其他组合与A的比值。本实施例中，非线性参量通过计算同一时刻造影信号和组织信号的比值或差值得出。对于一帧造影信息，其由M×N个像素值组成，每个像素值对应显示界面上的一个像素点，该像素值可能属于造影剂产生的信号，也可能属于组织残留信号。对于一帧造影信息而言，每个像素点都会有对应的一个非线性参量，如果该像素值是组织残留信息，则其对应的非线性参量值较小，如果该像素值是造影剂信息，则其对应的非线性参量值较大，因此，通过计算造影信息帧中各像素值的非线性参量，即可将属于组织残留信息的像素值和属于造影剂信息的像素值区分开。根据同一时刻的一帧非线性造影信息和一帧线性组织信息的非线性参量即可得出该造影信息帧中各像素点的非线性参量，对这些非线性参量进行统计，则可得到该造影信息帧的非线性参量分布图，得出的非线性参量分布图可以采用直方图方式表示，也可以采用二维图方式表示。
四、两极阈值划分
如图6所示为采用直方图方式表示的非线性参量分布图，包括一个大峰和一个小峰，直方图的横坐标就是计算出的非线性参量值，单位是dB，纵坐标表示造影信息帧中与该非线性参量值对应的点的数量。经试验证明，大峰区域中绝大部分区域所对应的造影信息帧中点反映的是造影剂微泡信息，小峰区域中绝大部分区域所对应的造影信息帧中点反映的是组织残留信息。
在确定两极阈值的另一实施例中，非线性参量分布图为二维图，如图7所示，其以色彩表示造影信息帧中各点非线性参量值，深蓝色的表示噪声区域，非线性参量值由大到小采用的颜色表示为红、黄、绿、蓝，因组织残留的非线性参量值较小，造影剂的非线性声参量值较大，故可以确定约小于-40dB的一定是组织残留，约大于-20dB的一定是造影剂，这里-40dB就是较小阈值，-20dB就是较大阈值。图7只是一个具体的例子，实际中，需要根据系统平台，通过较多非线性参量二维图，找出确定为组织残留的点的非线性参量值大小，其最大的值就是较小阈值，找出确定为造影剂的点的非线性参量值的大小，其最小值就是较大阈值。
本实施例中，通过两极阈值来划分非线性参量分布图，每次划分都分为造影剂区域和组织残留区域，对造影剂区域进行增强显示，对组织残留区域抑制显示。两极阈值就是极小、极大两个阈值，非线性参量不同，阈值的选择也不同，阈值的大小和系统的相消性有关，大阈值要大于系统的相消值，小阈值要小于系统的相消值，相消值简单讲就是注入造影剂最强时刻和没有注入造影剂时两个非线性成分大小的差值。具体如下：
在步骤29中，采用第一阈值对非线性参量分布图进行分割，如图6中，左边的直线阈值为第一阈值，该第一阈值使得非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧中对应的点属于组织残留信息，因此第一阈值也称为极小阈值。划分后可认为：在如图6所示的非线性参量分布图中，小于第一阈值的左边的区域为组织残留区域32，组织残留区域32中的点在造影信息帧中的对应点被认为属于组织残留信息。大于第一阈值的右边的区域为造影剂区域31，造影剂区域31中的点在造影信息帧中的对应点被认为可能属于造影剂产生的信息。此时的小于第一阈值的几乎都是组织残留，而大于第一阈值的可能属于造影剂。再结合噪声区域26，这样就把造影图像分成噪声、造影剂和组织残留三个不同的区域。
步骤35，为了达到增强造影剂、抑制组织残留的目的，对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理。例如对上述三个区域分别乘以或加上不同的权重系数P1、P2和P3，如果是乘以权重系数，则造影剂区域的系数P2大于1，组织残留区域的系数P3应小于1，噪声区域系数P1可以等于1，如果想抑制噪声，P1也可以小于1。同时，P2大于1，也增加了造影信号的动态范围。如果是加上权重系数，则P3<0<P2，这样能增强造影信号强度，减弱组织残留强度，P1≤0，P1<0也达到抑制噪声的目的。对上述三个区域分别乘以或加上不同的权重系数后再形成一个新的图像，称为第一图像数据。
在步骤30中，采用第二阈值对非线性参量分布图进行分割，如图6中，右边直线阈值为第二阈值，该第二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧中对应的点属于造影剂信息，因此，第二阈值也称为极大阈值，且第二阈值大于第一阈值。划分后可认为：在如图6所示的非线性参量分布图中，小于第二阈值的区域为组织残留区域34，组织残留区域34中的点在造影信息帧中的对应点被认为可能属于组织残留信息。大于第二阈值的区域为造影剂区域33，即造影剂区域33中的点在造影信息帧中的对应点被认为属于造影剂产生的信息。此时的大于第二阈值的几乎都是造影剂，而小于第二阈值的则可能属于组织残留。再结合噪声区域26，也将造影图像分成噪声、造影剂和组织残留三个不同的区域。
步骤36，为了达到增强造影剂、抑制组织残留的目的，对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理。例如对步骤30中的三个区域分别乘以或加上不同的权重系数P1、P4和P5。同样如果乘以权重系数，则造影剂区域的系数P4大于1，组织残留区域的系数P5应小于1，即P5<1<P4，而噪声区域系数P1等于或小于1。如果加上权重系数，则P5<0<P4，P1≤0。对上述三个区域分别乘以或加上不同的权重系数后再形成一个新的图像，称为第二图像数据。
为了达到增强造影图像的效果，对绝对属于造影剂区域的像素点增强的程度大，对绝对属于组织残留区域的像素点抑制的程度大，对介于造影剂和组织残留的区域增强和抑制的程度稍微弱些。如果乘以权重系数，增强和抑制的参数满足P4>P2>1，P3<P5<1；如果加上权重系数，增强和抑制的参数满足P4>P2>0，P3<P5<0。
五、图像融合
因为任何阈值都无法完全正确的区分造影剂和组织残留，所以第一图像数据和第二图像数据都有些极端，作为最终的增强图像都有些不正确。因此在接下来的步骤37中将两种图像进行融合处理，融合的方法有很多，本实施例提出了一种简单的融合方法，将第一图像数据乘以权重系数K，第二图像数据乘以权重系数1-K，0<K<1,再组合显示。为了达到增强造影剂的目的，1-K≥K，即K≤0.5。
本实施例最后通过对两个图像数据进行融合处理，相互补充，使对造影剂的增强、对造影剂的抑制，对介于造影剂和组织残留的保守处理，都能达到比较理想的造影图像增强效果。
本领域技术人员可以理解，上述实施方式中各种方法的全部或部分步骤可以通过程序来指令相关硬件完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器、随机存储器、磁盘或光盘等。
以上应用了具体个例对本发明进行阐述，只是用于帮助理解本发明，并不用以限制本发明。对于本领域的一般技术人员，依据本发明的思想，可以对上述具体实施方式进行变化。
Claims (17)
1.一种增强超声造影图像的方法,其特征在于包括：
获取信息，根据造影成像过程中的超声回波信号获得非线性的造影信息和线性的组织信息；
计算非线性参量，比较所述造影信息和组织信息得到一帧造影信息的非线性参量分布图，所述非线性参量是衡量超声波在媒质中传播时产生的非线性效应大小的参量；
第一阈值分割，采用第一阈值对非线性参量分布图进行分割，小于第一阈值的区域为组织残留区域，大于第一阈值的区域为造影剂区域，所述第一阈值使得非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧中对应的点属于组织残留信息；
形成第一图像，对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第一图像数据；
第二阈值分割，采用第二阈值对非线性参量分布图进行分割，小于第二阈值的区域为组织残留区域，大于第二阈值的区域为造影剂区域，所述第二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧中对应的点属于造影剂信息，所述第二阈值大于第一阈值；
形成第二图像，对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第二图像数据；
图像融合，将第一图像数据和第二图像数据进行融合，形成增强后的造影图像数据。
2.如权利要求1所述的方法，其特征在于，所述增强处理包括乘以大于1的权重系数或加上大于0的权重系数，所述抑制处理包括乘以小于1的权重系数或加上小于0的权重系数。
3.如权利要求2所述的方法，其特征在于，形成第一图像数据过程中进行抑制处理的权重系数小于形成第二图像数据过程中进行抑制处理的权重系数，形成第二图像数据过程中进行增强处理的权重系数大于形成第一图像数据过程中进行增强处理的权重系数。
4.如权利要求1所述的方法，其特征在于，在计算非线性参量之前还包括：判断造影信息中的噪声，将一帧造影信息分为噪声信号和除去噪声的非噪声造影信号；在计算非线性参量时采用非噪声造影信号计算非线性参量。
5.如权利要求1所述的方法，其特征在于，所述非线性参量分布图根据同一时刻的一帧非线性造影信息和一帧线性组织信息的比值或差值计算得出。
6.如权利要求1至5中任一项所述的方法，其特征在于，所述图像融合包括：将第一图像数据和第二图像数据分别乘以第一权重系数和第二权重系数，所述第一权重系数与第二权重系数的和为1且第二权重系数大于或等于第一权重系数。
7.一种超声造影图像增强系统,其特征在于包括：
信息获取单元，用于根据造影成像过程中的超声回波信号获得非线性的造影信息和线性的组织信息；
非线性参量计算单元，用于比较所述造影信息和组织信息得到一帧造影信息的非线性参量分布图，所述非线性参量是衡量超声波在媒质中传播时产生的非线性效应大小的参量；
第一阈值分割单元，用于采用第一阈值对非线性参量分布图进行分割，小于第一阈值的区域为组织残留区域，大于第一阈值的区域为造影剂区域，所述第一阈值使得非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧中对应的点属于组织残留信息；
第一图像形成单元，用于对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第一图像数据；
第二阈值分割单元，用于采用第二阈值对非线性参量分布图进行分割，小于第二阈值的区域为组织残留区域，大于第二阈值的区域为造影剂区域，所述第二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧中对应的点属于造影剂信息，所述第二阈值大于第一阈值；
第二图像形成单元，用于对该造影信息帧中属于组织残留区域的信息进行抑制处理，对该造影信息帧中属于造影剂区域的信息进行增强处理，并基于处理后的结果形成第二图像数据；
图像融合单元，用于将第一图像数据和第二图像数据进行融合，形成增强后的造影图像数据。
8.如权利要求7所述的系统，其特征在于，所述增强处理包括乘以大于1的权重系数或加上大于0的权重系数，所述抑制处理包括乘以小于1的权重系数或加上小于0的权重系数。
9.如权利要求8所述的系统，其特征在于，第一图像形成单元进行抑制处理的权重系数小于第二图像形成单元进行抑制处理的权重系数，第二图像形成单元进行增强处理的权重系数大于第一图像形成单元进行增强处理的权重系数。
10.如权利要求7所述的系统，其特征在于，还包括噪声识别单元，所述噪声识别单元用于在非线性参量计算单元计算非线性参量之前判断造影信息中的噪声，将一帧造影信息分为噪声造影和除去噪声的非噪声造影信号；非线性参量计算单元在计算非线性参量时采用非噪声造影信号计算非线性参量。
11.如权利要求7所述的系统，其特征在于，所述非线性参量计算单元根据同一时刻的一帧非线性造影信息和一帧线性组织信息的比值或差值计算该造影信息帧的非线性参量分布图。
12.如权利要求7至11中任一项所述的系统，其特征在于，所述图像融合单元将第一图像数据和第二图像数据分别乘以第一权重系数和第二权重系数后进行融合，所述第一权重系数与第二权重系数的和为1且第二权重系数大于或等于第一权重系数。
13.如权利要求7至11中任一项所述的系统，其特征在于，所述非线性参量分布图为直方图或二维图。
14.一种超声造影成像设备，其特征在于包括：
探头，用于向组织发射超声波并接收组织反射回的超声回波；
前信号处理模块，对超声回波信号进行处理，从超声回波信号中提取反映造影剂微泡信息的非线性成分和反映组织解剖特征的线性成分；
第一信号处理模块，用于对线性成分进行处理；
组织图像生成模块，用于将第一信号处理模块的输出处理成组织图像数据；
第二信号处理模块，用于对非线性成分进行处理；
造影图像生成模块，用于将第二信号处理模块的输出处理成造影图像数据；
如权利要求7至13中任一项所述的超声造影图像增强系统，用于对造影信息进行增强，以得到增强后的造影图像数据；
显示模块，用于对组织图像数据和增强后的造影图像数据进行可视化显示。
15.如权利要求14所述的超声造影成像设备，其特征在于，所述超声造影图像增强系统对前信号处理模块输出的射频域的造影信息进行增强处理。
16.如权利要求14所述的超声造影成像设备，其特征在于，所述超声造影图像增强系统对组织图像生成模块和造影图像生成模块输出的图像域的造影信息进行增强处理。
17.如权利要求14所述的超声造影成像设备，其特征在于，所述超声造影图像增强系统对第二信号处理模块处理过程中的造影信息进行增强处理。
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2016/087979 WO2018000359A1 (zh) | 2016-06-30 | 2016-06-30 | 一种增强超声造影图像的方法、系统及超声造影成像设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108135566A CN108135566A (zh) | 2018-06-08 |
| CN108135566B true CN108135566B (zh) | 2020-11-24 |
Family
ID=60785845
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201680039006.7A Active CN108135566B (zh) | 2016-06-30 | 2016-06-30 | 一种增强超声造影图像的方法、系统及超声造影成像设备 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN108135566B (zh) |
| WO (1) | WO2018000359A1 (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111200699B (zh) * | 2018-11-19 | 2022-04-26 | 瑞昱半导体股份有限公司 | 影像调整方法 |
| CN116433475A (zh) * | 2019-06-26 | 2023-07-14 | 图码思(成都)科技有限公司 | 一种基于特征点提取的图像拼接方法及终端 |
| CN113164160B (zh) * | 2020-05-18 | 2022-09-16 | 深圳迈瑞生物医疗电子股份有限公司 | 超声造影成像方法、装置和存储介质 |
| CN112837254A (zh) * | 2021-02-25 | 2021-05-25 | 普联技术有限公司 | 一种图像融合方法、装置、终端设备及存储介质 |
| CN113538299B (zh) * | 2021-09-13 | 2022-01-11 | 深圳瀚维智能医疗科技有限公司 | 超声图像去噪方法、装置、设备及计算机可读存储介质 |
| CN113768547B (zh) * | 2021-09-14 | 2024-03-22 | 南京超维景生物科技有限公司 | 冠状动脉成像方法及装置、存储介质及电子设备 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4662380A (en) * | 1985-10-31 | 1987-05-05 | General Electric Company | Adaptive time gain compensation system for ultrasound imaging |
| US6508767B2 (en) * | 2000-02-16 | 2003-01-21 | Koninklijke Philips Electronics N.V. | Ultrasonic harmonic image segmentation |
| US6371914B1 (en) * | 2000-04-13 | 2002-04-16 | Bracco Research S.A. | Single-shot phase cancellation ultrasound contrast imaging |
| US6592523B2 (en) * | 2001-11-21 | 2003-07-15 | Ge Medical Systems Global Technology Company, Llc | Computationally efficient noise reduction filter for enhancement of ultrasound images |
| US6679844B2 (en) * | 2002-06-20 | 2004-01-20 | Acuson Corporation | Automatic gain compensation for multiple mode or contrast agent imaging |
| CN101987023B (zh) * | 2009-07-31 | 2014-05-07 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像增益补偿及图像优化方法及其装置和系统 |
| CN103845077B (zh) * | 2012-12-05 | 2016-01-20 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像增益优化方法及超声成像增益自动优化装置 |
| CN104720850B (zh) * | 2013-12-23 | 2017-10-03 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声造影成像方法及造影图像的区域检测、显像方法 |
| CN105631867B (zh) * | 2015-12-25 | 2019-08-23 | 中国科学院深圳先进技术研究院 | 一种全自动超声造影图像分割方法 |
-
2016
- 2016-06-30 CN CN201680039006.7A patent/CN108135566B/zh active Active
- 2016-06-30 WO PCT/CN2016/087979 patent/WO2018000359A1/zh active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| CN108135566A (zh) | 2018-06-08 |
| WO2018000359A1 (zh) | 2018-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108135566B (zh) | 一种增强超声造影图像的方法、系统及超声造影成像设备 | |
| US10219786B2 (en) | Methods for optimizing gain of ultrasound images and automatic gain optimization apparatuses for ultrasound imaging | |
| Narayanan et al. | A view on despeckling in ultrasound imaging | |
| Aysal et al. | Rayleigh-maximum-likelihood filtering for speckle reduction of ultrasound images | |
| JP3432204B2 (ja) | 超音波診断装置 | |
| US11786216B2 (en) | Ultrasound contrast enhanced imaging method and ultrasound imaging system | |
| JP2020517363A (ja) | 改善されたクラッタ抑制を含むパワードップラー・イメージング・システム及び方法 | |
| Scorza et al. | A novel approach to a phantom based method for maximum depth of penetration measurement in diagnostic ultrasound: a preliminary study | |
| US8500646B2 (en) | Color Doppler mode image processing in an ultrasound system | |
| Park et al. | Combined chirp coded tissue harmonic and fundamental ultrasound imaging for intravascular ultrasound: 20–60 MHz phantom and ex vivo results | |
| WO2006015264A2 (en) | T-statistic method for suppressing artifacts in blood vessel ultrasonic imaging | |
| WO2015096352A1 (zh) | 一种超声造影成像方法及造影图像的区域检测、显像方法 | |
| US20080097212A1 (en) | Automatic identification of orientation in medical diagnostic ultrasound | |
| EP2757523A2 (en) | Method and apparatus for eliminating noise in magnetic resonance images | |
| WO2016161574A1 (zh) | 一种超声造影成像方法及装置 | |
| CN109363722B (zh) | 彩色血流成像中运动伪像的抑制方法及设备 | |
| Cowe et al. | Improving performance of pulse compression in a Doppler ultrasound system using amplitude modulated chirps and Wiener filtering | |
| KR101652728B1 (ko) | 초음파 영상 화질 개선 방법 및 이를 이용한 초음파 영상 장치 | |
| Ahmed et al. | Comparison of windowing effects on elastography images: Simulation, phantom and in vivo studies | |
| CN110599405B (zh) | 一种超声图像增强方法、装置和计算机设备 | |
| US11129598B2 (en) | Calibration for ARFI imaging | |
| Hourani et al. | Block-wise ultrasound image deconvolution from fundamental and harmonic images | |
| WO2020142760A1 (en) | Systems and methods for ultrasound attenuation coefficient estimation | |
| Hourani et al. | Joint deconvolution of fundamental and harmonic ultrasound images | |
| CN112263278B (zh) | 基于随机孔径的心脏超声杂波抑制方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |