JAVA图像处理系列(四)——噪声

图像噪声是指存在于图像数据中的不必要的或多余的干扰信息。为图像添加噪声,常见有以下几种方法:
(1)高斯噪声
(2)脉冲噪声
(3)泊松噪声

高斯噪声

高斯噪声是指它的概率密度函数服从高斯分布(即正态分布)的一类噪声。高斯分布,也称正态分布,又称常态分布,记为N(μ,σ2),其中μ,σ2为分布的参数,分别为高斯分布的期望和方差。当有确定值时,p(x)也就确定了,特别当μ=0,σ^2=1时,X的分布为标准正态分布。
给图像添加告诉高斯噪声的代码如下:

import java.awt.image.BufferedImage;
import java.util.Random;

public class GaussianNoise {
    private Random rand = new Random(System.currentTimeMillis());
    
    private int generateNoise(int color) {
        double alpha, beta, sigma, value;
        alpha = rand.nextDouble();
        if (alpha == 0.0)
            alpha = 1.0;

        double SigmaGaussian = 4.0;
        double TauGaussian = 20.0;

        double tau;
        beta = rand.nextDouble();
        sigma = Math.sqrt(-2.0 * Math.log(alpha)) * Math.cos(2.0 * Math.PI * beta);
        tau = Math.sqrt(-2.0 * Math.log(alpha)) * Math.sin(2.0 * Math.PI * beta);
        value = (double) color + Math.sqrt((double) color) * SigmaGaussian * sigma + TauGaussian * tau;

        if (value < 0.0)
            return 0;
        if (value > 255)
            return (int) 255;
        return (int) (value + 0.5);
    }

    public BufferedImage addNoiseImage(BufferedImage image) {
        BufferedImage bimg = new BufferedImage(image.getWidth(),image.getHeight(),BufferedImage.TYPE_INT_RGB);
        Pixel pixel = new Pixel();

        for (int y = 0; y < image.getHeight(); y++) {
            for (int x = 0; x < image.getWidth(); x++) {
                pixel.setRGB(image.getRGB(x, y));
                pixel.red = generateNoise(pixel.red);
                pixel.green = generateNoise(pixel.green);
                pixel.blue = generateNoise(pixel.blue);
                bimg.setRGB(x, y, pixel.getRGB());
            }
        }

        return bimg;
    }
}

可以通过改变SigmaGaussian,TauGaussian的值,改变噪声的分布。

脉冲噪声

脉冲噪声(pulse noise)在通信中出现的离散型噪声的统称。它由时间上无规则出现的突发性干扰组成。脉冲噪声(impulsive noise)是非连续的,由持续时间短和幅度大的不规则脉冲或噪声尖峰组成。

给图像添加告诉脉冲噪声的代码如下:

import java.awt.image.BufferedImage;
import java.util.Random;

public class ImpulseNoise {
    private Random rand = new Random(System.currentTimeMillis());
    
    private int generateNoise(int color) {
        double alpha, value;
        alpha = rand.nextDouble();
        if (alpha == 0.0)
            alpha = 1.0;

        double SigmaImpulse = 0.1;

        if (alpha < (SigmaImpulse / 2.0)) {
            value = 0;
        } else if (alpha >= (1.0 - (SigmaImpulse / 2.0))) {
            value = 255;
        } else {
            value = color;
        }

        if (value < 0.0)
            return 0;
        if (value > 255)
            return (int) 255;
        return (int) (value + 0.5);
    }

    public BufferedImage addNoiseImage(BufferedImage image) {
        BufferedImage bimg = new BufferedImage(image.getWidth(),image.getHeight(),BufferedImage.TYPE_INT_RGB);
        Pixel pixel = new Pixel();

        for (int y = 0; y < image.getHeight(); y++) {
            for (int x = 0; x < image.getWidth(); x++) {
                pixel.setRGB(image.getRGB(x, y));
                pixel.red = generateNoise(pixel.red);
                pixel.green = generateNoise(pixel.green);
                pixel.blue = generateNoise(pixel.blue);
                bimg.setRGB(x, y, pixel.getRGB());
            }
        }

        return bimg;
    }
}

可以通过调整SigmaImpulse的值改变噪声分布

泊松噪声

泊松噪声,就是噪声分布符合泊松分布模型
泊松分布的概率函数为:

\mathbb{P}\left ( X= k \right )= \frac{\lambda ^{k } }{k!}e^{-\lambda },k=0,1,...

