-
Notifications
You must be signed in to change notification settings - Fork 0
/
R_readcdss.R
108 lines (67 loc) · 3.31 KB
/
R_readcdss.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
library(here)
source(here::here('libraries.R'))
# Parques eólicos ---------------------------------------------------------
Hoy<- as.character(lubridate::today())
list_hoy<- readRDS(file = paste0(here::here('Data/Espana/'),"Espana_",Hoy,".RDS"))
#Lubian
Longitud_Parque=-6.84653
Latitud_Parque=42.04514
Lubian_list<- Cortar_datos(list_hoy = list_hoy,
Longitud_Parque = Longitud_Parque,
Latitud_Parque=Latitud_Parque)
if(!file.exists(paste0(here::here('Data/Parques/Lubian/'),"Lubian_",Hoy,".RDS"))){
saveRDS(Lubian_list, file = paste0(here::here('Data/Parques/Lubian/'),"/Lubian_",Hoy,".RDS"))
}else{print(paste0("Hoy ya se han guardado este archivo: ", "Lubian_",Hoy,".RDS"))}
#El Cerro
Longitud_Parque=-3.64
Latitud_Parque=42.87
ElCerro_list<- Cortar_datos(list_hoy = list_hoy,
Longitud_Parque = Longitud_Parque,
Latitud_Parque=Latitud_Parque)
if(!file.exists(paste0(here::here('Data/Parques/ElCerro/'),"/ElCerro_",Hoy,".RDS"))){
saveRDS(ElCerro_list, file = paste0(here::here('Data/Parques/ElCerro/'),"ElCerro_",Hoy,".RDS"))
}else{print(paste0("Hoy ya se han guardado este archivo: ", "ElCerro_",Hoy,".RDS"))}
#La Sia
Longitud_Parque=-3.57
Latitud_Parque=43.14
LaSia_list<- Cortar_datos(list_hoy = list_hoy,
Longitud_Parque = Longitud_Parque,
Latitud_Parque=Latitud_Parque)
if(!file.exists(paste0(here::here('Data/Parques/LaSia/'),"/LaSia_",Hoy,".RDS"))){
saveRDS(LaSia_list, file = paste0(here::here('Data/Parques/LaSia/'),"LaSia_",Hoy,".RDS"))
}else{print(paste0("Hoy ya se han guardado este archivo: ", "LaSia_",Hoy,".RDS"))}
#La Belesar
Longitud_Parque=-7.6854
Latitud_Parque=42.72343056
LaBelesar_list<- Cortar_datos(list_hoy = list_hoy,
Longitud_Parque = Longitud_Parque,
Latitud_Parque=Latitud_Parque)
if(!file.exists(paste0(here::here('Data/Parques/Belesar/'),"/Belesar_",Hoy,".RDS"))){
saveRDS(LaBelesar_list, file = paste0(here::here('Data/Parques/Belesar/'),"Belesar_",Hoy,".RDS"))
}else{print(paste0("Hoy ya se han guardado este archivo: ", "Belesar_",Hoy,".RDS"))}
## Convertir a listas por localizaciones y con toda la serie temporal
Lubian_lolat<- lon_lat_df_ls(Lubian_list)
ElCerro_lolat<- lon_lat_df_ls(ElCerro_list)
LaSia_lolat<- lon_lat_df_ls(LaSia_list)
LaBelesar_lolat<- lon_lat_df_ls(LaBelesar_list)
##Convetir componentes U y V a formato WS y WD
Lubian_lolat<- lapply(Lubian_lolat, uv_transformation)
ElCerro_lolat<- lapply(ElCerro_lolat, uv_transformation)
LaSia_lolat<- lapply(LaSia_lolat, uv_transformation)
LaBelesar_lolat<- lapply(LaBelesar_lolat, uv_transformation)
#Lo de los mapas esta a la espera de poder instalar JAVA
# Mapeos de puntos y demas ------------------------------------------------
#Seteamos el tamaño del mapa, para ello habrá que elegir
n=max(Datos2$lat)
s=min(Datos2$lat)
e=max(Datos2$lon)
w=min(Datos2$lon)
#Fijamos incremento para hacer más grande el mapa
incr<- 0.0215
if(n > 0){n<- n + incr}else{n<- n + incr}
if(s > 0){s<- s - incr}else{s<- s- incr}
if(e > 0){e<- e + incr}else{e<- e + incr}
if(w > 0){w<- w - incr}else{w<- w- incr}
ul <- round(c(n,w),digits = 3) #Upper Left
lr <- round(c(s,e), digits = 3) #Lower Right
download_maps(ul, lr, maptyp = "bing", res=15)