# aug 26

setwd("/Users/jameshanley/Documents/work/Espionage")

source("ExtractLinesFromLargeFileFunction.R")

result = ExtractLinesFromLargeFileFunction("jupiter.ps",20000,10000)


LENGTH=unlist(result[3]) ; LENGTH ; summary(LENGTH)

n.l = length(LENGTH) ; # n.l

X = result[1][[1]] ; Y = result[2][[1]] ; length(X)

x.min=unlist(result[4]) ; y.min=unlist(result[5])
x.max=unlist(result[6]) ; y.max=unlist(result[7])
c(x.min,y.min,x.max,y.max)

plot(x.range,y.range,col="white", 
 xlim=c(x.min,x.max),ylim=c(y.min,y.max) )
 
for(i in 1:length(X) ) lines(X[[i]],Y[[i]] )

x.range=c(400,490) ; y.range = c(605,680)

plot(x.range,y.range,col="white", xlim=x.range, ylim=y.range)
for(i in 1:length(X) ) lines(X[[i]],Y[[i]] )

# focus on tick marks 
plot(x.range,y.range,col="white", xlim=x.range, ylim=y.range)
for(i in 1:270 ) {
	L = LENGTH[i]	
	xa=0; if (X[[i]][1] > X[[i]][L]) xa=1
	ya=0; if (Y[[i]][1] > Y[[i]][L]) ya=1
	text(X[[i]][1],Y[[i]][1], toString(i), cex=0.5,adj=c(xa,ya) ) 
	lines(X[[i]],Y[[i]]) ; 
   text(X[[i]][L],Y[[i]][L], toString(i), cex=0.5,adj=c(xa,ya) ) 
}

X.sd = as.vector(sapply(X,function(x) sd(x)))
Y.sd = as.vector(sapply(Y,function(x) sd(x)))

dot = (Y.sd==0 & X.sd==0) ; sum(dot)

vertical.line =  (Y.sd > 0 & X.sd == 0 ) ; 
n.jumps=sum(vertical.line) ; n.jumps

i.jumps = (1:n.l)[vertical.line]

i.dot = i.jumps = (1:n.l)[dot]

t.dot =  ( unlist(X[i.dot]) - X[[i.y.0]][1] ) / x1

n.t = table(1+round(t.dot,1))