基于修正的retinex霧天圖像增強(qiáng)算法-資料下載頁
2025-06-19 12:45本頁面
【正文】 /2)^2)/(2*sigma*sigma))。 endendF = F./(sum(F(:)))。 Ffft=fft2(double(F))。 DB0=Bfft2.*Ffft。DB0=fftshift(DB0)。DB=ifft2(DB0)。 DBdouble=double(DB)。 u=mean2(DBdouble)。v=std2(DBdouble)。 u=mean2(DBdouble)。v=std2(DBdouble)。MIN = u3*v。MAX = u+3*v。for i=1:m for j=1:n if DBdouble(i,j)=MIN Rb(i,j)=0。 end if DBdouble(i,j)MINamp。amp。DBdouble(i,j)MAX Rb(i,j)=(DBdouble(i,j)MIN)./6*v。 end if DBdouble(i,j)=MAX Rb(i,j)=1。 end endendRl=(B0+)./(Rb+)。mi=min(min(Rl))。ma=max(max(Rl))。t=(mi+ma)./2。for i=1:m for j=1:n if Rl(i,j)miamp。amp。Rl(i,j)t Rl1(i,j)=mi+(tmi)*[(Rl(i,j)mi)./(tmi)]^2。 end if Rl(i,j)tamp。amp。Rl(i,j)ma Rl1(i,j)=ma(mat)*[(maRl(i,j))./(mat)]^2。 end endendIb_double=Rb.*Rl1。Final=cat(3,uint8(Ir_double),uint8(Ig_double),uint8(Ib_double))。 %將三通道圖像合并。 figure,imshow(Final)。 %顯示處理后的圖像。%%%MSR算法%%%%f=imread(39。39。)。fr=f(:, :, 1)。 fg=f(:, :, 2)。 fb=f(:, :, 3)。 mr=mat2gray(im2double(fr))。 mg=mat2gray(im2double(fg))。 mb=mat2gray(im2double(fb))。%%數(shù)據(jù)類型歸一化。alf1=1458。 %定義標(biāo)準(zhǔn)差alf=a^2/2 a=54n=161。%定義模板大小。n1=floor((n+1)/2)。%確定中心。 for i=1:n for j=1:n b(i,j) =exp(((in1)^2+(jn1)^2)/(4*alf1))/(pi*alf1)。 %高斯函數(shù)。 end end nr1 = imfilter(mr,b,39。conv39。, 39。replicate39。)。ng1 = imfilter(mg,b,39。conv39。, 39。replicate39。)。nb1 = imfilter(mb,b,39。conv39。, 39。replicate39。)。%卷積濾波。ur1=log(nr1)。 ug1=log(ng1)。 ub1=log(nb1)。tr1=log(mr)。tg1=log(mg)。tb1=log(mb)。yr1=(tr1ur1)/3。yg1=(tg1ug1)/3。yb1=(tb1ub1)/3。alf2=。 %定義標(biāo)準(zhǔn)差alf=a^2/2 a=x=31。%定義模板大小。x1=floor((n+1)/2)。%確定中心。for i=1:n for j=1:n a(i,j) =exp(((in1)^2+(jn1)^2)/(4*alf2))/(6*pi*alf2)。 %高斯函數(shù)。 end end nr2 = imfilter(mr,a,39。conv39。, 39。replicate39。)。ng2 = imfilter(mg,a,39。conv39。, 39。replicate39。)。nb2 = imfilter(mb,a,39。conv39。, 39。replicate39。)。%卷積濾波。ur2=log(nr2)。 ug2=log(ng2)。 ub2=log(nb2)。tr2=log(mr)。tg2=log(mg)。tb2=log(mb)。yr2=(tr2ur2)/3。yg2=(tg2ug2)/3。yb2=(tb2ub2)/3。alf3=。 %定義標(biāo)準(zhǔn)差alf=a^2/2 a=167l=501。%定義模板大小。l1=floor((n+1)/2)。%確定中心。for i=1:n for j=1:n e(i,j) =exp(((in1)^2+(jn1)^2)/(4*alf3))/(4*pi*alf3)。 %高斯函數(shù)。 end end nr3 = imfilter(mr,e,39。conv39。, 39。replicate39。)。ng3 = imfilter(mg,e,39。conv39。, 39。replicate39。)。nb3 = imfilter(mb,e,39。conv39。, 39。replicate39。)。%卷積濾波。ur3=log(nr3)。 ug3=log(ng3)。 ub3=log(nb3)。tr3=log(mr)。tg3=log(mg)。tb3=log(mb)。yr3=(tr3ur3)/3。yg3=(tg3ug3)/3。yb3=(tb3ub3)/3。dr=yr1+yr2+yr3。dg=yg1+yg2+yg3。db=yb1+yb2+yb3。 cr=im2uint8(dr)。 cg=im2uint8(dg)。 cb=im2uint8(db)。z=cat(3, cr, cg, cb)。 figure, imshow(z)%%%顯示處理后的圖像寧可累死在路上，也不能閑死在家里！寧可去碰壁，也不能面壁。是狼就要練好牙，是羊就要練好腿。什么是奮斗？奮斗就是每天很難，可一年一年卻越來越容易。不奮斗就是每天都很容易，可一年一年越來越難。能干的人，不在情緒上計(jì)較，只在做事上認(rèn)真；無能的人！不在做事上認(rèn)真，只在情緒上計(jì)較。拼一個(gè)春夏秋冬！贏一個(gè)無悔人生！早安！—————獻(xiàn)給所有努力的人.學(xué)習(xí)
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1