example9_4.m

来自「《matlab6.5辅助图像处理》源代码」· M 代码 · 共 23 行

%首先，显示原始图像
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)