Ejercicio 1

G174: Fotointerpretación y Teledetección

Activación del paquete raster

1 Apartado 1: generación de un rasterBrick.

Crea un objeto raster brick cuyo nombre será brick1. Este objeto consistirá en 4 bandas, denominadas B1, B2, B3, B4 y B5. Cada banda consistirá en 100 píxeles, cuyas características serán las siguientes:

• B1: números (enteros) aleatorios sin reemplazo entre 17 y 89.

• B2: números (decimales) entre 3 y 87.

• B3: números (enteros) aleatorios sin reemplazo entre 44 y 365.

• B4: números (enteros) aleatorios con reemplazo entre 83 y 148.

• B5: números (decimales) entre 1 y 63.

PASO 1: creación de vectores con esas condiciones

PASO 2. Construcción de un array a partir de esos vectores

PASO 3. Generación del rasterBrick e impresión en pantalla.

2 Apartado 2: cambio de nombre a las capas e indexado del rasterBrick.

¿Cuál es el nombre de la capa 5 y el valor correspondiente al pixel situado en las columnas 3 y 6?

¿Cuáles son los valores correspondientes a la banda 1 y a sus filas 7 a 9 y a las columnas 2 a 5.

3 Apartado 4: generación de un rasterLayer (ráster monobanda) a partir de un rasterBrick

A partir del objeto brick1 se extraerá la capa B4, sus valores se dividirán entre 0.73 y se formará con ella un objeto raster monobanda denominado raster1.

4 Apartado 5. Análisis estadístico de los valores del objeto ráster monobanda.

• ¿Hay valores repetidos en raster1?

• ¿Calcula los siguientes estadísticos correspondientes a raster1?:

• Promedio de todos los píxeles

• Suma de todos los píxeles

• Frecuencia de valores

• ¿Cuál es el valor del pixel 50? ¿y del pixel 85?

• ¿Cuáles son los valores de los píxeles que van del 17 al 25?

• ¿Cuáles son los valores de las siguientes filas y columnas

• Fila 5

• Columna 2

• Columnas 5 a 8

• ¿Cuál es el valor del píxel localizado en en la 5ª fila y la 3ª columna

• ¿Cuál es el valor de los píxeles localizados entre las filas 6 y 8 y las columnas 3 y 4

• ¿Cuáles son los valores que corresponden a los vértices de ráster1?

• ¿Cuál es la coordenada y de la celda 67? ¿y la coordenada X de la celda 99?

¿Y las coordenadas de las celdas 64,65,66 y 67?

• Y recíprocamente, a partir de coordenadas x e yse obtiene el número de celda.

• ¿A qué columna corresponde el pixel 83?

• ¿A qué fila corresponen los píxeles 21, 46 y 103?

• ¿Cuáles son las filas y columnas correspondientes a los píxeles 11, 63 y 89?

• ¿Cuál es el pixel correspondiente a la fila 9 y a la columna 3?

5 Apartado 5. Modificación manual de los valores de un objeto ráster.

En este apartado se van a modificar las bandas 1 y 5 del rasterBrick con valores introducidos manualmente. Para confirmar que se ha realizado correctamente, esa banda deberá ser representada gráficamente. Los cambios serán los siguientes:

• Banda 1: se sustituirá el valor original del pixel correspondiente a la fila 1 y la columna 6 por valores ausentes (\(NA\))

• Banda 5: se sustituirá el valor original del pixel de toda la fila 4 y de toda la columna 10 por (\(-9999\))