FR 3084192 A1 20200124 FR 1856732 A 20180719 La présente invention concerne un procédé de rendu d'une scène tridimensionnelle et des programme d'ordinateur et dispositif correspondants. L'invention s'applique plus particulièrement à un procédé de rendu d'une scène tridimensionnelle, comportant, pour chacun d'au moins un pixel d'une image représentant une vue de la scène tridimensionnelle : la détermination, à partir de la scène tridimensionnelle et de paramètres de la vue, d'une fonction de contribution donnant une contribution d'un chemin de lumière partant d'une source de lumière de la scène tridimensionnelle, à une valeur du pixel, une première sélection aléatoire de premiers chemins de lumière, et l'évaluation d'un estimateur d'une intégrale de la fonction de contribution, l'estimateur comportant un estimateur de Monte Carlo utilisant les premiers chemins de lumière. L'estimateur de Monte Carlo présente une variance qui diminue avec le nombre de chemins de lumière aléatoires utilisés. Cette variance est responsable d'instabilités spatiale et temporelle lors du rendu. En effet, pour l'instabilité spatiale, deux pixels dont les valeurs devraient être identiques ou proches (par exemple deux pixels adjacents dans l'image) sont estimés à partir de chemins de lumière aléatoires pouvant être différents, de sorte que les estimations peuvent être différentes, d'autant plus que la variance est élevée. Pour l'instabilité temporelle, dans le cas d'une vidéo, un même pixel devrait souvent avoir la même valeur ou des valeurs proches d'une image à la suivante. Or, comme sa valeur est à chaque fois estimée à partir de chemins de lumière aléatoires pouvant être différents d'une image à la suivante, sa valeur peut varier d'une image à la suivante, d'autant plus que la variance est élevée. Ainsi, pour limiter les instabilités spatiale et temporelle, il est connu d'augmenter le nombre de chemins de lumière aléatoires utilisés par l'estimateur de Monte Carlo. Cette solution présente néanmoins comme inconvénient d'augmenter la complexité des calculs des valeurs des pixels, et donc d'augmenter le temps de calcul et/ou la puissance du matériel nécessaire pour réaliser ces calculs. Il peut ainsi être souhaité de prévoir un procédé de rendu d'une scène tridimensionnelle qui permette de s'affranchir d'au moins une partie des problèmes et contraintes précités. Il est donc proposé un procédé de rendu d'une scène tridimensionnelle, comportant, pour chacun d'au moins un pixel d'une image représentant une vue de la scène tridimensionnelle : 2 la détermination, à partir de la scène tridimensionnelle et de paramètres de la vue, d'une fonction de contribution donnant une contribution d'un chemin de lumière partant d'une source de lumière de la scène tridimensionnelle, à une valeur du pixel, 5 une première sélection aléatoire de premiers chemins de lumière, et l'évaluation d'un estimateur d'une intégrale de la fonction de contribution, l'estimateur comportant un estimateur de Monte Carlo utilisant les premiers chemins de lumière, comportant en outre : 10 - la détermination d'au moins un paramètre d'une fonction de contribution approchée, de sorte que la fonction de contribution approchée donne une contribution d'un chemin de lumière à la valeur du pixel proche de la contribution donnée par la fonction de contribution, la fonction de contribution approchée comportant un réseau de neurones comportant : 15 - une couche d'entrée d'au moins un neurone, - une unique couche cachée d'au moins un neurone, - une couche de sortie d'au moins un neurone - des connexions entre le ou les neurones de la couche d'entrée et le ou les neurones de la couche cachée auxquelles des poids sont 20 respectivement associés, ces poids étant indépendants du ou des paramètres, et - des connexions entre le ou les neurones de la couche cachée et le ou les neurones de la couche de sortie auxquelles des poids dépendant du ou des paramètres sont respectivement associés, 25 et dans lequel l'estimateur comporte en outre une moyenne de contributions données par la fonction de contribution approchée appliquée à au moins une partie des premiers chemins de lumière. En effet, la présence dans l'estimateur de l'invention, noté Ê, des contributions approchées permet d'obtenir une variance réduite par rapport à la variance de 30 l'estimateur de Monte Carlo seul. L'utilisation de l'estimateur Ê permet donc d'obtenir une plus grande stabilité temporelle et spatiale lors du rendu que l'utilisation de l'estimateur de Monte Carlo seul. En outre, le choix de la fonction de contribution approchée sous la forme d'un réseau de neurones comportant une unique couche cachée et dans lequel seuls les poids des connexions de sortie (entre la couche 35 cachée et la couche de sortie) sont paramétrés, facilite la détermination du ou des 3 paramètres pour que la fonction de contribution approchée soit proche de la fonction de contribution. De façon optionnelle, le procédé comporte en outre une deuxième sélection aléatoire, indépendante de la première sélection, de deuxièmes chemins de lumière, 5 et le ou les paramètres de la fonction de contribution approchée sont déterminés à partir des deuxièmes chemins de lumière. De façon optionnelle également, la détermination du ou des paramètres est réalisée de manière à minimiser une distance entre la fonction de contribution approchée et la fonction de contribution. 10 De façon optionnelle également, la distance utilise les deuxièmes chemins de lumière. De façon optionnelle également, la distance comporte la somme des carrés des différences entre la fonction de contribution approchée et la fonction de contribution pour les deuxièmes chemins de lumière. 15 De façon optionnelle également, la distance comporte en outre un paramètre de régularisation et la détermination du ou des paramètres comporte la détermination du paramètre de régularisation, puis la détermination du ou des paramètres en utilisant le paramètre de régularisation déterminé. De façon optionnelle également, l'estimateur peut s'écrire matriciellement 20 sous la forme A[Im- M]B, où lm est une matrice identité de taille m, A et B sont des vecteurs et M est une matrice dépendant du paramètre de régularisation, et le paramètre de régularisation est déterminé de sorte que la valeur propre maximale de M présente au plus un écart prédéfini avec 1. De façon optionnelle également, les poids associés aux connexions entre le 25 ou les neurones de la couche d'entrée et le ou les neurones de la couche cachée sont déterminés de manière aléatoires. De façon optionnelle également, l'estimateur comporte en outre une espérance de la fonction de contribution approchée. Il est également proposé un programme d'ordinateur téléchargeable depuis un 30 réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions pour l'exécution des étapes d'un procédé de rendu d'une scène tridimensionnelle selon l'invention, lorsque ledit programme est exécuté sur un ordinateur. 4 Il est également proposé un dispositif de rendu d'une scène tridimensionnelle, comportant, pour chacun d'au moins un pixel d'une image représentant une vue de la scène tridimensionnelle : - des moyens de détermination, à partir de la scène tridimensionnelle et de 5 paramètres de la vue, d'une fonction de contribution donnant une contribution d'un chemin de lumière partant d'une source de lumière de la scène tridimensionnelle, à une valeur du pixel, des moyens de réalisation d'une première sélection aléatoire de premiers chemins de lumière, et 10 des moyens d'évaluation d'un estimateur d'une intégrale de la fonction de contribution, l'estimateur comportant un estimateur de Monte Carlo utilisant les premiers chemins de lumière, comportant en outre : des moyens de détermination d'au moins un paramètre d'une fonction de 15 contribution approchée, de sorte que la fonction de contribution approchée donne une contribution d'un chemin de lumière à la valeur du pixel (p) proche de la contribution donnée par la fonction de contribution, la fonction de contribution approchée comportant un réseau de neurones comportant : 20 - une couche d'entrée d'au moins un neurone, - une unique couche cachée d'au moins un neurone, - une couche de sortie d'au moins un neurone - des connexions entre le ou les neurones de la couche d'entrée et le ou les neurones de la couche cachée auxquelles des poids sont 25 respectivement associés, ces poids étant indépendants du ou des paramètres, et - des connexions entre le ou les neurones de la couche cachée et le ou les neurones de la couche de sortie auxquelles des poids dépendant du ou des paramètres sont respectivement associés, 30 et dans lequel l'estimateur comporte en outre une moyenne de contributions données par la fonction de contribution approchée appliquée à au moins une partie des premiers chemins de lumière. L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés 35 dans lesquels : 5 la figure 1 représente schématiquement la structure générale d'un dispositif de rendu de scène tridimensionnelle, selon un mode de réalisation de l'invention, la figure 2 illustre un exemple de scène tridimensionnelle, et 5 - la figure 3 illustre les étapes successives d'un procédé de rendu de scène tridimensionnelle, selon un mode de réalisation de l'invention. En référence à la figure 1, un dispositif 100 de rendu de scène tridimensionnelle va à présent être décrit. Dans l'exemple décrit, le dispositif 100 est un dispositif informatique 10 comportant une unité centrale de traitement 102, telle qu'un microprocesseur, et une mémoire 104 dans laquelle sont enregistrés un programme d'ordinateur 106, des données 108 décrivant une scène tridimensionnelle dans un espace tridimensionnel virtuel, et des données 110 représentant des paramètres d'une vue souhaitée de cette scène tridimensionnelle, cette vue étant destinée à être réalisée par une caméra 15 (virtuelle). Les paramètres de la vue incluent généralement des paramètres extrinsèques à la caméra (par exemple, une position et une orientation de la caméra dans l'espace tridimensionnel), ainsi que des paramètres intrinsèques à la caméra (par exemple, une focale, un pixel central, une distorsion optique, etc. de la caméra). Pour chacune des étapes du procédé de rendu de scène tridimensionnelle qui 20 sera décrit plus loin, le programme d'ordinateur 106 comporte des instructions pour la réalisation de cette étape par l'unité centrale de traitement 102. De manière générale, ces instructions peuvent être regroupées ou séparées selon toutes les combinaisons possibles en un ou plusieurs programmes d'ordinateur. Leur distinction est en outre purement fonctionnelle. Enfin, leurs fonctions pourraient 25 aussi être au moins en partie micro programmées ou micro câblées dans des circuits intégrés dédiés. Ainsi, en variante, le dispositif informatique 100 incluant l'unité de traitement 102 pourrait être remplacé par un dispositif électronique composé uniquement de circuits numériques (sans programme d'ordinateur) pour la réalisation des mêmes actions. 30 En référence à la figure 2, un exemple de scène tridimensionnelle 200 est illustré. Dans cet exemple, la scène 200 comporte un seul objet, un cube 202, et une seule source de lumière 204, ayant tous les deux des positions prédéfinies dans un espace tridimensionnel virtuel. 6 Une image 206 destinée à représenter une vue de la scène 200 est également illustrée sur la figure 2. L'image 206 comporte une pluralité de pixels, parmi lesquels un pixel p est représenté. De manière générale, chaque source de lumière de la scène émet de la 5 lumière qui peut atteindre l'image soit directement, soit indirectement après avoir interagi avec les objets de la scène, par exemple après réflexion sur un ou plusieurs objets ou bien après déviation par un ou plusieurs objets. Par exemple, sur la figure 2, des chemins de lumière 208, 210 partent de la source de lumière 204 et atteignent le pixel p. 10 En référence à la figure 3, un procédé 300 de rendu de scène tridimensionnelle va à présent être décrit. Dans l'exemple décrit, le procédé 300 réalise un rendu de la scène 200 de la figure 2. Cependant, le procédé 300 pourrait être utilisé pour réaliser le rendu d'autres scènes tridimensionnelles. En outre, le procédé 300 tel que décrit est mis en oeuvre par le dispositif 100 de la figure 1, mais il 15 pourrait être mis en oeuvre par d'autres moyens. Au cours d'une étape 302, le programme 106 récupère la scène 200 à rendre et les paramètres de la vue souhaitée, par exemple en lisant les données 108, 110 enregistrées dans la mémoire 104. Au cours d'une étape 304, le programme 106 réalise un rendu de la scène 20 200 afin de fournir l'image 206 représentant une vue de cette scène 200. Pour chaque pixel p de l'image 206, l'étape 304 comporte les étapes 306 et 308 suivantes. Au cours d'une étape 306, une fonction co(x) est déterminée à partir de la scène 200, des paramètres de la vue souhaitée et du pixel p considéré. La fonction 25 yo(x) est une fonction destinée à recevoir en entrée un chemin de lumière x partant d'une source de lumière de la scène 200 (par exemple, la source de lumière 204) et arrivant au pixel p, et à fournir une contribution de ce chemin de lumière x à la valeur du pixel p. Les chemins de lumière 208, 210 de la figure 2 sont des exemples de tels chemins de lumière. La fonction cp sera ainsi appelée par la suite « fonction de 30 contribution ». La valeur théorique du pixel p est égale à l'intégrale de la fonction de contribution co(x) sur tous les chemins de lumière x possibles arrivant au pixel p. Ainsi, l'intégrale Ip suivante donne la valeur théorique du pixel p : cp(x)p(dx) 7 où B est l'ensemble des chemins de lumière possibles arrivant au pixel p et pt(dx) une mesure associée à l'ensemble B. Théoriquement, l'ensemble B comporte un nombre infini de chemins de lumière. En outre, la fonction de contribution cp(x) n'est généralement pas intégrable 5 analytiquement. Ainsi, l'intégrale /p est très difficile à calculer de manière exacte et il faut passer par un estimateur de cette intégrale /p. Un exemple de fonction yo(x) pouvant être utilisée dans la présente invention est donné dans l'article « Sequential Monte Carlo Instant Radiosity », de Hedman et al., 2016, Proceeding 13D '16 Proceedings of the 20th ACM SIGGRAPH Symposium 10 on Interactive 3D Graphics and Games, pages 121-128. Au cours d'une étape 308, une valeur du pixel p est calculée à partir d'un estimateur Ê de l'intégrale /p. Cet estimateur Ê est égal à la somme, d'une part, d'un estimateur de Monte Carlo E de l'intégrale /p et, d'autre part, d'un terme stabilisateur T destiné, comme 15 cela sera expliqué par la suite à améliorer la stabilité temporelle et spatiale de l'estimation : Ê= E + T Comme cela est connu en soi, un estimateur de Monte Carlo est conçu pour estimer l'intégrale Ip en utilisant un nombre n prédéfini de chemins de lumière sélectionnés aléatoirement, par échantillonnage statistique, dans l'ensemble 20 B. Par exemple, un estimateur de Monte Carlo E réalise une moyenne des contributions de ces chemins de lumière {rd à la valeur du pixel p. L'échantillonnage statistique utilise par exemple la méthode, connue en soi, de tirage de rayons (« raytracing » en anglais). Dans l'exemple décrit, l'estimateur de Monte Carlo E réalise une moyenne 25 arithmétique des contributions des chemins de lumière {xi} et est ainsi donné par l'équation : n E({x(}) = t=1 Le terme stabilisateur T utilise au moins une partie des chemins de lumière (xi) utilisés par l'estimateur de Monte Carlo E, de sorte que le terme stabilisateur T présente une certaine corrélation c(T, E) avec l'estimateur de Monte Carlo E. 30 Dans l'exemple décrit, le terme stabilisateur T est donné de manière générale par l'équation suivante : 8 T = a(-T1 + m(Ti)) avec T1 une statistique corrélée à l'estimateur de Monte Carlo E, a proche de c(ETi) veril ' où c(E,T1) est la corrélation entre E et T1 et v(T1) est la variance de T1. Par exemple, (E,T, a est égal à c-v(Ti)) à 10% près, de préférence à 1% près, de préférence encore c()) exactement égal à ver1 5 Ce choix de a présente comme avantage de minimiser la variance de Ê, notée v(Ê), pour un Ti donné. En effet, en gardant a quelconque, v(Ê) est donnée par : v(Ê) = v(E) + a2v(T1) - 2a c(E,T1) qui est une équation du second degré de a admettant un minimum pour a = c(E'n) v(;) En outre, avec a est égal à ), la variance de Ê est donnée par : te, ) v(Ê) = v(E)(1 - c(E,T1))2 de sorte que plus la corrélation c(E,T1) est forte (c'est-à-dire proche de 1) plus la 10 variance v(Ê) est faible (proche de 0). En outre, dans l'exemple décrit, T1 est de la forme suivante : 1 T, = nIcp',(xi) i=1 spin est une fonction analytiquement intégrable destinée à être proche de la fonction de contribution cp de manière à donner une contribution approchée d'un chemin de lumière à la valeur du pixel p. La fonction Cpm sera ainsi appelée par la 15 suite « fonction de contribution approchée ». Par ailleurs, 140m) est l'espérance de la fonction de contribution approchée cpm. Ainsi, dans l'exemple décrit, le terme T1 est égal à une moyenne (moyenne arithmétique dans l'exemple décrit) des contributions approchées des chemins de lumière {xi}. 20 L'espérance du terme T1 est ainsi égale à l'espérance de la fonction de contribution approchée (Pm : 1 Pal/ = Li (Pm(Xi)) = 11(0m) i=1 Le terme stabilisateur T présente donc une espérance nulle : PM= u (a(-Ti +l(Ti))) = a (-P(Tl) ± Klei))) = a(-11(T1) ± Rai) = 9 Ainsi, l'espérance de l'estimateur Ê est égale à l'espérance de l'estimateur de Monte Carlo E, de sorte que l'estimation donnée par l'estimateur Ê est, en moyenne, une bonne estimation de l'intégrale /p. La fonction de contribution approchée çorp est une fonction paramétrée 5 comportant au moins un paramètre, ce ou ces paramètres étant destinés à être déterminés de manière à minimiser une distance entre la fonction de contribution approchée ify', et la fonction de contribution (p. En outre, la fonction de contribution approchée Cpin comporte un réseau de neurones, comportant des neurones et des connexions entre ces neurones, ces 10 connexions présentant des poids respectifs, dont au moins une partie dépend du ou des paramètres de la fonction de contribution approchée Com. Dans l'exemple décrit, la fonction de contribution approchée Om est formée d'un réseau de neurones dit « de projections de caractéristique aléatoire » (de l'anglais « random feature maps »). Un tel réseau de neurones comporte une couche 15 d'entrée de neurone(s), une unique couche cachée de neurone(s) et une couche de sortie de neurone(s). Les neurones de la couche d'entrée et ceux de la couche cachée sont connectés par des connexions auxquelles des poids W k aléatoires sont respectivement associés. Les neurones de la couche cachée et de la couche de sortie sont connectés par des connexions auxquelles des poids sont respectivement 20 associés, ces derniers dépendant du ou des paramètres de la fonction de contribution approchée cpm. Ainsi, dans l'exemple décrit, la fonction de contribution approchée Cpm est donnée par l'équation : 1 171\1 c (W T x Z A avec : Z =1a( 1/17Tg) 25 où X = icm] sont m chemins de lumière sélectionnés aléatoirement, rn étant aussi le nombre de paramètres, N est le nombre de neurones de la couche cachée, W = [(01 .-- tow] sont les poids des connexions entre la couche d'entrée et la couche cachée, A = [ai ... am] sont les paramètres, Z A sont les poids des connexions entre la couche cachée et la couche de sortie, et a étant la fonction 30 d'activation des neurones de la couche cachée. 10 La fonction d'activation a peut être n'importe quelle fonction dérivable. Dans l'exemple décrit, la fonction d'activation a est donnée par l'équation : = [cosi [sin] Le nombre N de neurones de la couche cachée est prédéfini et choisi en fonction de la performance visée. 5 Avec le choix précédent de çom , il a été trouvé que cv(E&T;) était égal à 1. Ainsi, a est aussi égal à 1, de sorte que, dans ce cas, l'estimateur Ê est donné par l'équation : Ê = E - T1+ p.(Ti) L'étape 308 comporte les étapes 310 à 320 suivantes. Au cours d'une étape 310, les poids cok sont calculés de manière aléatoire. 10 Dans l'exemple décrit, les poids wk suivent une loi normale de moyenne nulle et de matrice de covariance I = g! où g est un nombre prédéfini et I une matrice identité. Au cours d'une étape 312, les m chemins de lumière . fi'. }1...m sont aléatoirement sélectionnés, par échantillonnage statistique, dans l'ensemble B. Ainsi, les chemins de lumière Vcil sont des variables indépendantes et identiquement 15 distribuées. Les chemins de lumière {xi} sont par exemple sélectionnés par la méthode, connue en soi, de tirage de rayons (« ray-tracing » en anglais). Au cours d'une étape 314, les paramètres À1 sont déterminés de manière à minimiser une distance entre la fonction de contribution approchée Ço', et la fonction de contribution O. Or, il est possible d'exploiter le fait que la matrice M est de la forme M = k1ZTZ + k2/'.,, avec k1 et k2 deux constantes, pour inverser la matrice M avec une complexité en o(N3) où N est 10 le nombre de neurones de la couche cachée. Si la matrice M était de la forme M = kiK + k2t, où K est une matrice ne pouvant pas s'écrire sous la forme ZTZ, l'inversion de la matrice M présenterait une complexité en o(m3) où ni est le nombre de tirages. Ainsi, si N est inférieur à m, il est possible grâce à l'invention d'exploiter la présence de ZTZ pour diminuer la complexité de l'inversion de la matrice M, et donc 15 de faciliter la détermination des paramètres Ai. Au cours d'une étape 316, l'espérance [1(0,,,) est calculée. Dans l'exemple décrit, où les poids wk suivent une loi normale de moyenne nulle et de matrice de covariance le et où la fonction d'activation a vaut rsj, l'espérance p(CA.,) est donnée par l'équation suivante : 1 [1(0m) = Triv3Ze {Ou (WT exp(-iWTi0)1A 20 où 0/, est la fonction caractéristique associée à l'espérance p. Au cours d'une étape 318, n chemins de lumière sont aléatoirement sélectionnés par échantillonnage statistique, dans l'ensemble D, indépendamment de la sélection réalisée à l'étape 312. Les chemins de lumière sont par exemple sélectionnés par la méthode, connue en soi, de tirage de rayons (c ray-tracing » en 25 anglais). Au cours d'une étape 320, l'estimateur Ê est évalué à partir des chemins de lumière [xi} sélectionnés à l'étape 318 : 1 v, eux()) = 2j(v(xi) - (Pm (xi) + /40,.)) 1=1 et le résultat est assigné comme valeur du pixel p. n 13 Le procédé 300 retourne alors à l'étape 306 pour un nouveau pixel p de l'image 206, jusqu'à ce que les valeurs de tous les pixels de l'image 206 aient été calculées. Il apparaît clairement qu'un procédé de rendu de scène tridimensionnelle tel 5 que celui décrit précédemment permet d'améliorer la stabilité spatiale et temporelle des valeurs des pixels de l'image. En outre, il apparaît que les deux points suivants de l'invention sont indépendants l'un de l'autre : 1. L'estimateur Ê est de la forme Ê = E + a(-T, + µ(T1), avec a proche de 10 v(Tm) ' 2. T, est de la forme T1=7.1Ein_,Com(xi), avec Com un réseau de neurones comportant une unique couche cachée de neurone(s) et où seules les connexions de sortie sont paramétrées. On notera par ailleurs que l'invention n'est pas limitée au mode de réalisation 15 décrit précédemment. Il apparaîtra en effet à l'homme de l'art que diverses modifications peuvent être apportées au mode de réalisation décrit ci-dessus, à la lumière de l'enseignement qui vient de lui être divulgué. Dans la présentation détaillée de l'invention qui est faite précédemment, les termes utilisés ne doivent pas être interprétés comme limitant l'invention aux modes de réalisation exposés dans la 20 présente description, mais doivent être interprétés pour y inclure tous les équivalents dont la prévision est à la portée de l'homme de l'art en appliquant ses connaissances générales à la mise en oeuvre de l'enseignement qui vient de lui être divulgué. c(E,Ti) REVENDICATIONS1. Procédé (300) de rendu d'une scène tridimensionnelle (200), comportant, pour chacun d'au moins un pixel (p) d'une image (206) représentant une vue de la scène tridimensionnelle (200) : - la détermination (306), à partir de la scène tridimensionnelle (200) et de paramètres de la vue, d'une fonction de contribution donnant une contribution d'un chemin de lumière partant d'une source de lumière (204) de la scène tridimensionnelle (200), à une valeur du pixel (p), une première sélection (318) aléatoire de premiers chemins de lumière, et l'évaluation (320) d'un estimateur d'une intégrale de la fonction de contribution, l'estimateur comportant un estimateur de Monte Carlo utilisant les premiers chemins de lumière, caractérisé en ce qu'il comporte en outre : - la détermination (314) d'au moins un paramètre d'une fonction de contribution approchée, de sorte que la fonction de contribution approchée donne une contribution d'un chemin de lumière à la valeur du pixel (p) proche de la contribution donnée par la fonction de contribution, la fonction de contribution approchée comportant un réseau de neurones comportant : - une couche d'entrée d'au moins un neurone, - une unique couche cachée d'au moins un neurone, - une couche de sortie d'au moins un neurone - des connexions entre le ou les neurones de la couche d'entrée et le ou les neurones de la couche cachée auxquelles des poids sont respectivement associés, ces poids étant indépendants du ou des paramètres, et - des connexions entre le ou les neurones de la couche cachée et le ou les neurones de la couche de sortie auxquelles des poids dépendant du ou des paramètres sont respectivement associés, et en ce que l'estimateur comporte en outre une moyenne de contributions données par la fonction de contribution approchée appliquée à au moins une partie des premiers chemins de lumière. 15 2. Procédé (300) selon la revendication 1, comportant en outre une deuxième sélection (312) aléatoire, indépendante de la première sélection (318), de deuxièmes chemins de lumière, et dans lequel le ou les paramètres de la fonction de contribution approchée sont déterminés (314) à partir des deuxièmes chemins de 5 lumière. 3. Procédé (300) selon la revendication 1 ou 2, dans lequel la détermination (314) du ou des paramètres est réalisée de manière à minimiser une distance entre la fonction de contribution approchée et la fonction de contribution. 4. Procédé (300) selon les revendications 2 et 3 prises ensemble, dans 10 lequel la distance utilise les deuxièmes chemins de lumière. 5. Procédé (300) selon la revendication 4, dans lequel la distance comporte la somme des carrés des différences entre la fonction de contribution approchée et la fonction de contribution pour les deuxièmes chemins de lumière. 6. Procédé (300) selon la revendication 5, dans lequel la distance comporte 15 en outre un paramètre de régularisation et dans lequel la détermination (314) du ou des paramètres comporte la détermination du paramètre de régularisation, puis la détermination du ou des paramètres en utilisant le paramètre de régularisation déterminé. 7. Procédé (300) selon la revendication 6, dans lequel l'estimateur peut 20 s'écrire matriciellement sous la forme A[Im- M]l3, où 4, est une matrice identité de taille m, A et B sont des vecteurs et M est une matrice dépendant du paramètre de régularisation y, et dans lequel le paramètre de régularisation y est déterminé de sorte que la valeur propre maximale de M présente au plus un écart prédéfini avec 1. 8. Procédé (300) selon l'une quelconque des revendications 1 à 7, dans 25 lequel les poids associés aux connexions entre le ou les neurones de la couche d'entrée et le ou les neurones de la couche cachée sont déterminés (310) de manière aléatoires. 9. Procédé (300) selon l'une quelconque des revendications 1 à 8, dans lequel l'estimateur comporte en outre une espérance de la fonction de contribution 30 approchée. 10. Programme d'ordinateur (106) téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, caractérisé en ce qu'il comprend des instructions pour l'exécution des étapes d'un procédé de rendu d'une scène tridimensionnelle (200) selon l'une 16 quelconque des revendications 1 à 9, lorsque ledit programme est exécuté sur un ordinateur. 11. Dispositif (100) de rendu d'une scène tridimensionnelle (200), comportant, pour chacun d'au moins un pixel (p) d'une image (206) représentant une vue de la 5 scène tridimensionnelle (200) : - des moyens de détermination, à partir de la scène tridimensionnelle (200) et de paramètres de la vue, d'une fonction de contribution donnant une contribution d'un chemin de lumière partant d'une source de lumière (204) de la scène tridimensionnelle (200), à une 10 valeur du pixel (p), des moyens de réalisation d'une première sélection aléatoire de premiers chemins de lumière, et des moyens d'évaluation d'un estimateur d'une intégrale de la fonction de contribution, l'estimateur comportant un estimateur de 15 Monte Carlo utilisant les premiers chemins de lumière, caractérisé en ce qu'il comporte en outre : - des moyens de détermination d'au moins un paramètre d'une fonction de contribution approchée, de sorte que la fonction de contribution approchée donne une contribution d'un chemin de 20 lumière à la valeur du pixel (p) proche de la contribution donnée par la fonction de contribution, la fonction de contribution approchée comportant un réseau de neurones comportant : - une couche d'entrée d'au moins un neurone, - une unique couche cachée d'au moins un neurone, 25 - une couche de sortie d'au moins un neurone - des connexions entre le ou les neurones de la couche d'entrée et le ou les neurones de la couche cachée auxquelles des poids sont respectivement associés, ces poids étant indépendants du ou des paramètres, et 30 - des connexions entre le ou les neurones de la couche cachée et le ou les neurones de la couche de sortie auxquelles des poids dépendant du ou des paramètres sont respectivement associés, et en ce que l'estimateur comporte en outre une moyenne de contributions données par la fonction de contribution approchée appliquée à au moins une partie des premiers chemins de lumière.