lpdenhance_main.m

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

function lpdenhance_main
clc
clear all

in=imread('212.bmp');
im=rgb2gray(in);
figure,imshow(in);title('原始图像');axis on;
im = double(im) / 256;
n=5;
pflit='9-7';
lpy=lpd(im,pflit,n);
for i=4:6
lpy{i}=wdec(lpy{i});
end
imrec=lpr(lpy,pflit);
figure,imshow(imrec);title('增强后的图像');axis on;

function coef=wdec(x)
[m,n]=size(x);
[ca,ch,cv,cd]=dwt2(x,'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
coef=idwt2(ca,ch,cv,cd,'db1');