Procédé de démosaïquage appliqué à une image initiale présentant des pixels correspondants chacun à une bande spectrale initiale pour l’obtention d’une image démosaïquée à plusieurs bandes spectrales présentant des pixels correspondants chacun à ladite bande spectrale initiale ou à une bande spectrale interpolée, caractérisé en ce qu’il comprend : une transformation de stabilisation de la variance du bruit, appliquée aux pixels des bandes spectrales initiales, sur la base d’une pluralité de patchs de référence prédéterminés et pour chaque patch de référence parmi la pluralité de patchs de référence : une détermination (310) d’un ensemble de patchs similaires audit patch de référence, en exploitant les pixels d’au moins une des bandes spectrales initiales transformées une détermination (320) d’une matrice de covariance dudit ensemble de patchs similaires et pour chaque patch similaire dudit ensemble de patchs similaires, un calcul d’estimatifs des bandes spectrales interpolées, fonction d’une décomposition en bloc de la matrice de covariance (340), où la décomposition en quatre blocs est fonction du phasage de chaque patch similaire, une détermination des bandes spectrales interpolées par une agrégation (350) des estimatifs des bandes spectrales interpolées, où on agrège les estimatifs des bandes spectrales interpolées pour chaque pixel et une transformation inverse de la transformation de stabilisation, appliquée aux pixels des bandes spectrales initiales transformées. Procédé de démosaïquage selon la revendication 1, dans lequel la transformation de stabilisation de la variance du bruit est réalisée en fonction d’un modèle de bruit propre au détecteur. Procédé selon la revendication 1 ou 2, dans lequel les patches de référence prédéterminés sont répartis de façon à couvrir une proportion déterminée de l’image. Procédé selon l’une des revendications 1 à 3, dans lequel lors de la détermination de la matrice de covariance de l’ensemble des patchs similaires, les patchs similaires comprennent des bandes spectrales initiales transformées et des bandes spectrales intermédiaires, ces bandes spectrales intermédiaires étant issues d’un démosaïquage préalable par interpolation linéaire appliqué à l’image initiale stabilisée. Procédé selon la revendication précédente, dans lequel la détermination (310) dudit ensemble de patchs similaires à chaque patch de référence est déterminée par une fonction de similitude appliquée à une image intermédiaire résultant du démosaïquage préalable par interpolation linéaire de l’image initiale stabilisée, la fonction de similitude utilisant une ou plusieurs des bandes spectrales de l’image intermédiaire. Procédé selon l’une des revendications 1 à 3, dans lequel lors de la détermination de la matrice de covariance de l’ensemble des patchs similaires, les patchs similaires comprennent des bandes spectrales initiales transformées et des bandes spectrales intermédiaires, ces bandes spectrales intermédiaires étant issues d’un démosaïquage préalable par interpolation linéaire et d’un débruitage préalable appliqués à l’image initiale stabilisée. Procédé selon la revendication précédente, dans lequel la détermination (310) dudit ensemble de patchs similaires à chaque patch de référence est déterminée par une fonction de similitude appliquée à une image intermédiaire résultant du démosaïquage préalable par interpolation linéaire (100) et du débruitage préalable de l’image initiale stabilisée, la fonction de similitude utilisant une ou plusieurs des bandes spectrales de l’image intermédiaire. Procédé selon la revendication 6 ou 7, dans lequel le débruitage consiste en une fonction de débruitage local appliquée à chaque patch et fonction dudit ensemble de patches similaires au patch considéré. Procédé selon l’une des revendications 1 à 8, dans lequel les bandes spectrales correspondent à l’espace colorimétrique RGB ou YUV. Procédé selon l’une des revendications 1 à 9, dans lequel les bandes spectrales correspondent à l’espace colorimétrique RGB et les pixels de l’image initiale stabilisée sont organisés selon un motif de Bayer. Procédé selon l’une des revendications précédentes, dans lequel l’agrégation (350) des estimatifs des bandes spectrales interpolées est réalisée par un calcul d’une moyenne de ces estimatifs de bandes spectrales interpolées. Procédé selon l’une des revendications précédentes, dans lequel la décomposition en quatre blocs fonction du phasage dudit patch similaire, est fonction du couple (L,C) correspondant au L ième pixel et au C ième pixel dudit patch similaire, chaque valeur de la matrice à la ligne L et à la colonne C, étant, selon qu’il s’agisse de : - un pixel L correspondant à une bande spectrale initiale transformée et un pixel C correspondant à une bande spectrale initiale transformée, attribuée à un premier bloc - un pixel L ne correspondant pas à une bande spectrale initiale et un pixel C correspondant à une bande spectrale initiale transformée, attribuée à un deuxième bloc - un pixel L correspondant à une bande spectrale initiale transformée et un pixel C ne correspondant pas à une bande spectrale initiale, attribuée à un troisième bloc ou - un pixel L ne correspondant pas à une bande spectrale initiale et un pixel C ne correspondant pas à une bande spectrale initiale, attribuée à un quatrième bloc . Procédé selon la revendication précédente, dans lequel, pour chaque patch similaire dans ledit ensemble de patchs similaires, le calcul des estimatifs des bandes spectrales interpolées est réalisé comme suit : où : est le vecteur des estimatifs des bandes spectrales interpolées dudit patch similaire, est le vecteur des valeurs moyennes des pixels des patches dudit ensemble de patchs similaires, où on considère les pixels ne correspondant pas à une bande spectrale initiale, est le vecteur des pixels correspondant à une bande spectrale initiale transformée dans ledit patch similaire, est le vecteur des valeurs moyennes des pixels des patches dudit ensemble de patchs similaires, où on considère les pixels correspondant à une bande spectrale initiale transformée, est la matrice inverse du premier bloc et est le deuxième bloc. Procédé selon la revendication précédente, dans lequel lors du calcul des estimatifs des bandes spectrales interpolées pour des patchs similaires présentant un même phasage, intervient une même décomposition de la matrice de covariance en blocs qui sont mémorisés et réutilisés pour tous ces patchs similaires présentant le même phasage. Produit programme d’ordinateur, comprenant des instructions de code pour la mise en œuvre du procédé selon l’une des revendications précédentes, lorsqu’il est exécuté par un calculateur. Dispositif de traitement d’images comprenant un calculateur (11), caractérisé en ce qu’il est configuré pour mettre en œuvre le procédé selon l’une des revendications 1 à 14.