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

基于直方图均衡化法和小波分解的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); end end %利用双正交小波函数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); end end y=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); figure,imshow(E33); 经过前文介绍的算法，处理所给的图像，最终以一种算法过程，获得了弱小目标，便于观察弱小目标—镁溶液第一气泡的特征。