一种电缆瓷套终端红外图像自适应去噪方法.pdf

上传人：zhu****_FC

文档编号：6164314

上传时间：2019-05-15

格式：PDF

页数：12

大小：3.32MB

《一种电缆瓷套终端红外图像自适应去噪方法.pdf》由会员分享，可在线阅读，更多相关《一种电缆瓷套终端红外图像自适应去噪方法.pdf（12页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103839239 A (43)申请公布日 2014.06.04 CN 103839239 A (21)申请号 201410071698.6 (22)申请日 2014.02.28 G06T 5/00(2006.01) (71)申请人华南理工大学地址 510640 广东省广州市天河区五山路 381 号 (72)发明人吴炬卓牛海清吴倩叶开发 (74)专利代理机构广州市华学知识产权代理有限公司 44245 代理人蔡茂略 (54) 发明名称一种电缆瓷套终端红外图像自适应去噪方法 (57) 摘要本发明公开了一种电缆瓷套终端红外图像自适应去噪方法，包括以。

2、下步骤： (1) 输入待去噪的电缆瓷套终端红外图像； (2)将图像分解为R颜色分量图像、 G颜色分量图像和B颜色分量图像； (3) 得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数； (4) 得到图像的细节小波系数的Bayes最优估计； (5)得到去噪后的图像； (6)得到去噪后的结果图像。具有可以在不知红外图像先验信息的情况下，自适应逐层选择每一尺度上的能量最优基小波等优点。 (51)Int.Cl. 权利要求书 2 页说明书 6 页附图 3 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请权利要求书2页说明书6页附。

3、图3页 (10)申请公布号 CN 103839239 A CN 103839239 A 1/2 页 2 1. 一种电缆瓷套终端红外图像自适应去噪方法，其特征在于，包括以下步骤： (1) 输入待去噪的电缆瓷套终端红外图像； (2) 将图像分解为 R 颜色分量图像、 G 颜色分量图像和 B 颜色分量图像； (3) 对于 R 颜色分量图像、 G 颜色分量图像和 B 颜色分量图像均采用自适应小波选择算法逐层确定最优基小波并将确定的各层基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数； (4) 对于 R 颜色分量图像、 G 。

4、颜色分量图像和 B 颜色分量图像均保留步骤 (2) 得到的近似小波系数，对步骤(2)得到的细节小波系数进行自适应Bayes估计，得到图像的细节小波系数的 Bayes 最优估计； (5) 对于每个颜色分量图像，利用步骤 (2) 得到的最高分解层数上的近似小波系数和步骤 (4) 得到的不同分解层数上的细节小波系数的 Bayes 最优估计进行图像重构，得到去噪后的图像； (6) 利用步骤 (5) 得到的去噪后的 R 颜色分量图像、 G 颜色分量图像和 B 颜色分量图像合成为原来的彩色图像，得到去噪后的结果图像。 2. 根据权利要求 1 所述的电缆瓷套终端红外图像自适应去噪方法。

5、，其特征在于，所述步骤 (3) 包括以下步骤： (a) 给定一个小波函数库 i:i=1,2,.,N 和小波分解层数 K ； (b) 从小波函数库中依次选择小波 i，对红外图像信号 s 进行第一层小波分解，得到近似系数和细节系数l=1,2,3 ； (c) 依次计算小波 i对应的近似系数的能量百分比直到 i=N ； (d) 比较计算得到的如果为中最大值，其中， i=1,2,.,N，则对应的小波 p为第一层小波分解的最优小波； (e) 对第 k 层小波分解，其中， k 1，令：为最优小波 p在尺度 k-1 上的近似系数，从小波库中依次选择小波 i对新的信号 s 进行一。

6、层小波分解，得到近似系数 (f) 依次计算小波 i对应的近似系数的能量百分比直到 i=N ； (g)寻找出中最大的能量百分比值，其中， i=1,2,.,N，其对应的小波i为第j层小波分解的最优小波； (h) 重复步骤 (e) (g), 直到 k=K，从而得到所有分解层数上的最优基小波； (i) 将确定的各层最优基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数三个方向上的细节小波系数。 3. 根据权利要求 2 所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，所述的步骤 (i) 中，所述近似小波系数的能量百分比的计算公式为：权。

7、利要求书 CN 103839239 A 2 2/2 页 3 式中， k=1,2,K，为近似系数，为细节系数，其中 l 表示方向， l=1，代表水平方向； l=2，代表竖直方向； l=3，代表对角线方向， K 为小波分解层数。 4. 根据权利要求 1 所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，所述步骤 (4) 中，所述图像的细节小波系数的 Bayes 最优估计表达式为：式中， l 表示方向， l=1，代表水平方向； l=2，代表竖直方向； l=3，代表对角线方向， k 代表分解层数； (i,j) 代表小波系数的序号；代表第 k 层上，。

