function?[I_W??S]?=?func_DWT（I?level?Lo_D?Hi_D）;
%通過這個函數將I進行小波分解，并將分解后的一維向量轉換為矩陣形式
%?Matlab?implementation?of?SPIHT?（without?Arithmatic?coding?stage）
%
%?Wavelet?decomposition
%
%?input:????I?:?input?image
%???????????level?:?wavelet?decomposition?level
%???????????Lo_D?:?low-pass?decomposition?filter
%???????????Hi_D?:?high-pass?decomposition?filter
%
%?output:???I_W?:?decomposed?image?vector
%???????????S?:?corresponding?bookkeeping?matrix
%
%?please?refer?wavedec2?function?to?see?more
%
%?Jing?Tian
%?Contact?me?:?scuteejtian@hotmail.com
%?This?program?is?part?of?my?undergraduate?project?in?GuangZhou?P.?R.?China.
%?April?-?July?1999

[CS]?=?func_Mywavedec2（IlevelLo_DHi_D）;?

S（:3）?=?S（:1）.*S（:2）;????????%?dim?of?detail?coef?nmatrices?求低頻和每個尺度中高頻的元素個數
%st=S（13）+S（23）*3+S（33）*3;%%%%對前兩層加密
%C（1:st）=0;

L?=?length（S）;?%a求S的列數

I_W?=?zeros（S（L1）S（L2））;%設一個與原圖像大小相同的全零矩陣

%?approx?part
I_W（?1:S（11）??1:S（12）?）?=?reshape（C（1:S（13））S（11:2））;?%將LL層從C中還原為S（11）*S（22）的矩陣

for?k?=?2?:?L-1???%%%%%%%%%%%%%%%%%%%%%%將C向量中還原出HLHHLH?矩陣
????rows?=?[sum（S（1:k-11））+1:sum（S（1:k1））];
????columns?=?[sum（S（1:k-12））+1:sum（S（1:k2））];
????%?horizontal?part
????c_start?=?S（13）?+?3*sum（S（2:k-13））?+?1;
????c_stop?=?S（13）?+?3*sum（S（2:k-13））?+?S（k3）;
????I_W（?1:S（k1）??columns?）?=?reshape（?C（c_start:c_stop）??S（k1:2）?）;

????%?vertical?part
????c_start?=?S（13）?+?3*sum（S（2:k-13））?+?S（k3）?+?1;
????c_stop?=?S（13）?+?3*sum（S（2:k-13））?+?2*S（k3）;
????I_W（?rows??1:S（k2）?）?=?reshape（?C（c_start:c_stop）??S（k1:2）?）;

????%?diagonal?part
????c_start?=?S（13）?+?3*sum（S（2:k-13））?+?2*S（k3）?+?1;
????c_stop?=?S（13）?+?3*sum（S（2:k3））;
????I_W（?rows??columns?）?=?reshape（?C（c_start:c_stop）??S（k1:2）?）;

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1883??2015-05-31?13:33??func_DWT.m
?????文件????????1683??2015-05-31?13:33??func_Myappcoef2.m