Procédé de traitement d’images photoacoustiques, comprenant - une obtention (110) d’une succession temporelle d’images d’un échantillon acquises par un système d’imagerie photoacoustique pour M λ longueurs d’ondes d’impulsions d’excitations avec N images acquises par longueur d’onde; - un filtrage (120) spatio-temporel multispectral par décomposition en valeurs singulières appliqué à l’ensemble des N*M λ images acquises de sorte à obtenir N*M λ images filtrées ; - pour chaque longueur d’onde, un calcul (140) d’une image de variance filtrée, un pixel de coordonnée r dans l’image de variance filtrée étant égal à la variance de la distribution des valeurs de pixels de même coordonnée r dans les images filtrées obtenues pour cette longueur d’onde ; - pour chaque longueur d’onde, une correction (150A) de l’image de variance filtrée par suppression de fluctuations parasites de sorte à obtenir une image de fluctuations corrigée qui est proportionnelle en chaque point d’image à l’énergie absorbée en un point de l’espace correspondant au point d’image concerné. Procédé selon la revendication 1 dans lequel la correction (150A) de l’image de variance filtrée comprend - une correction de fluctuations dues au bruit électronique de capteurs d’ondes ultrasonores du système d’imagerie photoacoustique par soustraction à chaque pixel de coordonnée r dans l’image de variance filtrée d’une variance du bruit électronique pour obtenir une image de variance corrigée ; - une détermination de l’image de fluctuations corrigée dont chaque pixel est égal à la racine carrée du pixel correspondant de l’image de variance corrigée. Procédé selon la revendication 2, comprenant en outre une estimation d’une variance du bruit électronique produit par des capteurs d’ondes ultrasonores sur les images, la variance de bruit électronique produit par les capteurs d’ondes ultrasonores sur les images étant estimé en fonction d’une variance du bruit électronique produit dans les signaux photoacoustiques acquis en l’absence d’échantillon corrigée par une quantité de bruit éliminée par le filtrage spatio-temporel multispectral par décomposition en valeurs singulières, la quantité de bruit étant estimée sur la base des valeurs singulières correspondant aux composantes de plus basse énergie supprimées par le filtrage spatio-temporel multispectral par décomposition en valeurs singulières. Procédé selon l’une quelconque des revendications 1 à 3, comprenant pour chaque longueur d’onde considérée, une normalisation (150B) de l’image de fluctuations corrigée par une fonction de la fluence d’impulsions laser du système d’imagerie photoacoustique de sorte à obtenir une image de fluctuations d’absorption représentative des fluctuations d’absorption dues à l’échantillon. Procédé selon l’une quelconque des revendications 1 à 4, dans lequel le filtrage (120) spatio-temporel multispectral par décomposition en valeurs singulières comprend une sélection des composantes correspondant aux valeurs singulières de plus haute énergie à supprimer et une suppression de composantes sélectionnées, la sélection étant effectuée en choisissant parmi un ensemble de valeurs d’index un index identifiant la première composante à conserver pour lequel un rapport contraste à bruit est maximal, le rapport contraste à bruit déterminé pour un index étant déterminé pour des images de variance filtrées calculées par filtrage (120) spatio-temporel multispectral par décomposition en valeurs singulières appliquant cet index pour identifier la première composante à conserver puis normalisation par rapport à une fluence moyenne. Procédé selon la revendication 5 dans lequel le rapport contraste à bruit est déterminé en éliminant le contraste dû à la valeur moyenne des images acquises pour au moins une longueur d’onde . Procédé selon l’une quelconque des revendications 4 à 6 lorsqu’elle dépend de la revendication 4, le procédé comprenant un calcul (160) d’une image du taux de saturation en oxygène à partir d’au moins deux images de fluctuations d’absorption obtenues pour au moins deux longueurs d’onde correspondantes, sur la base d’un modèle exprimant pour chaque pixel de coordonnée r une relation entre une concentration totale en hémoglobine, un taux de saturation en oxygène et la valeur au pixel r de l’image de fluctuations d’absorption. Dispositif de traitement d’images photoacoustiques comprenant, au moins une mémoire de données comprenant des instructions de code de programme, au moins un processeur de données, le processeur de données étant configurées pour, lors que les instructions de code de programme sont exécutées par le processeur de données, causer l’exécution par le dispositif de traitement d’images photoacoustiques d’un procédé selon l’une quelconque des revendications 1 à 7. Support de données lisible par ordinateur comprenant des instructions de programme informatique qui, lorsqu’elles sont exécutées par un processeur, causent l’exécution d’un procédé selon l’une quelconque des revendications 1 à 7. Programme informatique comprenant des instructions de programme informatique qui, lorsqu’elles sont exécutées par un processeur, causent l’exécution d’un procédé selon l’une quelconque des revendications 1 à 7.