Skip to main content

Didacticiel : Créer des produits d’imagerie satellite via l’orthocartographie dans ArcGIS Pro

Disponible avec une licence Advanced.

Dans ArcGIS Pro, vous pouvez corriger les images satellite de manière photogrammétrique afin de supprimer toute distorsion géométrique induite par la plateforme et le déplacement de terrain. Une fois les distorsions supprimées, vous pouvez générer des produits d’orthocartographie, notamment une orthomosaïque, un modèle numérique de surface (MNS) et un maillage MNS. Dans ce didacticiel, vous allez générer un MNS et une orthomosaïque haute résolution.

Tout d’abord, vous configurerez un espace de travail d’orthocartographie pour gérer une collection d’images satellite. Ensuite, vous effectuerez un ajustement de bloc, puis un ajustement affiné à l’aide de points de contrôle au sol. Enfin, vous générerez une mosaïque orthorectifiée (orthomosaïque) et un MNS.

ArcGIS Pro peut traiter des images satellite à partir de plusieurs plateformes de capteur tant que l’orientation des images est décrite par un modèle de coefficients polynomiaux rationnels (RPC) ou un modèle de capteur rigoureux. Ce modèle est généralement incorporé dans le fichier image ou inclus comme fichier de métadonnées séparé.

Créer un espace de travail d’orthocartographie.

Un espace de travail d’orthocartographie est un sous-projet ArcGIS Pro dédié aux processus d’orthocartographie. Il s’agit d’un conteneur au sein d’un dossier de projet ArcGIS Pro qui contient les ressources et les fichiers dérivés appartenant à une collection d’images dans une tâche d’orthocartographie.

L’imagerie empaquetée pour ce didacticiel a été collectée et fournie par Maxar Technologies. Elle comprend une paire d’images multispectrales et panchromatiques, une table de points de contrôle au sol et un MNE.

Pour créer un espace de travail d’orthocartographie, procédez comme suit :

  1. Téléchargez le jeu de données du didacticiel et enregistrez-le dans C:\SampleData\orthomapping_satellite_tutorial.

  2. Dans ArcGIS Pro, créez un projet en utilisant le modèle Map (Carte).

  3. Dans l’onglet Imagery (Imagerie), dans le groupe Ortho Mapping (Orthocartographie), cliquez sur la liste déroulante New Workspace (Nouvel espace de travail) et sélectionnez New Workspace (Nouvel espace de travail).

  4. Dans la fenêtre Workspace Configuration (Configuration de l’espace de travail), saisissez un nom pour l’espace de travail.

  5. Dans le menu déroulant Type, sélectionnez Satellite.

  6. Si vous le souhaitez, vous pouvez choisir dans la liste déroulante Basemap (Fond de carte) un fond de carte comme arrière-plan de la collection d’images.

  7. Définissez la valeur Parallel Processing Factor (Facteur de traitement parallèle) de l’espace de travail sur 90 %.

    Définir la valeur Parallel Processing Factor (Facteur de traitement parallèle) sur 90 % signifie que 90 % du nombre total de cœurs d’UC sont utilisés pour prendre en charge le traitement d’orthocartographie. Il est possible également de changer la valeur Parallel Processing Factor (Facteur de traitement parallèle) sous Properties (Propriétés) dans le groupe Workspace (Espace de travail) de l’onglet Imagery (Imagerie), une fois l’espace de travail créé.

    Définir la valeur du paramètre Parallel Processing Factor (Facteur de traitement parallèle) sur une valeur supérieure à celle que le système peut prendre en charge entraîne un échec du traitement parallèle. La configuration couramment requise pour chaque processeur logique correspond à 2 Go de RAM. Si, par exemple, vous utilisez un système qui dispose de 6 cœurs, 12 processeurs logiques et 16 Go de RAM, définir le facteur de traitement parallèle sur 100 pour cent implique au minimum 24 Go de RAM pour le bon déroulement du traitement. Une valeur plus appropriée du paramètre Parallel Processing Factor (Facteur de traitement parallèle) basée sur cet exemple serait 50 pour cent, valeur qui nécessiterait environ 12 Go de RAM.

  8. Cliquez sur Next (Suivant).

    Assistant Nouvel espace de travail d’orthocartographie

  9. Dans la fenêtre Image Collection (Collection d’images), pour Sensor Type (Type de capteur), sélectionnez GeoEye-1.

  10. Sous Folder Containing Images (Dossier contenant les images), cliquez sur le bouton Browse (Parcourir) et accédez au dossier qui contient les données du didacticiel sur votre machine, sélectionnez le dossier d’images (GeoEye_NAD83_UTM10N), puis cliquez sur OK.

    Les informations appropriées sont automatiquement indiquées dans Workspace Spatial Reference (Référence spatiale de l’espace de travail).

    Fenêtre Collection d’images

  11. Cliquez sur Next (Suivant).

  12. Dans la fenêtre Data Loader (Chargeur de données), sous Elevation Source (Source d’élévation), sélectionnez DEM (MNE). Sous DEM (MNE), accédez au MNE fourni avec le jeu de données du didacticiel.

    Remarque :

    Ce MNE sera utilisé pour prendre en charge le traitement d’ajustement de bloc.

  13. Comme la plupart des données d’altitude utilisent des hauteurs orthométriques, vous devez appliquer une correction de géoïde. Sous Geoid correction (Correction géodésique), vérifiez que EGM96 est sélectionné.

  14. Sous Processing Template (Modèle de traitement), sélectionnez Panchomatic (Panchromatique).

  15. Développez Advanced Options (Options avancées).

  16. Vérifiez que l’option Estimate Statistics (Estimer les statistiques) est cochée.

  17. Acceptez tous les autres paramètres par défaut et cliquez sur Finish (Terminer).

    Chargeur de données satellite

    Une fois l’espace de travail créé, les images et les emprises des images sont affichées dans la carte. Une catégorie Ortho Mapping a également été ajoutée à la fenêtre Contents (Contenu), où les données d’imagerie source et les produits d’orthocartographie dérivés seront stockés.

    L’affichage initial de l’imagerie dans l’espace de travail confirme que toutes les images et les métadonnées nécessaires ont été fournies pour lancer l’espace de travail. Les images n’ayant pas été alignées ni ajustées, l’aspect de la mosaïque peut sembler incorrect du point de vue géométrique.

    Un nouvel onglet Ortho Mapping (Orthocartographie) est ajouté au menu principal ArcGIS Pro. Un clic sur cet onglet fait apparaître un ensemble d’outils et de processus dédiés à l’orthocartographie. Dans la catégorie Product (Produit), tous les boutons sont indisponibles car les images n’ont pas encore été ajustées.

    Espace de travail Satellite

    Image satellite © 2020 Maxar Technologies

