G Exercise 05
Last updated: 2021-01-16 22:29:08
G.1 Question 1
- Recreate the DEM of the Carmel area in UTM, as shown in the right panel in Figure 9.10.
- 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
NApixels in the image. (Hint: calculate the distances from the peak to all
NApixels, converted to points, then use
which.minto 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.
G.2 Question 2
- Read the
MOD13A3_2000_2019.tifraster, with monthly NDVI values, into a
rto UTM with
st_warp(r, crs = 32636).
- 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
forloop; see Section 6.1.3.
- The time series of “Nafa” where the highest and lowest average NDVI were observed, during the studied period, need to be highlighted in blue and red, respectively. The other “Nafa” should be displayed in grey. Don’t specific use “Nafa” indices (such as
6) or names (such as
"Golan"), instead calculate the required indices in your code.