Volume Illumination & Vector Field Visualisation

Transcription

1 Volume Illumination & Vector Field Visualisation Visualisation Lecture 11 Institute for Perception, Action & Behaviour School of Informatics Volume Illumination & Vector Vis. 1

2 Previously : Volume Rendering Image Order Volume Rendering Object Order Volume Rendering ray casting / intensity transfer function splatting / hybrid shear-warp rendering algorithm Volume Classification volume scalars colour & opacity feature isolation & classification transparency & clipping Volume Illumination & Vector Vis. 2

3 Light Propagation in Volumes Lighting in volume only transmission and emission considered (until now) can also: reflect light scatter light into different directions Volume Illumination & Vector Vis. 3

4 Global Illumination of Volumes For every voxel ray intersects, need to consider: Light absorbed. Light emitted. Light scattered out of the ray. Light scattered into the ray. Volume Illumination & Vector Vis. 4

5 Global Illumination of Volumes For every voxel ray intersects, need to consider: Light absorbed. Normally ignore scattering in volumetric illumination Light!Light emitted. scattered out of the ray. Light scattered into the ray. Why? : computational cost Volume Illumination & Vector Vis. 5

6 Example: single scatter Synthetic images of the Earth s atmosphere produced by Rayleigh scattering. Irwin 95. Light scattered to produce atmospheric haze effect Volume Illumination & Vector Vis. 6

7 Example : multiple scattering Scattering has scale dependence can see this is quite a small cloud. Light scattered multiple times to produce simulation of a cloud Volume Illumination & Vector Vis. 7

8 Example : sub-surface scattering Very large statue Medium sized Small statue. Varying how light is scattered inside a surface affects perception hence useful in visualisation why? - prior visual experience, perception of distance? Volume Illumination & Vector Vis. 8

9 Volume Illumination -? Anyway, scattering is to costly so we usually do not take them into account when doing volume rendering But we still can add slight shadows to the volume by illuminating them Volume Illumination & Vector Vis. 9

10 Volume Illumination Why do we want to illuminate volumes? illumination helps us to better understand 3D structure displays visual cues to surface orientation highlight significant gradients within volume Volume Illumination & Vector Vis. 10

11 What are we illuminating? embedded (iso-) surface sharp gradients in opacity Volume Illumination & Vector Vis. 11

12 Shading an Embedded isosurface classify volume with a step function use regular specular / diffuse surface shading Remember for lighting equations of lecture 2 require illumination direction camera model (position) surface orientation need to calculate and store surface normal (for our voxel in order to illuminate them) Volume Illumination & Vector Vis. 12

13 Estimating the surface normal from distance map Use distance map to the iso-surface value 1. Determine the threshold value 2. Determine the surface voxels based on the threshold 3. Compute the normal vectors based on centred difference method 1 zthe x = z x 1, y z x 1, y For example, if we sample 2 centre of the voxels, 1 z y = z x, y 1 z x, y 1 2 z x z y 1 N = 2, 2, z x z y 1 z x z y 1 z x z y 1 Volume Illumination & Vector Vis. 13

14 Result : illuminated iso-surface MIP technique Surface normals recovered from depth map of surface Shaded embedded isosurface. regions of similar curvature taken into consideration Volume Illumination & Vector Vis. 14

15 Normals are sensitive to step size Step size through the volume / over the depth map Step = 2.0 Step = 1.0 Step = 0.1 Artefacts with larger step sizes under standard lighting model Volume Illumination & Vector Vis. 15

16 Illuminating Opacity (Scalar) Gradient Illuminate scalar gradient instead of iso-surface requirement : estimate and store gradient at every voxel Composite Shaded Composite Shaded opacity gradient (shades changes in opacity) Volume Illumination & Vector Vis. 16

17 Estimating Opacity Gradient Use 3D centred difference operator I = I x, I y, I z = I, I, I x y z I x 1, y, z I x 1, y, z I= x 2 We can extract the normal vectors of the region where the scalar values are changing significantly, i.e. boundary of tissues Evaluate at each voxel and interpolate Volume Illumination & Vector Vis. 17

18 Illumination : storing normal vectors Visualisation is interactive compute normal vectors for surface/gradient once store normal perform interactive shading calculations Storage : 2563 data set of 1-byte scalars ~16Mb normal vector (stored as floating point(4-byte)) ~ 200Mb! Solution : quantise direction & magnitude as small number of bits Volume Illumination & Vector Vis. 18

