ndvi: Calculo del NDVI para cada seccion censal
In jsegu/greennessr: Espacios verdes en ciudades

Description Usage Arguments Details Value Examples

View source: R/calcular_ndvi.R

Extrae el valor del índice de vegetación normalizado (NDVI) para cada sección censal. Para ello utiliza los datos de una imagen del satelite Landsat, a la que se le realiza el cálculo del NDVI y una capa de polígonos correspondiente a las secciones censales. Por defecto el método de extracción es la media de píxeles de cada polígono, este método se puede modificar con el parámetro calculo.

ndvi(
  dir_img,
  sc,
  calculo = mean,
  dev_raster = FALSE,
  landsat8 = FALSE,
  puntos = FALSE,
  radio = 300
)

`dir_img`	Carácter: Directorio donde se encuentra la imagen Landsat, tiene que estar descomprimida.
`sc`	Objeto espacial `(SpatialPolygonsDataFrame / sf):` con las geometrías poligonales de las secciones censales.
`calculo`	Función para resumir los valores NDVI en la capa de polígonos, por defecto `mean`. Se puede modificar (ej. median, min o max).
`dev_raster`	Lógico: Si es TRUE en el resultado se incluye la capa ráster con los valores del NDVI.
`landsat8`	Lógico: Por defecto se calculan de manera automática las imagenes de landsat-5 y landsat-7. En caso de querer calcular una imagen landsat-8, cambiar el parámetro a `TRUE`.
`puntos`	Boleano `TRUE o FALSE:` Por defecto es `FALSE`, poner a `TRUE` si la base de datos `sc` es de puntos, dónde la propia función creará un buffer para convertir el dato puntual a polígono.
`radio`	Numérico: Se utiliza si el parámetro `puntos = TRUE`. Distancia en metros del radio del buffer, por defecto 300m.

El índice de vegetación normalizado (NDVI) se utiliza para estimar la cantidad, calidad y el desarrollo de la vegetación. Este estimador se calcula mediante los valores de intensidad de la radiación que las plantas emiten o reflejan en ciertos rangos del espectro electromagnético. Concretamente en las longitudes de onda del rojo y del infrarrojo cercano (nir). El cálculo es el siguiente:

NDVI = (nir - rojo) / (nir + rojo)

Para que la función pueda calcularlo es necesario una imagen satélite multiespectral, en este caso se utilizan las de los satélites Landsat-5, Landsat-7 y Landsat-8 (en caso de utilizar Landsat-8 cambiar el parámetro landsat8 = TRUE). Para descargar las imágenes es recomendable usar la página web earth explorer de la USGS. Es muy importante que al realizar la descarga se tenga en cuenta lo siguiente:

Que las imágenes seleccionadas sean de la colección 'Landsat Collection 1 Level-1 (L1TP)' ya que están radiometricamente calibradas y ortorectificadas.
Que la zona de estudio no este cubierta por nubes, en la página web de earth explorer podrá visualizar la imagen antes de descargarla.
Que la imagen cubra toda la superficie de la zona de estudio. En el caso de tener que descargar más de una imagen para cubrir su zona de estudio tendrá que juntar las imágenes previamente antes de introducirlas en la función. Esto lo puede realizar usando la función merge.

En este link encontrará más detalles a cerca de la descarga y los distintos niveles de procesamiento de las imágenes.

No hace falta tener en cuenta los sistemas de referencia de los datos geográficos, la propia función los cambiará si fuese necesario.

El resultado que devuelve es el mismo objeto espacial (SpatialPolygonsDataFrame / sf) introducido en el parámetro sc de la función, con una columna nueva llamada ndvi que tiene el resultado para cada polígono.

Si también se quiere obtener el mapa ráster con el resultado del ndvi, se tiene que cambiar el parámetro dev_raster = TRUE. Entonces el resultado devuelto es una lista con dos elementos: El primero es el objeto espacial (SpatialPolygonsDataFrame / sf) comentado anteriormente y el segundo una RasterLayer con el resultado del ndvi en forma de mapa de píxeles.

## Not run: 
library(greennessr)
sc_valencia <- ndvi(
  dir_img    = "~01_data/LE07_L1TP_199033_20000621_20170211_01_T1",
  sc         = sc_valencia,
  calculo    = median,
  dev_raster = FALSE,
  landsat8   = FALSE,
  puntos     = FALSE,
  radio      = 300
)
sc_valencia
## End(Not run)