3d pdf

draw3dPDF.m

+
+
+
+
+
+
+im = imread('Ellipse2.tif');
+imagesc(im)
+
+[r c]=find(im);
+obj=gmdistribution.fit([c,r],2);
+
+idx = [1:prod(size(im))];
+[r c]=ind2sub(size(im),idx);
+
+probXY=obj.pdf([c',r']);
+imP = reshape(probXY,size(im));
+
+figure;surf(imP);hold on
+
+% imrgb = ind2rgb(im,[0 0 0;1 1 1]);
+% imagesc(imrgb)
+[r c]=find(im);
+plot(c,r,'.r')
\ No newline at end of file