《基于方向加权中值滤波的超声图像处理方法.pdf》由会员分享，可在线阅读，更多相关《基于方向加权中值滤波的超声图像处理方法.pdf（16页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102663708 A (43)申请公布日 2012.09.12 CN 102663708 A *CN102663708A* (21)申请号 201210128319.3 (22)申请日 2012.04.27 G06T 5/00(2006.01) G06T 7/00(2006.01) (71)申请人 飞依诺科技 （苏州） 有限公司 地址 215123 江苏省苏州市工业园区星湖街 218 号 C8 座 5 楼 (72)发明人 马睿 (74)专利代理机构 苏州创元专利商标事务所有 限公司 32103 代理人 孙仿卫 赵艳 (54) 发明名称 基于方向加权中值滤波的超声图像。
2、处理方法 (57) 摘要 本发明公开了基于方向加权中值滤波的超声 图像处理方法， 该方法在中值滤波的基础上， 考虑 了图像的边缘方向和边缘幅度信息， 针对某个具 体的当前像素点， 只选取邻域窗口中与当前点边 缘方向相近的点做中值滤波， 而且中值的排序选 取是在用边缘幅度加权后的像素值上进行的， 本 发明的方法能够很好的区分超声图像的噪声和边 缘信息， 使得在去噪的同时有效地保护了超声图 像的边缘细节信息， 更加有利于后期的诊断与分 析。 而且本发明采用的基础算法是中值滤波， 算法 简单， 时间复杂度低， 处理耗时少， 更加符合超声 成像系统实时处理显示图像的要求。 (51)Int.Cl. 权利。
3、要求书 2 页 说明书 7 页 附图 6 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 7 页 附图 6 页 1/2 页 2 1. 基于方向加权中值滤波的超声图像处理方法， 其特征在于 ： 包括以下步骤 ： (1) 计算所述超声图像中每个像素点的边缘幅度和边缘方向 ； (2) 取所述图像中的一像素点作为当前点 ； (3) 根据所述当前点的边缘幅度或边缘方向判断该当前点是否为边缘点， 如果不是边 缘点， 则对该当前点进行自适应加权中值滤波处理 ； 如果是边缘点， 则对该当前点进行下面 步骤 (4)(5)(6) 的处理 ； (4)取该当前点的M*N邻域。
4、， 选择该M*N邻域中与该当前点的边缘方向相关的点作为处 理点 ； (5) 取所述处理点的像素值和边缘幅度值， 对所述处理点的像素值用相应的边缘幅度 值加权， 对加权后的值进行排序， 找出中间值 ； (6) 用所述中间值对应的像素值替换该当前点的像素值 ； (7) 取下一个像素点作为当前点进行步骤 (3) 的处理， 直到处理完所述图像中的所有 像素点， 最后输出处理后的图像。 2. 根据权利要求 1 所述的基于方向加权中值滤波的超声图像处理方法， 其特征在于 ： 在步骤 (1) 中， 所述像素点的边缘幅度和边缘方向的计算方法包括 ： (1a) 计算该像素点的水平方向和垂直方向的梯度， 设定该像。
5、素点的位置为 (i， j)， 该 像素点的像素值为 I(i， j)， 则该像素点的水平方向的梯度 Gx(i， j)、 垂直方向的梯度 Gy(i， j) 的计算方式如下 ： Gx(i， j) I(i， j+1)-I(i， j-1) Gy(i， j) I(i+1， j)-I(i-1， j) ； (1b) 计算该像素点的梯度幅度 GradientAm 和梯度方向 GradientAn， 计算方式如下 ： GradientAn atan(Gy/Gx) (1c) 判断该像素点的梯度幅度 GradientAm 是否大于预先设定的阈值， 如果该梯度幅 度GradientAm大于或等于所述阈值， 则取该梯度幅。
6、度GradientAm、 该梯度方向GradientAn 为该像素点的边缘幅度和边缘方向， 表示该像素点为边缘点 ； 如果该梯度幅度 GradientAm 小于所述阈值， 则设定该像素点的边缘幅度和边缘方向均为负值， 表示该像素点不是边缘 点。 3. 根据权利要求 1 所述的基于方向加权中值滤波的超声图像处理方法， 其特征在于 ： 在步骤 (1) 中， 所述像素点的边缘幅度和边缘方向的计算方法包括 ： (1a) 选取不同方向的边缘检测算子 ； (1b) 将每个方向的边缘检测算子中的数值与该像素点的邻域中的对应像素值相乘后 求和， 求和后取绝对值作为输出值 ； (1c) 在不同方向的输出值中， 。
7、将最大值减去最小值得差值， 如果该差值大于或等于预 先设定的阈值， 则取所述不同方向的输出值中的最大值为该像素点的边缘幅度值， 取该最 大值对应的边缘检测算子的方向为该像素点的边缘方向， 表示该像素点为边缘点 ； 如果该 差值小于所述阈值， 则设定该像素点的边缘幅度和边缘方向均为负值， 表示该像素点不是 边缘点。 4. 根据权利要求 1 所述的基于方向加权中值滤波的超声图像处理方法， 其特征在于 ： 在步骤 (4) 中， 从所述 M*N 邻域中选择与所述当前点的边缘方向相关的点作为处理点的方 法包括 ： 权 利 要 求 书 CN 102663708 A 2 2/2 页 3 (4a) 选取不同方。
8、向的取点模板 ； (4b) 从所述不同方向的取点模板中， 取与所述当前点的边缘方向最接近的方向对应的 取点模板 ； (4c) 在所述 M*N 邻域中， 取被所述最接近的方向对应的取点模板击中的点作为处理 点。 5. 根据权利要求 1 所述的基于方向加权中值滤波的超声图像处理方法， 其特征在于 ： 在步骤(4)中， 从所述M*N邻域中， 选择边缘方向与所述当前点的边缘方向的差值小于预先 设定的阈值的点作为处理点。 6. 根据权利要求 1 所述的基于方向加权中值滤波的超声图像处理方法， 其特征在于 ： 在步骤 (5) 中， 所述处理点的像素值用相应的边缘幅度值进行加权处理指得是该像素值与 相应的边。
9、缘幅度值进行相乘处理。 权 利 要 求 书 CN 102663708 A 3 1/7 页 4 基于方向加权中值滤波的超声图像处理方法 技术领域 0001 本发明涉及超声成像中的数据后处理技术， 特别涉及提高超声图像质量的图像处 理技术， 更具体地说， 涉及一种用于超声图像去噪和增强的方向加权中值滤波方法。 背景技术 0002 在超声图像形成过程中， 当超声波长与照射物体表面粗糙度相当时， 就会产生斑 点噪声， 这一现象可以用随机散射模型来解释。 这些噪声的存在， 使得超声图像的清晰度不 高， 这也是超声成像的主要缺点之一。超声图像特有的散斑噪声不仅使得超声图像的质量 较差， 尤其是掩盖和降低了。
10、图像某些细节信息， 也使得对图像细节的识别和分析变得更加 困难。为以后的病情诊断及定量分析， 图像特征提取和识别造成不利的影响。因此， 抑制这 些散斑噪声， 改善图像质量是超声图像分析和识别的重要预处理环节。 0003 超声图像的斑点噪声去除问题一直是国内外超声成像技术的重要课题之一。超 声图像去噪一般要求有效抑制散斑噪声， 同时要保留对后期分析和诊断有用的图像细节信 息。超声图像去噪的主要难点在于 ： 1) 散斑噪声可以大致看作为一种乘性噪声 ； 2) 噪声的 随机性质比较复杂 ； 3) 噪声易与图像细节相混， 而图像细节又复杂且多样。 0004 目前， 针对抑制超声图像散斑噪声， 改善图像。
11、质量， 人们已经提出了许多方法。1) 图像平均方法利用在不同时间、 不同频率或不同扫描方向得到同一目标的一系列图像， 将 它们平均形成一幅复合图像以提高图像的信噪比 (SNR)。尽管这种方法简单、 快速， 然而它 受到一些限制 ： 需要严格控制系列图像的形成， 而且需要图像的配准 ； 由于图像模糊的影 响， 一些小的细节 ( 例如小的血管、 纹理等 ) 会丢失， 因此降低了空间分辨能力。2) 自适 应加权中值滤波方法对图像中的每一个像素点的值， 用其局部邻域窗口的加权中值代替。 设定窗口中心点为 (in， in)， 窗口大小为 (2w+1)(2w+1)， 则窗口内各点加权系数计算为 ： wei。
12、ght(i， j) w(io， io)-a*d*2/m 0005 其中， 为取整运算， a 为系数， d 为点 (i， j) 到中心点 (in， in) 的距离， 2 和 m 分别为局部邻域窗口内的方差和均值。取邻域窗口内的加权中值代替中心点的像素值。 其中， 邻域窗口内各点的权值能够根据图像的局部统计特征来自动选取， 邻域窗口的大小 也能够根据其局部信噪比自动调节。相对于单纯的中值滤波， 这种自适应加权处理在保留 细节方面取得了一定的效果。 但是对窗口的选择很敏感， 限制了处理效果， 在去除噪声的同 时， 也会导致一些细微细节的损失。3) 小波阈值收缩是一类重要的去除散斑噪声的方法， 这种方。
13、法主要基于 Donoho 提出的小波软阈值 (soft-thresholding) 收缩去噪 ： 首先将图 像小波分解， 设定一个阈值， 对于小于阈值的小波系数用零代替， 而对于大于阈值的小波系 数用它减去阈值代替， 被处理后的小波系数作小波反变换， 就可以得到重建图像。同时， 还 有许多对于小波阈值的改进算法。然而， 上述方法主要是针对高斯分布的噪声， 解决其它 分布的噪声时， 效果不甚理想 ； 同时， 小波阈值的选取也是一个需要着重解决的问题， 既不 能选的太大而丢失一些边缘和局部细节， 也不能选的太小而对噪声的抑制不充分。如何选 择小波变换的尺度和阈值尚无确定方法。4) 基于各向异性扩散。
14、模型的去噪方法。各向异 说 明 书 CN 102663708 A 4 2/7 页 5 性扩散实际上是一个非线性的偏微分方程， 由图像的梯度来决定扩散速度， 能够同时兼顾 噪声消除和特征保持两个方面。目前， 以 Perona-Malik(PM) 模型为代表的各向异性扩散 方程 (Anisotropic Diffusion， AD) 得到了广泛的应用。Yu 等将各向异性扩散方程应用到 散斑抑制中， 提出了去除散斑噪声的各向异性扩散模型 (Speckle Reducing Anisotropic Diffusion， SRAD)。Yu 等修正了扩散系数， 使扩散方程能够根据图像噪声的情况而调整扩 散。
15、系数， 并且能够对图像的细节信息更加敏感。但 SRAD 模型也存在明显的缺陷 ： 模型中尺 度函数是由初始图像中尽可能大的均匀区域计算得到的， 模型的关键是如何选取图像中一 个尽可能大的， 合适的均匀区域， 这一区域的选取往往很大程度上影响着扩散结果， 给实验 结果带来较大的偶然性。此外， SRAD 模型所使用的局部统计信息实际上是各向同性的， 这 也背离了各向异性扩散算法的本质。 而且各向异性扩散需要多次迭代处理才能取得比较好 的效果， 迭代又是非常耗时的， 所以该算法很难满足超声实时处理的需求。 0006 总结一下， 现有的去除超声图像散斑噪声的方法有两个比较关键的不足需要提 高 ： 一是。
16、在对图像平滑去除噪声的同时， 对细节和边缘的模糊比较严重， 导致了一些边缘细 节的损失， 影响了后期的诊断与分析 ； 二是处理的时间复杂度比较高， 难以满足超声成像系 统实时处理显示图像的需求。 发明内容 0007 本发明的目的是提供基于方向加权中值滤波的超声图像处理方法， 该处理方法在 去噪的同时有效地保护了超声图像的边缘细节信息， 而且本处理方法采用的基础算法是中 值滤波， 算法简单， 时间复杂度低。 0008 为达到上述目的， 本发明采用了以下技术方案 ： 基于方向加权中值滤波的超声图 像处理方法， 包括以下步骤 ： 0009 (1) 计算所述超声图像中每个像素点的边缘幅度和边缘方向 ；。
17、 0010 (2) 取所述图像中的一像素点作为当前点 ； 0011 (3) 根据所述当前点的边缘幅度或边缘方向判断该当前点是否为边缘点， 如果不 是边缘点， 则对该当前点进行自适应加权中值滤波处理 ； 如果是边缘点， 则对该当前点进行 下面步骤 (4)(5)(6) 的处理 ； 0012 (4)取该当前点的M*N邻域， 选择该M*N邻域中与该当前点的边缘方向相关的点作 为处理点 ； 0013 (5) 取所述处理点的像素值和边缘幅度值， 对所述处理点的像素值用相应的边缘 幅度值加权， 对加权后的值进行排序， 找出中间值 ； 0014 (6) 用所述中间值对应的像素值替换该当前点的像素值 ； 001。
18、5 (7) 取下一个像素点作为当前点进行步骤 (3) 的处理， 直到处理完所述图像中的 所有像素点， 最后输出处理后的图像。 0016 步骤 (1) 的一种具体实施方式， 所述像素点的边缘幅度和边缘方向的计算方法包 括 ： 0017 (1a) 计算该像素点的水平方向和垂直方向的梯度， 设定该像素点的位置为 (i， j)， 该像素点的像素值为 I(i， j)， 则该像素点的水平方向的梯度 Gx(i， j)、 垂直方向的梯度 Gy(i， j) 的计算方式如下 ： 说 明 书 CN 102663708 A 5 3/7 页 6 0018 Gx(i， j) I(i， j+1)-I(i， j-1) Gy(。
19、i， j) I(i+1， j)-I(i-1， j) ； 0019 (1b) 计算该像素点的梯度幅度 GradientAm 和梯度方向 GradientAn， 计算方式如 下 ： 0020 GradientAn atan(Gy/Gx) 0021 (1c) 判断该像素点的梯度幅度 GradientAm 是否大于预先设定的阈值， 如果该 梯度幅度 GradientAm 大于或等于所述阈值， 则取该梯度幅度 GradientAm、 该梯度方向 GradientAn 为该像素点的边缘幅度和边缘方向， 表示该像素点为边缘点 ； 如果该梯度幅度 GradientAm 小于所述阈值， 则设定该像素点的边缘幅度。
20、和边缘方向均为负值， 表示该像素 点不是边缘点。 0022 步骤 (1) 的另一种具体实施方式， 所述像素点的边缘幅度和边缘方向的计算方法 包括 ： 0023 (1a) 选取不同方向的边缘检测算子 ； 0024 (1b) 将每个方向的边缘检测算子中的数值与该像素点的邻域中的对应像素值相 乘后求和， 求和后取绝对值作为输出值 ； 0025 (1c) 在不同方向的输出值中， 将最大值减去最小值得差值， 如果该差值大于或等 于预先设定的阈值， 则取所述不同方向的输出值中的最大值为该像素点的边缘幅度值， 取 该最大值对应的边缘检测算子的方向为该像素点的边缘方向， 表示该像素点为边缘点 ； 如 果该差值。
21、小于所述阈值， 则设定该像素点的边缘幅度和边缘方向均为负值， 表示该像素点 不是边缘点。 0026 步骤 (4) 的第一种具体实施方式， 从所述 M*N 邻域中选择与所述当前点的边缘方 向相关的点作为处理点的方法包括 ： 0027 (4a) 选取不同方向的取点模板 ； 0028 (4b) 从所述不同方向的取点模板中， 取与所述当前点的边缘方向最接近的方向对 应的取点模板 ； 0029 (4c) 在所述 M*N 邻域中， 取被所述最接近的方向对应的取点模板击中的点作为处 理点。 0030 步骤(4)的第二种具体实施方式， 从所述M*N邻域中， 选择边缘方向与所述当前点 的边缘方向的差值小于预先设。
22、定的阈值的点作为处理点。 0031 步骤 (5) 具体实施时， 所述处理点的像素值用相应的边缘幅度值进行加权处理指 得是该像素值与相应的边缘幅度值进行相乘处理。 0032 由于上述技术方案的运用， 本发明与现有技术相比具有下列优点 ： 本发明的基于 方向加权中值滤波的超声图像处理方法， 在中值滤波的基础上， 考虑了图像的边缘方向和 边缘幅度信息， 针对某个具体的当前像素点， 只选取邻域窗口中与当前点边缘方向相近的 点做中值滤波， 而且中值的排序选取是在用边缘幅度加权后的像素值上进行的， 本发明的 方法能够很好的区分超声图像的噪声和边缘信息， 使得在去噪的同时有效地保护了超声图 像的边缘细节信息。
23、， 更加有利于后期的诊断与分析。而且本发明采用的基础算法是中值滤 波， 算法简单， 时间复杂度低， 处理耗时少， 更加符合超声成像系统实时处理显示图像的要 求。 说 明 书 CN 102663708 A 6 4/7 页 7 附图说明 0033 附图 1 为本发明的流程示意图 ； 0034 附图 2 为本发明中通过梯度计算图像像素点的边缘幅度和边缘方向的流程示意 图 ； 0035 附图 3 为本发明中通过边缘检测算子计算图像像素点的边缘幅度和边缘方向的 流程示意图 ； 0036 附图 4 显示了边缘检测算子作用的 2*2 邻域的取法示意图 ； 0037 附图 5 显示了四个不同方向的边缘检测算子。
24、 ； 0038 附图 6 为本发明中通过当前点的边缘方向选择处理点的流程示意图 ； 0039 附图 7 显示了四个不同方向的取点模板 ； 0040 附图 8 为本发明中通过当前点与邻域中的点的边缘方向差来选择处理点的流程 示意图。 具体实施方式 0041 下面结合附图来进一步阐述本发明。 0042 本发明提出了一种方向加权中值滤波算法， 该算法的基本思路如图 1 所示， 以超 声图像为背景， 在中值滤波的基础上， 考虑了图像的边缘方向和边缘幅度信息， 针对某个具 体的当前点， 只选取滤波窗口中与当前点边缘方向相近的点做中值滤波， 而且中值的排序 选取是在用边缘幅度加权后的像素值上进行的。 方向。
25、加权中值滤波算法考虑到边缘不同于 背景， 其灰度分布存在明显的差异， 而且这种差异具有方向性， 因此采用边缘幅度和边缘方 向结合的方法。 0043 基于方向加权中值滤波算法， 用于超声图像的去噪和增强。该算法的流程图参见 图 1， 首先计算图像 I 中每个像素点的边缘幅度和边缘方向 ； 然后在此基础上， 逐点处理， 针 对每一像素点， 根据边缘方向选择其 M*N 邻域中方向相关的点作为处理点 ； 最后， 针对邻域 中选出的处理点， 用边缘幅度进行加权中值滤波处理。 在导入超声图像数据后， 主要包括以 下步骤 ： 0044 (1) 计算超声图像中每个像素点的边缘幅度 EdgeAm 和边缘方向 E。
26、dgeAn ； 0045 (2) 取图像中的一像素点作为当前点 ； 0046 (3) 根据当前点的边缘幅度或边缘方向判断该当前点是否为边缘点， 如果不是边 缘点， 则对该当前点进行自适应加权中值滤波处理， 该自适应加权中值滤波处理的方法采 用的是现有技术， 具体参见 “背景技术” ； 如果是边缘点， 则对该当前点进行下面步骤 (4)(5) (6) 的处理 ； 0047 (4)取该当前点的M*N邻域， 选择该M*N邻域中与该当前点的边缘方向相关的点作 为处理点 ； 0048 (5) 取上述处理点的像素值和边缘幅度值， 对处理点的像素值用相应的边缘幅度 值加权， 对加权后的值进行排序， 找出中间值。
27、 ； 0049 (6) 用上述中间值对应的像素值替换当前点的像素值 ； 0050 (7) 取下一个像素点作为当前点进行步骤 (3) 的处理， 直到处理完图像中的所有 像素点， 最后输出处理后的图像。 说 明 书 CN 102663708 A 7 5/7 页 8 0051 该算法中， 通过边缘幅度控制减少了噪声对中值滤波的影响， 边缘方向信息的加 入又使得该算法对边缘方向特别敏感， 保留细节能力强， 减少了边缘模糊的程度。因此， 这 种改进的中值滤波算法， 在去除散斑噪声能力、 保持边缘能力以及处理速度等指标上均取 得了比较好的效果。 0052 上述算法中， 步骤 (1) 计算图像 I 中每点的。
28、边缘幅度和方向， 有很多种方法能够 实现。这里给出两种具体的实施方法 ： 1. 通过图像的梯度幅度和方向得到图像的边缘幅度 和方向 ； 2. 用不同方向的边缘检测算子得到图像的边缘幅度和方向。但不仅限于这两种方 法。 0053 方法 1 ： 通过图像的梯度幅度和方向得到图像的边缘幅度和方向的方法， 其流程 图参见图 2。首先， 根据水平方向的梯度 Gx(i， j) 和垂直方向的梯度 Gy(i， j) 计算图像的 梯度幅度 GradientAm 和梯度方向 GraditenAn， 计算方式如下 ： 0054 设定图像像素点的位置为 (i， j)， 该像素点的像素值为 I(i， j)， 则该像素点。
29、的 Gx(i， j)、 Gy(i， j) ： 0055 Gx(i， j) I(i， j+1)-I(i， j-1) Gy(i， j) I(i+1， j)-I(i-1， j) 0056 则该像素点的梯度幅度 GradientAm 和梯度方向 GraditenAn ： 0057 GradientAn atan(Gy/Gx) 0058 其次， 判断该像素点的梯度幅度值是否大于预先设定的阈值 GradientAmThresh。 如果大于或等于阈值， 则表示该像素点为边缘点， 把梯度幅度看作是边缘幅度， 梯度方向看 作是边缘方向 ； 如果小于阈值， 则认为无边缘 ( 表示该像素点不是边缘点 )， 把边缘幅。
30、度和 方向都设成负值 ( 比如 -1)， 以便在后面处理中区分像素点是否为边缘点。即 ： 0059 0060 对图像 I 的所有像素点逐点进行上述处理， 则得到每个像素点的边缘幅度 EdgeAm 和边缘方向 EdgeAn。 0061 方法 2 ： 通过不同方向的边缘检测算子得到图像边缘幅度和方向的方法， 参见图 3， 首先选取一组有不同方向的边缘检测算子 ； 其次把这组算子分别作用于图像上的一像素 点的邻域上， 得到每个方向上算子的输出 ； 比较各个方向上的输出， 如果差别比较大， 选取 最大的输出作为当前像素点的边缘幅度， 最大输出对应的算子的方向为边缘方向 ； 如果差 别不大， 则认为当前。
31、像素点不是边缘点。 0062 举例说明， 选取如图 5 所示的四个方向 (90 度， 0 度， 45 度， 135 度 ) 上的边缘检测 算子， 使用该组边缘检测算子计算图像上每点的边缘幅度和方向， 参见图 3， 首先， 针对图像 中的一像素点(i， j)， 取如图4所示的2*2邻域上的点 ； 然后， 把图5所示的四个方向的边缘 检测算子分别作用于图 4 的 2*2 邻域， 取每个算子模板与邻域模板中的对应数值相乘后求 和的绝对值， 作为输出值 0peratorAm am90， am0， am45， am135， 相应的方向 0peratorAn 90， 0， 45， 135 ； 接着， 比较。
32、四个方向上的输出值， 如果差别比较大， 也就是四个方向的 输出值中最大值减去最小值的差大于或等于预先设定的阈值 0peratorAmThresh， 则表示该 像素点为边缘点， 当前像素点 (i， j) 的边缘幅度为输出的最大值， 边缘方向为最大值对应 的边缘检测算子的方向， 假设输出的最大值为am45， 则边缘幅度为am45， 边缘方向为45。 如 说 明 书 CN 102663708 A 8 6/7 页 9 果输出差别很小， 也就是最大值减去最小值的差小于阈值， 则认为无边缘 ( 表示当前像素 点不是边缘点 )， 当前像素点的边缘幅度和方向都设成负值 ( 比如 -1)。用公式表示如下 ： 0。
33、063 0064 其中， index 为 max(OperatorAm) 对应的序号。 0065 对图像的所有像素点都按照上述步骤处理， 就能得到每个像素点的边缘幅度 EdgeAm 和方向 EdgeAn。 0066 上述基于方向加权中值滤波算法的步骤(4)中， 如何选择M*N邻域中与当前点(i， j) 的边缘方向相关的点作为处理点， 这里也给出了两种具体实现方法 ： 1. 以当前点的边缘 方向为准， 选择邻域中这个方向上的点为处理点 ； 2. 选择与当前点的边缘方向差小于预设 阈值的点为处理点。但不仅限于这两种方法。 0067 方法 1 ： 以当前点的边缘方向为准， 选择邻域中这个方向上的点为。
34、处理点， 方法 1 的流程图如图 6 所示。首先选取不同方向的取点模板， 如图 7， 图 7 显示了四个方向 (90 度， 0 度， 45 度， 135 度 ) 的取点模板， 判断当前点 (i， j) 的边缘方向 EdgeAn(i， j) 与哪个方向 (90 度， 0 度， 45 度， 135 度 ) 最接近， 也就是说 EdgeAn(i， j) 落在如下哪个角度区间里面 ： 67.5， 112.5， 0， 22.5U157.5， 180， 22.5， 67.5， 112.5， 157.5 ； 然后， 按照如下规则 选取邻域中沿着这个方向上的点 ： 0068 如果与 90 度最接近， 也就是落。
35、在角度区间 67.5， 112.5 中， 选择图 7 中模板 2a) 或者 2b) 击中的点作为处理点 ； 0069 如果与 0 度最接近， 也就是落在区间 0， 22.5U157.5， 180 中， 选择图 7 中模板 3a) 或者 3b) 击中的点作为处理点 ； 0070 如果与 45 度最接近， 也就是落在区间 22.5， 67.5 中， 选择图 7 中模板 4a) 或者 4b) 击中的点作为处理点 ； 0071 如果与 135 度最接近， 也就是落在区间 112.5， 157.5 中， 选择图 7 中模板 5a) 或 者 5b) 击中的点作为处理点。 0072 假设选出的处理点的个数为。
36、 k， 边缘幅值为 0073 LocalEdgeAm ea1， ea2， ea3， eak， 像素值为 0074 LocalPixel p1， p2， p3， pk。 0075 需要说明的是， 为了直观起见， 图 7 中只给出了 5*5 大小的取点模板 2a)-5a) 或者 2b)-5b)， 实际上这里的模板大小应该是M*N大小的。 其中， 2a)-5a)是针对边缘宽度是单点 的情况下的取点模板， 2b)-5b)是针对边缘宽度比较宽的情况下的取点模板。 如果选择模板 2a)-5a)， 容易出现误把噪声点增强， 如果选择模板 2b)-5b)， 容易把一些细小的边缘抹掉， 所以要根据具体脏器的边缘。
37、宽窄情况， 选择合适的取点模板。 0076 方法 2 ： 选择与当前点 (i， j) 的边缘方向的差小于预设阈值的点为处理点， 方法 2 的流程图如图 8 所示。在当前点 (i， j) 的邻域中的 N*M 个点中， 选择边缘方向 EdgeAn 与当 前点 (i， j) 边缘方向 EdgeAn(i， j) 的差小于设定的阈值 EdgeAnThresh 的点为处理点。假 设选中的处理点有 k 个， 它们的边缘幅度值为 LocalEdgeAm ea1， ea2， ea3， eak， 像 素值为 LocalPixel p1， p2， p3， pk。 说 明 书 CN 102663708 A 9 7/7。
38、 页 10 0077 上述基于方向加权中值滤波算法中， 步骤 (5)(6) 中涉及到的， 在选出了处理点并 拿到了它们的边缘幅度值和像素值后， 进行加权中值滤波处理， 首先， 用边缘幅度对像素值 加权， 即 ： LEAP ea1*p1， ea2*p2， ea3*p3， eak*pk ； 其次， 把加权后的值 LEAP 排序， 取中间值 median(LEPA) 所对应的像素值 px 代替当前点 (i， j) 处的像素值， 如果中间值为 ea2*p2， 则用 p2 代替点 (i， j) 的像素值。 0078 本发明的基于方向加权中值滤波的超声图像处理方法， 能够很好的区分超声图像 的噪声和边缘信。
39、息， 使得在去噪的同时有效地保护了超声图像的边缘细节信息， 更加有利 于后期的诊断与分析。而且本发明采用的基础算法是中值滤波， 算法简单， 时间复杂度低， 处理耗时少， 更加符合超声成像系统实时处理显示图像的要求。 说 明 书 CN 102663708 A 10 1/6 页 11 图 1 说 明 书 附 图 CN 102663708 A 11 2/6 页 12 图 2 说 明 书 附 图 CN 102663708 A 12 3/6 页 13 图 3 说 明 书 附 图 CN 102663708 A 13 4/6 页 14 图 4 图 5 图 6 说 明 书 附 图 CN 102663708 A 14 5/6 页 15 图 7 说 明 书 附 图 CN 102663708 A 15 6/6 页 16 图 8 说 明 书 附 图 CN 102663708 A 16 。