leitorHidr.R

#### leitor pdo_term e pdo_oper

agoraa = now()
## funcao que transforma valores expressos no tipo caracter em tipo numérico
transf <- function(gf)
{
  aux = as.numeric(gf)
  
  if(is.numeric(aux) && length(aux))
  {
    
    for(i in 1:length(aux))
    {
      if(is.na(aux[i]))
      {
        a8 = strsplit(gf[i],split = "")
        if(is.integer(which(a8[[1]] == ",")) && length(which(a8[[1]] == ",")) == 1)
        {
          
          a8[[1]] = a8[[1]][-which(a8[[1]] == ",")]
          
        }
        if(is.integer(which(a8[[1]] == " ")) && length(which(a8[[1]] == " ")) == 1)
        {
          
          a8[[1]] = a8[[1]][-which(a8[[1]] == " ")]
          
        }
        a10 = paste0(a8[[1]],collapse = "")
        aux[i] = as.numeric(a10)
      }
      
    }
    
  }
  
  # caso nao haja valor algum em tal premio, por o número 0
  if(is.numeric(aux) && length(aux) == 0)
  {
    aux = 0
  }
  
  gf = aux
  
  return(gf)
  
}

tirarEspaco = function(aux)
{
  for(i in 1:length(aux))
  {
    a8 = strsplit(aux[i],split = "")
    a8[[1]] = a8[[1]][-which(a8[[1]] == " ")]
    aux[i] = paste0(a8[[1]],collapse = "")
  }
  
  return(aux)
  
}

library(readr)
library(tidyverse)
library(lubridate)
library(rapportools)

raiz = "C:\\Estágio\\dessem_semana"
setwd(raiz)

meses = c("Janeiro","Fevereiro","Marco","Abril","Maio","Junho","Julho","Agosto","Setembro","Outubro","Novembro","Dezembro")
earm_max = c(202691.938,19897.047,51602.07,15164.887) ## fonte ONS - energia agora - reservatórios

### Verificar quais os decks do dessem estão na pasta indicada

ds = list.files(getwd(), pattern = "DS_ONS")

a = strsplit(ds, split = "")
rv = NULL
data = NULL
for(i in 1:length(a))
{
  rv[i] = paste(a[[i]][15:17],sep = "",collapse = "")
  dia = paste(a[[i]][19:20],sep = "",collapse = "")
  mes = paste(a[[i]][8:9],sep = "",collapse = "")
  ano = paste(a[[i]][10:13],sep = "",collapse = "")
  if((rv[i]=="RV0")&&(as.numeric(dia)>23)&&(mes!="01"))
  {
    mes = as.numeric(mes)-1
  }
  if((rv[i]=="RV0")&&(as.numeric(dia)>23)&&(mes=="01"))
  {
    mes = 12
    ano = as.numeric(ano)-1
  }
  
  data[i] = as.character(as.Date(paste(ano,mes,dia,sep = "-")))
  
}
ds = cbind(ds,data,rv)

dadosHidr = NULL
t1 = now()
for(cont in 1:nrow(ds))   
{
  arq2 = list.files(getwd())
  unzip(ds[cont,1],junkpaths=TRUE)
  arq1 = list.files(getwd())
  arq = list.files(getwd(), pattern = ".dat")
  aux = setdiff(arq1,arq2)
  
  hidr = read_lines("pdo_hidr.dat")
  
  marc = grep("IPER",hidr)
  
  matriz = NULL
  
  for(i in marc:(length(hidr)))
  {
    matriz = rbind(matriz,strsplit(hidr[i],split = ";")[[1]])
  }
  
  matriz = cbind(matriz,rep(ds[cont,2],nrow(matriz)))
  
  dadosHidr = rbind(dadosHidr,matriz)
  
  
  file.remove(aux)
  
  print(ds[cont,2])
}
t2 = now() - t1

#load("dadosHidrDessem_18OUT_24OUT.RData")

dadosHidr = dadosHidr[c(1,which(dadosHidr[,6]=="  99 ")),]

#save(dadosHidr,file = "dadosHidrDessem_21OUT_18OUT.RData")