8、 l 方向上信号的小波系数方差；代表第 k 层上， l 方向上噪声的小波系数方差；代表第 k 层上， l 方向上带噪信号的小波系数；代表第 k 层上， l 方向上信号的小波估计。 5. 根据权利要求 4 所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，第所述的 k 层上， l 方向上的噪声估计均方差计算公式为：其中，为 k 层上， l 方向上带噪信号的小波系数中值。 6. 根据权利要求 4 所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，第所述的 k 层上， l 方向上的带噪信号的小波系数方差计算公式为：其中，为 k 层上， l 方向上带噪信号的小波系。

9、数的平均值。 7. 根据权利要求 4 所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，第所述的 k 层上， l 方向上图像信号的小波系数估计方差计算公式为：其中，为第 k 层上， l 方向上图像信号的小波系数估计方差，为第 k 层上， l 方向上的带噪信号的小波系数方差，为第 k 层上， l 方向上的噪声估计均方差。权利要求书 CN 103839239 A 3 1/6 页 4 一种电缆瓷套终端红外图像自适应去噪方法技术领域 0001 本发明涉及一种图像处理和在线监测技术，特别涉及一种电缆瓷套终端红外图像自适应去噪方法。背景技术 0002 基于红外图像的红外诊断。

10、技术已经被成功地运用于电气设备的状态检测和故障监测，并在发电机故障诊断、设备绝缘诊断等方面取得显著成效。电缆瓷套终端红外图像检测方法由于具有安全、准确、经济和实用等特点，近年来受到广泛关注。 0003 电缆瓷套终端红外图像是否清晰，图像细节是否完整，特征是否明显，是能否对其进行正确诊断的关键。红外图像在生成过程中受探测器噪声和各种电子器件噪声的影响，具有高噪声、低对比度的特点。因此采用高效率的滤波方法对电缆瓷套终端红外图像去噪，提升图像质量，是进行正确诊断的前提。 0004 针对红外图像的去噪，目前有很多算法，其中小波分析法由于其良好的时频局部化能力，。

11、具有较好的去噪效果。在小波去噪法中，基小波的选择对信号的畸变有着密切联系。常用的小波去噪方法是小波阈值法，即选定某个阈值对各尺度空间的小波系数进行截断处理，去除小于阈值的小波系数，保留大于阈值的小波系数并进行小波逆变换得到去噪后图像。小波阈值法虽然比常用的一些去噪方法对图像的模糊程度要小，但由于其是对小波系数进行截断式处理，这在一定程度上破坏了图像信息的完整性，造成部分有用信号的丢失。 0005 本发明要解决的技术问题是，克服现有技术中的不足，提供一种电缆瓷套终端红外图像自适应去噪方法，解决中压电缆电缆瓷套终端红外图像含有噪声的问题。该方法在自适应逐层选择最优。

12、基小波并对红外图像进行二维小波分解基础上，引入基于 Bayes 估计对细节小波系数进行处理。该方法能更有效地去除噪声，并完整保留图像细节。发明内容 0006 本发明的目的在于克服现有技术的缺点与不足，提供一种电缆瓷套终端红外图像自适应去噪方法，该自适应去噪方法有效解决了中压电缆电缆瓷套终端红外图像含有噪声的问题，并在自适应逐层选择最优基小波并对红外图像进行二维小波分解基础上，引入基于 Bayes 估计对细节小波系数进行处理，使该方法能更有效地去除噪声，并完整保留图像细节。 0007 本发明的目的通过下述技术方案实现：一种电缆瓷套终端红外图像自适应去噪方法，包。

13、括以下步骤： 0008 （1）输入待去噪的电缆瓷套终端红外图像； 0009 （2）将图像分解为将图像分解为 R 颜色分量图像、 G 颜色分量图像和 B 颜色分量图像； 0010 （3）对于每个颜色分量图像，采用自适应小波选择算法逐层确定最优基小波并将说明书 CN 103839239 A 4 2/6 页 5 确定的各层基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数； 0011 （4）对于每个颜色分量图像，保留近似小波系数，对不同分解层数上的三个方向上的细节小波系数进行自适应 Bayes 估计，得。

14、到图像的细节小波系数的 Bayes 最优估计； 0012 （5）对于每个颜色分量图像，利用最高分解层数上的近似小波系数和不同分解层数上的细节小波系数的 Bayes 最优估计进行图像重构，得到去噪后的图像； 0013 （6）利用去噪后的 R、 G、 B3 个颜色分量图像合成为原来的彩色图像，得到去噪后的结果图像。 0014 本发明中，所述步骤 (3) 包括以下步骤： 0015 (a) 给定一个小波函数库 i:i=1,2,.,N 和小波分解层数 K ； 0016 (b) 从小波函数库中依次选择小波 i，对红外图像信号 s 进行第一层小波分解，得到近似系数和细节系数l=1,。

