基于幻方编码的中值滤波方法技术

本发明专利技术公开了一种基于幻方编码的中值滤波方法，其包括以下步骤：步骤一、读取图像中各像素的坐标和灰度值；步骤二、遍历图像，采用基于幻方编码的方法判断各像素是否为噪声像素或者是否位于线性细节区；步骤三、对于非噪声像素，令其灰度值保持不变，对于噪声像素，通过修改所述噪声像素的灰度值去除噪声。采用本发明专利技术对可能作为噪声像素处理的像素进行了是否位于线性细节区的判断，不会造成线性细节区内的像素缺失，确保图像信息完整；对线性细节区内的像素采用幻方编码的方法进行判断，计算量小；根据灰度中值与灰度最大值和最小值的平均值之间的差异情况选择合适的值替换噪声像素，减少了因去除图像噪声处理造成的图像模糊的现象。

【技术实现步骤摘要】

本专利技术涉及一种滤除图像噪声的方法，尤其是一种基于幻方编码技术的中值滤波 去除图像噪声的方法。
技术介绍
研究表明，人类获取的视觉图像信息在人类接受的信息中的比重达到四分之三， 是人类获取信息、表达信息和传递信息的重要手段。但在图像的生成、传输或变换过程中， 由于图像采集系统、传输媒介或成像系统的不完善，会引入不同程度的噪声，从而使输出图 像的质量产生退化。由于图像的许多后续处理(如边缘检测、模式识别或者图像分割等) 的成败在很大程度上依赖于噪声去除的效果和质量，因此，在数字图像处理领域，图像噪声 的滤除一直是最重要、最基本的研究课题之一。对受到噪声污染的退化图像的复原可以采用线性滤波方法来处理，在许多情况下 是很有效的，但是多数线性滤波具有低通特性，在去除噪声的同时也使图像的边缘变模糊， 不能获得较好的复原效果。而中值滤波是一种去除噪声的非线性处理方法，在某些条件下 既可去除噪声又能保护图像边缘，获得了较好的图像复原效果。因此自TUrkeyl971年提 出中值滤波技术以来，该技术就得到了广泛的研究，并被应用到语音处理、图像恢复等领域 中。标准中值滤波的基本原理是把数字图像或数字序列中某一像素的值，用该像素的 一个邻域中各像素值的中值代替。虽然标准中值滤波技术在衰减噪声的同时能较好地保护 图像边缘，但由于其仅考虑滤波窗内输入数据的排序信息，故在图像处理中会破坏一些重 要的图像细节，如细线等线形细节。显然，构成细线等线形细节的像素是信号，不是噪声。而 标准中值滤波方法会把构成线形细节的像素作为噪声进行处理，使得图像线性细节区的信 息丢失。另外，标准的中值滤波方法包括如下步骤选择一个(2n+l) X (2n+l)的滑动窗口， 使其沿图像数据的行或者列方向逐像素滑动；每次滑动后，对窗口内的像素值进行排序； 用排序得到的中间值代替窗口中心位置像素值。这种方法存在计算量大的问题，需要消耗 大量的时间和资源。
技术实现思路
本专利技术要解决的技术问题是提供一种有效保留图像线性细节区信息且算法简单 的。为解决上述技术问题，本专利技术所采取的技术方案是按照以下步骤处理步骤一、读取图像中各像素的坐标和灰度值；步骤二、遍历图像，采用基于幻方编码的方法判断各像素是否为噪声像素a.取以待处理的像素为中心的(2N+1) X (2N+1)的图像块为处理单元，对所述处 理单元内的像素所处位置按照幻方进行编码，得到的幻方编码序列为{Ρμ,ΡΜ，…，P2N+1,2N， P2N+1,2N+J，所述幻方编码序列对应的灰度值序列为^ali2,…，a2N+1,2N，a2N+1,2N+1}，其中N为自然数，Pi,j为(i，j)位置处的幻方编码，a,,,.为(i，j)位置处的灰度值；b.将所述处理单元内像素的灰度值按照从小到大进行排序，设排序后的灰度值序列为■(气，A2，…, +1)/2，…， _,，％}，所述排序后的灰度值序列中各灰度值所对应的像素在所 述处理单元中所处的位置序列为U1, ifikD^id，ij，所述位置序列对应的排序后的在上述序列中，e= (2N+1) X (2N+1)S为最小的灰度值%为最大的灰度值flfWv2为灰度中值；c.当W ε (1)，则待处理的像素是非噪声像素，式⑴中ε的取值范围为,； 当式⑴不成立时，根据下述公式(2)至⑷同时成立或者(5)至(7)是否同时 成立来判断待处理的像素是否位于线性细节区 式中，aM为待处理的像素灰度值 ε i的取值范围为^^ , Q 为幻和，Q = (2N + l)x = (2N + 1)x(2N2 +2N + 1) ε 2的取值范围为」^· ,； 10028] 当公式(2)、(3)和(4)同时成立或者公式(5)、(6)和(7)同时成立时，则待处理 的像素位于线性细节区；步骤三、对于非噪声像素，令其灰度值保持不变；对于噪声像素，通过修改所述噪声像素的灰度值去除噪声，其方法为 当I巧 -一^( 8 )，则用所述灰度中值Λ 替换噪声像素； 当上式(8)不成立时，用 替换噪声像素，即用灰度中值及位于灰 度中值前后各N个灰度值的平均值替换噪声像素， 式(8)中ε 3的取值范围为 上述的技术方案中，本专利技术这种去除图像噪声的方法，将所述处理单元内的灰度 值按照从小到大进行排序，根据最大的灰度值与最小的灰度值之差与所设定的阈值ε的 关系，判断是否为噪声像素；将大于阈值ε的像素利用幻方编码判断是否位于线性细节 区，其基本原理为首先判断待处理的像素的灰度值与所述处理单元内最大的灰度值或最 小的灰度值之差与所设定的阈值ε工的关系，若小于所设定的阈值S1，说明待处理的像素 的灰度值与最小的灰度值或者最大的灰度值接近；然后判断待处理的像素的灰度值与排序 后的灰度序列中第(2Ν+1)个或者倒数第(2Ν+1)个灰度值之差与所设定的阈值ε 2的关系， 若小于所设定的阈值ε 2，说明则待处理的像素位于排序后的灰度值序列中前(2Ν+1)个或 者后(2Ν+1)个；最后判断排序后的幻方编码序列中的前(2Ν+1)或者后(2Ν+1)的幻方编码 之和与幻和的关系，若等于幻和，说明排序后的前(2Ν+1)或者后(2Ν+1)个像素恰好位于一 条过待处理的像素点的直线上，则该待处理的像素位于线形细节区。采用本专利技术所产生的有益效果是(1)对以待处理像素为中心的(2Ν+1) X (2Ν+1)的图像块进行幻方编码，根据待处 理的像素与其邻域像素幻方编码及灰度值的关系，判断其是否为噪声，这样，在降低噪声的 同时尽可能地保留了图像细节，而且加快了程序的执行速度。(2)由于本专利技术不是简单地使用噪声像素邻域内像素的灰度中值对噪声进行替 换，而是根据灰度中值与灰度最大值和最小值的平均值之间的差异情况选择合适的值替换 噪声像素，使得图像的去噪效果更好，减少了因去除图像噪声处理造成的图像模糊的现象。附图说明图1是本专利技术的处理流程图；图2是本专利技术替换噪声像素方法的流程图；图3是所述处理单元内像素位置的幻方编码序列；图4是图3所对应像素的灰度值序列；图5是3阶幻方的幻方编码序列样例；图6是5阶幻方的幻方编码序列样例。具体实施例方式的处理流程，参见图1和图2，其按照以下步骤处 理步骤101，读取图像中各像素的坐标和灰度值并存储。本步骤中，存储的方法可以与现有技术相同，可以存储到函数f(x，y)中，将各个 像素的灰度值与坐标值对应存储。步骤102，取以待处理的像素为中心的(2N+1) X (2N+1)的图像块为处理单元，对 所述处理单元内的像素所处位置按照幻方进行编码，得到的幻方编码序列为{PijPu，…， P2N+1.2N' Pmi,mJ，所述幻方编码序列对应的灰度值序列为Iaia，aI,2'…，a2N+l,2N' a2N+l,2N+lJ ‘其中N为自然数，Py为(i，j)位置处的幻方编码，为(i，j)位置处的灰度值，所述幻 方编码序列和对应的灰度值序列参见图3和图4。图5和图6是3阶幻方和5阶幻方的幻 方编码序列样例。步骤103，将所述处理单元内像素的灰度值按照从小到大进行排序，设排序后的灰 度值序列为{w.,\華,.·，,aie丨，所述排序后的灰度值序列中各灰度值所对应的像 素在所述处理单元中所处的位置序列为U1, ifi(e+1)/2…‘” 本文档来自技高网...

【技术保护点】
一种基于幻方编码的中值滤波方法，其特征在于，按照以下步骤处理：步骤一、读取图像中各像素的坐标和灰度值；步骤二、遍历图像，采用基于幻方编码的方法判断各像素是否为噪声像素ａ．取以待处理的像素为中心的（２Ｎ＋１）×（２Ｎ＋１）的图像块为处理单元，对所述处理单元内的像素所处位置按照幻方进行编码，得到的幻方编码序列为｛ｐ↓［１，１］，ｐ↓［１，２］，…，ｐ↓［２Ｎ＋１，２Ｎ］，ｐ↓［２Ｎ＋１，２Ｎ＋１］｝，所述幻方编码序列对应的灰度值序列为｛ａ↓［１，１］，ａ↓［１，２］，…，ａ↓［２；步骤三、对于非噪声像素，令其灰度值保持不变；对于噪声像素，通过修改所述噪声像素的灰度值去除噪声，其方法为：当｜ａ↓［ｉ↓［（ｅ＋１）／２］］－ａ↓［ｉ↓［１］］＋ａ↓［ｉ↓［ｅ］］／２｜＜ε↓［３］（＊），则用所述灰度中值ａ↓［ｉ↓［（ｅ＋１）／２］］替换噪声像素；当上式（８）不成立时，用１／２Ｎ＋１＊ａ↓［ｉ↓［ｊ］］替换噪声像素，即用灰度中值及位于灰度中值前后各Ｎ个灰度值的平均值替换噪声像素，式（８）中ε↓［３］的取值范围为：（「ａ↓［ｉ↓［ｅ］］－ａ↓［ｉ↓［１］］／３２」，「ａ↓［ｉ↓［ｅ］］－ａ↓［ｉ↓［１］］／１６」）。Ｎ＋１，２Ｎ］，ａ↓［２Ｎ＋１，２Ｎ＋１］｝，其中Ｎ为自然数，ｐ↓［ｉ，ｊ］为（ｉ，ｊ）位置处的幻方编码，ａ↓［ｉ，ｊ］为（ｉ，ｊ）位置处的灰度值；ｂ．将所述处理单元内像素的灰度值按照从小到大进行排序，设排序后的灰度值序列为｛ａ↓［ｉ↓［１］］，ａ↓［ｉ↓［２］］，…，ａ↓［ｉ↓［（ｅ＋１）／２］］，…，ａ↓［ｉ↓［ｅ－１］］，ａ↓［ｉ↓［ｅ］］｝，所述排序后的灰度值序列中各灰度值所对应的像素在所述处理单元中所处的位置序列为｛ｉ↓［１］，ｉ↓［２］…ｉ↓［（ｅ＋１）／２］…ｉ↓［ｅ－１］，ｉ↓［ｅ］｝，所述位置序列对应的排序后的幻方编码序列为｛ｐ↓［ｉ↓［１］］，ｐ↓［ｉ↓［２］］…ｐ↓［ｉ↓［（ｅ＋１）／２］］…ｐ↓［ｉ↓［ｅ－１］］，ｐ↓［ｉ↓［ｅ］］｝，在上述序列中，ｅ＝（２Ｎ＋１）×（２Ｎ＋１）ａ↓［ｉ↓［１］］为最小的灰度值ａ↓［ｉ↓［ｅ］］为最大的灰度值ａ↓［ｉ↓［（ｅ＋１）／２］］为灰度中值；ｃ．当ａ↓［ｉ↓［ｅ］］－ａ↓［ｉ↓［１］］＜ε（１），则待处理的像素是非噪声像素，式（１）中ε的取值范围为：（「ａ↓［ｉ↓［ｅ］］－ａ↓［ｉ↓［１］］／４」，「ａ↓［ｉ↓［ｅ］］－ａ↓［ｉ↓［１］］／２」）...

【技术特征摘要】

【专利技术属性】
技术研发人员：张有会，王志巍，赵金，霍利岭，杨会云，杨芳芳，
申请(专利权)人：河北师范大学，
类型：发明
国别省市：13[中国|河北]