Effectuer un ajustement de bloc

Une fois que vous avez créé l’espace de travail d’orthocartographie, l’étape suivante consiste à effectuer un ajustement de bloc à l’aide des outils des groupes Adjust (Ajuster) et Refine (Affiner). L’ajustement de bloc commence par calculer les points de rattachement, qui sont des points communs dans les zones de la superposition d’images. Les points de rattachement sont alors utilisés pour calculer l’orientation de chaque image, appelée "orientation extérieure" en photogrammétrie.

  1. Dans l’onglet Ortho Mapping (Orthocartographie), dans le groupe Adjust (Ajuster), cliquez sur Adjust (Ajuster) Ajuster.

  2. Dans la fenêtre Adjust (Ajuster), sous Transformation Type (Type de transformation), choisissez RPC. La transformation des coefficients polynomiaux rationnels (RPC) est appliquée à l’ajustement, qui est utilisé pour les images satellite qui contiennent des données RPC dans les métadonnées.

  3. Cochez la case en regard de Reproject Tie Points (Reprojeter les points de rattachement).

    Cela permet de garantir que les coordonnées de la carte des points de rattachement sont calculées.

  4. Développez la section Tie Point Matching (Appariement des points de rattachement) et vérifiez que la configuration des paramètres correspond à celle de l’exemple ci-dessous.

    Option d’ajustement de bloc pour l’imagerie satellite

  5. Cliquez sur Run (Exécuter) pour effectuer l’ajustement de bloc.

  6. Une fois l’ajustement effectué, activez la couche Points de rattachements dans la fenêtre Contents (Contenu) pour afficher la distribution des points de rattachement générés sur la carte. Votre distribution des points de rattachement peut être différente de celle illustrée ci-dessus.

    Points de rattachement distribués sur l’image satellite

  7. Les valeurs résiduelles et un rapport d’exactitude des points de rattachement peuvent être consultés dans le fichier des journaux. Dans l’onglet Ortho Mapping (Orthocartographie), dans le groupe Review (Examiner), cliquez sur Logs (Journaux) Affichage des journaux pour accéder au fichier Logs (Journaux). Les valeurs résiduelles des points de rattachement sont affichées sur la ligne intitulée RMSE_Tie_Image(x,y). L’EQM des points de rattachement est exprimée en pixels.

    Vos valeurs résiduelles peuvent différer légèrement de celles de l’exemple ci-dessus.

    Journal des points de rattachement

