CN104318520A - 一种像素局部区域方向检测方法 - Google Patents
一种像素局部区域方向检测方法 Download PDFInfo
- Publication number
- CN104318520A CN104318520A CN201410510938.8A CN201410510938A CN104318520A CN 104318520 A CN104318520 A CN 104318520A CN 201410510938 A CN201410510938 A CN 201410510938A CN 104318520 A CN104318520 A CN 104318520A
- Authority
- CN
- China
- Prior art keywords
- pixel
- filtering
- image
- detection method
- operator
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Image Processing (AREA)
Abstract
本发明提供一种像素局部区域方向检测方法,包括以下步骤:确定方向滤波算子;确定滤波算子的至少一个方向参数;依次采用滤波算子的方向参数对图像进行滤波,得到至少一个滤波值;选择图像滤波值中的最大值对应的方向参数所指示的方向为每个像素的方向。本发明有益效果是:在确定像素局部区域方向后,可以为多向中值滤波提供选择滤波模板的参考,以便保留中值滤波保留细节,同时进一步加强了保护边缘和纹理的能力。
Description
技术领域
本发明涉及一种像素局部区域方向检测方法,特别地,涉及一种能够在超声图像去噪过程中像素局部区域方向检测方法,该检测方法的应用,能够在去除超声图像中噪声的同时能够保护超声图像边缘和纹理的细节部分。
背景技术
超声检查是一种医疗诊断技术,由于其价格便宜、对人体安全以及较好的适应性,被广泛用于各种医疗检查。20世纪70年代广泛发展应用的超声诊断技术,已经经历了从静态向动态图像、从黑白向彩色图像、从二维图像向三维图像、从反射法向透射法的发展。特别是血管内超声、三维成像、新型声学造影剂的应用使超声诊断又上了一个新台阶,目前已成为临床多种疾病诊断的首选方法,并成为一种非常重要的多种参数的系列诊断技术。
超声成像的原理是利用回声与原声波的差异产生图像,超声波经物体反射后会产生变化,变化与物体的形状特性有关,所以可以根据反射波来确定物体的外形。临床应用上,超声射入体内,由表面到深部,将经过不同声阻抗和不同衰减特性的器官与组织,从而产生不同的反射与衰减。这种不同的反射与衰减是构成超声图像的基础。将接收到的回声,根据回声强弱,用明暗不同的光点依次显示在影屏上,则可显出人体的断面超声图像。人体器官表面有被膜包绕,被膜同其下方组织的声阻抗差大,形成良好界面反射,超声象图上出现完整而清晰的周边回声,从而显出器官的轮廓,根据周边回声能判断器官的形状与大小。与其他医学影像诊断相比,超声诊断具有无损、无痛、价廉、方便、实时等优点,因此超声诊断成为医学临床诊断的重要手段之一。
然而,超声检查获得的图像受到噪声的影响,导致了利用超声诊断的巨大挑战。超声图像对比度差、影像中表征组织特性的特征不明显,特别是由于超声成像的固有特点,相位不同的散射回波相互干涉产生斑点噪声噪声成为影响超声诊断质量的关键因素。斑点噪声既降低了图像的质量影响诊断结果,又给图像的边缘检测、特征提取等后续处理带来困难。
为了克服上述缺陷,超声图像去噪技术一直成为研究热点。图像去噪主要是去除均匀区域的噪声。一般情况下,区域内部值比较突兀的像素被认为是噪声。超声图像中斑点噪声的抑制方法有两类。第一类是混合方法,这类方法将某方式得到的一组同一目标的图像进行相干平均,以去除随机斑点噪声。具体方式包括在不同的时间、以不同的扫描频率或从不同的空间位置对组织进行扫描。此方法较为成熟,但实现过程较为繁琐。第二类方法是滤波方法。各种去噪方法被相继提出,如中值滤波。中值滤波是一种非线性滤波技术,在某些情况下该滤波技术可以在去除图像噪声的同时保留图像细节。由于边缘是视觉关注的重要部分,因此希望保留边缘区域;而均匀部分由于存在一些纹理信息,也不能像普通图像一样简单平滑。因而超声去噪的同时更希望能保留边缘和纹理区域的信息。虽然中值滤波技术能够保护部分细节,但是由于它对所有区域使用同一种方式操作造成了边缘和纹理的平滑过度问题。又例如,小波滤波方法,该方法在去除噪声的同时会不同程度地丢失了一些图像细节。又例如维纳滤波、小波域内细节抛弃法、小波软阈值去噪等非线性阈值算法,这些方法在抑制斑点噪声的同时不同程度地降低了图像的分辨率。
其他的去噪方法中代表性的有申请号为CN201410094044.5的专利文献中提到的技术,其公开了一种超声图像斑点噪声抑制方法,包括以下步骤:1.从多个空间维度拍摄提取含噪超声图像,对每个空间维度内的含噪图像进行点对点智能化识别,获得图像的组织信息;2.对每个空间维度内的组织信息进行像素点区域属性分类,利用局部几何结构将像素点划分为脉冲区域、边缘细节区域;3.根据像素点的分类,分别对每个空间维度内的含噪组织信息进行斑点噪声抑制,获得单维度去噪超声图像;4.将每个空间维度的单维度去噪超声图像合成为超声去噪图像。该方法针对不同的属性采用不同的斑点去噪方法,能滤除脉冲噪声而细节点保持不变,但从实际效果来看,其在去噪的同时,边缘保持度并不理想。
发明内容
为了在去噪同时更好地保留边缘和纹理部分,提高去噪后超声图像的视觉效果,本发明提供了一种像素局部区域方向检测方法,包括以下步骤:
确定方向滤波算子;
确定滤波算子的至少一个方向参数;
依次采用滤波算子的方向参数对图像进行滤波,得到至少一个滤波值;
选择图像滤波值中的最大值对应的方向参数所指示的方向为每个像素的方向。
优选地,所述滤波算子为
其中,(m,n)为二维滤波算子中的位置,σd为高斯函数的标准差,λd为sin函数的波长,θd为滤波算子的方向参数。所述依次采用滤波算子的方向参数对图像进行滤波的步骤为:使用方向参数θd依次对图像像素进行滤波,获得P个滤波值确定P个滤波值中最大值对应的θd为像素的方向。
所述滤波算子的方向参数θd为{0,π}中的P个数值,P≥2。更为优选地,所述滤波算子的方向参数θd∈{0,π/4,π/2,3π/4}。确定四个值中最大值Max{γ0(x,y),γπ/4(x,y),γπ/2(x,y),γ3π/4(x,y)}对应的θd为像素的方向。
更为优选地,σd=10,λd=10。
本发明有益效果是:在确定像素局部区域方向后,可以为多向中值滤波提供选择滤波模板的参考,以便保留中值滤波保留细节,同时进一步加强了保护边缘和纹理的能力。使得在去除超声图像噪声的同时,保留有助于诊断的边缘和纹理区域细节,增强图像的视觉效果,进一步提高通过超声诊断的准确程度。
附图说明
图1示出了本发明中σd=10,λd=10时,θd∈{0,π/4,π/2,3π/4}四个方向检测滤波器;
图2示意了一种典型的5像素对称结构的中值滤波模板;
图3其示出了本发明中θt∈{0,π/4,π/2,3π/4}四个方向参数所对应的非对称结构的中值滤波模板;
图4(a)为医学切片图像;
图4(b)为由图4(a)生成的超声模拟图像;
图4(c)为图4(b)添加斑点噪声后的图像;
图5为采用本发明的去噪方案与采用传统中值滤波方法、采用维纳滤波方法去噪后的图像对比效果图。
图6为本发明的技术方案、传统中值滤波、维纳滤波、不滤波之间的峰值信噪比PSNR分析结果图。
具体实施方式
以下详细描述本发明技术方案的具体实施方式,本领域的技术人员应当清楚,以下具体实施方式为本发明最优或者较优选的实现方法,而不是对本发明保护范围的限定。
方向是边缘与纹理的重要特征之一,由于噪声较大的挑战,对于具有强噪声、低对比度和纹理复杂的超声图像,方向检测是一个巨大挑战,相应的滤波器及滤波方法应满足下列要求:
1、不受噪声影响。特别地,应该具有从具有强噪声的区域中进行方向检测的能力。
2、滤波器应该适应低对比度条件。多数情况下,两个边缘的两侧差异特别小,从而使得方向检测变得困难。
3、该方法应该能够在复杂环境下检测不同的方向。
为此,在本发明的第一实施方式中,提供了一种图像像素局部区域方向检测方法,用于每个像素的局部区域方向,以选择多向去噪滤波模板。该方法包括以下步骤:
确定方向滤波算子;
确定滤波算子的至少一个方向参数;
依次采用滤波算子的方向参数对图像进行滤波,得到至少一个滤波值;
选择图像滤波值中的最大值对应的方向参数所指示的方向为每个像素的方向。
在本发明的第二实施方式中,提供了一种方向滤波算子,该算子可作为本发明的第一实施方式中所述的方向滤波算子。该算子如下:
其中,(m,n)为二维算子中的位置,σd为高斯函数的标准差,λd为sin函数的波长,θd为滤波算子的方向参数。
上述滤波算子的方向参数θd可以为{0,π}中的任意一个或多个数值。优选地,上述滤波算子的方向参数θd∈{0,π/4,π/2,3π/4},此时滤波器与局部区域方向基本一致。即意味着,式(1)包含了四个模板。
在此基础上,使用这四个模板依次对图像Id(x,y)进行滤波,获得滤波结果γ(x,y)=|Id(x,y)*gd|,其中Id(x,y)为原始图像,gd为式(1)中的滤波算子。得到每个像素(x,y)的四个滤波值γ0(x,y),γπ/4(x,y),γπ/2(x,y),γ3π/4(x,y)。另一方面,若方向参数θd为{0,π}中的P个数值,则根据式(1)可以确定P个模板,使用这P个模板依次对图像Id(x,y)进行滤波,得到P个滤波值
进一步,确定每四个值中最大值Max{γ0(x,y),γπ/4(x,y),γπ/2(x,y),γ3π/4(x,y)}对应的θd为像素的方向。若方向参数θd为{0,π}中的P个数值,则可得到P个滤波值则确定P个值中最大值对应的θd为每个像素的方向。
更进一步,σd=10,λd=10时,滤波效果最佳。附图1示出了σd=10,λd=10时,θd∈{0,π/4,π/2,3π/4}四个方向检测滤波器。本领域的技术人员应当清楚,上述σd和λd也可以是其他的取值,并不受取值为10的限制,可以根据图像的性质和特征进行调整。
在本发明的第三实施方式中,提供了一种中值滤波方法,优选地,该方法包括以下步骤:
获取至少一个带方向特性的中值滤波模板;
根据确定的像素方向,选择中值滤波模板;
使用选择的中值滤波模板对每个像素进行滤波。
中值滤波是一种典型的非线性数字滤波器,经常被用于图像去噪,因为其能够在去除噪声的同时能过很好地保护边缘,中值滤波的模板具有多样性,现有的中值滤波模板大多都是对称结构的,如附图2所示,其是一种典型的5像素对称结构的中值滤波模板,其具有方向性,对称结构的模板并不利于超声图像去噪,因为其不能有效地保护边缘和纹理。也就是说,附图2中的滤波模板呈对称分布,类似于圆形模板,较适用于均匀区域的噪声去除,然而在超声图像中的边缘和纹理区域,其效果并不明显。
本发明设计的带有方向性的中值滤波模板不仅考虑了均匀区域的去噪问题,同时也考虑了边缘和纹理区域,因此其保留了对称滤波模板(例如图2所示的圆形模板)的大部分,而不影响均匀区域的去噪,同时又采用了非对称分布,有利于去除边缘和纹理区域的噪声。
对此,本发明提供了非对称性中值滤波模板,优选地,如图3所示,其示出了θt∈{0,π/4,π/2,3π/4}四个方向参数所对应的非对称结构的中值滤波模板,在图3(a)所示的5像素中值滤波模板中,第一模板为:第一行和第五行像素值为0,其他的均为1;第二模板为:第一行的第一二像素、第二行的第一像素、第四行的第五像素、第五行的第四五像素值为0,其他的均为1;第三模板为:第一列和第五列像素值为0,其他的均为1;第四模板为:第一列的第四五像素、第二列的第五像素、第四列的第一像素、第五列的第一二像素值为0,其他的均为1。
在像素方向方面,如上所述,可根据P个滤波值中最大值对应的θd为每个像素的方向。确定了像素的方向,即可选择对应的中值滤波模板。优选地,当θt∈{0,π/4,π/2,3π/4}时,若
Max{γ0(x,y),γπ/4(x,y),γπ/2(x,y),γ3π/4(x,y)}=γ0(x,y)
则选择附图3(a)的中值滤波模板;若
Max{γ0(x,y),γπ/4(x,y),γπ/2(x,y),γ3π/4(x,y)}=γπ/4(x,y)
则选择附图3(b)的中值滤波模板;其他的情况与之类似。
确定好中值滤波模板后,对该像素进行中值滤波。
本发明的第四实施方式中提供了一种超声图像去噪方法,当然,本领域的技术人员应当理解,只要具有与超声图像类似的噪声、边缘、纹理分布的图像均可使用本发明提供的图像去噪方法进行去噪。
对整幅图像而言,去噪过程为一个遍历过程,也就是说,每一个像素都将被处理,每一个像素的处理步骤如下:
确定局部区域方向。在该步骤中,使用多个方向滤波器确定每一像素的局部区域方向,优选地,多个方向滤波器的滤波结果中最大值对应的θd被确定为每个像素的方向。
多向中值滤波。根据上一步骤中确定的像素方向θd选择中值滤波模板进行滤波,得到滤波结果其中,为根据像素方向θd选择滤波模板后的中值滤波函数。中值滤波函数为本领域公知的滤波函数,本领域的技术人员在知道滤波模板的前提下即能够确定本发明不再敖述。
在本发明的第五实施方式中,对第四实施方式提供的去噪方法进行了进一步优化。在第四实施方式多向中值滤波之后,还包括确定去噪结果的步骤。初步试验发现,上一步骤中,经过多向中值滤波后的像素不能作为去噪的最终结果,因为对于少数像素而言,其实真严重,为此,本实施方式中,对于最终的去噪结果采用如下确定方法:判断每个像素位置滤波后的值和原来值的差距,当相差不多时,最终的滤波结果保留原值,否则使用滤波后的值。即,
其中,F(x,y)为滤波后的像素值,Id(x,y)为原图像素值,t为阈值。
也就是说,如果滤波后的像素值F(x,y)与原图像素值Id(x,y)的差的绝对值超过了阈值t,则将原图像素Id(x,y)作为滤波后的像素值,也就意味着不进行去噪;反之,将滤波后的像素值F(x,y)作为最终的去噪结果。
最为优选地,t=5时,第五实施方式的去噪结果最好,当然,阈值t的选择并不局限于5。
附图4展示了用于说明本发明技术方案的实验图像,图4(a)为肾部切片图像,图4(b)为合成的实验图像,图4(c)为图4(b)添加噪声参数σ2=0.1后的图像。
图5(a)为采用本发明的去噪方案对图4(c)进行去噪后的图像,图5(b)为采用传统中值滤波方法去噪后的图像,图5(c)为采用维纳滤波方法去噪后的图像。通过三者对比可以发现,本发明的去噪效果最佳,其能够保留更多的边缘和纹理细节又能够去除斑点噪声,相比而言,维纳滤波过于平滑,而传统中值滤波则保留了太多的细节而不分均匀区域和边缘、纹理。
图6为本发明的技术方案、传统中值滤波、维纳滤波、不滤波之间的峰值信噪比PSNR分析。
其中I为原图,即图4(c),Id为去噪之后的图像,i,j为像素的横坐标和纵坐标,L为图像的灰阶,图4(b)中L=255。图6明显反映出本发明提供的技术方案具有更高的峰值信噪比。
Claims (8)
1.一种像素局部区域方向检测方法,包括以下步骤:
确定方向滤波算子;
确定方向滤波算子的至少一个方向参数;
依次采用方向滤波算子的方向参数对图像进行滤波,得到至少一个滤波值;
选择图像滤波值中的最大值对应的方向参数所指示的方向为每个像素的方向。
2.如权利要求1所述的像素局部区域方向检测方法,其特征在于:所述滤波算子为
其中,(m,n)为二维滤波算子中的位置,σd为高斯函数的标准差,λd为sin函数的波长,θd为滤波算子的方向参数。
3.如权利要求2所述的像素局部区域方向检测方法,其特征在于:所述滤波算子的方向参数θd为{0,π}中的P个数值,P≥2。
4.如权利要求2所述的像素局部区域方向检测方法,其特征在于:所述滤波算子的方向参数θd∈{0,π/4,π/2,3π/4}。
5.如权利要求3所述的像素局部区域方向检测方法,其特征在于:所述依次采用滤波算子的方向参数对图像进行滤波的步骤为:使用方向参数θd依次对图像像素进行滤波,获得P个滤波值
6.如权利要求5所述的像素局部区域方向检测方法,其特征在于:确定P个滤波值中最大值对应的θd为像素的方向。
7.如权利要求4所述的像素局部区域方向检测方法,其特征在于:确定四个值中最大值Max{γ0(x,y),γπ/4(x,y),γπ/2(x,y),γ3π/4(x,y)}对应的θd为像素的方向。
8.如权利要求2-7中任意一项权利要求所述的像素局部区域方向检测方法,其特征在于:σd=10,λd=10。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410510938.8A CN104318520A (zh) | 2014-09-28 | 2014-09-28 | 一种像素局部区域方向检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410510938.8A CN104318520A (zh) | 2014-09-28 | 2014-09-28 | 一种像素局部区域方向检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104318520A true CN104318520A (zh) | 2015-01-28 |
Family
ID=52373746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410510938.8A Pending CN104318520A (zh) | 2014-09-28 | 2014-09-28 | 一种像素局部区域方向检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104318520A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339989A (zh) * | 2016-08-13 | 2017-01-18 | 浙江莱达信息技术有限公司 | 医学x光影像适应性降噪方法 |
CN107103619A (zh) * | 2017-04-19 | 2017-08-29 | 腾讯科技(上海)有限公司 | 一种头发纹理方向的处理方法、装置及系统 |
CN107392086A (zh) * | 2017-05-26 | 2017-11-24 | 深圳奥比中光科技有限公司 | 人体姿态的评估装置、系统及存储装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520888A (zh) * | 2008-02-27 | 2009-09-02 | 中国科学院自动化研究所 | 基于方向场的视网膜图像血管的增强方法 |
US7664326B2 (en) * | 2004-07-09 | 2010-02-16 | Aloka Co., Ltd | Method and apparatus of image processing to detect and enhance edges |
CN102393905A (zh) * | 2011-07-13 | 2012-03-28 | 哈尔滨工程大学 | 手背静脉模式纹理提取方法 |
-
2014
- 2014-09-28 CN CN201410510938.8A patent/CN104318520A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7664326B2 (en) * | 2004-07-09 | 2010-02-16 | Aloka Co., Ltd | Method and apparatus of image processing to detect and enhance edges |
CN101520888A (zh) * | 2008-02-27 | 2009-09-02 | 中国科学院自动化研究所 | 基于方向场的视网膜图像血管的增强方法 |
CN102393905A (zh) * | 2011-07-13 | 2012-03-28 | 哈尔滨工程大学 | 手背静脉模式纹理提取方法 |
Non-Patent Citations (6)
Title |
---|
冯星奎等: "方向加权中值滤波算法", 《中国图象图形学报》 * |
刘尚平: "眼底视网膜图像自动分析中质量改善_特征检测及配准的算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
张婷婷等: "基于方向滤波的超声图像边缘增强", 《计算机应用与软件》 * |
李欢利: "虹膜特征表达与识别算法研究", 《中国博士学位论文全文数据库》 * |
薛帅 等: "基于Gabor滤波器的Sobel算子图像边缘检测算法", 《信息技术》 * |
陈桂友: "自动指纹识别系统中的关键算法研究及应用", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339989A (zh) * | 2016-08-13 | 2017-01-18 | 浙江莱达信息技术有限公司 | 医学x光影像适应性降噪方法 |
CN106339989B (zh) * | 2016-08-13 | 2018-12-18 | 浙江莱达信息技术有限公司 | 医学x光影像适应性降噪方法 |
CN107103619A (zh) * | 2017-04-19 | 2017-08-29 | 腾讯科技(上海)有限公司 | 一种头发纹理方向的处理方法、装置及系统 |
CN107103619B (zh) * | 2017-04-19 | 2022-03-08 | 腾讯科技(上海)有限公司 | 一种头发纹理方向的处理方法、装置及系统 |
CN107392086A (zh) * | 2017-05-26 | 2017-11-24 | 深圳奥比中光科技有限公司 | 人体姿态的评估装置、系统及存储装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104299191A (zh) | 一种超声图像去噪方法及中值滤波方法 | |
Jaybhay et al. | A study of speckle noise reduction filters | |
CN102306377B (zh) | 一种超声图像去噪的方法和装置 | |
Zhang et al. | Speckle filtering of medical ultrasonic images using wavelet and guided filter | |
US9934554B2 (en) | Ultrasound imaging method/technique for speckle reduction/suppression in an improved ultra sound imaging system | |
CN102663708B (zh) | 基于方向加权中值滤波的超声图像处理方法 | |
CN101459766B (zh) | 超声图像增强与噪声抑制的方法及其装置 | |
Guo et al. | A novel approach to speckle reduction in ultrasound imaging | |
US20060184021A1 (en) | Method of improving the quality of a three-dimensional ultrasound doppler image | |
CN104680495A (zh) | 超声图像的自适应去噪方法 | |
CN101980286B (zh) | 结合双树复小波与双变量模型的sar图像降斑方法 | |
Roomi et al. | Speckle noise removal in ultrasound images using particle swarm optimization technique | |
CN104318520A (zh) | 一种像素局部区域方向检测方法 | |
Vanithamani et al. | Performance analysis of filters for speckle reduction in medical ultrasound images | |
CN103927715A (zh) | 一种超声图像斑点噪声抑制方法 | |
Rajabi et al. | Non-local adaptive hysteresis despeckling approach for medical ultrasound images | |
Abd-Elmoniem et al. | Real time adaptive ultrasound speckle reduction and coherence enhancement | |
Goyal et al. | SGS-SRAD filter for denoising and edge preservation of ultrasound images | |
Garg et al. | Speckle noise reduction in medical ultrasound images using coefficient of dispersion | |
JP2007229283A (ja) | 超音波診断装置 | |
Gupta et al. | Despeckling of medical ultrasound images: a technical review | |
Gupta et al. | An efficient technique for speckle noise reduction in ultrasound images | |
Esmaeili et al. | 3D segmentation of retinal cysts from SD-OCT images by the use of three dimensional curvelet based K-SVD | |
CN102819832A (zh) | 基于超复数小波幅值软阈值的斑点噪声抑制方法 | |
Kongo et al. | Dual-tree complex wavelet in medical ultrasounds images restoration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150128 |