기술 요약: 얇은 기하 구조에서의 디폴레 확산 오차 (Dipole Diffusion Error in Thin Geometry)

문제 정의

컴퓨터 그래픽스의 하부 표면 산란(Subsurface Scattering, SSS)은 주로 디폴레 근사(dipole approximation) 및 그 파생 모델들(multipole, quantized diffusion 등)을 사용하여 모델링됩니다. 이러한 모델들은 **평평하고 반무한한 슬래브(flat, semi-infinite slab)**를 위해 유도된 방사형 반사 프로파일 $R_d(r)$을 가정합니다. 두꺼운 평면 구조에서는 효과적이지만, 이 가정은 얇거나 급격히 휘어진 물체(예: 잎의 가장자리, 보석의 모서리)에 적용될 때 체계적인 기하학적 오차를 유발합니다. 이러한 기하 구조에서는 빛이 디폴레 모델이 예측하는 것처럼 재방출되는 대신, 확산 거리 내에서 두 번째 경계와 마주치며 물체 밖으로 탈출하게 됩니다. 이는 얇은 특징 부분이 물리적으로 보여야 하는 것보다 더 밝게 보이는 "과도한 광휘(over-glow)" 아티팩트를 초래합니다.

멀티폴(multipole) 모델과 같은 기존 솔루션은 평면형 얇은 슬래브 문제는 해결하지만, 일반적인 곡면 기하 구조에는 실패합니다. 또한, 일반적인 기하 인식형 확산 추정기(geometry-aware diffusion estimator)는 내부 메쉬나 그리드를 필요로 하는데, 이는 계산 비용이 많이 들고 부호 거리 함수(Signed Distance Function, SDF) 기반 워크플로우에 통합하기 어렵습니다. 본 논문의 목적은 다음과 같습니다:

1. 디폴레 모델의 얇은 기하 구조에 대한 오차를 광학적 두께의 스칼라 함수로 특성화한다.
2. 내부 이산화 없이 SDF 상에서 직접 작동하는, 그리드 프리(grid-free) 방식의 기하학적 정밀 확산 추정기를 제공한다.
3. 렌더링, 오차 예측 및 역문제(inverse problems)에서의 유용성을 입증한다.

방법론

1. 슬래브 오차의 이론적 분석

저자들은 동일한 확산 상수와 경계 조건 하에서 표준 디폴레(반무한 슬래브)를 멀티폴(유한 슬래브)과 비교함으로써 기하학적 오차를 분리해 냅니다.

• 광학적 두께 법칙: 저자들은 상대적 확산 알베도(albedo)가 광학적 두께 $\tau = T/\ell_d$ (여기서 $T$는 물리적 두께, $\ell_d$는 확산 길이)에 따라 지수적으로 감소함을 유도했습니다.
• 폐쇄형 비율(Closed-Form Rate): 주요 오차 항은 $C e^{-2\tau}$이며, 여기서 비율 $2$는 슬래브를 통과하는 빛의 "왕복(round-trip)"을 의미하며, 전계 상수 $C$는 재질에 따라 달라집니다.
• 투과율(Transmittance): 얇은 슬래브에서 손실된 빛(디폴레가 고려하지 못한 부분)은 뒷면을 통해 투과됩니다. 이 투과율은 단일 패스 비율인 $e^{-\tau}$에 따라 감소합니다.
• 일반화: 공간적으로 변화하는 매질의 경우, 저자들은 최단 비용 경로를 나타내는 광학적 거리(optical distance, Agmon distance) $\tau^*(x)$ 개념을 사용하여 이를 확장했습니다. 이는 비균질 매질에 대해 슬래브 법칙을 일반화합니다.

2. 그리드 프리 기하학적 정밀 추정기

임의의 형상에 대해 메쉬 없이 확산 방정식을 평가하기 위해, 저자들은 스크린드 포아송 방정식(screened Poisson equation)을 풀기 위해 구체 위 걷기(Walk on Spheres, WoS) 기법을 사용합니다.
$$\nabla^2 \phi - c \phi = -f$$

• SDF 표현: 물체는 오직 부호 거리 함수(SDF)로 정의됩니다. WoS 알고리즘은 도메인 내부에 들어가는 가장 큰 구(ball) 내부의 점들을 샘플링하며, 거리 필드를 사용하여 구의 반지름과 경계 근접도를 결정합니다.
• 반평면 가정 배제: 국소 접평면 반평면을 가정하는 디폴레와 달리, WoS 워크는 뒷면과 타이트한 곡률을 포함한 실제 기하 구조를 존중합니다.
• 단일 산란(Single Scattering): 확산 모델은 다중 산란만을 다루므로, 저자들은 Henyey-Greenstein 위상 함수로부터 유도된 별도의 폐쇄형 단일 산란 항을 추가합니다. 이는 날카로운 근접장 피크를 포착하기 위해 짧은 레이 마칭(ray-march)을 통해 계산됩니다.
• 구현: 솔버는 CPU 참조 모델과 함께 순수 float Metal(GPU)로 구현되었으며, 균질 및 공간 가변 매질(null-scattering walk를 통해)을 모두 지원합니다.

3. 오차 예측 및 하이브리드 렌더링

**형상 직경 함수(Shape Diameter Function, 국소 두께)**를 $\ell_d$로 나눈 값을 사용하는 실용적인 예측기가 유도되었습니다.

• 보정(Correction): 이 예측기는 알베도 손실을 추정하고 표준 디폴레에 보정 계수 $\kappa(\tau)$를 적용합니다.
• 하이브리드 디스패치(Hybrid Dispatch): 예측기는 디폴레 오차가 유의미한 픽셀(예: 역광을 받는 얇은 영역)을 식별하는 데 사용됩니다. 이러한 픽셀은 비용이 많이 드는 그리드 프리 WoS 추정기로 렌더링되는 반면, 나머지는 보정된 디폴레를 사용합니다.

