In short, you can use thresholding, bwlabel, and regionprops to find the objects and their centroids. You can use rectangle() to place a circle in the overlay centered at the centroid of each object.
You might also look into the hough() demos - do a search on MATLAB Central, the newsgroup, and on the Answers forum.
