广告位联系
返回顶部
分享到

基于Java实现空间滤波的代码

java 来源:转载 作者:秩名 发布时间:2021-08-24 08:01:22 人浏览
摘要

空间滤波的定义 滤波的本义是指信号有各种频率的成分,滤掉不想要的成分,即为滤掉常说的噪声,留下想要的成分,这即是滤波的过程,也是目的。空间滤波是一种采用滤波处理的影像增强方法。其理论基础是空间卷积和空间相关。目的是改善影像质量,包括去除高频噪声

空间滤波的定义

滤波的本义是指信号有各种频率的成分,滤掉不想要的成分,即为滤掉常说的噪声,留下想要的成分,这即是滤波的过程,也是目的。空间滤波是一种采用滤波处理的影像增强方法。其理论基础是空间卷积和空间相关。目的是改善影像质量,包括去除高频噪声与干扰,及影像边缘增强、线性增强以及去模糊等。分为低通滤波(平滑化)、高通滤波(锐化)和带通滤波。

图像需要增强的原因

各类图像处理系统在图像的采集、获取、传送和转换(如成像、复制扫描、传输以及显示等)过程中,均处在复杂的环境中,光照、电磁多变,所有的图像均不同程度地被可见或不可见的噪声干扰。噪声源包括电子噪声、光子噪声、斑点噪声和量化噪声。如果信噪比低于一定的水平,噪声逐渐变成可见的颗粒形状,导致图像质量的下降。除了视觉上质量下降,噪声同样可能掩盖重要的图像细节,在对采集到的原始图像做进一步的分割处理时,我们发现有一些分布不规律的椒盐噪声,为此采取相应的对策就是对图像进行必要的滤波降噪处理。

(1) 中值滤波

PART/01

是将每个像元在以其为中心的M×N邻域内取中间亮度值来代替该像元值,以达到去尖锐“噪声”和平滑图像的目的。具体计算方法与模板卷积方法类似,仍采用活动窗口的扫描方法。取值时,将M×N窗口内所有像元按亮度值的大小排列,取中间值作为中间像元的值。所以M×N取奇数为好。一般来说,图像亮度为阶梯状变化时,取均值平滑比取中值滤波要明显得多,而对于突出亮点的“噪声”干扰,从去“噪声”后对原图的保留程度看取中值要优于取均值。

代码实现:

比如我们要对下面矩阵进行中值滤波运算,

4 4 3 7 6 8 8

4 4 3 7 6 8 8

2 2 15 8 9 9 9

5 5 8 9 13 10 10

7 7 9 12 15 11 11

8 8 11 10 14 13 13

8 8 11 10 14 13 13

输出结果如下图所示:


 

package NB;

import java.util.ArrayList;

import java.util.Arrays;

import java.util.Scanner;

public class zzlb {

public static void main(String[] args) {

    Scanner s=new Scanner(System.in);

    System.out.println("请输入矩阵的行数:");

    int x=s.nextInt();//获取键盘输入的数字

    System.out.println("请输入矩阵的列数:");

    int y=s.nextInt();//获取键盘输入的数字

    int[][]b=new int[x][y];//创建一个二维数组

    int[]p=new int[9];//创建一个可以存放9个元素的一维数组来获取#3*3窗口的像元值

    System.out.println("请输入"+x+"*"+y+"的矩阵:");

    for (int i = 0; i <y; i++) {

        for (int j = 0; j < x; j++) {

            b[i][j]=s.nextInt();//将键盘输入的矩阵存放到二维数组里面

        }

    }

    System.out.print("运算结果如下:\n");

    for (int i = 1; i <y-1; i++) {

        for (int j = 1; j < x-1; j++) {  //锁定到中心像元的位置,从(1,1)开始

            int g=0;

            for (int k =i-1; k <=i+1 ; k++) {

                for (int l =j-1; l <=j+1 ; l++) {

                    p[g++]=b[k][l]; //将3*3矩阵窗口存储到一维数组中

                }

            }

            for (int k = 0; k <p.length-1; k++) {

                for (int l = 0; l < p.length-1-k; l++) {

                    if(p[l]>p[l+1]) {

                        int temp = p[l];

                        p[l]=p[l+1];

                        p[l+1]=temp;//本次采用冒泡排序法对3*3窗口内的像元从小到大排序

                    }

                }

            }

            System.out.print(p[4]+" ");//输出每个3*3模板的中心值,下标都是4

        }

        System.out.println("\n");//输出完一行就换行

    }

}
}

