Intrduction

This is an R Markdown exercise on how to visualize geospatial data in R

Load the necessary packages

library(raster)
## Loading required package: sp
library(ggplot2)

Call tiff files for mapping

ndvi_2000 <- "C:/Tamirat_2024/IGSSA_PhD/Y1_S2/1_Geospatial_Modeling/Part_2_Geostatistical_Modelling/Project/TIFF_files/N1_SER_NDVI_2000_2.tif"
ndvi_2005 <- "C:/Tamirat_2024/IGSSA_PhD/Y1_S2/1_Geospatial_Modeling/Part_2_Geostatistical_Modelling/Project/TIFF_files/N1_SER_NDVI_2005_2.tif"
ndvi_2010 <- "C:/Tamirat_2024/IGSSA_PhD/Y1_S2/1_Geospatial_Modeling/Part_2_Geostatistical_Modelling/Project/TIFF_files/N1_SER_NDVI_2010_2.tif"
ndvi_2015 <- "C:/Tamirat_2024/IGSSA_PhD/Y1_S2/1_Geospatial_Modeling/Part_2_Geostatistical_Modelling/Project/TIFF_files/N1_SER_NDVI_2015_2.tif"
ndvi_2020 <- "C:/Tamirat_2024/IGSSA_PhD/Y1_S2/1_Geospatial_Modeling/Part_2_Geostatistical_Modelling/Project/TIFF_files/N1_SER_NDVI_2020_2.tif"
ndvi_2024 <- "C:/Tamirat_2024/IGSSA_PhD/Y1_S2/1_Geospatial_Modeling/Part_2_Geostatistical_Modelling/Project/TIFF_files/N1_SER_NDVI_2024_2.tif"

Convert them to raster format for R

raster_2000 <- raster(ndvi_2000)
raster_2005 <- raster(ndvi_2005)
raster_2010 <- raster(ndvi_2010)
raster_2015 <- raster(ndvi_2015)
raster_2020 <- raster(ndvi_2020)
raster_2024 <- raster(ndvi_2024)

Plot each of the rasters

plot(raster_2000, main = "NDVI of South Ethiopia Regional State in 2000", xlab = "Longitude", ylab = "Latitude")

plot(raster_2005, main = "NDVI of South Ethiopia Regional State in 2005", xlab = "Longitude", ylab = "Latitude")

plot(raster_2010, main = "NDVI of South Ethiopia Regional State in 2010", xlab = "Longitude", ylab = "Latitude")

plot(raster_2015, main = "NDVI of South Ethiopia Regional State in 2015", xlab = "Longitude", ylab = "Latitude")

plot(raster_2020, main = "NDVI of South Ethiopia Regional State in 2020", xlab = "Longitude", ylab = "Latitude")

plot(raster_2024, main = "NDVI of South Ethiopia Regional State in 2024", xlab = "Longitude", ylab = "Latitude")