Suite à l’ajustement initial, remarquez que tous les boutons dans la catégorie Product (Produit) sont devenus actifs. Ces boutons mettent en évidence les produits d’imagerie qui peuvent être générés. Avant de générer des produits, les points de contrôle au sol vont être utilisés pour améliorer la précision absolue des images.

Ajouter des points de contrôle au sol (PCS)

Les points de contrôle au sol (PCS) sont des points avec des coordonnées au sol x,y,z connues. Ils sont souvent obtenus à partir d’arpentage au sol ou de données existantes et permettent de garantir le géoréférencement précis des images dans le système de coordonnées au sol. Il est possible d’appliquer l’ajustement de bloc sans points de contrôle au sol et de garantir une précision relative, mais l’ajout de points de contrôle au sol augmente la précision absolue de l’imagerie ajustée. Si les points de contrôle au sol ne proviennent pas d’un arpentage au sol, mais que vous disposez d’une couche raster géoréférencée (jeu de données raster, jeu de données mosaïque ou service d’imagerie), vous pouvez l’ajouter en tant que référence pour le calcul des points de contrôle au sol.

Les points de contrôle au sol sélectionnés dans une couche raster et enregistrés dans un fichier texte peuvent également être importés et utilisés pour améliorer l’exactitude de l’ajustement. Dans le cadre de ce didacticiel, vous allez utiliser cette méthode d’importation pour ajouter des points de contrôle au sol au projet

Pour importer des points de contrôle au sol, procédez comme suit :

  1. Dans l’onglet Ortho Mapping (Orthocartographie), dans le groupe Refine (Affiner), cliquez sur Manage GCPs (Gérer les PCS) pour ouvrir la fenêtre GCP Manager (Gestionnaire de PCS).

  2. Dans la fenêtre GCP Manager (Gestionnaire de PCS), cliquez sur le bouton Import GCPs (Importer des PCS) Importer des PCS.

  3. Dans la fenêtre Import GCPs (Importer des PCS), accédez au fichier de PCS (Vancouver_NAD83-UTM10N.csv) et sélectionnez-le. Cliquez sur OK.

  4. Sous Set GCP Spatial Reference (Définir la référence spatiale des PCS), cliquez sur le bouton Browse (Parcourir) Parcourir la référence spatiale. Pour Horizontal systems (Systèmes horizontaux), développez Projected Coordinate System (Système de coordonnées projetées), North America (Amérique du Nord), Zone Series (Série de zones), UTM (NAD 1983), puis choisissez NAD 1983 UTM Zone 10N. Cliquez sur OK pour accepter les modifications et fermer la fenêtre Spatial Reference (Référence spatiale).

    Le système de coordonnées verticales n’a pas été défini, car il n’existe pas de système de ce type défini pour le modèle d’élévation numérique utilisé pour extraire les valeurs de hauteur des PCS. S’il avait existé un système de coordonnées verticales pour le MNE que vous avez utilisé, les systèmes verticaux auraient été définis avec les coordonnées d’appariement.

  5. Pour Geographic Transformations (Transformations géographiques), cliquez sur l’onglet Horizontal et sélectionnez WGS 1984 (ITRF00) To NAD 1983 (WGS 1984 [ITRF00] vers NAD 1983) dans la liste déroulante.

  6. Vérifiez que l’appariement des champs est correct.

  7. Cliquez sur le bouton Browse (Parcourir) en dessous de GCP Photo Location (Localisation des photos des PCS), accédez au dossier contenant les images des localisations des PCS et sélectionnez-le. Cliquez sur OK.

  8. Pour Photo File Extension (Extension du fichier photo), choisissez PNG dans la liste déroulante.

    Importer des PCS pour améliorer l’ajustement des images satellite

  9. Cliquez sur OK pour importer les PCS.

    Une fois les PCS importés, la table du Gestionnaire de PCS est alimentée.

    PCS ajoutés aux projets

    Dans la fenêtre GCP Manager (Gestionnaire de PCS), le premier PCS de la liste est sélectionné par défaut et un aperçu de l’image de sa localisation s’affiche dans la section Preview (Aperçu).

  10. Pour ajouter des points de rattachement pour le PCS sélectionné, dans la section Preview (Aperçu), cliquez sur le bouton View GCP Photo (Afficher la photo du PCS) Afficher la photo du PCS pour afficher le fragment d’image du PCS. Utilisez la molette de la souris pour zoomer légèrement sur le fragment d’image et afficher la localisation du PCS indiqué par une flèche rouge.

  11. Cliquez sur le menu déroulant à côté du bouton Add GCP or Tie Point (Ajouter un PCS ou un point de rattachement) Ajouter des points de rattachement dans la fenêtre GCP Manager (Gestionnaire de PCS) et sélectionnez Semi-Auto (Semi-Automatique) pour ajouter un point de rattachement dans la visionneuse de chaque image.

    Les points de rattachement des autres images sont automatiquement calculés par l’algorithme d’appariement d’image lorsque cela est possible, mais il convient de vérifier l’exactitude de chaque point de rattachement. Si le point de rattachement n’est pas identifié automatiquement, ajoutez-le manuellement en sélectionnant l’emplacement approprié dans l’image.

  12. Répétez les étapes 10 et 11 pour sélectionner et ajouter des points de rattachement pour les PCS restants.

  13. Une fois les PCS mesurés avec les points de rattachement, sélectionnez le point PP_GCP02 et cliquez dessus avec le bouton droit de la souris pour le modifier en Check Point (Point de vérification). Ce point sera exclu du traitement d’ajustement et servira à évaluer indépendamment la précision des résultats d’ajustement.

  14. Après avoir ajouté les PCS et les points de vérification, l’ajustement doit être réexécuté afin d’intégrer ces points. Cliquez sur Adjust (Ajuster).

    Gestionnaire de PCS

