problem with fixed.angle?

Fri, 19 Apr 2013 19:31:55 -0700


----- Forwarded message from Fabio Machado <macfa...@gmail.com> -----

Date: Thu, 18 Apr 2013 20:53:03 -0400
From: Fabio Machado <macfa...@gmail.com>
Reply-To: Fabio Machado <macfa...@gmail.com>
Subject: problem with fixed.angle?
To: morphmet@morphometrics.org

Dear List,

I was trying to use the fixed.angle function in geomorph, and there seams to be a problem. The function deforms the resulting shape, even in the example given in the package. The problem seams to be in the centering on the rotation landmark, as you can check with the following code:

data(plethspecies)
plot(plethspecies$land[,,1],asp=1, xlim=c(-1,1))
points(plethspecies$land[,,1]-plethspecies$land[1,,1],asp=1, col="red")

fixed<-fixed.angle(plethspecies$land,art.pt=1,angle.pts=c(5,6),rot.pts=c(2,3,4,5))
points(fixed[,,1],asp=1, col="blue", pch=4)


Apparently, the subtraction of the articulation point follows "downwards" through the column and recycle through the entire configuration. I've altered the code of the function as following

A[, , i] <- A[, , i] - A[art.pt, , i]

to

A[, , i] <- t(t(A[, , i]) - A[art.pt, , i])

And this seams to solve the issue, but I'm not sure if I'm altering anything essential for the function or even if the problem I'm having is not due to anything else.

Any input is appreciated.

Best,

-----------------------------------------------------
Fabio Andrade Machado
Laboratório de Evolução de Mamíferos 
Departamento de Genética e Biologia Evolutiva
Universidade de São Paulo
Rua do Matão, trav.14, nº277, 
Edificio Minas Gerais 
+55 11 3091-8758
+55 11 982-631-029
skype: fabio_a_machado




----- End forwarded message -----



Reply via email to