Discover MakerZone

MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn more

Discover what MATLAB® can do for your career.

Opportunities for recent engineering grads.

Apply Today

How can i implement this code?

Asked by Tomas on 17 Oct 2012

Hi

I'm looking for a matlab code which estimates the two principal directions (axes) using principal components analysis (PCA)?? Then any an other algorithm can be applied to estimate the asymmetry of the region of image in terms of the two principal axis. I tried several times to implement it and I looked everywhere to no avail. Now I ask my question on this forum hoping to get the right answer. thanks

0 Comments

Tomas

Tags

Products

No products are associated with this question.

3 Answers

Answer by Sachin on 17 Oct 2012

1 Comment

Tomas on 17 Oct 2012

thank you for this link. I visited it but since I'm a beginner in matlab I have not got to use these concepts to draw the axes of the region in the segmented image and estimate the asymmetry.

Sachin
Answer by Tomas on 17 Oct 2012
Edited by Tomas on 18 Oct 2012

I tried to do this code but i need your help please

I2=segmentation(I);
[n m] = size(I2);
AMean = mean(I2);
co=I2*I2'; % Compute the covariance matrix (co)
% Compute the eigen values and eigen vectors of the covariance matrix
[eigvector,eigvl]=eig(co);
eigvalue = diag(eigvl);
pc = eigvector * I2;
plot(pc(1,:),pc(2,:))

2 Comments

Matt Kindig on 18 Oct 2012

What is your question? What about this code doesn't work?

Tomas on 18 Oct 2012

Hi

I'm looking for a matlab code which estimates the two principal directions (axes) using principal components analysis (PCA)?? Then any an other algorithm can be applied to estimate the asymmetry of the region of interest of the segmented image in terms of the two principal axis. I tried to do this code but it doesn't work so i need your help.

Tomas
Answer by Tomas on 18 Oct 2012
Edited by Tomas on 18 Oct 2012

Hi, Can any one please correct this code

A=imread('aeroplane silhouette.png');
bw=~im2bw(A,0.5); %Threshold and invert
subplot(1,2,1), imshow(bw,[]);
[y,x]=find(bw>0.5); %Get coordinates of non zero pixels
centroid=mean([x y]); %Get (centroid) of data
hold on; plot(centroid(1),centroid(2),'rd'); %Plot shape centroid
C=cov([x y]); %Calculate covariance of coordinates
[U,S]=eig(C)
m=U(2,1)./U(1,1);
const=centroid(2)/m.*centroid(1);
xl=50:450; yl=m.*xl+const
subplot(1,2,2), imshow(bw,[]); h=line(xl,yl); %Display image and axes
set(h,'Color',[1 0 0],'LineWidth',2.0)
m2=U(2,2)./U(1,2);
const=centroid(2)/m2.*centroid(1);
x2=50:450; y2=m2.*x2+const
h=line(x2,y2);
set(h,'Color',[1 0 0], 'LineWidth',2.0)

you find the images here

http://img15.hostingpics.net/pics/679738aeroplanesilhouette.png

http://img15.hostingpics.net/pics/171229plane.png

2 Comments

Matt Kindig on 18 Oct 2012

Before asking us to correct the code, you should tell us what's wrong with the code. From the second image, it seems that the gray lines lie along what I would expect to be the principal components.

Tomas on 18 Oct 2012

These image explained what i should to do but this code didn't plot the two axes

Tomas

Contact us