Hello

Thanks for ur replies.....
I did that and now I get one more error. Basically I have wrritten
small code which are used in the final code. such as extbeamlts.m and
it is used here. Again this extbeamlts.m file has 2 different codes in
them named extractboundary.m, integral.m and weight4.m .

Please see this errror and let me know what to do.. is it because the
size of my image is too large ? How can I decrease t he size o an
image ? and how should I know what is size limitations ?

Thanks

code

im = imread('s.jpg');
im = rgb2gray(im);
im = im2double(im);
 g = im;
 im = imnoise(im,'speckle' , 0.6);
 im1 = ones(size(im));
 [m,n] = size(im);
 tic
 for i = 1:16:m-15
     for j = 1:16:n-15
         subim = im(i:i+15,j:j+15);
         [Energy wt sptx spty eptx epty]=extbeamlets(im)
         [selected xp0,yp0,xp1,yp1] = threshold(Energy,
sptx,spty,eptx,epty);
         result = drawbmlets(selected,xp0,yp0,xp1,yp1,subim);
         im1(i:i+15,j:j+15) = result;
     end
 end
 toc
 figure, imshow(im1);

 ??? Attempted to access im(95,85); index out of bounds because
size(im)=[94,150].

Error in ==> Integral at 25
        I=I+im(floor(y),floor(x));

Error in ==> weight4 at 35
             [I,im]=Integral(im,floor(x),floor(y),theta-pi/2,dir);

Error in ==> extbeamlets at 9
        wt(k)=weight4(im,xp(i),yp(i),xp(j),yp(j));

>>

Thanks!!!