CN104217401B - 超声成像方法及装置 - Google Patents
超声成像方法及装置 Download PDFInfo
- Publication number
- CN104217401B CN104217401B CN201310216537.7A CN201310216537A CN104217401B CN 104217401 B CN104217401 B CN 104217401B CN 201310216537 A CN201310216537 A CN 201310216537A CN 104217401 B CN104217401 B CN 104217401B
- Authority
- CN
- China
- Prior art keywords
- echo
- vector
- ulga
- image
- uhss
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 35
- 239000013598 vector Substances 0.000 claims abstract description 128
- 230000011218 segmentation Effects 0.000 claims abstract description 42
- 230000003044 adaptive effect Effects 0.000 claims abstract description 10
- 238000009826 distribution Methods 0.000 claims description 16
- 238000012937 correction Methods 0.000 claims description 12
- 238000005192 partition Methods 0.000 claims description 11
- 238000000605 extraction Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000011084 recovery Methods 0.000 claims description 8
- 230000005764 inhibitory process Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 15
- 238000012545 processing Methods 0.000 abstract description 7
- 150000001875 compounds Chemical class 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 9
- 238000013507 mapping Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005057 finger movement Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 125000001475 halogen functional group Chemical group 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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/52046—Techniques for image enhancement involving transmitter or receiver
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- 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/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
-
- 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
-
- 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/30196—Human being; Person
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Networks & Wireless Communication (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本申请公开了一种超声成像方法及装置,其基于对回声信号进行隶属函数的模糊分段,将输入图像按回声强度依次划分成若干模糊分段,据此将图像划分成若干层次,进而对若干层次内的隶属度取值进行空间邻域操作,从而得到空间模糊分段数组;另一方面通过基于图像内容的自适应计算或基于用户的输入而确定图像的操纵量向量,最后利用该操纵量向量对各个层次的空间模糊分段数组进行灵活的加权和拉伸处理,并使用得到的回声增量和亮度补偿对输入图像进行调整和复合,从而得到最终的输出图像。本实施例的方法及装置能够实现局部增益调节、高亮回声抑制、低回声恢复、和亮度均衡等功能,从而达到丰富图像信号细节分辨率等目的。
Description
技术领域
本申请涉及超声成像技术领域,尤其涉及一种超声成像方法及装置。
背景技术
随着超声成像技术发展的日新月异,超声成像设备采集的信号动态范围越来越大,高回声和低回声信号越来越丰富,而超声信号由于受到衰减、非均匀声束和探头指向性等的影响,信号和噪声的形态在整幅图像上表现并不统一,信号的动态范围在不同位置、不同时间各不相同。因此,为了将这些高回声和低回声信号同时显示在一起,并且显示更多的细节信号,相关领域中已经提出多种技术方案。
在图像处理领域,Chiu等人提出关于空间非均匀调整的方法(K Chiu,M Herf,PShirley;“用于高对比度图像的空间非均匀尺度函数(Spatially non-uniform scalingfunctions for high contrast images)”;《图形界面论文集》93年刊第245-253页),该方法利用图像本身的高斯平滑作为局部增益,使得不同动态范围的信号得以较好显示。该方案的不足之处在于,处理后的图像在强或弱的边缘周围易出现光晕(halo)现象。
在超声波成像领域,现有技术中还存在多种自适应调整方法,其主要思想都是利用噪声信号确定信号动态范围,再依据确定的该动态范围对超声信号进行映射。这些方案一个不足之处是噪声信号的确定难度较大,对于回波比较弱的超声信号,很难简单地将有用的组织信息和噪声区分开来,如果区分效果不好,就不能对弱的信号或噪声进行有效的控制;其另外一个不足在于这种简单的映射会带来不自然的图像局部对比度,对实际动态范围较小的区域进行过分拉伸,会给组织带来过多的黑洞,或将管腔内的斑点噪声凸显出来。
现有技术中还存在一种区分组织和噪声的对比度的方案,但该方案没有对图像进行伸缩操作,因此无法对弱信号进行显著增强;此外,这项技术还需要利用一个对比度分辨率较好的图像作为基础,从而对图像处理提出了额外要求,增大了处理难度。
综上所述,有必要提出一种应用于超声成像的更好的信号动态范围调整方案,从而在实现显示更多细节信号的同时避免上述相关技术方案存在的各种不足。
发明内容
针对上述问题,本申请的实施例旨在提供一种超声成像方法及装置。
为实现上述目的,本申请的实施例提供了一种超声成像方法,包括以下步骤:
以输入图像特定点处的回声强度作为函数输入,按照根据回声强度强弱分成的p个区间对所述函数输入进行模糊隶属度的划分,得到该特定点处的空间模糊分段数组V,其中V=[V1,V2,…,Vp],V1、V2、…Vp的值均介于0与1之间,且V1+V2+…+Vp≈1;
根据对所述输入图像的内容进行自适应计算得到的结果或者根据用户的输入来确定针对该特定点处的操纵量向量U,所述操纵量向量U为选自以下向量的任意组合之和:局部增益调节向量ULGA、高亮回声抑制向量UHSS、低回声恢复向量ULER、和亮度均衡向量UBE;其中,ULGA=[ULGA1,ULGA2,…,ULGAp],UHSS=[UHSS1,UHSS2,…,UHSSp],ULER=[ULER1,ULER2,…,ULERp],UBE=[UBE1,UBE2,…,UBEp];以及
根据所述空间模糊分段数组V和操纵量向量U求内积得到该特定点处的回声增量ΔI,并利用所述回声增量ΔI调节该特定点处的回声强度,并将经过调节的回声强度作为输出图像在对应点处的回声强度,其中ΔI=UT·V,UT表示U的转置矩阵。
本申请的实施例还提出一种超声成像装置,包括:
模糊分段确定单元,用于以输入图像特定点处的回声强度作为函数输入,按照根据回声强度强弱分成的p个区间对所述函数输入进行模糊隶属度的划分,得到该特定点处的空间模糊分段数组V,其中V=[V1,V2,…,Vp],V1、V2、…Vp的值均介于0与1之间,且V1+V2+…+Vp≈1;
操纵量确定单元,用于根据对所述输入图像的内容进行自适应计算得到的结果或者根据用户的输入来确定针对该特定点处的操纵量向量U,所述操纵量向量U为选自以下向量的任意组合之和:局部增益调节向量ULGA、高亮回声抑制向量UHSS、低回声恢复向量ULER、和亮度均衡向量UBE;其中,ULGA=[ULGA1,ULGA2,…,ULGAp],UHSS=[UHSS1,UHSS2,…,UHSSp],ULER=[ULER1,ULER2,…,ULERp],UBE=[UBE1,UBE2,…,UBEp];以及
图像输出单元,用于根据所述空间模糊分段数组V和操纵量向量U求内积得到该特定点处的回声增量ΔI,并利用所述回声增量ΔI调节该特定点处的回声强度,并将经过调节的回声强度作为输出图像在对应点处的回声强度,其中ΔI=UT·V,UT表示U的转置矩阵。
由上述技术方案可知,本申请的实施例基于回声强度的模糊分段,通过对各个分段使用不同的控制参数,由此实现局部增益调节、高亮回声抑制、低回声恢复、和亮度均衡等功能,能够达到丰富图像信号细节分辨率等目的。
附图说明
图1为本申请超声成像方法的实施例流程图;
图2为图1所示实施例中子块划分的示意图;
图3为图1所示实施例中使用的三角形隶属函数示意图;
图4为作为图3所示隶属函数的示意调整对象的斜坡函数示意图;
图5-图9分别为利用不同的参数组合对图4的斜坡函数信号进行调整的效果示意图;
图10-图12分别为图1所示实施例可使用的其他隶属函数示意图;
图13-图15分别为图1所述实施例中使用的由柱状图表示的回声信号分布特性示意图;
图16为本申请超声成像装置的实施例框图。
具体实施方式
下面将详细描述本申请的具体实施例。应当注意,这里描述的实施例只用于举例说明,并不用于限制本申请。
图1为本申请超声成像方法的实施例流程图,如图所示,本实施例的超声成像方法包括以下步骤S101-S105。
S101、对当前的超声图像进行子块划分;
超声图像除了受到探头、前端电路、信号处理的影响外,还受到待检查病人条件的影响,其中,不同病人的衰减程度、杂波、强反射界面、声窗条件等各不相同。因此,对图像进行检测和分析,实时地改变图像效果以适应不同的条件是一个重要的问题。在本步骤中,通过对图像进行子块划分,有助于掌握图像的局部信息,便于后续对图像内容进行调整。子块划分可按照预设的规则来进行,例如保持子块是矩形、圆形、或其他形状,又例如保持相邻子块之间有一定比例的重叠等等。在一个实施例中,如图2所示,超声图像S被划分成多个圆形子块S(1)、S(2)等等。在另一个实施例中,还可将大小为M×N的矩形超声图像划分成多个大小为m×n的矩形子块,这里,矩形超声图像的大小M×N和矩形子块的大小m×n的单位可以是实际物理长度单位,比如毫米等等,也可以是像素单位。
另外,本步骤中作为子块划分基础的超声图像可以是经过包络提取之后的图像,也即,本实施例的超声成像方法可以是在超声图像经过包络提取之后的某一环节进行,因为经过包络提取的信号可较好地体现回声的幅度,便于对图像进行统计分析。但本申请并不仅限于此,在包络提取之前进行本实施例的超声成像方法也是可行的。
S102、以输入图像特定点处的回声强度作为函数输入,按照根据回声强度强弱分成的p个区间对所述函数输入进行模糊隶属度的划分;
在经过步骤S101的子块划分之后,后续步骤的流程均以划分后的子块图像作为输入图像进行处理。
在本步骤中,为了更好地操纵输入图像的局部灰阶特性,首先可以定义p个强弱不同的回声强度区间。本发明的实施例中,定义的回声强度区间的数量p的值不做特别的限制。通常,p越大,即回声强度区间的数量越多,则可以实现更为精细的调整,但同时也会增加计算量。因此,p的大小可以根据实际情况的需要灵活选择。例如,在一个实施例中,可采取p为3-20左右的区间数进行模糊分段处理。为简便起见,在本文后续实施例中,均以超低回声、低回声、中等回声、高回声、和超高回声这5个区间(即p=5)为例进行详细说明,但本领域技术人员当可了解本申请并不仅限于此。
接续,用隶属函数将上述5个回声强度区间表示出来。下面以对称三角形隶属函数为例,如图3所示,分别用图中从左至右5个分段的三角形函数表示超低回声、低回声、中等回声、高回声、和超高回声这5个区间,图中横坐标表示回声强度,纵坐标表示隶属度,由图示可以看出,这5个函数的隶属度最高值对应的回声强度分别位于5db、25db、45db、65db、和85db处,单个分类的跨度为40db,并且在相邻区间之间具有20db的重叠。
按照各个隶属函数对输入图像进行模糊分段,确定超低回声、低回声、中等回声、高回声、和超高回声这5个区间,分别用VL(X)、L(X)、M(X)、H(X)、和VH(X)表示。这样,以输入图像每个特定点处的回声强度I作为函数输入(即X=I),则通过函数计算或查表便可输出一个长度为5的模糊分段数组:[VL(I),L(I),M(I),H(I),VH(I)],根据隶属函数的特性可知,VL(I)、L(I)、M(I)、H(I)、和VH(I)的取值均介于0到1之间,且它们的和为1。参见图3,以回声强度为5db为例,对应的输出为[1,0,0,0,0],以回声强度为10db为例,对应的输出大约为[0.5,0.5,0,0,0]。
经过本步骤的模糊隶属度划分,可得到强度不同的5个模糊分段,后续的流程即通过对这5个模糊分段施加不同的增益来实现图像细节的控制,在此之前,这里先将以斜坡函数(如图4所示)为例简要描述后续的分段调整方法对图像灰阶的影响。
假设对函数输入X处的各个模糊分段施加不同的增益分别为:GainVL、GainL、GainM、GainH、和GainVH,则X处的信号改变量可用如下计算式表示:
ΔX=GainVL·VL(X)+GainL·L(X)+GainM·M(X)+
GainH·H(X)+GainVH·VH(X)+ΔGainComp
其中,ΔGainComp为亮度补偿量,其大小与回声调整的分段有关,当调整中低回声的幅度时,对图像的整体亮度则不会有较大影响;相反,当调整中高回声时,则需要较大的幅度来调整图像的亮度。
首先参照图5和图6所示,其是通过GainM来调整信号,从而改变最终的图像灰阶。图5中设置GainM>0,其作用是增加中低回声的局部对比度,同时降低中高回声的局部对比度,GainM越大这种效果就越强。同理,图6中设置GainM<0,其作用是增加中高回声的对比度,同时降低中低回声的对比度。从结果上来看,通过改变GainM,可以在保持回声信号单调性的同时,实现图像高低回声层次的调整。
进一步,图7中设置GainVL、GainL、GainH和GainVH均大于0,而GainM=0,其效果类似于图6的效果,不同之处在于图像亮度有所增加,起到的灰阶调整效果则与图6完全一致的。
接续,为了提升中等偏高回声的对比,同时避免导致部分过强的高亮回声饱和,可以设置:GainM<0、GainH>0和GainVH<0,其效果如图8所示。
最后,对于动态范围过大的图像来说,将有效的图像信息有效地映射到显示设备是一个非常重要的问题。例如当信号动态范围超过显示动态范围时,高亮回声可能已出现大片饱和,或者低回声信息出现丢失,此时可分别设置GainVL>0和GainVH<0,由此将丢失的信号动态范围调整到人眼能够识别的灰阶范围内,效果如图9所示。
需要说明的是,上面的模糊隶属度划分是以对称的三角形隶属函数为例进行说明,但本申请显然不限于此,本领域技术人员可以理解还可采用其他形式的隶属函数。
例如,在一个实施例中,可以应用如图10所示三个顶点依次对应于横坐标a、b、c的一般三角形函数形式,其中f(x)=max(min((x-a)/(b-a),(c-x)/(c-b)),0)。其中a为图10中A点的横坐标,b为图10中B点的横坐标,c为图10中C点的横坐标。这里,A、B和C分别为图10中的三角形函数的顶点。
在另一个实施例中,为了改善信号梯度过渡的连续性,可使用光滑的隶属函数,例如图11所示的分布参数为、中心在C处高斯函数,其中
这里,c为图11中的高斯函数的顶点(中心点)C的横坐标。
在再一个实施例中,为了信号控制的灵活性,还可以将相邻两个或者多个隶属函数合并成一个,从而形成一个“平顶”的隶属函数,如图12所示的“钟形”函数,其中
这里,c为图12中的钟形隶属函数的左右对称中心C的横坐标,a和b是该钟形隶属函数的控制参数,用于控制该钟形隶属函数的形状和跨度。
S103、对模糊分段数组进行空间邻域的平滑操作;
如上所述,步骤S102是以输入图像特定点处的回声强度作为函数输入进行模糊隶属度的划分,以得到长度为p的模糊分段数组输出。然而,这种划分仅仅是根据回声强度进行分类,而并未考虑到图像的空间连续性。具体而言,考虑某一输入图像,如果存在离散的黑孔,当组织回声和这些黑孔的回声分属不同的模糊集合时,而两个集合的操纵量又相反时,这些图像上的黑孔将更加明显。
有鉴于此,为了改善最终图像的连续性,抑制组织上的黑孔和管腔内的噪点等局部突变,还可以经过本步骤对由输入图像所有点的模糊分段数组[VL(I),L(I),M(I),H(I),VH(I)]组成的信号矩阵进行空间邻域的平滑操作,使得输入图像中在空间上相邻近的点被划分到一个区间的隶属度较大,由此保留一定程度的“空间结构”,使图像能够在灰阶调整中保持空间上的连续性。在一个实施例中,上述空间邻域的平滑操作可采用线性或非线性的平滑来进行,例如形态学开/闭运算、非线性黑孔填充、低通滤波、旋转掩模、和中值滤波等等算法。经过上述平滑操作的信号矩阵能够保留局部灰阶层次的空间连续结构,而并不具有原图像的细节分辨率和对比度分辨率,由此对于输入图像的每个特定点处则均可得到一个对应的空间模糊分段数组V=[V1,V2,…,Vp]。注意到V1、V2、…Vp的取值仍介于0和1之间,表示当前回声图像相对于各个区间的隶属程度,而不再反映回声强度的大小,并且它们的和仍然近似等于1。
本步骤的空间邻域平滑操作是在模糊隶属度划分之后进行,但本申请并不仅限于此,在一个实施例中,为了改善图像的空间连续性,还可以先对输入图像进行平滑操作,然后再通过模糊隶属函数计算获得长度为5的空间模糊分段数组V=[V1,V2,V3,V4,V5]。
S104、确定针对该特定点处的操纵量向量U;
在经过上述步骤S101-S103得到输入图像特定点处对应的空间模糊分段数组V的基础上,接续便可以确定针对该特定点处的操纵量向量U。在一个实施例中,操纵量向量U是局部增益调节(Local Gain Adjustment)向量ULGA、高亮回声抑制(High SaturationSuppression)向量UHSS、低回声恢复(Low Echo Recovery)向量ULER、和亮度均衡(Brightness Equalization)向量UBE的任意组合之和,下面将分别予以详细阐述。
由于超声成像***的内部增益以及不同病人的情况等原因,不同子块输入图像的信号与噪声分布情况各不相同,这就有必要进行局部增益的校正,以改善整幅图像的一致性。局部增益参数的计算需要基于当前子块及邻近子块的回声强度分布信息,通过邻近区域的图像的直方图进行分析,由此计算出局部增益校正量。
有些图像由于信噪比的差异,仅通过局部增益无法改善图像的均匀性,这时就有必要配合进行亮度均衡的调整。与局部增益调整的方式不同,亮度均衡侧重于不同信噪比图像亮度的一致化。具体而言,信噪比较低的区域,如果仅通过局部增益来提升亮度,则增加信号的同时也会增加噪声,带来负面影响。因此,有必要根据当前信号分布情况,通过亮度均衡来调整各个回声模糊分段的比例,从而达到调整灰阶层次的目的,最终改善整体图像的一致性。
超声波遇到强反射界面,会产生较强的回波信号,由于显示动态范围有限,很容易产生较强的高亮回声饱和,体现在超声图像上就是出现一片没有层次的高亮灰阶,通过分析显示动态范围和回声强度,可以计算出现高亮回声的相对区间,通过调整映射曲线的形状,可以恢复部分高亮回声的层次。同理,低回声信号恢复是将容易向下溢出的超声信号调整到用户的显示动态范围之内,从而增加有效信号的信息量。高亮回声抑制和低回声信号恢复结合使用,配合用户动态范围,可以优化图像显示内容,改善图像质量。
为了实现上述的局部增益调节、高亮回声抑制、低回声恢复和亮度均衡,可以基于前述步骤的空间模糊分段来实现。具体而言,对于这4个功能,针对前述5个模糊分段可以用4组长度为5的控制向量来对应。在一个实施例中,这些控制向量可以根据图像内容进行自适应计算而得到,或者也可以通过与用户的交互而基于用户的输入得到,下面首先将对前一种情况(即图像内容的自适应计算)加以详细阐述。
(1)局部增益调节
局部增益调节向量的确定原则比较简单,只需要设置向量中针对超低回声、低回声、中等回声、高回声、和超高回声这5个区间的所有参数值均相等,即,ULGA1=ULGA2=…=ULGAp=ΔGain。其中ΔGain为当前点需要调整的局部增益,在一个实施例中,ΔGain可根据例如美国专利申请公开第US20120260736号、中国专利申请第200910109156中揭示的方案等现有技术来确定。
相对于局部增益调节而言,亮度均衡的参数设置灵活度更大。概略而言,对于信噪比较高的图像区域,针对中、高回声的参数应低于针对低回声的参数,反之,对于信噪比较低的图像,应根据情况适当调整中、低回声的参数增加中低信号的对比度。低回声信号恢复是将尽可能多的低回声信号映射到显示动态范围之内,比如前面结合图9提到的设置GainVL>0,并且,根据信号分布的特征,还可以适度设置GainL和GainM,从而得到更加适度的中低回声分布。同理,对于局部的高亮回声信号,可设置GainVH<0,同时设置GainH和GainM来保留中高回声的层次,抑制高亮回声抑制问题。进一步,对于亮度均衡、低回声信号恢复和高亮回声抑制来说,向量中各个参数可能各不相等,下面将对各向量的确定方法和原则分别加以说明。
(2)高亮回声抑制
按照上述步骤S102中结合图5-图9描述的调整示例可知,可以通过改变各操纵参数的符号和大小来控制图像的改变量。因此,在本步骤中,首先可以通过设置符号和大小各不相同的参数组合来预设不同级别的高亮回声抑制向量。例如,在一个实施例中,可以定义非常强、较强、中等、较小、不变这5个级别的高亮回声抑制向量。其中,对应每个级别的操纵参数的单位可以是db,大小可以根据临床图像标准的需要灵活设置。例如“非常强”这个级别的设置可以根据按照临床的要求能产生临床图像显著变化的要求,并且其余各个档位应有明显的依次变化,等等。
此外,实际操作参数也可以是介于两个档位之间的模糊值,比如介于“非常强”与“较强”之间的值,或者“不变”和“较小”之间的值,均根据实际临床图像的需要而定。
例如,一个实施例中,各个高回声抑制向量可以分别如下所示。
非常强:UHSS=[0,0,0,-5,-10]
较强:UHSS=[0,0,0,-3,-7]
中等:UHSS=[0,0,0,-2,-3]
较小:UHSS=[0,0,0,-1,-1]
不变:UHSS=[0,0,0,0,0]
接续,通过对输入图像进行分析获得回声信号的分布特性,并根据该分布特性确定高亮回声饱和的级别。具体可参见图13所示的示例,其柱状图即反映了当前输入图像的回声信号分布,在该示例中,假设超声设备的用户显示动态范围是50db(10db-60db),换言之,强度为10db的回声信号在屏幕上显示为灰阶0,强度为60db的信号显示为灰阶255,而小于10db或大于60db的信号都不会显示在屏幕上,分别对应为弱回声信号损失和高亮回声饱和。从图13可以看出,该输入图像中的大部分回声信号都超出了60db,这会引起较为严重的高亮回声饱和,由此从“非常严重、较严重、中等、轻微、无”这5个预设级别中选择当前输入图像的高亮回声饱和级别确定为“较严重”,或者可以按照模糊推理的原则(比如,可以使用如文献“Mamdani,E.H.,“Applications of fuzzy logic to approximate reasoningusing linguistic synthesis:,IEEE transactions on computers,vol.26,no.12,pp.1182-1191,1977.”中揭露的经典的Mamdani推理法)得到图像高亮饱和的级别介于“非常严重”和“较为严重”之间。需要说明的是,上述示例中是使用柱状图统计的方法来描述回声信号的分布特性,但本申请并不仅限于此,在其他实施例中,此处也可通过上述步骤S102中提及的隶属函数来描述回声分布特性,虽然隶属函数和柱状图的计算方法不同,但二者从模糊理论和统计学上的内涵是相通的,二者都能从图中直观得到高亮回声饱和的级别,是以此处对于隶属函数的表示方法不再赘述。
这之后,对应于确定的“较严重”的高亮回声饱和,便需要较强的高亮回声抑制功能,因此,对应即从上述预设级别的高亮回声抑制向量中选择“较强”一级的高亮回声抑制向量[0,0,0,-3,-7]作为最终确定的高亮回声抑制向量UHSS。
(3)低回声恢复
与上述高亮回声抑制向量的确定类似,在确定低回声恢复向量时,首先可以预设不同级别的低回声恢复向量。例如,在一个实施例中,可以定义非常强、较强、中等、较小、不变这5个级别的低回声恢复向量,其结果可以分别如下所示。
非常强:ULER=[10,5,0,0,0]
较强:ULER=[7,3,0,0,0]
中等:ULER=[3,2,0,0,0]
较小:ULER=[2,1,0,0,0]
不变:ULER=[0,0,0,0,0]
接续,通过对输入图像进行分析获得的回声信号分布特性来确定低回声丢失的级别。具体可参见图14所示的示例,从其柱状图中可以看出,该输入图像中的一部分回声信号都在10db以下,表面需要进行一定程度的低回声信号恢复,由此从“非常严重、较严重、中等、轻微、无”这5个预设级别中选择将当前输入图像的低回声丢失级别确定为“中等”。
这之后,对应于确定的“中等”的低回声丢失,便需要中等的低回声恢复功能,因此,对应即从上述预设级别的低回声恢复向量中选择“中等”一级的低回声恢复向量[3,2,0,0,0]作为最终确定的低回声恢复向量ULER。
(4)亮度均衡
与上述高亮回声抑制向量的确定类似,在亮度均衡向量时,首先可以预设不同级别的亮度均衡向量。例如,在一个实施例中,可以定义非常强、较强、中等、较小、不变这5个级别的亮度均衡向量,其结果可以分别如下所示。
非常强:UBE=[0,5,10,6,2]
较强:UBE=[0,3,7,3,1]
中等:UBE=[0,2,5,2,1]
较小:UBE=[0,1,2,1,0]
不变:UBE=[0,0,0,0,0]
接续,通过对输入图像进行分析获得的回声信号分布特性来确定亮度明暗的级别。具体可参见图15所示的示例,从其柱状图中可以看出,该输入图像的信号较弱,需要增加图像亮度,由此从“非常暗、比较暗、有些暗、轻微暗、正常”这5个预设级别中选择将当前输入图像的亮度明暗级别确定为“比较暗”。
这之后,对应于确定的“比较暗”的亮度明暗级别,便需要较强的亮度均衡功能,因此,对应即从上述预设级别的亮度均衡向量中选择“较强”一级的亮度均衡向量[0,3,7,3,1]作为最终确定的亮度均衡向量UBE。
由前文所述,本发明的实施例中,各个控制向量(例如,高亮回声抑制向量UHSS、低回声恢复向量ULER和亮度均衡向量UBE)可以是基于对输入图像进行分析获得。此时,本发明的实施例中,可以是对输入图像的一次分析中同时获得多个控制向量(例如,高亮回声抑制向量UHSS、低回声恢复向量ULER和亮度均衡向量UBE可以在对输入图像的分析中一次性同时获得),也可以是分别在对输入图像的不同的分析步骤中获得不同的控制向量。
在如上所述分别确定了局部增益调节向量、高亮回声抑制向量、低回声恢复向量、和亮度均衡向量之后,可对这些向量进行求和运算来得到最终的操纵量向量U=ULGA+UHSS+ULER+UBE。但本申请并不仅限于此,在一个实施例中,根据实际情况需要,还可以从局部增益调节向量、高亮回声抑制向量、低回声恢复向量、和亮度均衡向量中选取任意组合进行求和运算来得到操纵量向量U。
上述实施例中详细描述了根据图像内容进行自适应计算来确定操纵量向量U的示例,然而本申请并不仅限于此。在一个实施例中,还可通过设置用户输入装置,结合与用户的交互进行计算,根据预设规则对用户的输入进行解释和翻译来分别得到确定操纵量向量U。以触摸屏作为用户输入装置为例,用户(例如医生)可通过手指在触摸屏上产生输入,从而可对触摸点坐标及其随时间变化的规律、手指移动的开始点和终止点、移动速度和加速度、以及移动路径等等进行计算和分析,并根据计算和分析结果参照预设规则来得到局部增益调节向量、高亮回声抑制向量、低回声恢复向量、和亮度均衡向量等各个参数。
S105、根据空间模糊分段数组V和操纵量向量U求内积得到特定点处的回声增量,进而得到输出图像。
在经过步骤S103得到空间模糊分段数组V以及经过步骤S104得到操纵量向量U之后,本步骤即可进行图像的回声增量计算。具体而言,在所述特定点处,回声增量可表示为ΔI=UT·V。其中,“T”表示矩阵的转置符号,“·”表示求内积的运算符号。如上所述,空间模糊分段数组V和操纵量向量U的长度均为n,因此通过求内积运算可消掉长度为n的这一维,得到特定点处的回声增量ΔI。
在一个示例中,假设步骤S104中确定的局部增益调节向量、高亮回声抑制向量、低回声恢复向量、和亮度均衡向量分别如下:
ULGA=[3,3,3,3,3],其功能是提升局部图像整体亮度;
UHSS=[0,0,0,0,-2],其功能是抑制超高回声的高亮饱和;
ULER=[2,2,1,0,0],其功能是抬高超低、低回声信号;
UBE=[0,-1,2,5,0],其功能是增加中高回声层次。
则通过求和得到最终的操纵量向量U如下:
U=ULGA+UHSS+ULER+UBE=[5,4,6,8,1]
进而通过求内积得到回声增量ΔI如下:
ΔI=UT·V=[5,4,6,8,1]T·[V1,V2,V3,V4,V5]
=5*V1+4*V2+6*V3+8*V4+1*V5
最后,利用求得的回声增量ΔI对该特定点处的回声强度进行调节,并用经过调节后的回声强度作为输出图像特定点处的回声强度,由此可得到输出图像。
结合步骤S102中提及的信号改变量计算式可知,在利用回声增量ΔI进行回声强度的调节后,还可以对输入图像进行亮度校正,由于局部增益调节向量、高亮回声抑制向量、低回声恢复向量、和亮度均衡向量各参数均已确定,因此校正量ΔGainComp也是可计算的。
本发明的实施例中,ΔGainComp可以根据使校正后的图像的亮度与动态范围调节前的图像的亮度基本一致的要求,根据前述已经获得的局部增益调节向量、高亮回声抑制向量、低回声恢复向量、亮度均衡向量、最终的操纵量向量和/或前述获得的回声增量确定,具体的确定方法可以根据实际情况的需要而定,只要使校正后的图像的亮度与动态范围调节前的图像的亮度基本一致即可。
在利用得到的ΔGainComp进行亮度校正后,可进一步得到亮度均匀的输出图像。
与上述方法实施例对应,本申请的实施例还提供一种超声成像装置。图16为本申请超声成像装置的实施例框图,如图所示,本实施例的超声成像装置包括图像子块划分单元21、模糊分段确定单元22、空间邻域平滑单元23、操纵量确定单元24、图像输出单元25、和亮度校正单元26。其中:
图像子块划分单元21用于对经过包络提取后得到的图像进行子块划分,并以子块划分后得到的各个子块图像作为输入图像。在一个实施例中,图像子块划分单元21也可以在包络提取之前进行子块划分,并直接以划分后的子块图像作为输入图像。
模糊分段确定单元22用于以输入图像特定点处的回声强度作为函数输入,按照根据回声强度强弱分成的p个区间对函数输入进行模糊隶属度的划分,得到该特定点处的空间模糊分段数组V,其中V=[V1,V2,…,Vp],V1、V2、…Vp的值均介于0与1之间,且V1+V2+…+Vp≈1。
空间邻域平滑单元23用于在模糊分段确定单元22进行模糊隶属度的划分之后进行空间邻域的平滑操作,使得输入图像中在空间上相邻近的点被划分到一个区间的隶属度较大,由此得到最终的空间模糊分段数组V。
操纵量确定单元24用于根据对输入图像的内容进行自适应计算得到的结果或者根据用户的输入来确定针对该特定点处的操纵量向量U,该操纵量向量U为选自以下向量的任意组合之和:局部增益调节向量ULGA、高亮回声抑制向量UHSS、低回声恢复向量ULER、和亮度均衡向量UBE;其中,ULGA=[ULGA1,ULGA2,…,ULGAp],UHSS=[UHSS1,UHSS2,…,UHSSp],ULER=[ULER1,ULER2,…,ULERp],UBE=[UBE1,UBE2,…,UBEp]。
图像输出单元25用于根据模糊分段确定单元22得到的空间模糊分段数组V和操纵量确定单元24得到的操纵量向量U求内积来得到该特定点处的回声增量ΔI,并利用该回声增量ΔI调节所述特定点处的回声强度,且将经过调节的回声强度作为输出图像在对应点处的回声强度,其中ΔI=UT·V,UT表示U的转置矩阵。
亮度校正单元26则用于在图像输出单元25通过求内积得到该特定点处的回声增量ΔI之后,根据该回声增量ΔI确定亮度补偿量ΔGainComp,并利用该亮度补偿量ΔGainComp对图像进行亮度校正后得到最终的输出图像。
进一步,在一个实施例中,上述的操纵量确定单元24包括选自局部增益计算子单元241、高亮回声抑制计算子单元242、低回声恢复计算子单元243、和亮度均衡计算子单元244的任意组合。其中:
局部增益计算子单元241用于确定特定点处需要调整的局部增益ΔGain,并使得ULGA1=ULGA2=…=ULGAp=ΔGain,由此得到所述局部增益调节向量ULGA。
高亮回声抑制计算子单元242、低回声恢复计算子单元243、和亮度均衡计算子单元244则分别用于对输入图像进行分析以得到回声信号的分布特性,并根据该分布特性来分别确定高亮回声饱和的级别、低回声丢失的级别、和亮度明暗的级别,之后再根据确定的各级别分别查找得到对应级别的预先分级别设置的高亮回声抑制向量UHSS、低回声恢复向量ULER、和亮度均衡向量UBE。
综上所述,本实施例的超声成像方法及装置基于对回声信号进行隶属函数的模糊分段,例如可将噪声和弱回声信号划分为超低回声信号,将向上饱和的信号(在用户显示动态范围之外)划分为超高回声信号,然后依次划分成若干模糊分段,据此将图像划分成若干属于不同模糊分段的层次,进而对若干层次内的隶属度取值进行空间邻域操作,从而得到空间模糊分段数组。另一方面通过基于图像内容的自适应计算或基于用户的输入而确定图像的操纵量向量,最后利用该操纵量向量对各个层次的空间模糊分段数组进行灵活的加权和拉伸处理(求内积),并使用得到的回声增量和亮度补偿对输入图像进行调整和复合,从而得到最终的输出图像。本实施例的方法及装置避免了处理噪声本底和高对比度分辨率的图像,用灵活的空间模糊加权处理避免了直接映射带来的副作用。
由上述技术方案可知,本申请的实施例基于回声强度模糊分段的灰阶调整方法,通过对各个分段使用不同的控制参数,实现局部增益调节、高亮回声抑制、低回声恢复、和亮度均衡等功能,从而达到丰富图像信号细节分辨率等目的。
按照本申请实施例的超声成像方法及装置,可以通过硬件、软件、固件、或者其组合实现在超声成像***中,从而使得超声成像***可以采用基于本申请实施例的超声成像方法,或者包括基于本申请实施例的超声成像装置。按照本申请的上述教导,这种实现对于本领域普通技术人员而言是显而易见的,在此不再赘述。
虽然已参照几个典型实施例描述了本申请,但应当理解,所用的术语是说明和示例性、而非限制性的术语。由于本申请能够以多种形式具体实施而不脱离申请的精神或实质,所以应当理解,上述实施例不限于任何前述的细节,而应在随附权利要求所限定的精神和范围内广泛地解释,因此落入权利要求或其等效范围内的全部变化和改型都应为随附权利要求所涵盖。
Claims (10)
1.一种超声成像方法,包括以下步骤:
以输入图像特定点处的回声强度作为函数输入,按照根据回声强度强弱分成的p个区间对所述函数输入进行模糊隶属度的划分,得到该特定点处的空间模糊分段数组V,其中V=[V1,V2,…,Vp],V1、V2、…Vp的值均介于0与1之间,且V1+V2+…+Vp≈1;
根据对所述输入图像的内容进行自适应计算得到的结果或者根据用户的输入来确定针对该特定点处的操纵量向量U,所述操纵量向量U为选自以下向量的任意组合之和:局部增益调节向量ULGA、高亮回声抑制向量UHSS、低回声恢复向量ULER、和亮度均衡向量UBE;其中,ULGA=[ULGA1,ULGA2,…,ULGAp],UHSS=[UHSS1,UHSS2,…,UHSSp],ULER=[ULER1,ULER2,…,ULERp],UBE=[UBE1,UBE2,…,UBEp];以及
根据所述空间模糊分段数组V和操纵量向量U求内积得到该特定点处的回声增量ΔI,并利用所述回声增量ΔI调节该特定点处的回声强度,并将经过调节的回声强度作为输出图像在对应点处的回声强度,其中ΔI=UT·V,UT表示U的转置矩阵。
2.如权利要求1所述的超声成像方法,其中在进行所述模糊隶属度的划分之后还包括:进行空间邻域的平滑操作,使得所述输入图像中在空间上相邻近的点被划分到一个区间的隶属度较大,由此得到所述空间模糊分段数组V。
3.如权利要求1所述的超声成像方法,其中在进行所述模糊隶属度的划分之前还包括:对经过包络提取后得到的图像进行子块划分,以子块划分后得到的各个子块图像作为所述输入图像。
4.如权利要求1所述的超声成像方法,其中在通过所述求内积得到该特定点处的回声增量ΔI之后还包括根据所述局部增益调节向量ULGA、所述高亮回声抑制向量UHSS、所述低回声恢复向量ULER、所述亮度均衡向量UBE和/或所述回声增量ΔI确定亮度补偿量ΔGainComp,之后则利用该亮度补偿量ΔGainComp对所述输出图像进行亮度校正。
5.如权利要求1至4任一项所述的超声成像方法,其中所述对输入图像的内容进行自适应计算包括以下步骤的任意组合:
确定所述特定点处需要调整的局部增益ΔGain,使ULGA1=ULGA2=…=ULGAp=ΔGain,由此得到所述局部增益调节向量ULGA;以及
对所述输入图像进行分析以得到回声信号的分布特性,根据该分布特性确定高亮回声饱和的级别、和/或低回声丢失的级别、和/或亮度明暗的级别,之后根据确定的各级别分别查找得到对应级别的预先分级别设置的高亮回声抑制向量UHSS、和/或低回声恢复向量ULER、和/或亮度均衡向量UBE。
6.一种超声成像装置,包括:
模糊分段确定单元,用于以输入图像特定点处的回声强度作为函数输入,按照根据回声强度强弱分成的p个区间对所述函数输入进行模糊隶属度的划分,得到该特定点处的空间模糊分段数组V,其中V=[V1,V2,…,Vp],V1、V2、…Vp的值均介于0与1之间,且V1+V2+…+Vp≈1;
操纵量确定单元,用于根据对所述输入图像的内容进行自适应计算得到的结果或者根据用户的输入来确定针对该特定点处的操纵量向量U,所述操纵量向量U为选自以下向量的任意组合之和:局部增益调节向量ULGA、高亮回声抑制向量UHSS、低回声恢复向量ULER、和亮度均衡向量UBE;其中,ULGA=[ULGA1,ULGA2,…,ULGAp],UHSS=[UHSS1,UHSS2,…,UHSSp],ULER=[ULER1,ULER2,…,ULERp],UBE=[UBE1,UBE2,…,UBEp];以及
图像输出单元,用于根据所述空间模糊分段数组V和操纵量向量U求内积得到该特定点处的回声增量ΔI,并利用所述回声增量ΔI调节该特定点处的回声强度,并将经过调节的回声强度作为输出图像在对应点处的回声强度,其中ΔI=UT·V,UT表示U的转置矩阵。
7.如权利要求6所述的超声成像装置,其中还包括空间邻域平滑单元,用于在所述模糊分段确定单元进行模糊隶属度的划分之后进行空间邻域的平滑操作,使得所述输入图像中在空间上相邻近的点被划分到一个区间的隶属度较大,由此得到所述空间模糊分段数组V。
8.如权利要求6所述的超声成像装置,其中还包括图像子块划分单元,用于在所述模糊分段确定单元进行模糊隶属度的划分之前,对经过包络提取后得到的图像进行子块划分,并以子块划分后得到的各个子块图像作为所述输入图像。
9.如权利要求6所述的超声成像装置,其中还包括亮度校正单元,用于在所述图像输出单元通过求内积得到该特定点处的回声增量ΔI之后,根据所述局部增益调节向量ULGA、所述高亮回声抑制向量UHSS、所述低回声恢复向量ULER、所述亮度均衡向量UBE和/或所述回声增量ΔI确定亮度补偿量ΔGainComp,并利用该亮度补偿量ΔGainComp对所述输出图像进行亮度校正。
10.如权利要求6-9任一项所述的超声成像装置,其中所述操纵量确定单元包括以下子单元的任意组合:
局部增益计算子单元,用于确定所述特定点处需要调整的局部增益ΔGain,使ULGA1=ULGA2=…=ULGAp=ΔGain,由此得到所述局部增益调节向量ULGA;以及
高亮回声抑制计算子单元、低回声恢复计算子单元、和亮度均衡计算子单元,分别用于对所述输入图像进行分析以得到回声信号的分布特性,并根据该分布特性确定高亮回声饱和的级别、低回声丢失的级别、和亮度明暗的级别,之后根据确定的各级别分别查找得到对应级别的预先分级别设置的高亮回声抑制向量UHSS、低回声恢复向量ULER、和亮度均衡向量UBE。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310216537.7A CN104217401B (zh) | 2013-06-03 | 2013-06-03 | 超声成像方法及装置 |
PCT/CN2014/074326 WO2014194710A1 (zh) | 2013-06-03 | 2014-03-28 | 超声成像方法及装置 |
US14/958,781 US9818185B2 (en) | 2013-06-03 | 2015-12-03 | Ultrasonic imaging method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310216537.7A CN104217401B (zh) | 2013-06-03 | 2013-06-03 | 超声成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104217401A CN104217401A (zh) | 2014-12-17 |
CN104217401B true CN104217401B (zh) | 2018-09-11 |
Family
ID=52007510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310216537.7A Active CN104217401B (zh) | 2013-06-03 | 2013-06-03 | 超声成像方法及装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9818185B2 (zh) |
CN (1) | CN104217401B (zh) |
WO (1) | WO2014194710A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110322413A (zh) * | 2019-07-05 | 2019-10-11 | 深圳开立生物医疗科技股份有限公司 | 超声血流图像的增益调节方法、装置、设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1098616A (zh) * | 1994-05-15 | 1995-02-15 | 重庆医科大学附属第二医院 | 超声定量诊断***及其回声图象定量化方法 |
CN1646934A (zh) * | 2002-04-01 | 2005-07-27 | 皇家飞利浦电子股份有限公司 | 具有自动受控对比度和亮度的超声波诊断成像*** |
CN1718164A (zh) * | 2004-07-07 | 2006-01-11 | 株式会社东芝 | 超声波诊断装置、图像处理装置和图像处理方法 |
CN1798523A (zh) * | 2003-06-03 | 2006-07-05 | 株式会社日立医药 | 超声波诊断装置 |
CN101987023A (zh) * | 2009-07-31 | 2011-03-23 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像增益补偿及图像优化方法及其装置和*** |
CN102130660A (zh) * | 2009-11-13 | 2011-07-20 | 超级科技公司 | 低噪声二进制编码增益放大器和医疗超声成像中时间增益补偿的方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69416584T2 (de) * | 1994-08-25 | 1999-10-07 | Co.Ri.M.Me. Consorzio Per La Ricerca Sulla Microelettronica Nel Mezzogiorno, Catania | Unscharfes Verfahren zur Bildrauschverringerung |
US7689016B2 (en) * | 2005-05-27 | 2010-03-30 | Stoecker & Associates, A Subsidiary Of The Dermatology Center, Llc | Automatic detection of critical dermoscopy features for malignant melanoma diagnosis |
CN100595789C (zh) * | 2006-09-21 | 2010-03-24 | 复旦大学 | 一种基于二维模糊聚类的超声图像分割方法 |
CN101207697B (zh) * | 2006-12-22 | 2011-06-22 | 昆山杰得微电子有限公司 | 一种融合多模糊算子的数字图像模糊增强方法 |
US8478042B2 (en) * | 2009-10-26 | 2013-07-02 | Texas Instruments Incorporated | Method and apparatus for enhancing image or video quality using an exposure aware scene adaptive global brightness contrast |
CN101719277B (zh) * | 2009-12-31 | 2011-11-30 | 华中科技大学 | 一种遗传模糊聚类图像分割方法 |
WO2012154260A2 (en) * | 2011-02-17 | 2012-11-15 | The Johns Hopkins University | Multiparametric non-linear dimension reduction methods and systems related thereto |
US8795179B2 (en) * | 2011-04-12 | 2014-08-05 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Methods, modules, and systems for gain control in B-mode ultrasonic imaging |
US8311973B1 (en) * | 2011-09-24 | 2012-11-13 | Zadeh Lotfi A | Methods and systems for applications for Z-numbers |
-
2013
- 2013-06-03 CN CN201310216537.7A patent/CN104217401B/zh active Active
-
2014
- 2014-03-28 WO PCT/CN2014/074326 patent/WO2014194710A1/zh active Application Filing
-
2015
- 2015-12-03 US US14/958,781 patent/US9818185B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1098616A (zh) * | 1994-05-15 | 1995-02-15 | 重庆医科大学附属第二医院 | 超声定量诊断***及其回声图象定量化方法 |
CN1646934A (zh) * | 2002-04-01 | 2005-07-27 | 皇家飞利浦电子股份有限公司 | 具有自动受控对比度和亮度的超声波诊断成像*** |
CN1798523A (zh) * | 2003-06-03 | 2006-07-05 | 株式会社日立医药 | 超声波诊断装置 |
CN1718164A (zh) * | 2004-07-07 | 2006-01-11 | 株式会社东芝 | 超声波诊断装置、图像处理装置和图像处理方法 |
CN101987023A (zh) * | 2009-07-31 | 2011-03-23 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像增益补偿及图像优化方法及其装置和*** |
CN102130660A (zh) * | 2009-11-13 | 2011-07-20 | 超级科技公司 | 低噪声二进制编码增益放大器和医疗超声成像中时间增益补偿的方法 |
Non-Patent Citations (1)
Title |
---|
基于模糊数模型的医学超声图像分割方法;杨先凤 等;《计算机工程与设计》;20111231;第32卷(第5期);第1742-1744页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2014194710A1 (zh) | 2014-12-11 |
US20160196646A1 (en) | 2016-07-07 |
CN104217401A (zh) | 2014-12-17 |
US9818185B2 (en) | 2017-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Santhi et al. | Adaptive contrast enhancement using modified histogram equalization | |
CN110706826B (zh) | 一种基于视频图像的非接触式的实时多人心率和血压测量方法 | |
Hasikin et al. | Adaptive fuzzy contrast factor enhancement technique for low contrast and nonuniform illumination images | |
WO2020177434A1 (zh) | 图像处理方法及装置、图像设备及存储介质 | |
Hassan1&2 et al. | A new approach for contrast enhancement using sigmoid function | |
Mohan et al. | Particle Swarm Optimization based Contrast Limited enhancement for mammogram images | |
CN105046663A (zh) | 一种模拟人类视觉感知的自适应低照度图像增强方法 | |
CN111999731B (zh) | 一种基于感知生成对抗网络的电磁逆散射成像方法 | |
CN111127476A (zh) | 一种图像处理方法、装置、设备及存储介质 | |
US20120069020A1 (en) | Lighting Control for Occlusion-based Volume Illumination of Medical Data | |
CN112560785A (zh) | 基于人工智能的面部追踪调整多屏亮度的控制方法 | |
CN105266849A (zh) | 实时超声弹性成像方法和*** | |
Verma et al. | Modified sigmoid function based gray scale image contrast enhancement using particle swarm optimization | |
CN104217401B (zh) | 超声成像方法及装置 | |
Al-Ameen | Contrast Enhancement of Digital Images Using an Improved Type-II Fuzzy Set-Based Algorithm. | |
Fan et al. | Cervical lesion image enhancement based on conditional entropy generative adversarial network framework | |
CN103345730A (zh) | 基于侧抑制网络的红外图像处理方法 | |
CN111127344A (zh) | 一种基于bp神经网络的自适应双边滤波的超声图像降噪的方法 | |
Liu et al. | Illumination compensation and feedback of illumination feature in face detection | |
Chang | A flexible contrast enhancement method with visual effects and brightness preservation: Histogram planting | |
JPH0991421A (ja) | 画像処理方法および装置 | |
CN112513673B (zh) | 具有改进的动态范围控件的超声成像*** | |
CN114648463A (zh) | 一种灰度图像动态范围调整方法、装置、设备及介质 | |
Huang et al. | An efficient single image enhancement approach using luminance perception transformation | |
Ramlan et al. | Image Enhancement Methods for Anterior Segment Photographed Images |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20141217 Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd. Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd. Contract record no.: X2022440020009 Denomination of invention: Ultrasound imaging method and device Granted publication date: 20180911 License type: Common License Record date: 20220804 |