Examiner les résultats de l’ajustement

La qualité des résultats de l’ajustement peut être consultée dans la fenêtre GCP Manager (Gestionnaire de PCS) en analysant les résidus de chaque PCS. Les valeurs résiduelles représentent la différence entre la position mesurée et la position calculée d’un point. Elles sont mesurées dans les unités du système de référencement spatial du projet. Après l’ajustement avec des PCS, trois nouveaux champs (dX, dY et dZ) affichant les valeurs résiduelles de chaque PCS sont ajoutés à la table GCP Manager (Gestionnaire de PCS). La qualité de l’adéquation entre le bloc ajusté et le système de coordonnées de la carte peut être évaluée à l’aide de ces valeurs. L’erreur quadratique moyenne (EQM) des résidus peut être affichée en développant la section Residual Overview (Vue d’ensemble des résidus) de la fenêtre GCP Manager (Gestionnaire de PCS).

Examiner les résultats de l’ajustement satellite

Des statistiques d’ajustement supplémentaires sont fournies dans le rapport de l’ajustement. Pour générer ce rapport, dans l’onglet Ortho Mapping (Ortho-cartographie), dans le groupe Review (Examiner), cliquez sur Adjustment Report (Rapport d’ajustement).

Générer un modèle numérique de surface (MNS)

Une fois l’ajustement de bloc effectué, des produits d’imagerie 2D peuvent être générés à l’aide des outils du groupe Product (Produit) de l’onglet Reality Mapping. Les paires d’images stéréo d’une collecte d’images permettent de générer un nuage de points (points 3D) pour lequel des données d’élévation peuvent être calculées. Les données d’élévation dérivées sont classées soit comme un modèle numérique de terrain (MNT) qui inclut la surface du sol, soit comme un modèle numérique de surface (MNS) qui inclut l’élévation des arbres, des bâtiments et d’autres entités au-dessus du sol.

Remarque :

Si une surface est fortement boisée ou présente une autre couverture végétale dense, il n’est pas possible de déduire la surface au sol du MNT car le sol n’est pas visible. Ce problème peut également survenir dans les zones urbaines denses, où les bâtiments cachent le sol. Dans ce cas, la surface d’élévation la plus appropriée est un MNS, qui crée spécifiquement une surface représentant la partie supérieure de l’environnement urbain.