泊松分布的参数λ是单位时间(或单位面积)内随机事件的平均发生次数。 泊松分布适合于描述单位时间内随机事件发生的次数。

泊松分布的期望和方差均为

\lambda

特征函数为

\Psi \left ( t \right )= exp\left \{ \lambda \left ( e^{it} - 1\right ) \right \}

给图像添加泊松噪声代码如下:

import java.awt.image.BufferedImage;
import java.util.Random;

public class PoissonNoise {
    private Random rand = new Random(System.currentTimeMillis());
    
    private int generateNoise(int color) {
        double alpha, beta, value;
        alpha = rand.nextDouble();
        if (alpha == 0.0)
            alpha = 1.0;

        double SigmaPoisson = 0.05;

        int i = 0;
        for (i = 0; alpha > Math.exp(-SigmaPoisson * color); i++) {
            beta = rand.nextDouble();
            alpha = alpha * beta;
        }
        value = i / SigmaPoisson;

        if (value < 0.0)
            return 0;
        if (value > 255)
            return (int) 255;
        return (int) (value + 0.5);
    }

    public BufferedImage addNoiseImage(BufferedImage image) {
        BufferedImage bimg = new BufferedImage(image.getWidth(),image.getHeight(),BufferedImage.TYPE_INT_RGB);
        Pixel pixel = new Pixel();

        for (int y = 0; y < image.getHeight(); y++) {
            for (int x = 0; x < image.getWidth(); x++) {
                pixel.setRGB(image.getRGB(x, y));
                pixel.red = generateNoise(pixel.red);
                pixel.green = generateNoise(pixel.green);
                pixel.blue = generateNoise(pixel.blue);
                bimg.setRGB(x, y, pixel.getRGB());
            }
        }

        return bimg;
    }
}

可以通过改变SigmaPoisson的值,改变造成分布。

辅助类Pixel实现

辅助类Pixel代码实现如下:

import java.awt.image.*;

public class Pixel {
  public int red;
  public int green;
  public int blue;
  public int alpha=0xFF;
  public double hue;
  public double saturation;
  public double luminosity;
  private int rgb;

  public Pixel() {
  }

  public void setRGB(int rgb) {
    red = (rgb & 0x00FF0000) / 0x00010000;
    green = (rgb & 0x0000FF00) / 0x00000100;
    blue = rgb & 0x000000FF;
    alpha = (rgb >> 24)&0x0ff;
    this.rgb = rgb;
  }

  public int getRGB() {
    rgb =  alpha<<24 | 0x00010000 * red | 0x00000100 * green | blue;
    return this.rgb;
  }

  public static void setRgb(BufferedImage image, int x, int y, int red, int green, int blue) {
    int rgb = 0xFF000000 | red * 0x00010000 | green * 0x00000100 | blue;
    image.setRGB(x, y, rgb);
  }

  public static int pixelIntensity(int rgb) {
    int red = (rgb&0x00FF0000)/0x00010000;
    int green = (rgb&0x0000FF00)/0x00000100;
    int blue = rgb&0x000000FF;
    return (int) (0.299 * red + 0.587 * green + 0.114 * blue + 0.5);
  }
}

测试及效果

使用以下代码进行测试:

public class ImageEffectTest  {
   public static void main(String[] argv) throws IOException {
   BufferedImage img = read(new File("girl.jpg"));

   GaussianNoise gauss = new GaussianNoise();
   ImpulseNoise impulse = new ImpulseNoise();
   PoissonNoise poisson = new PoissonNoise();

   BufferedImage img2 = gauss.addNoiseImage(img);
   ImageIO.write(img2, "jpeg", new File("noise-girl-gaussian.jpg"));

   img2 = impulse.addNoiseImage(img);
   ImageIO.write(img2, "jpeg", new File("noise-girl-impulse.jpg"));

   img2 = poisson.addNoiseImage(img);
   ImageIO.write(img2, "jpeg", new File("noise-girl-poisson.jpg"));
 }
}

产生的图像效果如下:(第一张为原图,第二张为高斯噪声效果,第三张为脉冲噪声效果,第四张为泊松噪声效果)


girl3.jpg

noise-girl-gaussian.jpg

noise-girl-impulse.jpg

noise-girl-poisson.jpg
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 202,905评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,140评论 2 379
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 149,791评论 0 335
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,483评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,476评论 5 364
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,516评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,905评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,560评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,778评论 1 296
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,557评论 2 319
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,635评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,338评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,925评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,898评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,142评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,818评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,347评论 2 342

推荐阅读更多精彩内容