Créer un composite géomédian

Cet exercice fait suite à la section précédente. Dans la section précédente, nous avons chargé les données Sentinel-2 en utilisant load_ard(). Cela a permis d’appliquer un masque de nuage aux données en identifiant et en supprimant les points de données nuageux. Nous pouvons maintenant utiliser ces données pour créer un composite géomédiatique.

Dans cette section, nous allons réexécuter le notebook de la section précédente, et l’étendre pour créer un composite géomédien. Tout d’abord, nous allons ouvrir et exécuter notre notebook Geomedian_composite.ipynb. Cela permettra de s’assurer qu’il charge les données Sentinel-2 masquées par les nuages. Ensuite, nous utiliserons une fonction appelée xr_geomedian pour calculer les géomédians. Cela nous permettra de tracer une image composite géomédiane sans nuage.

Note

Nous allons utiliser le notebook que nous avons créé dans la section précédente, Cloud masking with load_ard(). Si vous n’avez pas encore mis en place un notebook appelé Geomedian_composite.ipynb avec les paquets et fonctions nécessaires, suivez les instructions de la section précédente. Assurez-vous d’avoir effectué toutes les étapes, y compris le chargement du jeu de données Sentinel-2 pour la côte ghanéenne.

Ouvrir et exécuter le notebook

Si vous suivez directement à partir de la dernière section, vous pouvez ignorer cette étape. Si vous avez fermé l’onglet de votre navigateur Sandbox ou vous êtes déconnecté d’Internet entre les exercices, suivez ces étapes pour vous assurer que les importations de packages et la connexion au datacube sont correctes.

  1. Accédez au dossier Training.

  2. Double-cliquez sur Geomedian_composite.ipynb. Il s’ouvrira dans le Lanceur.

  3. Select Kernel -> Restart Kernel and Run All Cells…

  4. Lorsque vous y êtes invité, sélectionnez Restart.

L’exécution du notebook peut prendre un peu de temps. Vérifiez que toutes les cellules ont été exécutées avec succès sans message d’erreur.

Vérifiez le jeu de données Sentinel-2

Confirmons que nous avons chargé avec succès notre jeu de données Sentinel-2. Rappelez-vous que nous l’avons nommé sentinel_2_ds. Faites défiler jusqu’au bas du carnet de notes et exécutez le nom du jeu de données dans la prochaine cellule vide.

image7

La sortie devrait montrer que nous avons 4 pas de temps.

Note

Si vous rencontrez un message d’erreur NameError, vous devrez peut-être recharger vos paquets et fonctions. Suivez les étapes ci-dessus, sous Open and run notebook.

Créer un composite géomédian

Maintenant nous allons créer un composite géomédian en utilisant la fonction xr_geomedian. Dans la prochaine cellule vide ci-dessous, entrez le code suivant, puis exécutez-le.

image8

Cela génère un jeu de données appelé geomedian_composite qui contient une valeur géomédienne pour chaque pixel.

Créer une image RVB du composite

Dans la nouvelle cellule ci-dessous, entrez le code suivant, puis exécutez-le pour générer une image RVB du composite.

image9

L’image RVB ressemblera à ceci:

image10

Excellent! Il n’y a pas de nuages ou d’ombres de nuages dans le composite géomédien.

Il y a quelques anomalies visuelles, mais la majeure partie de l’image représente fidèlement la zone. Notez que cette image est une représentation de ce à quoi ressemblait la zone sur la plage de temps chargée, donc du 1er février 2018 au 15 mars 2018.

Note

Enregistrez votre travail en cliquant sur Fichier -> Enregistrer le carnet ou en appuyant sur Ctrl + S.

Conclusion

Nous avons créé un composite géomédien de données Sentinel-2 et vu qu’il s’agit d’une représentation précise d’une zone. Il remplit les données manquantes ou masquées, telles que les nuages et les ombres des nuages.

C’est la fin de la session 3. Dans cette session, vous avez été initié au concept d’utilisation d’ensembles de données et d’images composites pour combler les lacunes de données. Nous avons également appris que la meilleure méthode statistique pour combiner des ensembles de données est le géomédian. Bon travail!