详析数字图像中高斯模糊理论及实现

高斯模糊(英语:Gaussian Blur),也叫高斯平滑。是在Adobe Photoshop、GIMP以及Paint.NET等图像处理软件中广泛使用的处理效果,通经常使用它来降低图像噪声以及降低细节层次。

简单介绍:高斯模糊(Gaussian Blur)是美国Adobe图像软件公司开发的一个图像处理软件:Adobe Photoshop(系列)中的一个滤镜,详细的位置在:滤镜—模糊——高斯模糊!

高斯模糊的原理中,它是依据高斯曲线调节像素色值,它是有选择地模糊图像。说得直白一点,就是高斯模糊能够把某一点周围的像素色值按高斯曲线统计起来,採用数学上加权平均的计算方法得到这条曲线的色值,最后能够留下人物的轮廓,即曲线.是指当 Adobe Photoshop 将加权平均应用于像素时生成的钟形曲线。

原理:

1 周边像素的平均值

所谓”模糊”。能够理解成每一个像素都取周边像素的平均值。

上图中,2是中间点。周边点都是1。

”中间点”取”周围点”的平均值。就会变成1。在数值上,这是一种”平滑化”。在图形上。就相当于产生”模糊”效果,”中间点”失去细节。

显然,计算平均值时,取值范围越大,”模糊效果”越强烈。

上图各自是原图、模糊半径3像素、模糊半径10像素的效果。模糊半径越大,图像就越模糊。从数值角度看,就是数值越平滑。

接下来的问题就是,既然每一个点都要取周边像素的平均值。那么应该怎样分配权重呢?

假设使用简单平均。显然不是非常合理,由于图像都是连续的,越靠近的点关系越密切,越远离的点关系越疏远。

因此,加权平均更合理,距离越近的点权重越大。距离越远的点权重越小。

2 正态分布权重

正态分布显然是一种可取的权重分配模式。

在图形上。正态分布是一种钟形曲线,越接近中心,取值越大,越远离中心,取值越小。计算平均值的时候,我们仅仅须要将”中心点”作为原点,其它点依照其在正态曲线上的位置,分配权重。就能够得到一个加权平均值。

3 高斯函数

上面的正态分布是一维的,图像都是二维的。所以我们须要二维的正态分布。

正态分布的密度函数叫做”高斯函数”(Gaussian function)。它的一维形式是:

当中,μ是x的均值,σ是x的方差。由于计算平均值的时候。中心点就是原公式进一步推导,由于计算平均值的时候,中心点就是原点,所以μ等于0。

依据一维高斯函数。能够推导得到二维高斯函数:

有了这个函数 ,就能够计算每一个点的权重了。

4 权重矩阵

假定中心点的坐标是(0,0),那么距离它近期的8个点的坐标例如以下,则权重前的矩阵为:

为了计算权重矩阵,须要设定σ的值。假定σ=1.5,则模糊半径为1的权重矩阵例如以下:

这9个点的权重总和等于0.4787147。假设仅仅计算这9个点的加权平均。还必须让它们的权重之和等于1,因此上面9个值还要分别除以0.4787147,得到终于的权重矩阵。

5 计算高斯模糊

有了权重矩阵,就能够计算高斯模糊的值了。

假设现有9个像素点。灰度值(0-255)例如以下:

每一个点乘以自己的权重值:

将这9个值加起来,就是中心点的高斯模糊的值。

对全部点反复这个过程,就得到了高斯模糊后的图像。假设原图是彩色图片,能够对RGB三个通道分别做高斯模糊。

6 实现代码

此代码为简单的高斯模糊生成代码,核心思想为取中心点像素值,依照步长为1 进行模糊。

(java代码)

package filter;
import java.awt.Color;
import java.awt.image.BufferedImage;
import java.io.File;
import java.io.IOException;
import javax.imageio.ImageIO;

