Volume Rendering - Introduction. Markus Hadwiger Visual Computing Center King Abdullah University of Science and Technology

Transcription

1 Volume Rendering - Introduction Markus Hadwiger Visual Computing Center King Abdullah University of Science and Technology

2 Volume Visualization 2D visualization: slice images (or multi-planar reformation: MPR) Indirect 3D visualization: isosurfaces (or surface-shaded display: SSD) Direct 3D visualization: direct volume rendering: DVR

3 Motivation

4 Motivation

5 Motivation

6 Direct Volume Rendering

7 Ray-Casting Image order approach: Image Plane Data Set Eye For each pixel { calculate color of the pixel }

8 Decrease Increase Physical Model of Radiative Transfer true emission in-scattering true absorption out-scattering

9 Volume Rendering Integral How do we determine the radiant energy along the ray? Physical model: emission and absorption, no scattering Initial intensity at s 0 Absorption along the ray segment s 0 - s viewing ray Optical depth τ Absorption к Without absorption all the initial radiant energy would reach the point s.

10 Volume Rendering Integral How do we determine the radiant energy along the ray? Physical model: emission and absorption, no scattering viewing ray Every One point along the the viewing ray emits additional radiant energy. Active emission at point s ~ Absorption along the distance s - ~ s

11 Numerical Solution

12 Numerical Solution can be computed recursively Radiant energy observed at position i Radiant energy emitted at position i Absorption at position i Radiant energy observed at position i 1

13 Numerical Solution back-to-front compositing Early Ray Termination: Stop the calculation when iterate from i=0 (back) to i=max (front): i increases can be computed recursively front-to-back compositing (used in ray-casting) iterate from i=max (front) to i=0 (back): i decreases

14 Ray-Casting: Implementation Ray setup Loop over ray Resample scalar value Classification Shading Compositing

15 Ray-Casting: Implementation Ray setup Loop over ray Resample scalar value Classification Shading Compositing

16 Rasterization-Based Ray Setup Pixel/fragment == ray Need ray start position and direction vector Rasterize bounding box - = Identical for orthogonal and perspective projection!

17 Object-Order Empty Space Skipping Modify initial rasterization step rasterize bounding box rasterize tight" bounding geometry

18 Ray-Casting: Implementation Ray setup Loop over ray Resample scalar value Classification Shading Compositing

19 Classification: Transfer Functions Assign optical properties: define the look of the data Which parts are transparent? Which parts have what color?

20 1D Transfer Functions texture = scalar field scalar value S resampling T(S) transferfunction texture = [Emission RGB, Absorption A] S RGBA

21 1D Transfer Functions

22 Pre-Classification / Post-Classification Pre-interpolative vs. post-interpolative classification pre-classification post-classification same TF, same resolution, same sampling rate

23 Pre-Classification / Post-Classification Pre-interpolative vs. post-interpolative classification pre-classification post-classification same TF, same resolution, same sampling rate

24 Pre-Integrated Classification front scalar back scalar pre-integrate all possible combinations in the TF s f s b store integral into table d constant sampling distance d for pre-integration of RGB, A s f s b s b s f

25 Pre-Integrated Classification

26 2D Transfer Functions 1D transfer function Horizontal axis: scalar value Vertical axis: number of voxels 1D Histogram 2D transfer function Horizontal axis: scalar value Vertical axis: gradient magnitude 2D Scatterplot [Kniss et al. 2002]

27 2D Transfer Functions 1D transfer function Horizontal axis: scalar value Vertical axis: number of voxels 2D transfer function Horizontal axis: scalar value Vertical axis: gradient magnitude [Kniss et al. 2002]

28 Ray-Casting: Implementation Ray setup Loop over ray Resample scalar value Classification Shading Compositing

29 Volume Shading Local illumination vs. global illumination Gradient-based or gradient-less Shadows, (multiple) scattering,

30 Ray-Casting: Implementation Ray setup Loop over ray Resample scalar value Classification Shading Compositing

31 Compositing

32 Compositing

33 Large Volumes: Bricking Subdivide volume into bricks Octrees, page tables, fixed grid, Combine memory management and empty space skipping Memory management: large bricks Empty space skipping: small bricks Cull using min-max of bricks Only store needed bricks in cache texture: swap in/out on demand

34 Thank You! Thank you: Christof Rezk-Salama, Klaus Engel, Johanna Beyer, Henning Scharsach, Stefan Wolfsberger, Christian Sigg, Katja Bühler, Joe Kniss