Méthode de cartographie multi-espèces d’une zone comprenant les étapes de : Acquérir (401) un ensemble de données spectrales associant pour chaque point de la zone un spectre, Appliquer (402) aux données spectrales une transformée en Ondelettes de manière à produire un ensemble de coefficients comprenant des coefficients de détail caractéristiques du bruit et des coefficients d’approximation caractéristiques du signal, Mettre à zéro (403) les coefficients de détail dont la valeur de la densité spectrale d’énergie est inférieure à un seuil prédéterminé, Appliquer (404) une transformée en Ondelettes inverse à l’ensemble de coefficients pour obtenir un nouvel ensemble de données spectrales parcimonieuses Appliquer (405) une analyse en composantes principales au nouvel ensemble de données spectrales parcimonieuses de manière à : Déterminer une base de projection L dans le domaine spectral, chaque composante de la base de projection correspondant à une signature spectrale d’une espèce, Projeter les données spectrales dans la base de projection L pour obtenir, pour chaque composante, une carte de la zone indiquant, pour chaque point de la zone, la proportion de l’espèce. Méthode de cartographie multi-espèces selon la revendication 1 dans laquelle l’étape (405) d’analyse en composantes principales est réalisée à partir d’une matrice de covariance filtrée à l’aide d’une fonction noyau déterminée de manière à augmenter le rapport signal à bruit des données spectrales sur lesquelles est calculée la matrice de covariance. Méthode de cartographie multi- espèces selon la revendication 2 dans laquelle la fonction noyau est du type , avec α un paramètre prédéterminé positif. Méthode de cartographie multi- espèces selon la revendication 3 dans laquelle la valeur du paramètre α dépend du rapport signal à bruit des données spectrales. Méthode de cartographie multi- espèces selon l’une quelconque des revendications précédentes dans laquelle la transformée en Ondelettes est une transformée en Ondelettes discrète, par exemple une transformée en Ondelettes bi-orthogonale. Méthode de cartographie multi- espèces selon l’une quelconque des revendications précédentes dans laquelle le seuil prédéterminé est compris entre 90% et 99% de la valeur maximale de la densité spectrale d’énergie des coefficients de détail. Méthode de cartographie multi- espèces selon l’une quelconque des revendications précédentes comprenant en outre une étape d’association de chaque composante de la base de projection L à une espèce en fonction d’une comparaison avec une signature spectrale de l’espèce. Méthode de cartographie multi- espèces selon l’une quelconque des revendications précédentes dans laquelle les espèces sont des éléments chimiques ayant chacun une signature spectrale donnée. Méthode de cartographie multi- espèces selon l’une quelconque des revendications précédentes dans laquelle les données spectrales sont des données hyperspectrales. Dispositif de cartographie multi- espèces comprenant un appareil d’acquisition d’un ensemble de données spectrales et une unité de calcul configurée pour exécuter les étapes de la méthode selon l’une quelconque des revendications précédentes. Programme d’ordinateur comprenant des instructions qui conduisent le dispositif selon la revendication 10 à exécuter les étapes de la méthode selon l’une quelconque des revendications 1 à 9. Support lisible par ordinateur sur lequel est enregistré le programme d’ordinateur selon la revendication 11.