%首先，顯示原始圖像
I1=imread('rice.tif');
figure,imshow(I1)
%原始圖像如圖9-2所示
%然后，選取閾值為0.2，對原始圖像進行四叉樹分解，并以灰度圖的形式顯示分解所得的稀疏矩陣。
S = qtdecomp(I1,0.2);
S2=full(S);
figure,imshow(S2)
%稀疏矩陣如圖9-3所示
%最后，通過查看S2來了解四叉樹分解的具體結果
[vals,r,c] = qtgetblk(I1,S2,2);
size(vals)
[vals2,r,c] = qtgetblk(I1,S2,4);
size(vals2)
[vals3,r,c] = qtgetblk(I1,S2,8);
size(vals3)
[vals4,r,c] = qtgetblk(I1,S2,16);
size(vals4)
[vals5,r,c] = qtgetblk(I1,S2,32);
size(vals5)
[vals6,r,c] = qtgetblk(I1,S2,1);
size(vals6)