G Exercise 05

Last updated: 2025-06-20 13:02:14

G.1 Materials

The materials for this exercise are:

  • Geometric operations with rasters (Chapter 9)
  • Combining rasters and vector layers (Chapter 10)

G.2 Question 1

  • Recreate the DEM of the Carmel area in UTM, as shown in the right panel in Figure 9.10, using the following code section:
library(stars)
dem1 = read_stars('data/srtm_43_06.tif')
dem2 = read_stars('data/srtm_44_06.tif')
pol = st_read('data/israel_borders.shp', quiet = TRUE)
dem = st_mosaic(dem1, dem2)
dem = dem[, 5687:6287, 2348:2948]
names(dem) = 'elevation'
dem = st_warp(src = dem, crs = 32636, method = 'near', cellsize = 90)
  • Calculate a point geometry of the peak of the Carmel mountain, by finding the centroid of the pixel that has the highest elevation in the image.
  • Calculate the nearest point in the Mediterranean sea from the peak, assuming the sea area is all of the NA pixels in the image. (Hint: calculate the distances from the peak to all NA pixels, converted to points, then use which.min to find the nearest one.)
  • Plot the DEM, the peak, the nearest “sea” point, the line connecting them, and the distance, as shown in Figure G.1.
Length of the shortest line between Carmel mountain peak and the Mediterranean sea

Figure G.1: Length of the shortest line between Carmel mountain peak and the Mediterranean sea

(50 points)

G.3 Question 2

  • Read the 'MOD13A3_2000_2019.tif' raster, with monthly NDVI values, into a stars object named r, then reproject it to UTM with st_warp(r,crs=32636).
r = read_stars('data/MOD13A3_2000_2019.tif')
r = st_warp(r, crs = 32636)
  • Use 'MOD13A3_2000_2019_dates.csv' to assign date values to raster bands (see Section 6.3.2).
  • Read the Shapefile named 'nafot.shp', which includes polygons of “Nafa” administrative regions in Israel, and reproject it to UTM as well.
  • Calculate the average NDVI for each polygon in each date.
  • Plot the resulting time series of average NDVI for each of the “Nafot”, as shown in Figure G.2. Hint: you can start from an empty plot using plot(..., col=NA), then add the NDVI times series using lines in for loop; see Section 6.1.3.
  • Mark the time series of “Nafa” where the lowest average NDVI was observed, during the studied period, in red. The other “Nafa” should be displayed in grey. Don’t use “Nafa” indices or names, or any other specific values; instead calculate the required indices in your code.
Average NDVI over time in "Nafa" polygons, the "Nafa" with the lowest observed average NDVI is marked in red

Figure G.2: Average NDVI over time in “Nafa” polygons, the “Nafa” with the lowest observed average NDVI is marked in red

(50 points)