wavelet decomposition

I have an image 480*640 and i want to use wavelet decomposition in order to produce 4 subimages within a small neighborhood window of size 9*9. how should i do that? do i have to use blocproc?

the code is:

I = imread('4.bmp');

[F1,F2] = wfilters('haar','d') [C,S] = wavedec2(I,2,F1,F2);

fs1 = appcoef2(C,S,'haar',1); fh1 = detcoef2('h',C,S,1); fv1 = detcoef2('v',C,S,1); fd1 = detcoef2('d',C,S,1);

fs2 = appcoef2(C,S,'haar',2); fh2 = detcoef2('h',C,S,2); fv2 = detcoef2('v',C,S,2); fd2 = detcoef2('d',C,S,2);

i want to do that in blocks 9*9

thanks in advance