文章详情

短信预约-IT技能 免费直播动态提醒

请输入下面的图形验证码

提交验证

短信预约提醒成功

opencv中如何调整图像亮度和对比度,以及亮度和对比度调整的原理

2023-09-22 22:27

关注

文章目录

1. 调整图像亮度和对比度的原理

通过设计一个映射曲线就可完成对比度、亮度调整,具体过程如下图所示,
a) 为原图;
b) 把亮度调高,像素强度+固定值;
c) 把亮度调低,像素强度-固定值;
d) 增大像素强度75附近的对比度;增大附近的斜率
e) 增大像素强度150附近的对比度;
f) 增大像素强度75和225附近的对比度。

曲线斜率大于45度角的区域灰度被拉伸、精度上升、对比度变高;
曲线斜率小于45度角的区域灰度被压缩、精度下降、对比度变低。

下面图像是 0-255 映射到 0-255, a原图,b原图加上一个数,c原图减去一个数,d在value=75处拉伸,e在value=150出拉伸,f在75和225出进行拉伸。拉伸操作可以提升该区域的对比度
在这里插入图片描述

opencv关于对比度和亮度的误解

2. python实现的一个代码如下,可以直接运行

首先减去均值,然后调整对比度(乘上大于1的系数进行拉伸),然后加上均值,最后调整亮度。
aa是拉伸系数,bb是亮度变化
img_a = aa * (img - bri_mean) + bb + bri_mean

import cv2import numpy as npfrom matplotlib import pyplot as pltif __name__ == "__main__":    file = r'D:\dataset\data_gain1_2\cap_frame_0002_raw.png'    img = cv2.imread(file).astype(np.float32)    bri_mean = np.mean(img)    a = np.arange(5, 16, 5) / 10    b = np.arange(-30, 31, 30)    a_len = len(a)    b_len = len(b)    print(a_len, b_len)    plt.figure()    for i in range(a_len):        for j in range(b_len):            aa = a[i]            bb = b[j]            img_a = aa * (img-bri_mean) + bb + bri_mean            print(i, j, aa, bb)            img_a = np.clip(img_a,0,255).astype(np.uint8)            plt.subplot(a_len+1, b_len, (j + b_len * i + 1))            plt.imshow(img_a, cmap='gray')    plt.subplot(a_len + 1, b_len, a_len*b_len+1)    plt.imshow(img.astype(np.uint8), cmap='gray')    plt.show()

得到不同亮度和对比度的结果
在这里插入图片描述

3. 在opencv种调节图像和对比度的方法

可以参考下面链接的demo code:
opencv demo code
运行程序,对应直方图的变化:
原图
在这里插入图片描述
增加亮度,直方图向右平移
在这里插入图片描述
增加对比度, 直方图更加均衡
在这里插入图片描述

a. opencv中的实现方法:

首先认识 convertTo函数

image.convertTo(dst, CV_8U, a, b);

converTo 中的a,b是相乘和相加系数

opencv里是这样得到 a,b系数的在这里插入图片描述
opencv中时如何求 alpha和beta的?如下:

if (contrast > 0){    delta = 127. * contrast / 100;    a = 255. / (255. - delta * 2);    b = a * (brightness - delta);}else{    delta = -128. * contrast / 100;    a = (256. - delta * 2) / 255.;    b = a * brightness + delta;}

b. 根据咱们的公式,我们也可以得到自己的a,b系数,实现结果表明和opencv的方法具有类似的效果

img_a = aa * (img - bri_mean) + bb + bri_mean
bri_mean锚定为127的时候
img_a = aa * (img - 127) + bb + 127 = aa * img - 127*aa + bb + 127 = aa * img + (1-aa) * 127 + bb

opencv中_brightness 和 _contrast 的范围是 0-200, 等于100时不进行调整, 参考和运行下面的demo程序

a = _contrast / 100.0f;b = (1 - a) * 127 + _brightness - 100;

通过对比直方图的变化效果 可知上面两个方法 作用是类似的。
本来想找opencv的方法对应的原理链接,但是链接已经失效

c.两种方案的完整demo,c++实现:

#include "opencv2/core/utility.hpp"#include "opencv2/imgproc.hpp"#include "opencv2/imgcodecs.hpp"#include "opencv2/highgui.hpp"#include using namespace cv;using namespace std;int _brightness = 100;int _contrast = 100;Mat image;static void updateBrightnessContrast(int , void*){    int histSize = 64;    int brightness = _brightness - 100;    int contrast = _contrast - 100;        double a, b;        double delta=0;    //opencv中的方法    if (contrast > 0)    {        delta = 127. * contrast / 100;        a = 255. / (255. - delta * 2);        b = a * (brightness - delta);    }    else    {        delta = -128. * contrast / 100;        a = (256. - delta * 2) / 255.;        b = a * brightness + delta;    }         //a = _contrast / 100.0f;    //b = (1 - a) * 127 + _brightness - 100;    //printf("ab value :%.4lf, %.4lf, %.4lf\n",delta, a, b);    Mat dst, hist;    image.convertTo(dst, CV_8U, a, b);    imshow("image", dst);    calcHist(&dst, 1, 0, Mat(), hist, 1, &histSize, 0);    Mat histImage = Mat::ones(200, 320, CV_8U) * 255;    normalize(hist, hist, 0, histImage.rows, NORM_MINMAX, CV_32F);    histImage = Scalar::all(255);    int binW = cvRound((double)histImage.cols / histSize);    for (int i = 0; i < histSize; i++)        rectangle(histImage, Point(i * binW, histImage.rows),            Point((i + 1) * binW, histImage.rows - cvRound(hist.at<float>(i))),            Scalar::all(0), -1, 8, 0);    imshow("histogram", histImage);}const char* keys ={    "{help h||}{@image|baboon.jpg|input image file}"};int main(int argc, const char** argv){    CommandLineParser parser(argc, argv, keys);    parser.about("\nThis program demonstrates the use of calcHist() -- histogram creation.\n");    if (parser.has("help"))    {        parser.printMessage();        return 0;    }    string inputImage = parser.get<string>(0);    // Load the source image. HighGUI use.    image = imread(samples::findFile(inputImage), IMREAD_GRAYSCALE);    if (image.empty())    {        std::cerr << "Cannot read image file: " << inputImage << std::endl;        return -1;    }    namedWindow("image", 0);    namedWindow("histogram", 0);    createTrackbar("brightness", "image", &_brightness, 200, updateBrightnessContrast);    createTrackbar("contrast", "image", &_contrast, 200, updateBrightnessContrast);    updateBrightnessContrast(0, 0);    waitKey();    return 0;}

来源地址:https://blog.csdn.net/tywwwww/article/details/126626804

阅读原文内容投诉

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

软考中级精品资料免费领

  • 历年真题答案解析
  • 备考技巧名师总结
  • 高频考点精准押题
  • 2024年上半年信息系统项目管理师第二批次真题及答案解析(完整版)

    难度     813人已做
    查看
  • 【考后总结】2024年5月26日信息系统项目管理师第2批次考情分析

    难度     354人已做
    查看
  • 【考后总结】2024年5月25日信息系统项目管理师第1批次考情分析

    难度     318人已做
    查看
  • 2024年上半年软考高项第一、二批次真题考点汇总(完整版)

    难度     435人已做
    查看
  • 2024年上半年系统架构设计师考试综合知识真题

    难度     224人已做
    查看

相关文章

发现更多好内容

猜你喜欢

AI推送时光机
位置:首页-资讯-后端开发
咦!没有更多了?去看看其它编程学习网 内容吧
首页课程
资料下载
问答资讯