15、2,3 ； 0017 (c) 依次计算小波 i对应的近似系数的能量百分比直到 i=N ； 0018 (d) 比较计算得到的如果为中最大值 (i=1,2,.,N)，则对应的小波 p为第一层小波分解的最优小波； 0019 (e) 对第 k(k 1) 层小波分解，令为最优小波 p在尺度 k-1 上的近似系数。从小波库中依次选择小波i对新的信号s进行一层小波分解，得到近似系数 0020 (f) 依次计算小波 i对应的近似系数的能量百分比直到 i=N ； 0021 (g) 寻找出中最大的能量百分比值 (i=1,2,.,N)，其对应的小波 i为第 j 层小波分解的最优小波； 0022 (。

16、h) 重复步骤 (e) (g), 直到 k=K，从而得到所有分解层数上的最优基小波； 0023 (i) 将确定的各层最优基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数。 0024 本发明中，所述的步骤 (i) 中，所述近似小波系数的能量百分比计算公式为： 0025 0026 式中， k=1,2,.,K，为近似系数，为细节系数，其中 l 表示方向。l=1，代表水平方向； l=2，代表竖直方向； l=3，代表对角线方向。K 为小波分解层数。 0027 本发明中，所述步骤 (4) 中，所述图像的细节小波系数。

17、的 Bayes 最优估计表达式为： 0028 0029 式中， l 表示方向， l=1，代表水平方向； l=2，代表竖直方向； l=3，代表对角线方向，说明书 CN 103839239 A 5 3/6 页 6 k 代表分解层数； (i,j) 代表小波系数的序号；代表第 k 层上， l 方向上信号的小波系数方差；代表第 k 层上， l 方向上噪声的小波系数方差；代表第 k 层上， l 方向上带噪信号的小波系数；代表第 k 层上， l 方向上信号的小波估计。 0030 本发明中，第所述的 k 层上， l 方向上的噪声估计均方差计算公式为： 0031 0032 。

18、其中，为 k 层上， l 方向上带噪信号的小波系数中值。 0033 本发明中，第所述的 k 层上， l 方向上的带噪信号的小波系数方差计算公式为： 0034 0035 其中，为 k 层上， l 方向上带噪信号的小波系数的平均值。 0036 本发明中，第所述的 k 层上， l 方向上图像信号的小波系数估计方差计算公式为： 0037 0038 其中，为第 k 层上， l 方向上图像信号的小波系数估计方差，为第 k 层上， l 方向上的带噪信号的小波系数方差，为第 k 层上， l 方向上的噪声估计均方差。 0039 本发明相对于现有技术具有如下的优点及效果： 0040 1、本发明的。

19、基于小波分解尺度系数能量最大原则的基小波选择方法具有自适应性，可以在不知红外图像先验信息的情况下，自适应逐层选择每一尺度上的能量最优基小波。 0041 2、本发明运用 Bayes 最优估计准则对带噪声的电缆瓷套终端红外图像在各个尺度空间，各个方向的小波系数进行估计，并利用估计的小波系数进行重构得到去噪后的图像。该方法在对电缆瓷套终端红外图像去噪的同时，能最大程度保留图像细节。附图说明 0042 图 1 为本发明的一种电缆瓷套终端红外图像自适应去噪方法的总体框图。 0043 图 2 为基于每层小波分解尺度系数能量最大原则，逐层自适应选择最优基小波步骤图。 0044 图。

20、3 为含噪的电缆瓷套终端红外图像。 0045 图 4 为用本发明方法去噪后得到的电缆瓷套终端红外图像。具体实施方式说明书 CN 103839239 A 6 4/6 页 7 0046 下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。 0047 实施例 0048 一种电缆瓷套终端红外图像自适应去噪方法，包括以下步骤： 0049 1. 将真彩色电缆瓷套终端图像分解为 R、 G、 B3 个颜色分量图像。 0050 2. 给定小波函数库为 sym1-sym30, 小波分解层数为 3 0051 3. 对于每个颜色分量图像，采用自适应小波选择算法逐层确定最优基小。

21、波并将确定的各层基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数。 0052 4.对于每个颜色分量图像，保留步3得到的近似小波系数，对步3得到的三个方向上的细节小波系数进行自适应 Bayes 估计，得到图像的细节小波系数的 Bayes 最优估计。 0053 5. 对于每个颜色分量图像，利用步 3 得到的最高分解层数上的近似小波系数和步 4 得到的不同分解层数上的细节小波系数的 Bayes 最优估计进行图像重构，得到去噪后的图像。 0054 6. 利用步骤 5 得到的去噪后的 R、 G、 B3 个颜色分量图像合成为。

