绿平衡（Green Balance）

ISPer_yy

于 2023-03-11 18:39:56 发布

阅读量1.5k

点赞数 1

分类专栏： ISP-Pipeline 文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/m0_73188795/article/details/129467006

版权

ISP-Pipeline 专栏收录该内容

5 篇文章 3 订阅

订阅专栏

图像信号处理芯片设计原理----07 绿平衡+去马赛克
 Green Imbalance
专利CN105530500的基本思路：
关键：在平衡Gr和Gb的同时，要尽量保持颜色边缘的清晰。

分别计算某范围内Gr和Gb的平均值，记为Gr_mean 和 Gb_mean，计算平均值的差值的一半，记为diff_mean:
$diff\_mean = \frac{|Gr\_mean-Gb\_mean|}{2}$
计算各个2x2块中Gr和Gb的差值：
$d i ff = ∣ G r - G b ∣$
判断diff_mean是否大于第一阈值thed1或者diff是否大于第二阈值thed2，二者满足其一则不修改Gr和Gb的值，否则执行下一步。
修改
$Gb)/2\\ Gb = (Gr + Gb)/2$

复现的MATLAB代码

function out_raw = GIC(img_raw,params,thed1,thed2)

H                   = [1 0 1 0 1; 0 0 0 0 0; 1 0 1 0 1; 0 0 0 0 0; 1 0 1 0 1]/9;
raw_mean            = imfilter(img_raw,H,"symmetric");
out_raw             = img_raw;
switch params.pattern
    case 'rggb'
        Gr          = img_raw(1:2:end,2:2:end);
        Gr_mean     = raw_mean(1:2:end,2:2:end);
        Gb          = img_raw(2:2:end,1:2:end);
        Gb_mean     = raw_mean(2:2:end,1:2:end);
    case 'bggr'
        Gr          = img_raw(2:2:end,1:2:end);
        Gr_mean     = raw_mean(2:2:end,1:2:end);
        Gb          = img_raw(1:2:end,2:2:end);
        Gb_mean     = raw_mean(1:2:end,2:2:end);
    case 'grbg'
        Gr          = img_raw(1:2:end,1:2:end);
        Gr_mean     = raw_mean(1:2:end,1:2:end);
        Gb          = img_raw(2:2:end,2:2:end);
        Gb_mean     = raw_mean(2:2:end,2:2:end);
    case 'gbrg'
        Gr          = img_raw(2:2:end,2:2:end);
        Gr_mean     = raw_mean(2:2:end,2:2:end);
        Gb          = img_raw(1:2:end,2:2:end);
        Gb_mean     = raw_mean(1:2:end,2:2:end);
end

mean_diff           = abs(Gr_mean - Gb_mean) / 2;
diff                = abs(Gr - Gb);
flag1               = (mean_diff > thed1);
flag2               = (diff > thed2);
flag                = min(1,flag1+flag2);

% Gr_out              = max(0,min(255,Gr + mean_diff));
% Gb_out              = max(0,min(255,Gb - mean_diff));

Gr_out              = max(0,min(255,(Gr + Gb)/2));
Gb_out              = max(0,min(255,(Gr + Gb)/2));

Gr_out              = Gr.*flag + Gr_out.*(1-flag);
Gb_out              = Gb.*flag + Gb_out.*(1-flag);

switch params.pattern
    case 'rggb'
        out_raw(1:2:end,2:2:end) = Gr_out;
        out_raw(2:2:end,1:2:end) = Gb_out;
    case 'bggr'
        out_raw(1:2:end,2:2:end) = Gb_out;
        out_raw(2:2:end,1:2:end) = Gr_out;
    case 'grbg'
        out_raw(1:2:end,1:2:end) = Gr_out;
        out_raw(2:2:end,2:2:end) = Gb_out;
    case 'gbrg'
        out_raw(1:2:end,1:2:end) = Gb_out;
        out_raw(2:2:end,2:2:end) = Gr_out;
end