1.参考文献
论文:Night video enhancement using improved dark channel prior
论文作者:Xuesong Jiang,Hongxun Yao,Shengping Zhang,Xiusheng Lu and Wei Zeng
DOI:10.1109/ICIP.2013.6738114
略有改动,只用于处理单张图像,去掉了金字塔放缩过程,有视频处理需要可根据原论文在本代码上进行简单修改即可
2.matlab代码
%论文:Night video enhancement using improved dark channel prior
%论文作者:Xuesong Jiang,Hongxun Yao,Shengping Zhang,Xiusheng Lu and Wei Zeng
%DOI:10.1109/ICIP.2013.6738114
%代码作者:slsl97
%略有改动,只用于处理单张图像,去掉了金字塔放缩过程,有视频处理需要可根据原论文在本代码上进行简单修改即可
clear;
clc;
Input=imread('Input.jpg');
w=0.98;
miu=0.95;
%%
%Calculate A
I_inv=255-Input;
I_min=min(I_inv,[],3);
I_min_med=medfilt2(I_min);
A=max(max(I_min_med));
%%
%Calculate Ac
I_inv_r=I_inv(:,:,1);
I_r_med=medfilt2(I_inv_r);
A_r=max(max(I_r_med));
I_inv_g=I_inv(:,:,2);
I_g_med=medfilt2(I_inv_g);
A_g=max(max(I_g_med));
I_inv_b=I_inv(:,:,3);
I_b_med=medfilt2(I_inv_b);
A_b=max(max(I_b_med));
I_inv_A(:,:,1)=double(I_inv_r)/double(A_r);
I_inv_A(:,:,2)=double(I_inv_r)/double(A_g);
I_inv_A(:,:,3)=double(I_inv_r)/double(A_b);
%%
I_dark_til=min(I_inv_A,[],3);
I_med=medfilt2(I_dark_til);
I_detail=medfilt2(abs(I_med-I_dark_til));
I_smooth=I_med-I_detail;
I_dark_cal(:,:,1)=miu*I_dark_til;I_dark_cal(:,:,2)=I_smooth;
I_dark=min(I_dark_cal,[],3);
t=1.0-w*I_dark;
J_inv=(double(I_inv)-double(A))./t+double(A);
J=255.0-J_inv;
J=uint8(J);imwrite(J,'output.png');
3.实验结果
输入图片从原论文中截取,可能会有部分失真,效果跟原论文中略有差异(可能因为只有单张图片,也可能有其他原因,欢迎指正)。
输入图片:
输出图片:
matlab自带直方图均衡效果: