Résumé Technique : Erreur de Diffusion Dipolaire en Géométrie Fine

Énoncé du Problème

La diffusion de la lumière sous-surfacique (SSS) est principalement modélisée dans l'infographie par l'approximation dipolaire et ses descendants (multipolaire, diffusion quantifiée, etc.). Ces modèles supposent un profil de réflectance radiale $R_d(r)$ dérivé pour une plaque semi-infinie et plane. Bien qu'efficace pour des surfaces planes et épaisses, cette hypothèse introduit une erreur géométrique systématique lorsqu'elle est appliquée à des objets fins ou fortement courbés (ex: bords de feuilles, arêtes de gemmes). Dans de telles géométries, la lumière rencontre une seconde frontière avant une longueur de diffusion, s'échappant de l'objet plutôt que d'être réémise comme le prédit le modèle dipolaire. Cela entraîne des artefacts de « sur-éclat » (over-glow) où les éléments fins paraissent plus brillants qu'ils ne devraient physiquement l'être.

Les solutions existantes, telles que le modèle multipolaire, traitent les plaques fines planes mais échouent sur les géométries courbes générales. De plus, les estimateurs de diffusion tenant compte de la géométrie générale nécessitent généralement un maillage ou une grille interne, ce qui est coûteux en calcul et difficile à intégrer dans les flux de travail basés sur les fonctions de distance signée (SDF). L'article vise à :

1. Caractériser l'erreur du modèle dipolaire sur la géométrie fine comme une fonction scalaire de l'épaisseur optique.
2. Fournir un estimateur de diffusion sans grille, géométriquement exact, opérant directement sur les SDF sans discrétisation intérieure.
3. Démontrer l'utilité de cet estimateur pour le rendu, la prédiction d'erreur et les problèmes inverses.

Méthodologie

1. Analyse Théorique de l'Erreur de Plaque

Les auteurs isolent l'erreur géométrique en comparant le dipôle standard (plaque semi-infinie) au multipôle (plaque finie) sous des constantes de diffusion et des conditions aux limites identiques.