22、原来的彩色图像，得到去噪后的结果图像。 0055 下面结合附图对本发明作进一步的详细说明。 0056 本发明提出一种电缆瓷套终端红外图像自适应去噪方法，该方法首先将真彩色红外图像分解为 R、 G、 B 颜色分量图像。然后基于每层小波分解尺度系数能量最大原则，逐层自适应选择最优基小波并分别对每个颜色分量图像进行二维小波分解。基于 Bayes 估计对细节小波系数进行处理，重构新的细节小波系数分量和最高尺度的低频小波系数分量获得去噪后的 R、 G、 B 颜色分量图像。将去噪后的 R、 G、 B 颜色分量图像合成为原来的彩色图像，得到结果的去噪图像。该方法不仅去除了电缆瓷套终端红外图。

23、像中的噪声，还完整保留了图像细节。 0057 如图 1 所示，为本发明的一种电缆瓷套终端红外图像自适应去噪方法的总体框图。 0058 1. 将真彩色电缆瓷套终端红外图像分解为 R、 G、 B 这 3 个颜色分量图像； 0059 2. 逐层自适应选择各个颜色分量图像最优小波基； 0060 如图 2 所示，为基于每层小波分解尺度系数能量最大原则，逐层自适应选择最优基小波步骤图； 0061 给定小波函数库为 sym1-sym30, 小波分解层数为 3， 0062 从小波函数库中依次选择小波 i，分别对每个颜色分量图像进行第一层小波分解，得到近似系数和细节系数l=1,2,3，。

24、依次计算小波 i对应的近似系数的能量百分比近似系数的能量百分比计算公式为： 0063 说明书 CN 103839239 A 7 5/6 页 8 0064 其中， k=1,2,K，为近似系数，为细节系数，其中 l 表示方向。l=1，代表水平方向； l=2，代表竖直方向； l=3，代表对角线方向。K 为小波分解层数。 0065 比较计算得到的如果为中最大值(i=1,2,.,N)，则对应的小波 p为第一层小波分解的最优小波。 0066 对第 k(k 1) 层小波分解，令为最优小波 p在尺度 k-1 上的近似系数。从小波库中依次选择小波 i对新的信号 s 进行一层小波分解，。

25、得到近似系数 0067 依次计算小波 i对应的近似系数的能量百分比直到 i=N。 0068 寻找出中最大的能量百分比值 (i=1,2,.,N)，其对应的小波 i为第 j 层小波分解的最优小波。 0069 得到所有分解层数上的最优基小波。 0070 对每个颜色分量图像，用确定的各层最优基小波对图像进行二维小波分解，得到最高分解层数上的近似小波系数和不同分解层数上的三个方向上的细节小波系数。 0071 3. 对每个颜色分量图像，采用基于 Bayes 估计对不同分解层数上的三个方向上的细节小波系数进行处理。图像的细节小波系数的 Bayes 最优估计表达式为： 0072 0073 其。

26、中， l 表示方向， l=1，代表水平方向； l=2，代表竖直方向； l=3，代表对角线方向， k 代表分解层数； (i,j) 代表小波系数的序号；代表第 k 层上， l 方向上信号的小波系数方差；代表第 k 层上， l 方向上噪声的小波系数方差；代表第 k 层上， l 方向上带噪信号的小波系数；代表第 k 层上， l 方向上信号的小波估计。 0074 第 k 层上， l 方向上的噪声估计均方差计算公式为： 0075 0076 其中，为 k 层上， l 方向上带噪信号的小波系数中值。 0077 第 k 层上， l 方向上的带噪信号的小波系数方差计算公式为： 0078。

27、 0079 其中，为 k 层上， l 方向上带噪信号的小波系数的平均值。 0080 第 k 层上， l 方向上图像信号的小波系数估计方差计算公式为： 0081 说明书 CN 103839239 A 8 6/6 页 9 0082 其中，为第k层上， l方向上图像信号的小波系数估计方差，为第k层上， l 方向上的带噪信号的小波系数方差，为第 k 层上， l 方向上的噪声估计均方差。 0083 4. 对于每个颜色分量图像，利用最高分解层数上的近似小波系数和不同分解层数上的细节小波系数的 Bayes 最优估计进行图像重构，得到去噪后的图像。 0084 5. 利用去噪后的 R、 G、 B3。

