图像噪声--添加噪声

椒盐噪声

        椒盐噪声就是给图片添加黑白噪点，椒指的是黑色的噪点(0,0,0),盐指的是白色的噪点(255,255,255)，通过num来控制噪声多少，值越大添加的噪声越多，图像损坏的更加严重。

void add_salt_pepper_noise(Mat& src,Mat& dst,int num)
{
    dst = src.clone();

    unsigned seed = std::chrono::system_clock::now().time_since_epoch().count();
    std::default_random_engine generator (seed);

    std::uniform_int_distribution distribute_x(0, src.cols - 1);//包含max
    std::uniform_int_distribution distribute_y(0, src.rows - 1);//包含max

    //black
    for(int i = 0;i < num;i ++){
        int x = distribute_x(generator);
        int y = distribute_y(generator);
        dst.at(y,x) = Vec3b(0,0,0);
    }

    //white
    for(int i = 0;i < num;i ++){
        int x = distribute_x(generator);
        int y = distribute_y(generator);
        dst.at(y,x) = Vec3b(255,255,255);
    }
}

下图是设置num=1000时产生的噪声图像： 

 

高斯噪声

        高斯噪声就是给图片添加一个服从高斯分布的噪声，可以通过调节高斯分布标准差sigma的大小来控制添加噪声程度，sigma越大添加的噪声越多图片损坏的越厉害.

void add_gassion_noise(Mat& src,Mat& dst,int mean,int sigma)
{
    dst = src.clone();

    Mat _gassion = Mat::zeros(src.rows,src.cols,CV_8UC3);
    cv::RNG rng;
    rng.fill(_gassion, RNG::NORMAL, mean, sigma);

    dst = dst + _gassion;
}

下图是mean=50,sigma=50时产生的噪声图像： 

随机形状斑点噪声

        随机斑点噪声是通过生成随机的多变形来实现的。

注意以下代码中image是QImage，num表示产生20个斑点。

    Mat _msk = Mat::zeros(image->height(),image->width(),CV_8UC1);
    num = 20;
    Mat kernal = Mat::ones(3,3,CV_8UC1);
    kernal.at(0,0) = 0;
    kernal.at(0,2) = 0;
    kernal.at(2,0) = 0;
    kernal.at(2,2) = 0;

    for(int i = 0;i < num;i ++){
        quint32 value = QRandomGenerator::global()->generate();
        int cx = value % (this->width() / 2 - 1);

        value = QRandomGenerator::global()->generate();
        int cy = value % (this->height() - 1);

        std::vector _polygon;
        int _line_num = 20,_min_radius = 2,_max_radius = 7;
        double _angle_step = 2 * 3.1415926 / _line_num;
        Point center(cx,cy);
        for (int i = 0;i < _line_num;i ++){
            quint32 value = QRandomGenerator::global()->generate();
            double k = value % 100 / 100.0;
            double _angle = _angle_step * i;
            _angle = _angle + (k - 0.5) * _angle_step * 0.25;

            value = QRandomGenerator::global()->generate();
            k = value % 100 / 100.0;
            double radius = _min_radius + k * (_max_radius - _min_radius);

            int x = radius * cos(_angle);
            int y = radius * sin(_angle);

            _polygon.push_back(Point(center.x + x,center.y + y));
        }

        value = QRandomGenerator::global()->generate();
        int c = value % 255;

        cv::fillPoly(_msk,_polygon,Scalar(c));

        //cv::morphologyEx(_msk,_msk,MORPH_DILATE,kernal,Point(-1,-1),1);
    }

    cv::GaussianBlur(_msk,_msk,Size(5,5),0);

    for (int h = 0;h < _msk.rows;h ++)
    {
        uchar *d1 = image->scanLine(h);
        uchar *d2 = _msk.ptr(h);
        for (int w = 0;w < _msk.cols;w ++)
        {
            int w4 = 4*w;

            double v = d1[w4 + 0] * (1 - d2[w] / 255.0);
            if (v > 255)v = 255;
            if (v < 0)v = 0;
            d1[w4 + 0] = v;

            v = d1[w4 + 1] * (1 - d2[w] / 255.0);
            if (v > 255)v = 255;
            if (v < 0)v = 0;
            d1[w4 + 1] = v;

            v = d1[w4 + 2] * (1 - d2[w] / 255.0);
            if (v > 255)v = 255;
            if (v < 0)v = 0;
            d1[w4 + 2] = v;
        }
    }

以上代码只产生黑色斑点， 白色斑点只要将合并方式改成加即可：double v = d1[w4 + 0] * (1 + d2[w] / 255.0);