#Importação PREPS


#Set workspace-----------
setwd('/home/gandra/preps')


rm(list=ls()) ## Removendo as variáveis
gc()
packages<-c('reshape','gdata','rgdal','lubridate')

package.check <- lapply(packages, FUN = function(x) {
  if (!require(x, character.only = TRUE)) {
    install.packages(x, dependencies = TRUE)
    library(x, character.only = TRUE)
  }
})


#Embarcacoes------------------
load('embarcacoes_final.Rda')
e<-efim
# ef<-e[which(e$modalidade=='Arrasto de Fundo' | 
#               e$modalidade=='Espinhel Pelágico de Superfície' |
#               e$modalidade=='Rede de Espera' |
#               e$modalidade=='Cerco'
# ),]
# nrow(e)-nrow(ef)


#Apenas as novas artes de pesca
emb<-unique(efim$cod_embarcacao)

#Tracks----------------
l<-list.files(recursive = T,pattern = ".txt")

d<-read.delim2(l[90], header=T,sep=',',quote="")
d<-d[0,]

i=100
for(i in 1:length(l)){
  time<-Sys.time()
  df<-read.delim2(l[i], header=T,sep=',',quote="")
  #Convertendo tipos de dados
  df$nom_embarcacao<-as.character(df$nom_embarcacao)
  df$utc<-as.POSIXct(df$data_hora_us, "%Y-%m-%d %H:%M:%S" )
  df$velocidade<-as.numeric(as.character(df$velocidade))
  df$lat_dec<-as.numeric(as.character(df$lat_dec))
  df$long_dec<-as.numeric(as.character(df$long_dec))
  df$distancia<-as.numeric(as.character(df$distancia))
  df$velocidade<-as.numeric(as.character(df$velocidade))
  colnames(df)[5] <- "lat"
  colnames(df)[6] <- "lon"
  #Apenas algumas colunas
  df<-df[c(3:6,10:12,17)]
  #somente embarcacoes das modalidades selecionadas
  df<-df[df$cod_embarcacao %in% emb,]
  #somente na area de estudo
  # df<-df[df$lat<=-23 & df$lat>=-38 & df$lon<=-43 & df$lon>=-55,]
  df<-merge(df,e[,c('cod_embarcacao','modalidade')],by=c('cod_embarcacao'))
  d<-rbind(d,df)
  rm(df)
  print(i)
  print(Sys.time()-time)
}

save(d,efim,file = 'preps_brasil.Rda')