28、个颜色分量图像合成为原来的彩色图像，得到去噪后的结果图像。 0085 如图 3 所示，为含噪的电缆瓷套终端红外图像，如图 4 所示，为本发明的方法去噪后得到的电缆瓷套终端红外图像，经对比，本发明采用的方法能有效地去除噪声，并且能最大程度保留图像细节。 0086 上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。说明书 CN 103839239 A 9 1/3 页 10 图 1 说明书附图 CN 103839239 A 10 2/3 页 11 图 2 说明书附图 CN 103839239 A 11 3/3 页 12 图 3 图 4 说明书附图 CN 103839239 A 12 。

摘要
申请专利号：	CN201410071698.6	申请日：	2014.02.28
公开号：	CN103839239A	公开日：	2014.06.04
当前法律状态：	驳回	有效性：	无权
法律详情：	发明专利申请公布后的驳回IPC(主分类):G06T 5/00申请公布日:20140604\|\|\|实质审查的生效IPC(主分类):G06T 5/00申请日:20140228\|\|\|公开
IPC分类号：	G06T5/00	主分类号：	G06T5/00
申请人：	华南理工大学
发明人：	吴炬卓; 牛海清; 吴倩; 叶开发
地址：	510640 广东省广州市天河区五山路381号
优先权：
专利代理机构：	广州市华学知识产权代理有限公司 44245	代理人：	蔡茂略
PDF完整版下载：	PDF下载

内容摘要

本发明公开了一种电缆瓷套终端红外图像自适应去噪方法，包括以下步骤：(1)输入待去噪的电缆瓷套终端红外图像；(2)将图像分解为R颜色分量图像、G颜色分量图像和B颜色分量图像；(3)得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数；(4)得到图像的细节小波系数的Bayes最优估计；(5)得到去噪后的图像；(6)得到去噪后的结果图像。具有可以在不知红外图像先验信息的情况下，自适应逐层选择每一尺度上的能量最优基小波等优点。

权利要求书

权利要求书
1.  一种电缆瓷套终端红外图像自适应去噪方法，其特征在于，包括以下步骤：
(1)输入待去噪的电缆瓷套终端红外图像；
(2)将图像分解为R颜色分量图像、G颜色分量图像和B颜色分量图像；
(3)对于R颜色分量图像、G颜色分量图像和B颜色分量图像均采用自适应小波选择算法逐层确定最优基小波并将确定的各层基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数；
(4)对于R颜色分量图像、G颜色分量图像和B颜色分量图像均保留步骤(2)得到的近似小波系数，对步骤(2)得到的细节小波系数进行自适应Bayes估计，得到图像的细节小波系数的Bayes最优估计；
(5)对于每个颜色分量图像，利用步骤(2)得到的最高分解层数上的近似小波系数和步骤(4)得到的不同分解层数上的细节小波系数的Bayes最优估计进行图像重构，得到去噪后的图像；
(6)利用步骤(5)得到的去噪后的R颜色分量图像、G颜色分量图像和B颜色分量图像合成为原来的彩色图像，得到去噪后的结果图像。

2.  根据权利要求1所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，所述步骤(3)包括以下步骤：
(a)给定一个小波函数库{ψi:i=1,2,...,N}和小波分解层数K；
(b)从小波函数库中依次选择小波ψi，对红外图像信号s进行第一层小波分解，得到近似系数和细节系数l=1,2,3；
(c)依次计算小波ψi对应的近似系数的能量百分比直到i=N；
(d)比较计算得到的如果为中最大值，其中，i=1,2,...,N，则对应的小波ψp为第一层小波分解的最优小波；
(e)对第k层小波分解，其中，k＞1，令：为最优小波ψp在尺度k-1上的近似系数，从小波库中依次选择小波ψi对新的信号s进行一层小波分解，得到近似系数
(f)依次计算小波ψi对应的近似系数的能量百分比直到i=N；
(g)寻找出中最大的能量百分比值，其中，i=1,2,...,N，其对应的小波ψi为第j层小波分解的最优小波；
(h)重复步骤(e)～(g),直到k=K，从而得到所有分解层数上的最优基小波；
(i)将确定的各层最优基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数三个方向上的细节小波系数。

3.  根据权利要求2所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，所述的步骤(i)中，所述近似小波系数的能量百分比的计算公式为：
Ea=ΣiΣj(a(i,j)K)2ΣiΣj(a(i,j)K)2+ΣlΣkΣiΣj(W(i,j)k(l))2,]]>
式中，k=1,2,…,K，为近似系数，为细节系数，其中l表示方向，l=1，代表水平方向；l=2，代表竖直方向；l=3，代表对角线方向，K为小波分解层数。

