[R] off-topic: affine transformation matrix
dylan.beaudette at gmail.com
Tue May 29 21:56:57 CEST 2007
Thanks for the prompt and clear reply! The simplicity of the solution may have
been why I initially overlooked this approach...
The results look convincing (http://188.8.131.52/~dylan/temp/affine.png),
now I just need to verify that the output from coef() is in the format that I
need it in.
l <- lm(cbind(nx,ny) ~ x + y, data=g)
(Intercept) 6.87938629 5.515261158
x 1.01158806 -0.005449152
y -0.04481893 0.996895878
## convert to format needed for affine() function in postGIS?
(Intercept) x y
nx 6.879386 1.011588063 -0.04481893
ny 5.515261 -0.005449152 0.99689588
note that the format that I am looking for looks something like the matrix
defined on this page:
On Monday 28 May 2007 15:18, Prof Brian Ripley wrote:
> Isn't this just a regression (hopefully with a near-zero error).
> coef(lm(cbind(xnew, ynew) ~ xold + yold))
> should do what I think you are asking for. (I am not clear which
> direction you want the transformation, so choose 'old' and 'new'
> On Mon, 28 May 2007, Dylan Beaudette wrote:
> > This may sound like a very naive question, but...
> > give two lists of coordinate pairs (x,y - Cartesian space) is there any
> > simple way to compute the affine transformation matrix in R.
> > I have a set of data which is offset from where i know it should be. I
> > have coordinates of the current data, and matching coordinates of where
> > the data should be. I need to compute the composition of the affine
> > transformation matrix, so that I can apply an affine transform the entire
> > dataset.
> > any ideas?
> > thanks in advance!
Soils and Biogeochemistry Graduate Group
University of California at Davis
More information about the R-help