专利名称：一种去除实时超声图像斑点噪声的方法
技术领域：
本发明涉及超声诊断技术领域，具体涉及一种去除实时超声图像斑 点噪声的方法。
背景技术：
在超声成像中，当人体组织的结构尺寸与入射超声波波长相近或小 于波长时，超声束发生散射，相位不同的散射回波相互千涉产生斑点噪 声，它降低了超声图像的质量，使对比度较低的软组织中的正常组织和 肿瘤病变组织不易分别。对斑点噪声的统计特性的研究表明，常见的斑点噪声服从Rayleigh 分布，其均值与标准差成正比，这说明斑点噪声是乘性的。Jain在1989 年提出了一个乘性与加性噪声相組合的超声图像噪声模型式中 和""分别是乘性噪声和加性噪声，x表示未被噪声污染的信号，y 表示被噪声污染后观察到的信号。在超声医学图像中，加性噪声的作用相对于乘性噪声来说很小，因 此在实际应用中，可以忽略加性噪声的影响，所以式(l)可以简化为(2)目前，超声图像中斑点噪声的抑制方法有两类。第一类是复合方法。 这类方法将某方式得到的一组同 一 目标的图像进行相干平均，以去除随 机斑点噪声。具体方式包括在不同的时间、以不同的扫描频率或从不同
的空间位置对组织进行扫描。此方法较为成熟，但实现过程较为繁瑣。 第二类方法是滤波方法。从现有文献看，目前许多去除斑点噪声的技术比4交著名的包括应用在雷达凄t据处理中的Lee、 Kuan和Frost滤波器。 与此同时，不是由斑点统计模型推导的其他滤波器也已提出并应用到斑 点噪声抑制中。诸如均值率波器、中值滤波器、几何滤波器、小波变换 滤波器、形态率波器、维纳滤波，近几年还有^f艮多利用前几种算法和别 的算法的结合发展的新的斑点噪声去除方法。这些方法在抑制斑点噪声 的同时不同程度地降低了图像的分辨率，并且很难保证图像的实时性。目前的中低档B超中，大多利用FPGA采用帧相关处理方法来达到去 除斑点噪声的目的。帧相关处理是最常见也是最古老的一种消除斑点的 时间复合。帧相关处理一般是帧间的一阶递归滤波，可以表示为如下的 形式><")=qK"-i)+(i-"W")0< <i (3)式(2)中"")表示当前帧某线某采样点的实测值，yw表示^w的滤波值，而^"-D表示前一帧同样位置采样点的滤波值，滤波器系数"一般 是可选的。"越接近l，速归滤波器的带宽越窄，抑制斑点噪声的作用越 大，但图象的实时跟踪性越差。滤波器系数"越接近l,递归滤波器的带 宽越窄，抑制斑点噪声的作用越大。但是我们现有B超系统的帧相关处 理算法中，"一般取为0.5，再增大"，图像的实时跟踪性就变得很差， 而且采样FPGA来进行去噪成本相对来说比较高。发明内容本发明要解决的技术问题是提供一种去除实时超声图像斑点噪声 的方法，克服现有技术去除实时超声图像斑点噪声的方法，采用硬件实 现，图像实时跟踪性差，去噪声成本高的缺陷。本发明去除实时超声图 像斑点噪声的方法在PC上用软件实现，既保存图像细节，又保证图像实时处理速度。本发明为解决上述技术问题所采用的技术方案为一种去除实时超声图像斑点噪声的方法，包括步骤Al、读入源图像象素数据；A2、对所述源图像泉素数据进行同态变换；A3 、对同态变换后的数据进行多级非线性加权平均中值滤波计算；A4、对滤波计算后的数据进行同态反变换；A5、将向态反变换后的数据进行数字扫描变换后显示。所述步骤A2中的同态变换是对源图像中的每一点取以自然数e为 底的对数。所述步骤A3包括步骤定义中心在(i， j),大小为(2K+l)x(2K+l)的滤波窗W的四个子集 『1<formula>formula see original document page 6</formula>其中K取正整数，M-K-l2kK, N匿K-l^j^K。所述步骤A3还包括步骤 计算滤波.窗W的四个子集的中值2力，力<formula>formula see original document page 6</formula>其中M-K画l^i》K， N-K-l^j^K。
所述步骤A3还包括步骤 计算四个中值的平均值'Z(/, /) = * (, ，力+ Z2 (/,力+ Z3 (,',力+ Z4 (/， 乂))，其中M-K-l^ikK， N國S-l^j2K。所述步骤A3还包括步骤 4姿照下列数学式进4亍多级非线性加权平均中值滤波计算邻,力=1 (IM * z"") +^(") + #4 * &(" ) + #4*&(")) 4 z0,力 Z(g力 Z(z,力 Z(z,力其中M-K-l》i^K， N-K-l^j^K。所述步骤A4中按照下列数学式进行同态反变换 y(/，y') = ez(")_l .其中Y(i, j)为经过同态反变换的图像元素，其中，i=0、 1、、、M-1, j-O、 1、、、、 N-1，M、 N如图2所示。本发明的有益效果为本发明利用软件方法实现斑点噪声的去除， 代替传统的硬件实现方式，节省了成本；去除斑点噪声的同时，使图像 的边缘和细节得以很好的保存；保证了图像的实时跟踪性。
图1为本发明滤波窗及其子集示意图； 图2为本发明具体实施过程示意图。
具体实施方式
下面根据附图和实施例对本发明作进一 步详细说明:本发明基于Rayleigh噪声模型，由于斑点噪声服从负指数统计规律， 因此其为乘性噪声，表达式为<formula>formula see original document page 8</formula> (4)这里X表示未被噪声污染的信号，Y表示观测到的已被噪声污染的信号， n表示零均值方差为1的高斯白噪声。本发明方法主要分为以下几个步 骤同态变换即把图像中的每一点取以自然数e为底的对数。对于乘性噪声的同态变 换为G(Y)=ln(Y+l)( 5 )经过同态变换后，原来的乘性噪声变成了加性噪声，降低了处理的复杂 性。多级非线性加权平均中值滤波计算设X(.，.)表示图像元素序列，W表示中心在(i， j),大小为(2K+l)x(2K+l)的正方形滤波窗。那么滤波窗W的四个子集可定义为<formula>formula see original document page 8</formula> (7 - a)<formula>formula see original document page 8</formula> (7 - b)<formula>formula see original document page 8</formula> (7 - C)<formula>formula see original document page 8</formula>7 - d)如图1所示,W1(i，j)、 W2(i，j)、 W3(i，j)、 W4(i，j)、 分别表示以(i， j) 为中心，滤波窗在水平方向、右对角线方向、垂直方向、左对角线方向 所覆盖的图像元素所组成的序列。假定Zs(i，j)、 (s=l, 2, 3, 4)是四个子集元素的中值，即，<formula>formula see original document page 8</formula> (8 - a)<formula>formula see original document page 8</formula> (8 — b)<formula>formula see original document page 8</formula> (8 - c) <formula>formula see original document page 9</formula>)
然后求四个中值的平均值<formula>formula see original document page 9</formula>
(9) 多级非线性加;f又平均中值滤波的输出为<formula>formula see original document page 9</formula>(10)
通过对上述四个子集图像元素的中值进行加权平均，给与其中较大的值 以较大的权重，以期弥补图像灰度均值降低和噪声抑制能力弱的缺陷， 较好的恢复噪声图像。其中Kd, ^ffi""^表示求图像元素序列的中值。 同态反变4夹最后图像的输出由(6)式计算得出 y<formula>formula see original document page 9</formula>如图2所示，中X[i][j], Data[i][j], G[i][j]， Z(i, j)， Y(i, j)都对应图像 中的同一个点，其中，i = 0、 1、、、 M-l， j = 0、 1、、、、 N-l。 过程①表示把源图像X[M][N]中的每一个元素X[i][j]存入数组Data[i][j] 中，其中，i = 0、 1、、、 M-l， j = 0、 1、、、、 N-l;过程②表示把数组Data[M][N]中的每一个元素Data[i][j]取以自然数e为 底的对数，计算后的值存入数组G[i][j],即G[i][j] = ln(Data[i][j])，其中， i = 0、 1、、、 M-l, j = 0、 1、、、、 N-l;过程③表示把数组G[M][N]中的每一个元素G[i][j]经过多级非线性加权 平均中值滤波计算后的值存入数组Z[i][j], Z[i][j]的值由式(10)计算得 出，其中，i-K+l、 K+2、、、 M-l-K， j = K+l、 K+2、、、 M-l-K，即我们 对边缘的点不做处理；过程④表示把数组Z[M][N]中的每一个元素Z[i][j]经过同态反变换的值 存入数组Y[i][j]中，即}10[_/] = ￡'，-1,其中，i = 0、 1、、、 M國l， j = 0、 1、、、、 N-l;过程⑤表示把数组Y[M][N]中的值经过数字扫描变换后显示。 本发明的具体实施步骤读入源图像X，假设图像横向有N个采样点，纵向有M个釆样点。 在内存中申请四个大小为MxN的数組Data[M][N]、 G[M][N]、 Z[M][N]、 Y[M][N]分别用来存放源图像数据、同态变换后图像数据、多级非线性 加权平均中值滤波后图像数据、同态反变换后图像数据。 把源图像中处于第i行第.j列的点存入Data的第i行第j歹'〗，即Data[i][j] =X[i][j]，其中，i-0、 1、、、 M-l, j = 0、 1、、、、 N-l。 同态变换把数组Data[M][N]中的所有元素取以自然数为底的对数，所 得结果存入数组G[M][N]中，即G〖i][j] = ln(Data[i][j]+l),其中，i = 0、 l…M-l， j = 0、 1、、、、 N-l。多级非线性加权平均中值滤波以G[M][N]中的每一个元素G[i][j]为中 心，把其大小为(2K+l)x(2K+l)的正方形邻域内的元素分别按水平方向、 垂直方向、右对角线方向、左对角线方向分为四组，分别用『力',/)、 W，力、『3(")、『4(!'J)表示，如图1所示；然后分别求序列W，力、w，y)、 『3(z',力、『4(/,力的中值，分別用2,(,',力、Z2(/,/)、 Z3(/，y)、 24(/,力表示; 再计算四个中值的平均值郭,力，最后由式(IO),计算多级非线性加权 平均中值滤波的输出，把输出但存入到Z[i][j]中，其中，i = K+l、 K+2、、、 M-l-K, j=K+l、 K+2、、、 M-l-K,即我们对边缘的点不做处理。 同态反变换由式(11)把以Z[M][N]中的每一个元素Z[i][j]为指数， 以自然数e为底数计算所得的值存入数组Y[i][j]中，其中，i = 0、 1、、、 M-l， j = 0、 1、、、、 N-l。 把Y[M][N]中的数据进行数字扫描变换后显示。 本领域技术人员不脱离本发明的实质和精神，可以有多种变形方案 实现本发明，.以上所述仅为本发明较佳可行的实施例而已，并非因此局 限本发明的权利范围，凡运用本发明说明书及附图内容所作的等效结构 变化，均包含于本发明的权利范围之内。
权利要求
1、一种去除实时超声图像斑点噪声的方法，其特征在于，包括步骤A1、读入源图像象素数据；A2、对所述源图像象素数据进行同态变换；A3、对同态变换后的数据进行多级非线性加权平均中值滤波计算；A4、对滤波计算后的数据进行同态反变换；A5、将同态反变换后的数据进行数字扫描变换后显示。
2、 根据权利要求1所述的去除实时超声图像斑点噪声的方法，其 特征在于所述步骤A2中的同态变换是对源图像中的每一点取以自然 数e为底的对数。
3、 根据权利要求1所述的去除实时超声图像斑点噪声的方法，其 特征在于，所述步骤A3包括步骤定义中心在(i, j)，大小为(2K+l)x(2K+l)的滤波窗W的四个子集 <formula>formula see original document page 2</formula>其中K取正整数，i^K, j^K。
4、 根据权利要求3所述的去除实时超声图像斑点噪声的方法，其 特征在于，所述步骤A3还包括步骤计算滤波窗W的四个子集的中值Z,0',力(s=l ，2,3,4)，Z2(〖，/) = AffiX)/^V(『2 (/,_/))， Z3 (i,力=M五D崖(『3 (z.，力)， Z4(!',/) = 7^Z)/^V(『4 (! ，_/))，其中i^0， j^0。
5、 根据权利要求4.所述的去除实时超声图像斑点噪声的方法，其 特征在于，所述步骤A3还包括步骤 计算四个中值的平均值郭,_/) = *(Z, (/,力+ Z2 (/， /) +'Z3 (! ，力+ Z4 力)， 其中i^O, j^0。
6、 根据权利要求5所述的去除实时超声图像斑点噪声的方法，其 特征在于，所述步骤A3包括步骤按照下列数学式进4亍多级非线性加权平均中值滤波计算邻，力=1(製*， +酒*， + #4*， +辦，,力)其中kO， j、0。
7、根据权利要求1至6任一所述的去除实时超声图像斑点噪声的 方法，其特征在于，所述步骤A4中按照下列数学式进行同态反变换其中Y(i, j)为经过同态反变换的图像元素，kO， >0
全文摘要
一种去除实时超声图像斑点噪声的方法，包括步骤A1.读入源图像象素数据；A2.对所述源图像象素数据进行同态变换；A3.对同态变换后的数据进行多级非线性加权平均中值滤波计算；A4.对滤波计算后的数据进行同态反变换；A5.将同态反变换后的数据进行数字扫描变换后显示。本发明利用软件方法实现斑点噪声的去除，代替传统的硬件实现方式，节省了成本；去除斑点噪声的同时，使图像的边缘和细节得以很好的保存；保证了图像的实时跟踪性。
文档编号A61B8/00GK101396279SQ200710123680
公开日2009年4月1日 申请日期2007年9月29日 优先权日2007年9月29日
发明者刘明宇, 齐丽芸 申请人:深圳市蓝韵实业有限公司