基于直方图均衡化法和小波分解的matlab图像处理的一种镁溶液第一气泡的检测方案（含源代码）.doc

基于直方图均衡化法和小波分解的matlab图像处理的一种镁溶液第一气泡的检测方案（含源代码）面对复杂背景的弱小目标检测，我们首先采用直方图均衡化法，提高图像的对比度，再通过利用构造的小波基进行小波分解重构处理，下面仔细介绍下程序的算法过程：%目标检测完整程序 %完整程序%读入图像image_1=imread(00228.bmp);%直方图均衡化图像image_2=histeq(image_1);%显示原始灰度图像和均衡化后的图像figure,imshow(uint8(image_1);figure,imshow(uint8(image_2);image_3=imread(00230.bmp);image_4=histeq(image_3);figure,imshow(uint8(image_3);figure,imshow(uint8(image_4);%将处理后的两图像作差分运算I=imsubtract(uint8(image_4),uint8(image_2);figure,imshow(uint8(I);%利用双正交小波函数bint53.9进行小波分解c,s=wavedec2(I,4,bint53.9);sizec=size(c);%在频域里对图像进行增强for i=1:sizec(2) if(c(i)350) c(i)=2*c(i); else c(i)=0.5*c(i); endend%利用双正交小波函数bint53.9进行小波重构C=waverec2(c,s,bint53.9);%显示重构图像figure,imshow(uint8(C); f=double(C);%二维离散傅里叶变换k=fft2(f);%将零点平移到频谱中心g=fftshift(k);%高频增强N1,N2=size(g);n=3;d0=20;%朝最近的整数取整u0=round(N1/2);v0=round(N2/2);for i=1:N1 for j=1:N2 d=sqrt(i-u0)2+(j-v0)2); h=1/(1+0.414*(d/d0)(2*n); y(i,j)=h*g(i,j); endendy=ifftshift(y);E1=ifft2(y);E2=uint8(real(E1);%显示图像figure,imshow(uint8(E2); %增强后图像和原小波重构后的图像作差I=imsubtract(uint8(C),uint8(E2);figure,imshow(I); %二值化处理，将图像黑白明显D=im2bw(I,0.06);figure,imshow(D); %数学形态学开运算se=strel(disk,1);E=imopen(D,se);figure,imshow(E); %发现所要求的目标相比另一噪点不明显，重新二值化及数学形态学处理，使得噪点明显，最好只有噪点，与E作差获得目标点即可%修改二值化系数D1=im2bw(I,0.07);figure,imshow(D1); se=strel(disk,1);E11=imopen(D1,se);figure,imshow(E11); %将E11和E1作差I33=imsubtract(E,E11);figure,imshow(I33); se=strel(disk,1);E33=imopen(I33,se);f