clear;
close?all;
clc;

%%?基于暗通道的去霧算法
I?=?imread（‘../test?images/haze2.jpg‘）;

PR=I（::1）;????????????????????????????????????????%提取紅色分量?
PG=I（::2）;????????????????????????????????????????%提取綠色分量?
PB=I（::3）;????????????????????????????????????????%提取藍色分量?
gray?=?min（min（PRPG）PB）;
Jdark=ordfilt2（gray1ones（1515））;?????????????????%最小值濾波，得到暗通道

%從暗通道圖中按照亮度的大小取前0.1最亮的像素????????????????
[t?index]?=?sort（Jdark（:）‘?‘descend‘）;???????????????????
n?=?ceil（length（Jdark）*0.001）????????????????????????%前0.1的像素個數
threshold?=?max（t（1:n））;????????????????????????????%暗通道中最亮的閾值
index（1:n）;?????????????????????????????????????????%暗通道閾值點集的位置
%根據暗通道所取閾值和位置，找到原圖對應位置里的最亮通道值，構造全局大氣光A矩陣
R?=?PR（:）‘;
G?=?PG（:）‘;?
B?=?PB（:）‘;?
for?i?=?1:n?????????????????????????????????????????%結合閾值求原圖通道最大值
????a?=?min（threshold?max（max（R（index（i））?G（index（i）））?B（index（i））））;
end
A?=?max（a）
%求透射率t
w?=?0.95;
A?=?I*0?+?A;
t?=?1?-?w*double（gray）./double（A（::1））;???????????%存在負數
tmean?=?mean（mean（t））
%求目標圖像J
t?=?max（t0.1）;?????????????????????????????????????%除去負數，防止白場過渡

J（::1）?=?（double（PR）?-?double（A（::1）））./t?+?double（A（::1））;
J（::2）?=?（double（PG）?-?double（A（::2）））./t?+?double（A（::2））;
J（::3）?=?（double（PB）?-?double（A（::3）））./t?+?double（A（::3））;
J?=?uint8（J）;???????????????????????????????????????%一定要轉回與I相同的類型，不然顯示不出（用round無效）
%直方圖均衡化
J（::1）?=?histeq（J（::1））;
J（::2）?=?histeq（J（::2））;
J（::3）?=?histeq（J（::3））;
figure;?
subplot（221）;imshow（I）;title（‘原圖像‘）;?
subplot（222）;imshow（Jdark）;title（‘暗通道‘）;?
subplot（223）;imshow（t）;title（‘透射率圖‘）;?
subplot（224）;imshow（J）;title（‘目標圖‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2020-04-02?00:41??圖像去霧實驗代碼\
?????文件???????75059??2020-03-29?23:27??圖像去霧實驗代碼\haze1_w0.95_Amax_histeqNO_filt15.png
?????文件???????79505??2020-03-29?23:26??圖像去霧實驗代碼\haze1_w0.95_Amin_histeqNO_filt15.png
?????文件???????85204??2020-03-29?23:06??圖像去霧實驗代碼\haze1_w0.95_Amin_histeq_filt15.png
?????文件???????86647??2020-03-29?23:43??圖像去霧實驗代碼\haze2_w0.95_Amin_histeqNO_filt15.png
?????文件???????95030??2020-03-29?23:07??圖像去霧實驗代碼\haze2_w0.95_Amin_histeq_filt15.png
?????文件???????70235??2020-03-29?23:07??圖像去霧實驗代碼\haze3_w0.95_Amin_histeq_filt15.png
?????文件???????90023??2020-03-29?23:37??圖像去霧實驗代碼\haze4_w0.6_Amin_histeqNO_filt15.png
?????文件???????92984??2020-03-29?23:36??圖像去霧實驗代碼\haze4_w0.8_Amin_histeqNO_filt15.png
?????文件???????87098??2020-03-29?23:36??圖像去霧實驗代碼\haze4_w0.95_Amin_histeqNO_filt15.png
?????文件???????92268??2020-03-29?23:08??圖像去霧實驗代碼\haze4_w0.95_Amin_histeq_filt15.png
?????文件???????88269??2020-03-29?23:08??圖像去霧實驗代碼\haze5_w0.95_Amin_histeq_filt15.png
?????文件????????1862??2020-03-29?23:59??圖像去霧實驗代碼\HazeRemoval.m
?????文件?????4389504??2020-03-20?12:13??圖像去霧實驗代碼\Single?Image?Haze?Removal?Using?Dark?Channel?Prior_cvpr_09.pdf
?????目錄???????????0??2020-04-02?00:40??圖像去霧實驗代碼\test?images\
?????文件???????17105??2020-03-17?22:52??圖像去霧實驗代碼\test?images\haze1.jpg
?????文件???????17666??2020-03-17?22:52??圖像去霧實驗代碼\test?images\haze2.jpg
?????文件???????14261??2020-03-17?22:52??圖像去霧實驗代碼\test?images\haze3.jpg
?????文件???????20033??2020-03-17?22:52??圖像去霧實驗代碼\test?images\haze4.jpg
?????文件???????21264??2020-03-17?22:52??圖像去霧實驗代碼\test?images\haze5.jpg
?????文件??????358533??2020-03-29?20:54??圖像去霧實驗代碼\test?images\haze6.jpg