Logiciel :

• QGIS 3.28
  
• Plugins :
  • Terrain Profile
  • Qgis2threejs Exporter

Dossier SVN :

• ANS-PROC\Documents\CARTO\Templates\Scripts\DATA_PROCESSING\AMSR

Vecteur (.shp) :

• AMSR (entités polygonales), champs requis respectivement : MOC | SECTOR | BUFFER | NAME

⇔ minimum obstacle clearance (m) | Altitude max Obstacle pour le secteur de base (m) | distance tampon (emprise secondaire, NM) | nom du secteur

Plugin QGIS :

• amsr_amg_rasterisation (dossier SVN)

Méthode manuelle (caduc)

Feuille de calcul :

• MODEL_PARAMETERS.xlsx (dossier SVN)

Modèle de script :

• AMSR_LIMITES_POLY.model3 (dossier SVN)

Script Python :

• AMSR_SURFACE_POLY.py (dossier SVN)

Rastérisation d'une surface d'Altitudes Minimales de Sécurité Radar (AMSR) dans le but d'effectuer des calculs (Calculatrice Raster) entre l'AMSR et un modèle numérique de terrain (MNT).

Cette modélisation s'effectue en plusieurs étapes à l'aide d'un plugin QGIS.

La partie ci-dessous concerne une méthode manuelle qui est désormais désuète excepté la partie pour les résultats secondaires…

Cette modélisation s'effectue en plusieurs étapes :

1. AMSR_LIMITES_POLY.model3 :
   • découplage vectoriel des secteurs et génération d'un tampon (emprise secondaire) de distance pouvant varier pour chaque secteur
2. AMSR_SURFACE_POLY.py :
   • génération d'une surface d'interpolation continue 2.5D (raster) pour chacun des secteurs de base en incluant la Marge de Franchissement d'Obstacles (MFO)
   • découpage de chaque raster selon leur emprise secondaire respective
   • génération du raster de la surface AMSR (valeur minimale conservée pour chaque empilement de pixels)
   • calcul des résidus de l'ajustement de la surface modélisée
     

1. Calcul des correctifs de température (référent : Loïc GIROUD)
2. Explorateur Windows :
   • créer un DOSSIER nommé “code-OACI”_AMSR
   • créer un sous dossier nommé LIMITES_EMPRISE
   • copier les chemins du dossier et du sous dossier
3. Calcul des paramètres de l'emprise géographique pour la modélisation (MODEL_PARAMETERS.xlsx) :

• NORME fixe
• NORME variable
• MODEL_PARAMETERS

Sous QGIS

1. AMSR_LIMITES_POLY.model3
2. AMSR_SURFACE_POLY.py
   
3. Vérification de l’ajustement de la surface raster :

• Résidus
• Transect

1. Sous dossier LIMITES_EMPRISE : le modèle de script génère des fichiers .shp ET un raster de l'emprise des aires secondaires
2. Dossier “code-OACI_AMSR”, trois sous dossiers ainsi que le raster de la surface AMSR sont créés, à la racine du dossier principal, lors de l'exécution du script Python :
   • RASTERS_TINs : rasters d'interpolation linéaire, des secteurs de bases à leur emprise secondaire, en prenant en compte la MFO
   • CLIPS : découpage des RASTERS_TINs selon l'emprise vectorielle des aires secondaires
   • RSD_clipped_TINs : calcul des résidus d'ajustement de la surface AMSR par rapport aux différents secteurs, les résidus doivent être ≤ 0

• AMSR_SURFACE.tif

• HEIGHT DIFFERENCE AMSR-MNT : la différence de hauteur entre la surface rasterisée de l'AMSR et un MNT peut être calculée sous QGIS, ou tout autre SIG, à l'aide de la Calculatrice Raster

"AMSR_SURFACE@1" - "LFPM_MNT_75M_LAMB93_RGF93_IGN69@1"

• Statistiques de couche raster ('MIN', 'MAX' VALUES) : sous SIG les Statistiques de la couche raster HEIGHT DIFFERENCE AMSR-MNT peuvent être calculées et sont nécessaires afin d'établir une table de reclassification
  

• Reclassification par table :

Limites de plages : min < = valeur < max

• Fichier .prj (WGS 84) : nécessaire afin d'importer les rasters dans GTT
• Polygoniser : après reclassification du raster HEIGHT DIFFERENCE AMSR-MNT, il convient d'effectuer une conversion Polygoniser (raster vers vecteur) afin de calculer les superficies de chaque classe
• Calculatrice de champ (fonction Géométrie : $area), il est nécessaire d'ajouter un champ (e.g. AREA_km², nombre réel) via la calculatrice de champ selon la syntaxe suivante :

$area / 10^6 

• Créer une Nouvelle couche virtuelle… et effectuer une Requête SQL sur AREA :

SELECT DN, SUM(AREA_km) FROM TABLE_VIRTUAL_LAYER_NAME GROUP BY "DN";