• Loi de l'Épaisseur Optique : Ils démontrent que l'erreur relative de l'albédo diffusique décroît exponentiellement avec l'épaisseur optique $\tau = T/\ell_d$ (où $T$ est l'épaisseur physique et $\ell_d$ la longueur de diffusion).
• Taux en Forme Fermée : Le terme d'erreur principal est $C e^{-2\tau}$, où le taux $2$ correspond à un « aller-retour » de la lumière à travers la plaque, et le préfacteur $C$ dépend du matériau.
• Transmittance : La lumière perdue par la plaque fine (que le dipôle ne parvient pas à comptabiliser) est transmise par la face arrière. Cette transmittance décroît à un taux de passage unique de $e^{-\tau}$.
• Généralisation : Pour les milieux spatialement variables, les auteurs étendent cela en utilisant le concept de distance optique (distance d'Agmon), $\tau^*(x)$, qui représente le chemin de coût minimal vers la frontière lointaine. Cela généralise la loi de la plaque aux milieux non homogènes.

2. Estimateur Sans Grille Géométriquement Exact

Pour évaluer l'équation de diffusion sur des formes arbitraires sans maillage, les auteurs utilisent la méthode Walk on Spheres (WoS) pour résoudre l'équation de Poisson écrantée :
$$\nabla^2 \phi - c \phi = -f$$

• Représentation SDF : L'objet est défini uniquement par une fonction de distance signée. L'algorithme WoS échantillonne des points à l'intérieur de la plus grande sphère s'inscrivant dans le domaine, en utilisant le champ de distance pour déterminer le rayon de la sphère et la proximité de la frontière.
• Absence d'Hypothèse de Demi-Espace : Contrairement au dipôle, qui suppose un demi-espace tangent local, le WoS respecte la géométrie réelle, incluant les faces arrière et les courbures serrées.
• Diffusion Simple : Puisque les modèles de diffusion ne traitent que la diffusion multiple, les auteurs ajoutent un terme de diffusion simple distinct, dérivé de la fonction de phase de Henyey-Greenstein. Celui-ci est calculé via un court ray-marching pour capturer le pic étroit du champ proche.
• Implémentation : Le solveur est implémenté en Metal (GPU) en flottants simples avec une référence CPU, supportant les milieux homogènes et spatialement variables (via des marches à diffusion nulle).

3. Prédiction d'Erreur et Rendu Hybride

Un prédicteur pratique est dérivé en utilisant la Fonction de Diamètre de Forme (épaisseur locale) divisée par $\ell_d$.

• Correction : Ce prédicteur estime la perte d'albédo et applique un facteur de correction $\kappa(\tau)$ au dipôle standard.
• Dispatch Hybride : Le prédicteur est utilisé pour identifier les pixels où l'erreur du dipôle est significative (ex: régions fines rétroéclairées). Ces pixels sont rendus avec l'estimateur WoS coûteux, tandis que les autres utilisent le dipôle corrigé.

4. Problème Inverse (Tomographie Optique Diffuse)

Parce que le solveur direct est différentiable, les auteurs démontrent un problème inverse sans grille : la récupération de l'emplacement et du contraste d'une inclusion absorbante cachée à partir de mesures de bordure diffuses.

Contributions Clés

1. Loi de la Plaque et Prédicteur d'Épaisseur :

   • Ont prouvé que l'erreur d'albédo du dipôle suit un taux exponentiel indépendant du matériau $e^{-2\tau}$ dans le cas planaire.
   • Ont fourni une correction en forme fermée $\kappa(\tau)$ qui réduit significativement l'erreur d'albédo pour les éléments fins.
   • Ont validé le taux (1,99 contre 2,0 prédit) en utilisant un traceur de chemin volumétrique de force brute.

2. Estimateur Sans Grille Géométriquement Exact :

   • Implémentation d'un solveur WoS pour l'équation de Poisson écrantée sur des SDF, ne nécessitant aucun maillage intérieur.
   • Ont atteint une erreur de 0,75 % par rapport aux solutions de forme fermée et ont réduit l'erreur par rapport aux références de path-tracing par rapport au dipôle.
   • Gèrent avec succès les milieux spatialement variables et les conditions aux limites complexes.

3. Validation Complète :

   • Comparaison de la méthode avec une référence de path-tracing sur quatre cas (rétroéclairage, éclairage frontal, lentilles latérales et un tore).
   • Ont constaté que si l'estimateur sans grille améliore les cas de rétroéclairage régis par l'épaisseur, il ne surpasse pas universellement le dipôle dans les scénarios d'éclairage frontal ou de forte courbure où la diffusion simple ou le transport angulaire prédominent.

4. Solveur Inverse Différentiable :

   • Démonstration de la récupération d'une inclusion gaussienne cachée avec une haute précision (centre à moins de 5,4 % du rayon de l'objet) grâce à la carte directe différentiable.

Résultats

• Précision : L'estimateur sans grille correspond aux solutions de plaque de forme fermée à 0,75 % et reproduit exactement la solution multipolaire. Par rapport à une référence de path-tracing, il réduit l'erreur brute moyenne de 72 % à 45 % dans les scénarios de lentilles rétroéclairées.
• Performance : Sur un Apple M4 Pro, le noyau GPU soutient 4,1 $\times$ 10$^8$ marches de sphères par seconde.
• Efficacité Hybride : Une approche hybride (utilisant le solveur sans grille uniquement pour les pixels signalés comme rétroéclairés) atteint l'erreur de ton moyen la plus faible (0,132) comparée au dipôle pur (0,180) ou au sans-grille pur (0,143).
• Limites Identifiées :
  • Courbure : Sur une géométrie générale, le préfacteur d'erreur dépend de la courbure moyenne ($\ell_d H$), et non seulement de l'épaisseur. L'épaisseur seule n'explique que 38 % de la variance de l'erreur ; l'ajout de la courbure porte ce chiffre à 94 %.
  • Erreur de Modèle : L'erreur restante dans les comparaisons de path-tracing est dominée par l'approximation de diffusion elle-même (vs le transport radiatif), plutôt que par l'erreur géométrique du solveur, particulièrement pour $\tau < 0,4$.
  • Coût : Les milieux spatialement variables à haut contraste augmentent le bruit de Monte Carlo, et les SDF non exactes (champs de distance lâches) augmentent considérablement le nombre d'étapes de marche ($\alpha^{-2,1}$).

Signification et Revendications

L'article affirme quantifier et atténuer l'erreur géométrique spécifique inhérente à l'approximation dipolaire pour les éléments fins. Il ne prétend pas remplacer le transport radiatif ni résoudre tous les défis de la SSS.

• Portée Modeste : Les auteurs précisent explicitement que la méthode réduit l'erreur géométrique au sein du régime de diffusion, mais ne remplace pas le transport pour les objets très fins ou optiquement minces où la diffusion simple domine.
• Utilité Pratique : Le prédicteur d'épaisseur dérivé sert d'heuristique utile pour le dispatch de ressources de calcul (rendu hybride) ou pour corriger l'albédo dans les pipelines de production sans nécessendre un path-tracing complet.
• Aperçu Fondamental : Ce travail établit que l'« échec du dipôle » n'est pas un artefact aléatoire mais une fonction prévisible de la distance optique, régie par la probabilité que les marches aléatoires frappent une seconde frontière.
• Capacité Inverse : La différentiabilité du solveur sans grille ouvre la voie à la tomographie optique diffuse sans grille, une capacité non disponible avec les modèles basés sur des profils.

L'article conclut que bien que le dipôle reste un standard pour son efficacité, l'erreur géométrique est systématique et corrigeable via l'épaisseur optique, et que les solveurs de diffusion sans grille offrent une alternative viable et sans maillage pour un rendu de sous-surface géométriquement précis dans des régimes spécifiques.