19 Illumination : storing normal vectors Quantize vector direction into one of N directions on a sub-divided sphere Subdivide an octahedron into a sphere. Number the vertices. Encode the direction according the nearest vertex that the vector passes through. For infinite light sources, only need to calculate the shading values once and store these in a table. Volume Illumination & Vector Vis. 19

20 Volume Illumination Summary Illumination can aid visual understanding of volumes scattering gives clue to object scale / effects perception too expensive at current time for interactive visualisation Illumination of: embedded isosurface gradient magnitude finite difference estimation Normal vector storage for each cell need to estimate surface normal use depth map use index into triangulated sphere Volume Illumination & Vector Vis. 20

21 Visualising Vectors Examples of vector data: meteorological analyses / simulation medical blood flow measurement Computational simulation of flow over aircraft, ships, submarines etc. visualisation of derivatives Why is visualising these difficult? not just of flow itself 2 or 3 components per data point, temporal aspects of vector flow, vector density Volume Illumination & Vector Vis. 21

22 Insight in Vector Fields Two properties of vector fields to visualise : local view of the flow global view of the flow e.g a meteorological wind forecast Local : for given location, what is the current wind strength and direction Global: a given location, where has the wind flow come from, and where will it go to. Volume Illumination & Vector Vis. 22

23 Two Methods of Flow Visualisation Visualise Flow wrt fixed point e.g. plot flow glyphs to show local direction and magnitude local view of vector field Visualise flow as the trajectory of a particles transported by the flow e.g plot particle traces, streamlines etc. global view of vector field require integrating the flow equation Volume Illumination & Vector Vis. 23

24 State of Flow : Steady / Unsteady Steady flow remains constant over time state of equilibrium or snapshot use particle traces known as streamlines Unsteady flow varies with time implications to tracing massless particles particle traces known as streaklines show little information about flow direction or magnitude Volume Illumination & Vector Vis. 24

25 Vectors : local visualisation Set of basic methods for showing local view: oriented lines, hedgehogs & glyphs colour mapping vector components warping animation Volume Illumination & Vector Vis. 25

26 Local vector visualisation : lines Draw line at data point indicating vector direction scale according to magnitude indicate direction as vector orientation problems non-uniform spacing showing large dynamic range field e.g. speed Option : use barbs to show speed Volume Illumination & Vector Vis. 26

27 Example : meteorology NOAA/FSL Lines are drawn with constant length, barbs indicate wind speed. Also colour mapped scalar field of wind speed. Volume Illumination & Vector Vis. 27

28 Example : lines in 3D Problems : Difficult to understand position and orientation in projection to 2D image. Clutter is also a problem. Volume Illumination & Vector Vis. 28

29 Local vector visualisation : colour map Dense Normal Vectors in 3D capture of large scale environment X component = Red. Y component = Green Z component = Blue. Returned laser power Distance to object (darker is closer black is no data) Volume Illumination & Vector Vis. 29

30 Local vector visualisation : Glyphs 2D or 3D objects inserted at data point, oriented with vector flow problem : scaling scaling glyph results in non-linear change in appearance surface area changes with square of size problem : clutter e.g. blood flow (reduced data) colourmap shows magnitude in addition to glyph scale Volume Illumination & Vector Vis. 30

31 Local vector visualisation : Warping deform geometry according to the vector field vector fields often associated with motion, or displacement. e.g vibration of a beam. Volume Illumination & Vector Vis. 31

32 Example : warping Insert slice planes into the data volume Displace surface according to flow momentum take care with scaling to avoid excessive geometric distortion surfaces may intersect, or even turn inside-out Volume Illumination & Vector Vis. 32

33 Local vector visualisation : animation Animation to enhance lines or glyphs improved clarity of magnitude and/or direction draw lines or glyphs & animate over time removes ambiguity in line or glyph direction also move glyphs along a streamline to visualise transport Volume Illumination & Vector Vis. 33

34 Vector Visualisation Summary Vector visualisation: local view / global view steady / unsteady flow methods of vector visualisation : Local Vector Visualisation: lines, hedgehogs & glyphs colour mapping, warping & animation Next lecture : more vector fields! Volume Illumination & Vector Vis. 34