function [rc,cc] = singspotcent(A) %SINGSPOTCENT returns centroid location for first image in input 3D matrix MB=double(A(:,:,1)); [rn cn]=size(MB); M=MB-mean(mean(min(MB))); r=1; c=1; sumIc=0; sumIr=0; while c<(cn+1) sumpIr=0; i=1; while i<(rn+1) sumpIr=sumpIr+i*M(i,c)^2; i=i+1; end sumIr = sumIr + sumpIr; c=c+1; end c=1; while r<(rn+1) sumpIc=0; i=1; while i<(cn+1) sumpIc=sumpIc+i*M(r,i)^2; i=i+1; end sumIc = sumIc + sumpIc; r=r+1; end sumI=sum(sum(M.^2)); rc=sumIr/sumI; cc=sumIc/sumI; end