for(j in 1:(ncol(dadosHidr)-1))
{
  dadosHidr[,j] = tirarEspaco(dadosHidr[,j])
}


#filtros
usinas = as.matrix(dadosHidr[-which(duplicated(dadosHidr[,4])),4])
usinas = usinas[-1]
subsis = as.matrix(dadosHidr[-which(duplicated(dadosHidr[,5])),5])
subsis = subsis[-1]
periodo = as.matrix(dadosHidr[-which(duplicated(dadosHidr[,1])),1])
periodo = periodo[-1]
datas = as.matrix(dadosHidr[-which(duplicated(dadosHidr[,ncol(dadosHidr)])),ncol(dadosHidr)])

a =intersect(dadosHidr[,1],periodo[49:length(periodo)])

for(i in 1:length(a))
{
  dadosHidr = dadosHidr[-which(dadosHidr[,1]==a[i]),]
}
periodo = as.matrix(dadosHidr[-which(duplicated(dadosHidr[,1])),1])
periodo = periodo[-1]

marcs = list(usinas, subsis, periodo)

## geracao termeletrica por subs individualizada
gerSubs = list()
for(i in 1:length(marcs[[2]]))
{
  gerSubs[[i]] = dadosHidr[which(dadosHidr[,5]==marcs[[2]][i]),]
}

## geracao termeletrica por subs a cada instante
gerTot = list()
for(i in 1:length(gerSubs))
{
  matriz = matrix(nrow = length(datas), ncol = length(periodo))
  
  for(idata in 1:length(datas))
  {
    mm = which(gerSubs[[i]][,ncol(gerSubs[[i]])]==datas[idata])
    if(sum(is.empty(mm))==0)
    {
      aux = gerSubs[[i]][mm,]
      a = NULL
      for(per in 1:length(periodo))
      {
        matriz[idata,per] = sum(as.numeric(as.matrix(aux[which((aux[,1]==periodo[per])&(aux[,6]=="99")),33])))
      }
    }
  }
  
  row.names(matriz) = datas
  gerTot[[i]] = t(matriz) ## linha = semi-hora, col = dia
}


## geracao individualizada
gerInd = list()
for(i in 1:length(usinas))
{
  xx = dadosHidr[which(dadosHidr[,4]==marcs[[1]][i]),]
  
  y = as.matrix(xx[-which(duplicated(xx[,ncol(xx)])),ncol(xx)])
  x = as.matrix(xx[-which(duplicated(xx[,1])),1])
  
  matriz = NULL
  for(idata in y)
  {
    a = NULL
    for(per in x)
    {
      a = c(a,sum(as.numeric(as.matrix(xx[which((xx[,1]==per)&(xx[,ncol(xx)]==idata)&(xx[,6]=="99")),33]))))
    }
    
    matriz = rbind(matriz,a)
  }
  
  row.names(matriz) = y
  gerInd[[i]] = t(matriz) ## linha = semi-hora, col = dia
}



## plots

#sh = read.csv("entrada_semihora.csv")
semih_7d=rbind("00:00:00",as.matrix(read.csv("entrada_semihora.csv")))
colnames(semih_7d)="semi-hora"
semih_7d=as.vector(semih_7d[1:48,])
datarep=rep(ds[,2],each=length(semih_7d))
Semihora=rep(semih_7d,length(datarep)/48)



ttm <- paste(datarep,Semihora,sep = " ")
tts <- strptime(ttm,format="%Y-%m-%d %H:%M:%S")

#ts = 1:length(tts)
#tts = tts[which(as.numeric(ts)%%2==1)]


for(i in 1:length(gerTot))
{
  plot(tts,as.numeric(gerTot[[i]]),col="white", main = paste0("GH - ",subsis[[i]]),
       ylab = "MWm", xlab = "Semi-hora")
  lines(tts,as.numeric(gerTot[[i]]))
  grid()
}

for(i in 1:length(gerInd))
{
  plot(tts,as.numeric(gerInd[[i]]),col="white", main = paste0("GH - ",usinas[[i]]),
       ylab = "MWm", xlab = "Semi-hora")
  lines(tts,as.numeric(gerInd[[i]]))
  grid()
}



agoraa = now() - agoraa ## 12 minutos para execucao