图像分裂合并的matlab实现.docxVIP

图像分裂合并的matlab实现 包含四段程序 CODE1: clear;I=imread(xingshi32.bmp);if(isgray(I)==0)disp(请输入灰度图像,本程序用来处理128 *128的灰度图像!);else?? if (size(I)~=[128,128])?? disp(图像的大小不合程序要求!);?? elseH.color=[1 1 1];?? %设置白的画布figure(H);imshow(I);title(原图像);??zeroImage=repmat(uint8(0),[128 128]);figure(H); ?? ??? %为分裂合并后显示的图设置画布meansImageHandle=imshow(zeroImage);title(块均值图像);%%%%%设置分裂后图像的大小由于本图采用了128像素的图blockSize=[128 64 32 16 8 4 2];%%设置一个S稀疏矩阵用于四叉树分解后存诸数据S=uint8(128);S(128,128)=0;??threshold=input(请输入分裂的阈值(0--1):);%阈值threshold=round(255*threshold);M=128;dim=128;tic%%%%%%%%%%%%%%%%% 分裂主程序%%%%%%%%%%%while (dim1)[M,N] = size(I);Sind = find(S == dim);numBlocks = length(Sind);if (numBlocks == 0)%已完成break;end?? rows = (0:dim-1);?? cols = 0:M:(dim-1)*M;?? rows = rows(:,ones(1,dim));?? cols = cols(ones(dim,1),:);?? ind = rows + cols;?? ind = ind(:);tmp = repmat(Sind, length(ind), 1);ind = ind(:, ones(1,numBlocks));ind = ind + tmp;blockValues= I(ind);blockValues = reshape(blockValues, [dim dim numBlocks]);if(isempty(Sind))?? ??? %已完成?? ??? break;end[i,j]=find(S);set(meansImageHandle,CData,ComputeMeans(I,S));maxValues=max(max(blockValues,[],1),[],2);minValues=min(min(blockValues,[],1),[],2);doSplit=(double(maxValues)-double(minValues))threshold;dim=dim/2;Sind=Sind(doSplit);Sind=[Sind;Sind+dim;(Sind+M*dim);(Sind+(M+1)*dim)];S(Sind)=dim;end[i,j]=find(S);?? % 用来寻找四叉机分解结果中大小为S的块的位置set(meansImageHandle,CData,ComputeMeans(I,S)); % 显示分解结果块均值图像Numberofbloks=length(i);?? %计算块数%sizev=size(v);endendtoc CODE2: function means = ComputeMeans(I, S)% 用来计算给定图像和稀疏矩阵的块均值%?? I: 为给定的图像%?? S: 为稀疏矩阵?? means = I;?? for dim = [128 64 32 16 8 4 2 1];values = getblk(I, S, dim);??if (~isempty(values))?? %%%%%以下的句子是将小块的平均值来代替原图像中相应的块处的像素%%%%?? if (min(min(values))=60)???? ?? ??? means = setblk(means, S, dim, 0);?? %用于合并时的阈值?? else??? ?? %means = setblk(means, S, dim, sum(sum(values,1),2) ./ dim^2+std2(values));?? %means = setblk(means, S, dim, sum(s