MATLAB: Problem adding components of a vector using a while loop

summing vectorswhile loop

Can anyone tell me what's wrong with this code?

x = [1 23 43 72 87 56 98 33];

b = 0;

k = 1;

while k<=length(x)

    b = b+k(x)
    k = k+1

end

I get index exceeds matrix dimensions. Thank you very much for your help

Best Answer

You should be using

b = b + x(k)

rather than

b = b + k(x)

x is your vector, k is the index into it.

Related Solutions

MATLAB: Write program by For loop, how can I stop

A=[40 47 50 60 70 75 80 80 90 55];
B(1)=70;
C(1)=60;
for k=1:9;
    C(k+1)= B(k)+C(k)-A(k);
    if C(k+1)<0,break,end  % Stops the loop.
    B(k+1)= A(k);
    if C(k+1)>=40
        B(k+1)=0
        k=k+1
    end   
end

MATLAB: Connect bwlabeled components

Using bwconncomp on a three dimensional image will group them in 3d. You'll have to call it on each individual slice to see how many objects per slice there are:

sliceObjs = zeros(size(stack,3),1);
for ii = 1:size(stack,3)
 CC = bwconncomp(stack(:,:,ii));
 sliceObjs(ii) = CC.NumObjects;
end

EDIT This is a script, with the engine completely generalized and ready to go. I saved you above image; cropped it and translated it for visualization error checking purposes. The translation was done using my function FEX:imtranslate

%%Data and stuff
I = imread('ans627.png'); %your image
I = I(35:276,86:381,1)==255; %keep relevant part; convert to binary
I(:,:,2) = imtranslate(I,[10 -14]); %artificially translate for testing
I(5:11,270:275,2) = true;  %add an eleventh object with no friends
szI = size(I); 
Ds = cell(szI(3)-1); %preallocate place to store stuff
CCold = bwconncomp(I(:,:,1));  %cc of first slice
RPold = regionprops(CCold,'centroid');
centsA = vertcat(RPold(:).Centroid); %extract centroids
for ii = 1:(szI(3)-1)
  CCnew = bwconncomp(I(:,:,ii+1));  %cc of next slice
  RPnew = regionprops(CCnew,'centroid');
  centsB = vertcat(RPnew(:).Centroid);
    dim = 2; %which dimension will control in the bsxfun expression?
    if CCold.NumObjects>CCnew.NumObjects
        dim = 1;
    end    
    %Engine:
    xyDiff = bsxfun(@(x,y)abs(x-y),reshape(centsA,[],1,2),reshape(centsB,1,[],2)); %Centroids distances in each dimension
    [~,idx] = min(hypot(xyDiff(:,:,1),xyDiff(:,:,2)),[],dim);
    if dim==2;
        displacements = centsB(idx,:)-centsA;
        H = quiver(centsA(:,1),centsA(:,2),displacements(:,1),displacements(:,2));
    else
        displacements = centsB-centsA(idx,:);
        H = quiver(centsB(:,1),centsB(:,2),displacements(:,1),displacements(:,2));
    end
    RPold = RPnew; %get ready to move on
    CCold = CCnew; 
    centsA = centsB;
    Ds{ii} = displacements;  %can store other stuff too.
end
%SCd 06/27/2011

If you now look at displacements, you'll see that the artificial translation was recovered. I didn't place a threshold on distance - just used the minimum one. This could be added easily enough. (replace objects without a match to nan to maintain order)

Best Answer

Related Solutions

MATLAB: Write program by For loop, how can I stop

MATLAB: Connect bwlabeled components

Related Question