descarga de información espacial desde R

[Jesicapaz89]

@lidefi87
Hola! estoy intentando acceder a la información de temperatura superficial del mar desde el servidor en R-Studio, por ejemplo quisiera poder acceder a la info del datasets id nasa_jpl_d686_3623_a91b // https://upwell.pfeg.noaa.gov/erddap/search/advanced.html?page=1&itemsPerPage=1000&searchFor=&protocol=griddap&cdm_data_type=%28ANY%29&institution=%28ANY%29&ioos_category=%28ANY%29&keywords=%28ANY%29&long_name=%28ANY%29&standard_name=%28ANY%29&variableName=sst&maxLat=-28&minLon=-69&maxLon=-51&minLat=-43&minTime=&maxTime=)

sería MODIS AQUA L3 SST THERMAL MONTHLY 9KM NIGHTTIME V2019.0 (MODISA L3 SMI, MODIS AQUA L3 SST THERMAL MONTHLY 9KM NIGHTTIME V2019.0) [time][lat][lon], 0.08334°, 2002-present

copio y pego el script de lo que hice hasta el momento, aunque lo intento no puedo lograr que me muestre esta información (la cual sí está en la página)

noaa_url <- servidores %>%
filter(short_name == "UAF") %>%
#Ahora que sabemos el nombre corto, podemos utilizar esto para filtrar los servidores
pull(url)#extrae la información que está en una celda

#Veamos el resultado
noaa_url #ya tengo el servidor

?ed_search
datos=ed_datasets('grid')

#Podemos hacer una busqueda sencilla de datos de temperatura de la superficie del mar (SST por sus siglas en ingles) usando la funcion ed_search.

temp_sup_res_ <- ed_search(query = "sst",
which = "griddap",
url = noaa_url)
head(temp_sup_res_$info)

temp_sup_res_2 <- ed_search_adv(query = "sst",
#Datos en grilla
protocol = "griddap",
#Servidor erddap
url = noaa_url,
#Datos mensuales
keywords = "monthly",
#Limites espaciales
#Limites temporales
minTime = "2010",
maxTime = "2020")

también quisiera saber si en la carpeta donde se descarga se puede acceder a cada archivo particular, por ejemplo vimos este ejemplo dado que solo habia un archivo en esa carpeta : archivo <- list.files("Data/", pattern = ".nc", full.names = T)

[lidefi87]

Hola Jessica. No me queda muy claro cuál es tu pregunta. ¿Tienes un vínculo a tu script? ¿Qué error te sale y en qué sección del código? Yo acabo de buscar el set de datos que te interesa (nasa_jpl_d686_3623_a91b) y me apareció como una opción en la búsqueda avanzada. ¿Me podrías dar más detalles de lo que necesitas?

En cuanto a tu segunda pregunta, list.files() es una función que te da una lista de todos los archivos que se encuentran dentro de una carpeta en particular. Si utilizas el argumento pattern, entonces puedes filtrar archivos que contengan palabras específicas. En el ejemplo que dí, enlista todos los archivos dentro de la carpeta Data que contengan .nc en su nombre.

[Jesicapaz89]

Todavía no me queda bien en claro como editar desde el jupiter lab por eso lo hago desde mi computadora. ami ese código cuando lo busco en la línea 42 del script que subí acá no me figura. es decir sí figura en la página ERDDAP pero no cuando lo hago desde en R, de hecho no me aparece la información de ningún satélite MODIS

Jesicapaz89/proy5-regiones-comparacion@517757a

[lidefi87]

Acabo de resolver tus dudas, mira este script en el repositorio de tu proyecto: https://github.com/Intercoonecta/proy5-regiones-comparacion/blob/main/rdatosespacialesjesi_zonagaviotines.R Busca mis iniciales DFA y te explico lo que estoy haciendo

[Jesicapaz89]

Genia total! era eso, muchas gracias!