Procédé de tatouage numérique (10) d’une image (1) par un filigrane (2), le procédé mis en œuvre par un système de calcul dédiée ou un ordinateur comprenant les étapes suivantes : Fournir (101) l’image I m à tatouer, l’image I m présentant une taille N×M, Appliquer (102), à l’image I m , une transformation fréquentielle FT pour obtenir (103) une matrice FT _I m d’éléments FT_Im (i,j) et de même taille que l’image I m fournie, Fournir (104) un code utilisateur Wi , mettre en œuvre (105) un générateur de nombres (pseudo)aléatoires basé sur un algorithme déterministe en l’initialisant avec le code utilisateur Wi , de sorte à générer (106) une suite de nombres (pseudo)aléatoires AWMi(j) avec j = 1, …, Li, où la taille Li de la suite est strictement inférieure à la taille N×M de l’image I m , à partir de la suite de nombres (pseudo)aléatoires AWMi(j), générer (107) le filigrane WMi à tatouer sur l’image I m , le filigrane WMi présentant la taille Li de la suite de nombres (pseudo)aléatoires AWMi(j), intégrer (108) le filigrane WMi dans la matrice FT _I m : en définissant (1081), dans la matrice FT _I m , une zone d’intégration EA d’éléments EA(i,j) et de taille P×Q inférieure à la taille N×M de la matrice F T _I m et à partir d’un élément FT_Im(ni, mi) déterminé de la matrice FT_Im tel que (P-mi+1) en convertissant (1082) la zone d’intégration EA en un vecteur LV de longueur K égale ou inférieure à la taille P×Q de la zone d’intégration EA , la conversion consistant à : attribuer, à un premier élément LV(1)= EA(1,1) du vecteur LV , la valeur de l’élément FT_Im(ni, mi) déterminé de la matrice FT _I m , puis attribuer respectivement, à chacun des éléments suivants du vecteur LV , la valeur de chacun des éléments FT_Im(i,j) suivants l’élément FT_Im(ni, mi) déterminé de la matrice F T _I m en parcourant la zone d’intégration EA selon un schéma dit d’extraction prédéterminé, en combinant (1083), avec le filigrane WMi , une partie du vecteur LV qui est comprise entre un élément de départ LV(kStart) et un élément d’arrivée LV(kEnd), avec kStart>1, kEnd L+2(kStart-1), de sorte à générer (1084) un vecteur WMV de taille K composé des kStart-1 premiers éléments du vecteur LV , d’éléments WMV(kStart) à WMV(kEnd) issus de la combinaison des éléments de LV(kStart) à LV(kEnd) avec les éléments de WMi(1) à WMi(Li) du filigrane WMi et des derniers K-kEnd éléments du vecteur LV , puis en introduisant (1085) le vecteur WMV dans la matrice F T _I m en lieu et place de la zone d’intégration EA , en parcourant la zone d’intégration EA selon un schéma dit d’insertion identique audit schéma d’extraction prédéterminé, pour générer (1086) une matrice F T _ WM I m qui représente, dans le domaine fréquentiel, l’image I m combinée au filigrane WMi , appliquer (109), à la matrice F T _ WM I m , une transformation inverse IT à la transformation fréquentielle FT, pour obtenir une image tatouée WMIm . Procédé de tatouage numérique (10) selon la revendication précédente, dans lequel : l’image I m étant une image en noir et blanc, la transformation fréquentielle FT est appliquée à l’image I m , ou l’image I m étant une image en couleurs, la transformation fréquentielle FT est appliquée : soit à une luminance Y de l’image I m en une représentation luminance-chrominance ( Y, Cb, Cr ), soit à au moins une des trois composantes rouge, vert et bleu de l’image I m en une représentation couleur rouge-vert-bleu. Procédé de tatouage numérique (10) selon l’une quelconque des revendications précédentes, dans lequel la génération (106) du filigrane WMi à tatouer sur l’image I m comprend un filtrage (1061) de la suite de nombres (pseudo)aléatoires AWMi(j) par un filtre passe-bas. Procédé de tatouage numérique (10) selon l’une quelconque des revendications précédentes, dans lequel la génération (106) du filigrane WMi à tatouer sur l’image I m comprend une multiplication (1062) par une fonction de pondération F(j) de chaque nombre de la suite de nombres (pseudo)aléatoires AWMi(j). Procédé de tatouage numérique (10) selon la revendication précédente, dans lequel la fonction de pondération F(j) est donnée par l’équation suivante : avec j=1, …, Li, où Li est la taille de la suite de nombres (pseudo)aléatoires AWMi(j). Procédé de tatouage numérique (10) selon l’une quelconque des revendications précédentes, dans lequel le schéma d’extraction est défini par un parcours en zig-zag (30) des éléments de la matrice F T _I m compris dans la zone d’intégration EA, en partant de l’élément FT_Im(ni, mi) et en finissant par l’élément FT_Im(ni+P-1, mi+Q-1). Procédé de tatouage numérique (10) selon l’une quelconque des revendications précédentes, dans lequel l’intégration (108) du filigrane WMi dans la matrice FT_I m est telle que les éléments WMV(kStart) à WMV(kEnd) du vecteur WMV sont donnés par l’équation suivante : WMV(k)=LV(k)+α*DLV*WMi(j), avec j=1, …, Li, et k = kStart, … kEnd, où DLV est l’écart-type des éléments du vecteur LV(k), WMi(j) est l’élément j du filigrane WMi , et où α est un coefficient de pondération prenant une valeur prédéterminée. Procédé de tatouage numérique (10) selon la revendication précédente, dans lequel la valeur du coefficient de pondération α est choisie en fonction d’un compromis entre visibilité du filigrane WMi dans l’image I m tatouée et fiabilité de détection du filigrane WMi , par exemple la valeur du coefficient de pondération α est sensiblement comprise entre 0,01 et 2,5, de préférence sensiblement comprise entre 0,05 et 0,5, et par exemple sensiblement égale à 0,1. Procédé de tatouage numérique (10) selon l’une quelconque des revendications précédentes, dans lequel l’indice kStart est supérieur à 10, de préférence supérieure à 50 et encore plus préférentiellement sensiblement égale à 100, et/ou la taille Li de la suite de nombres (pseudo)aléatoires AWMi(j) est supérieure à 500, de préférence égale ou supérieure à 1000. Procédé de tatouage numérique (10) selon l’une quelconque des revendications précédentes, dans lequel la même image I m est tatouée avec une pluralité de filigranes WMi , WMj , etc. Procédé de tatouage numérique (10) selon la revendication précédente, dans lequel chaque filigrane WMi étant associé à une zone d’intégration EA i qui lui est propre, la zone d’intégration EA i d’au moins un filigrane WMi recouvre au moins partiellement la zone d’intégration EA j d’au moins un autre filigrane WMj . Procédé de détection (20) d’au moins un filigrane (2) dans une image (3), le procédé étant mis en œuvre par un système de calcul dédié ou un ordinateur mémorisant, en association avec chaque filigrane WMi à détecter, un générateur de nombres (pseudo)aléatoires basé sur un algorithme déterministe, un schéma dit d’extraction prédéterminé et des données d’identification de chaque filigrane WMi à détecter, ces dernières comprenant, voire consistant en, un code utilisateur Wi , une taille Li du filigrane WMi et des coordonnées (ni, mi) d’un élément WMIm(ni, mi) de l’image WMIm à partir duquel le filigrane WMi est à rechercher, le procédé comprenant les étapes suivantes : Fournir (201) l’image WM I m dans laquelle le filigrane WMi est à détecter, l’image WMI m présentant une taille N×M strictement supérieure à la taille Li du filigrane WMi à détecter, Appliquer (202), à l’image WM I m , une transformation fréquentielle FT pour obtenir (203) la matrice FT_ WM I m d’éléments FT_WMIm(i,j) et de même taille que l’image WM I m fournie, en fonction de la taille Li du filigrane WMi et du code utilisateur Wi fourni, mettre en œuvre (205) le générateur de nombres (pseudo)aléatoires en l’initialisant (204) avec le code utilisateur Wi , pour générer (206) une suite de nombres (pseudo)aléatoires AWMi(j) présentant la taille Li, à partir de la suite de nombres (pseudo)aléatoires AWMi(j), générer (206) le filigrane WMi , ce dernier présentant la taille Li de la suite de nombres (pseudo)aléatoires AWMi(j), à partir de la taille Li du filigrane WMi , des coordonnées (ni, mi) de l’élément WMIm(ni, mi) de l’image WMIm à partir duquel le filigrane WMi est à rechercher et du schéma d’extraction prédéterminé, extraire (207), de la matrice FT_WMI m , un vecteur LVm de longueur K=P*Q strictement supérieure à la taille Li du filigrane WMi et strictement inférieure à la taille N×M de l’image WMIm , l’extraction consistant à : attribuer, à un premier élément LVm(1) du vecteur LVm , la valeur de l’élément FT_WMIm(ni, mi) de la matrice FT_WMI m , puis attribuer respectivement, à chacun des éléments suivants du vecteur LVm , la valeur de chacun des éléments FT_WMIm(i,j) suivants l’élément FT_WMIm(ni, mi) déterminé de la matrice FT_WMI m en la parcourant en partie selon le schéma dit d’extraction prédéterminé, puis en fonction du filigrane WMi généré et d’une partie du vecteur LVm extrait qui est comprise entre un élément de départ LVm(kStart) et un élément d’arrivée LVm(kEnd), où kEnd-kStart=Li -1 et K>Li+2(kStart-1), calculer (208) les valeurs Cxy(k) d’une fonction de corrélation croisée Cxy entre le filigrane WMi généré et le vecteur LVm extrait, pour une pluralité de valeurs k de l’indice de départ kStart comprises entre 1 et K-Li+1, et vérifier (209) que les valeurs Cxy(k) calculées de la fonction de corrélation croisée Cxy comprennent une valeur maximale qui est représentative d’une corrélation effective entre le filigrane WMi généré et le vecteur LVm extrait, de sorte à détecter, en cas de vérification positive, le filigrane WMi dans l’image WMIm . Procédé de détection (20) selon la revendication précédente, dans lequel les valeurs Cxy(k) de la fonction de corrélation croisée Cxy sont définies par l’équation suivante : avec et Procédé de détection (20) selon l’une quelconque des revendications 12 et 13, dans lequel, le système de calcul dédié ou l’ordinateur mémorisant en outre une valeur de l’indice kStart en association avec le filigrane WMi à détecter, les valeurs Cxy(k) de la fonction de corrélation croisée Cxy entre le filigrane WMi généré et le vecteur LVm extrait sont calculées (208) pour les seules valeurs de l’indice k qui sont comprises entre kStart-d et kStart+d, où d est un entier prédéterminé, strictement inférieur à kStart et compris entre 1 et 20, de préférence entre 2 et 10, et encore plus préférentiellement égal à 4. Procédé de détection (20) selon l’une quelconque des revendications 12 à 14, dans lequel l’image WM I m a été tatouée en mettant en œuvre le procédé de tatouage numérique (10) selon l’une quelconque des revendications 1 à 11. Procédé de détection (20) selon la revendication précédente, dans lequel l’image WMI m telle que tatouée a été modifiée (1020), suite à son tatouage, par filtrage et/ou ajout de bruit et/ou modification de l’intensité des pixels qui la composent, avant d’être fournie (201). Procédé d’authentification et/ou d’identification d’une image I m et/ou de l’auteur, du propriétaire ou de l’utilisateur autorisé d’une image I m mettant en œuvre le procédé de tatouage numérique (10) selon l’une quelconque des revendications 1 à 11 et le procédé de détection (20) selon l’une quelconque des revendications 12 à 16. Produit programme d’ordinateur comprenant des instructions qui, lorsqu’elles sont implémentées par un système de calcul dédié ou un ordinateur, exécute les étapes d’au moins l’un parmi le procédé de tatouage numérique (10) selon l’une quelconque des revendications 1 à 11 et le procédé de détection (20) selon l’une quelconque des revendications 12 à 16.