将二进小波变换用于图像增强的MATLAB代码

将二进小波变换用于图像增强的MATLAB代码。

function [ALLALHAHLAHH] = MyDWT2（Ahg）
%   MyDWT2 Single-level discrete 2-D wavelet transform.
%   DWT2 performs a single-level 2-D wavelet decomposition
%   with respect to wavelet filters （h g）.
%
%   [CACHCVCD] = DWT2（XLo_DHi_D） computes the 2-D wavelet
%   decomposition as above given these filters as input:
%   h is the decomposition low-pass filter.
%   g is the decomposition high-pass filter.
%   h and g must be the same length.
%
%   size（ALL） = size（ALH） = size（AHL） = size（AHH） = 1/2 Size（A）


lh = length（h）;lg = length（g）;
sA = size（A）;s1 = sA（11）; s2 = sA（12）;
sizeEXTh = （lh-1）/2; sizeEXTg=（lg-1）/2;
sizeKEPT = 2*ceil（sA/2）;

%对A作周期延拓：边界两边分别增加sizeEXTh行和sizeEXTg列对原图象分别用h，g作行滤波，得到AL和AH

y = wextend（‘2D‘‘per‘A[0sizeEXTh]）;
AL = wconv（‘row‘yh）; 
AL = wkeep（ALsA）;%figureimshow（AL/256）;
y = wextend（‘2D‘‘per‘A[0sizeEXTg]）;
AH = wconv（‘row‘yg）; 
AH = wkeep（AHsA）; 

%对行变换的结果AL和AH再分别用hg作列变换，得到
y = wextend（‘2D‘‘per‘AL[sizeEXTh0]）;
ALL = wconv（‘col‘yh）;ALL = wkeep（ALLsizeKEPT）;
y = wextend（‘2D‘‘per‘AL[sizeEXTg0]）;
ALH = wconv（‘col‘yg）;ALH = wkeep（ALHsizeKEPT）;% x方向的高频细节
y = wextend（‘2D‘‘per‘AH[sizeEXTh0]）;
AHL = wconv（‘col‘yh）;AHL = wkeep（AHLsizeKEPT）;% y方向的高频细节
y = wextend（‘2D‘‘per‘AH[sizeEXTg0]）;
AHH =wconv（‘col‘yg）;AHH = wkeep（AHHsizeKEPT）;

ALL = dyaddown（ALL‘c‘1）;%取偶数列
ALL = dyaddown（ALL‘r‘1）;%取偶数行
ALH = dyaddown（ALH‘c‘1）;
ALH = dyaddown（ALH‘r‘1）;%
AHL = dyaddown（AHL‘c‘1）;%
AHL = dyaddown（AHL‘r‘1）;
AHH = dyaddown（AHH‘c‘1）;%
AHH = dyaddown（AHH‘r‘1）;%

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    ......R      1672  2003-05-29 22:13  二进小波变换的图像增强\MyDWT2.m

    ......R      2729  2003-06-06 09:32  二进小波变换的图像增强\NonlinearEnhancement.m

    ......R      1869  2003-06-06 09:39  二进小波变换的图像增强\WaveletEnhance_IE3_Global.m

    ...D..R         0  2004-04-23 20:13  二进小波变换的图像增强

----------- ---------  ---------- -----  ----

                 6270                    4