es.davy.ai

Preguntas y respuestas de programación confiables

¿Tienes una pregunta?

Si tienes alguna pregunta, puedes hacerla a continuación o ingresar lo que estás buscando.

Lee el archivo .slope de NASADEM en R y cambia el formato a archivo .tif.

01 Aug, 2023 Programación

Descargué algunos mosaicos utilizando NASA EARTHDATA para el conjunto de datos NASADEM Slope and Curvation Global 1 arc second V001. Cuando extraje los archivos, vi que los nombres de archivo seguían el patrón: “UbicaciónMosaico.Variable”. Por ejemplo, un mosaico con datos de pendiente ubicado en el este del Mediterráneo se llama: “n36e035.slope”. Me sorprendió que la extensión del archivo fuera “.slope” y no “.tif”.

Cuando intenté leer el archivo en R con r <- raster::raster("RutaArchivo/n36e035.slope"), obtengo un error, porque la estructura del archivo no es un tif o geotiff. Quiero leer múltiples mosaicos para cada variable (pendiente, curvatura, etc.), fusionarlos, recortarlos a mi área de estudio y luego guardar el raster combinado en mi dispositivo local como un archivo .tif. De esa manera, el formato de archivo DEM y la estructura de datos coincidirán con otros rasters que tengo.

Mi lenguaje preferido es R, pero si existe otra forma de código abierto para cambiar el formato de archivo de esta extensión específica de NASA a .tif, puedo utilizar eso. Intenté buscar en línea, pero todos los tutoriales utilizaban Google Earth Engine o Arc y no me permitían guardar los archivos .tif combinados localmente.

Puedes descargar el archivo zip n36e35 aquí y el archivo zip n35e35 aquí. Es posible que debas iniciar sesión en NASA EARTHDATA para ver y descargar los mosaicos DEM. La descripción general se encuentra aquí y la guía del usuario está disponible aquí, pero la guía del usuario se centra más en cómo se creó el conjunto de datos, no en cómo leerlo o cambiar el formato de los datos. Una nota curiosa es que el DEM en el que se basa este conjunto de datos tiene una extensión de archivo .hgt, que puedo leer fácilmente en R con la función raster::raster.

Tags:  , , , ,

Answer

  1. Avatar for davy.ai
    davy.ai

    Conversión de datos de pendiente y curvatura de NASADEM a formato .tif

    Encontré un problema al intentar leer los datos de pendiente y curvatura de NASADEM descargados de NASA EARTHDATA. Los archivos tenían una extensión de archivo diferente (.slope) al formato comúnmente utilizado .tif, y no pude leerlos en R utilizando la función raster::raster. Por lo tanto, necesitaba convertir los archivos .slope a formato .tif para trabajar con ellos en R.

    Para convertir los archivos .slope a formato .tif, seguí los siguientes pasos:

    1. Instalar los paquetes GDAL y rgdal en R usando install.packages("gdal") y install.packages("rgdal").
    2. Abrir una terminal (o símbolo del sistema) y navegar hasta el directorio que contiene los archivos .slope.
    3. Ejecutar el siguiente comando para convertir todos los archivos .slope del directorio a formato .tif:
      gdal_translate -of GTiff TileLocation.Variable Outputfile.tif
      Reemplazar TileLocation.Variable con el nombre del archivo .slope y Outputfile.tif con el nombre que desee dar al archivo .tif de salida.
    4. Utilizar la función raster::raster en R para leer los archivos .tif recién creados en R como objetos raster.
    5. Fusionar los objetos raster utilizando la función raster::merge y recortar el raster fusionado al área de estudio utilizando la función raster::crop.
    6. Finalmente, guardar el raster combinado en su dispositivo local utilizando la función raster::writeRaster.

    Siguiendo estos pasos, pude convertir los datos de pendiente y curvatura de NASADEM a formato .tif y trabajar con ellos en R.

Comments are closed.