public class GaussianBlur {
/*   简单的高斯模糊算法
  @param args
  @return void
  @see[类成员,类方法]*/
    public static void main(String[] args) throws IOException {
        //图片读取路径
        String pathname="d:\\my\\1.jpg";
        BufferedImage img = ImageIO.read(new File(pathname));
        System.out.println(img);
        int height =img.getHeight();
        int width = img.getWidth();
        int[][]matrix = new int [3][3];
        int []values  = new int [9];
        for (int i = 0; i < width; i++) {
            for (int j = 0; j < height; j++) {
                readPixel(img,i,j,values);//读取图像的像素值
                fillMatrix(matrix,values);//将图像的像素值填充到矩阵中
                img.setRGB(i, j, avgMatrix(matrix));//又一次依照平均值来生成新的img,将依照RGB模型整数像素数组设置为图像数据
            }
        }
        //图片保存路径
        String pathname1="d:\\my\\2.jpg";
        ImageIO.write(img, "jpeg", new File(pathname1));

    }
    //此函数的功能是将BuffedImage矩阵的像素值取出来
    public static void readPixel(BufferedImage img, int x, int y, int[] pixels) {
        int xStart = x-1;
        int yStart = y-1;
        int current = 0;
        for (int i = xStart; i <3+xStart ; i++) {
            for (int j = yStart; j <3+yStart;j++) {
                int tx=i;
                if (tx<0){
                    tx =-tx;
                }else if (tx>=img.getWidth()){
                    tx=x;
                }
                int ty=j;
                if (ty<0){
                    ty =-ty;
                }else if (ty>=img.getHeight()){
                    ty=y;
                }
            pixels[current++]= img.getRGB(tx, ty);
            }
            }
        }

//此函数功能为给定像素值。将像素值填充入矩阵中去
    public static void fillMatrix(int[][] matrix, int ... values) {
     int filled = 0;
     for (int i = 0; i < matrix.length; i++) {
        int []x =matrix[i]  ;
        for (int j = 0; j < x.length; j++) {
            x[j]=values[filled++];
        }
    }
  }
//此函数的功能是将像素的矩阵进行高斯模糊,又一次生成一个模糊后的矩阵
    public static int avgMatrix(int[][] matrix) {
        int r=0;
        int g=0;
        int b=0;
        for (int i = 0; i < matrix.length; i++) {
            int[] x =matrix[i];
            for (int j = 0; j < x.length; j++) {
                if (j==1) {
                    continue;
                }
            Color c = new Color(x[i]);
            r+= c.getRed();
            g+= c.getGreen();
            b+= c.getGreen();
        }
        }
        return new Color(r/8,g/8,b/8).getRGB();
  }
}

效果例如以下:高斯模糊前后对照图

我们再用MATLAB 进行实现一下(用MATLAB 处理图像既专业又方便,哈哈)

模糊又能够称为退化。我们能够用直接从库中调取。

无论是高斯模糊还是运动模糊非常easy实现地。。。

function ifbl = degrade(im, LEN, THETA, noisetype, M, V)
%Converting the image to grayshade by eliminating hue and
%saturation将图像转化为灰色来消除色相饱和度
if size(im, 3) == 3,
    im = rgb2gray(im);
end

%Converting to double
imd = im2double(im);

%Converting the image to frequency domain将图像转化到频率域中
imf = fft2(imd);

%Create PSF of degradation,滤波器卷积。逆时针角度THETA
PSF = fspecial(‘gussian‘,LEN,THETA);

%Convert psf to otf of desired size
%OTF is Optical Transfer Function
%fbl is blurred image in frequency domain,psf2otf是将点扩散函数转换成光学转换函数
OTF = psf2otf(PSF,size(im));

%Blurring the image
fbl = OTF.*imf;

%Converting back to spatial domainf返回空间域
ifbl = abs(ifft2(fbl));

%Checking if the image is in the range 0 to 1检查图像的范围
for i = 1:size(ifbl, 1)
    for j = 1:size(ifbl, 2)
        if ifbl(i, j) >= 1
            ifbl(i, j) = 0.999999;
        end
        if ifbl(i, j) <= 0
            ifbl(i, j) = 0.000001;
        end
    end
end        

%Adding noise
if nargin>3
    if nargin==4      %using default values
        ifbl = imnoise(ifbl, noisetype);
    elseif nargin==5  %specifying additional parameters explicitly
        ifbl = imnoise(ifbl, noisetype, M);
    elseif nargin==6  %specifying additional parameters explicitly
        ifbl = imnoise(ifbl, noisetype, M, V);
    end
end

高斯模糊。默认的高斯噪声均值m为0,方差var为0.01(參数能够自己设置)

好了,基本能够了,就是这个feel。。。

时间: 2024-10-29 19:11:35

详析数字图像中高斯模糊理论及实现的相关文章

Android 布局之LinearLayout 子控件weight权重的作用详析

