【图像增强】基于频域+密度分割+灰度级-彩色变换实现伪彩色图像增强含Matlab源码 .zip

Image=imread('qiguan.bmp'); [height,width]=size(Image); NewImage=zeros(height,width,3); G=fft2(Image); n=2; %阶数n为2 RedD0=150; GreenD0=200; %低通、高通滤波器参数设置 Bluecenter=150; Bluewidth=100; Blueu0=10; Bluev0=10; %带通滤波器参数设置 for u=1:height for v=1:width D(u,v)=sqrt(u^2+v^2); RedH(u,v)=1/(1+(sqrt(2)-1)*(D(u,v)/RedD0)^(2*n)); %低通滤波器转移函数 GreenH(u,v)=1/(1+(sqrt(2)-1)*(GreenD0/D(u,v))^(2*n)); %高通滤波器转移函数 BlueD(u,v)=sqrt((u-Blueu0)^2+(v-Bluev0)^2); BlueH(u,v)=1-1/(1+BlueD(u,v)*Bluewidth/((BlueD(u,v))^2-(Bluecenter)^2)^(2*n)); %带通滤波器转移函数 end end Red=RedH.*G; RedC=ifft2(Red); Green=GreenH.*G; GreenC=ifft2(Green); Blue=BlueH.*G; BlueC=ifft2(Blue); NewImage(:,:,1)=uint8(real(RedC)); NewImage(:,:,2)=uint8(real(GreenC)); NewImage(:,:,3)=uint8(real(BlueC)); figure;imshow(NewImage);title('频域伪彩色增强'); imwrite((NewImage),'pinwei.bmp');