有小伙伴看到就会说:我求的是55的矩阵,为什么输入的是77的矩阵?原因在于我们运用的33矩阵窗口放在矩阵4个角的像元时,需要额外添加邻近的像元来构成33的窗口,对于这个添加的临时像元值一般和最近的像元值保持一致。

(2)罗伯特锐化

PART/02

图像锐化是为了突出图像上地物的边缘、轮廓,或某些线性目标要素的特征。这种滤波方法提高了地物边缘与周围像元之间的反差,因此也被称为边缘增强。锐化的方法很多,在此只介绍罗伯特梯度。梯度反映了相邻像元的亮度变化率,也就是说,图像中如果存在边缘,如湖泊、河流的边界,山脉和道路等,则边缘处有较大的梯度值。对于亮度值较平滑的部分,亮度梯度值较小。因此,找到梯度较大的位置,也就找到边缘,然后再用不同的梯度计算值代替边缘处像元的值,也就突出了边缘,实现了图像的锐化。

 

不过在讲解之前小编先给大家引入一个新的概念——图像卷积运算, 卷积运算:可看作是加权求和的过程,使图像区域中的每个像素分别与卷积核(权矩阵)的每个元素对应相乘,所有乘积之和作为区域中心像素的新值。它是在空间域上对图像作局部检测的运算,以实现平滑和锐化的目的。具体作法是选定一个卷积函数,又称为“模板”,实际上是一个M×N图像。二维的卷积运算是在图像中使用模板来实现运算的。运算方法从图像左上角开始,选定与模板同样大小的矩阵元素窗口,图像窗口与模板像元对应的亮度值相乘后再相加,最后一般将计算结果放在窗口中心位置(当M和N都是奇数时),代替原来的像元灰度值。然后活动窗口向右移动一个像元再以同样的方法进行卷积运算,逐行扫描,直到全幅影像都扫描一遍,最后生成新图像。罗伯特锐化方法使用的两个模板如下:

F=|aidxi|+|aidyi|,其中a是矩阵中的2*2个元素,i是第i个元素,i<=4,计算结果放在左上角像元,代替之前的灰度值,但是右下角的像元模板范围内没有其他像元了,无法进行计算,为此使用罗伯特锐化方法都要在原始矩阵的最下边和最右边添加一行和一列,值都和最临近的像元灰度值一样。为了直白了断,小编直接上例子吧,比如有某个矩阵,如下图(1)所示,我要用上面两个模板进行罗伯特锐化,首先得在右边和下边分别添加一行和一列,如下图(2)所示:

接下来对左上角像元进行运算:|2*(-1)+30+40+51|+|20+3*(-1)+41+50|=4,

接下来对右上角像元进行运算:|3*(-1)+30+50+51|+|30+3*(-1)+51+50|=4,

接下来对左下角像元进行运算:|4*(-1)+50+40+51|+|40+5*(-1)+41+50|=2,

 

接下来对右下角像元进行运算:|5*(-1)+50+50+51|+|50+5*(-1)+51+50|=0,最终输出矩阵如下所示:

这种算法的意义在于用交叉的方法检测出像与其领域在上下之间或左右之间或斜方向之间的差异,最终产生一个梯度影像,达到提取边缘信息的目的。

代码实现:

比如我们要对下面矩阵进行中值滤波运算,

2 2 10 10 10 10

2 2 10 10 10 10

2 2 10 10 10 10

2 2 2 2 2 2

2 2 2 2 2 2

2 2 2 2 2 2

输出结果如下图所示:

 
package NB;

import java.util.Scanner;

