图像除雾序列——非参数图像增强5

5. 基于迭代的Retinex算法

参考文献[3 、 14] 和 http://www.cs.sfu.ca/~colour/publications/IST-2000/ ，图像增强效果如图 8 所示。

 

图8-(a) 原始图像

 

图8-(b) McCann99 增强效果

 

图8-(b) Frankle-McCann增强效果

Frankle-McCann算法的代码如下：

function Retinex = retinex_frankle_mccann(L, nIterations)

% RETINEX_FRANKLE_McCANN: 
%         Computes the raw Retinex output from an intensity image, based on the
%         original model described in:
%         Frankle, J. and McCann, J., "Method and Apparatus for Lightness Imaging"
%         US Patent #4,384,336, May 17, 1983
%
% INPUT:  L           - logarithmic single-channel intensity image to be processed
%         nIterations - number of Retinex iterations
%
% OUTPUT: Retinex     - raw Retinex output
%
% NOTES:  - The input image is assumed to be logarithmic and in the range [0..1]
%         - To obtain the retinex "sensation" prediction, a look-up-table needs to
%         be applied to the raw retinex output
%