FR 3036197 A1 20161118 FR 1554363 A 20150513 1 Procédé de génération numérique d'un hologramme, dispositif et programme d'ordinateur associés 1. Domaine de l'invention Le domaine de l'invention est celui de la génération numérique d'un champ lumineux d'hologramme. 2. Présentation de l'art antérieur L'holographie est souvent considérée comme la technologie de visualisation 3D la plus prometteuse puisqu'elle elle fournit une illusion du relief la plus naturelle et la plus réaliste possible, et ceci sans la nécessité de porter de lunettes et sans causer de fatigue visuelle. En effet, cette technologie permet de fournir tous les indices de perception de la profondeur du système visuel humain. Ces dernières années, plusieurs méthodes de génération numérique d'hologrammes ont été proposées. Grâce à ces méthodes, il est possible d'obtenir l'hologramme d'une scène réelle ou synthétique en simulant la propagation de la lumière provenant de la scène vers le plan de l'hologramme. On connait du document de génération de J. S. Underkoffler, intitulé « Occlusion Processing and smooth surface shading for fully connputed synthetic holography », publié dans les proceedings de la conférence « Practical Holography XI and Holographic Materials II, en avril 1997, Proc. SPIE 3011' pp 19-30, 00048, une technique de génération numérique d'hologrammes consistant à décomposer la scène en un nuage de points lumineux et à calculer l'onde lumineuse issue de la scène comme la somme des ondes sphériques émises par chacun des points. Dans la suite de cet article, nous appellerons cette approche "basée points". Cette approche basée points est très flexible n'impose aucune restriction sur la géométrie de la scène. Cependant, la complexité algorithmique de cette approche est très grande puisqu'elle nécessite un calcul par point de la scène et par pixel de l'hologramme. De plus, pour que les surfaces apparaissent continues, la scène doit être échantillonnée à l'aide d'un très grand nombre de points, rendant la génération de l'hologramme beaucoup trop lente. Afin de réduire la complexité algorithmique de cette approche, plusieurs méthodes ont été proposées, basées sur l'utilisation de tables de correspondance pré-calculées, la réduction des redondances spatiales, l'utilisation de formules de récurrence, l'utilisation de plans intermédiaires, ou encore le recours à des cartes graphiques ou à du matériel spécifique. De plus, il est possible de réduire le nombre de tests de visibilité nécessaires pour la 3036197 2 gestion des occultations en groupant les pixels de l'hologramme ou encore en groupant les points de la scène. On connait en outre du document de K. Matsushinna, intitulé « Exact hidden-surface removal in digitally synthetic full-parallax holograms », publié dans les Proceedings de la 5 conférence « Practical holography XIX : Material and applications », qui s'est tenue à San Jose en Californie, en 2005, une autre approche consistant à décomposer la scène en un ensemble de plans parallèles au plan de l'hologramme et à calculer la propagation de l'onde lumineuse issue de la scène d'un plan à un autre en utilisant des formules de propagation de la lumière telles que la formule du spectre angulaire. Dans la suite, nous appellerons cette approche "basée champ". Le 10 calcul du spectre angulaire nécessitant deux Transformées de Fourier Rapide (FFT), sa complexité algorithmique est plus importante que celle du calcul de l'onde sphérique issue d'un point. En revanche, l'onde émise par chacun des points compris dans un même plan est calculée en une seule fois grâce au spectre angulaire. Ainsi, cette approche est plus efficace que l'approche basée points lorsque les objets de la scène sont constitués de grands plans parallèles au plan de 15 l'hologramme. En revanche, lorsque la scène comprend des formes plus complexes, un nombre important de plans contenant peu de points est nécessaire pour l'échantillonner, rendant cette approche moins efficace que l'approche basée points. Quelle que soit la méthode utilisée, une contrainte importante en holographie est la capacité à restituer les occultations entre les objets de la scène. En effet, les occultations dans la 20 scène constituent l'un des indices de perception de la profondeur les plus importants. La prise en compte des occultations en génération numérique d'hologrammes est l'équivalent de la détermination des surfaces cachées en synthèse d'image. Cependant, contrairement aux images synthétiques, un hologramme généré numériquement fournit la parallaxe de mouvement, qui est liée aux occultations. Ainsi, la prise en compte des occultations en génération numérique 25 d'hologrammes implique que les objets visibles changent en fonction du point de vue de l'utilisateur. Les techniques de gestion des occultations dépendent de l'approche utilisée. Dans l'approche basée points, les occultations entre objets de la scène sont habituellement prises en compte en effectuant un test de visibilité pour vérifier la présence ou non d'obstacles entre chaque point de la scène et chaque pixel de l'hologramme. 30 Dans l'approche basée champ, les occultations entre objets de la scène sont prises en compte à chaque étape de diffraction en multipliant le champ lumineux incident sur un plan donné par le masque binaire de ce plan, supprimant ainsi la nécessité d'effectuer de test de visibilité. 3036197 3 3. Inconvénients de l'art antérieur La méthode basée points est très calculatoire et la méthode basée champ n'est pas adaptée aux scènes réelles comprenant des formes complexes. 5 4. Objectifs de l'invention L'invention vient améliorer la situation. L'invention a notamment pour objectif de pallier ces inconvénients de l'art antérieur. Plus précisément, un objectif de l'invention est de proposer une solution qui permette de générer un hologramme à partir d'une scène réelle pour un coût moindre en termes de ressources 10 de calcul que les techniques existantes. Un autre objectif de l'invention est de proposer une solution qui soit aussi capable de gérer efficacement les occultations de la scène. 5. Exposé de l'invention Ces objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints à l'aide d'un 15 procédé de génération numérique d'un hologramme dans un plan, dit plan de l'hologramme, à partir d'une scène tridimensionnelle, comprenant une étape de découpage de ladite scène en une pluralité de plans parallèles au plan de l'hologramme. Le procédé selon l'invention est particulier en ce qu'il comprend les étapes suivantes, mises en oeuvre pour un plan de la pluralité de plans, dit plan courant, les plans étant parcourus du plus 20 éloigné au plus proche du plan de l'hologramme : - comptage d'un nombre de points de la scène d'amplitude non nulle dans le plan courant ; - choix d'une première ou d'une deuxième technique de propagation d'une onde lumineuse émise par le plan courant en fonction d'un nombre de points d'amplitude non nulle de la portion de scène comprise dans le plan courant et d'une valeur de seuil prédéterminée, la 25 première technique, dite basée points, destinées à calculer la propagation d'une somme d'ondes lumineuses émises par des sources ponctuelles constituées par les points de la portion de scène d'amplitude non nulle du plan courant sur un plan suivant, la deuxième technique, dite basée champ, destinée à calculer globalement la propagation d'un champ lumineux émis par la portion de scène comprise dans le plan courant sur un plan suivant ; 3036197 4 - traitement du plan courant en fonction de la technique de propagation choisie, destiné à calculer une onde lumineuse émise par le plan courant. Selon l'invention, lesdites étapes sont répétées pour la pluralité de plans et le procédé comprend en outre une étape d'obtention de l'hologramme à partir des ondes lumineuses calculées pour la 5 pluralité de plans. Alors que les méthodes de l'état de l'art se basent sur une décomposition de la scène en une pluralité de plans qu'elles traitent soit globalement selon un technique basée champ ou soit localement selon une technique basée points, l'invention utilise les deux approches conjointement et tire ainsi parti des avantages de chacune d'entre elles. 10 Tout d'abord, la scène est décomposée en un ensemble de plans parallèles au plan de l'hologramme. Ensuite, l'onde lumineuse provenant de la scène est propagée d'un plan à un autre en partant du plan le plus distant. Selon le nombre de points de la scène compris dans un plan l'onde lumineuse émise par les points de la scène est calculée soit en considérant chaque point indépendamment, soit en considérant le plan dans son ensemble. Enfin, l'onde lumineuse est 15 propagée jusqu'au plan de l'hologramme afin d'obtenir l'hologramme final. Avec l'invention, on choisit la méthode la plus adaptée à la scène. Les plans comprenant un faible nombre de points d'amplitude non nulle sont avantageusement traités par la méthode basée points, pour un coût raisonnable, tandis que les plans comprenant un nombre de points supérieur à la valeur d'un seuil prédéterminé sont traités globalement par la méthode basée champ. 20 Selon un aspect de l'invention, ladite étape de traitement du plan courant comprend les sous- étapes suivantes : - identification du dernier plan traité selon la deuxième technique, dit dernier plan basé champ; 25 - calcul de l'onde lumineuse émise par la portion de scène comprise dans le plan courant ; - lorsque la deuxième technique a été choisie pour le plan courant, dit plan basé champ, l'étape de traitement comprend les sous-étapes suivantes : o calcul d'une propagation de l'onde lumineuse émise par au moins le dernier plan basé champ sur le plan courant, selon la technique choisie pour ledit plan; 3036197 5 o calcul de l'onde lumineuse émise par le plan courant par sommation de l'onde lumineuse émise par le plan courant et au moins la propagation de l'onde lumineuse émise par au moins le dernier plan basé champ ; et l'étape d'obtention de l'hologramme comprend le calcul de la propagation de l'onde lumineuse 5 émise par le plan basé champ le plus proche de l'hologramme sur le plan de l'hologramme. Avantageusement, l'invention traite la propagation des ondes lumineuses de proche en proche, d'un plan basé champ à l'autre. L'hologramme est obtenu par simple propagation du champ lumineux calculé par le plan basé champ le plus proche du plan de l'hologramme, selon la deuxième technique. 10 Selon encore un autre aspect de l'invention, lorsque ledit au moins un plan précédent comprend en outre au moins un plan basé points, intercalé entre le dernier plan basé champ et le plan courant, le calcul de la propagation de l'onde lumineuse émise par au moins le dernier plan basé champ sur le plan courant comprend en outre le calcul d'une propagation de l'onde lumineuse émise par les points d'amplitude non nulle de la portion de scène comprise dans ledit au moins un 15 plan intercalaire sur le plan courant selon la première technique. Les ondes lumineuses émises par des plans basés points intercalaires sont elles aussi propagées directement jusqu'au plan basé champ courant. De cette manière, l'invention ne gère que la propagation d'une onde lumineuse issue d'un plan basé points sur un plan basé champ, dont la complexité dépend du nombre de points d'amplitude non nulle initialement présents dans la 20 portion de scène comprise dans le plan basé points. Un avantage est de maîtriser la complexité calculatoire de la méthode proposée. Selon une caractéristique avantageuse de l'invention, l'étape de traitement du plan courant comprend une sous-étape de gestion d'une occultation de l'onde lumineuse propagée par au moins le dernier plan basé champ dans le plan courant par application d'un masque binaire occultant à 25 ladite onde lumineuse, ledit masque prenant une valeur nulle aux points de la portion de scène comprise dans le plan courant et une valeur non nulle en dehors de cette portion et en ce que la sous-étape de calcul de l'onde lumineuse totale émise par le plan courant comprend une sommation de l'onde lumineuse émise par le plan courant à l'onde lumineuse obtenue après application du masque binaire occultant. 30 De cette manière, pour un plan courant basé champ, on ne propage pas l'onde lumineuse issue des plans précédents aux points non nuls de la portion de scène comprise dans le plan courant. 3036197 6 Selon encore un autre aspect de l'invention, lorsque la première technique a été choisie pour traiter le plan courant, dit plan dit basé points, l'étape de traitement comprend en outre les sous-étapes suivantes : - calcul d'une onde lumineuse propagée par au moins le dernier plan basé champ selon 5 la première technique et bloquée par les points d'amplitude non nulle de la portion de scène comprise dans le plan courant ; - calcul de l'onde lumineuse émise par le plan courant par soustraction de l'onde lumineuse bloquée par les points d'amplitude non nulle de la portion de scène comprise dans le plan courant à l'onde lumineuse émise par lesdits points. 10 L'invention permet de gérer les occultations dues aux plans basés point intercalaires, tout en maîtrisant la complexité calculatoire de la méthode. En effet, elle évite de créer de nouveaux points dans les plans basés points en propageant le champ lumineux d'un plan précédent basé point ou les ondes lumineuses d'un plan précédent basé point, de manière à ne pas illuminer de points supplémentaires dans le plan basé points courant et à ne pas augmenter le temps de calcul dû au 15 traitement du plan basé points courant. A la place, elle propose de calculer au niveau de chaque plan intercalaire basé points, l'amplitude de l'onde lumineuse bloquée par la scène dans ce plan et de la propager jusqu'au plan basé champ suivant. De cette manière, elle transmet au plan basé champ suivant l'information relative à la quantité d'onde lumineuse issue des plans précédents qui est occultée par ce plan intermédiaire. L'amplitude de l'onde lumineuse bloquée est ensuite 20 soustraite à celle de l'onde lumineuse issue des plans précédents qui avait été propagée directement vers le plan basé champ suivant, comme s'il n'y avait pas d'occultation par les plans intercalaires. Ainsi, l'invention propose une solution astucieuse au problème de la gestion des occultations entre deux plans de nature différente, qui, pour un plan donnée, ne remet pas en cause le principe de choix d'une technique de propagation adaptée à la complexité de la portion de 25 scène comprise dans ce plan. Selon encore un autre aspect de l'invention, la sous-étape de calcul de l'onde lumineuse bloquée comprend le calcul des ondes lumineuses propagées au moins par le dernier plan basé champ selon la première technique et l'application d'un masque binaire occultant inversé à l'onde 30 lumineuse calculée, ledit masque prenant une valeur non nulle aux points de la portion de scène comprise dans le plan courant et une valeur nulle en dehors de cette portion. 3036197 7 Un premier avantage de cette technique de masque inversé est qu'elle permet d'obtenir simplement l'amplitude des ondes lumineuses bloquées par le plan basé points courant. Un deuxième avantage est qu'il permet de limiter le calcul de l'onde bloquée aux points d'amplitude non nulle de la scène, puisqu'ailleurs le masque binaire inversé prend une valeur nulle. 5 Selon un autre aspect de l'invention, la sous-étape de calcul de l'onde lumineuse bloquée comprend un découpage en une pluralités de zones du plan courant, une détermination des points d'un plan précédent à partir du dernier plan basé champ, contribuant à la propagation d'une onde lumineuse dans une zone de la pluralité de zones comprenant des points d'amplitude non nulle de la portion de scène comprise dans le plan courant, le calcul de l'onde lumineuse bloquée n'étant 10 effectué que pour les points déterminés. Avantageusement, on limite le calcul aux points du plan précédents qui contribuent à illuminer les points non nuls de la scène comprise dans le plan courant. Selon un aspect de l'invention, le procédé comprend une étape d'évaluation du seuil au moins en fonction d'un temps de traitement estimé pour la première technique de propagation, 15 d'un temps de traitement estimé pour la deuxième technique de propagation pour le plan courant, et d'un nombre de points d'amplitude non nulle des portions de scène comprises dans les plans précédents. Avantageusement, le seuil est évalué spécifiquement pour chaque plan de la scène, ce qui garantit de choisir, pour ce plan, la méthode la moins coûteuse en temps de calcul. Les temps de 20 traitement estimés pour chacune des techniques de propagation peuvent aussi prendre en compte la gestion des occultations selon l'invention. Le procédé qui vient d'être décrit dans ses différents modes de réalisation est avantageusement mis en oeuvre par un dispositif de génération numérique d'un hologramme dans un plan, dit plan de l'hologramme, à partir d'une scène tridimensionnelle, comprenant une unité de 25 découpage de ladite scène en une pluralité de plans parallèles au plan de l'hologramme. Un tel dispositif comprend les unités suivantes, aptes à être mises en oeuvre pour un plan de la pluralité de plans, dit plan courant, les plans étant parcourus du plus éloigné au plus proche du plan de l'hologramme : - comptage d'un nombre de points de la scène d'amplitude non nulle dans le plan courant ; 30 - choix d'une première ou d'une deuxième technique de propagation d'une onde lumineuse émise par le plan courant en fonction d'un nombre de points d'amplitude non nulle de la 3036197 8 portion de scène comprise dans le plan courant et d'une valeur de seuil prédéterminée, la première technique, dite basée points, destinée à calculer la propagation d'une somme d'ondes lumineuses émises par des sources ponctuelles constituées par les points de la portion de scène d'amplitude non nulle du plan courant sur un plan suivant, la deuxième 5 technique, dite basée champ, destinée à calculer globalement la propagation d'une onde lumineuse émise par la portion de scène comprise dans le plan courant sur un plan suivant ; - traitement du plan courant selon la technique de propagation choisie, destiné à calculer une onde lumineuse émise par le plan courant. 10 Selon l'invention, la mise en oeuvre desdites unités est répétée pour la pluralité de plans et le dispositif comprend en outre une unité d'obtention de l'hologramme à partir des ondes lumineuses calculées pour la pluralité de plans. Corrélativement, l'invention concerne aussi un terminal d'utilisateur comprenant un module de commande d'un dispositif de restitution d'un hologramme. Ledit terminal comprend en outre un 15 dispositif de génération numérique dudit hologramme selon l'invention. L'invention concerne encore un programme d'ordinateur comportant des instructions pour la mise en oeuvre des étapes d'un procédé de génération d'un hologramme tel que décrit précédemment, lorsque ce programme est exécuté par un processeur. Ce programme peut utiliser n'importe quel langage de programmation. Il peut être 20 téléchargé depuis un réseau de communication et/ou enregistrés sur un support lisible par ordinateur. L'invention se rapporte enfin à un support d'enregistrement, lisible par un processeur, intégré ou non au dispositif de génération numérique d'un hologramme, éventuellement amovible, mémorisant respectivement un programme d'ordinateur mettant en oeuvre un procédé de 25 génération numérique d'un hologramme, tels que décrit précédemment. 6. Liste des figures D'autres avantages et caractéristiques de l'invention apparaîtront plus clairement à la lecture de la description suivante d'un mode de réalisation particulier de l'invention, donné à titre 30 de simple exemple illustratif et non limitatif, et des dessins annexés, parmi lesquels : 3036197 9 - la figure 1 illustre un exemple de scène tridimensionnelle découpée en une pluralité de plans parallèles au plan de l'hologramme ; - la figure 2 présente de façon schématique les étapes du procédé de génération numérique d'un hologramme selon un mode de réalisation de l'invention ; 5 - la figure 3 présente de façon schématique le découpage d'une scène tridimensionnelle en une pluralité de plans parallèles au plan de l'hologramme, et la propagation des ondes lumineuses entre deux plans basés champ successifs selon un mode de réalisation de l'invention ; - la figure 4 détaille l'étape de traitement d'un plan courant selon un mode de réalisation de 10 l'invention ; - la figure 5 illustre de façon schématique une région de contribution d'un point lumineux prise en compte par la technique de propagation d'une onde lumineuse basée points ; la figure 6 illustre de façon schématique le traitement des ondes lumineuses propagées par les plans précédents sur un plan basé champ selon l'invention ; 15 - la figure 7 présente de façon schématique un exemple de masque occultant mis en oeuvre par l'invention pour la portion de scène comprise dans un plan basé champ ; - la figure 8 présente de façon schématique le masque occultant inverse mis en oeuvre par l'invention pour un plan basé points, selon un mode de réalisation de l'invention ; - la figure 9 illustre de façon schématique le principe d'un mode de réalisation de l'invention 20 permettant de simplifier le calcul de l'occultation des ondes lumineuses la portion de scène comprise dans un plan basé points; - la figure 10 présente de façon schématique la structure matérielle d'un dispositif de génération numérique d'un champ d'hologramme selon l'invention ; et - la figure 11 illustre un autre exemple d'usage de l'invention. 25 7. Description d'un mode de réalisation particulier de l'invention Le principe général de l'invention repose sur le découpage d'une scène tridimensionnelle en une pluralité de plans parallèles au plan de l'hologramme et sur une solution hybride de génération numérique d'hologramme, mettant en oeuvre une première technique basée points pour les plans 30 comprenant un nombre limité de points de la scène d'amplitude non nulle et une deuxième 3036197 10 technique basée champ pour les plans comprenant un nombre plus important de points de la scène. En relation avec la Figure 1, on considère la géométrie d'une scène tridimensionnelle Sc comprenant une bouilloire, un plan PH, dit plan de l'hologramme sur lequel un hologramme 5 représentatif de la scène Sc doit être généré numériquement et une pluralité de plans Po, P - avec Nzentier non nul supérieur ou égal à 2, parallèles au plan de l'hologramme PH. On considère en outre un référentiel (O, x, y, z), tel que l'hologramme soit inscrit dans le plan (O, x, y) et que l'axe Oz soit perpendiculaire à PH et à la pluralité de plans P. Par exemple, les plans Pi sont situés à égale distance d, les uns des autres. Po est le plan le plus 10 éloigné du plan de l'hologramme, à la distance zmax de PH et PNz_i est le plan le plus proche du plan de l'hologramme, à la distance zff',, de PH, avec zff',, Mimriax, on décide de mettre en oeuvre la deuxième méthode, dite basée champ pour traiter le plan courant Pi. Le plan Pi est alors désigné par plan « basé champ ». En relation avec la Figure 3, on présente un exemple de découpage en plans successifsPo à PNz_i. Les plans P0, P3 et P, sont des plans basés champ. Les plans P1, p P p et P6 sont des -1/ - 2/ - 4, - 5 - 6 10 plans dits « basés points ». En E4, le plan courant Pi est traité en fonction de la technique de propagation choisie. Cette étape fournit l'onde lumineuse totale ui émise par ce plan. De façon avantageuse, elle est stockée dans une mémoire Ml. En E5, on teste si le plan courant est le dernier plan, c'est-à-dire si i = Ar, -1. Si ce n'est 15 pas le cas, on incrémente i et on répète les étapes précédemment décrites pour le plan suivant. Si i = Ar, -1, on passe à l'étape E7. Au cours de cette étape, on calcule la propagation de l'onde lumineuse UN z_1 émise par le dernier plan P',z_1 sur le plan de l'hologramme, selon la technique de propagation basée champ et on obtient l'hologramme H. L'étape de traitement E4 va maintenant être détaillée en relation avec la Figure 4, pour un 20 plan basé champ, par exemple le plan P3, puis pour un plan basé points, par exemple P2. En E41 on identifie parmi les plans précédents, déjà traités, le plan basé champ le plus proche ou dernier plan basé champ. Dans l'exemple de la Figure 3, il s'agit du plan P0. Dans la suite, on lui affecte l'indice I, avec I entier compris entre 0 et Nz-1. En E42, on calcule l'onde lumineuse émise par la portion de scène comprise dans le plan P. 25 Chaque plan Pi fonctionne comme une source de lumière qui émet une onde complexe oi donnée par : oi (x, y) = ai (x, y)exp [j (X, y)] (2) où ai(x, y) est l'amplitude du point aux coordonnées (x, y) dans le plan Pi et cpi(x, y) est sa phase, initialisée aléatoirement afin de restituer une scène diffuse. 3036197 12 L'amplitude ai(x,y) est calculée en utilisant un modèle d'illumination, par exemple celui décrit par Phong dans l'article « Illumination for Computer Generated Pictures » publié dans le volume 18, issue 6 du journal « Communications of the ACM », pages 311-317, en 1975. En E44, on calcule la propagation des ondes lumineuses complexes émises par les plans précédents 5 à partir du plan basé champ précédent PI à Pi-1, à savoir dans l'exemple de la Figure 3 P0, P1 et P2, jusqu'au plan courant Pi selon la technique de propagation T1 ou T2 choisie pour chaque plan précédent. Pour le plan PI = Po basé champ, la propagation de la lumière jusqu'au plan Pi est calculée à l'aide de la formule du spectre angulaire connue de l'homme de métier et décrite par Goodman dans le 10 livre intitulé « Introduction to Fourier Optics », pages 55 à 58 et édité par Roberts and Company Publishers en 2005. Le spectre angulaire exprime la propagation d'une onde lumineuse entre deux plans parallèles séparés d'une distance z par : PZ {o}(x, y) = F-1 {r {o} exp (-j2n-V2.-2 - fX - f z)) (x, y) (3) 15 où À est la longueur d'onde de la lumière, f' et fy sont les fréquences spatiales, et .7 et.-1 sont respectivement la Transformée de Fourier et la Transformée de Fourier Inverse. Ces transformées peuvent être calculées en utilisant un algorithme de Transformée de Fourier Rapide de type FFT (pour «Fast Fourier Transform » en anglais. La complexité algorithmique de l'opérateur Pw est donc de O(N2log(N)). 20 Dans la suite, on désignera par Pzie, l'opérateur de propagation d'une onde lumineuse selon la technique basée champs sur une distance z. La propagation de l'onde lumineuse u, =0, du plan Po sur le plan Pi peut donc s'exprimer comme suit : 25 uW (x, y) = Pzw.._idztuol (x, y) = F-1 fFtitol exp (- j2n-12.-2 - - fy2 (zmax - iclz))} (x, y) (4) On considère maintenant la propagation des ondes lumineuses émises par les plans intercalaires Pi avec j entier compris entre 1+1 et i-1. Dans l'exemple de la Figure 3, les plans P1 et P2 sont basés points. On désigne par u1 l'onde lumineuse émise par le plan intercalaire P. L'onde lumineuse up u2 est obtenue en mettant en oeuvre l'étape de traitement E4 selon l'invention dans le cas d'un plan 30 « basé points ». Elle sera détaillée dans la suite de la description. La propagation de l'onde lumineuse u1 du plan P1 basé points sur le plan Pi, i = 3 est donc calculée en utilisant la première technique dite basé points. 3036197 13 La propagation de l'onde lumineuse uj d'un plan Pj basé points sur le plan Pi est donc calculée en utilisant la première technique dite basé points. Suivant l'approche basée points, les points de la scène situés dans chaque plan sont considérés 5 comme des sources lumineuses sphériques et la propagation de la lumière provenant de ce plan jusqu'au plan Pi est calculée comme la somme des ondes sphériques émises par chacun des points du plan. L'onde sphérique émise par un point lumineux k de coordonnées (xk, Yk, z) est donnée par une formule de spectre angulaire connue de l'homme de métier et par exemple décrite dans le livre de Goodman, intitulé « Introduction to Fourier Optics », pages 57 à 61 et 10 édité par Roberts and Company Publishers en 2005, de la façon suivante : ( 1 -j2n\IÀ.-2-n-f,z} Wk (X, 37) = 0 (Xk, Y k)F e 0 8(x - xk, y - Yk) (5) où o(xk, yk) est l'amplitude complexe du point et est l'opérateur de convolution. La convolution d'une fonction avec une impulsion de Dirac la recentre autour de l'impulsion. Ainsi, si l'on connait le terme de la transformée de Fourier dans l'équation (5) à l'avance, wk peut être calculée 15 simplement en factorisant ce terme avec l'amplitude du point, puis en la décalant spatialement. Afin d'accélérer le calcul, il est possible d'utiliser une table de correspondance pré-calculée T donnée par h(x, y, z) (6) (x, y, z) = F-1 région de contribution d'un point donné, h est une fonction fenêtre permettant (T e--127T \I,1-2-n-f,z de restreindre la 20 égale à un dans la région de contribution du point et à zéro ailleurs. Cette fonction limite les fréquences spatiales de l'onde complexe afin d'éviter le recouvrement de spectre dans l'hologramme. Suivant le critère de Nyquist, la fréquence spatiale maximale pouvant être représentée avec un pas d'échantillonnage de p est donnée par furax 25 En relation avec la Figure 5, on considère maintenant la région de contribution d'un point lumineux (x,y,z) sur un plan parallèle au plan de l'hologramme, par exemple le plan Pi courant. L'équation des réseaux de diffraction, connue de l'homme de métier, et énoncée par exemple dans le document de Goodman déjà cité, donne la relation entre la fréquence spatiale maximale furax et l'angle de diffraction maximum 0 comme sin(0) Àfmax. 30 Ainsi, comme illustré par la Figure 5, la région de contribution d'un point lumineux à la profondeur z est donnée par son rayon maximum Rmax tel que Rmax = z tan(0) = z tan (sin-1 (M) 2p (7) 3036197 14 La fonction fenêtre h peut ainsi être définie comme (8) h(x, y) = si ix2 + y2 i, comme illustré par la Figure 8. 30 L'étape de traitement E4 met en oeuvre la sous-étape E42 déjà décrite de calcul de l'onde lumineuse oi émise par les points de la portion de scène comprise dans le plan Pi qui ont une amplitude non nulle selon l'équation (2) déjà décrite. 3036197 16 En E47, on calcule l'amplitude des ondes lumineuses émises par des plans précédents à partir du plan basé champ P1 jusqu'au plan et bloquées par des obstacles occultants de la scène comprise dans P. Il s'agit encore une fois des points de la scène qui présentent une amplitude non nulle. 5 Il s'agit donc de gérer les occultations des ondes lumineuses provenant des plans précédents, par un plan basé points. Les inventeurs ont constaté que, bien que parfaitement adapté aux plans basés champ, dont l'onde lumineuse est calculée en utilisant l'opérateur Pw, la technique de prise en compte des occultations décrite en relation avec les Figures 6 et 7 pour les plans basés champs n'est pas du 10 tout adaptée aux plans basés points. En effet, cette technique nécessite de calculer préalablement la propagation des ondes lumineuses de proche en proche, c'est-à-dire d'un plan précédent sur le plan courant basé points. Or, à chaque étape de diffraction, l'onde complexe émise par la scène s'étend graduellement d'un plan à un autre, et l'onde incidente sur le plan Pi peut s'étendre sur un grand nombre d'échantillons. Ainsi, même si le plan Pi ne comporte qu'un petit nombre de points 15 de la scène pour lesquels l'amplitude de l'onde lumineuse oi diffère de zéro, l'onde lumineuse totale ui' émise par ce plan peut avoir un grand nombre de valeurs non nulles. Chaque valeur non nulle de l'onde Li; est alors considérée comme une source lumineuse sphérique par l'opérateur de propagation Ps précédemment décrit. Autrement dit, l'onde incidente va illuminer un certain nombres de points supplémentaires dans le plan basé points courant et complexifier le traitement 20 ultérieur de ce plan. Or, comme on l'a vu précédemment, la complexité algorithmique de Ps dépend du nombre de sources lumineuses sphériques situées dans le plan de source. L'utilisation de la technique précédente de prise en compte des occultations risquerait de rendre complètement inefficace le principe même de l'invention, à savoir de choisir la méthode de traitement adaptée à un plan en 25 fonction du nombre de points d'amplitude non nulle dans la portion de scène qu'il contient. L'invention propose une approche différente, adaptée aux plans basés points. Elle consiste à ne calculer l'amplitude de l'onde lumineuse provenant des plans précédents que pour les points d'amplitude non nulle de la scène du plan courant Pi'. Avantageusement, ce calcul s'appuie sur un masque binaire inversé. Il fournit l'amplitude de l'onde lumineuse bloquée par le plan courant. 30 En relation avec la Figure 8, on considère un obstacle occultant situé dans le plan Pd+1 et on utilise un masque binaire inversé défini par ad+1(x, y)- 1 - md+1(x, Y) (17) 3036197 17 En remplaçant ma' dans l'équation (14) par l'expression de ad+1 selon l'équation (17), ud+2 devient : ud+2(x,Y)=Pdz{od+i-F(1-ad+i) Pdz{ud}}(x,Y) ud+2(x,Y)= P2dz{ud}(x,Y) +Pdz{od+i-ad+i Pdz{ud}}(x,Y) 5 ud+2(x,Y)- P2dz{ud}(x,Y) +Pdz{0d+1}(x,Y) (18) Le calcul de ud+2 en utilisant le masque binaire inversé nécessite l'utilisation de l'opérateur de propagation P trois fois, alors qu'il n'est nécessaire de l'utiliser que deux fois en utilisant le 10 masque binaire m. Cependant, la propagation numérique du plan Pd au plan Pd+1 n'a besoin d'être calculée que dans la région définie par le masque binaire inversé ad+i, qui correspond aux coordonnées des points de la scène situés dans le plan Pd+1- C'est pour cette raison qu'on utilise l'opérateur de propagation basé points Ps, même pour le plan précédent basé champ PI. Ainsi, le nombre de valeurs non nulles de Cld+1 est le même que le nombre de points de la scène pour 15 lesquels l'amplitude de l'onde od+1 est différente de zéro. L'utilisation du masque binaire inversé présente l'avantage de ne pas augmenter le nombre de sources lumineuses sphériques dont l'onde doit être calculée à chaque pas de diffraction. Cette technique est ainsi extrêmement efficace pour la prise en compte des occultations par les plans basés points. Cette technique de gestion des occultations se généralise pour le plan courant basé points Pi par la 20 formule suivante : ub (x , y) = (x , y) Daz{ui} (x , y) (19) où ubi est l'onde lumineuse propagée depuis les plans précédents et bloquée par les points de la scène d'amplitude non nulle du plan basé points P. En E48, on calcule l'onde lumineuse totale ui émise par le plan Pi en soustrayant l'onde obi 25 bloquée à l'onde lumineuse émise par la portion de scène comprise dans le plan P. ui(x, y) = oi(x, y) - ub (x , y) = o i(x , y) - ai(x , y) z{ui} (x , y) (20) On comprend que le fait de retrancher de l'amplitude de l'onde lumineuse émise par le plan basé points courant Pi celle de l'onde lumineuse provenant des plans précédents et bloquée par le 30 plan courant permet de transmettre au plan basé champ suivant l'information de ce qui a été occulté dans le plan intercalaire P. 3036197 18 La valeur de l'onde lumineuse Li; obtenue pour le plan basé point Pi contribue au calcul de l'onde lumineuse ud émise par les plans basés points intercalaires Pi, avec j = l+1,...d-1 et propagée sur le plan Pd . Elle est ensuite exploitée pour calculer l'onde lumineuse totale émise ud par le plan basé champ suivant Pd, selon l'équation (16) précédemment décrite : 5 ud (x, y) = od (x, y) + md (x, y) [ud (x, y) (x, y)] (16) De cette manière, on compense au niveau du plan basé champ Pd, le fait d'avoir propagé les ondes lumineuses émises par les plans précédents sans prendre en compte les occultations par les plans intercalaires. 10 En relation avec la Figure 2, on revient maintenant sur l'étape E2 d'évaluation d'un seuil pour le plan courant Pi selon un mode de réalisation de l'invention. Afin de pouvoir implémenter la méthode proposée, il faut déterminer expérimentalement la valeur de seuil M -d,max- On note Tis le temps nécessaire au calcul de la propagation et de l'occultation de l'onde lumineuse issue de la scène par le plan Pi comportant Mi points lumineux en utilisant 15 rapproche basée points et Tile le temps nécessaire à ce même calcul en utilisant l'approche basée champ. Tis correspond à la somme du temps nécessaire au calcul de l'occultation de l'onde lumineuse par le plan Pi en utilisant la technique du masque binaire inversé (étape E47) et du temps nécessaire au calcul de la propagation de l'onde lumineuse émise par le plan Pi en utilisant l'opérateur Ps (étape E44): 20 Td =aMdM+(3MdN2 (21) où a et )3 sont des coefficients constants, M est le nombre total de points de la scène situés dans les plans d'indice j compris entre I et i-1 et N2 le nombre total de pixels de l'hologramme H : m = lm; 1=1 Tile est donné par Tiw = yN21og (N) (22) 25 où y est un coefficient constant. a, )3 et y doivent être déterminés de manière expérimentale et dépendent des capacités et performances du système sur lequel la méthode est exécutée, notamment en termes de puissance de calcul. A titre d'exemple, on donne ci-dessous des valeurs adaptées à notre système : a = 1, )3 = 2 et y = 5,35.10-8 30 Afin de maximiser l'efficacité de notre méthode, la valeur de seuil M -d,max est fixée telle que Ts- Tw - "21°g (N) (23) d - d Md,max aM+131V^2 3036197 19 Un avantage de ce mode de réalisation de l'invention est qu'il permet de déterminer la technique la plus adaptée au traitement du plan courant sur la base d'une comparaison réaliste des temps de calcul associés au deux méthodes mises en compétition. On notera que l'évaluation du seuil peut avantageusement être mise en oeuvre au cours 5 d'une première passe de prétraitement de la pluralité de plans. En effet, selon le mode de réalisation de l'invention qui vient d'être décrit, le nombre de points par plan basé points pris en compte dans le procédé reste inchangé et égal au nombre de points d'amplitude non nulle de la scène Sc d'origine. On notera également qu'en variante on pourrait tout aussi bien évaluer le seuil M -d,max à 10 partir d'estimations de complexité de calcul de la première et de la deuxième méthode pour le plan courant plutôt qu'à partir d'estimations de temps de calcul, de façon à rester indépendant d'éventuelles optimisations d'implémentation de la méthode sur un système particulier. En relation avec Figure 9, on décrit maintenant un mode de réalisation de l'invention permettant de simplifier l'étape E47 de calcul de l'amplitude obi bloquée par les points de la portion 15 de scène comprise dans un plan intercalaire basé points P. On considère un plan basé point Pi et un plan précédent Pi_i.Le plan occultant Pi est divisé en zones de A à H. Par exemple, ces zones présentent des formes géométriques simples, typiquement rectangulaires et ont toutes les mêmes dimensions. Dans l'exemple de la Figure 9, on considère les points q1, q2, q3 et q4 du plan Pi, disposés 20 de telle sorte que : - q2 E A ; - qi E B ; - q3, q4 e C. On remarque que les zones D, E, F, G et H sont vides. 25 Dans le plan précédent Pi-1, on considère les points qk, avec k entier compris entre 1 et 5. Pour chaque point pn, avec n entier compris entre 1 et 4, du plan Pi_i, on calcule sur quelles zones du planPi sa lumière émise s'étend. On exploite pour cela la région de contribution d'une onde lumineuse sphérique décrite précédemment en relation avec la Figure 5. Avantageusement, on détermine les zones du plan Pi qui ont une intersection non nulle avec le disque de rayon Rmax 3036197 20 projeté du plan précédent Pi-1 sur le plan courant Pi. Par exemple, on voit que le point Pi émet sur les zones B et D. Pour chaque zone du plan Pi, on définit une liste I comprenant les points pk du plan Pi_1 qui l'illuminent: 5 - 1(A) = {132}; - 1(B) = {pl} ; - 1(C) = {132,134} ; - 1(D) = P2, P3} ; - 1(E) = {134,P5} ; 10 - l(F) = {133,P5}; - l(G)= tP51; - 1(1)= {P5}. Chaque point q; du plan Pi n'occulte donc que la lumière provenant des points pk du plan 15 qui illuminent la zone dans laquelle il se trouve. Par exemple, le point qi E B n'occulte que la lumière provenant des points compris dans la liste 1(B), c'est-à-dire que la lumière provenant de pi. De cette manière, on ne calcule l'amplitude ob de l'onde lumineuse bloquée par la portion de scène comprise dans le plan Pi que pour les points de Pi_1 qui contribuent effectivement à illuminer les zones de Pi comprenant des points de la scène d'amplitude non nulle. La complexité 20 calculatoire de l'étape E47 se trouve donc réduite. On notera que l'invention qui vient d'être décrite, peut être mise en oeuvre au moyen de composants logiciels et/ou matériels. Dans cette optique, les termes « module » et « entité », utilisés dans ce document, peuvent correspondre soit à un composant logiciel, soit à un composant matériel, soit encore à un ensemble de composants matériels et/ou logiciels, aptes à mettre en 25 oeuvre la ou les fonctions décrites pour le module ou l'entité concerné(e). 3036197 21 En relation avec la figure 10, on présente maintenant un exemple de structure simplifiée d'un dispositif 100 de génération numérique d'un hologramme à partir d'une scène tridimensionnelle selon l'invention. Le dispositif 100 met en oeuvre le procédé de génération selon l'invention qui vient d'être décrit en relation avec les Figures 2 et 4. 5 Par exemple, le dispositif 100 comprend une unité de traitement 110, équipée d'un processeur pl, et pilotée par un programme d'ordinateur Pgl 120, stocké dans une mémoire 130 et mettant en oeuvre le procédé de selon l'invention. A l'initialisation, les instructions de code du programme d'ordinateur Pg1 120 sont par exemple chargées dans une mémoire RAM avant d'être exécutées par le processeur de l'unité de 10 traitement 110. Le processeur de l'unité de traitement 110 met en oeuvre les étapes du procédé décrit précédemment, dans ces différents modes de réalisation, selon les instructions du programme d'ordinateur 120. Dans cet exemple de réalisation de l'invention, le dispositif 100 comprend au moins une unité CARD de comptage d'un nombre de points de la scène d'amplitude non nulle dans le plan 15 courant, une unité DEC de choix d'une première ou d'une deuxième technique de propagation d'une onde lumineuse émise par le plan courant en fonction d'un nombre de points d'amplitude non nulle de la portion de scène comprise dans le plan courant et d'une valeur de seuil prédéterminée, la première technique, dite basée points, destinées à calculer une somme d'ondes lumineuses émises par des sources ponctuelles constituées par les points de la portion de scène d'amplitude 20 non nulle du plan courant sur un plan suivant, la deuxième technique, dite basée champ, apte à calculer globalement une onde lumineuse émise par la portion de scène comprise dans le plan courant sur un plan suivant , une unité de traitement du plan courant selon la technique de propagation choisie, apte à calculer une onde lumineuse émise par le plan courant et une unité GET H d'obtention de l'hologramme à partir des ondes lumineuses calculées pour la pluralité de plans. 25 Avantageusement, l'unité PROC de traitement du plan courant comprend des sous-unités d'identification du dernier plan basé champ traité, de calcul de l'onde lumineuse émise par la portion de scène comprise dans le plan courant et éventuellement, selon la technique choisie, une sous-unité de calcul de l'onde lumineuse propagée des plans précédents, une sous unité de calcul de l'onde lumineuse occultée par le plan courant et une sous unité de calcul de l'onde lumineuse 30 totale Li; émise par le plan courant. 3036197 22 Le dispositif 100 comprend en outre une unité M1 de stockage des ondes lumineuses calculées pour la pluralité de plans. Ces unités sont pilotées par le processeur p..1 de l'unité de traitement 110. De façon avantageuse, un tel dispositif 100 peut être intégré à un terminal d'utilisateur TU. A 5 titre d'exemple, un tel terminal peut être un ordinateur personnel (pour « personal computer », en anglais), une tablette, un téléphone intelligent (pour « smartphone », en anglais) ou tout système matériel et/ou logiciel comprenant des ressources de calcul suffisantes pour mettre en oeuvre le procédé selon l'invention. Un tel terminal comprend en outre un module de commande par l'intermédiaire de moyens de connexion filaires ou non, d'un dispositif de restitution DR, de type 10 téléviseur 2D, 3D ou 3D holographique, apte à restituer l'hologramme généré par le procédé selon l'invention. Bien sûr, pour un dispositif de restitution non holographique, il est nécessaire de prévoir, par exemple au sein du module de commande MC, un traitement intermédiaire d'adaptation du champ lumineux généreux afin de le rendre visualisable sur le dispositif de restitution DR. 15 L'invention qui vient d'être présentée permet de générer numériquement un champ lumineux d'hologramme à partir d'une scène tridimensionnelle. En relation avec la Figure 11, on présente maintenant un autre exemple d'usage possible de l'invention. On considère qu'on dispose d'un premier hologramme H qui a été généré à l'aide du procédé de l'invention ou d'un autre procédé. 20 On ne connait pas a priori la première scène à l'origine de l'hologramme H. Dans le scénario dit d'holographie augmentée que l'on considère, on souhaite intégrer une deuxième scène comprenant au moins un objet supplémentaire devant la première scène. Dans l'exemple de la figure 11, la deuxième scène se limite à un objet, en l'espèce un lapin. Pour pouvoir s'intégrer correctement dans la scène, le lapin doit correctement occulter les éléments 25 de la première scène se trouvant derrière lui. Comme nous ne connaissons pas la première scène correspondant à l'hologramme H, il est impossible d'utiliser des tests de visibilité entre le lapin et les autres éléments de la scène. On propose de mettre en oeuvre le procédé de génération d'un champ lumineux d'hologramme selon l'invention, en découpant la deuxième scène en une pluralité de plans équidistants à partir du 30 plan de l'hologramme H, tel que Po =H. La pluralité de plans comprend Nz plans qui sont traités selon l'invention comme des plans basés champs ou basés points, le plan Po = H et le plan P - Nz-1 étant des plans basés champ. 3036197 23 Le nouvel hologramme H' est obtenu par propagation de l'onde lumineuse uNz_i calculée pour le dernier plan basé champ PNz_i sur le plan de l'hologramme H'. Ensuite, on peut éventuellement propager le champ lumineux H' obtenu dans le plan H pour conserver la position de l'hologramme original et voir apparaître la deuxième scène devant la 5 première. Il va de soi que les modes de réalisation qui ont été décrits ci-dessus ont été donnés à titre purement indicatif et nullement limitatif, et que de nombreuses modifications peuvent être facilement apportées par l'homme de l'art sans pour autant sortir du cadre de l'invention. 10 REVENDICATIONS1. Procédé de génération numérique d'un hologramme (H) dans un plan (PH), dit plan de l'hologramme, à partir d'une scène tridimensionnelle (Sc), comprenant une étape (EO) de Nz-1, découpage de ladite scène en une pluralité de plans (Po, P parallèles au plan de l'hologramme, caractérisé en ce que le procédé comprend les étapes suivantes, mises en oeuvre pour un plan (pi) de la pluralité de plans, dit plan courant, les plans étant parcourus du plus éloigné au plus proche du plan de l'hologramme : - comptage (E1) d'un nombre de points (Mi) de la scène d'amplitude non nulle dans le plan courant ; - choix (E2) d'une première ou d'une deuxième technique (T1,T2) de propagation d'une onde lumineuse émise par le plan courant en fonction du nombre de points d'amplitude non nulle de la portion de scène comprise dans le plan courant et d'une valeur de seuil prédéterminée (Mi,max), la première technique (T1), dite basée points, étant destinée à calculer la propagation d'une somme d'ondes lumineuses émises par des sources ponctuelles constituées par les points de la portion de scène d'amplitude non nulle du plan courant sur un plan suivant, la deuxième technique (T2), dite basée champ, destinée à calculer globalement la propagation d'un champ lumineux émis par la portion de scène comprise dans le plan courant sur un plan suivant ; - traitement (E4) du plan courant en fonction de la technique de propagation choisie, destiné à calculer une onde lumineuse (ui) émise par le plan courant; lesdites étapes étant répétées pour la pluralité de plans, et en ce que le procédé comprend en outre une étape (E6) d'obtention de l'hologramme à partir des ondes lumineuses calculées pour la pluralité de plans. 2. Procédé de génération numérique d'un hologramme selon la revendication 1, caractérisé en ce que ladite étape de traitement (E4) du plan courant comprend les sous-étapes suivantes : - identification (E41) du dernier plan traité selon la deuxième technique, dit dernier plan basé champ (PI); - calcul (E42) de l'onde lumineuse émise par la portion de scène comprise dans le plan courant ; 3036197 25 - lorsque la deuxième technique a été choisie pour le plan courant, dit plan basé champ, l'étape de traitement comprend en outre les sous-étapes suivantes : o calcul (E44) d'une propagation de l'onde lumineuse émise par au moins le dernier plan basé champ (PI) sur le plan courant, selon la technique choisie pour ledit plan; 5 o calcul (E46) de l'onde lumineuse émise par le plan courant par sommation de l'onde lumineuse émise par le plan courant et au moins la propagation de l'onde lumineuse émise par au moins le dernier plan basé champ ; et en ce que l'étape (E6) d'obtention de l'hologramme comprend le calcul de la propagation de l'onde lumineuse émise par le plan basé champ le plus proche de l'hologramme (PN,_1) sur le 10 plan de l'hologramme (PH). 3. Procédé de génération numérique d'un hologramme selon la revendication 2, caractérisé en ce que, lorsque ledit au moins un plan précédent comprend en outre au moins un plan basé points (Pi), intercalé entre le dernier plan basé champ et le plan courant, le calcul de la propagation de l'onde lumineuse émise par au moins le dernier plan basé champ sur le plan 15 courant (Pi) comprend en outre le calcul d'une propagation de l'onde lumineuse émise par les points d'amplitude non nulle de la portion de scène comprise dans ledit au moins un plan intercalaire sur le plan courant selon la première technique. 4. Procédé de génération numérique d'un hologramme selon l'une des revendications 2 et 3, caractérisé en ce que l'étape de traitement (E4) du plan courant comprend une sous-étape 20 (E45) de gestion d'une occultation de l'onde lumineuse propagée par au moins le dernier plan basé champ (PI) dans le plan courant par application d'un masque binaire occultant (nni) à ladite onde lumineuse, ledit masque prenant une valeur nulle aux points de la portion de scène comprise dans le plan courant (Pi) et une valeur non nulle en dehors de cette portion et en ce que la sous-étape (E46)de calcul de l'onde lumineuse (ui) émise par le plan courant 25 comprend une sommation de l'onde lumineuse (oi) émise par le plan courant à l'onde lumineuse obtenue après application du masque binaire occultant (mi). 5. Procédé de génération numérique d'un hologramme selon l'une des revendications 3 et 4, caractérisé en ce lorsque la première technique a été choisie pour traiter le plan courant, dit 30 plan dit basé points, l'étape de traitement (E4) comprend en outre les sous-étapes suivantes : 3036197 26 - calcul (E47) d'une onde lumineuse (obi) propagée par au moins le dernier plan basé champ (PI) selon la première technique et bloquée par les points d'amplitude non nulle de la portion de scène comprise dans le plan courant ; - calcul (E48) de l'onde lumineuse (ui) émise par le plan courant par soustraction de 5 l'onde lumineuse (obi) bloquée par les points d'amplitude non nulle de la portion de scène comprise dans le plan courant à l'onde lumineuse (oi) émise par lesdits points. 6. Procédé de génération numérique d'un hologramme selon la revendication 5, caractérisé en ce que la sous-étape (E47) de calcul de l'onde lumineuse bloquée comprend le calcul de l'onde lumineuse propagée au moins par le dernier plan basé champ (PI) selon la première technique 10 et l'application d'un masque binaire occultant inversé (ai) à l'onde lumineuse calculée, ledit masque prenant une valeur non nulle aux points de la portion de scène comprise dans le plan courant et une valeur nulle en dehors de cette portion. 7. Procédé de génération numérique d'un hologramme selon la revendication 6, caractérisé en ce que la sous-étape de calcul (E47) de l'onde lumineuse bloquée comprend un découpage en 15 une pluralités de zones du plan courant, une détermination des points d'un plan précédent à partir du dernier plan basé champ, contribuant à la propagation d'une onde lumineuse dans au moins une zone de la pluralité de zones comprenant des points d'amplitude non nulle de la portion de scène comprise dans le plan courant, le calcul de l'onde lumineuse bloquée n'étant effectué que pour les points déterminés. 20 8. Procédé de génération numérique d'un hologramme selon l'une des revendications précédentes, caractérisé en ce qu'il comprend une étape (E3) d'évaluation du seuil (Mi,max) au moins en fonction d'un de temps de traitement estimé pour la première technique de propagation, d'un temps de traitement estimé par la deuxième technique de propagation pour le plan courant et un nombre de points d'amplitude non nulle de la portion de scène 25 comprise dans les plans précédents. 9. Dispositif (100) de génération numérique d'un hologramme (H) dans un plan (PH), dit plan de l'hologramme, à partir d'une scène tridimensionnelle (Sc), comprenant une unité (CUT) de découpage de ladite scène en une pluralité de plans (Po, PNz-1) parallèles au plan de l'hologramme, caractérisé en ce que le dispositif comprend les unités suivantes, aptes à être 30 mises en oeuvre pour un plan (Pi) de la pluralité de plans, dit plan courant, les plans étant parcourus du plus éloigné au plus proche du plan de l'hologramme : 3036197 27 - comptage (CARD) d'un nombre de points de la scène d'amplitude non nulle dans le plan courant ; - choix (DEC) d'une première ou d'une deuxième technique de propagation d'une onde lumineuse émise par le plan courant en fonction d'un nombre de points d'amplitude non 5 nulle de la portion de scène comprise dans le plan courant et d'une valeur de seuil prédéterminée, la première technique (T1), dite basée points, destinée à calculer une somme d'ondes lumineuses émises par des sources ponctuelles constituées par les points de la portion de scène d'amplitude non nulle du plan courant sur un plan suivant, la deuxième technique (T2), dite basée champ, destinée à calculer globalement une 10 onde lumineuse émise par la portion de scène comprise dans le plan courant sur un plan suivant ; - traitement (PROC) du plan courant selon la technique de propagation choisie, destiné à calculer une onde lumineuse émise par le plan courant; la mise en oeuvre desdites unités étant répétée pour la pluralité de plans, 15 et en ce que le dispositif comprend en outre une unité (GET H) d'obtention de l'hologramme à partir des ondes lumineuses calculées pour la pluralité de plans. 10. Terminal d'utilisateur (TU) comprenant un module de commande (MC) d'un dispositif de restitution (DR) d'un hologramme (H), caractérisé en ce qu'il comprend en outre un dispositif (100) de génération numérique dudit hologramme selon la revendication 9. 20 11. Programme d'ordinateur (Pg1) comprenant des instructions pour la mise en oeuvre du procédé de traitement selon l'une des revendications 1 à 8, lorsqu'il est exécuté par un processeur. 12. Support d'enregistrement lisible par un processeur, caractérisé en ce qu'il est apte à mémoriser le programme d'ordinateur selon la revendication 11.