clc
close?all
clear?
%%實現論文??暗通道先驗的大壩水下裂縫圖像增強算法?中勻光算法
srcImage=imread（‘5.png‘）;
imshow（srcImage）
title（‘彩色圖像‘）;
src_gray=rgb2gray（srcImage）;%彩色圖像轉灰度圖像
figure;imshow（src_gray）;
title（‘灰度圖像‘）
%圖像勻光處理
num?=?91;%設置勻光處理的像素帶個數，必須為奇數。
[mnk]?=?size（srcImage）;
pixelNum?=?m*n;%圖像像素個數
m_m?=?floor（m/4）;%掩模的尺寸
n_n?=?floor（n/4）;
X_mask?=?（1/（m_m*n_n））*ones（m_mn_n）;%矩形濾波掩模

line_filting=imfilter（src_grayX_mask‘replicate‘）;%通過對原圖灰度圖進行模糊處理，擬合出光照強弱分布圖
figure;imshow（line_filting）;title（‘光照分布圖‘）;
line_filting=im2double（line_filting）;

%??光照分布圖的最大值與最小值
mask_max=max（max（line_filting））;
mask_min=min（min（line_filting））;
pix_diff=mask_max-mask_min+1.0/255;?%像素值差

gray_distance?=?double（pix_diff）/num;%每個區域的灰度間隔
normal?=?ceil（num/2）;%正常區域
mask=zeros（mnnum）;
program_mean1=zeros（1num）;%每個區域的均值
program_std1=zeros（1num）;%每個區域的均方差

for?k=1:num
????src_gray1=im2double（src_gray）;
????mask（::k）=（mask_min+gray_distance*（k-1）-0.0001）????src_gray1（~mask（::k））=0;
????if（k==normal）
???????[normal_meankk]=mean1（src_gray1）;%求正常區域的均值跟標準差
???????program_mean1（k）=normal_mean;
???????normal_std=std1（src_gray1normal_meankk）;
???????program_std1（k）=normal_std;
????else
????????[program_meane]=mean1（src_gray1）;%問題區域均值
????????program_mean1（k）=program_mean;
????????program_std=std1（src_gray1program_meane）;%問題區域標準差
????????program_std1（k）=program_std;
????end
end
%
out_gray=zeros（size（src_gray））;%定義勻光之后的數組
for?i=1:num
????src_gray2=im2double（src_gray）;
?????if（i==normal）
?????????src_gray2（mask（::i）==0）=0;
?????????out_gray=src_gray2+out_gray;
?????else
?????????src_gray2（mask（::i）==0）=0;
?????????out_gray=（Uniformlight（?src_gray2out_graynormal_meannormal_stdprogram_mean1（i）program_std1（i）））;
?????end
end
figure;imshow（out_gray[]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2061??2018-01-27?10:47??勻光算法\dodging.m

?????文件????????519??2018-01-26?16:43??勻光算法\mean1.m

?????文件????????591??2017-10-14?14:51??勻光算法\std1.m

?????文件????????470??2017-10-14?19:22??勻光算法\Uniformlight.m

?????目錄??????????0??2018-01-27?10:46??勻光算法

-----------?---------??----------?-----??----

?????????????????3641????????????????????5