Pour générer un MNS à l’aide de l’assistant, procédez selon les étapes suivantes :

  1. Dans l’onglet Ortho Mapping (Orthocartographie), cliquez sur le bouton DSM (MNS) dans le groupe Product (Produit).

    La fenêtre Ortho Mapping Products Wizard (Assistant de produits d’orthocartographie) s’affiche.

  2. Cliquez sur Next (Suivant) pour faire passer l’assistant à la fenêtre Point Cloud Settings (Paramètres du nuage de points).

  3. Dans la fenêtre Point Cloud Settings (Paramètres du nuage de points), pour Matching Method (Méthode d’appariement), choisissez Semiglobal Matching (Appariement semi-global) dans le menu déroulant.

    En règle générale, cette méthode est utilisée pour les images de zones urbaines et permet de capturer des informations de terrain plus détaillées.

  4. Acceptez la valeur Point Ground Spacing (Espacement au sol des points) par défaut.

    Il définit l’espacement, en mètres, auquel sont générés les points 3D. La valeur par défaut est égale à trois fois la résolution de l’imagerie source.

  5. Acceptez tous les autres paramètres par défaut, puis cliquez sur Next (Suivant).

    Pour obtenir des informations sur Advanced Settings (Paramètres avancés), voir Créer des données d’élévation à l’aide de l’assistant MNT d’ortho-cartographie.

    Paramètres du nuage de points

  6. Dans la fenêtre DEM Settings (Paramètres MNE), pour Cell Size (Taille de cellule), utilisez la valeur par défaut égale à 3 x GSD.

    Cela détermine la résolution du MNS qui correspond à trois fois la résolution de l’imagerie, dans notre cas de figure.

  7. Acceptez les autres paramètres par défaut, puis cliquez sur Finish (Terminer).

    Le MNS est généré.

Générer une orthomosaïque

Vous allez à présent générer une orthomosaïque. Une orthomosaïque est un produit d’imagerie orthorectifié et mosaïqué à partir d’un ensemble d’images. La distorsion géométrique a été corrigée et les couleurs des images ont été équilibrées de manière à obtenir une mosaïque.

  1. Dans l’onglet Ortho Mapping (Orthocartographie), dans le groupe Product (Produit), cliquez sur Orthomosaic (orthomosaïque) pour lancer l’assistant d’orthomosaïque.

  2. Vérifiez que les options Color Balance (Équilibrer les couleurs) et Generate Seamlines (Générer les lignes de raccord) sont sélectionnées.

  3. Cliquez sur Next (Suivant).

  4. Dans la fenêtre Color Balance Settings (Paramètres d’équilibrage des couleurs), pour Balance Method (Méthode d’équilibrage), sélectionnez Dodging (Éclaircissement) dans la liste déroulante et acceptez toutes les autres valeurs par défaut.

  5. Cliquez sur Next (Suivant).

  6. Dans la fenêtre Seamline Settings (Paramètres de ligne de raccord), pour Computation Method (Méthode de calcul), sélectionnez Voronoi dans la liste déroulante.

  7. Développez la section Processing (Traitement), puis entrez 20 pour Blend Width (Largeur de fusion).

    Paramètres des lignes de raccord d’orthocartographie

  8. Cliquez sur Next (Suivant).

    Le processus guidé de l’assistant passe à la fenêtre suivante, Orthomosaic Settings (Paramètres de l’orthomosaïque).

  9. Acceptez tous les paramètres par défaut dans la fenêtre Orthomosaic Settings (Paramètres de l’orthomosaïque) et cliquez sur Finish (Terminer).

    L’orthomosaïque est générée, répertoriée dans la fenêtre Contents (Contenu) et chargée dans l’affichage cartographique.

    Remarque :

    Plusieurs produits peuvent être générés simultanément à l’aide de l’assistant Plusieurs produits.

Résumé

Dans ce didacticiel, vous avez créé un espace de travail d’orthocartographie pour l’imagerie satellit et utilisé les outils de l’onglet Ortho Mapping (Orthocartographie) pour appliquer un ajustement photogrammétrique avec des points de contrôle au sol. Vous avez ensuite utilisé l’assistant d’orthocartographie MNS pour générer un MNE et une orthomosaïque.

Les images satellite utilisées dans ce didacticiel ont été acquises et fournies par Maxar Technologies.

Rubriques connexes