polarize.r <- function(x,y,x0=0,y0=0){
	r <- vector()
	for (i in 1:length(x))
		r[i] <- sqrt((x[i]-x0)^2+(y[i]-y0)^2)
r}

polarize.t <- function(x,y,x0=0,y0=0){
	t <- vector()
	for (i in 1:length(x)){
		t[i] <- atan((y[i]-y0)/(x[i]-x0))*180/pi
		if (x[i]-x0<0) t[i] <- t[i]+180
		if (x[i]-x0<0 & y[i]-y0<0) t[i] <- t[i]-360}
t}

library(vowels)
setwd("~/Fabricius")
dat <- read.table("rp_data.txt",header=T)

normed <- norm.wattfabricius(dat)
normed$r <- polarize.r(-normed[,5],-normed[,4],-1,-1)
normed$t <- polarize.t(-normed[,5],-normed[,4],-1,-1)

dat[,"F1/S(F1)"] <- normed[,"F1/S(F1)"]
dat[,"F2/S(F2)"] <- normed[,"F2/S(F2)"]

dat$r <- normed$r
dat$t <- normed$t

write.csv(dat,"rp_dat_wf_polar.csv")