public class robet {
public static void main(String[] args) {

    Scanner s1 = new Scanner(System.in);

    System.out.println("请输入矩阵的行数:");

    int x = s1.nextInt();//获取键盘输入的数字

    System.out.println("请输入矩阵的列数:");

    int y = s1.nextInt();//获取键盘输入的数字

    int[][] b = new int[x][y];//创建一个二维数组

    System.out.println("请输入" + x + "*" + y + "的矩阵:");

    for (int i = 0; i < y; i++) {

        for (int j = 0; j < x; j++) {

            b[i][j] = s1.nextInt();

        }///将键盘输入的矩阵存放到二维数组里面

    }

    for (int i =0; i <y-1; i++) {

        for (int j = 0; j <x-1; j++) {

            b[i][j]=Math.abs(b[i][j]-b[i+1][j+1])+Math.abs(b[i+1][j]-b[i][j+1]);

            System.out.print(b[i][j]+" ");

        }//对两个模板进行图像卷积运算后取绝对值相加,并且把结果放到2*2窗口的左上方

        System.out.println("\n");//输出完一行就换行

    }

}
}


版权声明 : 本文内容来源于互联网或用户自行发布贡献,该文观点仅代表原作者本人。本站仅提供信息存储空间服务和不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权, 违法违规的内容, 请发送邮件至2530232025#qq.cn(#换@)举报,一经查实,本站将立刻删除。
原文链接 : https://blog.csdn.net/qq_43173805/article/details/119868866
相关文章
  • SpringBoot自定义错误处理逻辑介绍

    SpringBoot自定义错误处理逻辑介绍
    1. 自定义错误页面 将自定义错误页面放在 templates 的 error 文件夹下,SpringBoot 精确匹配错误信息,使用 4xx.html 或者 5xx.html 页面可以打印错误
  • Java实现手写一个线程池的代码

    Java实现手写一个线程池的代码
    线程池技术想必大家都不陌生把,相信在平时的工作中没有少用,而且这也是面试频率非常高的一个知识点,那么大家知道它的实现原理和
  • Java实现断点续传功能的代码

    Java实现断点续传功能的代码
    题目实现:网络资源的断点续传功能。 二、解题思路 获取要下载的资源网址 显示网络资源的大小 上次读取到的字节位置以及未读取的字节
  • 你可知HashMap为什么是线程不安全的
    HashMap 的线程不安全 HashMap 的线程不安全主要体现在下面两个方面 在 jdk 1.7 中,当并发执行扩容操作时会造成环形链和数据丢失的情况 在
  • ArrayList的动态扩容机制的介绍

    ArrayList的动态扩容机制的介绍
    对于 ArrayList 的动态扩容机制想必大家都听说过,之前的文章中也谈到过,不过由于时间久远,早已忘却。 所以利用这篇文章做做笔记,加
  • JVM基础之字节码的增强技术介绍

    JVM基础之字节码的增强技术介绍
    字节码增强技术 在上文中,着重介绍了字节码的结构,这为我们了解字节码增强技术的实现打下了基础。字节码增强技术就是一类对现有字
  • Java中的字节码增强技术

    Java中的字节码增强技术
    1.字节码增强技术 字节码增强技术就是一类对现有字节码进行修改或者动态生成全新字节码文件的技术。 参考地址 2.常见技术 技术分类 类
  • Redis BloomFilter布隆过滤器原理与实现

    Redis BloomFilter布隆过滤器原理与实现
    Bloom Filter 概念 布隆过滤器(英语:Bloom Filter)是1970年由一个叫布隆的小伙子提出的。它实际上是一个很长的二进制向量和一系列随机映射
  • Java C++算法题解leetcode801使序列递增的最小交换次

    Java C++算法题解leetcode801使序列递增的最小交换次
    题目要求 思路:状态机DP 实现一:状态机 Java 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 class Solution { public int minSwap(int[] nums1, int[] nums2) { int n
  • Mybatis结果集映射与生命周期介绍

    Mybatis结果集映射与生命周期介绍
    一、ResultMap结果集映射 1、设计思想 对简单的语句做到零配置,对于复杂一点的语句,只需要描述语句之间的关系就行了 2、resultMap的应用场
  • 本站所有内容来源于互联网或用户自行发布,本站仅提供信息存储空间服务,不拥有版权,不承担法律责任。如有侵犯您的权益,请您联系站长处理!
  • Copyright © 2017-2022 F11.CN All Rights Reserved. F11站长开发者网 版权所有 | 苏ICP备2022031554号-1 | 51LA统计