Let's say I have some coplanar points in R^3. I want to describe these points by coordinates on the plane they share, in R^2. This way I can do some tests faster. How should I do this?

**1**

# Project points to subspace

###
#2
Crossbones+ - Reputation: **2140**

Posted 15 January 2014 - 09:53 AM

Let's say I have some coplanar points in R^3. I want to describe these points by coordinates on the plane they share, in R^2. This way I can do some tests faster. How should I do this?

Take any 3 distinct points A, B, and C. Take the cross product of (B - A) and (C - A) to get the normal, N, of the plane. Let (B - A) be one of your axes (normalize it). Take the cross product of N and (B - A) to get another orthogonal vector in the plane as the other axes (normalize it). Then project your points (relative to A) in your new axes.

-Josh

**Edited by jjd, 15 January 2014 - 09:59 AM.**

--www.physicaluncertainty.com

--linkedin

--irc.freenode.net#gdnet

###
#4
Crossbones+ - Reputation: **2140**

Posted 15 January 2014 - 10:33 AM

Thanks! Am I right when I think that the final projection could be a matrix multiplication by the 2x3 vectors we obtained?

Yes, once you have subtracted A from all of the points.

-Josh

--www.physicaluncertainty.com

--linkedin

--irc.freenode.net#gdnet

###
#5
Crossbones+ - Reputation: **19960**

Posted 15 January 2014 - 10:36 AM

If you have an orthonormal basis of your subspace, you can find the coefficients of the projection by simply computing the inner product of the vector to be projected and the basis vectors. So yes, the computation boils down to a matrix multiplication that turns three numbers into two.

In order to get an orthonormal basis for the sub-space, you can use the Gram-Schmidt process.

[EDIT: I am not saying anything new, just rephrasing what others have said, in the hopes that my description might be easier to understand to some.]

**Edited by Álvaro, 15 January 2014 - 10:38 AM.**

###
#7
Crossbones+ - Reputation: **19960**

Posted 06 February 2014 - 04:54 PM

I'm wondering if there's a way to avoid having to subtract the subspace's origin before doing projection, by encoding the origin in the matrix. Is it possible using homogenous coordinates?

Yes, you can encode a translation into the matrix by adding a row and a column and by appending a "1" to the coordinates of the point. This is related to homogeneous coordinates, but I am not going to try to describe the details here.

###
#9
Crossbones+ - Reputation: **19960**

Posted 06 February 2014 - 05:46 PM

Projection(P - O) = Projection(P) - Projection(O) (because the projection is linear)