4.  根据权利要求1所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，所述步骤(4)中，所述图像的细节小波系数的Bayes最优估计表达式为：
W^S(i,j)k(l)=σS,k2(l)σS,k2(l)+σn,k2(l)W(i,j)k(l),]]>
式中，l表示方向，l=1，代表水平方向；l=2，代表竖直方向；l=3，代表对角线方向，k代表分解层数；(i,j)代表小波系数的序号；代表第k层上，l方向上信号的小波系数方差；代表第k层上，l方向上噪声的小波系数方差；代表第k层上，l方向上带噪信号的小波系数；代表第k层上，l方向上信号的小波估计。

5.  根据权利要求4所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，第所述的k层上，l方向上的噪声估计均方差计算公式为：
σ^n,k=median|(W(i,j)k(l))0.6745,]]>
其中，为k层上，l方向上带噪信号的小波系数中值。

6.  根据权利要求4所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，第所述的k层上，l方向上的带噪信号的小波系数方差计算公式为：
σk2(l)=E{[W(i,j)k(l)-E](W(i,j)k(l))2},]]>
其中，为k层上，l方向上带噪信号的小波系数的平均值。

7.  根据权利要求4所述的电缆瓷套终端红外图像自适应去噪方法，其特征在于，第所述的k层上，l方向上图像信号的小波系数估计方差计算公式为：
σ^S,k2(l)=σk2(l)-σ^n,k2(l),]]>
其中，为第k层上，l方向上图像信号的小波系数估计方差，为第k层上，l方向上的带噪信号的小波系数方差，为第k层上，l方向上的噪声估计均方差。

说明书

