一种图像增强技术的实现

function [IE,Mean]=DerivationEnhancer(Img) % 函数功能：获得增强子。 % 输入参数：输入图像 Img % 输出参数：增强子图像 IE % 时间：2017-8-31 %clear %Img=rgb2gray(imread('分割图像.JPG')); %Img=rgb2gray(imread('lena.jpg')); [H,L]=size(Img); %计算输入图像Img的像素的平均值 sum=0; for i=1:H for j=1:L sum=sum+double(Img(i,j)); end end Mean=sum/(H*L); Median=median(reshape(Img,1,H*L)); %将输入图像划分成亮暗两个区域,命名为BT，DT BT=zeros(H,L); for i=1:H for j=1:L if Img(i,j)>=Mean BT(i,j)=Img(i,j); else DT(i,j)=Img(i,j); end end end %中值滤波的方法得到模糊图像BI BI=medfilt2(Img,[3,3]); %小波滤波的方法得到模糊图像BI %BI=Blurred(Img); %亮增强子 for i=1:H for j=1:L if BI(i,j)>=Mean BE(i,j)=double((BI(i,j)-Mean))/double((256-Mean)); else BE(i,j)=0; end end end %暗增强子 for i=1:H for j=1:L if BI(i,j)<=Mean DE(i,j)=(Mean-double(BI(i,j)))/Mean; else DE(i,j)=0; end end end %增强子IE %计算暗区域非零元素的平均值 dsum=0; dcount=0; %暗区域的像素个数 for i=1:H for j=1:L dsum=dsum+double(DT(i,j)); if (DT(i,j)>0) dcount=dcount+1; end end end dMean=double(dsum/dcount); %暗区域的平均值 DTZ=nonzeros(DT)'; %取出暗区域的非零元素 dMedian=median(DTZ); %计算中值 bcount=H*L-dcount; %亮区域的像素个数 DR=(double(dMean)/double(dMedian)*double((128-dMean))/128+double(Mean)/double(Median)*double(dcount)/double(bcount))/4;%暗比率 bsum=0; for i=1:H for j=1:L bsum=bsum+double(BT(i,j)); end end bMean=double(bsum/bcount); %亮区域非零元素的均值 BR=1-(Mean/bMean)^0.5; %增强图像IE for i=1:H for j=1:L IE(i,j)=Img(i,j)*(1+DR*DE(i,j)-BR*BE(i,j)); end end % subplot(2,1,1) % imshow(Img); % subplot(2,1,2) % imshow(IE); imwrite(IE,'增强子.JPG');