关于Android开发中的LinearLayout子控件权重android:layout_weigh参数的作用,网上关于其用法有两种截然相反说法: 说法一:值越大,重要性越高,所占用的空间越大: 说法二:值越大,重要性越低,所占用的空间越小. 到底哪个正确?哪个错误?抑或还有其他解释?请点击查看关于weight 权重参数作用的详分析: 其实这两种情况都不太准确: 准确的解释是,weight 权限 是用于分配父控件某一方向上尺寸-所有子控件在该方向上设定尺寸和 所得值的一个参数,把这个相减得到的结

JS的预编译和执行顺序 详析

原文:JS的预编译和执行顺序 详析 最近在复习javascript的事件处理时发现了一个问题,然后也是我来写javascript的预编译和执行顺序的问题   代码:   复制代码 代码一 <html>   <head>     <title>事件处理</title>     <meta http-equiv="content-type" content="text/html;charset=utf-8"/>

利用QJM实现HDFS自动主从切换(HA Automatic Failover)源码详析

最近研究了下NameNode HA Automatic Failover方面的东西,当Active NN因为异常或其他原因不能正常提供服务时,处于Standby状态的NN就可以自动切换为Active状态,从而到达真正的高可用 NN HA Automatic Failover架构图 为了实现自动切换,需要依赖ZooKeeper和ZKFC组件,ZooKeeper主要用来记录NN的相关状态信息,zkfc组件以单独的JVM进程的形式运行在NN所在的节点上.下面首先分析下NN的启动流程,NN对象在实例化过

详解Android中那些酷炫返回方式的实现

Android手机都会有返回键,不管是实体键,还是虚拟键.Android用户主要也都是通过这个返回键操控页面返回方式的,不比IOS逼格甚高的只保留一个操作键.这种方式是最普遍的返回方式,还有一种也是比较常见的,那就是页面内部自己响应.绝大多数APP每个页面的设计图顶部左侧都会有一个返回键图标,偶尔也有奇葩的设计放在底部左侧,点击这个图标即finish掉当前页面.简单的介绍完了最常见的两种方式,下面为大家介绍两种更友好的交互方式. 拿大家比较常用的三款社交软件的交互来说.腾讯微博的返回方式除去上述

CSS学习笔记(9)--详解CSS中:nth-child的用法

详解CSS中:nth-child的用法 前端的哥们想必都接触过css中一个神奇的玩意,可以轻松选取你想要的标签并给与修改添加样式,是不是很给力,它就是":nth-child". 下面我将用几个典型的实例来给大家讲解:nth-child的实际用途: Tips:还用低版本的IE浏览器的哥们请绕过! :nth-child(2)选取第几个标签,"2可以是你想要的数字" .demo01 li:nth-child(2){background:#090} :nth-child(n

详解HTTP中get和post的区别

1.详解HTTP中GET和POST的区别 http://www.jellythink.com/archives/806 2.HTTP 方法:GET 对比 POST http://www.cnblogs.com/liu-ke/p/4198815.html

详解mysql中的Using与On的用法

多用才可以体会各个关键字的用法啊... 原文来自[http://bbs.php100.com/read-htm-tid-148469.html] 在用Join进行多表联合查询时,我们通常使用On来建立两个表的关系.其实还有一个更方便的关键字,那就是Using.那么这两个关键字在使用上有啥区别呢?往下看.假设有如下两张表:    mysql> select * from pets; +---------+---------+--------+-----------+ | pets_id | ani

详解Webwork中Action 调用的方法

详解Webwork中Action 调用的方法 从三方面介绍webwork action调用相关知识: 1.Webwork 获取和包装 web 参数 2.这部分框架类关系 3.DefaultActionProxyFactory.DefaultActionProxy.DefaultActionInvocation 终于要开始 webwork 核心业务类的总结,webwork 通过对客户端传递的 web 参数重新包装,进行执行业务 Action 类,并反馈执行结果,本篇源码分析对应下图 WebWork

实例详解 EJB 中的六大事务传播属性--转

前言 事务 (Transaction) 是访问并可能更新数据库中各种数据项的一个程序执行单元 (unit).在关系数据库中,一个事务可以是一条或一组 SQL 语句,甚至整个程序.它有通常被称为 ACID 的原子性(Atomicity).一致性(Consistency).隔离性(Isolation).持续性(Durability)四大特性: 原子性(Atomicity):一个事务是一个不可分割的工作单位,事务中包括的诸操作要么都做,要么都不做. 一致性(Consistency):事务必须是使数据库