function centroiddata=centroidplot(M,N)
%a function to read the image matrix M and plot the centroids of each plot
%on the image
H=M(:,:,N);
cd /opt/EDTpdv/hartmann/src/
centroiddata = centroid_image(H);
cd /opt/EDTpdv/JKmatlab/

v=centroiddata.current_centroids;
r=v(:,1);
c=v(:,2);
display_dalsa_image(H);
hold
scatter(r,c,'CData',[1 0 0],'Marker','o')
hold
end