# 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))