当前位置:   article > 正文

OpenCV-白平衡(完美反射算法)_opencv 白平衡

opencv 白平衡

作者:翟天保Steven
版权声明:著作权归作者所有,商业转载请联系作者获得授权,非商业转载请注明出处

实现原理

       白平衡的意义在于,对在特定光源下拍摄时出现的偏色现象,通过加强对应的补色来进行补偿,使白色物体能还原为白色。

       完美反射算法是白平衡各种算法中较常见的一种,比灰度世界算法更优。它假设图像世界中最亮的白点是一个镜面,能完美反射光照;基于白点,将三通道的数值进行适当地调整,以达到白平衡效果;除此之外,还需要统计最亮的一定区间的三通道均值,该均值与该通道最大值的差距决定了该通道调整的力度。

       通俗的讲,若图像中绿色分量最大值是255,但是绿色最亮的前百分之10个点的平均值只有80,说明原图的绿色分量整体较低,需要对其加强;若最大值只有100,那么加强的系数就较低,白平衡的效果就不达预期。这就是完美反射算法比较依赖图像中存在白点的原因,白点的三通道灰度值接近【255,255,255】。最下方将用实际图像作进一步说明,以帮助读者理解。

       完美反射算法的实现流程如下:   

       1.计算图像RGB三通道各自的灰度最大值Rmax、Gmax、Bmax。

       2.利用三通道数值和,确定图像最亮区间的下限T。

       3.计算图像三通道数值和大于T的点的三通道均值Rm、Gm、Bm。

       4.计算三通道的补偿系数,即单通道最大值除以单通道亮区平均值。

功能函数代码

  1. // 白平衡-完美反射
  2. cv::Mat WhiteBalcane_PRA(cv::Mat src)
  3. {
  4. cv::Mat result = src.clone();
  5. if (src.channels() != 3)
  6. {
  7. cout << "The number of image channels is not 3." << endl;
  8. return result;
  9. }
  10. // 通道分离
  11. vector<cv::Mat> Channel;
  12. cv::split(src, Channel);
  13. // 定义参数
  14. int row = src.rows;
  15. int col = src.cols;
  16. int RGBSum[766] = { 0 };
  17. uchar maxR, maxG, maxB;
  18. // 计算单通道最大值
  19. for (int i = 0; i < row; ++i)
  20. {
  21. uchar *b = Channel[0].ptr<uchar>(i);
  22. uchar *g = Channel[1].ptr<uchar>(i);
  23. uchar *r = Channel[2].ptr<uchar>(i);
  24. for (int j = 0; j < col; ++j)
  25. {
  26. int sum = b[j] + g[j] + r[j];
  27. RGBSum[sum]++;
  28. maxB = max(maxB, b[j]);
  29. maxG = max(maxG, g[j]);
  30. maxR = max(maxR, r[j]);
  31. }
  32. }
  33. // 计算最亮区间下限T
  34. int T = 0;
  35. int num = 0;
  36. int K = static_cast<int>(row * col * 0.1);
  37. for (int i = 765; i >= 0; --i)
  38. {
  39. num += RGBSum[i];
  40. if (num > K)
  41. {
  42. T = i;
  43. break;
  44. }
  45. }
  46. // 计算单通道亮区平均值
  47. double Bm = 0.0, Gm = 0.0, Rm = 0.0;
  48. int count = 0;
  49. for (int i = 0; i < row; ++i)
  50. {
  51. uchar *b = Channel[0].ptr<uchar>(i);
  52. uchar *g = Channel[1].ptr<uchar>(i);
  53. uchar *r = Channel[2].ptr<uchar>(i);
  54. for (int j = 0; j < col; ++j)
  55. {
  56. int sum = b[j] + g[j] + r[j];
  57. if (sum > T)
  58. {
  59. Bm += b[j];
  60. Gm += g[j];
  61. Rm += r[j];
  62. count++;
  63. }
  64. }
  65. }
  66. Bm /= count;
  67. Gm /= count;
  68. Rm /= count;
  69. // 通道调整
  70. Channel[0] *= maxB / Bm;
  71. Channel[1] *= maxG / Gm;
  72. Channel[2] *= maxR / Rm;
  73. // 合并通道
  74. cv::merge(Channel, result);
  75. return result;
  76. }

