[f1,f2] = freqspace(size(J), 'meshgrid' ); D = sqrt(f1.^2 + f2.^2); H_b=1./((1+0.1./D).^2); %Butterworth high-pass filter H_em=0.5+0.75*H_b; %High frequency emphasis filter H_em=ifftshift(H_em); % Perform filtering in frequency domain
I_f=fft2(J); I_f1=I_f.*(H_em); I2=uint8(ifft2(I_f1)); imshow(I2); histeq(I2);
Best Answer