Procédé et dispositif de codage d’un ensemble de coefficients, procédé et dispositif de décodage d’un ensemble de coefficients, flux de données et programme d’ordinateur associés Domaine technique de l'invention La présente invention concerne le domaine technique de l’holographie numérique. Elle concerne plus particulièrement un procédé et un dispositif de codage d’un ensemble de coefficients participant à une représentation d’un hologramme numérique, ainsi qu’un procédé et un dispositif de décodage d’un ensemble de coefficients participant à une représentation d’un hologramme numérique, un flux de données et un programme d’ordinateur associés. Etat de la technique Un hologramme numérique est l’enregistrement sous forme numérique d’une section, au niveau d’un plan de référence, d’un champ lumineux se propageant dans un espace tridimensionnel, en vue de pouvoir restituer ultérieurement ce champ lumineux à un utilisateur. Les propriétés physiques des champs lumineux ont conduit à représenter les hologrammes numériques au moyen de coefficients respectivement associés à des éléments d’un domaine spatio-fréquentiel défini par deux dimensions spatiales (qui correspondent au plan de référence) et par les deux dimensions fréquentielles associées, c’est-à-dire à des multiplets formés chacun par un premier paramètre spatial, un second paramètre spatial, un premier paramètre fréquentiel et un second paramètre fréquentiel. Présentation de l'invention L’invention propose un procédé de codage d’un ensemble de coefficients participant à une représentation d’un hologramme numérique en étant respectivement associés à des multiplets formés chacun par un premier paramètre spatial, un second paramètre spatial, un premier paramètre fréquentiel et un second paramètre fréquentiel, le procédé comprenant les étapes suivantes, effectuées pour au moins une paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel : - codage de données indicatives de positions associées respectivement à des coefficients non-nuls parmi un ensemble des positions associées de manière prédéfinie aux coefficients associés respectivement aux différents multiplets définis par ladite paire de valeurs du premier paramètre fréquentiel et du second paramètre fréquentiel ; - codage des coefficients associés respectivement aux positions indiquées par lesdites données indicatives. Le regroupement des coefficients relatifs à des multiplets définis par la même paire de valeurs des paramètres fréquentiels en vue de leur codage permet de coder des informations mieux localisées, en général au niveau de certains coefficients seulement, et ainsi d’augmenter l’efficacité du codage. Selon une possibilité, au moins une donnée indicative parmi lesdites données indicatives peut représenter une différence entre une position associée à un coefficient non-nul et une autre position associée à un autre coefficient non-nul. Selon une autre possibilité, lesdites données indicatives peuvent désigner au moins une position associée à un coefficient non-nul sans désigner au moins une autre position associée à un autre coefficient non-nul et indiquée par des données indicatives codées pour une autre paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel. On code ainsi la position des coefficients non-nuls relativement à celles rencontrées pour une autre paire de valeurs, ce qui augmente l’efficacité du codage. On peut prévoir que les coefficients associés aux différents multiplets définis par ladite paire de valeurs du premier paramètre fréquentiel et du second paramètre fréquentiel soient ordonnés de manière prédéfinie en une séquence ordonnée de coefficients. La donnée indicative d’une position associée à un coefficient non-nul peut alors indiquer la position de ce coefficient non-nul dans ladite séquence ordonnée de coefficients, ce qui permet de limiter la quantité de données nécessaires au codage des positions. On peut prévoir par ailleurs qu’à l’étape de codage des coefficients, les coefficients non-nuls soient codés dans l’ordre défini par la séquence ordonnée. L’invention propose également un procédé de décodage d’un ensemble de coefficients participant à une représentation d’un hologramme numérique en étant respectivement associés à des multiplets formés chacun par un premier paramètre spatial, un second paramètre spatial, un premier paramètre fréquentiel et un second paramètre fréquentiel, le procédé comprenant les étapes suivantes, effectuées pour au moins une paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel : - décodage de données indicatives de positions associées respectivement à des coefficients non-nuls parmi un ensemble des positions associées de manière prédéfinie aux coefficients associés respectivement aux différents multiplets comprenant ladite paire de valeurs du premier paramètre fréquentiel et du second paramètre fréquentiel ; - décodage des coefficients associés respectivement aux positions indiquées par lesdites données indicatives. Comme déjà indiqué pour le procédé de codage, on peut prévoir qu’au moins une donnée indicative parmi lesdites données indicatives représente une différence entre une position associée à un coefficient non-nul et une autre position associée à un autre coefficient non-nul. Selon une autre possibilité, lesdites données indicatives peuvent désigner au moins une position associée à un coefficient non-nul sans désigner au moins une autre position associée à un autre coefficient non-nul et indiquée par des données indicatives décodées pour une autre paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel. L’invention propose par ailleurs un dispositif de codage d’un ensemble de coefficients participant à une représentation d’un hologramme numérique en étant respectivement associés à des multiplets formés chacun par un premier paramètre spatial, un second paramètre spatial, un premier paramètre fréquentiel et un second paramètre fréquentiel, comprenant : - un premier module de codage conçu pour coder, pour au moins une paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel, des données indicatives de positions associées respectivement à des coefficients non-nuls parmi un ensemble des positions associées de manière prédéfinie aux coefficients associés respectivement aux différents multiplets définis par ladite paire de valeurs du premier paramètre fréquentiel et du second paramètre fréquentiel ; - un second module de codage conçu pour coder les coefficients associés respectivement aux positions indiquées par lesdites données indicatives. L’invention propose en outre un dispositif de décodage d’un un ensemble de coefficients participant à une représentation d’un hologramme numérique en étant respectivement associés à des multiplets formés chacun par un premier paramètre spatial, un second paramètre spatial, un premier paramètre fréquentiel et un second paramètre fréquentiel, comprenant : - un premier module de décodage conçu pour décoder, pour au moins une paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel, des données indicatives de positions associées respectivement à des coefficients non-nul parmi un ensemble des positions associées de manière prédéfinie aux coefficients associés respectivement aux différents multiplets comprenant ladite paire de valeurs du premier paramètre fréquentiel et du second paramètre fréquentiel ; - un second module de décodage conçu pour décoder les coefficients associés respectivement aux positions indiquées par lesdites données indicatives. L’invention propose aussi un flux de données représentant un ensemble de coefficients participant à une représentation d’un hologramme numérique en étant respectivement associés à des multiplets formés chacun par un premier paramètre spatial, un second paramètre spatial, un premier paramètre fréquentiel et un second paramètre fréquentiel, le flux de données comprenant, pour au moins une paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel : - des données indicatives de positions associées respectivement à des coefficients non-nuls parmi un ensemble des positions associées de manière prédéfinie aux coefficients associés respectivement aux différents multiplets définis par ladite paire de valeurs du premier paramètre fréquentiel et du second paramètre fréquentiel ; - des données représentatives des coefficients associés respectivement aux positions indiquées par lesdites données indicatives. L’invention propose par ailleurs un programme d’ordinateur comprenant des instructions conçues pour mettre en œuvre un procédé tel que proposé ci-dessus lorsque ces instructions sont exécutées par un processeur. Les procédés proposés ci-dessus peuvent ainsi être mis en œuvre par ce processeur, par exemple lorsque ce processeur équipe un ordinateur ou un dispositif de codage ou de décodage tel que mentionné ci-dessus et décrit plus loin. L’invention propose enfin un support d’informations (par exemple non-transitoire) lisible par un ordinateur, sur lequel est enregistré un tel programme d'ordinateur. Bien entendu, les différentes caractéristiques, variantes et formes de réalisation de l'invention peuvent être associées les unes avec les autres selon diverses combinaisons dans la mesure où elles ne sont pas incompatibles ou exclusives les unes des autres. Description détaillée de l'invention De plus, diverses autres caractéristiques de l'invention ressortent de la description annexée effectuée en référence aux dessins qui illustrent des formes, non limitatives, de réalisation de l'invention et où : - la figure 1 représente schématiquement un hologramme numérique ; - la figure 2 représente un exemple de dispositif d’holographie numérique ; - la figure 3 représente un exemple de procédé mis en œuvre par un module de traitement du dispositif d’holographie numérique de la figure 2 ; - la figure 4 représente un procédé de décodage d’un hologramme numérique à partir de données reçues dans un flux de données ; - la figure 5 représente un exemple de dispositif de codage conforme à l’invention ; - la figure 6 représente un exemple de dispositif de décodage conforme à l’invention ; - la figure 7 représente un procédé d’échange de données entre le dispositif de codage de la figure 5 et le dispositif de décodage de la figure 6 ; - la figure 8 représente une variante de réalisation du procédé de la figure 7 ; - la figure 9 représente un procédé de partitionnement du domaine spatio- fréquentiel D ; - la figure 10 représente un procédé de codage de coefficients associés respectivement à des éléments du domaine spatio-fréquentiel D ; - la figure 11 représente un procédé de décodage correspondant au procédé de codage de la figure 10 ; - la figure 12 représente un exemple de dispositif de codage conforme aux enseignements de l’invention ; et - la figure 13 représente un exemple de dispositif de décodage conforme aux enseignements de l’invention. La figure 1 représente schématiquement un hologramme numérique défini au moyen d’une fonction ψ0 dans un plan de référence Pref d’équation (z=0) dans un espace de propagation E où les coordonnées tridimensionnelles d’un point sont définies dans un repère orthonormé (O,x,y,z). Autrement dit, les axes (Ox) et (Oy) (ou plus simplement les axes x et y) sont compris dans le plan de référence Pref et/ou définissent le plan de référence Pref. La direction de propagation (avant application de la transformation mentionnée plus bas) correspond ainsi ici à l’axe (Oz) ; dans cette situation, la direction de propagation (z) est perpendiculaire (ou orthogonale) au plan de référence Pref. Le champ d’onde (représentant ici le champ lumineux monochromatique de longueur d’onde λ) associé à cet hologramme en un point de l’espace de propagation E de coordonnées r=(x,y,z) s’écrit alors (avant application de la transformation mentionnée plus bas) : [Math.1] où la transformée de Fourier de la fonction ψ ₀ est définie (dans le plan de référence Pref) par [Math.2] λ est la longueur d’onde de la lumière (monochromatique) concernée, η = (ηx, η _y , η _z ) avec η _z = SQRT(1 - η _x ² - η _y ² ), SQRT() étant l’opérateur racine carrée, <r,η> est le produit scalaire de r et η, et D _λ est le domaine défini par {η _x ² + η _y ² < 1}. Le champ d’onde défini par l’hologramme numérique peut ainsi être vu comme une somme (en théorie continue ou infinie) d’ondes planes de longueur d’onde λ. Afin de tenir compte des particularités des hologrammes, on utilise des représentations de ces hologrammes sous forme de distribution dans un domaine spatio-fréquentiel D (ou espace des phases) à 4 dimensions, dont les 2 premières dimensions correspondent aux coordonnées spatiales dans le plan P de l’hologramme et les 2 autres dimensions correspondent aux fréquences spatiales dans ce même plan. Une telle distribution comprend par exemple les valeurs d’énergie lumineuse respectivement associées aux coordonnées spatiales et aux fréquences spatiales du domaine spatio-fréquentiel D. Selon un premier exemple, on peut appliquer une transformée de Fourier fenêtrée (en anglais : "Windowed Fourier Transform") au champ d’onde ψ ; la distribution Wψ ^Φ s’écrit alors : [Math.3] où Φ est une fonction dite fenêtre, ψ* est le conjugué (complexe) de ψ, X représente les coordonnées spatiales (x, y) et Ξ les fréquences spatiales (ξ _x , ξ _y ). Selon un second exemple, on peut utiliser une distribution de Wigner-Ville Wψ obtenue en appliquant une transformée de Wigner-Ville au champ d’onde ψ : [Math.4] avec comme précédemment X représentant les coordonnées spatiales (x, y) et Ξ les fréquences spatiales (ξ _x , ξ _y ). Il s’agit donc d’une transformée de Fourier fenêtrée particulière, dans laquelle la fonction fenêtre utilisée est la fonction ψ. Selon un troisième exemple, on peut utiliser une décomposition en trames de Gabor (en anglais : "Gabor frames"), qui correspond à une discrétisation de la transformée de Fourier fenêtrée. Pour ce faire, on définit un maillage Λ du domaine spatio-fréquentiel D comme suit : [Math.5] où A est une matrice inversible de dimensions 4x4. La décomposition en trames de Gabor peut alors s’écrire : [Math.6] ψ = ^ ^^ ^{^} ^ (^)^^(^)^ ^∈^ où Φ ^D est la fenêtre duale de la fonction Φ, π(ζ) est la translation dans le domaine spatio-fréquentiel D de vecteur ζ et avec l’opérateur de modulation/translation suivant : [Math.7] où X, Ξ sont respectivement la composante spatiale (à deux dimensions) et la composante fréquentielle (à deux dimensions) du vecteur ζ, soit : ζ = (X , Ξ), et où < , > représente le produit scalaire (ici dans le plan). On étudie ici le champ ψ _T obtenu lorsqu’on applique une transformation T de l’espace tridimensionnel E à l’hologramme numérique. D’après Kyoji Matsushima, Hagen Schimmel & Frank Wyrowski, "Fast calculation method for optical diffraction on tilted planes by use of the angular spectrum of plane waves", in J. Opt. Soc. Am. A 20, 1755-1762 (2003), le champ d’onde (ou champ lumineux) modifié par la transformation T s’écrit comme suit (au niveau d’un plan parallèle au plan de référence P _ref et de coordonnée z) : [Math.8] Cette expression implique notamment le produit scalaire du vecteur de coordonnées (x y z) ^t de l’espace de propagation tridimensionnel E et du vecteur Tη résultant de l’application de la transformation T au vecteur d’onde (η _x , η _y , η _z ) de norme prédéfinie C et dont les deux premières coordonnées sont les variables fréquentielles (ηx, ηy) respectivement associées aux coordonnées spatiales (x,y) dans le plan de référence Pref de l’hologramme numérique. La norme prédéfinie C dépend de la représentation utilisée et vaut ici 1 (représentation normalisée). En variante, la norme prédéfinie pourrait être égale à la longueur d’onde λ. En écrivant : [Math.9] ^(^, ^) =< ^, ^^ >, l’expression de la fonction ψ _z ^T donnée ci-dessus peut être vue comme un opérateur intégral de Fourier utilisé en analyse semi-classique dans le domaine de la mécanique quantique, où la fonction φ définie ci-dessus est la fonction génératrice (en anglais "generating function") associée à cet opérateur. En s’intéressant à un plan particulier Pz parallèle au plan de référence, de coordonnée z selon la direction perpendiculaire au plan de référence Pref, on peut écrire la fonction génératrice comme une fonction à quatre variables (c’est-à-dire une fonction à valeurs sur ℝ ⁴ ) comme suit : [Math.10] La fonction génératrice φz est donc le produit scalaire : - d’un premier vecteur (de coordonnées tridimensionnelles), dont les deux premières coordonnées correspondent aux deux premières variables des quatre variables et dont la troisième coordonnée correspond à la coordonnée z du plan Pz concerné, et - d’un vecteur résultant de l’application de la transformation T à un second vecteur, dont les deux premières coordonnées correspondent aux deux dernières variables des quatre variables et, ici, dont la norme C est une constante prédéfinie (cette norme C valant ici 1 du fait de la normalisation déjà mentionnée ; en variante, dans d’autres représentations, cette norme C pourrait par exemple être égale à la longueur d’onde λ). Les explications sont données dans la suite pour un plan P _z de coordonnée z donné (plan parallèle au plan de référence P _ref ), mais sont valables pour tout plan parallèle au plan de référence Pref (i.e. pour z quelconque). Par souci de simplification, on note dans la suite φ la fonction génératrice φ _z associée à ce plan P _z . Dans le contexte défini ci-dessus, si on note χ le symplectomorphisme du domaine spatio-fréquentiel D vers lui-même qui correspond à la transformation T, le sous-ensemble (ou graphe) Γ de ℝ ⁴ x ℝ ⁴ comprenant les éléments de la forme (e, χ(e)), pour e variant dans ℝ ⁴ , peut être également défini comme l’ensemble des éléments de la forme : [Math.11] pour (x,y,η _u , η _v ) variant dans ℝ ⁴ . On pourra se référer à l’ouvrage "Semi-classical Analysis", de Victor Guillemin et Shlomo Sternberg, International Press, 2013, ISBN 1571462767, 9781571462763 pour plus d’explications à ce sujet. On rappelle qu’un symplectomorphisme est un difféomorphisme qui conserve la forme symplectique, c’est-à-dire que dont la matrice jacobienne M vérifie : J=M ^t JM, avec : [Math.12] En notant [Math.13] on peut donner l’expression des variables η _u et η _v en fonction de e = (x,y,ξ _x ,ξ _y ) : [Math.14] Les 5 ^ème et 6 ^ème coordonnées (u, v) des éléments du sous-ensemble Γ peuvent donc quant à elles s’écrire : [Math.15] On obtient ainsi une expression paramétrique du sous-ensemble Γ en fonction de ses quatre premières coordonnées (x,y,ξ _x ,ξ _y ), ce qui permet d’obtenir l’expression suivante du symplectomorphisme χ (étant rappelé que le sous-ensemble Γ correspond aux éléments de la forme (e, χ(e)) pour e variant dans ℝ ⁴ ) : [Math.16] En considérant à présent que la transformation T étudiée est une transformation rigide de l’espace de propagation E, cette transformation T est donc la combinaison d’une translation de vecteur (t _x ,t _y ,t _z ) et d’une rotation R qui peut-être décomposée en 3 rotations à l’aide des angles d’Euler : R = R ^z x z αRθRβ, où R ^z _α et Rz _β sont des rotations autour de l’axe z d’angles respectives α et β, et R ^x θ est une rotation autour de l’axe x d’angle θ. On peut ainsi définir des symplectomorphismes χt, χα, χθ, χβ respectivement associés à la translation et à ces 3 rotations R ^z x z α, Rθ, Rβ, et dont la composition est le symplectomorphisme χ associé à la transformation T : χ = χβ o χθ o χα o χt, où o est l’opérateur de composition. Les symplectomorphismes χ _t , χ _α , χ _β associés à la translation et aux rotations Rz _α , R ^x _θ , Rz _β autour de l’axe Z s’écrivent simplement : χt(x,y,ξx,ξy) = (x+tx,y+ty,ξx,ξy) χα(x,y,ξx,ξy) = (x.cos α–y.sin α,x.sin α+y.cos α, ξx.cos α–ξy.sin α, ξx.sin α+ξy.cos α) χ _β (x,y,ξ _x ,ξ _y ) = (x.cos β –y.sin β,x.sin β +y.cos β, ξ _x .cos β –ξ _y .sin β, ξ _x .sin β +ξ _y .cos β) On note que la troisième composante t _z de la translation est quant à elle prise en compte dans la propagation, c’est-à-dire en considérant le plan Pz+tz. On propose d’utiliser les développements exposés ci-dessus pour déterminer le symplectomorphisme χ x θ associé à la rotation Rθ d’angle θ autour de l’axe x. Dans ce cas, la fonction génératrice φ s’écrit : [Math.17] soit : [Math.18] ^ ⁽ ^, ^ ⁾ = ^. ^ _^ + ^. (^ _^ . cos ^ − ^ _^ . sin ^) + ^. (^ _^ . sin ^ + ^ _^ . cos ^) En utilisant le fait que : [Math.19] (voir ci-dessus la définition du sous-ensemble Γ) et en remplaçant ηz par son expression en fonction de η _x et de η _y donnée plus haut, on obtient : ξx = ηx (ηy.cos θ) – SQRT(1 – ηx ² – ηy ² ).sin En prenant le carré des deux termes de cette équation, on obtient une équation du second degré dont la seule solution ayant un sens physique est : ηy = (ξy.cos θ) + SQRT(1 – ξx ² – ξy ² ).sin θ En utilisant à nouveau la définition du sous-ensemble Γ, on peut finaliser l’expression du symplectomorphisme χθ comme suit : [Math.21] ^ = = ^. (cos (sin ^^ _^ On note que les variables η _x et η _y dans ces dernières expressions sont exprimées en fonction de et de ξy quelques lignes plus haut. On s’intéresse à présent aux conséquences de la transformation T de l’espace tridimensionnel de propagation E sur la représentation W _ψ ^Φ , W _ψ de l’hologramme numérique (également dénommée distribution). Les représentations Wψ ^Φ , Wψ de l’hologramme dans le domaine spatio- fréquentiel D peuvent être vues comme le symbole P d’un opérateur pseudo- différentiel Op(P). En poursuivant le parallèle avec les théories de la mécanique quantique, on pourra se référence à l’article "How Wigner Functions Transform Under Symplectic Maps", de Alex J. Dragt et Salman Habib, in Proceedings of the Advanced Beam Dynamics Workshop on Quantum Aspects of Beam Physics (Monterey, CA) 1998. Dans le premier exemple donné ci-dessus (transformée de Fourier fenêtrée), l’opérateur différentiel de symbole W _ψ ^Φ est Op(Wψ ^Φ ) = | Φ >< ψ | où on utilise la notation usuelle en mécanique quantique introduite dans l’article "A new notation for quantum mechanics", de P. A. M. Dirac, in Mathematical Proceedings of the Cambridge Philosophical Society, juillet 1939, p.416-418. Dans le second exemple (distribution de Wigner-Ville), l’opérateur différentiel de symbole W _ψ est : Op(Wψ) = | ψ >< ψ |. Concernant le troisième exemple (trames de Gabor), on remarque que la décomposition utilise des fonctions du même type que celle utilisées dans le premier exemple (transformées de Fourier fenêtrées), au sein d’une somme sur les éléments du maillage, et les développements ci-dessous s’appliquent donc du fait qu’ils s’appliquent aux différents termes de la somme. En continuant le parallèle avec les théories de l’analyse semi-classique (on pourra par exemple se référer à l’article "Uniform semiclassical estimates for the propagation of quantum observables" de A. Bouzouina et D. Robert in Duke Mathematical Journal, 111(2) 223-252, 1 ^er février 2002), les effets de la transformation T sur l’hologramme peuvent alors être rendus en appliquant l’opérateur de propagation Op(χ) (associé au symplectomorphisme χ) à l’opérateur pseudo- différentiel Op(P) de sorte que le symbole P’ de l’opérateur différentiel propagé Op(χ)Op(P)[Op(χ)] ^-1 forme une représentation de l’hologramme résultant de la transformation T. D’après le théorème d’Egorov (voir par exemple l’article de A. Bouzouina et D. Robert précité), le symbole P’ de l’opérateur différentiel propagé Op(χ)Op(P)[Op(χ)] ^-1 peut s’écrire : [Math.22] où χ est le symplectomorphisme associé à la transformation T comme expliqué plus haut, λ est la longueur d’onde de la lumière concernée comme déjà mentionné, et où les fonctions fn sont définies comme une somme d’intégrales d’opérateurs différentiels bi-différentiels ; par exemple, dans la quantification de Weyl utilisée, la fonction f ₁ est nulle et les fonctions f _n peuvent être définies comme suit pour n > 1 : [Math.23] où H _t et χ _t sont respectivement le hamiltonien et le flot hamiltonien à l’instant t (voir ci-dessous), i est ici comme précédemment le nombre imaginaire tel que i²=-1, et où le crochet de Poisson étendu est défini comme suit : [Math.24] où g et h sont des fonctions scalaires de (q, p) sur le domaine spatio-fréquentiel D. Ainsi, lorsque la fonction génératrice φ est non-quadratique, c’est-à-dire que le symplectomorphisme χ est non-linéaire comme expliqué dans l’article "How Wigner Functions Transform Under Symplectic Maps" précité, la représentation P’ de l’hologramme déformé par la transformation T ne se limite pas à une déformation du domaine spatio-fréquentiel D au moyen du symplectomorphisme χ (terme P(χ) dans la formule ci-dessus), mais comprend des termes correctifs additionnels qui dépendent chacun de la longueur d’onde λ, précisément qui incluent chacun un facteur égal à une puissance de la longueur d’onde λ. En pratique, comme indiqué plus bas, on peut utiliser un seul terme correctif additionnel qui inclut un facteur égal à la longueur d’onde λ ou, ici, égal au carré λ² de la longueur d’onde λ. Le hamiltonien Ht et le flot hamiltonien χt sont liés par les équations de Hamilton : [Math.25] En pratique, on peut décomposer le flux hamiltonien χt en deux flux hamiltoniens respectivement liés à une translation et à une rotation, et donc considérer un hamiltonien H ^trans lié à la translation et un hamiltonien H ^rot lié à la rotation. Le hamiltonien lié à la translation H ^trans correspond simplement à la propagation libre dans un milieu homogène d’indice de réfraction n ₀ (on peut prendre en pratique n ₀ = 1 pour une propagation dans l’air) : H ^trans = SQRT(n ₀ ²-|p|²). Concernant le hamiltonien lié à la rotation H ^rot , on peut obtenir l’expression de ses dérivées (qui seules sont nécessaires pour évaluer les fonctions f _n définies ci- dessus) en reportant l’expression littérale du symplectomorphisme χθ dans les formules de Hamilton rappelées ci-dessus : [Math.26] [Math.27] [Math.28] [Math.29] avec R=SQRT(1 – ξ _x ² – ξ _y ²). Les dérivées d’ordre supérieur peuvent être obtenues sur cette base par des méthodes numériques (telle que celle de différences finies). La figure 2 représente un exemple de dispositif d’holographie numérique 10 utilisant les enseignements de l’invention. Ce dispositif d’holographie numérique 10 comprend un module de traitement 2, un module de mémorisation 4, un module d’acquisition 6 et un module de restitution holographique 8. Le dispositif d’holographie numérique 10 est par exemple réalisé sous la forme d’un casque de visualisation ou de lunettes de réalité augmentée. Le module de traitement 2 est par exemple un processeur, ici un microprocesseur. Le module d’acquisition 6 comprend au moins un capteur de mouvement (par exemple de type accéléromètre ou gyroscope) qui transmet des signaux de mesure (sous forme analogique ou numérique) au module de traitement 2, sur la base desquels le module de traitement 2 peut élaborer des données représentatives de la position et/ou de l’orientation du dispositif d’holographie numérique 10. Le module de mémorisation 4 comprend par exemple une mémoire et/ou un disque dur. Le module de mémorisation 4 mémorise notamment des instructions de programme d’ordinateur conçues de sorte que le module de traitement 2 mette en œuvre certaines au moins des étapes du procédé décrit ci-dessous en référence à la figure 3 lorsque ces instructions sont exécutées par le module de traitement 2. Comme expliqué ci-dessous, le module de mémorisation 4 mémorise en outre une représentation W d’un hologramme numérique dans le domaine spatio-fréquentiel D déjà mentionné. En pratique, le module de mémorisation mémorise des valeurs W(L) respectivement associées à différents éléments L du domaine spatio-fréquentiel D (c’est-à-dire en pratique à différents éléments de ℝ ⁴ ). Le module de restitution holographique 8 comprend ici un modulateur de lumière et une source de lumière de longueur d’onde λ, et permet la restitution de l’hologramme numérique représenté par la représentation W à un utilisateur. On pourra par exemple se référer à la demande de brevet publiée en tant que WO2019/001968 à ce sujet. La figure 3 représente un exemple de procédé mis en œuvre par le module de traitement 2 du dispositif d’holographie numérique 10. Ce procédé débute par une étape E2 de réception de signaux de mesure en provenance du module d’acquisition 6. Le module de traitement 2 détermine alors des données représentatives de la position et/ou de l’orientation du dispositif d’holographie numérique 10 sur la base des signaux de mesure reçus (étape E4). Afin de présenter l’hologramme numérique à l’utilisateur en tenant compte du mouvement du dispositif d’holographie numérique 10 (et de simuler ainsi le mouvement réciproque de l’hologramme numérique par rapport à l’utilisateur), le module de traitement 2 applique à l’hologramme numérique (représenté par la représentation W comme mentionné ci-dessus) une transformation T définie dans l’espace de propagation tridimensionnel E en fonction des données représentatives précitées. Pour ce faire, le module de traitement 2 détermine à l’étape E6 le symplectomorphisme χ associé à cette transformation T grâce à la méthode proposée ci-dessus. Comme déjà indiqué, dans un plan Pz parallèle au plan de référence Pref, la fonction χ (symplectomorphisme) est telle que que le sous-ensemble des éléments de l’espace produit D x D de forme (e, χ(e)) est le sous-ensemble défini par les éléments de la forme (x, y, ηx, ηy) pour (x, y, ηx, ηy) variant dans ℝ ⁴ , où φ est une fonction génératrice non quadratique à quatre variables définie comme produit scalaire : - d’un premier vecteur, dont les deux premières coordonnées correspondent respectivement aux deux premières variables des quatre variables et dont la troisième coordonnée correspond à la coordonnée dudit plan parallèle P _z selon la direction z perpendiculaire au plan de référence Pref, et - d’un vecteur résultant de l’application de la transformation à un second vecteur, dont les deux premières coordonnées correspondent respectivement aux deux dernières variables des quatre variables et, ici, dont la norme est prédéfinie. Le module de traitement peut ainsi déterminer à l’étape E8 la représentation W’ de l’hologramme numérique transformé (ou hologramme numérique modifié par application de la transformation T) sur la base de la représentation W mémorisée dans le module de mémorisation 4 et de la fonction χ (symplectomorphisme) déterminé à l’étape E6. En variante, la représentation W (c’est-à-dire en pratique des données représentatives des valeurs W(L) précitées) pourrait être reçue par le dispositif d’holographie numérique 10 au moyen d’un module de communication (non représenté). Conformément à ce qui a été exposé précédemment, la représentation W’ de l’hologramme numérique transformé peut être déterminé comme suit : [Math.30] pour les différents L du domaine spatio-fréquentiel D utilisées dans la représentation concernée, où N est une valeur prédéterminée pour un dispositif d’holographie numérique donné (qui correspond au degré de précision souhaitée), λ est comme déjà indiqué la longueur d’onde associée à (et donc utilisée pour restituer) l’hologramme numérique, o est comme déjà indiqué l’opérateur de composition et où les fonctions fn sont définies plus haut. En pratique, on peut utiliser seulement le premier terme (dit "terme correctif") non-nul de la somme dans la formule ci-dessus ; le module de traitement détermine dans ce cas la représentation W’ de l’hologramme numérique modifié comme suit (le facteur f ₁ étant nul avec la quantification de Weyl utilisée ici) : [Math.31] avec [Math.32] 1 ^ ^ _^ = − ^{ ^ _^ , ^^^ _^ ^} ^^ ^^ 4 ^{^} _{^ ^^^} ^ où {.} ₃ est le crochet de Poisson défini plus haut. Le module de traitement 2 est ainsi conçu pour appliquer la fonction χ lors de l’application de la transformation T à l’hologramme numérique. La fonction χ est notamment appliquée lors de la composition de cette fonction avec la représentation W (terme W o χ). Par ailleurs, comme indiqué ci-dessus, on prévoit ici que la distribution W’ représentant l’hologramme numérique transformé par la transformation T soit déterminée par le module de traitement 2 en utilisant au moins un terme correctif qui dépend de la longueur d’onde λ. Ce terme correctif est non-nul du fait que la fonction χ est non-linéaire (la fonction génératrice φ étant non-quadratique comme expliqué plus haut). Le module de traitement 2 peut alors commander la restitution de l’hologramme numérique transformé (représenté par la représentation ou distribution W’) au moyen du module de restitution holographique 8. Pour ce faire, le module de traitement 2 reconstruit par exemple à l’étape E10 l’hologramme numérique transformé à partir de la représentation W’ obtenue à l’étape E8. Cette étape de reconstruction peut comprendre en pratique une sommation de contributions lumineuses respectivement associées aux valeurs W’(L) définissant la représentation W’ pour différents éléments L du domaine spatio-fréquentiel D. Le module de traitement 2 envoie alors à l’étape E12 des signaux représentatifs de l’hologramme numérique reconstruit à destination du module de restitution holographique 8 (en pratique au modulateur de lumière équipant ce module de restitution holographique 8) afin de permettre la restitution de l’hologramme numérique transformé à l’utilisateur, ici en illuminant, avec la source de lumière de longueur d’onde λ, le modulateur de lumière configuré sur la base des signaux représentatifs de l’hologramme numérique reconstruit. En variante, la restitution de l’étape E12 pourrait être réalisée sur un système stéréoscopique. L’étape E12 comprend alors une sous-étape de conversion stéréoscopique pour convertir l’hologramme reconstruit en signaux stéréoscopiques et une sous-étape d’affichage sur un écran du système stéréoscopique à partir des signaux stéréoscopiques. La figure 4 représente un procédé mis en œuvre dans un autre contexte d’utilisation de l’invention. Il s’agit d’un procédé de décodage d’un hologramme numérique Hn+1 à partir de données reçues dans un flux de données. Ce procédé est mis en œuvre dans un dispositif de décodage. Un tel dispositif de décodage comprend par exemple un processeur et une mémoire mémorisant des instructions de programme d’ordinateur conçues pour mettre en œuvre certaines étapes au moins du procédé de la figure 4 lorsque ces instructions sont exécutées par le processeur. Ce dispositif de décodage peut être un ordinateur (par exemple un ordinateur personnel) ou un dispositif d’holographie numérique tel que décrit ci-dessus en référence à la figure 2. Le procédé de la figure 4 débute par une étape E20 de réception d’un flux de données. Le dispositif de décodage reçoit par exemple ce flux de données en provenance d’un serveur distant, typiquement au moyen d’un module de communication du dispositif de décodage. Le procédé se poursuit par une étape E22 d’obtention d’un hologramme numérique de référence H _n . Le dispositif de décodage peut obtenir l’hologramme numérique de référence par décodage de données reçues dans le flux de données précité. Ces données sont par exemple représentatives d’une représentation Wn de l’hologramme numérique Hn (sans référence à un autre hologramme numérique). En variante, l’hologramme numérique H _n peut être obtenu au moyen d’un procédé du même type que celui décrit ici pour obtenir l’hologramme numérique Hn+1. Le procédé comprend également une étape E24 de décodage d’un résidu en à partir d’autres données reçues dans le flux de données. Le procédé comprend également une étape E26 d’obtention d’un prédicteur PRED, par exemple par lecture dans une mémoire du dispositif de décodage. En variante, le prédicteur PRED peut être obtenue par décodage de données contenues dans le flux de données reçu. Le procédé comprend alors une étape E28 de prédiction d’un hologramme numérique prédit H’ _n par application à l’hologramme numérique de référence H _n d’une transformation T définie par le prédicteur PRED, en utilisant les techniques décrites plus haut pour transformer la représentation P en représentation P’, qui comprennent donc l’application d’une fonction χ déterminée comme décrit ci-dessus. Comme dans le cadre du procédé de la figure 3, l’hologramme numérique de référence Hn est par exemple représenté par une représentation W dans le domaine spatio-fréquentiel D (c’est-à-dire par des valeurs W(L) respectivement associées à différents éléments L du domaine spatio-fréquentiel D). Dans ce cas, l’étape E28 de prédiction peut inclure une détermination d’une distribution W’ représentant l’hologramme numérique prédit H’n sur la base de la représentation W de l’hologramme numérique de référence H _n conformément à ce qui est exposé ci-dessus pour l’étape E8 de la figure 3, c’est-à-dire comme suit : [Math.33] pour les différents L du domaine spatio-fréquentiel D utilisés dans la représentation concernée, où N est une valeur prédéterminée pour un dispositif de décodage donné (qui correspond au degré de précision souhaitée, avec par exemple N = 1, auquel cas la formule contient un seul terme correctif λ.f _n (L)), λ est comme déjà indiqué la longueur d’onde associée aux hologrammes numériques, o est comme déjà indiqué l’opérateur de composition et où les fonctions fn sont définies plus haut. La fonction χ est ainsi notamment appliquée lors de la composition de cette fonction avec la représentation W (terme W o χ). Par ailleurs, comme indiqué ci-dessus, on prévoit ici que la distribution (ou représentation) W’ représentant l’hologramme numérique prédit H’n soit déterminée par le dispositif de décodage en utilisant au moins un terme correctif qui dépend de la longueur d’onde λ. Ce terme correctif est non-nul du fait que la fonction χ est non- linéaire comme expliqué plus haut. Le procédé peut alors comprendre une étape E30 d’obtention de l’hologramme décodé H _n+1 par combinaison du résidu e _n et de l’hologramme numérique prédit H’ _n . Par exemple, pour les différents L du domaine spatio-fréquentiel D utilisés dans la représentation concernée, la représentation (ou distribution) W’’ représentant l’hologramme numérique décodé Hn+1 est obtenue comme suit : W’’(L) = W’(L) + E(L) où E est la représentation du résidu en dans le domaine spatio-fréquentiel D. Lorsque le dispositif de décodage comprend un module de restitution holographique (par exemple lorsque le dispositif de décodage est un dispositif d’holographie numérique comme celui de la figure 2), le procédé peut en outre comprendre une étape de reconstruction de l’hologramme à restituer sur la base de la représentation W’’ de l’hologramme décodé Hn+1 et/ou une étape de E32 de restitution de l’hologramme numérique décodé (et reconstruit) au moyen du module de restitution holographique. La figure 5 représente un exemple de dispositif de codage 20 conforme à l’invention. Ce dispositif de codage est par exemple situé au sein d’un serveur de manière à délivrer des données représentatives d’une séquence d’hologrammes numériques à un dispositif électronique client, comme décrit dans la suite. Le dispositif de codage 20 comprend ici une première unité de codage 22, une unité de prédiction 24, une unité de détermination de résidu 26, une seconde unité de codage 28 et une unité de communication 30. Le dispositif de codage 20 est par exemple réalisé sous forme d’un dispositif électronique comprenant un processeur (tel qu’un microprocesseur) ; dans ce cas, chacune des unités précitées peut être mise en œuvre du fait de l’exécution, par ce processeur, d’instructions de programme d’ordinateurs mémorisées au sein du dispositif électronique formant le dispositif de codage 20 (par exemple dans une mémoire de ce dispositif électronique). Toutefois, en variante, une au moins des unités précitées peut être réalisée au moyen d’un circuit électronique dédié, par exemple un circuit intégré à application spécifique. La première unité de codage 22 est conçue pour produire des données D1 qui sont représentatives d’un hologramme numérique Hk, ici d’un hologramme numérique de référence, sans référence à des données relatives à un autre hologramme numérique. On décrit dans la suite en référence à la figure 10 un exemple de réalisation envisageable pour le codage effectué par la première unité de codage 22. L’unité de prédiction 24 est conçue pour prédire un hologramme numérique H’k, dit "hologramme numérique prédit", par transformation de l’hologramme numérique de référence Hk, en utilisant les techniques qui ont été présentées ci-dessus en référence aux figures 1 à 4, notamment en appliquant la fonction χ définie ci-dessus pour la détermination de l’hologramme numérique prédit. L’unité de détermination de résidu 26 est conçue pour déterminer un résidu ek+1 entre l’hologramme numérique prédit H’k produit par l’unité de prédiction 24 et l’hologramme numérique courant H _k+1 . La seconde unité de codage 28 est conçue pour produire des données D3, D4 représentatives du résidu ek+1 produit par l’unité de détermination de résidu 26. Des éléments complémentaires relatifs au fonctionnement de ces différentes unités sont données plus loin en référence notamment à la figure 7. Enfin, l’unité de communication 30 est conçue pour transmettre les données représentatives produites par la première unité de codage 22 et/ou par la seconde unité de codage 28 à destination du dispositif électronique client. La figure 6 représente un exemple de dispositif de décodage 40 conforme à l’invention. Ce dispositif de décodage 40 est par exemple un dispositif électronique (tel que le dispositif électronique client susmentionné) conçu pour recevoir des données représentatives d’une séquence d’hologramme numériques en provenance d’un serveur (tel que le serveur susmentionné). Le dispositif de décodage 40 comprend une première unité de décodage 42, une seconde unité de décodage 44, une unité de prédiction 46, une unité de combinaison 48, une unité de reconstruction 49 et une unité de communication 50. Le dispositif de décodage 40 est par exemple réalisé sous forme d’un dispositif électronique comprenant un processeur (tel qu’un microprocesseur) ; dans ce cas, chacune des unités précitées peut être mise en œuvre du fait de l’exécution, par ce processeur, d’instructions de programme d’ordinateurs mémorisées au sein du dispositif électronique formant le dispositif de décodage 40 (par exemple dans une mémoire de ce dispositif électronique). Toutefois, en variante, une au moins des unités précitées peut être réalisée au moyen d’un circuit électronique dédié, par exemple un circuit intégré à application spécifique. La première unité de décodage 42 est conçue pour obtenir un hologramme numérique Hk, ici un hologramme numérique de référence, à partir de données D1 représentatives de cet hologramme numérique. Ces données représentatives D1 sont par exemple reçues par l’unité de communication 50 (en provenance du dispositif de codage 20 dans les exemples décrits dans la suite). Un exemple de procédé de décodage qui peut être mis en œuvre par la première unité de décodage est décrit plus loin en référence à la figure 11. La seconde unité de décodage 44 est conçue pour décoder un résidu ek+1 à partir de données D3, D4 représentatives de ce résidu e _k+1 . Ces données D3, D4 représentatives du résidu sont ici reçues par l’unité de communication 50 (par exemple en provenance du dispositif de codage 20). L’unité de prédiction 46 est conçue pour prédire un hologramme numérique H’ _k , dit "hologramme numérique prédit", par transformation de l’hologramme numérique de référence Hk, en utilisant les techniques qui ont été présentées ci-dessus en référence aux figures 1 à 4, notamment en appliquant la fonction χ définie ci-dessus pour la détermination de l’hologramme numérique prédit. Son fonctionnement est semblable à celui de l’unité de prédiction 24 du dispositif de codage 20. L’unité de combinaison 48 est conçue pour obtenir un hologramme numérique H _k+1 (dit "hologramme numérique décodé") par combinaison du résidu e _k+1 décodé par la seconde unité de décodage 44, et de l’hologramme prédit H’k produit en sortie de l’unité de prédiction 46. L’unité de reconstruction 49 est conçue pour reconstruire l’hologramme numérique décodé H _k+1 sur la base d’une représentation de cet hologramme numérique, cette représentation utilisant ici une décomposition en trames de Gabor. Des éléments complémentaires relatifs au fonctionnement de ces différentes unités sont données plus loin en référence notamment à la figure 7. La figure 7 représente un procédé d’échange de données entre le dispositif de codage 20 et le dispositif de décodage 40 en vue par exemple de la restitution, sur un dispositif de restitution holographique, d’une séquence d’hologrammes numériques codés par le dispositif de codage 20, transmis du dispositif de codage 20 au dispositif de décodage 40 (par exemple via un réseau de communications), et décodés par le dispositif de décodage 40. Ce procédé débute par une étape E2 d’initialisation de l’indice k qui désigne l’hologramme numérique courant Hk dans la séquence d’hologrammes numériques traitée. L’indice k est ici initialisé à 0. Le procédé se poursuit par une étape E4 de décomposition de l’hologramme numérique courant Hk en trames de Gabor Wk. Comme précédemment indiqué, l’hologramme numérique courant H _k est alors représenté par un ensemble de premiers coefficients (ici complexes) W _k (ζ) respectivement associés à des éléments ζ du domaine spatio-fréquentiel D. On rappelle que chacun de ces éléments du domaine spatio-fréquentiel D est formé par un premier paramètre spatial x, un second paramètre spatial y, un premier paramètre fréquentiel η et un second paramètre fréquentiel ξ. Les coefficients Wk(ζ) peuvent par exemple être calculés par la formule suivante : Wk(ζ) = < ^^(ζ)Φ ^D , Hk > en reprenant les notations de la formule [Math.6] ci-dessus et avec < f , g > le produit scalaire entre les fonctions f et g. En pratique, la représentation de l’hologramme numérique courant Hk utilisée est par exemple une représentation parcimonieuse, c’est-à-dire que les coefficients Wk(ζ) sont choisis pour représenter au mieux l’hologramme numérique courant H _k tout en étant en nombre limité. Pour ce faire, on utilise par exemple un algorithme de poursuite de base (ou "basis pursuit" selon la terminologie anglo-saxonne couramment utilisée) pour sélectionner les coefficients W _k (ζ) issus de la décomposition de l’hologramme numérique courant H _k en trames de Gabor à utiliser dans la représentation parcimonieuse. On pourra se référer par exemple à l’article "Basis pursuit", de S. Chen et D. Donoho in Proceedings of 199428th Asilomar Conference on Signals, Systems and Computers (vol.1, pp.41-44), octobre 1994 en ce qui concerne l’algorithme de poursuite de base. Le procédé de la figure 7 se poursuit alors par une étape E6 de codage de l’hologramme numérique courant H _k , c’est-à-dire ici de codage des coefficients W _k (ζ) représentant l’hologramme numérique courant Hk. Cette étape de codage est ici effectuée par la première unité de codage 22 mentionnée ci-dessus. Un exemple de mise en œuvre envisageable pour l’étape E6 est décrit plus loin en référence à la figure 10. L’étape E6 permet de générer des données codées D1 représentatives de l’hologramme numérique courant H _k . Ces données codées D1 sont émises, par l’unité de communication 30 du dispositif de codage 20 et via un réseau de communication, à destination du dispositif de décodage 40 (étape E8). En variante, les données codées D1 pourraient être mémorisées au sein du dispositif de codage 20 (par exemple pour utilisation ultérieur dans le cas où le dispositif de codage 20 et le dispositif de décodage 40 sont un même et unique dispositif électronique). Dans les étapes E10 à E14 décrites à présent, l’hologramme numérique H _k d’indice k est utilisé comme hologramme numérique de référence pour le codage différentiel (c’est-à-dire sous forme d’un résidu) de l’hologramme numérique H _k+1 qui suit l’hologramme numérique Hk dans la séquence d’hologrammes numériques. Lors de l’étape E10, l’unité de prédiction 24 détermine une transformation T _k+1 (ici une transformation rigide comme déjà indiqué) sur la base de l’hologramme numérique de référence Hk et de l’hologramme numérique suivant Hk+1, puis détermine un hologramme numérique prédit H’k par transformation de l’hologramme numérique de référence H _k en utilisant la fonction (ou symplectomorphisme) χ associée comme expliqué plus haut à la transformation Tk+1 déterminée. La transformation T _k+1 est par exemple déterminée par analyse de mouvement entre la scène représentée par l’hologramme numérique de référence H _k et l’hologramme numérique suivant Hk+1, comme proposé dans la demande de brevet WO 2021/004797 La représentation W’ _k de l’hologramme numérique prédit H’ _k (au moyen de seconds coefficients W’ _k (ε)) est alors déterminée à partir de la représentation W _k de l’hologramme numérique de référence Hk (au moyen des premiers coefficients Wk(ε)) en utilisant la fonction χ associée à la transformation T _k+1 conformément à ce qui a été précédemment décrit, c’est-à-dire qu’on a selon une notation simplifiée : où o est l’opérateur de composition de fonctions, λ est la longueur d’onde d’enregistrement des hologrammes numériques et Pn des fonctions complexes définies sur le domaine spatio-fréquentiel D comme expliqué précédemment. Un coefficient associé à un élément (x ₀ , y ₀ , η ₀ , ξ ₀ ) et participant à la représentation de l’hologramme numérique prédit H’ _k est donc déterminé en utilisant la formule : Wk ( χ(x0, y0, η0, ξ0) ) + ^∑ ^ ^ _^^ ^ ^{^} ^ _^ (x0, y0, η0, ξ0). Ainsi, pour la détermination d’un second coefficient W’k(ε0) associé à un élément ε ₀ formé par une première valeur x ₀ relative au premier paramètre spatial, une seconde valeur y ₀ relative au second paramètre spatial, une troisième valeur η ₀ relative au premier paramètre fréquentiel et une quatrième valeur ξ0 relative au second paramètre fréquentiel, on utilise notamment le premier coefficient Wk associé à l’image, par la fonction χ, de l’élément ε ₀ formé par la première valeur x ₀ , la seconde valeur y ₀ , la troisième valeur η0 et la quatrième valeur ξ0. Selon la formule proposée ci-dessus pour calculer les coefficients W’k(ε) représentant l’hologramme numérique prédit H’ _k , on utilise un terme correctif qui dépend de la longueur d’onde λ d’acquisition de l’hologramme numérique. Ce terme correctif étant une somme de fonctions, on peut selon une variante envisageable utiliser un seul terme de cette somme en tant que terme correctif lors de la détermination des coefficients W’k(ε). Selon une autre variante, pour alléger les traitements, on peut ne pas tenir compte du terme correctif ; dans ce cas, tout second coefficient W’ _k (ε ₀ ) déterminé à l’étape E10 et associé à un élément ε0 formé par une première valeur x0 relative au premier paramètre spatial, une seconde valeur y ₀ relative au second paramètre spatial, une troisième valeur η ₀ relative au premier paramètre fréquentiel et une quatrième valeur ξ0 relative au second paramètre fréquentiel, est égal au premier coefficient Wk associé à l’image, par la fonction χ, de l’élément ε0 formé par la première valeur x0, la seconde valeur y ₀ , la troisième valeur η ₀ et la quatrième valeur ξ ₀ . Par ailleurs, il est également possible d’obtenir les coefficients W’ _k (ε) en multipliant l’expression donnée ci-dessus (avec ou sans utilisation du terme correctif), qui forme alors une valeur intermédiaire, par une fonction de distribution spatio- fréquentielle Δ à valeurs réelles sur le domaine spatio-fréquentiel D et variable en fonction seulement selon les deux dimensions fréquentielles (c’est-à-dire qu’il existe une fonction f à deux variables telle que Δ(x,y,η,ξ)=f(η,ξ) pour tout multiplet (x,y,η,ξ)), par exemple pour obtenir une atténuation de l’hologramme numérique dans une certaine direction. Cette fonction de distribution spatio-fréquentielle est par exemple de la forme : [Math.34] avec θ une valeur prédéterminée liée à la direction d’atténuation, et la formule générale qui donne l’expression des coefficients W’k représentant l’hologramme prédit H’ _k est alors : W’ _k = [W _k o χ + ∑ ^ ^ _^^ ^ ^{^} ^ _^ ].Δ. Le procédé comprend alors une étape E12 de décomposition de l’hologramme numérique H _k+1 en trames de Gabor W _k+1 . Cette étape est similaire à l’étape E4 réalisée pour l’hologramme numérique Hk et permet d’obtenir, pour représenter l’hologramme numérique H _k+1 , un ensemble de coefficients (ici complexes) Wk+1(ζ) respectivement associés à des éléments ζ du domaine spatio-fréquentiel D. Le procédé de la figure 7 se poursuit par une étape E14 de détermination d’un résidu ek+1 entre l’hologramme numérique Hk+1 et l’hologramme numérique prédit H’k, ici par différence entre l’hologramme numérique Hk+1 et l’hologramme numérique prédit H’k. Cette étape est ici réalisée par l’unité de détermination de résidu 26. En pratique, cette différence est faite entre les coefficients W _k+1 (ζ) représentant l’hologramme numérique Hk+1 et les coefficients W’k(ζ) représentant l’hologramme numérique prédit H’k, et on note ainsi de manière simplifiée : e _k+1 = W _k+1 – W’ _k Plus précisément, pour tous les éléments ε du domaine spatio-fréquentiel D pour lesquels l’un au moins des coefficients Wk+1(ε) et W’k(ε) est non-nul, l’unité de détermination de résidu 26 détermine un coefficient (complexe) e _k+1 (ε) du résidu e _k+1 par l’opération : ek+1(ε) = Wk+1(ε) – W’k(ε) L’unité de communication 30 peut alors transmettre à l’étape E16 des données représentatives D2, D3, D4 de la transformation Tk+1 et du résidu ek+1 à destination du dispositif de décodage 40 (via le réseau de communication précité). En variante, ces données représentatives pourraient être mémorisées dans une mémoire du dispositif de codage 20 pour utilisation ultérieure. Les données D2 représentatives de la transformation Tk+1 comprennent des valeurs qui caractérisent la transformation Tk+1, par exemple les valeurs des composantes du vecteur définissant la translation impliquée et les valeurs d’angle définissant les différentes rotations impliquées. On propose de coder comme suit le résidu ek+1 pour sa transmission au dispositif de décodage 40. En variante, on pourrait utiliser la méthode de codage décrite plus bas en référence à la figure 10. On traite ici séparément la partie réelle et la partie imaginaire de chacun des coefficients e _k+1 (ε). On forme des images I _i,j qui regroupent chacune la partie réelle des coefficients ek+1(ε) associés à des éléments dont le premier paramètre fréquentiel et le second paramètre fréquentiel sont fixes. Autrement dit, un pixel en ligne p et colonne q d’une image I _i,j a pour valeur : Ii,j(p,q) = Re[ ek+1 ( (xp, yp, ηi, ξj) ) ] où Re[γ] désigne la partie réelle du nombre complexe γ et η _i et ξ _j sont respectivement le premier paramètre fréquentiel et le second paramètre fréquentiel associés à l’image Ii,j. La séquence d’images formée par les images Ii,j (lorsque i et j varient de manière à balayer tous les couples premier paramètre fréquentiel η _i - second paramètre fréquentiel utilisés dans la représentation) est alors codée au moyen d’un codeur vidéo (par exemple conformément à une norme de type AVC ou HEVC). La séquence d’images précitée peut être formée en prenant les images I _i,j dans l’ordre suivant : les différentes images I1,j avec j croissant, puis les différentes images I2,j avec j décroissant, puis les différentes images I3,j avec j croissant, etc. Les données codées D3 produites par le codage de la séquence d’images I _i,j sont envoyées au dispositif de décodage 40 par le circuit de communication 30. De la même manière, on forme des images Ji,j qui regroupent chacune la partie imaginaire des coefficients e _k+1 (ε) associés à des éléments dont le premier paramètre fréquentiel et le second paramètre fréquentiel sont fixes. Autrement dit, un pixel en ligne p et colonne q d’une image Ji,j a pour valeur : Ji,j(p,q) = Im[ ek+1 ( (xp, yp, ηi, ξj) ) ] où Im[γ] désigne la partie imaginaire du nombre complexe γ et η _i et sont respectivement le premier paramètre fréquentiel et le second paramètre fréquentiel associés à l’image Ji,j. La séquence d’images formée par les images J _i,j (lorsque i et j varient de manière à balayer tous les couples premier paramètre ηi - second paramètre ξj utilisés dans la représentation) est alors codée au moyen d’un codeur vidéo (par exemple conformément à une norme de type AVC ou HEVC). La séquence d’images précitée peut être formée en prenant les images J _i,j dans l’ordre suivant : les différentes images J1,j avec j croissant, puis les différentes images J _2,j avec j décroissant, puis les différentes images J _3,j avec j croissant, etc. Les données codées D4 produites par le codage de la séquence d’images J _i,j sont envoyées au dispositif de décodage 40 par le circuit de communication 30. Le procédé de la figure 7 se poursuit par une étape E18 d’incrémentation (d’une unité) de l’indice k pour traitement de l’hologramme numérique suivant. Le procédé comprend alors une étape E20 où le dispositif de codage 20 détermine si une réinitialisation du codage différentiel est nécessaire c’est-à-dire si le nouvel hologramme numérique H _k+1 à coder doit être codé entièrement (comme à l’étape E6) ou peut encore être codé par codage différentiel. Le critère pour déterminer si une réinitialisation est nécessaire est par exemple un critère de fréquence (de manière à obtenir un nombre prédéterminé d’hologramme numériques successifs codés par codage différentiel) ou un critère de limitation de l’erreur (en limitant à un seuil prédéterminé une mesure d’erreur entre l’hologramme numérique d’origine et l’hologramme numérique reconstruit après décodage comme expliqué dans la suite). Si le dispositif de codage 20 détermine à l’étape E20 qu’une réinitialisation est nécessaire, le procédé se poursuit (flèche P) à l’étape E22 pour une nouvelle incrémentation de l’indice k (d’une unité) puis boucle à l’étape E4 pour le traitement du nouvel hologramme numérique, à présent désigné H _k . Si le dispositif de codage 20 détermine à l’étape E20 qu’une réinitialisation n’est pas nécessaire, le procédé se poursuit (flèche N) à l’étape E10 pour codage différentiel du nouvel hologramme numérique Hk+1. On décrit à présent (toujours en référence à la figure 7) le décodage des données transmises par le dispositif de codage 20 comme expliqué ci-dessus en vue de la reconstruction de la séquence d’hologramme numériques. À réception des données codées D1 représentatives de l’hologramme numérique Hk, la première unité de décodage 42 décode les données codées D1 (étape E30) de manière à obtenir les coefficients W _k (ζ) formant la représentation de l’hologramme numérique Hk. Un exemple de procédé de décodage envisageable est décrit plus loin en référence à la figure 11. L’unité de reconstruction 49 peut alors reconstruire l’hologramme numérique H _k à l’étape E32 sur la base des coefficients Wk(ζ) décodés à l’étape E30. On utilise pour ce faire la transformation inverse de celle utilisée à l’étape E4 pour obtenir la représentation de l’hologramme numérique. Dans l’exemple décrit ici, on utilise donc la transformation inverse à celle utilisée pour obtenir la décomposition en trames de Gabor. L’hologramme numérique H _k est donc reconstruit ici selon la formule [Math.6] présenté plus haut (en utilisant les coefficients décodés Wk(ζ) en tant que valeurs de ^^ ^{^} ^ dans cette équation). L’hologramme numérique H _k peut alors être restitué au moyen d’un dispositif de restitution holographique (étape E34). Un tel dispositif de restitution holographique comprend par exemple un modulateur de lumière et une source de lumière de longueur d’onde λ. En appliquant respectivement les valeurs des pixels de l’hologramme numérique H _k aux pixels du modulateur de lumière et en illuminant le modulateur de lumière au moyen de la source de lumière, on restitue l’hologramme numérique H _k . Ultérieurement, à réception des données codées D2, D3, D4, la seconde unité de décodage 44 décode ces données à l’étape E40 de manière à obtenir la transformation Tk+1 et le résidu ek+1. La seconde unité de décodage 44 détermine la transformation T _k+1 à partir des données codées D2, ici à partir des valeurs des composantes du vecteur définissant la translation caractéristique de la transformation Tk+1 et des valeurs d’angle définissant les différentes rotations caractéristiques de la transformation T _k+1 . La seconde unité de décodage 44 reconstruit par ailleurs la séquence d’images Ii,j à partir des données codées D3 et peut ainsi déterminer les parties réelles respectives des coefficients ek+1(ε) représentant le résidu ek+1. De même, la seconde unité de décodage 44 reconstruit la séquence d’images Ji,j à partir des données codées D4 et peut ainsi déterminer les parties imaginaires respectives des coefficients ek+1(ε) représentant le résidu ek+1. La seconde unité de décodage 44 peut ainsi déterminer les coefficients (complexes) ek+1(ε), c’est-à-dire décoder le résidu ek+1, à partir des données D3, D4 reçues du dispositif de codage 20. L’unité de prédiction 46 peut par ailleurs prédire l’hologramme prédit H’k à l’étape E42 par transformation d’un hologramme numérique de référence, par exemple l’hologramme numérique Hk décodé à l’étape E30 ou un hologramme numérique obtenu lors d’un précédent passage à l’étape E44 décrite ci-dessous. Lors de cette étape E42, les coefficients W’ _k (ε) de l’hologramme numérique prédit H’k sont déterminés sur la base des coefficients Wk(ε) représentant l’hologramme numérique de référence Hk et de la transformation Tk+1 (décodée à l’étape E40 sur la base des données D2 reçues du dispositif de codage 20), selon un processus identique à celui employé à l’étape E10 décrite plus haut et qui ne sera donc pas décrit à nouveau. La détermination de l’hologramme numérique prédit H’k lors du décodage utilise donc ici aussi la fonction χ associée comme expliqué plus haut à la transformation Tk+1. A l’étape E44, l’unité de combinaison 48 peut ainsi obtenir l’hologramme numérique décodé Hk+1 par combinaison du résidu ek+1 (décodé à l’étape E40) et de l’hologramme numérique prédit H’ _k , ici par sommation (coefficient à coefficient) du résidu ek+1 et de l’hologramme numérique prédit H’k : H _k+1 = H’ _k + e _k+1 . Autrement dit, pour tout élément ε du domaine spatio-fréquentiel D correspondant à un coefficient dans la représentation du résidu ek+1 ou de l’hologramme numérique prédit H’k, l’unité de combinaison 48 détermine le coefficient Wk+1(ε) associé à l’élément ε dans la représentation de l’hologramme numérique décodé H _k+1 comme suit : Hk+1(ε) = H’k(ε) + ek+1(ε) . L’unité de reconstruction 49 peut alors reconstruire l’hologramme numérique Hk+1 à l’étape E46 sur la base des coefficients Wk+1(ε) obtenus à l’étape E44, de la même manière que cela a été fait pour l’hologramme numérique Hk à l’étape E32. L’hologramme numérique Hk+1 peut alors être restitué au moyen du dispositif de restitution holographique (étape E48), ici en appliquant respectivement les valeurs des pixels de l’hologramme numérique Hk+1 aux pixels du modulateur de lumière et en illuminant le modulateur de lumière au moyen de la source de lumière. La figure 8 représente une variante de réalisation du procédé de la figure 7 selon laquelle on considère plusieurs objets Om (ici M objets O1, O2, …, OM) dans la scène tridimensionnelle représentée par les hologrammes numériques. Le procédé de la figure 8 peut lui aussi être mis en œuvre par la coopération du dispositif de codage 20 et du dispositif de décodage 40. Le procédé de la figure 8 commence également par l’initialisation (par exemple à la valeur 0) d’un indice k désignant l’hologramme numérique en cours de traitement (étape E50). Le dispositif de codage 20 effectue alors une étape E54 de décomposition de l’hologramme numérique courant H _k en trames de Gabor W _k . Cette étape E54 est identique à l’étape E4 décrite plus haut en référence à la figure 7 et ne sera donc pas décrite en détail. L’étape E54 permet ainsi d’obtenir des coefficients (ici complexes) W _k (ζ) respectivement associés à des éléments ζ du domaine spatio-fréquentiel D et représentant l’hologramme numérique courant Hk. Le dispositif de codage 20 procède alors à une étape E56 de codage des coefficients W _k (ζ) représentant l’hologramme numérique H _k . Cette étape de codage est identique à l’étape E6 décrit plus haut et ici effectuée par la première unité de codage 22 mentionnée ci-dessus. Cette étape de codage utilise par exemple un codage tel qu’exposé plus loin en référence à la figure 10. L’étape de codage E56 permet de produire des données codées D1 représentatives des coefficients Wk(ζ). L’unité de communication 30 peut alors envoyer à l’étape E58 les données D1’ représentatives des coefficients W _k (ζ) à destination du dispositif de décodage 40. L’étape E58 est ici identique à l’étape E8 décrite plus haut. Le procédé comprend alors une étape E60 (identique à l’étape E12 du procédé de la figure 7) de décomposition de l’hologramme numérique H _k+1 en trames de Gabor Wk+1. Cette étape E60 permet d’obtenir, pour représenter l’hologramme numérique Hk+1, un ensemble de coefficients (ici complexes) Wk+1(ζ) respectivement associés à des éléments ζ du domaine spatio-fréquentiel D. Le procédé de la figure 8 se poursuit par une étape E61 de partition du domaine spatio-fréquentiel D en une pluralité de parties Pm disjointes respectivement associées aux objets Om. Un exemple de procédé permettant d’obtenir une telle partition est donné ci- dessous en référence à la figure 9. Le dispositif de codage 20 peut alors former à l’étape E62 des ensembles disjoints Em en regroupant dans un ensemble Em donné les coefficients Wk+1(ζ) associés à des éléments ζ appartenant à une partie P _m du domaine spatio-fréquentiel D associée à cet ensemble Em donné. Autrement dit, un coefficient W _k+1 (ζ) appartient à un ensemble E _m si et seulement si l’élément ζ appartient à la partie P _m . Les coefficients Wk+1(ζ) d’un ensemble Em représentent ainsi un hologramme numérique partiel Hk+1,m qui correspond à l’objet Om. Dans l’exemple décrit ici (et comme expliqué plus bas lors de la description de la figure 9), on propose de déterminer comme suit à quel ensemble Em appartient un coefficient Wk+1(ζ) associé à un élément ζ = (x,y,η,ξ) du domaine spatio-fréquentiel D. On convertit les coordonnées fréquentielles (η,ξ) en coordonnées spatiales (px,py) avec : px = 1/tan(arcsin(λ.η)) et py = 1/tan(arcsin(λ.ξ)), puis on détermine à quel sous-ensemble Gx,y,m appartient le couple (px, py), par exemple en calculant de quel côté de chaque droite (s _x,y,m ,i, s _x,y,m,i+1 ) se trouve le point (px, py). (On rappelle que les sous-ensembles Gx,y,m, délimités par les segments [sx,y,m,i, s _x,y,m,i+1 ] sont définis par la mise en œuvre du procédé décrit ci-dessous en référence à la figure 9 lors de l’étape E61.) Le coefficient Wk+1(ζ) associé à l’élément ζ = (x,y,η,ξ) sera alors attribué à l’ensemble Em associé au sous-ensemble Gx,y,m. Le procédé se poursuit par une étape E63 d’initialisation (ici à la valeur 1) d’un indice m désignant l’objet O _m (et donc également l’ensemble E _m ) en cours de traitement. L’unité de prédiction 24 prédit alors à l’étape E64 un hologramme numérique prédit partiel H’k,m par transformation d’un hologramme numérique de référence partiel Hk,m défini par les coefficients Wk(χm(ζ)) associés aux éléments χm(ζ) pour ζ variant dans la partie Pm associée à l’objet Om. Comme à l’étape E10 du mode de réalisation de la figure 7, cette étape E64 utilise une fonction χm associée à la transformation Tk+1,m associée à l’objet Om conformément à ce qui est décrit plus haut pour le symplectomorphisme χ. La transformation T _k+1,m est donc ici déterminée au moyen d’un algorithme d’estimation du mouvement de l’objet Om. Comme déjà indiqué, on pourra utiliser la solution proposée la demande de brevet WO 2021/004797. Ici l’unité de prédiction 24 détermine les coefficients W’k,m(ε) de l’hologramme prédit partiel H’ _k,m comme suit : pour les éléments ε appartenant la partie P _m du domaine spatio-fréquentiel D, c’est-à-dire de manière à prédire les coefficients Wk+1(ε) appartenant à l’ensemble Em (ces coefficients sont notés Wk+1,m(ε) dans la suite). L’unité de détermination de résidu 26 peut alors déterminer un résidu partiel ek+1,m entre l’hologramme numérique partiel Hk+1,m et l’hologramme numérique prédit partiel H’k,m, ici par différence entre l’hologramme numérique partiel Hk+1,m et l’hologramme numérique prédit partiel H’ _k,m ., soit pour tout élément ε appartenant la partie P _m du domaine spatio-fréquentiel D (ici pour tout coefficient W _k+1,m (ε) appartenant à l’ensemble Em) : ek+1,m(ε) = Wk+1,m(ε) - W’k,m(ε). L’unité de communication 30 transmet alors à l’étape E68 des données D2’ représentatives de la transformation T _k+1,m et du résidu partiel e _k+1,m . Pour obtenir ces données D2’, on utilise par exemple un codage identique à celui décrit plus haut pour le codage de la transformation T _k et du résidu e _k+1 . En particulier, selon une possibilité de réalisation, le résidu partiel e _k+1,m peut être codé au moyen de la méthode de codage décrite ci-dessous en référence à la figure 10. Le dispositif de codage 20 détermine alors à l’étape E70 si tous les objets Om ont été traités, c’est-à-dire ici si l’on a m=M. Dans la négative (flèche N), l’indice m est incrémenté de 1 (étape E72) et le procédé boucle à l’étape E64 décrite ci-dessus. Dans l’affirmative (flèche P), le procédé se poursuit à l’étape E74 décrite à présent. A l’étape E74, l’indice k est incrémenté de 1. Le procédé comprend alors une étape E76 (identique à l’étape E20 décrite plus haut) où le dispositif de codage 20 détermine si une réinitialisation du codage différentiel est nécessaire, c’est-à-dire si le nouvel hologramme numérique Hk+1 à coder doit être codé entièrement ou peut encore être codé par codage différentiel. On pourra se référer à la description de l’étape E20 en ce qui concerne le critère pour déterminer si une réinitialisation est nécessaire. Si le dispositif de codage 20 détermine à l’étape E76 qu’une réinitialisation est nécessaire, le procédé se poursuit (flèche P) à l’étape E78 pour une nouvelle incrémentation de l’indice k (d’une unité) puis boucle à l’étape E54 pour le traitement du nouvel hologramme numérique, à présent désigné Hk. Si le dispositif de codage 20 détermine à l’étape E76 qu’une réinitialisation n’est pas nécessaire, le procédé se poursuit (flèche N) à l’étape E60 pour codage différentiel du nouvel hologramme numérique Hk+1. On décrit à présent toujours en référence à la figure 8 le décodage des données D1’, D2’ transmises par le dispositif de codage 20 comme expliqué ci-dessus en vue de la reconstruction de la séquence d’hologramme numériques. Ce procédé de décodage comprend une étape E80 de décodage des données D1’ telles que reçues par l’unité de communication 50 de manière à obtenir les ensembles disjoints Em et les coefficients Wk(ζ) contenus dans chacun de ces ensembles disjoints Em. L’unité de reconstruction 49 peut alors reconstruire l’hologramme numérique H _k à l’étape E82 sur la base de tous les coefficients Wk(ζ) décodés à l’étape E80 (c’est- à-dire des coefficients contenus dans tous les ensembles E _m ). On remarque que les informations relatives aux ensembles disjoints E _m ne sont pas utilisées à cette étape, mais seront utilisées dans la suite. Cette étape E82 est similaire à l’étape E32 décrite dans le cadre de la description de la figure 8 et ne sera donc pas décrite en détail. L’hologramme numérique H _k peut alors être restitué à l’étape E84 au moyen d’un dispositif de restitution holographique du même type que celui utilisé à l’étape E34 décrite plus haut. On décrit à présent le décodage d’un hologramme numérique Hk+1 codé de manière différentielle par rapport à un hologramme numérique de référence Hk (tel que par exemple l’hologramme numérique décodé à l’étape E80), ici au moyen des données D2’ mentionnées plus haut. Le dispositif de codage 20 transmet par exemple des données de signalisation au dispositif de décodage 40 (au sein du flux de données émis par le module de communication 30 et reçu par le module de communication 50) pour indiquer ce codage différentiel du prochain hologramme numérique de la séquence. Ces données de signalisation peuvent par exemple indiquer le nombre M d’objets présents dans la scène tridimensionnelle. Le dispositif de décodage 40 initialise alors un indice m à 1 (étape E86). Le dispositif de décodage 40 décode alors à l’étape E88 les données D2’ reçues de manière à obtenir la transformation T _k+1,m et le résidu partiel e _k+1,m . L’unité de prédiction 46 peut alors prédire à l’étape E90 l’hologramme numérique prédit partiel H’k,m sur la base de la transformation Tk+1,m et de l’hologramme numérique de référence partiel Hk,m défini par les coefficients Wk(χm(ε)) tels que l’élément ε appartient à la partie P _m du domaine spatio-fréquentiel D associée à l’objet courant Om. Les traitements effectués lors de cette étape E90 sont identiques à ceux effectués lors de l’étape E64 au cours du codage et ne seront donc détaillés à nouveau ici. Ces traitements permettent d’obtenir les coefficients W’ _k,m (ε) représentant l’hologramme numérique prédit partiel H’k,m. L’unité de combinaison 48 obtient alors l’hologramme numérique partiel Hk+1,m par combinaison du résidu partiel ek+1,m (décodé à l’étape E88) et de l’hologramme numérique prédit partiel H’ _k,m , ici par sommation (coefficient à coefficient) du résidu partiel ek+1,m et de l’hologramme numérique prédit partiel H’k,m : H _k+1,m = H’ _k,m + e _k+1,m . Autrement dit, pour tout élément ε appartenant à la partie P _m du domaine spatio- fréquentiel D, l’unité de combinaison 48 détermine le coefficient Wk+1,m(ε) associé à l’élément ε dans la représentation de l’hologramme numérique partiel Hk+1,m comme suit : W _k+1,m (ε) = W’ _k,m (ε) + e _k+1,m (ε) . Le dispositif de décodage 40 détermine alors à l’étape E94 si tous les objets Om ont été traités, c’est-à-dire ici si l’on a m=M. Dans la négative (flèche N), l’indice m est incrémenté de 1 (étape E96) et le procédé boucle à l’étape E88 décrite ci-dessus. Dans l’affirmative (flèche P), le procédé se poursuit à l’étape E74 par une étape de fusion de tous les coefficients W _k+1,m (ε) obtenus au cours des différents passages à l’étape E92 et qui représentent (ensemble) l’hologramme numérique décodé Hk+1. L’unité de reconstruction 49 peut alors reconstruire l’hologramme numérique H _k+1 à l’étape E100 sur la base des coefficients fusionnés à l’étape E98, de la même manière que cela a été fait pour l’hologramme numérique Hk à l’étape E82. L’hologramme numérique Hk+1 peut alors être restitué au moyen du dispositif de restitution holographique (étape E102), ici en appliquant respectivement les valeurs des pixels de l’hologramme numérique H _k+1 aux pixels du modulateur de lumière et en illuminant le modulateur de lumière au moyen de la source de lumière. On décrit à présent en référence à la figure 9 un procédé de partitionnement du domaine spatio-fréquentiel D. Dans l’exemple décrit ici, on partitionne le domaine spatio-fréquentiel D en parties Pm en partitionnant, pour chaque point de l’hologramme numérique de coordonnées (x,y), un ensemble M _x,y contenant tous les couples (u,v) de variables fréquentielles utilisés dans la représentation de l’hologramme numérique concerné, c’est-à-dire toutes les directions associées à ces couples (u,v), ou, ce qui est équivalent, en partitionnant une partie d’un plan PPN parallèle au plan de l’hologramme numérique, comme décrit ci-dessous. Le procédé débute ainsi en initiant une boucle qui parcourt l’ensembles des points X de coordonnées (x,y) de l’hologramme numérique (étape E122). Pour chaque point X, le procédé initie à l’étape 124 une boucle qui parcourt l’ensemble des objets Om de l’objet le plus éloigné du plan de l’hologramme à l’objet le plus proche. Dans certains cas d’application, les objets, leurs positions et leurs formes sont connus. C’est le cas par exemple lorsque l’hologramme numérique concerné a été construit à partir d’une représentation de la scène tridimensionnelle comprenant ces objets. Dans d’autres cas d’application, les objets sont reconstruits à partir de l’hologramme numérique concerné, en utilisant par exemple les techniques décrites dans l’article "Comparative analysis of autofocus functions in digital in-line phase- shifting holography", E. S. R. Fonseca, P. T. Fiadeiro, M. Pereira et A. Pinheiro, , Appl. Opt., AO, vol.55, n° 27, p.7663-7674, septembre 2016. Pour chaque objet Om, le procédé initie à l’étape E126 une boucle qui parcourt tous les points P de l’objet Om. On suppose que l’objet est discrétisé un nombre fini de points et on note (P _x , P _y , P _z ) les coordonnées du point P courant. Pour chaque point P de l’objet Om, on détermine à l’étape E128 le point P’ intersection du segment [PX] et du plan PPN parallèle au plan (de référence) de l’hologramme numérique. Ici on prend pour plan PPN le plan d’équation z=1 (étant rappelé que le plan de référence est le plan d’équation z=0). Ici, les coordonnées du point P’ sont donc ((Px-x)/Pz,(Py-y)/Pz,1). On détermine ensuite à l’étape E130 un point P’’ de l’hologramme numérique de coordonnées (P’’ _x ,P’’ _y ) tel qu’aucun autre point de l’hologramme numérique ne soit plus proche du point P’ que le point P’’. On peut alors remplir à l’étape E132 un élément d’une matrice Gx,y de mêmes dimensions que l’hologramme numérique : on fixe à la valeur m (correspondant à l’objet O _m courant) l’élément de coordonnées (P’’ _x ,P’’ _y ) dans la matrice G _x,y . Si cet élément a déjà été rempli lors d’un précédent passage à l’étape E132, la valeur précédente est écrasée (c’est-à-dire remplacée par la valeur courante de la variable m). Une fois les boucles sur les objets Om et leurs points P mentionnées ci-dessus terminées, on procède à l’étape E134 à laquelle on définit des ensembles Gx,y,m : Gm = {(x’,y’) : Gx,y(x’,y’) = m}. Autrement dit, un ensemble Gx,y,m donné est l’ensemble des couples (x’,y’) pour lesquels l’élément de coordonnées (x’,y’) dans la matrice Gx,y vaut m. Les ensembles G _x,y,m sont dès lors des ensembles disjoints et correspondent chacun (grâce à la projection perspective de centre (x,y,0) effectuée à l’étape E128) à l’ensemble des directions où est situé l’objet O _m concerné (depuis le point de coordonnées (x,y,0)), c’est-à-dire à l’ensemble des couples de coordonnées fréquentielles (u,v) de l’ensemble Mx,y associés à ces directions. La définition des ensembles Gx,y,m revient ainsi à un partitionnement de l’ensemble M _x,y et ainsi, lorsque cela est réalisé pour tous les couples (x,y), à un partitionnement du domaine spatio-fréquentiel D en parties P _m , avec l’équivalence suivante : (x,y,u,v) appartient à Pm si et seulement si le couple (1/tan(arcsin(λ.u)) , 1/tan(arcsin(λ.v))) appartient à G _x,y,m . Le procédé se poursuit par une étape E136 de détermination, pour chaque ensemble Gx,y,m, d’un ensemble de segments [sx,y,m,i,sx,y,m,i+1] délimitant l’ensemble Gx,y,m. Par exemple, on parcourt l’ensemble des points bordant l’ensemble Gx,y,m (c’est-à-dire voisins avec un point qui n’appartient pas à l’ensemble G _x,y,m ) et on utilise chacun de ces points comme l’un des points sx,y,m,i définissant une borne d’un segment [sx,y,m,i,sx,y,m,i+1]. Après l’étape E136, le procédé boucle éventuellement à l’étape E124 pour traitement d’un nouveau point X de l’hologramme numérique, jusqu’à ce que tous les points de l’hologramme numérique aient été traités. On décrit à présent en référence à la figure 10 un procédé de codage de coefficients associés respectivement à des éléments du domaine spatio-fréquentiel D et utilisés pour représenter un hologramme numérique. Dans l’exemple décrit ici, il s’agit du codage des coefficients W _k (ζ) représentant un hologramme numérique H _k de manière à obtenir des données D1 représentatives de cet hologramme numérique H _k . Dans le contexte décrit plus haut, ce procédé est mis en œuvre par la première unité de codage 22 du dispositif de codage 20. En variante toutefois, ce procédé pourrait être utilisé pour coder d’autres types de coefficients, par exemple les coefficients formant le résidu ek+1 mentionné ci-dessus. Comme indiqué plus haut, l’hologramme numérique H _k est ici représenté par une pluralité de coefficients Wk(ζ) associés respectivement à des éléments ζ qui font partie d’un maillage Λ du domaine spatio-fréquentiel D (voir l’équation [Math 5] ci-dessus pour la définition du maillage Λ). On rappelle que les coefficients Wk(ζ) sont des nombres complexes. Dans la suite, on note ^ _^ (^ _^ entier variant de 1 à Nx) les indices respectivement associés aux valeurs ^ _^^ du premier paramètre spatial utilisées dans la représentation de l’hologramme numérique Hk, ^ _^ ( ^ _^ entier variant de 1 à Ny) les indices respectivement associés aux valeurs ^ _^^ du second paramètre spatial utilisées dans cette même représentation, ^ _^ ( ^ _^ entier variant de 1 à Mx) les indices respectivement associés aux valeurs du premier paramètre fréquentiel utilisées dans la représentation et ^ _^ (^ _^ entier variant de 1 à M _y ) les indices respectivement associés aux valeurs ^ _^^ du second paramètre fréquentiel utilisées dans la représentation. Pour faciliter l’exposé qui suit, on définit des blocs ^ _^^,^^ comme suit (pour tous les quadruplets (^ _^ , ^ _^ , ^ _^ , ^ _^ ) utilisés dans la représentation) : ) si la représentation de H _k comprend un coefficient associé ^ _^^,^^ (^ _^ , ^ _^ ) = 0 sinon. Autrement dit, chaque bloc ^ _^^,^^ comprend (tous) les coefficients associés aux éléments ayant ^ _^^ pour valeur du premier paramètre fréquentiel et ^ _^^ pour valeur du second paramètre fréquentiel, organisés sous forme de matrice en fonction des valeurs respectives du premier paramètre spatial et du second paramètre spatial. L’élément défini par la position (^ _^ , ^ _^ ) dans le bloc ^ _^^,^^ est nul lorsqu’aucun coefficient n’est associé à l’élément ^ _^^ , ^ _^^ , ^ _^^ ) dans la représentation de l’hologramme numérique H _k . Le procédé de la figure 10 débute par une étape E150 de codage des positions des coefficients non-nuls dans un premier bloc ^ _^^,^^ avec ici ^ _^ = 1 et ^ _^ = 1. La définition d’un bloc ^ _^^,^^ (telle que proposée ici) est une manière particulière d’associer de manière prédéfinie un ensemble des positions (dans le bloc) aux coefficients associés respectivement aux différents multiplets de la forme En pratique, la position utilisée pour le codage peut être définie au sein d’une séquence ordonnée des coefficients au sein du bloc ^ _^^,^^ concerné (ici le bloc ^ _^^,^^ Pour ce faire, pour le bloc concerné ^ _^^,^^ , on ordonne par exemple les coefficients du bloc par ordre de balayage de trame (couramment dénommé "raster scan order") : on ordonne les coefficients ^ _^^,^^ (1, ^ _^ ) par ordre croissant d’indice ^ _^ relatif au second paramètre spatial (^ _^ croissant de 1 à Ny), puis les coefficients ^ _^ ) par ordre décroissant d’indice ^ _^ (^ _^ décroissant de Ny à 1), puis à nouveau par ordre croissant d’indice ^ _^ pour les coefficients ^ _^^,^^ (3, ^ _^ ), etc. On propose ici de coder les positions des coefficients ^ _^^,^^ (^ _^ , ^ _^ ) non-nuls au moyen de données indicatives représentant la différence entre les positions de deux coefficients non-nuls successifs dans la séquence ordonnée de coefficients (sauf pour le premier coefficient dont la position sera par exemple codée en valeur absolue). Ainsi par exemple, les positions des coefficients non-nuls c1, c2, c3, c4 dans une séquence (0 c100 c2 c3000 c4000) seront codées par les chiffres : 2, 3, 1, 4, où 2 représente la position du premier coefficient non-nul c ₁ , 3 la différence entre la position du second coefficient non-nul c2 et la position du premier coefficient non-nul c1, 1 la différence entre la position du troisième coefficient non-nul c3 et la position du second coefficient non-nul c ₂ et 4 la différence entre la position du quatrième coefficient non-nul c4 et la position du troisième coefficient non-nul c3. Dans la suite, on note D11 les données indicatives des positions des coefficients non-nuls dans le bloc Le procédé initie alors une boucle à l’étape E152 de manière à traiter tous les blocs ^ _^^,^^ autres que le bloc Les différents blocs ^ _^^,^^ peuvent par exemple être traités dans un ordre de balayage de trame en considérant l’ensemble des blocs ^ _^^,^^ comme éléments d’une trame : on traite les blocs ^ _^,^^ par ordre d’indice ^ _^ croissant (de 2 à My dans l’exemple décrit ici), puis les blocs ^ _^,^^ par ordre d’indice ^ _^ décroissant (ici de My à 1), puis les blocs ^ _^,^^ par ordre d’indice ^ _^ croissant (ici de 1 à My), puis les blocs ^ _^,^^ par ordre d’indice ^ _^ croissant (ici de 1 à My), etc. Quel que soit l’ordre de traitement des blocs, on note dans la suite le bloc traité immédiatement avant le bloc courant ^ _^^,^^ (ce bloc pouvant être le bloc ^ _^^,^^ lors du traitement du premier bloc ^ _^^,^^ désigné par la boucle initiée à l’étape 152, ou un bloc précédemment traité à l’étape E152 décrite ci-dessous). Pour tous les blocs ^ _^^,^^ autres que le bloc ^ _^^,^^ (et ici dans l’ordre de traitement défini ci-dessus), le procédé effectue à l’étape E154 le codage des positions des coefficients non-nuls de ce bloc ^ _^^,^^ (bloc courant). Dans l’exemple décrit ici, on propose de coder d’une part (au sein de données D12) les positions des coefficients non-nuls dans le bloc précédent ^′ _^^,^^ pour lesquelles le coefficient correspondant dans le bloc courant ^ _^^,^^ est nul, et d’autre part (au sein de données D13) les positions des coefficients non-nuls dans le bloc courant ^ _^^,^^ pour lesquelles le coefficient correspondant dans le bloc précédent ^′ _^^,^^ est nul. Autrement dit, les données D12 représentent les positions des coefficients qui "disparaissent" de la représentation lorsque l’on passe du bloc précédent ^′ _^^,^^ au bloc courant ^ _^^,^^ , tandis que les données D13 représentent les positions des coefficients qui "apparaissent" dans la représentation lorsque l’on passe du bloc précédent au bloc courant ^ _^^,^^ . Ainsi, les données D13 désignent au moins une position associée à un coefficient non-nul dans le bloc courant ^ _^^,^^ , sans désigner en revanche au moins une autre position associée à un autre coefficient non-nul lorsque cette autre position est déjà indiquée par des données indicatives codées (dans le cadre du codage des positions dans le bloc précédent ^′ _^^,^^ ) pour une autre paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel. Les données D12 désignent quant à elles des positions correspondant (chacune) à un coefficient nul dans le bloc courant ^ _^^,^^ , mais pour lesquelles le coefficient correspondant était non-nul dans le bloc ^′ _^^,^^ précédemment traité, ces positions étant donc indiquées dans des données indicatives codées pour une autre paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel. Aucune donnée n’est donc codée dans le flux de données lorsque, pour une position donnée, les coefficients du bloc courant et du bloc précédent sont tous deux nuls, ou sont tous deux non-nuls. Les positions codées à l’étape E154 sont par exemple les positions absolues des coefficients concernés au sein de la séquence ordonnée (par ordre de balayage de trame) définie plus haut au sein de chaque bloc ^ _^^,^^ . Le procédé comprend alors une étape E156 de codage des coefficients (non- nuls) associés aux positions indiquées par les données indicatives D11, D12, D13. Ces coefficients sont par exemple codés dans l’ordre des traitements effectués aux étapes E150 et E154, c’est-à-dire en commençant par le bloc ^ _^^,^^ (ici le bloc ^ _^,^ ) puis par ordre de balayage de trame au sein des blocs et, au sein de chaque bloc, dans l’ordre de balayage de trame (décrit ci-dessus dans le cadre de l’étape E150). Pour le codage d’un coefficient W _k ( ) donné, on procède d’une part à la quantification de la partie réelle de ce coefficient et d’autre part à la quantification de la partie imaginaire du coefficient, la partie réelle quantifiée et la partie imaginaire quantifiée étant ensuite concaténée. Le quantificateur utilisé pour la quantification est par exemple un quantificateur scalaire uniforme ou un quantificateur tel que décrit dans l’article "Exact global motion compensation for holographic video compression", de Raees Kizhakkumkara Muhamad et al. in Applied Optics58, G204-G217 (2019). Dans l’exemple décrit ici, le codage des coefficients de l’étape E156 est réalisé postérieurement au codage des positions des étapes E150 et E154. En variante, le codage des coefficients pourrait être intercalé entre des étapes de codage des positions : par exemple, le codage des coefficients (non-nuls) d’un bloc ^ _^^,^^ pourrait être réalisé à la suite du codage des positions des coefficients non-nuls dans ce bloc ^ _^^,^^ tel que réalisé à l’étape E150 ou à l’étape E154. L’ensemble des données obtenues par les différentes étapes de codage mentionnées ci-dessus (étape E150, E154, E156) peuvent en outre être codées par un codeur entropique, par exemple un codeur arithmétique adaptatif, afin d’obtenir les données D1 à transmettre au dispositif de décodage 40 (ou en variante à mémoriser pour décodage ultérieur). Enfin, des marqueurs peuvent être insérés entre les différents types de données (ici entre les données D11, les données D12, les données D13 et les parties réelles et imaginaires de coefficients quantifiées) pour permettre de retrouver ces données au niveau du décodeur. D’autres solutions peuvent toutefois être utilisées, telles que l’utilisation d’une longueur fixe pour les données représentatives des parties réelles et des parties imaginaires des coefficients, par exemple. Le flux de données D1 transmis (du dispositif de codage 20 au dispositif de décodage 40 dans le contexte décrit plus haut) pour représenter l’hologramme numérique H _k comprend donc : - des données D11 indicatives des positions associées respectivement à des coefficients non-nuls au sein du bloc ^ _^^,^^ (c’est-à-dire parmi un ensemble des positions associées de manière prédéfinie aux coefficients associés respectivement aux différents multiplets définis par la valeur pour le premier paramètre fréquentiel et par la valeur pour le second paramètre fréquentiel) ; - puis successivement pour chacun des autres blocs ^ _^^,^^ , des données D12, D13 indicatives des positions associées respectivement à des coefficients non-nuls au sein du bloc ^ _^^,^^ concerné (c’est-à-dire parmi un ensemble des positions associées de manière prédéfinie aux coefficients associés respectivement aux différents multiplets définis par la valeur ^ _^^ pour le premier paramètre fréquentiel et par la valeur ^ _^^ pour le second paramètre fréquentiel) ; - pour chaque bloc, des données représentatives des coefficients associés respectivement aux positions indiquées par lesdites données indicatives relatives au bloc concerné. Comme déjà indiqué, pour chaque bloc, la position des coefficients utilisée pour le codage peut être définie au sein d’une séquence ordonnée des coefficients au sein du bloc ^ _^^,^^ concerné, par exemple l’ordre de balayage de trame déjà présenté. On décrit ci-dessous en référence à la figure 11 un procédé de décodage d’un tel flux de données, tel qu’il est mis en œuvre par exemple par le dispositif de décodage 40 à l’étape E30 décrite ci-dessus. Pour chaque bloc ^ _^^,^^ défini ci-dessus dans le cadre du procédé de codage de la figure 10, le dispositif de décodage (précisément la première unité de décodage 42 dans le contexte décrit en référence à la figure 6) décode à l’étape E160 les données indicatives des positions associées à des coefficients non-nuls dans le bloc ^ _^^,^^ concerné, c’est-à-dire parmi un ensemble des positions associées de manière prédéfinie (ici par définition du bloc) aux coefficients associés respectivement aux différents multiplets dont la première coordonnée fréquentielle est la valeur associée à l’indice ^ _^ et dont la seconde coordonnée fréquentielle est la valeur ^ _^^ associée à l’indice ^ _^ . Les différents blocs sont par exemple transmis (et donc présentés au sein du flux de données) dans un ordre prédéfini par convention entre le dispositif de codage 20 et le dispositif de décodage 40, tel que l’ordre déjà mentionné de balayage de trame au sein des blocs agencés sous forme d’une matrice de blocs. En variante, des données complémentaires pourraient être placées en en-tête des données relatives à un bloc ^ _^^,^^ et indiquer des paramètres désignant le bloc ^ _^^,^^ concerné (ces paramètres étant par exemple les indices ^ _^ et ^ _^ du bloc concerné). Dans le cas du bloc ^ _^^,^^ (ici le bloc ^ _^,^ ), les données D11 indicatives des positions des coefficients non-nuls (ici au sein d’une séquence ordonnée de coefficients du bloc représentent une différence entre une position associée à un coefficient non-nul et une autre position associée à un autre coefficient non-nul, comme exposé lors de la description du procédé de codage. Dans le cas des blocs ^ _^^,^^ autres que le bloc ^ _^^,^^ , les données D12, D13 indicatives des positions des coefficients non-nuls (ici au sein d’une séquence ordonnée de coefficients du bloc ^ _^^,^^ ) désignent soit les positions des coefficients nuls dans le bloc courant ^ _^^,^^ alors que le coefficient à la même position était non- nul dans le bloc ^′ _^^,^^ précédemment décodé (données D12), soit les positions des coefficients non-nuls dans le bloc courant ^ _^^,^^ alors que le coefficient à la même position était nul dans le bloc ^′ _^^,^^ précédemment décodé (données D13). Ainsi, les données D13 désignent au moins une position associée à un coefficient non-nul sans désigner au moins une autre position associée à un autre coefficient non-nul lorsque cette autre position a été indiquée par des données indicatives décodées en relation avec un autre bloc, c’est-à-dire pour une autre paire de valeurs respectives du premier paramètre fréquentiel et du second paramètre fréquentiel (cette paire de valeurs étant associé aux coefficients formant l’autre bloc). Lorsque les données indicatives des positions des coefficients non-nuls ont été décodées pour un bloc donné ^ _^^,^^ , le dispositif de décodage 40 (précisément ici la première unité de décodage 42) peut décoder les coefficients associés respectivement aux positions indiquées par ces données indicatives. Le dispositif de décodage 40 retrouve ainsi les coefficients non-nuls ^ _^^,^^ (^ _^ , ^ _^ ), c’est-à-dire les coefficients non-nuls Wk( (^ _^^ ^ _^^ , ^ _^^ , ^ _^^ ) ). En parcourant l’ensemble des blocs ^ _^^,^^ (pour lesquels les données sont successivement transmises dans le flux de données), le dispositif de décodage 40 obtient ainsi l’ensemble des coefficients Wk non-nuls de la représentation de l’hologramme numérique H _k . On remarque que, comme déjà indiqué dans le cadre du procédé de codage de la figure 10, le dispositif de décodage 40 peut recevoir toutes les données indicatives des positions des coefficients non-nuls pour l’ensemble de blocs ^ _^^,^^ , puis les coefficients associés à ces positions, ou, en variante, en considérant successivement les différents blocs ^ _^^,^^ , les données indicatives des positions des coefficients non- nuls de ce bloc et les coefficients associés à ces positions. On décrit à présent en référence à la figure 12 un exemple de dispositif de codage 60 conforme aux enseignements de l’invention. Il peut s’agir d’un dispositif de codage tel que le dispositif de codage 20 décrit plus haut en référence à la figure 5, auquel cas les modules décrits ci-dessous font partie de la première unité de codage 22 mentionnée plus haut. Le dispositif de codage 60 est par exemple réalisé sous forme d’un dispositif électronique comprenant un processeur (tel qu’un microprocesseur) ; dans ce cas, chacun des modules décrits ci-dessous peut être mis en œuvre du fait de l’exécution, par ce processeur, d’instructions de programme d’ordinateur mémorisées au sein du dispositif de codage (par exemple dans une mémoire 70 de ce dispositif de codage). Toutefois, en variante, au moins un des modules décrits ci-dessous peut être réalisé au moyen d’un circuit électronique dédié, par exemple un circuit intégré à application spécifique. Le dispositif de codage 60 comprend un module de traitement de bloc initial 62, un module de traitement de bloc ultérieur 64, un module de traitement des coefficients 66 et un module de finalisation 68. Le dispositif de codage 60 est ainsi apte à mettre en œuvre le procédé décrit ci- dessus en référence à la figure 10, ici lorsque le processeur du dispositif de codage 60 exécute certaines des instructions de programme d’ordinateur mémorisées dans la mémoire 70 du dispositif de codage et déjà mentionnées ci-dessus. Comme expliqué plus haut à propos de la figure 10, les coefficients à coder par le dispositif de codage 60 sont des coefficients respectivement associés à des multiplets formés chacun par un premier paramètre spatial ^ _^^ , un second paramètre spatial ^ _^^ , un premier paramètre fréquentiel et un second paramètre fréquentiel ^ _^^ . Comme déjà indiqué, on décrit ici le cas où ces coefficients sont des coefficients W _k qui forment (à eux seuls) une représentation d’un hologramme numérique. En variante, ces coefficients pourraient représenter un hologramme numérique en combinaison avec d’autres éléments, par exemple dans le cas où les coefficients définissent un résidu, tel que le résidu e _k+1 mentionné plus haut. Les modules 62, 64, 66, 68 peuvent accéder à la mémoire 70 du dispositif de codage 60 où sont mémorisés les coefficients Wk à coder et peuvent donc traiter ces coefficients sous forme de blocs ^ _^^,^^ , ^ _^^,^^ , comme expliqué à présent. Le module de traitement de bloc initial 62 est conçu pour coder des données D11 indicatives des positions des coefficients non-nuls au sein du bloc ^ _^^,^^ de coefficients W _k défini plus haut, par codage de différence entre deux positions successives associées à des coefficients non-nuls. Ce bloc ^ _^^,^^ regroupe les coefficients associés aux multiplets dont le premier paramètre fréquentiel est égal à et donc le second paramètre fréquentiel est égal Le module de traitement de bloc ultérieur 64 est conçu pour coder des données D12, D13 indicatives des positions des coefficients non-nuls au sein de chacun des blocs ^ _^^,^^ de coefficients W _k définis plus haut, par codage des positions des coefficients non-nuls dans le bloc précédent ^′ _^^,^^ pour lesquelles le coefficient correspondant dans le bloc courant ^ _^^,^^ est nul (données D12), et par codage des positions des coefficients non-nuls dans le bloc courant ^ _^^,^^ pour lesquelles le coefficient correspondant dans le bloc précédent ^′ _^^,^^ est nul (données D13). Comme expliqué plus haut, pour chaque couple ( ^ _^ , ^ _^ ), le bloc ^ _^^,^^ regroupe les coefficients Wk associés aux multiplets dont le premier paramètre fréquentiel est égal et donc le second paramètre fréquentiel est égal à ^ _^^ . Le module de traitement des coefficients 66 est conçu pour coder les coefficients W _k non-nuls, c’est-à-dire les coefficients W _k dont les positions sont indiquées dans les données indicatives D11, D12, D13, selon un ordre prédéfini, par exemple par ordre de balayage de trame au sein de tous les blocs ^ _^^,^^ , et, au sein de chaque bloc, par ordre de balayage de trame parmi les coefficients Wk du bloc. On note D14 les coefficients codés. Le module de finalisation 68 procède au codage entropique des données D11, D12, D13, D14 et/ou insère des marqueurs entre les différents types de données afin d’obtenir les données D1 à émettre. On décrit à présent en référence à la figure 13 un exemple de dispositif de décodage 80 conforme aux enseignements de l’invention. Il peut s’agir d’un dispositif de décodage tel que le dispositif de décodage 40 décrit plus haut en référence à la figure 6, auquel cas les modules décrits ci-dessous font partie de la première unité de décodage 42 mentionnée plus haut. Le dispositif de décodage 80 est par exemple réalisé sous forme d’un dispositif électronique comprenant un processeur (tel qu’un microprocesseur) ; dans ce cas, chacun des modules décrits ci-dessous peut être mis en œuvre du fait de l’exécution, par ce processeur, d’instructions de programme d’ordinateur mémorisées au sein du dispositif de décodage (par exemple dans une mémoire de ce dispositif de décodage). Toutefois, en variante, au moins un des modules décrits ci-dessous peut être réalisé au moyen d’un circuit électronique dédié, par exemple un circuit intégré à application spécifique. Le dispositif de décodage 80 comprend un module d’analyse 82, un module de construction de bloc initial 82, un module de construction de bloc ultérieur 84 et un module de décodage des coefficients 88. Le dispositif de décodage 80 est ainsi apte à mettre en œuvre le procédé décrit ci-dessus en référence à la figure 11, ici lorsque le processeur du dispositif de décodage 80 exécute certaines des instructions de programme d’ordinateur mémorisées dans la mémoire du dispositif de décodage et déjà mentionnées ci- dessus. Le module d’analyse 82 reçoit les données D1 et y retrouve les différents données D11, D12, D13, D14, par exemple à l’aide des marqueurs susmentionnés. Le module d’analyse 82 peut éventuellement effectuer une étape de décodage entropique (lorsqu’un codage entropique a été effectué au sein du module de finalisation 68 décrit ci-dessus). Le module de construction de bloc initial 84 est conçu pour décoder les données D11 représentatives des différences entre deux positions successives associées à des coefficients non-nuls au sein du bloc de coefficients Wk, de manière à obtenir les positions des coefficients non-nuls au sein de ce bloc ^ _^^,^^ . On rappelle à cet égard que les coefficients du bloc ^ _^^,^^ sont ordonnés de manière prédéfinie, par exemple un ordre de balayage de trame. Diverses solutions sont envisageables pour mémoriser les positions de ces coefficients non-nuls. Par exemple, le module de construction de bloc initial 84 initialise à la valeur zéro les éléments d’une table correspondant aux coefficients de la forme Wk(x,y, ^ _^ , ^ _^ ), puis met à une valeur non nulle prédéfinie (par exemple la valeur 1) les éléments correspondant aux coefficients indiqués comme non-nuls d’après les données D11. (Lorsque les coefficients sont des nombres complexes, on peut procéder ainsi d’une part pour la partie réelle des coefficients et d’autre part pour la partie imaginaire des coefficients.) Pour chaque couple (^ _^ , ^ _^ ) (différent de (^ _^ , ^ _^ )), le module de construction de bloc ultérieur 86 est conçu pour obtenir les positions des coefficients Wk non-nuls au sein de ce bloc ^ _^^,^^ par décodage d’une part des données D12 indicatives des positions des coefficients non-nuls dans le bloc précédemment traité ^′ _^^,^^ pour lesquelles le coefficient correspondant dans le bloc courant ^ _^^,^^ est nul et des données D13 indicatives des positions des coefficients non-nuls dans le bloc courant ^ _^^,^^ pour lesquelles le coefficient correspondant dans le bloc précédemment traité ^′ _^^,^^ est nul. Ici encore, diverses solutions sont envisageables pour mémoriser les positions de ces coefficients non-nuls. Par exemple, le module de construction de bloc ultérieur 86 initialise à la valeur zéro les éléments d’une table correspondant aux coefficients de la forme Wk(x,y, ^ _^ , ^ _^ ), puis met à une valeur non nulle prédéfinie (par exemple la valeur 1) les éléments correspondant aux coefficients indiqués comme non-nuls d’après les données D12 et D13. (Lorsque les coefficients sont des nombres complexes, on peut procéder ainsi d’une part pour la partie réelle des coefficients et d’autre part pour la partie imaginaire des coefficients.) Le module de décodage des coefficients 88 est conçu pour décoder les données D14 représentant les coefficients W _k (non-nuls) associés respectivement aux positions indiquées par les données D11, D12, D13. Comme indiqué ci-dessous, les coefficients Wk sont codés au sein des données D14 dans un ordre prédéfini, de sorte que le module de décodage des coefficients 88 peut associer chaque coefficient décodé à une position d’un coefficient non nul, et ainsi retrouver la valeur d’un coefficient Wk( ). Par exemple, dans le cadre de la possibilité de mise en œuvre décrite ci-dessous, le module de décodage des coefficients 88 parcourt les tables susmentionnées dans un ordre prédéfini correspondant à l’ordre de codage des coefficients, puis, lorsqu’un coefficient est décodé, la valeur du coefficient décodé est attribuée au prochain élément de valeur prédéfinie non-nulle dans la table en cours de traitement.