Rapport hebdomadaire 11

Author

Alix Biton

Published

March 27, 2026

Semaine 11 (23-27/03)

Objectifs

  • Etablir des petites zones d’études
  • Faire le masque de forêt
  • Faire le modèle type Kurt McLaren

Masque de forêt

Les masques de forêts disponibles grâce aux travaux de Thani permettent d’avoir une petite zone d’étude de “forêt publique” ou de “parc ou réserve”. J’ai en premier lieu fait avec seulement “forêt publique” et, couplé aux données NDVI disponibles en décembre et janvier, il ne reste seulement quelques zones disponibles.

J’ai pu en extraire quelques données quantitatives récapitulatives : - Altitude - Longitude - Latitude - Vitesse Maximale au cours du cyclone - Exposition topographique (pour différents angles allant de 6 à 20°) - Pente - Différence de NDVI (janvier - décembre)

Grâce à cela, quelques métriques et premiers modèles émergent.

Modèle type Kurt McLaren

Le modèle utilisé par Kurt McLaren (https://doi.org/10.1016/j.agrformet.2019.107621) est comme un indice de vulnérabilité d’exposition décrit comme ceci :\((sum(vitesse_i * exposure_i))/n\) avec n le nombre total de lieux d’intérêt et i l’indice d’un de ces lieux.

J’ai donc calculé cet indice tel que i soit un pixel et n le nombre de pixels disponibles.

De plus, dans cet article, les auteurs mentionnent un angle de 20° (angle d’ouverture de canopée) qui serait utilisé aussi en Jamaique par Boose et al. (1994). J’ai donc regardé le % de zones exposées, de zones protégées ainsi que le % de dégâts en fonction de la zone dans laquelle elle se situe (exposées ou protégée)

# compute EV as Kurt McLaren
n <- nrow(dta)
ev <- sum(dta$msw_max * dta$ang6, na.rm = TRUE) / n
print(ev)
[1] 6.689356
ev/mean(dta$ang6)
[1] 53.40212
mean(dta$msw_max)
[1] 53.11845

Cela suggère que la vitesse maximale moyenne du cyclone (MSW) sur la zone est d’environ 53.02 m/s. La moyenne des vitesses maximale est de 52.69 m/s donc globalement, l’exposition a accentué l’impact du cyclone sur la zone.

# model kurt mclaren 
mkurt <- lm(perte_ndvi ~ msw_max * ang6, data = dta)
summary(mkurt)

Call:
lm(formula = perte_ndvi ~ msw_max * ang6, data = dta)

Residuals:
     Min       1Q   Median       3Q      Max 
-0.63211 -0.10092  0.00021  0.09086  0.75558 

Coefficients:
               Estimate Std. Error t value Pr(>|t|)    
(Intercept)   1.0559631  0.0195205  54.095   <2e-16 ***
msw_max      -0.0225362  0.0003673 -61.353   <2e-16 ***
ang6          0.7259099  0.0844258   8.598   <2e-16 ***
msw_max:ang6 -0.0167696  0.0015775 -10.630   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.1612 on 21767 degrees of freedom
Multiple R-squared:  0.2432,    Adjusted R-squared:  0.2431 
F-statistic:  2331 on 3 and 21767 DF,  p-value: < 2.2e-16
  • effet de la vitesse seule : À chaque fois que le vent augmente de 1 m/s , le NDVI diminue de 0.03 unités.
  • effet de l’exposition topo seul : pas trop interprétable car considère un vent nul ce qui n’arrive pas lors du cyclone.
  • effet de l’interaction : L’effet destructeur du vent est accentué sur les versants exposés. Sur un versant face au vent (ang6 = 1), la dégradation du NDVI n’est plus de −0.030, mais elle s’additionne avec l’interaction : −0.030+(−0.026)=−0.056 donc à vitesse de vent égale, la perte de végétation est deux fois plus forte sur un versant exposé que sur un terrain plat.

Zones d’études possibles ?

Il faut que les zones d’études soient une zone forestière avec une différence de NDVI marquée c’est-à-dire significative une exposition topographique plutôt marquée (altitude élevée, plus que pour Mbouzi).