C++测试代码

  1. #include <iostream>
  2. #include <opencv.hpp>
  3. using namespace std;
  4. // 白平衡-完美反射
  5. cv::Mat WhiteBalcane_PRA(cv::Mat src)
  6. {
  7. cv::Mat result = src.clone();
  8. if (src.channels() != 3)
  9. {
  10. cout << "The number of image channels is not 3." << endl;
  11. return result;
  12. }
  13. // 通道分离
  14. vector<cv::Mat> Channel;
  15. cv::split(src, Channel);
  16. // 定义参数
  17. int row = src.rows;
  18. int col = src.cols;
  19. int RGBSum[766] = { 0 };
  20. uchar maxR, maxG, maxB;
  21. // 计算单通道最大值
  22. for (int i = 0; i < row; ++i)
  23. {
  24. uchar *b = Channel[0].ptr<uchar>(i);
  25. uchar *g = Channel[1].ptr<uchar>(i);
  26. uchar *r = Channel[2].ptr<uchar>(i);
  27. for (int j = 0; j < col; ++j)
  28. {
  29. int sum = b[j] + g[j] + r[j];
  30. RGBSum[sum]++;
  31. maxB = max(maxB, b[j]);
  32. maxG = max(maxG, g[j]);
  33. maxR = max(maxR, r[j]);
  34. }
  35. }
  36. // 计算最亮区间下限T
  37. int T = 0;
  38. int num = 0;
  39. int K = static_cast<int>(row * col * 0.1);
  40. for (int i = 765; i >= 0; --i)
  41. {
  42. num += RGBSum[i];
  43. if (num > K)
  44. {
  45. T = i;
  46. break;
  47. }
  48. }
  49. // 计算单通道亮区平均值
  50. double Bm = 0.0, Gm = 0.0, Rm = 0.0;
  51. int count = 0;
  52. for (int i = 0; i < row; ++i)
  53. {
  54. uchar *b = Channel[0].ptr<uchar>(i);
  55. uchar *g = Channel[1].ptr<uchar>(i);
  56. uchar *r = Channel[2].ptr<uchar>(i);
  57. for (int j = 0; j < col; ++j)
  58. {
  59. int sum = b[j] + g[j] + r[j];
  60. if (sum > T)
  61. {
  62. Bm += b[j];
  63. Gm += g[j];
  64. Rm += r[j];
  65. count++;
  66. }
  67. }
  68. }
  69. Bm /= count;
  70. Gm /= count;
  71. Rm /= count;
  72. // 通道调整
  73. Channel[0] *= maxB / Bm;
  74. Channel[1] *= maxG / Gm;
  75. Channel[2] *= maxR / Rm;
  76. // 合并通道
  77. cv::merge(Channel, result);
  78. return result;
  79. }
  80. int main()
  81. {
  82. // 载入原图
  83. cv::Mat src = cv::imread("test21.jpg");
  84. // 白平衡-完美反射
  85. cv::Mat result = WhiteBalcane_PRA(src);
  86. // 显示
  87. cv::imshow("src", src);
  88. cv::imshow("result", result);
  89. cv::waitKey(0);
  90. return 0;
  91. }

测试效果

图1 原图
图2 白平衡后图像

       如图1所示,是傍晚的一张图像,众所周知,傍晚的色温是较低的,此时采用高于傍晚色温的色温值拍照,就会得到一张暖色系的图片,偏黄;对其进行白平衡,使图片颜色回归真实的环境色温,就得到如图2的效果。

       如果你用过灰度世界算法,你会发现完美反射算法的效果更亮些;这是因为蓝通道最大值和蓝通道亮区平均值的差距较大,因而补偿的强度很强;若原图中不存在白色区,那蓝通道的补偿就会较弱,达不到较好的预期,因此该算法所处理的图像中最好有较白的点。

图3 单色原图
图4 单色图白平衡效果

        如图3所示,是一张色彩相对一致的图像,整体呈粉色系。灰度世界法将三通道数值进行平均再补偿,就会使三通道的数值趋近一致,进而呈现灰色;而完美反射算法,计算的是三通道各自数值最大值和其亮区平均值的差距,再进行补偿,因此三通道的数值都会适当加强,使光感更强,即图片更亮。

        接下来做个有趣的测试,将原本粉色的墙纸设为较纯的绿色。

图5 调色后的图像
图6 白平衡效果图

        如图5所示,因为图像中存在色调相冲的两个部分,完美反射算法的优势在这种场景下体现的很明显。因为绿色区域较大,灰度世界算法中单纯的计算均值再平衡,会无脑地将整图的绿色分量降低,红色分量提高,使得结果异常滑稽;而完美反射算法,因为图中有白值,三通道的最大值均在250左右,对绿色分量而言,其最亮区的均值接近于230,而红色分量和蓝色分量而言,其最亮区的均值也接近于200多,因此三通道的平衡结果就是整体提亮,而不是被无脑平均。

        灰度世界算法效果图见:

OpenCV-白平衡(灰度世界算法)_翟天保Steven的博客-CSDN博客

       如果函数有什么可以改进完善的地方,非常欢迎大家指出,一同进步何乐而不为呢~

       如果文章帮助到你了,可以点个赞让我知道,我会很快乐~加油!

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/小小林熬夜学编程/article/detail/287477
推荐阅读
相关标签
  

闽ICP备14008679号