%-----------------------------5/3整數小波3級分解函數---------------------------------
??
function?[S11d11]=decompose53（imMatdimn）
%?imMat?輸入圖像的數據變量
%??dim??輸入圖像的維數
%???n???表示第n次分解
%??S11??表示分解的低頻部分
%??d11??表示分解的高頻部分

%?imMat=imread（‘cameraman.tif‘）;
%?dim=size（imMat1）;
%?n=1;
imMat=imMat（1:dim/n1:dim/n）;
for?m=1:dim/n
????for?i=1:2:dim/n-1
????????O1（m（（i+1）/2））=imMat（mi）;%得到第一級奇序列
????end
end
for?m=1:dim/n
????for?i=2:2:dim/n
????????E1（mi/2）=imMat（mi）;%得到第一級偶序列
????end
end
E1=[E1zeros（dim/n2）];%再加兩列相等行數的零元素（對于行變換）
for?m=1:dim/n
????for?i=1:dim/（n*2）
????????d1（mi）=O1（mi）-E1（mi）/2-E1（mi+1）/2;
????end
end
d_1=[d1zeros（dim/n2）];%再加兩列相等行數的零元素（行）
for?m=1:dim/n
????for?i=1:dim/（2*n）
????????S1（mi）=E1（mi）+（d_1（mi）+d_1（mi+2）+2）/4;
????end
end
L1=[S1d1];
for?m=1:dim/n
????for?i=1:2:dim/n-1
????????O11（（（i+1）/2）m）=L1（im）;%得到列變換的第一級奇序列
????end
end
for?m=1:dim/n
????for?i=2:2:dim/n
????????E11（i/2m）=L1（im）;%得到列變換的第一級偶序列
????end
end
E11=[E11;zeros（2dim/n）];%再加兩行相等列數的零元素（列）
for?m=1:dim/n
????for?i=1:dim/（2*n）
?????????d11（im）=O11（im）-E11（im）/2-E11（i+1m）/2;%各行差值數據；不對細節數據再做變換，只對S1再預測
????end
end
d_11=[d11;zeros（2dim/n）];%再加兩行相等列數的零元素（列）
for?m=1:dim/n
????for?i=1:dim/（2*n）
??????????S11（im）=E11（im）+（d_11（im）+d_11（i+2m）+2）/4;?
???end
end?

%?imshow（S11[]）
%?imwrite（uint8（S11）‘cameraman_1.bmp‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1477??2017-05-28?20:14??recompose53.m

?????文件???????1564??2017-11-23?19:59??decompose53.m

-----------?---------??----------?-----??----

?????????????????3041????????????????????2