How about thresholding followed by regionprops() to ask for the centroid?
binaryImage = real(fftImage) > threshold;
labeledImage = bwlabel(binaryImage);
measurements = regionprops(labeledImage, 'Centroid');
How about imregionalmax() followed by regionprops() to ask for the centroid?
binaryImage = imregionalmax(real(fftImage));
labeledImage = bwlabel(binaryImage);
measurements = regionprops(labeledImage, 'Centroid');
You forgot to post your image by the way.
