ACTIVIDADES PRELIMINARES

Crea dentro de ese directorio de trabajo una carpeta denominada ejercicio1.

if(!dir.exists("./ejercicio3"))                                                 # Si no existe la carpeta `ejercicio1` dentro del directorio de trabajo
  dir.create("./ejercicio3")                                                    # Crea la carpeta `ejercicio1`

Descarga directamente el fichero comprimido con los datos de la actividad.

download.file("https://personales.unican.es/rasillad/docencia/g174/6_manejo_raster/datos/actividad_evaluacion_continua_3.zip",    # URL fichero a descargar.
              "./ejercicio3/actividad_evaluación_continua_3.zip",         # Nombre que recibirá el fichero y su localización.
              method = "auto")                                                  # Método usado por defecto para descargar los ficheros.

Descomprime el fichero *.rar*

unzip("./ejercicio3/Datos_actividad_evaluación_continua_3.zip",                 # Fichero que se debe descomprimir
      exdir = "./ejercicio3")                                                   # Directorio/carpeta donde se ubicarán los ficheros descompromidos

Convierte esa carpeta en tu directorio de trabajo

#setwd("..../ejercicio3")

Por último, activa el paquete raster

EJERCICIOS DE EVALUACIÓN

PREGUNTA 1: Con los ficheros de imagen descargados y descomprimidos, lista los ficheros que existen en la carpeta dónde has alojado los ficheros de la actividad

PREGUNTA 2. Crea un objeto RasterBrick con nombre brick que contendrá los ficheros correspondientes a las bandas 1, 2, 3 y 4.

Para ello deberás crear primero un objeto RasterStack y luego convertirlo en RasterBrick.

PREGUNTA 3: cambia los nombres de las bandas del objeto brick, de tal manera que cada una de ellas siga la sintaxis “B1”, “B2” etc…

PREGUNTA 4: Graba el objeto brick con el nombre bandas1234 con formato *.tif. Tras ello, elimina todos los objetos existentes en el Global Environment.

PREGUNTA 5: Importa el fichero bandas1234.tif y denomínalo raster1.

PREGUNTA 6: Determina los siguientes atributos del objeto raster1:

  • Número de celdas

  • Número de capas

  • Resolución espacial (dimensión x).

  • Sistema de Coordenadas de Referencia (CRS).

PREGUNTA 7: Comprueba qué nombre tienen las bandas 1, 2, 3, y 4 de Landsat 8 y cambia en consecuencia las 4 bandas del objeto raster1.

PREGUNTA 8. Extrae la capa correspondiente a la banda “azul” y denominala banda_azul.

PREGUNTA 9. ¿Qué tipo de objeto raster se ha creado al importer ese fichero?.

PREGUNTA 10. El objeto raster asocia internamente una serie de atributos y propiedades que en R se denominan slots. Lista los slots que componen el objeto que acabas de importer.

PREGUNTA 11. Comprueba el tamaño del objeto creado en función del número de filas, columnas y de pixels del objeto que acabas de importar.

PREGUNTA 12. Extrae la resolución especial del objeto y el Sistema de Referencia de Coordenadas del objeto que acabas de importar.

PREGUNTA 13. Calcula los siguientes valores estadísticos del objeto que acabas de importar

PREGUNTA 14. ¿Cuál es el valor del pixel correspondiente a?

  • Fila 43

  • Columna 68

  • Fila 87 y columna 34

PREGUNTA 15. ¿Cuáles son las coordenadas de los siguientes píxeles?

  • Coordenada y del píxel 232

  • Coordenada x del píxel 444

  • Coordenadas x e y del píxel 623

PREGUNTA 16. Crea un objeto RasterStack con nombre stack1 que contendrá los ficheros correspondientes a las bandas 5, 6 y 7

Creación del stack

PREGUNTA 17: Comprueba qué nombre tienen las bandas 5, 6 y 7 de Landsat 8 y cambia en consecuencia las 3 bandas del objeto stack.

PREGUNTA 18: Graba los datos del stack usando el nombre bandas567 y el formato *.grd.

PREGUNTA 19: importa el fichero *.grd que has salvado anteriormente con el nombre raster2.

PREGUNTA 20: Contesta a las siguiente preguntas acerca del objeto raster2.

  • Número de filas.

  • Número de columnas.

  • Dimensiones.

PREGUNTA 21. Identifica cuál es la segunda capa del objeto raster2.

PREGUNTA 22. Cambia los nombres de las capas a su denominación original: azul, verde, rojo, nir.

PREGUNTA 23. Comprueba si ambas bandas tienen igual extensión, número de filas, columnas etc…

PREGUNTA 23: Visualiza de forma conjunta todas capas que conforman el objeto raster2.