彩色图像增强的同态滤波的MATLAB实现

%function homofilter（imdrcn） 
%HOMOFILTER Summary of this function goes here
%   Detailed explanation goes here
%%%%%%%%%%Butterworth high pass filter %%%%%%%%%%%%%% 
clc 
close all 
clear all 
histgram=zeros（1256）;             %生成直方图数组并置0
cdf=zeros（1256）;
d=1; 
n=2;
%img=（imread（‘hofi.bmp‘））; 
img=（imread（‘E:\MATLAB project\photos\water\cai1.png‘））; 
 
[r c]=size（img（::1））; 
A=zeros（rc）; 
H=zeros（rc）;
for i=1:r 
    for j=1:c 
        R=（（（i-r/2）.^2+（j-c/2）.^2））.^（.5）; 
        H（ij）=1/（1+（d/R）^（2*n））; 
    end 
end 
%%%%%%%%%%%%%Using it for my application as homomorphic filtering is 
%%%%%%%%%%%%%application specific taking the value of alphaL and alphaH 
%%%%%%%%%%%%%values accordingly. 
alphaL=0.1; 
aplhaH=1.01; 
H=（（aplhaH-alphaL）.*H）+alphaL; 
H=1-H; 
im_e=img;

for k=1:3
%%%%%log of image 
im_l=log2（1+double（img（::k）））; 
%%%%%DFT of logged image 
im_f=fft2（im_l）; 
%%%%%Filter Applying DFT image 
im_nf=H.*im_f;   
%%%%Inverse DFT of filtered image 
im_n=abs（ifft2（im