setting up the problem set

DrawPrincipalEllipse.m

-function DrawPrincipalEllipse(c,r,szColor)
-
-x = c-mean(c);
-y = r-mean(r);
-
-[U,S,V] = svd([x,y],'econ');
-xyp = [x y]*V;
-
-
-r1 = (max((xyp(:,1)))-min((xyp(:,1))))/2;
-r2 = (max((xyp(:,2)))-min((xyp(:,2))))/2;
-
-theta=linspace(0,2*pi,1000);
-cXp = r1*sin(theta)' ;
-cYp=  r2*cos(theta)';
-% plot(xyp(:,1),xyp(:,2),'.g');hold on
-% plot(cXp,cYp,'.k')
-cXY = [cXp cYp]/V;
-cXY = [cXY(:,1)+mean(c) cXY(:,2)+mean(r)];
-
-plot(cXY(:,1),cXY(:,2),szColor);
-

GetRandomEllipse.m

-% randomEllipses
-function [pts ptsEdge] = GetRandomEllipse(mu)
-
-
-r1 = 10+40*rand();
-r2 = 10+10*rand();
-
-theta = 2*pi*rand();
-
-cv = [r1^2 0;0 r2^2];
-
-[pts]=genEllipse(mu,cv);
-ptsEdge = genEllipse(mu,cv,1,0);
-
-pts=pts-repmat(mu,[size(pts,1) 1]);
-ptsEdge = ptsEdge-repmat(mu,[size(ptsEdge,1) 1]);
-rotate = [cos(theta) -sin(theta);sin(theta) cos(theta)];
-
-pts=pts*rotate;
-ptsEdge=ptsEdge*rotate;
-
-pts=pts+repmat(mu,[size(pts,1) 1]);
-ptsEdge=ptsEdge+repmat(mu,[size(ptsEdge,1) 1]);
-
-bw=false(1000);
-idx=sub2ind(size(bw),round(ptsEdge(:,2)),round(ptsEdge(:,1)));
-
-bw(idx)=1;
-
-bw = imfill(bw,'holes');
-[r c]=find(bw);
-pts=[c r];

GetRandomEllipseImage.m

@@ -150,3 +150,40 @@ for i=1:length(idxPtsOverlap)
         end
     end
 end
+
+
+function [pts_xy, bw]=genEllipse(mu_xy,cv,c,bFill)
+
+if nargin<3
+    c=1;
+end
+
+if nargin<4
+    bFill=1;
+end
+bw=false(1000);
+    
+theta=linspace(0,2*pi,1000);
+[evec, eval]=eig(cv);
+[eval, order] = sort(diag(eval),'descend');
+evec=evec(:,order);
+
+a1 = sqrt(eval(1))*c*evec(:,1);
+a2 = sqrt(eval(end))*c*evec(:,2);
+
+x = (mu_xy(1) + a1(1)*sin(theta) + a2(1)*cos(theta));
+y = (mu_xy(2) + a1(2)*sin(theta) + a2(2)*cos(theta));
+x = max(x, 1);
+y = max(y, 1);
+x = min(x, 1000);
+y = min(y, 1000);
+
+pts_xy=[x' y'];
+idx=sub2ind(size(bw),round(y),round(x));
+bw(idx)=1;
+
+if bFill
+    bw = imfill(bw,'holes');
+    [r c]=find(bw);
+    pts_xy=[c r];
+end
\ No newline at end of file

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));
-% imP=imP-eps;
-figure;surf(imP);hold on
-

genEllipse.m

-%
-% genEllipse(mu,cv,c,bFill)
-% return the ellipse with center mu, covariance cv and isocontour c
-% if bFill is 0, only returns the ellipse boundary points
-%
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%
-% Copyright (c) 2016, Drexel University
-% All rights reserved.
-% 
-% Redistribution and use in source and binary forms, with or without
-% modification, are permitted provided that the following conditions are met:
-% 
-% * Redistributions of source code must retain the above copyright notice, this
-%   list of conditions and the following disclaimer.
-% 
-% * Redistributions in binary form must reproduce the above copyright notice,
-%   this list of conditions and the following disclaimer in the documentation
-%   and/or other materials provided with the distribution.
-% 
-% * Neither the name of PixelRep nor the names of its
-%   contributors may be used to endorse or promote products derived from
-%   this software without specific prior written permission.
-% 
-% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
-% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
-% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
-% DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
-% FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
-% DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
-% SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-% CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
-% OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
-% OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
-%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-function [pts_xy, bw]=genEllipse(mu_xy,cv,c,bFill)
-
-if nargin<3
-    c=1;
-end
-
-if nargin<4
-    bFill=1;
-end
-bw=false(1000);
-    
-theta=linspace(0,2*pi,1000);
-[evec, eval]=eig(cv);
-[eval, order] = sort(diag(eval),'descend');
-evec=evec(:,order);
-
-a1 = sqrt(eval(1))*c*evec(:,1);
-a2 = sqrt(eval(end))*c*evec(:,2);
-
-x = (mu_xy(1) + a1(1)*sin(theta) + a2(1)*cos(theta));
-y = (mu_xy(2) + a1(2)*sin(theta) + a2(2)*cos(theta));
-x = max(x, 1);
-y = max(y, 1);
-x = min(x, 1000);
-y = min(y, 1000);
-
-pts_xy=[x' y'];
-idx=sub2ind(size(bw),round(y),round(x));
-bw(idx)=1;
-
-if bFill
-    bw = imfill(bw,'holes');
-    [r c]=find(bw);
-    pts_xy=[c r];
-end
\ No newline at end of file

go1Ellipse.m

-
-
-im = GetRandomEllipseImage(1);
-figure;imagesc(im)
-hold on
-[r c]=find(im);
-DrawPrincipalEllipse(c,r,'.m')
\ No newline at end of file

goCluster.m

-im = imread('Ellipse2.tif');
-% imagesc(im)
-
-[r c]=find(im);
-obj=gmdistribution.fit([c,r],2);
-
-idx = obj.cluster([c,r]);
-idx1=find(1==idx);
-idx2=find(2==idx);
-hold on
-plot(c(idx1),r(idx1),'.r');
-plot(c(idx2),r(idx2),'.g');
-
-DrawPrincipalEllipse(c(idx1),r(idx1),'.k');
-DrawPrincipalEllipse(c(idx2),r(idx2),'.m');
-
-

goMakeGMM_PS.m

-load('Experiment\SWgeo\2.mat')
-[mm idx]=max([rgExperiment.scoreGMM]);
-Experiment=rgExperiment(idx);
-save('Best2GM.mat','Experiment')
-
-ptsCombined=vertcat(Experiment.pts{:});
-idxCombined=sub2ind(size(im),1000-ptsCombined(:,2),ptsCombined(:,1));
-bw=zeros(1000,1000);
-bw(idxCombined)=1;
-bw=bw(470:550,450:540);
-figure;imshow(bw)
\ No newline at end of file