library('geomorph')
rm(list=ls())

source('ShapeRotator.R')


#SIMPLE.ROTATION
radioulna <- readland.tps("XXX.tps", specID = "ID")
humerus <- readland.tps("XXX.tps", specID = "ID")



land.a=XX
land.b=XX
land.c=XX 

land.d=XX
land.e=XX 
land.f=XX 

angle= 90

radioulna_t <- translate(radioulna, land.a)
humerus_t <- translate(humerus, land.d)

data.1 <- radioulna_t
data.2 <- humerus_t


arms90 <- simple.rotation(data.1= radioulna_t, data.2= humerus_t, land.a= 9, land.b= 1, land.c= 17, land.d= 52,land.e= 19, land.f= 107, angle= 90)

plot.rotation.3D(arms90, radioulna_t, humerus_t, 1) 


writeland.tps(arms90, file="Arms_90degrees", scale = NULL, specID = TRUE)