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