说明书一种电缆瓷套终端红外图像自适应去噪方法
技术领域
本发明涉及一种图像处理和在线监测技术，特别涉及一种电缆瓷套终端红外图像自适应去噪方法。
背景技术
基于红外图像的红外诊断技术已经被成功地运用于电气设备的状态检测和故障监测，并在发电机故障诊断、设备绝缘诊断等方面取得显著成效。电缆瓷套终端红外图像检测方法由于具有安全、准确、经济和实用等特点，近年来受到广泛关注。
电缆瓷套终端红外图像是否清晰，图像细节是否完整，特征是否明显，是能否对其进行正确诊断的关键。红外图像在生成过程中受探测器噪声和各种电子器件噪声的影响，具有高噪声、低对比度的特点。因此采用高效率的滤波方法对电缆瓷套终端红外图像去噪，提升图像质量，是进行正确诊断的前提。
针对红外图像的去噪，目前有很多算法，其中小波分析法由于其良好的时频局部化能力，具有较好的去噪效果。在小波去噪法中，基小波的选择对信号的畸变有着密切联系。常用的小波去噪方法是小波阈值法，即选定某个阈值对各尺度空间的小波系数进行截断处理，去除小于阈值的小波系数，保留大于阈值的小波系数并进行小波逆变换得到去噪后图像。小波阈值法虽然比常用的一些去噪方法对图像的模糊程度要小，但由于其是对小波系数进行截断式处理，这在一定程度上破坏了图像信息的完整性，造成部分有用信号的丢失。
本发明要解决的技术问题是，克服现有技术中的不足，提供一种电缆瓷套终端红外图像自适应去噪方法，解决中压电缆电缆瓷套终端红外图像含有噪声的问题。该方法在自适应逐层选择最优基小波并对红外图像进行二维小波分解基础上，引入基于Bayes估计对细节小波系数进行处理。该方法能更有效地去除噪声，并完整保留图像细节。
发明内容
本发明的目的在于克服现有技术的缺点与不足，提供一种电缆瓷套终端红外图像自适应去噪方法，该自适应去噪方法有效解决了中压电缆电缆瓷套终端红外图像含有噪声的问题，并在自适应逐层选择最优基小波并对红外图像进行二维小波分解基础上，引入基于Bayes估计对细节小波系数进行处理，使该方法能更有效地去除噪声，并完整保留图像细节。
本发明的目的通过下述技术方案实现：一种电缆瓷套终端红外图像自适应去噪方法，包括以下步骤：
（1）输入待去噪的电缆瓷套终端红外图像；
（2）将图像分解为将图像分解为R颜色分量图像、G颜色分量图像和B颜色分量图像；
（3）对于每个颜色分量图像，采用自适应小波选择算法逐层确定最优基小波并将确定的各层基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数；
（4）对于每个颜色分量图像，保留近似小波系数，对不同分解层数上的三个方向上的细节小波系数进行自适应Bayes估计，得到图像的细节小波系数的Bayes最优估计；
（5）对于每个颜色分量图像，利用最高分解层数上的近似小波系数和不同分解层数上的细节小波系数的Bayes最优估计进行图像重构，得到去噪后的图像；
（6）利用去噪后的R、G、B3个颜色分量图像合成为原来的彩色图像，得到去噪后的结果图像。
本发明中，所述步骤(3)包括以下步骤：
(a)给定一个小波函数库{ψi:i=1,2,...,N}和小波分解层数K；
(b)从小波函数库中依次选择小波ψi，对红外图像信号s进行第一层小波分解，得到近似系数和细节系数l=1,2,3；
(c)依次计算小波ψi对应的近似系数的能量百分比直到i=N；
(d)比较计算得到的如果为中最大值(i=1,2,...,N)，则对应的小波ψp为第一层小波分解的最优小波；
(e)对第k(k＞1)层小波分解，令为最优小波ψp在尺度k-1上的近似系数。从小波库中依次选择小波ψi对新的信号s进行一层小波分解，得到近似系数
(f)依次计算小波ψi对应的近似系数的能量百分比直到i=N；
(g)寻找出中最大的能量百分比值(i=1,2,...,N)，其对应的小波ψi为第j层小波分解的最优小波；
(h)重复步骤(e)～(g),直到k=K，从而得到所有分解层数上的最优基小波；
(i)将确定的各层最优基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数。
本发明中，所述的步骤(i)中，所述近似小波系数的能量百分比计算公式为：
Ea=ΣiΣj(a(i,j)K)2ΣiΣj(a(i,j)K)2+ΣlΣkΣiΣj(W(i,j)k(l))2,]]>
式中，k=1,2,...,K，为近似系数，为细节系数，其中l表示方向。l=1，代表水平方向；l=2，代表竖直方向；l=3，代表对角线方向。K为小波分解层数。
本发明中，所述步骤(4)中，所述图像的细节小波系数的Bayes最优估计表达式为：
W^S(i,j)k(l)=σS,k2(l)σS,k2(l)+σn,k2(l)W(i,j)k(l),]]>
式中，l表示方向，l=1，代表水平方向；l=2，代表竖直方向；l=3，代表对角线方向，k代表分解层数；(i,j)代表小波系数的序号；代表第k层上，l方向上信号的小波系数方差；代表第k层上，l方向上噪声的小波系数方差；代表第k层上，l方向上带噪信号的小波系数；代表第k层上，l方向上信号的小波估计。
本发明中，第所述的k层上，l方向上的噪声估计均方差计算公式为：
σ^n,k=median|(W(i,j)k(l))0.6745,]]>
其中，为k层上，l方向上带噪信号的小波系数中值。
本发明中，第所述的k层上，l方向上的带噪信号的小波系数方差计算公式为：
σk2(l)=E{[W(i,j)k(l)-E](W(i,j)k(l))2},]]>
其中，为k层上，l方向上带噪信号的小波系数的平均值。
本发明中，第所述的k层上，l方向上图像信号的小波系数估计方差计算公式为：
σ^S,k2(l)=σk2(l)-σ^n,k2(l),]]>
其中，为第k层上，l方向上图像信号的小波系数估计方差，为第k层上，l方向上的带噪信号的小波系数方差，为第k层上，l方向上的噪声估计均方差。
本发明相对于现有技术具有如下的优点及效果：
1、本发明的基于小波分解尺度系数能量最大原则的基小波选择方法具有自适应性，可以在不知红外图像先验信息的情况下，自适应逐层选择每一尺度上的能量最优基小波。
2、本发明运用Bayes最优估计准则对带噪声的电缆瓷套终端红外图像在各个尺度空间，各个方向的小波系数进行估计，并利用估计的小波系数进行重构得到去噪后的图像。该方法在对电缆瓷套终端红外图像去噪的同时，能最大程度保留图像细节。
附图说明
图1为本发明的一种电缆瓷套终端红外图像自适应去噪方法的总体框图。
图2为基于每层小波分解尺度系数能量最大原则，逐层自适应选择最优基小波步骤图。
图3为含噪的电缆瓷套终端红外图像。
图4为用本发明方法去噪后得到的电缆瓷套终端红外图像。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
实施例
一种电缆瓷套终端红外图像自适应去噪方法，包括以下步骤：
1.将真彩色电缆瓷套终端图像分解为R、G、B3个颜色分量图像。
2.给定小波函数库为sym1-sym30,小波分解层数为3
3.对于每个颜色分量图像，采用自适应小波选择算法逐层确定最优基小波并将确定的各层基小波对图像逐层进行二维一次小波分解，得到最高分解层数上的近似小波系数和不同分解层数上三个方向上的细节小波系数。
4.对于每个颜色分量图像，保留步3得到的近似小波系数，对步3得到的三个方向上的细节小波系数进行自适应Bayes估计，得到图像的细节小波系数的Bayes最优估计。
5.对于每个颜色分量图像，利用步3得到的最高分解层数上的近似小波系数和步4得到的不同分解层数上的细节小波系数的Bayes最优估计进行图像重构，得到去噪后的图像。
6.利用步骤5得到的去噪后的R、G、B3个颜色分量图像合成为原来的彩色图像，得到去噪后的结果图像。
下面结合附图对本发明作进一步的详细说明。
本发明提出一种电缆瓷套终端红外图像自适应去噪方法，该方法首先将真彩色红外图像分解为R、G、B颜色分量图像。然后基于每层小波分解尺度系数能量最大原则，逐层自适应选择最优基小波并分别对每个颜色分量图像进行二维小波分解。基于Bayes估计对细节小波系数进行处理，重构新的细节小波系数分量和最高尺度的低频小波系数分量获得去噪后的R、G、B颜色分量图像。将去噪后的R、G、B颜色分量图像合成为原来的彩色图像，得到结果的去噪图像。该方法不仅去除了电缆瓷套终端红外图像中的噪声，还完整保留了图像细节。
如图1所示，为本发明的一种电缆瓷套终端红外图像自适应去噪方法的总体框图。
1.将真彩色电缆瓷套终端红外图像分解为R、G、B这3个颜色分量图像；
2.逐层自适应选择各个颜色分量图像最优小波基；
如图2所示，为基于每层小波分解尺度系数能量最大原则，逐层自适应选择最优基小波步骤图；
给定小波函数库为sym1-sym30,小波分解层数为3，
从小波函数库中依次选择小波ψi，分别对每个颜色分量图像进行第一层小波分解，得到近似系数和细节系数l=1,2,3，依次计算小波ψi对应的近似系数的能量百分比近似系数的能量百分比计算公式为：
Ea=ΣiΣj(a(i,j)K)2ΣiΣj(a(i,j)K)2+ΣlΣkΣiΣj(W(i,j)k(l))2,]]>
其中，k=1,2,…,K，为近似系数，为细节系数，其中l表示方向。l=1，代表水平方向；l=2，代表竖直方向；l=3，代表对角线方向。K为小波分解层数。
比较计算得到的如果为中最大值(i=1,2,...,N)，则对应的小波ψp为第一层小波分解的最优小波。
对第k(k＞1)层小波分解，令为最优小波ψp在尺度k-1上的近似系数。从小波库中依次选择小波ψi对新的信号s进行一层小波分解，得到近似系数
依次计算小波ψi对应的近似系数的能量百分比直到i=N。
寻找出中最大的能量百分比值(i=1,2,...,N)，其对应的小波ψi为第j层小波分解的最优小波。
得到所有分解层数上的最优基小波。
对每个颜色分量图像，用确定的各层最优基小波对图像进行二维小波分解，得到最高分解层数上的近似小波系数和不同分解层数上的三个方向上的细节小波系数。
3.对每个颜色分量图像，采用基于Bayes估计对不同分解层数上的三个方向上的细节小波系数进行处理。图像的细节小波系数的Bayes最优估计表达式为：
W^S(i,j)k(l)=σS,k2(l)σS,k2(l)+σn,k2(l)W(i,j)k(l),]]>
其中，l表示方向，l=1，代表水平方向；l=2，代表竖直方向；l=3，代表对角线方向，k代表分解层数；(i,j)代表小波系数的序号；代表第k层上， l方向上信号的小波系数方差；代表第k层上，l方向上噪声的小波系数方差；代表第k层上，l方向上带噪信号的小波系数；代表第k层上，l方向上信号的小波估计。
第k层上，l方向上的噪声估计均方差计算公式为：
σ^n,k=median|(W(i,j)k(l))0.6745,]]>
其中，为k层上，l方向上带噪信号的小波系数中值。
第k层上，l方向上的带噪信号的小波系数方差计算公式为：
σk2(l)=E{[W(i,j)k(l)-E](W(i,j)k(l))2},]]>
其中，为k层上，l方向上带噪信号的小波系数的平均值。
第k层上，l方向上图像信号的小波系数估计方差计算公式为：
σ^S,k2(l)=σk2(l)-σ^n,k2(l),]]>
其中，为第k层上，l方向上图像信号的小波系数估计方差，为第k层上，l方向上的带噪信号的小波系数方差，为第k层上，l方向上的噪声估计均方差。
4.对于每个颜色分量图像，利用最高分解层数上的近似小波系数和不同分解层数上的细节小波系数的Bayes最优估计进行图像重构，得到去噪后的图像。
5.利用去噪后的R、G、B3个颜色分量图像合成为原来的彩色图像，得到去噪后的结果图像。
如图3所示，为含噪的电缆瓷套终端红外图像，如图4所示，为本发明的方法去噪后得到的电缆瓷套终端红外图像，经对比，本发明采用的方法能有效地去除噪声，并且能最大程度保留图像细节。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。