AHD算法的matlab实现Bayer转RGB算法

Adaptive Homogeneity-Directed Demosaic Algorithm的matlab实现。

function [RGB] = MNartifact（RGBiteration）
% ****************************************************
% * Adaptive Homogeneity-Directed Demosaic Algorithm *
% ****************************************************
%   Designed by:  Keigo Hirakawa
%                 khirakawa1@udayton.edu
%
%   [RGB] = MNartifact（RGBiteration）
%         R         red image
%         G         green image
%         B         blue image
%         iteration number of iteration to run the G-R/B 
%                   correction stage
%
%   MNartifact reduces interpolation （color） artifacts
%
%   This algorithm was developed according to Hirakawa‘s master‘s 
%   thesis.
%

% adjust R/B then adjust G
for i=1:iteration
    R = G+median（cat（3conv2（R-G[1 0 0;0 0 0;0 0 0]‘same‘）conv2（R-G[0 1 0;0 0 0;0 0 0]‘same‘）conv2（R-G[0 0 1;0 0 0;0 0 0]‘same‘）conv2（R-G[0 0 0;1 0 0;0 0 0]‘same‘）conv2（R-G[0 0 0;0 0 1;0 0 0]‘same‘）conv2（R-G[0 0 0;0 0 0;1 0 0]‘same‘）conv2（R-G[0 0 0;0 0 0;0 1 0]‘same‘）conv2（R-G[0 0 0;0 0 0;0 0 1]‘same‘））3）;
    B = G+median（cat（3conv2（B-G[1 0 0;0 0 0;0 0 0]‘same‘）conv2（B-G[0 1 0;0 0 0;0 0 0]‘same‘）conv2（B-G[0 0 1;0 0 0;0 0 0]‘same‘）conv2（B-G[0 0 0;1 0 0;0 0 0]‘same‘）conv2（B-G[0 0 0;0 0 1;0 0 0]‘same‘）conv2（B-G[0 0 0;0 0 0;1 0 0]‘same‘）conv2（B-G[0 0 0;0 0 0;0 1 0]‘same‘）conv2（B-G[0 0 0;0 0 0;0 0 1]‘same‘））3）; 
    Gr = R+median（cat（3conv2（G-R[0 1 0;0 0 0;0 0 0]‘same‘）conv2（G-R[0 0 0;1 0 0;0 0 0]‘same‘）conv2（G-R[0 0 0;0 0 1;0 0 0]‘same‘）conv2（G-R[0 0 0;0 0 0;0 1 0]‘same‘））3）;
    Gb = B+median（cat（3conv2（G-B[0 1 0;0 0 0;0 0 0]‘same‘）conv2（G-B[0 0 0;1 0 0;0 0 0]‘same‘）conv2（G-B[0 0 0;0 0 1;0 0 0]‘same‘）conv2（G-B[0 0 0;0 0 0;0 1 0]‘same‘））3）;
    G = （Gr+Gb）/2;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3504  2012-03-16 11:46  MNdemosaic.m
     文件        1783  2012-03-16 11:47  MNartifact.m
     文件         988  2012-03-16 11:48  MNballset.m
     文件        1420  2012-03-16 11:47  MNdemosaicX.m
     文件         871  2012-03-16 11:47  MNdemosaicY.m
     文件        1465  2012-03-16 11:47  MNhomogeneity.m
     文件         506  2006-08-28 14:40  MNpad.m
     文件        1142  2012-03-16 11:48  MNparamA.m
     文件        1187  2012-03-16 11:48  MNrgb2lab.m
     文件         576  2002-11-05 21:41  mosaic.m
     文件        1819  2012-03-16 11:45  readme.txt