暗通道去雾，实现图像增强处理

%这个程序用于图片去雾 clear;clc;%close all; tic; a=imread('q.jpg'); %暗通道去雾技术 m=20;n=20; % for i=1:1 % a0=padarray(a(:,:,i),[m,n],'replicate'); % g(:,:,i)=colfilt(a0,[m,n],'sliding',@gmean); % end; % imshow(uint8(g)); toc [g,ga]=antongdao(a,[m n]);toc [M,N]=size(a(:,:,1)); % g=zeros(M,N,3)+0.5; %imshow(g1); % ga=double(max(g(:))); % ga=5*g; ga=double(ga); ep=180; %为缓解图片中天空去雾效果不佳的现象而设置的最高阈值 % ga(ga>=ep)=ep; t0=0.1; %为缓解图片中因透射率值过低 %而导致图片向白场过度的现象而设置的低阈值 gt=1-double(g)./ga; %适当乘以0.95以保留一定程度的雾 % gt((gt<t0))=t0; w=fspecial('average',[10*m,10*n]);%加窗 %均值滤波器 % [M,N]=size(a(:,:,1)); % [I,J]=meshgrid(1:N,1:M); bb=zeros(M,N,3); %预分配内存 %%向量化------------------------------------ i=1:3; [gt(:,:,i),ga(:,:,i)]=guidefilter2(a(:,:,i),gt(:,:,i),ga(:,:,i),50,0.001); gt(gt<=t0)=t0; % ga(:,:,i)=guidefilter(a(:,:,i),ga(:,:,i),50,0.001); % ga(ga>=ep)=ep; %导向滤波获取透射率图 % bb(:,:,i)=(a00-ga(:,:,1).*(1-gt(:,:,1)))./gt(:,:,1); bb(:,:,i)=(double(a(:,:,i))-ga(:,:,i))./gt(:,:,i)+ga(:,:,i); % bb(:,:,i)=bb(:,:,i)+0.1*(255-bb(:,:,i));%适当增加亮度 % bb(:,:,i)=mat2gray(bb(:,:,i)); %%----------------------------------------- % figure,imshow([gt,mat2gray(ga)]); toc;%与tic来记录程序运行时间，toc记录程序完成时间，tic用来保存当前时间 figure();imshow([a bb]); figure();imshow(uint8(bb));xlabel('暗通道处理后的图像'); % figure();subplot(121);imshow(a);title('原图'); % subplot(122);imshow(bb);title('处理后的图像');