lpwenhance_main.m

来自「包内提供了Contourlet和Wavelet变换进行图像增强、去噪的各种算法和」· M 代码 · 共 39 行

function lpwenhance_main
clc
clear all

in=imread('212.bmp');
im=rgb2gray(in);
figure,imshow(in);title('原始图像');axis on;
im = double(im) / 256;
pfilt='9-7';                % choose LP decomposition filter
[h,g]=pfilters(pfilt);
[co,de]=lpdec(im,h,g);
%figure,imshow(co);title('粗糙图像'),axis on;
%figure,imshow(de);title('细节图像'),axis on;

[m,n]=size(de);
[ca,ch,cv,cd]=dwt2(de,'db1');
g=sqrt(ch.^2+cv.^2);
mj=200;
c=3;
p=0.5;
for i=1:m/2
    for j=1:n/2
        if abs(g(i,j)<c)
            y=(mj/c)^p;
        elseif abs(g(i,j))>=mj
            y=1;
        else y=(mj/abs(g(i,j)))^p;
        end
        ch(i,j)=y*ch(i,j);
        cv(i,j)=y*cv(i,j);
    end
end
de=idwt2(ca,ch,cv,cd,'db1');
[h,g]=pfilters(pfilt);
imrec=lprec(co,de,h,g);
g=1.1;
newim=adjgamma(imrec,g);
figure,imshow(newim);title('增强后的图像');axis on;