%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