4. 역문제 (확산 광학 단층 촬영)

포워드 솔버가 미분 가능하기 때문에, 저자들은 미분 가능한 포워드 맵을 사용하여 은닉된 흡수 포함물(absorbing inclusion)의 위치와 대비를 복구하는 그리드 프리 역문제를 시연합니다.

주요 기여

1. 슬래브 법칙 및 두께 예측기:

   • 디폴레의 알베도 오차가 평면 케이스에서 재질과 무관한 지수적 비율 $e^{-2\tau}$를 따른다는 것을 증명했습니다.
   • 얇은 특징에 대해 알베도 오차를 크게 줄이는 폐쇄형 보정 $\kappa(\tau)$를 제공했습니다.
   • 브루트 포스 볼륨 패스 트레이싱을 사용하여 예측된 비율(1.99 vs 2.0)을 검증했습니다.

2. 기하학적 정밀 그리드 프리 추정기:

   • 내부 메쉬 없이 SDF 상의 스크린드 포아송 방정식을 위한 WoS 솔버를 구현했습니다.
   • 폐쇄형 슬래브 솔루션에 대해 0.75%의 오차를 달성했으며, 패스 트레이싱 참조 모델과 비교했을 때 디폴레보다 오차를 크게 줄였습니다.
   • 공간 가변 매질과 복잡한 경계 조건을 성공적으로 처리합니다.

3. 종합 검증:

   • 역광, 전광, 측면 조명 렌즈, 토러스의 네 가지 케이스에 대해 패스 트레이싱 참조 모델과 비교했습니다.
   • 얇은 두께가 지배적인 역광 케이스에서는 그리드 프리 추정기가 개선 효과를 보이지만, 단일 산란이나 각운동 수송이 지배적인 전광 또는 고곡률 시나리오에서는 디폴레보다 반드시 우월하지는 않음을 발견했습니다.

4. 미분 가능한 역 솔버:

   • 미분 가능한 포워드 맵을 사용하여 물체 반지름의 5.4% 이내의 중심점을 찾아내며 은닉된 가우시안 포함물을 높은 정밀도로 복구하는 것을 입증했습니다.

결과

• 정확도: 그리드 프리 추정기는 폐쇄형 슬래브 솔루션을 **0.75%**까지 일치시키며 멀티폴 솔루션을 정확히 재현합니다. 패스 트레이싱 참조 모델에 대해서는 역광 렌즈 시나리오에서 평균 원시 오차를 72%에서 45%로 줄였습니다.
• 성능: Apple M4 Pro에서 GPU 커널은 초당 4.1 $\times$ 10$^8$ 개의 볼 워크(ball walks)를 유지합니다.
• 하이브리드 효율성: (그리드 프리 솔버를 플래그된 역광 픽셀에만 사용하는) 하이브리드 접근 방식은 순수 디폴레(0.180)나 순수 그리드 프리(0.143)보다 낮은 평균 톤맵 오차(0.132)를 달성했습니다.
• 한계점 식별:
  • 곡률(Curvature): 일반적인 기하 구조에서 오차 전계는 두께뿐만 아니라 평균 곡률($\ell_d H$)에 의존합니다. 두께만으로는 오차 분산의 38%만을 설명하지만, 곡률을 추가하면 이를 94%까지 높입니다.
  • 모델 오차: 패스 트레이싱 비교에서의 나머지 오차는 디폴레 자체의 확산 근사(vs 복사 수송)에 의해 지배되며, 특히 $\tau < 0.4$인 경우에 그러합니다.
  • 비용: 높은 대비를 가진 공간 가변 매질은 몬테카를로 노이즈를 증가시키며, 비정밀한 SDF(느슨한 거리 필드)는 워크 스텝 수를 크게 증가시킵니다($\alpha^{-2.1}$).

의의 및 주장

본 논문은 얇은 특징에 대한 디폴레 근사의 특정 기하학적 오차를 정량화하고 완화한다고 주장합니다. 이는 복사 수송(radiative transport)을 대체하거나 모든 SSS 문제를 해결하겠다는 것이 아닙니다.

• 제한된 범위: 저자들은 이 방법이 확산 영역 내에서의 기하학적 오차를 줄이지만, 단일 산란이 지배적인 매우 얇거나 광학적으로 얇은 물체에 대한 수송 모델을 대체하지는 않는다고 명시합니다.
• 실용적 유용성: 유도된 두께 예측기는 (하이브리드 렌더링을 위한) 컴퓨팅 자원 할당을 위한 유용한 휴리스틱 역할을 하며, 전체 패스 트레이싱 없이도 프로덕션 파이프라인에서 알베도를 보정할 수 있게 합니다.
• 기초적 통찰: 본 연구는 "디폴레 실패"가 무작위적인 아티팩트가 아니라, 무작위 워크가 두 번째 경계에 부딪힐 확률에 의해 제어되는 광학적 거리에 따른 예측 가능한 함수임을 확립합니다.
• 역 역량: 그리드 프리 솔버의 미분 가능성은 그리드 프리 확산 광학 단층 촬영으로 가는 길을 열어주며, 이는 프로파일 기반 모델에서는 불가능한 기능입니다.

결론적으로, 디폴레가 효율성 면에서 표준으로 남아있지만, 기하학적 오차는 체계적이며 광학적 두께를 통해 보정 가능하다는 것을 보여줍니다. 또한 그리드 프리 확산 솔버가 특정 영역에서 기하학적으로 정확한 하부 표면 렌더링을 위한 실행 가능한 대안이 될 수 있음을 제시합니다.