(00:02) muy bien nuevamente buenos días a todos y todos Muchas gracias Por quienes se conectaron muy puntual estamos ya con los dos ponentes afortunadamente así que bueno vamos a arrancar esta segunda sesión eh de los talleres de formación intermedios y vamos a dejar en el uso de la palabra a Héctor Villalobos adelante Héctor Gracias qué tal Muy buenos días a todos este vamos a empezar hablando hoy de de un poquito de análisis de datos temporales
(00:32) nr entonces voy a compartir mi pantalla para que me vayan siguiendo un segundito a ver nada más me confirman por favor si se ve Sí Héctor se ve bien okay Bueno recordarán ustedes que lo que nos platicó ayer Emilio las instrucciones para acceder a esta plataforma del Jupiter Hub entonces este nos vamos a ir metiendo ahí el enlace por si no lo ven es este de aquí que dice.
(01:13) 212c Cloud Entonces voy a ir entrando ahí nos aparece la pantalla de login Y en este caso vamos a escoger como opción de servidor el r no con r Studio otra cosa que también es importante recordar ahorita en lo que arranca esto es que ahí mismo en donde donde estaba hace un momento Déjenme entrar de nuevo en lo que abre inter conecta Ah ah vamos al material de Emilio de ayer otra vez aquí mismo estaban las instrucciones para clonar Eh Pues el el repositorio de todas las
(02:15) lecciones que vamos a ver nada más lo lo pongo muy rápido pero Yo supongo que ya lo hicieron No había más que presionar el enlace que está aquí no para para jalar todo el repositorio donde vamos a a ocupar ahorita el el material que está entonces voy a esperar un momentito para que me digan Si todos están ahí no sé si ya les terminó de abrir a mí todavía no estamos esperando pero si alguien tiene problemas en estas en esta etapa por favor Avísenme para poder
(02:52) este hacerlo no Bueno entonces nadie me dice nada Supongo yo que todo está bien Voy a hctor Perdona Héctor te voy leyendo algunos mensajes del chat que no sé si tenga que ver con el ingreso directo como preguntadas pero pregunta Cristian tolosa que dónde puede acceder a estos o consultar esos tutoriales bueno es que esas instrucciones se les pasaron ya les voy a copiar aquí el el enlace en el chat a ver
(03:30) lo veo comenta que ya lo logró ya lo logró Cristian Gracias Héctor s Emilio también está Emilio también está respondiendo Gracias em en mi caso también ya ya está realizado Okay Okay gracias gracias a todos Eh bueno Entonces miren ustedes pueden ya sea desde desde su página web bueno su su navegador de internet meterse a la en la parte de talleres intermedios al el día 17 que es hoy y en el enlace que dice material de datos temporales enr ahí pueden ver si
(04:06) ustedes seleccionan el archivo datos temporales en rmd pueden ver qué es lo que vamos a hacer ahorita es lo mismo que vamos a ver en r Studio no pero aquí está digamos la versión pues que todo el mundo puede consultar sin tener que tener acceso al al Jupiter Hub como ustedes en este momento lo están haciendo solamente eso como referencia vamos a dejarlo por ahí pero lo vamos a trabajar desde acá entonces y ya todos
(04:32) clonaron su repositorio les debe aparecer en La pestaña de acá abajo del R Studio ya V Bueno ahí tienen ustedes diferentes pestañas pero si seleccionan la que dice files les debe aparecer la carpeta de talleres intermedios y hicieron el proceso de clonado a ver creo que Hans tiene la mano levantada a ver adelante Sí disculpe este la clase anterior hicimos un Clon de los es intermedios pero he visto de que han hecho una actualización necesitaríamos hacer un Pull o algo para actualizar la
(05:05) carpeta o eh Sí yo creo que habría que hacerlo de nuevo no Emilio si nos puedes apoyar ahí con que vuelvan a presionar el enlace efectivamente solo tienen que presionar el enlace de nuevo bien entonces eh les decía yo en la pestaña de files si se van a talleres intermedios van a ver esta exactamente el mismo contenido que Les acabo de mostrar en el navegador de internet entes ahorita vayan hasta la carpeta dector y el archivo que vamos a trabajar
(05:39) es el rmd no el md que hace un momento les mostré Ese md es markdown es para visualizarlo exclusivamente en el navegador aquí ocupamos el otro el rmd ya sea que lo lo hagan de aquí o como normalmente lo hacemos en el menú de file del R Studio iríamos a open file vemos las mismas carpetas de las que estoy hablando ahorita 17 de octubre Héctor y seleccionamos el archivo rmd es el que yo ya tengo abierto en este momento igual me espera aquí un momentito para asegurarme de que todos
(06:14) están ahí bueno no hay no veo manos levantadas en el chat hay algo a ver qué dice Ah pero es un comentario de Emilio bueno el chat lo voy a dejar ahí pues ya les pediría por favor apoyo a Alejandra Emilio si me si me avisan si hay algo por ahí relevante no okay Bueno entonces no parece ser que no hay problemas hasta acá eh vamos a empezar entonces con con con esta sesión de datos temporales miren algo de lo de lo más importante que para mí
(06:54) es relacionado con este tema que creo que hay que abordar es entender el Cómo se almacena esta información no las fechas y las horas eh Y eso es sumamente importante porque es bien común yo que trabajo aquí con estudiantes de de maesta y doctorado que cuando vienen con sus datos después traen una mezcla de cómo se o sea inconsistencias en Cómo capturaron la información no algunos por ejemplo ponen una columna para el año otra para el mes otra para
(07:25) el día pero a veces lo invierten a veces ponen el mes con letra a veces con número Entonces eso es un verdadero caos no porque cuando lo importamos a r y lo queremos ya pasar a un formato formal de fecha pues tenemos muchos problemas Entonces algo de lo que está explicado aquí en esta breve introducción es justamente eso no que tenemos que ser consistentes básicamente los formatos utilizados pues son como este
(07:52) que está acá no dos dígitos para el día dos dígitos para el mes dos dígitos para el año y como decía yo en una misma columna no no mezclar no voltear el orden poner primero el mes y luego el día etcétera no el separador que le puse Aquí un guion Pues también puede ser una una diagonal como acá ahí por ejemplo está con el mes abreviado eso también es válido pero de nuevo consistencia si ustedes ponen el mes con el nombre
(08:20) abreviado completo pues no se vale que a veces lo pongan en mayúsculas luego en minúsculas completo etcétera no entonces mucho cuidado con esos con esos de detalles antes de importar los nr y pues de los formatos más estándar digamos perdónenme el chiste Pero no ofender a nadie eh más utilizados son Día mes y año o bien año mes y día la forma que se usa en algunos países de habla inglés al poner primero el mes luego el día y luego el año es muy confusa y genera muchos problemas Entonces yo sugiero que
(08:52) la traten de evitar supongan que ya tenemos por ahí eh nuestra información capturada y lo lo que sea no cuando importamos datos de ese tipo en r eh Hay dos formas digamos tenemos imaginen que tenemos un archivo as un archivo de texto que ya fue No sé procesado tal vez en Excel o en otro o en una base de datos real y eh cuando lo importamos con rable que es la función que se utiliza para importar datos de archivos sasi esa columna va a pasar
(09:26) como una clase de de texto no va a aparecer como carácter eh En cambio cuando lo hacemos con Excel es probable que lo reconozca ya como fecha siempre y cuando el formato haya sido consistente no entonces esos son dos diferencias importantes en la manera en la que ustedes pueden importar los datos no archivos sasi o archivos Excel directamente eh Y bueno pues una vez dentro de r vamos a tener clases no en r los objetos que se que se importan las función
(10:00) los datos etcétera se les asignan clases dependiendo de cuál es su contenido entonces en el caso de las fechas la más eh sencilla digamos de todas es la clase fecha no la clase day y en esto para ilustrarlo Pues aquí les estoy poniendo un ejempl generando algunas fechas así al azar voy a presionar el triangulito por acá eh el triangulito verde de este lado espero que ya todos sepan esa parte si no por favor me interrumpen y me dicen yo aquí generé con la función pa
(10:35) una cadena de cuatro fechas pegando el año 2024 con 09 10 11 12 para indicar el mes y después los días así no continuos no 13 18 25 31 etcétera separados por un guion y pues aquí podemos ver el resultado No eso pues ya para nosotros parece una fecha no efectivamente el es que si yo le pregunto a r ver de qué clase es ese vector me va a decir Ah pues ese vector es de clase carácter Entonces no lo puedo usar yo en gráficos porque no no significa nada en términos
(11:11) de una fecha no entonces necesito convertirlo a a ese formato que les decía y lo vamos a hacer A través de la función as date esta función lo que va a hacer es tratar de de reconocer qué parte de ese texto es el año cuál parte es el mes cuál parte es el día eh pero imagínense No si ustedes intercambian el lugar del mes y el día en diferentes posiciones en ese vector pues la función va a fallar aquí pues como todo es consistente no hay problema generamos el
(11:43) el el la conversión sigue viéndose igualito que antes no no cambió nada esto es exactamente lo mismo que tenía acá arriba pero ahora al preguntarle la clase me dice Ah bueno la clase de ese objeto es fecha y y eso ya me da muchas ventajas no ya lo puedo utilizar para procesar datos Entonces cómo es que el r va almacenar esta información si yo le presiono aquí o le utilizo la función oncl me va decir esos
(12:15) numerotes no Y esos números son el número de días transcurridos desde el primero de enero de 1970 no solamente para que sepan Cómo se almacena Entonces lo podemos comprobar miren yo voy aquí a dar una fecha esta fecha del origen eh Que estoy mencionando ahorita el primero de enero del 1970 y si lo resto de mis fechas pues me da exactamente los valores que me dio onclass con eso compruebo que está almacenado en r como el número de días
(12:46) transcurridos no y bueno vamos a avanzar por aquí un poquito eh yo con la función sis Time puedo consultar Cuál es la la hora de la computadora en la que trabajando en este caso fíjense yo hice el ejemplo en mi máquina aquí en en en México en la paz y ahora que lo lo estamos corriendo acá pues va a tomar la hora del servidor donde estamos trabajando no Entonces si se fijan aquí me dice pues que estoy el día 17 a las 14:18 en horario utc no que
(13:21) es el tiempo universal coordinado que no es el horario de mi de mi computadora local yo en el en lo que les puse acá adelante es que eh y vuelvo a utilizar la función as date con este objeto que es otra clase de fecha que ahorita voy a mencionar adelante me lo va a me va a descartar la la parte que se refiere a la hora no ahí lo lo pueden ver ustedes lo está descartando nada más que aquí yo para el ejemplo que estaba manejando antes pues le puse mi zona de tiempo no en este
(13:54) caso como es utc pues podíamos podríamos haberlo puesto así no como utc eso es importante porque eh si ya estuviéramos digamos en un horario aquí local en el que la hora hubiera estado Ya pasado el día 17 en utc Pues en lugar de darme acá el más más bien me hubiera dado el día 17 cuando ya en utc sería el 18 por eso es importante eso Y ustedes pueden conocer las zonas de tiempo si usan esta función que se llama hson Names no Y ustedes la presionan les va a dar todos los códigos de las zonas de tiempo de
(14:33) todo el mundo ento Sería cuestión nada más de Buscar el de ustedes y pues ya lo pueden poner por acá Okay entonces esa clase date solamente es para fecha hay otras dos que son yo creo que las más utilizadas que tiene que ver con un estándar internacional se llama posix no Entonces está el formato posix Lt y eh vamos a ver cómo es que se almacena este miren Yo aquí estoy generando con la función strp Time le estoy pasando
(15:03) una cadena de texto Esta es una cadena de texto no est entrecomillado es texto eso que no quede duda Pero bueno está en un formato que no es muy muy ambiguo hasta el año mes día y una hora no hora minutos y segundos separado por dos puntos Entonces yo aquí en format le puedo indicar Cómo es que está eso es decir yo pude haberlo puesto en en el otro formato que les decía hace un momento puedo cambiar aquí le pongo 17 y acá le pongo 10 ahí evidentemente el 17 no se
(15:34) refiere al al mes verdad porque no hay mes 17 entonces para que esta conversión funcione Miren la voy a hacer ahorita y me dice na porque porque la conversión falló porque yo le dije que en la segunda posición estaba el mes y no es cierto no puedo haber mes 17 entonces Tendría que haber hecho el cambio Aquí y ahora si tiene que funcionar ven entonces realmente puede estar capturado Como sea pero siempre de
(16:04) manera consistente y que ustedes conozcan estos códigos que yo estoy usando acá para decirle qué parte corresponde a qué cosa entonces eh Muy brevemente el porcentaje y mayúscula significa el año con toda la centuria no es decir 20 24 o 1985 si fuera el año sin la centuria O sea si yo tuviera aquí solamente de 24 entonces tendría que poner una y minúscula y ese sería mi año sin la centura y de todas maneras la conversión
(16:38) funciona Okay el porcentaje d evidentemente es para el día el porcentaje M es para el mes pero el mes como número también puedo poner el mes como texto no como abreviado o completo como sea eh Y la hora pues porcentaje H porcentaje m para minutos porcentaje s para segundos no Y la zona de tiempo Entonces ustedes pueden ver todos esos códigos que acabo yo de de de poner consultando la ayuda de la función
(17:08) stpr Time eso ya lo saben estamos en un taller intermedio Entonces no me no me detengo mucho pueden ustedes ver ahí en la ayuda todos los códigos que se utilizan no porcentaje a minúscula mayúscula por ejemplo es el porcentaje B Es para el mes abreviado mayúscula es para el nombre completo en el local de su computadora etcétera no todo eso lo pueden ustedes ver ahí para las diferentes posibilidades entonces este esta clase posix Lt cuando nosotros eh vemos cómo
(17:44) está guardada de nuevo usando la función onclass me va a dar un montón de cosas miren por ejemplo esto que tengo aquí es una lista no Y me dice Estos son segundos los minutos la hora El el mes perdón el día del mes no este es el el mes pero le puse octubre Y por qué me da 9 Ah pues porque la cuenta empieza desde cero no el enero es el mes cero y así sucesivamente y en el caso del año me da 124 Bueno eso es 124 porque en este caso el origen del año se considera 1900
(18:21) Entonces si suman 1900 má 124 pues les da el año no que yo le puse acá y el día de la semana el del año o día juliano también se le llama eh una banderita para ver si tenemos el el el horario eh Cómo se le llama el el horario de verano para el ahorro de energía no en este caso en México Ya no hacemos cambio casi en todo el país no eh la zona de tiempo etcétera no todo eso está es como se está almacenando esta información
(18:51) aquí muy bien eh También podríamos eh usar la función en lugar de esta que usé acá para convertir para generar digamos o convertir una cadena de texto a Pos Lt existe una función que se llama aspix Lt y aquí pues va de nuevo el formato fíjense como no le estoy indicando nada de de de zona de tiempo de formato porque no es ambigua la forma en la que esto está escrito entonces la función puede muy bien eh identificar eh las
(19:31) diferentes componentes de esta de esta cadena de texto y lo sabe poner bien no de nuevo Pues habría que eh tener cuidado en cómo está Cómo se le está pasando no Entonces fíjense como pues me vuelve a dar un un un objeto de clase posix Lt No aquí con Class le preguntamos la clase y el otro el otro formato existe eh o la otra clase es posic c es muy parecido el la manera de operar miren aquí la función es muy parecida la la voy a ejecutar ya no no
(20:05) Me detengo mucho ya es lo mismo que hace rato lo que sí tenemos que saber acá es cómo está guardado a diferencia del anterior en este me va a dar un numerote Miren el número que me da esa fecha Y eso porque esto se almacena como el número de segundos eh No años segundos transcurridos desde el primero de enero de 1970 comprobación muy fácil me voy a crear un objeto que le llamo orig de origen donde le pongo pues ese no el primero de enero de 1970 a la c horas y
(20:41) a ver ejecutamos y me dice no la diferencia de tiempo entre el el objeto psx que es el que hice yo acá Arribita este de acá pxz y el origen en segundos pues es este numerote que es justamente este de acá entonces todo esto es importante porque a veces nos pasan datos en donde el tiempo bien puede venir Así guardado e como ese número y nosotros tenemos que saber cómo convertirlo a una fecha porque si no no va a tener sentido entonces creo que me
(21:18) tengo que dar un poquito de prisa porque ya son 26 minutos creo que s alcanzamos vamos a a proceder con un ejemplo de datos de información de de de Copérnico no sé si han escuchado de esa plataforma creo que lo vamos a ver en otros tutoriales más adelante eh No me voy a tener mucho ahí ahorita pues un archivo que ya está en en las carpetas de del repositorio eh Y ustedes pueden aquí con este código leer no exactamente leer Digamos
(21:51) como que tener acceso a un archivo net cdf que está ahí en en esa carpeta de datos de del tutorial que estamos siguiendo ahorita y pueden consultar por ejemplo con esto eh Cuáles son las horas que Ay perdón Es que aquí se me vinieron todas y por eso tenía acá abajo Head para solamente ver algunas me borrar lo vuelvo a ejecutar entonces fíjense esos numerotes que tenemos ahí esas son fechas no pero tenemos que conocer el formato en el que están pues para poderlo consultar entonces si yo tecleo acá en la la
(22:28) consola ncf que es donde estoy digamos creando la conexión a ese archivo puedo ver todos los detalles de lo que ese archivo contiene eh de nuevo no Me detengo mucho son datos de temperatura para diferentes latitudes y longitudes y la parte que me interesa ahorita es esta no donde dice el tiempo y me dice Ah bueno las unidades del tiempo son segundos desde el primero de enero del 70 Entonces con eso ya sé cómo procesar
(22:55) esta información no la la versión correcta sería la función aspos ct esa ese vector de de valores con esos números Ay insisto yo en jalarlos todos si son muchos y se me va perdón eh Y eh le doy Ah pues el origen es este en horario utc y ahora sí ya los tengo convertidos a fecha no O en este archivo que estamos trabajando ahorita tengo datos diarios desde el primero de enero del 2019 hasta el 31 de diciembre de 2020 no son son dos años de
(23:36) datos diarios Pero bueno eso ahorita nada más en lo que respecta a la conversión digo esto Aunque parece un poco laborioso podemos fácilmente codificarlo nosotros en funciones y es el ejemplo que vamos a hacer ahorita eh miren tuvimos algunos detalles con con la instalación de algunos paquetes Entonces les voy a pedir que en este blo en el que estamos ahorita le quiten los comentarios a esas dos líneas y los ejecuten por favor
(24:04) estas dos no el que dice library dep Tools y install github H villala Team por favor hagan esas dos cosas yo ya lo hice No lo voy a hacer y los voy a volver a comentar pero por favor para que lo demás nos funcione ustedes ejecuten esas dos líneas y me espero un momentito ahí a ver si alguien me puede decir por favor si está funcionando o no Bueno no veo manos levantadas No veo nada en el chat Supongo que está
(24:54) funcionando Héctor Perdona tengo un mensaje directo de Ronald abrego dice yo no he podido entrar a Ocean hack week 2 y2c Cloud en mi usuario rona brego o del github entonces no sé si te lo mando a ti directamente o si se lo mando igual a Emilio para que qui nos está apoyando también Sí pues yo creo que ahí Si pudiera Emilio ayudarnos a ver qué a orientarlo a ver cómo puede acceder Sí claro eh voy a a confirmar
(25:40) pero básicamente todo el mundo está todo menos una persona está ya inscrita e Bueno gracias está la información un momento eh En lo que pegó no sé si sea eso pero dice ghob ron abrego @g coma ail.com no sé si haya habido un error con su correo voy a confirmar también e sugiero que Don Héctor continúe y yo voy a a investigarlo y comunicarme por el chat Así que Hans por favor eh ponerle atención al Ah no perdón Ronald al al
(26:19) chat Gracias Emilio bueno no no me dijeron si hubo problemas o no Yo espero que no los haya habido ayer yo hice la instalación y funcionó rápido y sin problema no entonces este paquete sa ya tiene digamos les podría enseñar ahorita el código de la función pero hace lo básicamente lo que dijimos identifica el origen de de los datos el formato en el que están las fechas y lo convierte ya un formato apropiado además de importar evidentemente los datos de temperatura no Entonces con esta función que dice
(26:54) r me va a importar esos datos podemos ver preguntarle a ver de qué clase de datos son me dices clase satín y tecleando su nombre lo voy a hacer aquí directamente en la consola puedo ver qué es lo que hay en esos datos no son datos de temperatura la variable se llama t o porque es temperatura potencial del mar en grados centígrados diarios Ya lo habíamos visto eso una resolución espacial de 9.
(27:20) 2 kmr y ahí están las dimensiones de los datos Qué son esas dimensiones el 61 se refiere al número de latitudes que tenemos en esos en esa información 60 es el número de latitudes es el tamaño de la matriz y 731 sería el número de días que tengo no porque Repito está aquí lo pueden ver en Data ranges está desde el primero de enero del 19 hasta el 31 de diciembre del 2020 y además está para dos profundidades diferentes no 0.49 m y 1.
(27:54) 54 entonces ya está todo importado ya las fechas está como deben estar Y pues vamos a ver entonces un poquito Cómo están almacenados esos datos el código este que dice str me muestra la estructura Este es un objeto de clase S4 donde yo puedo con el con el símbolo este de arroba extraer partes de lo que está ahí no veo que tengo un vector de longitudes uno de latitudes un arreglo de datos no que las dimensiones que ya comentamos la información que me dijo de la variable etcétera y el que me importa
(28:29) ahorita no el el periodo digamos de datos que están guardados ahí entonces yo para extraerlo decía uso la función ara se llaman slots cada una de esas partes no entonces del objeto este que importé tetao con @ period voy a extraer parte de los periodos pero en esta información Siempre aunque sea la misma en ocasiones tengo un tiempo de inicio y un tiempo final no sé por qué pero así lo tienen aunque sean datos diarios así vi eh entonces con Head puedo ver cuáles
(29:00) son los primeros datos ya veo que eso ya es un formato de fecha formal y con Tail puedo ver pues los los finales no ahí para que vean que se convirtió correctamente y bueno un plot de ese información pues me daría no rápido una visualización del tipo de datos que tenemos no repito datos de temperatura en grados centígrados para esta región del Golfo de California en México Y bueno pues ahora vamos a extraer datos de un de un pixxel no porque aquí estoy viendo solamente un día pero yo
(29:34) tengo 731 Entonces con el siguiente código le digo a ver Quiero un dataframe para extraer el Pixel que está en la longitud 110 y la latitud 26 no ent con esto aquí la función es extract Point ex pts me extrae esta información de acá no vean me dice ahí las dimensiones de lo que acabo de sacar es un renglón y 1468 columnas por qué es esto tengo 731 días pero tengo dos profundidades Entonces sería 731 por 2 pero aún así no me da el número completo no Ah bueno lo
(30:15) que pasa es que en este objeto digamos que va a guardar lo voy a hacer acá mejor para que se vea completo te guarda digamos identificador del punto que yo seleccioné el que yo elegí con las coordenadas x y donde yo dije quiero el dato aquí Pero puede ser que en ese en ese punto en particular sobre todo si hago clic en el punto en lugar de darle un valor como tal como yo lo hice ahorita no haya o sea el el dato está referido
(30:46) al centro del pixxel entonces en donde dice Long y lat me está diciendo En qué en qué punto Sí hay dato el punto más cercano donde hay dato y la de sería un control nada más para que vea yo que no está muy lejos de donde yo seleccioné y donde hay dato no en este caso no podría exceder 9.2 km porque es la es la resolución de los píxeles no entonces Bueno pues ahí a partir de de ahí tengo todos estos que serían ya los valores de temperatura para cada uno de los días
(31:20) terminaría primero los de la parte superficial y luego los de la siguiente profundidad entonces me interesan a mí los primeros y es lo que voy a hacer yo ordenándolos acá pero lo voy a hacer ahora tomando digamos del objeto original el periodo de inicio que es un formato de fecha formal y las temperaturas lo voy a eh combinar con las temperaturas extraídas pero a partir de la séptima columna para brincarme estas primeras seis que son informativas y me tengo que ir hasta el final de la
(31:52) primera profundidad es decir 731 días más se que son las columnas que me que tengo de exceso al principio no entonces ya lo tengo aquí arreglado como una tablita como un dataframe de fechas y temperaturas que podemos graficar y pues ahí tendría yo mi ciclo de en dos años no los dos ciclos anuales eh de temperaturas diarias en esta zona no con esto yo puedo trabajar este Híjole ya estoy en a las 7:37 y creo que tengo 40
(32:23) minutos verdad Bueno voy a voy a darme prisa con esta parte ya eh podemos explotar esta información de la fecha gracias al paquete lub en el cual yo puedo extraer de mi fecha la parte que corresponde al mes al día o al año que es lo que hace este código de acá no y lo estoy combinando en la misma tablita que ya tenía No aquí voy a repetir esto de y ahí tendría extraído digamos en columnas aparte el mes y el año porque esas las puedo yo usar con otras funciones de r como t apply para sacar
(33:03) por ejemplo el promedio por año no ahí lo ven ustedes eh le digo de la temperatura me interesa que me la eh digamos eh ordenes por año y quiero la media No ahí tengo medias por año o las medias por mes sin importar el año o las medias por año mes no también ahí ya tengo combinado Y si quisiera yo digamos graficar solamente el promedio mensual en lugar de usar la función t apply que me devuelve una matriz o una Ray un
(33:35) arreglo usaría la función aggregate la sintaxis es prácticamente la misma eh Y el arreglo es un poquito diferente aquí no vemos todos pero estarían digamos todos los meses del 2019 y luego serían los del 2020 con esto Pues podría tener un gráfico muy simplificado del otro porque ya no tendría los valores diarios tendría los promedios mensuales no no no Me detengo mucho aquí porque bueno el código ahí lo tienen y lo pueden ejecutar si
(34:06) ustedes lo desean no debe funcionar sin problema bien aquí es donde tengo un problemita porque el paquete Smooth no se no se logró instalar Entonces les voy a copiar ahorita en el chat un pedazo de código que es mi plan b y yo creo que con eso me voy a quedar porque no nos va a dar tiempo de terminarlo lo otro eh eh voy a copiar todo esto de aquí Se los voy a pegar aquí en el chat a híjole Creo que va a quedar como respuesta a algo pero no sé dónde hacer
(34:51) una nuevo lo voy a pegar ahí es lo que está hasta abajo Entonces eso que acabo de pegar ahí por favor en la parte donde dice library Smooth pueden ustedes seleccionar antes del plot y pegan ahí lo que acabo de copiar antes del plot lo del plot no lo vayan a borrar por favor entonces Bueno eso es porque la el paquete smoo no se puedo instalar no sé por qué entonces Bueno lo que hacía ese paquete pues era generar
(35:24) un promedio móvil que espero que todo el mundo sepa lo que es un promedio móvil qué me permite suavizar digamos No es simplemente aquí le estoy diciendo Yo que cada 30 valores cada 30 días me dé un promedio entonces va a tomar del día un al día 30 promedio del día 2s al día 31 promedio Y así sucesivamente hasta que termina y eso pues permite quitar un poquito de la variabilidad diaria que hay en los datos de temperatura y me
(35:50) muestra ya un suavizado de esos datos no eh Miren la la siguiente parte pues la voy a dejar el código está si instalaron lo que Ah bueno aquí también habr que hacer una instalación como la que hicieron con satín ustedes podrían correrla Y si hay preguntas pues ya las ponen por ahí en las diferentes opciones que tienen No ya sea el discord o bien en la plataforma de También estamos usando y podamos platicar por ahí si existieran dudas No
(36:29) Quisiera yo comer más tiempo Pues porque tenemos tres tres temas diferentes que ver esta mañana no Entonces lo voy a dejar ahí y pues estaré atento ahí a sus preguntas por si las hay Gracias Perdona Héctor no sé si me oyes Pero tenemos un tiempo para preguntas también eh o sea que ahora si quieren preguntar en directo pueden hacerlo tranquilamente y si no por discord como prefieran okay Bueno entonces este pues escuchemos a ver si hay
(37:03) pregunta se me apagó un monitor aquí alguien tiene alguna pregunta Antonio una pregunta era Era sobre cargar eh datos directa algún enlace para cargar directamente de copernicus eh No sin cargar en local sin cargar en local o sea directo Pues yo al menos no conozco ninguna Bueno sí hay por ahí algunas opciones pero el problema o al menos donde yo trabajo luego donde platico estos temas es que el internet No es muy bueno entonces No podemos
(37:49) confiarnos a estar jalando los datos en el mismo momento desde un servidor porque a veces pasamos minutos valiosos descargando datos y no podemos avanzar en lo demás entonces la manera en la que yo trabajo es descargo los datos en mi computadora y ya sobre eso trabajo no Entonces eso es más bien el enfoque que yo que yo prefiero no por esta cuestión les digo no en ocasiones por ejemplo con tengo alumnos que están trabajando con modelos de distribución de especies y tienen áreas enormes de de trabajo y cuando
(38:22) quieren descargarlos pues no se puede no hay que bajar los archivos individuales más bien y ya hacer la nosotros la la el procesamiento aquí en local no vale Vale gracias no hay alguien más bueno Yo nada más les les comentaría si si ustedes descargan o ya tienen el repositorio en su máquina local eh podrían ejecutar en local Y ahí sí van a poder instalar el paquete Smooth y no habr necesidad del Pequeño cambio que yo introduje ahorita eso eso fue Simplemente para salir del del paso con
(39:07) esta parte no Espero que les haya funcionado que haya corrido nadie dijo que haya tenido problema a ver Hans la lantada este esa función del paquete satín solo te funciona para los archivos en cdf de temperatura puede ser para cualquiera siempre que venga cnicus ese sirve para todos normalmente debería funcionar para todos los datos que están en Copérnico no las diferentes bases aunque la que hemos más hemos probado es la de El global Ocean reanálisis
(39:42) no también puede leer los de datos de la otra que es la de datos de que trae planton que trae oxígeno a veces los formatos cambian un poquito de los net cdf pero en general funciona bien para todas y hay en el paquete otras funciones para bajar o importar datos de Ocean color de productividad primaria de la Universidad de Oregon y los de Ultra alta resolución de temperatura superficial del mar de un km ya genial
(40:19) graciasa alguien más Cuánto es el tiempo para preguntas tenemos unos unos 15 minutos Héctor Ah sí sí sí así que bueno les invitamos a todas y todos aprovechando el tiempo de Héctor y su presencia que puedan plantear sus preguntas y también abrir su cámara por favor para que podamos interactuar de mejor forma todos pudieron instalar el el paquete satín nadie tuvo problemas con eso No creo que no si no desp yo sí tuve problemas eh Para
(41:09) instalar el paquete Ah s no sabemos Qué pasó porque Emilio lo había instalado en el el Jupiter Hub este ese y el de Smooth pero por alguna razón no dice que no están no no sabemos por qué s voy a investigar eso pero desgraciadamente esta semana continuará el problema hoy pues hoy Que estaremos trabajando en r esto no debería ocurrir están instalados eh No no se presentó ningún error durante la instalación pero en la práctica no funcionan no los reconoce entonces para espero que la
(41:50) mayoría de ustedes continúen el evento delat en noviembre y espero que parece entonces y habremos resuelto esto o en realidad para los que piensen regresar a este tutorial ahora y elat volver a la plataforma en línea si lo resolvemos digamos la semana que viene Entonces ya estará ya tendrán acceso a eso Bueno pues si no hay más preguntas pues tal vez podríamos pasar a a Marina si les parece eh Héctor un momento si no te importa para que yo pueda retomar
(42:42) podrías enseñarnos tu gráfica con las temperaturas y los años que yo he tratado de capturarla y no me ha dado tiempo y digo Mira pues empiezo por ahí que puede ser buen El problema es que no sé por qué se meag El Monitor secundario Ah y entonces no puedo ver lo que está ahí ni regresarme nada lo siento no no sé qué pasó Estoy aquí tratando de ver si ay ay ya ya creo que ahí va ahí va Dame un segundito creo
(43:14) que se Se aflojó acá el cable ya ahí va Ahí va un segundito estaba yo aquí intrigado porque qué pasa porque ya no veo nada y iba a tener que salir para dejar de compartir pero ahí va ahí viene ya está es la gracia del directo Sí a ver entonces a qué parte me di me regreso e esos gráficos tenías otro gráfico arriba con con los puntitos este este gráfico sí que nos quedemos un poco con la copla de que o sea con la idea de que Héctor nos está enseñando el promedio de temperatura en un punto en
(43:46) concreto para muchos años desde para para para diferentes meses y diferentes años vale quedemos un poco ahí y ahora yo retomo con esa con esa imagen y con ese gráfico retomo yo la parte de datos espaciales que al final está muy relacionado con todo lo que nos ha presentado Héctor porque la parte espacial y la temporal no podemos separarla entonces eh Bueno pues yo puedo empezar compartiendo pantalla vamos a empezar con a ver cómo hago yo
(44:18) esto a ver vale Yo voy a empezar con una presentación a ver compartir aquí a ver lo que no sé si verán el modo presentación si lo pongo modo presentación a ver están viéndolo ahora mismo la presentación Sí Marina ya perfecto genial pues entonces ahora en esta parte vamos a ver los datos espaciales que sería mediante una serie de paquetes de de gis de sistemas de información geográfica nr y posteriormente Cómo podemos usar este tipo de herramientas para cálculos de
(45:11) métricas climáticas y comoo esto por ejemplo se puede aplicar en temas de cambio climático para medio Marino en ecosistemas marinos que pues yo ahora mismo la verdad es que no sé muy bien Qué nivel tiene cada uno y entonces es difícil hacer es posible que las cosas que diga ya las sepáis O puede que no Entonces He decidido empezar por por lo más básico para que las personas que no que desconozcan Cómo funcionan los sistemas de información geográfica por lo menos puedan tener acceso a esta
(45:41) información y si ya lo sabéis Pues estupendo ratificar que que son conocimientos que ya conocéis entonces eh en sistemas de computación la información geográfica se puede usar para muchas cosas para combinar para manipular almacenar analizar datos producir otra serie de productos y todo ello se combina con sistemas de información geográfica y esto se puede computar a través de por ejemplo un lenguaje como r si utilizamos los paquetes adecuados en el caso de Héctor
(46:15) nos ha estado enseñando cómo trabajar con datos temporales de variables climáticas para utilizando r básico que suele ser un tipo de lenguaje dentro de r Y utilizando determinados paquetes por ejemplo él nos ha presentado lubridate que es una que es una es una función ampliamente usada para combinar fechas y poder y poder trabajar con ellas porque a fin de cuentas necesitamos poder convertir unos productos en otros y poder modificarlos es una de las varias
(46:49) cosas que nos ha enseñado Y en este caso utilizaremos otro tipo de paquetes en los datos espaciales tanto en en r como en argis eh realmente hace unos años se trabajaba en argis que es un programa de de sistemas de información geográfica de pago y ahora ya en tenemos en en abierto en sistemas libres como r se han creado paquetes que nos ayudan a trabajar sistemas de información geográfica de la misma manera simulando argis no Exactamente igual pero simulando este tipo de funciones Entonces cuando
(47:25) empezamos a trabajar con datos espaciales vemos que tenemos eh Por ejemplo puntos líneas y polígonos entonces eh Por ejemplo puntos pueden ser eh puertos ciudades puntos de muestreo donde hemos ido a muestrear las coordenadas de una localización geográfica concreta las líneas serían líneas que unen diferentes puntos pues por ejemplo la ruta de un barco o la trayectoria de un animal que se ha desplazado a veces se les pone un GPS a diferentes animales y y o a los barcos
(48:00) que también llevan un sistema de información geográfica donde nos indica En qué posición está ese barco pesquero por ejemplo y podemos ver la trayectoria de ese barco a lo largo de muchos días si vamos uniendo todos los todas las coordenadas o por ejemplo polígonos si enlazamos diferentes líneas podemos tener un área Marina protegida con diferentes niveles de protección o la zona económica exclusiva de un país tenemos eh diferentes eh diferentes
(48:31) H objetos Por así decirlo diferentes datos espaciales esto serían un ejemplo de puntos en un en un mapa mundial diferentes líneas que por ejemplo indicarían los ríos de una los ríos del mundo el anterior indicaba las capitales del mundo eh un polígono por ejemplo una serie de polígonos que indican todos los países a nivel mundial y cada punto cada línea y cada polígono está definido por sus coordenadas que en
(49:02) este caso las coordenadas son X e Y pero cuando trabajamos con sistemas de información geográfica son latitud y longitud entonces pues aquí tenemos puntos que tenemos la la la x y la i si si os fijáis entre paréntesis aparece 1 2 La la el uno sería la x el dos es la y de la misma manera para las líneas y para Los polígonos eh del mismo modo necesitamos todas esas coordenadas que conforman los que conformarían un polígono Qué sucede a nivel a nivel en
(49:41) oceanografía por ejemplo y en en investigación biológica no solamente tanto en Marino como en terrestre podemos utilizar lo que llamamos rasters que son eh el raster raster son una serie de ficheros que acumulan información que manejan información de muchos tipos y que la conforman de una manera que nos hace muy fácil manejarla por ejemplo en este caso tenemos un mapa eh raster de una determinada variable
(50:07) que en este caso eh yo creo que es temperatura los raster tienen los raster los vectores tienen diferente resolución en este caso un vector serían las coordenadas del punto las coordenadas eh para cada vector y el raster lo que hace es pasarlas a a un grid de malla entonces imaginamos que tenemos el mapa del mundo en un grid de malla y todo el mapa todo tu mapa está conformado por cuadrículas estas estas cuadrículas serían esto que vemos aquí a la derecha
(50:41) de cada cada cada mapa con cuadrículas está formado por un por un raster que se puede almacenar de muchas maneras o como o como imagen o como grd G que serían ficheros georreferenciados para poder abrir y exportarlos luego en por ejemplo en kugis en programas que que nos lo plote rápidamente y que podemos ver la localización de estas zonas e aquí no voy a entrar en más detalle porque luego cuando abramos el el código de r podremos entrar un poco más en
(51:14) detalle es muy importante tener en cuenta la diferencia entre extensión y resolución la extensión hace la función extent extensión hace referencia a todo tu mapa Cómo grande es el mapa que tú quieres eh analizar y la resolución hace referencia a cómo de Cuál es el tamaño de tu de tu grid Cuál es el tamaño de tu malla de información en este caso pues tenemos un un grid muy amplio y esta resolución se va haciendo O sea este grid se va haciendo más pequeño lo cual nos aporta más
(51:50) información de la de la zona por ejemplo el primer grid sería de grado uno que se llama y conforme se va haciendo más pequeñito ese ese tamaño de celda se va reduciendo el grid y el siguiente será de unos Yo diría que unos 0 0,5 porque se divide cada grid se va dividiendo en sub celdas todo esto es un poco a groso modo Porque si me extiendo mucho más me puedo tirar 5 horas hablando de este asunto Entonces esto es eh resumido Qué
(52:21) es la resolución de un fichero no voy a entrar en las eh proyecciones espaciales porque ya complej aría mucho esta esta esta sesión Pero bueno aquí podemos ver que cada cada píxel tiene un valor numérico con su localización y con el valor asociado de una variable por ejemplo entonces cada cada línea cada layer que se llama en inglés eh se puede acumular en un raster y acumular temporalmente y espacialmente en este
(52:57) caso por ejemplo lo que y hao acumular eh una serie de mapas de 2017 a 2020 que indicarían cada capa de de nuestro raster entonces eh una misma capa un mismo mapa lo podemos acumular en en en capas con la función stack y tendremos cada cada cada mapa tendrá una serie de valores por ejemplo para temperatura y cada una de esas capas referencia a una fecha en el caso de los datos que vamos a ver son datos yo creo que son mensuales Y luego podemos trabajar con ellos a nivel anual en el caso de los
(53:37) datos por ejemplo que nos enseñaba Héctor eran datos mensuales y en lo que hacía era hacer una gráfica para mostrarnos el promedio de temperatura para cada año o para cada mes en este caso es algo parecido lo que pasa es que lo tenemos todo almacenado en un mismo documento que sería lo que llamamos un brick que lo que haces con un brick es almacenar una serie de mapas en un mismo fichero y cada capa hace referencia a una fecha y a unas coordenadas en este
(54:11) caso pues tenemos aquí Estados Unidos porque para variar siempre lo tenemos que tener en todos los sitios pero Eh bueno son Bueno es es un ejemplo cualquiera pero cada año eh tenemos eh cada mapa sería un año respondería con un año y con una variable ahora veremos Cómo esto se aplica al abrirlo en en en r en este caso lo que haremos será abrir un fichero en sí que se llama un fichero netcdf que tiene que tú te lo descargas de copernicus Julia mañana nos enseñará
(54:44) Cómo hay que descargarse ficheros de copernicus que creo que lo enseñará directamente Cómo descargarlos directamente desde python que antes ha habido Alguien que ha preguntado por justamente por esto mañana podremos aprender Cómo descargar datos lo que no sé si nos lo enseñará en r creo que solo lo enseñará en python en r también es posible hacerlo de nuevo Yo por ejemplo en mi caso la mayoría de las veces lo que hago es descargármelo en mi
(55:10) ordenador en en local como indicaba Héctor que también hacía lo mismo y trabajar con sus ficheros locales en esto cada uno pues hace un poco lo que va pudiendo y lo que necesita en cada momento y eh en en este caso lo que haremos será abrir un un fichero de temperatura y lo convertiremos a brick Y entonces al pasar netcdf es un formato en sí es un formato que está ampliamente extendido cuando cuando trabajamos con datos de satélite la gran mayoría de
(55:40) ficheros vienen como ndf o en Sí porque son son ficheros muy Versátiles muy útiles ocupan poco espacio casi todo el mundo son compatibles con la gran mayoría de sistemas de programación entonces todo el mundo es importante saber cómo abrir un fichero cdf Y cómo trabajar con él en mi caso yo me he dado cuenta que abrirlo como brick eh Y abrirlo con con sistemas de información geográfica con este tipo de paquetes
(56:07) resulta muy sencillo y resulta muy muy fácil de visualizar y de entender qué hay dentro de un de un fichero eh NC que no son no no siempre son sencillos bueno Hay muchísimas fuentes de datos Online para para para información geográfica que necesitemos por ejemplo tenemos un mapa de temperatura de todo el Mediterráneo Yo yo trabajo con el Mediterráneo vosotros y vosotras no sé dónde estaréis pero Eh pues queréis toda la costa eh de Centroamérica o la costa de de un país
(56:44) concreto la costa de Guatemala pues podemos entrar en emodnet y Buscar la batimetría de tu región y descargarte la batimetría y las profundidades a las cuales está está hay que ver que hay disponible que no hay disponible pero a través de este tipo de servidores gratuitos podemos saber cuál es la profundidad Y si tú tienes una serie de puntos de muestreo eh de cualquier análisis por ejemplo o un puerto si combinas la localización
(57:14) del puerto que te interesa con una capa de batimetría podrás saber qué profundidades hay alrededor de ese puerto esto En el caso de batimetría pero lo tenemos para nivel terrestre NASA tiene también una serie de servidores la FAO para temas de pesca O sea hay hay hay muchísimos datos geográficos Y georreferenciados que ahora mismo están en abierto y si buscamos datos gratuitos de de sistema de información los encontraremos sin sin problemas no todo y no todo al máximo detalle pero hay
(57:52) mucho disponible ahora mismo qué paquetes Con qué paquetes vamos a trabajar ahora pues vamos a trabajar con paquetes raster SP y plot y raster ahora ya ha sido actualizado por Terra y Terra es un paquete que está siendo muy versátil y muy muy ampliamente utilizado también otro como sf que ahora mismo Está más eh se aplica más ggplot todo esto Si no conocéis los nombres no os preocupéis porque son al final los paquetes que que que combinan las
(58:23) funciones que nos interesan Y entonces Entonces al final es es importante conocer los paquetes porque son los que nos van a dar las herramientas para tener acceso y hacer todos estos cálculos que que necesitamos a ver y bueno ahora Nos pasamos a a r y vamos a ver eh Cómo cómo continuamos vale aquí a ver no sé si a ver a ver a ver ahora tendré que dejar de compartir a ver y compartiros Júpiter Júpiter Hub Vale ahora mismo estáis viendo mi mi ventana de a ver Uy
(59:06) no esto no ay que esto no sé cómo hacerlo a ver vale ahora mismo estáis viendo eh la el Júpiter a ver espera sí estaréis viendo el Júpiter al que tenéis todos acceso si entramos en el enlace que que ya hemos entrado el que ya entramos el otro día y volvemos a entrar hoy el de Ocean Hub week 2i 2c iCloud en este Hub os lo pongo por aquí por si acaso ya lo ha puesto antes Héctor Pero bueno lo pongo por aquí para que para que podáis entrar y aquí en lugar de python El otro
(59:45) día entramos en python en esta ocasión entraremos en r Studio al igual que hemos hecho durante la sesión de Héctor Uy Y ahora no me deja a ver A ver por qué no me deja bueno supongamos que me dejase vale si suponemos que me deja una una vez entramos tenéis que entrar a rstudio no sé si si queréis ir diciéndome si tenéis acceso o no aquí entraríamos en rstudio y aquí nos nos dirigiría Se os abriría eh rstudio en blanco y aquí
(1:00:20) tenéis que fichar tenéis que entrar clicar en Open file y ahora está pensando lo que normalmente debería salir es eh talleres intermedios y en talleres intermedios entrar en el 17 de octubre 17 de octubre Marina y ahí hay una carpeta que pone métricas climáticas y datos espaciales y abrir ese rmd que es un rmd un fichero rmd eso es lo que a mí no sé no sé por qué no me está dando me está dando error y no os lo puedo enseñar pero exactamente si entráis Open file no me
(1:01:00) está dejando Ahora si me contáis si os si alguien puede entrar me vendría muy bien saber si estáis pudiendo entrar veis aquí tenemos datos temporales de r que era el el fichero que habíamos abierto antes con Héctor y ahora abrirías el que pone sí que habéis podido vale lo ideal es que abráis el que pone Def 2 pone métricas climáticas eh Def y hay otro que es Def de definitivo dos que sería el más reciente pero yo no he podido abrirlo en el día de hoy Así que utilizaré el que pone
(1:01:33) solo Def es este de aquí si abrís eh Bueno aquí entramos en el hay un fichero también que es html que si lo abrís saldrá todo el todo el Script ya funcionando y podréis Ver todas las imágenes de todo lo que de todo lo que hay en este contenido pero ahora vamos a ir paso a paso En caso de que tengáis problema y por algún motivo no funcione si abrís el html podréis ver exactamente el resultado que da cada una de las funciones entonces aquí tenemos si le damos al
(1:02:07) Play Uy ahora no me está dejando Ayer me funcionaba perfecta no Marina a ver aquí sí decirme eh No sé veo que ti varias pestañas abiertas y y me parece confuso tal vez en realidad está desconectado lo que te sugiero es que nos va a tomar unos minutos estoy mandando un enlace al al chat que vayas a ese enlace y salirte y volverte a meter vale estaba ahí puedes ver que estaba desconectada my server le digo stop my Sí eso y ahora de nuevo estáis viendo
(1:02:50) para eliminar la confusión te sugiero que cerres las demás pestañas por lo menos las otras tres que dice Ah vale bor cierro todo lo anterior vale Ay lo acabo de cerrar todo me caches en la m demasiado literal al cerrar todo entonces le digo stop my server Start my server ahora r Studio Start te saliste no estamos viendo tu pantalla Ah no esto mi pantalla vale discp a ver la comparto de nuevo aquí vale dle Ok a eso como dije puedes confirmar ahora que funcion
(1:03:36) efectivamente cualquier cosa que puedas hacer para para confirmar que parce que se ve bien lo estáis viendo ahora mismo s vale a ver entonces aquí le damos a Open file ahora tenemos veremos talleres intermedios entramos en talleres intermedios entramos en el 17 de octubre Marina datos espaciales y aquí tenemos el el métricas climáticas Ocean hawick global Def Esto es lo que os indicaba que había hay otro colgado que es def2 pero vamos a trabajar ahora con Def vale con el fichero definitivo entonces aquí
(1:04:16) entramos y eh podemos ir dándole le damos al Play directamente Estos son chun cada uno de los huecos grises ese es un chun trmd Y si le damos al Play no tenemos que hacer nada más si le das al Play lo que hace es eh hace funcionar cada una de las líneas de código y en este caso se nos abrirían todas las eh todas las todos los todas las librerías en este caso si no os funcionase que no sé si os funcionará eh en principio funciona no hay problema
(1:04:53) si os da problema la librería de de boc me avisáis pero por el la única librería que os podría dar problemas es BC pero yo creo que el que que no está dando problemas el otro día a mí me entró dice que no está instalada Vale pues en este caso aquí he dicho otra opción si las anteriores no funciona en este caso le quitamos la la la almohadilla quitamos esta este símbolo de aquí y eh Y corremos aquí le hacemos funcionar el install package
(1:05:23) remote y no le decimos que no nos lo que no reempe porque si no tardaremos mucho y remotes y que nos estas dos líneas de fichero son las que hay que hacer funcionar vale que ahora mismo tienen una almohadilla hay que quitar la almohadilla y hacerlo funcionar entonces una vez quitamos esas dos almohadillas y hacemos funcionar estas dos No deberíamos de de tener problema lo que podemos hacer es poner aquí almohadilla en en esta en esta de aquí que pone library boook y ya no nos dará problemas
(1:05:58) y Le podremos dar al Play y funcionará entero entonces eh A ver seguimos por aquí vale no debería haber problemas en seguir haciendo funcionar el resto el resto de librerías que ya se han que ya están cargadas con lo cual aquí no vamos a perder tiempo pero cuando estéis trabajando en vuestro ordenador será importante que eh veréis que las librerías no las tenéis en en el ordenador Y entonces habrá que instalar cada una de esas librerías descargarlas
(1:06:31) instalarlas para que os funcionen esto es importante porque de pronto te encuentras con que no te funciona una librería y es porque hay que instalarla primero en este caso utilizaremos datos de copernicus pero no vamos a entrar como Julia nos va a enseñar mañana cómo se hace todo esto a través de copernicus en este caso no vamos a entrar a ver cómo se obtendrán los dat desde copernicus pero eh allí podremos acceder a datos de
(1:07:00) temperatura superficial y eh Y y Bueno aquí Tenemos aquí yo he dejado escrito como un tipo de Cuál es el fichero que vamos a utilizar y cuáles son algunos de los artículos que aplican este tipo de cosas que vamos a que vamos a ver ahora pero en este caso e a ver bajamos hacia abajo y le damos a esta a esta línea de código a la 61 que nos dice source Data que lo que hace es Descargar datos que ya están subidos a a jth datos de temperatura superficial que
(1:07:38) del satélite hadi Y entonces aquí tenemos un fichero que se llama Had sst y le damos a ver qué hay aquí dentro de Had sst aquí podemos ver como tenemos datos de 11 que en este caso sería un grado de resolución que era lo que veíamos anteriormente que es eh una resolución típica para mapas grandes mapas de mapas muy amplios por ejemplo mapas a nivel mundial si quisiéramos trabajar a nivel de kilómetros a escala mucho más pequeña trabajaríamos con resoluciones más
(1:08:16) Menores en este caso un grado hace referencia a un grado de de de la malla de latitud longitud del sistema de coordenadas qué corresponde con 111 km que está por aquí escrito en el Script corresponde con 111 km de de tamaño en un mapa vale en este caso eh A ver aquí podemos ver que que tenemos el mapa que la extensión del mapa es de 180 a 90 que eso quiere decir que es de todo a nivel mundial es un mapa mundial y
(1:08:50) vemos Que en en la función nombres nos dice que tenemos datos desde 1000 1955 y que los como vemos va cambiando el número del mes y tenemos datos mensuales a nivel mundial si le preguntamos qué clase es este este documento nos dice es un raster porque lo hemos abierto directamente como raster si si descargá demos un fichero en sí que también eh lo podemos ver si nos da tiempo y lo abrimos directamente como raster trabajaré del mismo modo en en este tipo
(1:09:28) de de objetos objetos rast aquí le preguntamos Cuántos nombres Cuáles son los nombres Dentro de este Dentro de este objeto y nos dice que van desde 1900 subimos hacia arriba y nos dice que va desde 1995 desde el 1 de enero del del 16 del 1 de enero de 1955 hasta eh diciembre de en este caso diciembre de 2010 que aquí tendríamos la información temporal de la que nos hablaba Héctor previamente si hacemos un cálculo rápido esto son 56 años y
(1:10:11) eh si tenemos por 12 meses tendremos 672 capas una cada una de ellas hace referencia a un mes entonces Si volvemos de nuevo al objeto Had J sst en las dimensiones tenemos 672 capas Entonces esto quiere decir que tenemos un raster stack apilado y cada una de las capas corresponde a un a un a un mes de información en este caso de temperatura superficial si le decimos ploteos la primera capa del fichero ji Pues nos enseñará eh A ver nos enseñará aquí tenemos eh la temperatura
(1:11:01) almacenada a nivel mundial en este mapa del mundo en grados centígrados si si queremos entonces aquí eh la información que tenemos es un paquete aquí podemos verlo un un dentro de un mismo objeto dentro de un único objeto Tenemos un montón de meses georreferenciados porque tenemos aquí la la latitud la longitud Aquí vemos eh la extensión de latitud y longitud 180 a 90 crs sería la proyección en el wg84 del
(1:11:35) pside 84 esta sería la proyección de coordenadas que utilizamos y esto quiere decir que tenemos un montón de capas almacenadas eh con el con la forma de un mapa mundial qué Qué sucede por ejemplo eh tenemos un digamos que tenemos un montón de tenemos este mapa de aquí a ver a ver a ver a ver [Música] qué tenemos este mapa de aquí este mapa hace referencia a la Gráfica eh a la Gráfica que nos había enseñado previamente Héctor Héctor nos
(1:12:15) enseñaba una gráfica en la que para cada fecha había un promedio de temperatura en este caso cada mapa está asociado a una fecha de manera que si cogiese una sola celda de cada uno de estos mapas y empezar a plotear y a decirle Por ejemplo si le decimos ploteos en vez de plotear noos O sea si empezás demos a jugar con este número de aquí podremos ir viendo podríamos plotear todos los días directamente Por
(1:13:12) ejemplo si hacemos solo plot eh Y no le decimos exactamente Cuál es la la la línea de código que queremos que nos pinte tendremos un mapa con todo todos los todas las capas de nuestro raster entonces aquí tenemos que esto es enero febrero marzo vale Para que veamos Aquí acaba en en abril porque no te lo puede dibujar todo te tendría que dibujar 672 672 mapas que están almacenados en esta documentación de aquí si cogemos estos mapas cada uno de estas capas y
(1:13:51) queremos que nos diga Cómo ha cambiado la Cuál ha sido la el cambio eh anual Cómo ha cambiado la temperatura a lo largo de toda nuestra serie de datos al fin de cuentas esto lo que es es una serie temporal de datos de temperatura desde 1955 hasta el año 2010 si plote o sea si dibujamos todas todos esos datos y trazamos una recta de regresión para ver cómo ha cambiado temporalmente y espacialmente la temperatura superficial
(1:14:24) del mar podemos evaluar Cuánto ha cambiado ha habido un aumento de la temperatura del mar o no sabemos que hay cambio climático Pero cómo vemos este cambio climático en en los datos que nosotros con los que nosotros trabajamos con los que nosotras Trabajamos para ello hay una función que a mí me ha parecido muy útil que es procede del paquete boc velocidad del cambio climático que lo ha desarrollado Jorge García molino un investigador y esta función se llama Sam series y
(1:14:56) esta función lo que hace es calcular el promedio anual pasando los datos de de días o meses a años en este caso por ejemplo esta esta sería la función y aquí le decimos nuestro objeto inicial que era el el fichero raster que teníamos que es Had sst aquí le decimos la fecha el año y el mes inicial de nuestra de nuestra base de datos no no de nuestra base de datos de los datos de los que nosotras tenemos interés en este caso por ejemplo estamos interesadas en trabajar con datos de
(1:15:34) 1960 hasta 2009 y queremos que comience a analizar los datos desde el 1 de enero hasta diciembre le decimos 19601 que corresponde con el 1 de enero 2009 a 12 a diciembre y tenemos que decirle también exactamente Cuál es la fecha de la serie de datos por ejemplo si si aquí lo estamos explicando de manera muy rápida pero si le decimos help que sería la la ayuda el interrogante Sam series lo he dejado por
(1:16:09) ahí en el fichero Pero bueno lo repasamos de nuevo se nos Abre una ventana donde nos indica nos da toda la información en inglés de Qué hace la función Sam series y aquí nos Explica qué tenemos que poner en cada uno de los en cada uno en cada coma cada objeto separado por una coma en este caso nos dice r Sam series r y r es un raster stack que en este caso es el fichero que tenemos Had sst que contenga series temporales de una
(1:16:42) variable climática en nuestro caso temperatura la p que sería el siguiente el siguiente objeto que tenemos que poner la siguiente información que le tenemos que dar a esta función y nos no dice nos pide el periodo para el cual queremos calcular serie de datos eh Gear 0 año inicial nos pide exactamente qué día empieza nuestra serie de datos si os acordaba Si volvemos al fichero Had sst nuestra serie de datos aquí en names empieza en [Música] 195.
(1:17:10) 09 166 a la x no le tenemos que hacer caso en este caso eh empieza el 1 de diciembre el 1 de enero de 1955 Mirad aquí me he dado cuenta que lo he puesto mal yo le he dicho uno y en realidad sería 16 luego le decimos el número de capas que tiene y para ello ponemos directamente en leyers que son el número de capas dentro de nuestro propio fichero la función x que sería la esto ya está dentro de la de esta propia función que es para calcularte las tasas de cambio de
(1:17:51) temperatura y queremos que la frecuencia que estamos introduciendo en nuestro fichero es mensual porque tenemos datos de enero febrero marzo y queremos que nos nos devuelva los datos en años entonces lo que estamos haciendo es hacer un promedio de datos mensuales a anuales para unos mapas espaciales Entonces estamos combinando temporal con con con espacial y Aquí hacemos funcionar esta esta línea de código la que está en
(1:18:25) 133 134 y le decimos Vale ahora dibujan con la función plot dibújalos eh dibújalos Qué hay en Dentro de este fichero y ahora en lugar de tener ficheros y y mapas de cada mes para cada año pasamos a tener eh cada un un mapa por cada año porque es el promedio anual Entonces pasamos de 1990 a 1900 1960 61 62 veis aquí tenemos eh cada uno de los mapas pero con esta información seguimos sin saber exactamente cómo ha cambiado la temperatura en los últimos eh en los
(1:19:10) últimos 60 años cuál en los últimos 50 50 años eh para para extraer más información de de todo esto si si alguien tiene algún problema voy rápido porque es que si no no hay mucho tiempo y no me da tiempo a a presentarlo todo pero si alguien va teniendo problemas a la hora de hacerlo funcionar o de entenderlo puede encender su micrófono entender su cámara y preguntarme o bien preguntar por el chat que tenemos a Emilio y a Marian que nos
(1:19:42) pueden que nos pueden y al resto de tutores que nos pueden ayudar a resolver dudas eh justo hay un mensaje de que dice por alguna razón no grafica error inl New no abierto bien parece que es el el p h el fichero Sí si el problema es el fichero hay que volver a la a la función source es importante que en la función source Data en esta línea de aquí de código de la 61 esta nos funcione bien que podamos Acceder al objeto Had y J sst que podamos plotear para
(1:20:27) luego seguir trabajando si si mira ahora ya se te ha cargado estupendo Este fichero sí que está cargado Entonces en ese caso deberías de poder decirle plot H sst un para plotear hasta ahí todo va bien genial genial genial Abigail y entonces ahora eh el problema Entonces sí que sí que has podido plotear cada uno de los de los días Eh os funciona si eso te ha funcionado debería funcionarte el el la función Sam series Sería bueno
(1:21:06) saber si os ha funcionado la función hadi la función de de plotear los datos que hay para hadi del fichero del del del objeto hadi os ha funcionado estas dos de aquí la línea de código 79 y 80 no justo lo que no estoy pudiendo es hacer el plot o sea todo lo demás funciona bien pero marca un error como a la hora de crear el eh la Gráfica e No ni ni el plot 1 todo lo demás me manda bien la clase los nombres Had aparec con todos sus datos solo plot plot es el que
(1:21:45) no Marina A veces lo que pasa con r Studio Supongo que pasa También aquí en la plataforma es que Si cierras mucho la parte derecha se puede mover el el tamaño de la parte derecha y la parte izquierda entonces si en la parte derecha no tiene suficiente espacio para plotear a veces da error como está haciendo ella Entonces si está muy cerrado aquí a la derecha A veces intenta mover hacia la izquierda esto a ver y hay que bajar hacia abajo porque
(1:22:11) en este fichero esto es un fichero RM down y en muchas ocasiones no se te dibuja en la sección en la cuadrícula de la derecha sino que se te dibuja más abajo entonces aquí hay un montón de información que yo he dejado para otro tipo por si nos daba tiempo a verlo pero si después de J sst si deslizamos hacia abajo aquí en la línea 118 es donde tenemos que ver los los los mapas puede ser que bajando hacia la 118 veáis el
(1:22:46) mapa pero en ese caso no le daría error no al darle error plot New por eso me ha recordado más a esto que tenga muy pequeño la parte de la derecha y has podido probar abig No no es eso tampoco Marina mar yo pienso Perdón aquí aquí en mi r Studio me dice que no tengo instalado map Tools Entonces ya le comenté milio porque es extraño que en tu caso no apareció y acá sí probablemente se deba a eso porque es
(1:23:18) posible que map Tools lo esté utilizando para los mapas y por eso no le aparezca claro eso sí puede ser eso sí puede ser sí eso sí que puede ser que haya un problema a la hora de las librerías Eh no os preocupéis Porque ahora mismo lo importante es dice que sí que tiene instalado map Tools comprueba que tengas todos los paquetes instalados eh Por favor todo el mundo si ha tenido problemas intentad comprobar si todas
(1:23:44) las librerías os os funcionan es no se me ocurre otro motivo porque la función plot es muy sencilla a veces puede ser ser que de problemas funciones de raster o de grid extra hay paquetes más complejos pero plot no suele no suele dar problemas prob va hacer como dice Marian abrir esta parte de aquí desplazando este este este eje para poder para que se nos dibujen aquí en la sección plot los mapas en lugar de que se dibujen en el propio en el propio Script si os parece Vamos a
(1:24:21) continuar para no parar toda la sesión pero en la próxima semana podréis volver a revisitar estas métricas mira ya le funcionó a Abigail estupendo estup tiene su mano levantada si es de ahorita o se quedó desde hace rato Ah perfecto Hans una consulta más que todo de concepto no me qued completamente claras trajar con estos archivos Debería ser un raster eh stack casi no se raster stack verdad sí no sé si alguien lo Oye yo
(1:25:05) prácticamente no lo oigo Ah no no se me escucha bien está un poquito mejor ahora sí este mencionaba que lo que me venías diciendo era que se tenía que trabajar con un raster stack Pero bueno la clase del objeto con el que vienes trabajando es un raster brick en este caso vendrían a ser conceptos es como que r los está trabajando de manera similar o es solamente por el paquete porque sabemos que raster brick y raster stack como objetos como tal r los toma de manera diferente los procesa de manera distinta
(1:25:35) pero cuando leí la ayuda también decía debe ser un raster stack pero ahí la clase del objeto que has medido metiendo es un raster brick Entonces no hay ningún problema con que sea un raster brick o un raster stack y en qué momento recomendarías trabajar con un raster bak con un raster stack propiamente eh sí es un es un objeto raster en este caso funciona bien si es un brick al por mi experiencia yo lo que veo es que en algunas ocasiones tenemos que usar raster stack eh de manera
(1:26:09) restrictiva pero en muchas otras nos permite usar un raster brick eh En ocasiones son restrictivos pero en este caso por ejemplo la función eh Sam series funciona bien con un con un raster brick y con un raster stack a fin de cuentas son son si haces brick lo que haces es un tanto la función raster como la función brick lo que hacen es acumular mapas que que tengan una una zona similar y y que tengan la misma proyección geográfica y la la misma
(1:26:44) Estructura de datos si queréis seguimos vale para ver un poco el resultado uno de los resultados finales que podemos obtener Entonces en este caso retomamos la función Sam series A ver no sé si he respondido a tu pregunta eh creo que sí Entonces continuamos el entonces aquí le metemos datos de de mensuales los obtenemos como año y lo que hacemos es obtener un fichero que son años los plote y podemos verlo podemos verlos aquí debajo que tardará un poquito
(1:27:28) e y eh A ver ahora está tardando en plotear pero bueno como ya lo hemos visto antes sí aquí tendremos cada uno de los mapas para cada año si le preguntamos slot names nos dirá Cuál es Qué información hay en cada uno de sus mapas Cuántas capas tiene en este caso tenemos 50 capas que corresponden con 50 años cuáles son las dimensiones en este caso nos dice eh 180 360 que son los grados de nuestro mapa global aquí
(1:28:03) podríamos eh dentro del dentro del Script he dejado una una opción B pero he puesto opción B no usar porque es para trabajar con datos de ficheros en sí descargados directamente de copernicus pero para que veáis cómo se trabajaría si te descargas los datos y los tienes como fichero en sí que es eh que es bastante sencillo de convertir y luego podríamos agregar estos datos en este caso como estamos trabajando con un grado no los agregamos pero podrían agregarse si nuestros datos
(1:28:38) fuesen de una resolución mayor Se podrían agregar y unir todos en en en celdas de nuestro grid mayores y aquí por ejemplo podemos eh calcular el promedio de todos nuestros mapas si aquí le decimos min y nuestro objeto aquí podemos y nuestro objeto que yo creo que es un raster brick si le preguntamos a ver qué objeto es grst que sería como el promedio anual es un raster en este caso es un raster stack nos ha pasado de un raster brick a
(1:29:13) un raster stack usando la función Sam series aquí le decimos promedio para que nos calcule el promedio de 50 años temperatura Y si le preguntamos que nos lo dibuje aquí tenemos cuál sería el promedio de temperatura desde eh desde el año aquí la fecha está mal puesta sería desde el año eh desde el 20160 al 2010 sería lo correcto que no sé por qué está así y ahora una función que a mí me parece muy útil es la función temp Trend que es tendencia temporal nos Calcula Cuál es la tendencia climática en este caso
(1:29:55) usando temperatura superficial de cada celda local entonces lo que hace es dibujar una recta de regresión linear para cada punto de nuestro grid de malla si le preguntamos eh Qué es la función tem Trend aquí en la en la función de ayuda vemos que lo que nos pide es un raster stack que sería el fichero que ya hemos obtenido y eh Y y bueno directamente le decimos podemos decirle el th es un es es el mínimo de observaciones que
(1:30:30) necesita tu serie de datos en nuestro caso tenemos 50 años de datos con lo cual no hay problema superamos los los superamos 10 los 10 años por ejemplo aquí le decimos que necesitamos 10 años de datos pero en este caso no es no es relevante este th y aquí le decimos tem Trend para que nos haga eh Para que nos calcule la tendencia temporal y nos produce eh un raster stack de tres capas eh Y una o sea de tres layers cada una
(1:31:04) de las layers si si abrimos este este objeto que es temporal Trend vemos que es un raster stack y que tenemos las tendencias temporales en la primera columna que sería spls trends el error estándar de nuestras tendencias temporales y la significancia de cada de de esas tendencias temporales por ejemplo vamos a verlo porque así nadie nadie lo entiende no es un poco difícil entender Este fichero raster stack pero
(1:31:33) si le decimos plote danoslo dibujan nolo aquí podemos ver que tenemos las tendencias temporales de temperatura que a nivel mundial en la mayoría del mundo están aumentando si por qué porque porque están en muchas zonas del planeta están a 0.06 que esto sería 0.06 gr al año Cuál es la la la la tasa de aumento de temperatura superficial en cada uno de nuestros puntos de nuestros de nuestros grids de malla de nuestros de nuestros eh datos Geo eh georreferenciados Y temporales y
(1:32:16) aquí le podemos nos daría Cuál es el error estándar de esta eh de esta tendencia temporal de cada uno de ellos y cuál es el p valor aquí se ha ploteado un poco eh se ha ploteado muy pequeño no nos da los datos de de error estándar pero si lo ploteos sin el r mar down o haciendo más grande esto no debería darnos problema o la otra opción es dibujarlo uno a uno para que nos Aquí vemos que nos llega al segundo decimal tenemos el podemos ver que hay algunas zonas
(1:32:51) donde las tendencias de temperatura son as pero en la gran mayoría dominan las tendencias positivas esto que la pendiente es es positiva está en aumento si le decimos tem tren de nuestro objeto en vez de uno le decimos dos nos dibujará nos dará la información del error estándar que lo tenemos aquí si lo vemos individualmente podemos ver hasta que Cuál es la resolución de nuestro error estándar y eh cuando le
(1:33:23) decimos entren tres obtendremos los datos de p valor esto quiere decir que en las zonas donde está blanco eh el p valor es cercano a cero Habría que ver habría que no entiendo por qué nos da no nos da Mayor Mayor información en este gráfico porque normalmente si entramos solo en tem tren 3 aquí podemos ver eh los p valores a los los p valores van desde 0.9 elevado a la 17 O sea que tendríamos información muy específica y muy detallada sobre el p valor de nuestra la significancia de nuestras regresiones
(1:34:08) eh todo esto se podría por ejemplo cortar y reducir para para para una región más pequeña por ejemplo la zona del Mediterráneo aquí en el Script he dejado Cómo cortarlo a nuestro mapa por ejemplo no nos interesa hacer esto mismo para una región más pequeña pues por ahí he dejado la función crop al final del documento he dejado la función crop que si tú tienes un polígono espacial si tú tienes las coordenadas del polígono que tú quieres crear y le dices que es un
(1:34:41) fichero un spaal polígono un polígono espacial con este con esta con esta parte de aquí en la que en la que pone cortar polígono haces un zoom haces una profundización en una zona específica usando este tipo de datos podemos llegar a calcular cosas como por ejemplo la velocidad del cambio climático no voy a entrar a a a profundizar en todo eso Pero si tenemos la tendencia temporal de de cambio de temperatura podemos calcular también cómo cambia espacialmente la temperatura
(1:35:15) Y utilizando el gradiente espacial y la tendencia temporal si dividimos la tendencia temporal por el gradiente espacial podemos calcular la velocidad del cambio climático no entraré en en estas funciones las he dejado aquí escritas pero eh yo creo que eh que no es necesario entrar a tanta a tanta profundidad por lo menos en el día de hoy no obstante si alguien le interesa Y quiere hacer funcionar esas esas funciones ver cómo funciona y tiene
(1:35:44) dudas me puede escribir por el por el el discord por la aplicación y ahí podemos resolver dudas eh No sabía muy bien Me he dejado un poco la parte del gradiente espacial y de la velocidad climática a la espera de ver cómo íbamos de tiempo qué tal vamos de tiempo yo creo que ya ya me he pasado no puede ser estamos bien Marina te comento todavía tenemos unos 19 minutitos para que puedas terminar y hacer una sesión de preguntas y respuestas vale muy bien quedan 20
(1:36:18) minutos verdad no hay nadie más que presente después de mí o sea que no hay problema sí correcto Cómo cómo cómo habéis funcionado hasta ahora os ha os ha funcionado todo o no contadme Si habéis logrado sacar los mapas de tendencia Temporal y si si habéis entendido qué es lo que estamos haciendo o no Porque eso es importante saber eh que me vayáis contando que nos contéis al resto si os da si os funciona o no os funciona Cómo vais
(1:37:07) todo ha funcionado bien vale entonces si nadie dice nada entiendo que a todo el mundo le funciona bien si nadie tiene dudas continúo pero si alguien tiene dudas No me importa parar ahora es un buen momento para parar y resolver las dudas que no hay problema no os doy un minutillo más para que os lo penséis pero pero si no hay problema Pues continúo a ver Mira Bernardo Leal eh claro Este es un tema que no
(1:37:40) encuentra la función sus Data me generó este problema efectivamente esto puede ser e esto no sé muy bien por qué puede darnos este error porque debería de funcionar a mí me funcionaba sin problemas estáis seguros que tenéis esta dirección bien puesta que no habéis borrado nada que no habéis cambiado nada a ver true no estoy viendo que está todo bien puesto esto le pasa a todo el mundo o no sé muy bien por qué puede por qué puede pasar esto Mirad si entramos aquí
(1:38:17) en este enlace de dónde están los datos lo que hacemos Es entrar directamente al del autor de este paquete de Jorge García molinos que él ha subido una serie de datos estos datos globales de temperatura y están disponibles no tenemos que entrar a descargarlos en en copernicus no sé muy bien ahora a mí no me está funcionando por ejemplo github No ahora mismo no sé por qué no está pero parece que lo que le está dando error es la función no encuentra No es que no encuentre el fichero sino la
(1:38:55) función ya sí no os iba a enseñar Cuál es el cuál es el cuál es el dónde nos lleva este fichero por lo menos para que supiera dónde debería llevarle el fichero Y por qué no y por lo menos saber qué es lo que qué es lo que debería conseguir pero no sé muy bien cómo cómo solucionar eso la verdad si alguien más puede echarnos una mano pero una opción es Buscar cuando algo no te das te está dando error pues una opción
(1:39:29) es Buscar a ver eh una pregunta Marina esa función source Data es algo que viene en con r o es parte de un paquete Eh Pues yo creo que es Yo creo que es es una pregunta basada en la ignorancia es r base tiene que estar en la instalación normal de r sí es de revés a meaba que no definitivamente n nadie debería tener ningún problema con eso no a m me sale que el source Data es de un paquete que se llama redm sí del paquete redm es el que me sale el source Data y los datos están en
(1:40:16) el repositorio de jthf sí os funciona este este paquete el de r Miss a ver a ver Bernardo a ver lo mismo lo mismo Bernardo tiene el problema con el paquete redmis a ver si corremos aquí justo si hacemos funcionar el el 30 la línea 38 veremos a ver si el problema es con con uno de los paquetes o no Marina Disculpa te puedo hacer una consulta problemas es posible realizar un mapa o algo ya representando las promedios de temperatura en base a las tendencias que estás generando porque son cambios de
(1:41:05) velocidad cambios de tendencia en los mapas que estabas presentando pero no es como que una proyección de las temperaturas como tal No sé si me dejo entender o sea por ejemplo puedes mostrar el gráfico que estabas mostrando hace un momento de temperatura Sí el de temperatura y el de la tendencia que estás a ver por ejemplo a ver el de la tendencia aquí Sí claro ese es un gráfico de tendencia pero por ejemplo de de del año actual al año 10 tengo entendido que pusiste 10 años verdad como el cambio de
(1:41:48) temperatura o no en este caso en este caso lo pusimos desde el año 1960 a 2009 aquí en la función Sam series es donde le indicamos de de cuándo a cuándo queremos que nos calcule que nos calcule exactamente eh los los años que nos interesan Entonces en este caso vamos de 1960 a 2009 y le indicamos eh que le que le estamos aportando los datos en mensuales para que nos proporcione el promedio anual entonces una vez tenemos el promedio anual ahí es cuando con la
(1:42:34) función temporal Trend nos calcula la tendencia temporal la recta de regresión linear en cada píxel de nuestro mapa no sé si es resuelto tu solo es para ese periodo pero no tiene nada que ver con ese 10 que pusiste en la función Trend más abajito no que te está calculando la tendencia de ese en base a ese modelo a 10 años no en este caso claro no en este caso en la función help le preguntamos vale le preguntamos tem Trend qué nos quiere decir ese th que no entendemos Pues aquí
(1:43:11) le decimos y nos dice mínimo número mínimo de observaciones que necesitas para calcular la tendencia en cada celda por ejemplo imagina que tú tuvieras un mapa donde no tienes años donde no tienes información todos los años en este caso tenemos mucha suerte porque estamos trabajando con tendencias con con temperatura superficial del mar de satélite pero si fuera por ejemplo estamos en el Ártico en el Ártico o en una zona de muchas nubes en el que no conseguimos tener
(1:43:44) datos para todas para todas las regiones Y de pronto Tenemos un mapa con con nubes que no han conseguido aportar la tendencia la temperatura superficial entonces lo que hacemos Es decirle pues mínimo necesito 10 años para que calcules necesito datos de los de los 50 años que me estás proporcionando mínimo tienes que tener 10 años para calcular la recta de regresión si no tienes 100 años no si no tienes 10 años no me lo calcules Y entonces tendríamos huecos vacíos en las
(1:44:18) celdas dur de las que no tenemos por lo menos 10 años Y en este caso tenemos la suerte de que son datos de 50 años entonces está está estupendo porque no todas nuestras celdas todos nuestros grid de de 100 km por 100 km Es decir de 1 grado por 1 grado contienen información de temperatura superficie Sí ya me qued Sid ya me quedó más claro digo entonces y en base a esa tendencia que tenemos ya para cada una de las grillas es posible hacer algún tipo de proyección o algo porque ya
(1:44:54) A qué te refieres con un tipo de proyección proyección a futuro ejemplo ya tengo Exacto Exacto eh Sí tú podrías en realidad lo que se suele hacer es eh descargarte o sea en este caso estamos calculando Cómo ha cambiado la temperatura superficial del mar en en un tiempo contemporáneo que se llama o o presente o sea desde el año que no es tanto presente en realidad es pasado pero bueno depende de dónde pongas el
(1:45:20) límite si es son Uy aquí yo puse No desde el año 1000 desde el año Sí desde desde el 1000 Esto está mal desde el 1960 al 2010 Pues sería Cómo ha cambiado en pasado si queremos si lo que queremos saber es cómo va a cambiar en los próximos 100 años lo que tenemos que hacer es descargarnos datos que ya están proyectados entonces entraríamos en buscaríamos eh los modelos cemip que se llaman y y en los modelos temip hay hay una serie de de ficheros de datos igual
(1:45:59) que estos de temperatura superficial del mar pero que ya han sido proyectados entonces tú puedes descargarte Cuál es la proyección Cuál es la previsión de de temperaturas superficiales del mar para esos años con sus correspondientes incertidumbres Por supuesto porque todo eso conlleva incertidumbre y podrías calcular Cuál es el la tendencia que se prevé y que se proyecta para los próximos 50 años o 100 años con con con este tipo de datos aquí no
(1:46:34) llegamos a tanto pero se puede hacer perfectamente con este tipo est esa sería una de las aplicaciones que tendría este este tipo de de Script y este tipo de análisis ya me quedó más claro Muchas gracias hay alguna pregunta más por ahí si no yo os sigo contando Cómo se hace el gradiente espacial y cómo se calcula la la velocidad Sí yo tengo una pregunta Sí mira en el caso de de cambio climático y de eventos el niño no
(1:47:12) e Cómo se podría hacer o o Med no un Script para hacer cómo te digo a ver tengo datos de 30 años continuos pero lo que quiero hacer es quitar los eventos de niños ahí que son digamos que lo voy a asumir que son anuales es decir supongamos que tengo desde el año dos 2000 bueno 1900 o 1950 a fecha actual Vale entonces ya en esto Quiero quitar los los años Donde ha habido algún evento climático de de de lenzo no Cómo podría ser para
(1:47:56) quitarlos pues eh yo creo que lo que se podría hacer es seleccionar las capas seleccionar hacer un seleccionar las capas de los años que que tú quieres eliminar y restárselo a tu raster stack eso podrías hacerlo puedes o crear uno paralelo o restarle esa serie de capas y quitarlas por ejemplo otra de las aplicaciones que que no he mostrado porque eh Por limitación de tiempo es seleccionar Qué meses quieres por ejemplo tú dices es que no quiero los meses de de primavera de de otoño y
(1:48:38) invierno Yo solo quiero saber el Cómo cambia la temperatura cada año durante el mes de verano porque es cuando suceden olas de calor por ejemplo o o porque son épocas eh meses especialmente relevantes o porque mi muestreo se hace en esos meses por ejemplo hay hay muchos motivos por los que puedes estar interesado en acotar la serie temporal puedes seleccionar los meses para los que quieres eh Para los que quieres que se calcule y
(1:49:08) del mismo modo puedes hacerlo Con los años que tú quieres que se que se calcule si alguien está interesado en eso yo puedo facilitarle un scp para para los meses y podría mirar también cómo hacerlo para cada año para el por ejemplo Vale pues pues yo sí estoy interesado tendrá noticias mías muy bien claro que sí gracias por ha puesto de ayuda nada bueno lo lo que os comentaba en el en el Script también podemos calcular el
(1:49:49) gradiente espacial que sería Cómo cambia la temperatura en función del espacio porque siempre hablamos de los cambios de eh de Cómo afecta el cambio climático a la temperatura a nivel temporal pero una de las grandes preguntas ahora mismo es Cómo afecta a nivel espacial y podríamos usar la función spaal gradient gradiente espacial para ver cómo son esos cambios a nivel espacial eh Y si si tenéis ganas de
(1:50:17) profundizar y dedicarle más tiempo a ese Script pues eh le dais Caña esta semana y y en algún momento a través de de la de la aplicación resolvemos resolvemos dudas No sé si alguien más tiene tiene preguntas en este caso es lo lo bueno creo que el tanto la clase de Héctor como la mía en no son exactamente compatibles porque utilizamos paquetes diferentes pero a fin de cuentas eh Es es muy importante poder poder entender cómo
(1:50:56) está codificada la la fecha en este tipo de ficheros y cómo se combinan porque en otras ocasiones no queremos trabajar con paquete raster y Terra sino que lo que queremos es extraer esos datos Y utilizarlos como un fichero csv por ejemplo como una matriz de datos que eso también se podría hacer perfectamente o sea podríamos convertir nuestro stack O Nuestro brick que tiene datos acumulados en una en una matriz de
(1:51:26) datos o en una en un Data frame en un en un fichero de de datos en tabla para trabajar con ellos de otra manera con otro tipo de sin la referencia o sea sin la referencia geográfica manteniendo la referencia geográfica como latitud y longitud en columnas individuales hay alguna pregunta más veo que no no sé si finalmente Bernardo lean siento que no hayamos podido resolver tu duda podemos podemos intentar Trata de hacerlo funcionar todo y y a lo mejor puede ser buena idea que
(1:52:13) os descarguéis los los scripts y los hagáis funcionar en vuestro propio en vuestro propio ordenador para ver Bueno a mí me parece Útil para para hacerlo funcionar Pero bueno cada uno trabaja como en Júpiter también es una buena manera en este servidor de aquí pero no sé muy bien Por qué a Bernardo no le ha funcionado Ah muy bien instalaste los paquetes y ya te funcionó dejó dejaste de tener problemas con el sus Data genial genial pues eso es lo que suele pasar muy a menudo una una de los
(1:52:51) errores frecuentes es que hay que volver revisar si han funcionado las librerías y comprobar que está todo bien instalado que no nos hemos dejado ninguna fila por ahí pendiente y volver a volver a hacerlo todo funcionar estupendo Bernardo me alegro alguien más tiene dudas o alguien más tiene alguna cosa que incluir os ha servido pensáis que os servir para algo o no Sí exactamente Es muy es muy aplicativo lo que estás mostrando lo que has enseñado Me parece
(1:53:39) genial tengo una pregunta ya que es ustedes manejan mucho información satelital Cuál es eh la entre la las plataformas que que hay de de libre acceso a los datos de temperatura Cuáles nos recomendarías Bueno me recomendarías a mí datos de temperatura cerca de la costa te digo cerca de la costa 5 km de de la orilla No qué datos de dónde podría sacarlo más reales No pues esto Julia mañ os explicar Cómo descargar los datos en
(1:54:23) copernicus por ejemplo copernicus sé que por ejemplo para la zona del Yo creo que es a nivel mundial copernicus que uno de los uno de los productos que tiene copernicus tienen datos de hasta 0.04 grados que estos son 4 km de de grid de malla O sea que es muy posible conseguir datos de una muy buena resolución de alta a nivel a nivel costero para toda la para toda la zona no no debería de suponer problema eh incluso de mayor
(1:55:02) resolución a 1 km Héctor Villalobos propone aquí por ejemplo nos ha mandado un enlace de NASA yo los utilizo de copernicus en Europa usamos mucho copernicus porque es el servidor que se usa pero muchísima gente trabaja con con NASA por ejemplo mira si entráis en el en el el enlace que nos ha mandado que nos ha mandado Héctor si queréis os lo os lo comparto sí Emilio también Comparte el produc otro de los productos e a ver mirar os comparto pantalla si
(1:55:39) queréis a [Música] ver dónde Bueno ahora no me sale Pero bueno no comparto pantalla pero si entráis en el enlace que que ha mandado Héctor podéis ver que NASA tiene tiene ficheros con con muy buena resolución y la gran mayoría de Ellos están en formato eh el formato os lo dejo por el chat escrito net cdf net cdf 4 y en este caso en el Script os lo puedo enseñar ya ya no creo que ya deberíamos acabar pero
(1:56:17) bueno os enseño a ver a ver dónde tenía yo el Júpiter será que lo he cerrado un momento HM a ver aquí tengo el servidor Aquí está aquí estáis viendo la pantalla verdad de r estudio Sí aquí en una de las partes del Script lo que yo he dejado si si lo buscáis en opción B la una de las eh opción a sería utilizando los datos que os que os que ya os he facilitado que está en que están en github pero si usamos la opción b
(1:57:06) m Aquí hay una manera de a ver dónde lo puse estará arriba aquí en opción B hay una manera en la que tú puedes abrir los datos directamente eh tú te los descargas de copernicus que los que es lo que nos o de Nasa que lo que nos va a enseñar Julia mañana los guardas en una en tu carpeta aquí le indicas le dices set working directory dónde tienes esos datos de temperatura eh lo abres como fichero en sí y
(1:57:35) directamente le dices brick función brick y tu tu netcdf Entonces en este caso ya rápidamente te te puedes trabajar con ello y ver los mapas directamente es una manera muy rápida de de verlo visualmente que que muchas veces nos ayuda a comprender Qué tipo de datos tenemos a nivel a nivel visual es una manera de eso abrirlo rápidamente Sí hay alguna pregunta más Bueno si no surgen más preguntas os quería comentar que bueno he puesto por el discord no sé si lo habéis visto o no
(1:58:31) pero el Jorge Cornejo que es el tercer ponente de hoy de nr Eh no ha podido asistir hoy en vivo Pero ha hecho una grabación de su presentación de su tutorial entonces Bueno al principio habíamos pensado ponerlo a continuación pero la verdad es que quizás es mejor que lo colguemos directamente en la plataforma de decir y en el y en la web de interconecta y que lo veáis a vuestro ritmo podéis ir parando y y haciendo así siguiendo el Script también en la
(1:59:03) plataforma y bueno si os parece bien a todos Quizás es una mejor forma de actuar en este caso ya que Jorge no está aquí y os puede responder a cualquier pregunta por por discord no sé si qué os parece Así que teníamos do horas y media planeadas para hoy Ho Pero quizás lo podemos dejar aquí con estas dos horas y y como digo que vosotros vayáis trabajando más a vuestro ritmo eh esta tercera este tercer tutorial que es sobre sobre informes en
(1:59:37) r hacemos así si te parece Alejandra Sí De acuerdo marián muchas gracias No sé si alguien tiene alguna otra pregunta o algún otro comentario también relacionado con la plataforma interconecta aprovechando y y bueno eh mientras ustedes se deciden si si tienen alguna pregunta o demás solo recordarles por favor que completen la lista de asistencia que estuve compartiendo en el en el chat la voy a volver a enviar es muy importante para nosotros tener la
(2:00:06) lista de asistencia no solo para registrar que ustedes han asistido sino Para futuras convocatorias e informaciones que compartimos vía correo electrónico Así que agradezco que se tomen unos minutitos para que ahora que terminemos puedan llenar esa lista de asistencia Bueno un que yo voy a bombardear con preguntas Lo siento Eh bueno tenía una pregunta para para vosotros era sabéis si hay algún lugar bueno un sitio donde podréis
(2:00:41) Descargar eh datos no así de oceanografía Pero dentro de bahía para dentro de las bahías en el el caso de de latinoamérica la zona de Perú de bahía Independencia hay datos fuera de la Bahía la que da cara al mar pero no dentro sabéis si Existe algún algún lugar donde Descargar datos si hay este imágenes de satélite o no sé no sé si me explicado bien yo creo que eso depende mucho de la resolución espacial por ejemplo en algunos de los datos de copernicus que están a 9.
(2:01:23) 2 km pues no vas a tener datos en bahías pequeñas no tendrías que irte a los de 4 km pero esos están Digamos como variables independientes no no en esos productos que agregan por ejemplo copernicus hace mucho eso no que dice a ver variables físicas variables biológicas en fin yo sugeriría que aras por ejemplo el sitio de NASA Ocean color lo puedes poner así en el en el Google NASA Ocean color que ahí están digamos las variables separadas encuentras
(2:02:04) clorofil encuentras temperatura y verás que en muchas de esas tienes al menos hasta Bueno más bien lo más detallado es 4 km entonces Esas probablemente información al menos de algunos píxeles dentro de las bahías y si tú contrast eso con los que mencioné hace rato bueno el link que puse para los de Ultra alta resolución de 1 km verás que ese sí te da en las bahías sin problemas no eh Hay muchos datos ahí vale Vale gracias si
(2:02:36) puedo Añadir también eh Como dice Héctor para esto Generalmente eh No lo vas a encontrar en productos agregados eh Y Tenés que ir más acerca de la Fuente de los datos eh Héctor mencionó producto de 1 km el el satélite o el sensor de la nasta o los Estados Unidos por lo menos Bears tiene resolución hasta 250 m eh pero probablemente para acceder los tendrás que ir directamente a la fente y ver que Qué tipo de producto ya existe o sea que es un poco más difícil eh
(2:03:12) eh encontrar Dónde están y investigar eh Qué tipo de producto satisfacen tus necesidades me imagino que no estoy tan tan familiarizado con el los satélites sentinel de de Europa Pero me imagino que hay uno también de alta resolución pero igual supondría que no no hay un producto tan fácil de acceder no fácil tan fácil de encontrar tal vez en la pregun la manera de describirlo Eh entonces tenés que eh hacerte todo un proyecto de investigación ir directamente a la fuente y y de ahí si
(2:03:47) queremos aún extender más el el eh lo que estás buscando depende de la naturaleza de tu pregunta Supongo por ejemplo hay datos tanto de lanza como de sentinel eh de resolución muy muy alta de 30 Met o algo así pero el lo que perdés es la resolución temporal Pero eso depende de la pregunta que estás haciendo Si en tu caso tener una resolución super alta e pero una vez cada una semana si eso satisface tus necesidades entonces perfecto vale Ya te
(2:04:28) entiendo Bueno lo que quería saber era si existe datos históricos temporal históricos No desde 1950 actual Si estás buscando saté no existen datos de saté antes de 75 8 muy útiles entonces esa esa es una respuesta fácil Vale y de los últimos 20 años quizás no Definitivamente sí pero de nuevo depende de que tus requisitos eh si son de Así que si ests extremadamente ambicioso y queres una resolución super alta eh
(2:05:08) espacial y una resolución temporal super alta vas a tener que investigarlo bien vale pues decir ya tengo trabajito Mira por ejemplo en este caso yo que trabajo en el Mediterráneo de del de desde copernicus aquí por ejemplo sí que os puedo enseñar en este enlace que he pasado por el chat eh Cómo entrar en este enlace de copernicus y te dice mediterránea en sí alta resolución eh temperatura superficial del mar reprocesado entonces aquí tenemos datos del Mediterráneo y de
(2:05:44) la costa atlántica hasta hasta y toda la costa del continente africano y ahí estos datos sí que son de de 0.05 que son 5 km de 5 km por 5 km de muy alta no existen más de mayor alta resolución en el Mediterráneo Yo diría que más de que mayor de esta resolución espacial Yo diría que no existe puede que me equivoque pero diría que no de fuentes de de satélites y de modelos reprocesadas zonas como indica Emilio hay que buscarlo como un exacto de esa
(2:06:22) región es probable que sí que exista eh pero yo ahora mismo no se me ocurre Héctor nos ha mandado un enlace que yo no lo no lo conozco es este enlace que ha mandado Héctor ha sido este verdad sí que en el Ocean virtual laboratory ahí se podría se podría buscar también también más cosas Pero efectivamente hay que hacer búsquedas muy concretas de de esa zona vale en la la ventana anterior que me mostraste del
(2:06:56) Mediterráneo hay datos antes o sea se tomarlos en cuenta de la zona Delta elbro s sin problema El delta delro Mira El delta debro está aquí esto Esto es la zona del delta del Y entonces esto quiere decir que hay datosa es bastante buena es es bastante buena y de la Bahía el alfax que está cuál que va ahí del alfax está atrás del tele de la Es que yo trabajo ahí est ha estoy haciendo pesquerías ahí Ah de al directamente sería difícil para cualquiera de nosotros
(2:07:48) darte una respuesta tan concreta estás preguntando lo que te recomendaría es que vayas a ese sistema y lo no es que er Aprovechando que tenemos un poco de tiempo era un poquito investig bueno no si me voy a poner a investigar Gracias Marina cómo era al alfas del pi teer no al fax El Alfa laa del Alfa que está del celebro dond está el puntero que está señalando más abajo no este es estás yote más hacia arriba el hebro Ahí está río río Delta El hebro
(2:08:26) ahí sí como dice ex si te acercas si te acercas más es verdad que luego conforme vas aumentando la resolución a nivel costero perdemos hay menos datos Sí al final las cosas globales yo creo que a nivel a nivel de investigación están bastante abarcadas y conforme vamos reduciendo la la localización y queremos cosas más locales cada vez todo se hace más difícil para casi todo muy bien pues muchas gracias algunos compañeros tendrán
(2:09:15) preguntas Sí mira Abigail nos dice en el sentinel 3 genera datos a 1 kilómetro y ha mandado un enlace de de NASA efectivamente ahí se pueden se pueden buscar es muy probable que haya que haya genial pues muchas gracias important información eh Adrián tiene la mano levantada buenas se logra logra buenas les pido disculpas que to entr saliendo tenido problemas de internet en la universidad si quiereser una pregunta es
(2:10:23) para poder descargar datos del Grace el Grace foo creo que también de altimetría y y de radar insar etcétera Cómo se se descargaría y hay librerías en python que permite el procesamiento de de esa Data eh Cuál sería el repositorio esa Data temporal eh de satélites como el Grace y similares ya no me digas que estoy ah no s est en se me logró Escuchar Sí sí excelente voy a comentar algo pues ya que nadie está comentando eh En primer lugar dijiste Grace y en la
(2:11:19) misma oración radares entonces lo que diría es que son completamente diferentes no tienen absolutamente nada que ver entonces no van a encontrar un sitio que tendría los dos de una manera pues posible pero son tecnologías super diferente la y no hay otra cosa como Grace sea no hay así como para temperatura hay una gama de sensores y bueno para los demás es un satélite tema especial eh que mido los cambios muy pequeños en gravedad a la medida que va pasando y basado en esos
(2:11:59) cambios de gravedad e e estima puede estimar cosas como el cambio en en la masa de agua el efecto de esa masa de agua enc suelos en el mar en en en la gravedad que detecta satélite la resolución me imagino que Grace ya sabes que la resolución es bastante cruda la de Grace eh Ahorita no recuerdo Perdón que sí ese le sigo estoy confundido un poquito con con los nombres de los de los satélites eh sé que hay algunos que están enfocados en en gravimetría
(2:12:42) eh pero está buscando también algunos que estuvieran con enfoques eh altimétricos eh precisamente para relacionar con información del nivel del mar por ejemplo y amarrarlo a otros a otros parámetros entonces La pregunta era repositorio donde uno pueda conseguir esos datos s y librerías de python o de r eh que permitan procesar digamos por ejemplo cuando es insar el procesamiento se suele trabajar Que si en en envi o alguno de estos programas
(2:13:18) con licencia entonces a ver si existía el equ equivalente que se pudiera utilizar en python o en r eh Y pues eso los repositorios de datos temporales de altimetría y similares voy a dice mar que Laura pod hablar de eso mañana fantástico lo voy a decir que de dejemos Grace a un lado altimetría ya Esto me lo pone en una categoría más clara eh Y sí ha habido una gama de sensores Que hacen altimetría Eh y existen productos ahorita no me acuerdo te recomiendo que nos hagas la
(2:13:54) pregunta a través del Canal general en discord porque así podemos seguir con eh con más tiempo e sé que existen librerías para esto pero no se me vienen ahorita eh algo que puedo mencionar y algo sumamente emocionante es que desde hace un año eh año y medio la NASA creo que en colaboración con e con Francia los Estados Unidos y la nas pues y la agencia espacial francesa tienen una colaboración que yo creo que Data de casi 30 años de la tecnología de punta
(2:14:28) en altimetría Marina eh Y este nuevo satélite que se llama swot swot es revolucionario en términos de la resolución que tiene e no sé Exacto Gracias eh No sé el estado actual del acceso Pero me imagino que sé que ya está operacional y me imagino que los datos ya se pueden acceder de manera fácil eh pero ahorita no se me viene de nuevo si por favor hacer la pregunta en discord y ahí le daremos seguimiento y luego mañana tal vez con laur podemos
(2:15:00) podes hacer la pregunta de nuevo y podemos hablar un poco de eso porque estoy seguro que Swat es algo que sería de interés los productos de altimetría de SW e dada la alta resolución tanto espacial como temporal creo serían de gran interés lo que ya dado que este es un sensor nuevo lo que no te va a dar por supuesto es historia no te va a dar nada que ocurrió antes de B supuesto pero para eso estoy seguro
(2:15:23) que hay producto Sencillamente no se me vienen ahita perfecto muchas gracias Bueno adelantó voy a pedir disculpas creo que mañana no voy a poder asistir pero justamente a las 10:30 terminó una defensa de tesis de chicos de pregrado a nivel de licenciatura De hecho mi tesista es la que va va a defender entonces pido disculpas De antemano pero el resto del equipo que también estamos inscritos probablemente estén entonces cualquier cosita cuadro con ellos de Cuál es la
(2:15:57) las preguntas de seguimiento Muchas gracias de cualquier modo realmente te recomiendo mand las preguntas Sí por supuesto por supuesto Gracias nada Bueno veo que no hay más preguntas e Sí veo que Alejandra ya iba hablar también damos por finalizada la sesión de hoy Y de nuevo animaros a ver revisar todo el material entrar en la plataforma revisar los vídeos que os iremos colgando y ya continuamos mañana con Alejandra si quieres decir alguna cosa más no solamente María Muchas
(2:17:04) gracias Bueno agradecerte también a ti Marina a Héctor a Emilio y a todo el grupo de expertos que están detrás no solo en este momento en las sesiones síncronas sino permanentemente respondiendo todos los planteamientos dudas y demás y bueno como decía marián nos encontramos mañana a la misma hora yo les Volveré el enlace para que lo tengan todos a mano Y pues muchas gracias por su asistencia que tengan feliz
(2:17:28) tarde Gracias gracias a todos esta mañana buen día y buenas noches allá en Europa
Usamos cookies para asegurar que te damos la mejor experiencia en nuestra web. Si continúas usando este sitio, asumiremos que estás de acuerdo con ello.
