Eurographics 2012, Cagliari, Italy
Tunç O. Aydın Disney Research Zurich, Switzerland
Piotr Didyk MPI Informatik, Germany
Diego Gutierrez Universidad de Zaragoza, Spain
Tobias Ritschel Télécom ParisTech/
MPI Informatik, Germany Alessandro Artusi
CaSToRC Cyprus Institute Francesco Banterle Visual Computing Lab ISTI-CNR, Italy
Elmar Eisemann Télécom ParisTech / CNRS-LTCI, France
Rafał Mantiuk University of Bangor, UK
Special thanks to:
Karol Myszkowski, MPI Informatik, Germany
Mapping images to target devices:
spatial, temporal, stereo, tone, and color
Eurographics 2012, Cagliari, Italy
Introduction: The Problem
Eurographics 2012, Cagliari, Italy
Introduction: Outline
• Dynamic Range and Color Retargeting (~70 mins):
•Rafał Mantiuk, Tobias Ritschel, and Alessandro Artusi
• Reverse/Inverse Tone Mapping (~45 mins) :
•Francesco Banterle
• Image Spatial Resolution Retargeting (~45 mins) :
•Diego Gutierrez
• Temporal Image Retargeting (~70 mins) :
•Tobias Ritschel and Elmar Eisemann
• Image and Video Quality Assessment(~70 mins) :
•Tunç O. Aydın
• Stereo Content Retargeting (~50 mins):
•Piotr Didyk
Eurographics 2012, Cagliari, Italy
1
Rafal Mantiuk
http://www.bangor.ac.uk/mantiuk/
Bangor University, UK
Research Institute of Visual Computing
Tone Mapping
Multidimensional image retargeting
Eurographics 2012, Cagliari, Italy
2
Learning outcomes
• What is tone-mapping?
• What problem(s) does it solve?
• Why is the problem so difficult?
• How do we perceive high dynamic range images?
• What are the major approaches to tone- mapping?
• How to choose a tone-mapping for a particular application?
Eurographics 2012, Cagliari, Italy
4
Tone-mapping problem
luminance range [cd/m2]
conventional display
simultaneously human vision
adapted
Tone mapping
Eurographics 2012, Cagliari, Italy
5
Question to the audience
• Who has never used a tone-mapping operator?
Each camera needs to tone-map a real-world captured light before it can be stored as a JPEG.
This is essentially the same process as tone- mapping, although knows as ‘color reproduction’
or ‘color processing’.
Eurographics 2012, Cagliari, Italy
7
Color space retargeting problem
7
Real-world
Display
Goal: map colors to a restricted color space
Eurographics 2012, Cagliari, Italy
8
Perceptual retargeting problem
8
Real-world
Display The eye adapted to
the real-world viewing conditions
The eye adapted to the display viewing conditions
Goal: match color appearance
Eurographics 2012, Cagliari, Italy
9
Tone Mapping?
• HDR ?
• Or something else ?
Eurographics 2012, Cagliari, Italy
10
What is tone-mapping?
Although tone-mapping may have different meanings, this course is about:
A) Transformation of an image from an unrestricted color gamut of real world or an abstract scene to the restricted color gamut of a device
B) Retargeting the perceptual appearance from one viewing conditions to another
Eurographics 2012, Cagliari, Italy
12
Input and output
• HDR
• (approximate) physical units
• luminance
• linear RGB
• scene-referred
• LDR (SDR)
• pixel values
• luma
• gamma corrected R’G’B’
• display referred
Tone mapping Tone mapping
Eurographics 2012, Cagliari, Italy
13
Luminous efficiency function
(weighting) Light spectrum (radiance)
Luminance
• Luminance – perceived brightness of light, adjusted for the sensitivity of the visual system to wavelengths
L
V
0L() V()d
Luminance
Eurographics 2012, Cagliari, Italy
14
Do HDR images contain luminance values?
• Not exactly, because:
• the combination of camera red, green and blue spectral sensitivity curves will not match the luminous efficiency function
• But they contain a good-enough approximation for most applications
• For multi-exposure camera capture the error in luminance measurements is 10-15%
Eurographics 2012, Cagliari, Italy
15
Sensitivity to luminance
• Weber-law – the just-noticeable difference is proportional to the magnitude of a stimulus
The smallest detectable luminance difference Background
(adapting) luminance
Constant
L
Typical stimuli:
ΔLErnst Heinrich Weber
[From wikipedia]Eurographics 2012, Cagliari, Italy
16
Consequence of the Weber-law
• Smallest detectable difference in luminance
• Adding or subtracting luminance will have different visual impact depending on the background luminance
• Unlike LDR luma values, HDR luminance values are not perceptually uniform!
L ΔL
100 cd/m
21 cd/m
21 cd/m
20.01 cd/m
2Eurographics 2012, Cagliari, Italy
17
How to make luminance (more) perceptually uniform?
• Using Fechnerian integration
R(L) 1
L(l) dl
0
Lluminance - L
response - R
1
ΔL
) ( ) 1
( L L L
dR
Derivative of response Derivative of
response
Detection threshold Detection threshold
Luminance transducer:
Eurographics 2012, Cagliari, Italy
18
Assuming the Weber law
• and given the luminance transducer
• the response of the visual system to light is:
R(L) 1
L(l) dl
0
LEurographics 2012, Cagliari, Italy
19
Fechner law
• Practical insight from the Fechner law:
• The easiest way to adopt image processing algorithms to HDR images is to convert luminance (radiance) values to the logarithmic domain
Gustav Fechner
[From Wikipedia]
R(L) aln(L)
Eurographics 2012, Cagliari, Italy
20
But…the Fechner law does not hold for the full luminance range
• Because the Weber law does not hold either
• Threshold vs. intensity function:
L ΔL
The Weber law region The Weber law
region
Eurographics 2012, Cagliari, Italy
21
Weber-law revisited
• If we allow detection threshold to vary with luminance according to the t.v.i. function:
• we can get more accurate estimate of the
“response”:
R(L) 1
L(l) dl
0
LL
ΔL
Eurographics 2012, Cagliari, Italy
22
Fechnerian integration and Steven’s law
22 R(L) - function
derived from the t.v.i. function R(L) - function derived from the
t.v.i. function
R(L) 1
L(l) dl
0
LEurographics 2012, Cagliari, Italy
24
Major approaches to tone-mapping
• Illumination & reflectance separation
• Forward visual model
• Forward & inverse visual models
• Constraint mapping problem
• This is not a crisp categorization
• Some operators combine several approaches
Eurographics 2012, Cagliari, Italy
25
Major approaches to tone-mapping
• Illumination & reflectance separation
• Forward visual model
• Forward & inverse visual model
• Constraint mapping problem
Eurographics 2012, Cagliari, Italy
26
Illumination &
reflectance separation
Input
Illumination
Reflectance
Eurographics 2012, Cagliari, Italy
27
Illumination and reflectance
Illumination
• Sun ≈ 10
9cd/m
2• Lowest perceivable luminance ≈ 10
-6cd/m
2• Dynamic range 10,000:1 or more
• Visual system partially discounts illumination
Reflectance
• White ≈ 90%
• Black ≈ 3%
• Dynamic range < 100:1
• Reflectance critical for object & shape detection
Eurographics 2012, Cagliari, Italy
28
Reflectance & Illumination TMO
• Distortions in reflectance are more apparent than the distortions in illumination.
• Tone mapping could preserve reflectance but compress illumination
• for example:
I d R T(I)
Tone-mapped image Tone-mapped
image
Reflectance Reflectance
Illumination Illumination
Tone-mapping Tone-mapping
I d R L 1/
Eurographics 2012, Cagliari, Italy
29
How to separate the two?
• (Incoming) illumination – slowly changing
• except very abrupt transitions on shadow boundaries
• Reflectance – low contrast and high frequency variations
Eurographics 2012, Cagliari, Italy
30
Gaussian filter
• First order approximation
• Blurs sharp boundaries
• Causes halos
Tone mapping result
Eurographics 2012, Cagliari, Italy
31
Bilateral filter
• Better preserves sharp edges
• Still some blurring on the edges
• Reflectance is not perfectly separated from illumination near edges
Tone mapping result
[Durand & Dorsey, SIGGRAPH 2002]
Eurographics 2012, Cagliari, Italy
32
WLS filter
• Weighted-least-squares optimization
• [Farbman et al., SIGGRAPH 2008]
-> min Make reconstructed image u
possibly close to input g Make reconstructed image u
possibly close to input g
Smooth out the image by making partial derivatives close to 0 Smooth out the image by making
partial derivatives close to 0
Spatially varying smoothing – less smoothing near the edges Spatially varying smoothing – less
smoothing near the edges
Eurographics 2012, Cagliari, Italy
33
WLS filter
• Stronger smoothing and still distinct edges
• Can produce stronger effects with fewer artifacts
Tone mapping result
Eurographics 2012, Cagliari, Italy
34
Retinex
• Retinex algorithm was initially intended to separate reflectance from illumination [Land 1964]
• There are many variations of Retinex, but the general principle is to eliminate from an image small gradients, which are attributed to the illumination
1 step: compute gradients in log domain
2
ndstep: set to 0 gradients less than the threshold
t
G
in
G
out3
rdstep: reconstruct an image from the vector field
For example by solving the Poisson equation
Eurographics 2012, Cagliari, Italy
35
Retinex examples
35
From: http://dragon.larc.nasa.gov/retinex/757/Original After Retinex
From:http://www.ipol.im/pub/algo/lmps_retinex_poisson_equation/#ref_1
Eurographics 2012, Cagliari, Italy
36
Gradient domain HDR compression
• Similarly to Retinex, it operates on log-gradients
• But the function amplifies small contrast instead of removing it
36 [Fattal et al., SIGGRAPH 2002]
Retinex Gradient domain
Contrast compression achieved by global contrast reduction
Enhance reflectance, then compress everything
Eurographics 2012, Cagliari, Italy
37
Contrast domain image processing
Original Image
Modified Image
Perceived contrast representation
Contrast enhancement
Rationale: Human eye is more sensitive to contrast than luminance Perceived contrast
representation
[Mantiuk et al., ACM Trans. Applied Perception, 2006]
Eurographics 2012, Cagliari, Italy
38
Contrast domain image processing
1
stlevel 2
ndlevel
Wavelets
Gradients
Contrast pyramid
1
stlevel 2
ndlevel
Image transform: Multi-scale contrast pyramid
Eurographics 2012, Cagliari, Italy
39
Contrast transducer function
Goal: Transform contrast to the representation that is possibly perceptually uniform.
Input Input
Output – approximates perceived contrastOutput – approximates perceived contrastEurographics 2012, Cagliari, Italy
40
Contrast Equalization: Examples
Log-Linear Scaling Contrast mapping
Contrast equalization
Eurographics 2012, Cagliari, Italy
41
Contrast Equalization: Examples
Log-Linear Scaling Contrast mapping Contrast equalization
Eurographics 2012, Cagliari, Italy
42
Tone mapping in photography
• Dodging and burning
• Darken on brighten image parts by occluding photographic paper during exposure
• Ansel Adams, The print, 1995
• Photoshop tool
• Essentially – attenuate low-pass frequencies associated to illumination
Eurographics 2012, Cagliari, Italy
43
Automatic dodging and burning
• Reinhard et al., Photographic tone reproduction for digital images.€ • SIGGRAPH 2002
• Choose dodging an burning kernel size adaptively
• depending on the response of the center-surround filter
• thus avoid halo artifacts
43
Eurographics 2012, Cagliari, Italy
44
Major approaches to tone-mapping
• Illumination & reflectance separation
• Forward visual model
• Forward & inverse visual model
• Constraint mapping problem
Eurographics 2012, Cagliari, Italy
45
Forward visual model
• Mimic the processing in the human visual system
• Assumption: what is displayed is brightness or abstract response of the visual system
Original image Original
image Visual model Visual model Displayed image Displayed
image Luminance,
radiance
Brightness, abstract response
Eurographics 2012, Cagliari, Italy
46
Forward visual model: Retinex
• Remove illumination component from an image
• Because the visual system also discounts illuminant
• Display ‘reflectance’ image on the screen
• Assumption:
• The abstract ‘reflectance’ contains most important visual information
• Illumination is a distraction for object recognition and
scene understanding
Eurographics 2012, Cagliari, Italy
47
Photoreceptor response
• Dynamic range reduction inspired by photoreceptor physiology
• [Reinhard & Devlin ‘05]
• From gamma to sigmoidal response:
Input luminance
Pixel value
Eurographics 2012, Cagliari, Italy
48
Results: photoreceptor TMO
48
Eurographics 2012, Cagliari, Italy
49
Photoreceptor models
• Naka-Rushton equation:
• Response of the photoreceptor to a short flicker of light - less applicable to viewing static images
Experiment:
time σ
L
Detection
Eurographics 2012, Cagliari, Italy
50
Sigmoidal tone-curves
• Very common in digital cameras
• Mimic the response of analog film
• Analog film has been engineered for many years to produce
optimum tone-reproduction (given that he tone curve must not change)
• Effectively the most commonly used tone- mapping!
Eurographics 2012, Cagliari, Italy
51
Why sigmoidal tone-curves work
• Because they mimic photoreceptor response
• Unlikely, because photoreceptor response to steady light is not sigmoidal
• Because they preserve contrast in mid-tones, which usually contains skin color
• We are very sensitive to variation in skin color
• Because an image on average has Gaussian distribution of log-luminance
• S-shape function is the result of histogram equalization of an image with a Gaussian-shape histogram
Eurographics 2012, Cagliari, Italy
52
Lightness perception
• Lightness perception in tone-reproduction for high dynamic range images [Krawczyk et al. ‘05]
• Based on Gilchrist lightness perception theory
• Perceived lightness is anchored to several
frameworks
Eurographics 2012, Cagliari, Italy
53
Gilchrist lightness perception theory
• Frameworks – areas of common illumination
• Anchoring – the tendency of
• highest luminance
• largest area
• to appear white
• Tone-mapping
• Rescale luminance in each framework to its anchor
Eurographics 2012, Cagliari, Italy
54
Results – lightness perception TMO
Eurographics 2012, Cagliari, Italy
55
Major approaches to tone-mapping
• Illumination & reflectance separation
• Forward visual model
• Forward & inverse visual model
• Constraint mapping problem
Eurographics 2012, Cagliari, Italy
56
Forward and inverse visual model
Original image Original
image Visual model Visual model
Inverse visual model Inverse visual
model Luminance, radiance
abstract response
Displayed image Displayed
image
Luminance, radiance
Editing (optional)
Editing (optional) World viewing
conditions
Display viewing conditions
Eurographics 2012, Cagliari, Italy
57
Contrast domain image processing
Original Image
Modified Image
Perceived contrast representation
Contrast enhancement
Rationale: Human eye is more sensitive to contrast than luminance Perceived contrast
representation [Mantiuk et al., ACM Trans. Applied Perception, 2006]
Eurographics 2012, Cagliari, Italy
58
Multi-scale model
• Multi-scale model of adaptation and spatial vision and color appearance
• [Pattanaik et al. ‘98]
• Combines
• psychophysical threshold and superthreshold visual models
• light & dark adaptation models
• Hunt’s color appearance model
• One of the most sophisticated visual models
Forward visual model Inverse visual model
Eurographics 2012, Cagliari, Italy
59
Results: multiscale model …
Eurographics 2012, Cagliari, Italy
60
Forward and inverse visual model
• Advantages of F&I visual models
• Can render images for different viewing conditions
• Different state of chromatic or luminance adaptation
• Physically plausible
• output in the units of luminance or radiance
• Shortcomings F&I visual models
• Assume that a standard display can reproduce the impression of viewing much brighter or darker scenes
• Cannot ensure that the resulting image is within the dynamic range of the display
• Not necessary meant to reduce the dynamic range
• Visual models are difficult to invert
Eurographics 2012, Cagliari, Italy
61
Major approaches to tone-mapping
• Illumination & reflectance separation
• Forward visual model
• Forward & inverse visual model
• Constraint mapping problem
Eurographics 2012, Cagliari, Italy
62
Constraint mapping problem
• Goal: to restrict the range of values while reducing inflicted damage
Eurographics 2012, Cagliari, Italy
63
Global tone mapping operator
Image histogram Image histogram
Best tone- mapping is the one which does not do anything, i.e.
slope of the tone-mapping curves is equal
to 1.
Best tone- mapping is the one which does not do anything, i.e.
slope of the tone-mapping curves is equal
to 1.
Eurographics 2012, Cagliari, Italy
64
Display limitations
But in practice contrast (slope) must be limited due to display limitations.
But in practice
contrast (slope)
must be limited
due to display
limitations.
Eurographics 2012, Cagliari, Italy
65
Tone mapping
Global tone- mapping is a compromise between clipping and
contrast compression.
Global tone- mapping is a compromise between clipping and
contrast compression.
Eurographics 2012, Cagliari, Italy
66
Histogram equalization
• 1. Compute cumulative distribution function:
• 2. Use that function to assign new pixel values
) (
inout
c Y
Y
Y
in YoutEurographics 2012, Cagliari, Italy
67
Histogram equalization
• Steepest slope for strongly represented bins
• Enhance contrast, if many pixels
• Reduce contrast, if few pixels
• HE distributes contrast distortions relative to the “importance” of a brightness level
Input log intensity
Input log intensity
FrequencyOutput log intensityEurographics 2012, Cagliari, Italy
68
Histogram adjustment with a linear ceiling
• [Larson et al. 1997, IEEE TVCG]
Linear mapping Histogram equalization
Histogram equalization with ceiling
Eurographics 2012, Cagliari, Italy
69
Histogram adjustment with a linear ceiling
• Truncate the bins that exceed the ceiling
• Recompute the ceiling based on the truncated histogram
• Repeat until converges
Input log intensity
Frequency
Ceiling, based on the detection thresholds of the
visual system Ceiling, based on
the detection thresholds of the
visual system
Eurographics 2012, Cagliari, Italy
70
Display adaptive tone-mapping
Goal: Minimize the visual difference between the input and displayed images
[Mantiuk et al., SIGGRAPH 2008]
Eurographics 2012, Cagliari, Italy
71
Display adaptive tone-mapping
[Mantiuk et al., SIGGRAPH 2008]
Forward-inverse visual model Forward-inverse
visual model
Eurographics 2012, Cagliari, Italy
72
Display adaptive TMO Non-adaptive TMO
10 300 10 000
lux
Results: ambient illumination compensation
Eurographics 2012, Cagliari, Italy
73
Display adaptive TMO Non-adaptive TMO
10 300 10 000
lux
Results: ambient illumination compensation
Eurographics 2012, Cagliari, Italy
74
Display adaptive TMO Non-adaptive TMO
10 300 10 000
lux
Results: ambient illumination compensation
Eurographics 2012, Cagliari, Italy
75
Results: display contrast
ePaper standard LCD HDR display
Eurographics 2012, Cagliari, Italy
76
Tone-mapping for video compression
• Find the tone-curve that minimizes distortion in a backward-compatible HDR video encoding
76
[Mai et al., IEEE TIP 2010]
Closed-form
solution:
Eurographics 2012, Cagliari, Italy
77
Which tone-mapping to choose?
• Illumination & reflectance separation
• Forward visual model
• Forward & inverse visual model
• Constraint mapping problem
77
1. Think what is the target application
- and thus the goal of your tone-mapping
2. Consider which tone-mapping approach(es) will deliver that goal
Eurographics 2012, Cagliari, Italy
78
Future of tone-mapping
Tone-mapping of today
• Built into cameras
• Assumes that all displays are the same
Tone-mapping of tomorrow
• Display tone-maps content on demand
• Depending on viewing conditions, viewer, its capabilities
• Content recorded, stored and transmitted in an HDR format
HDR- HDMI
Eurographics 2012, Cagliari, Italy
79
Thank you
Eurographics 2012, Cagliari, Italy
Apparent Contrast and Brightness Enhancement
Tobias Ritschel MPI Informatik Ramon Cajal Fellow
Eurographics 2012, Cagliari, Italy
Motivation
• Image display
• Dynamic range of existing displays is limited
• No reproduction real-world contrast/brightness
• Good image appearance doesn’t require that
• Modern tone mapping operators good at optimizing the physical contrast and luminance use
Eurographics 2012, Cagliari, Italy
Motivation
• Human preference
• Enhanced contrast and brightness improve image appearance
• Can we still boost the contrast and
brightness impression?
Eurographics 2012, Cagliari, Italy
Human perception
• Spatial vision
• Cornsweet illusion Apparent contrast boost
• Glare illusion
Apparent brightness boost
Eurographics 2012, Cagliari, Italy
Contrast Enhancement: Motivation
• Usual contrast enhancement techniques
• either enhance everything
• or require manual intervention
• change image appearance
• Tone mapping often gives numerically optimal solution
• no dynamic range left for enhancement
tone mapping result
HDR image (reference)
restore missing contrast
Krawczyk et al. EG2007
Eurographics 2012, Cagliari, Italy
Overview
Reference HDR Image Tone Mapped Image
Measure Lost Contrast
at Several Feature Scales
Enhance Lost Contrast in Tone Mapped Image
Enhanced TM Image communicate lost
image contents
maintain image appearance
Krawczyk et al. EG2007
Eurographics 2012, Cagliari, Italy
Cornsweet Illusion
Create apparent contrast based on Cornsweet illusion
Countershading
– gradual darkening / brightening towards a contrasting edge – contrast appears with ‘economic’ use of dynamic range
Enhanced Image
Krawczyk et al. EG2007
Eurographics 2012, Cagliari, Italy
Details of Contrast Illusion
WHAT YOU SEE ACTUAL SIGNAL
Krawczyk et al. EG2007
Eurographics 2012, Cagliari, Italy
Details of Contrast Illusion
• Luminance profiles cause contrast
• Properties:
• Shape matches shape of the enhanced feature
• Amplitude defines the perceived contrast
• Noise (texture) does not cancel the illusion
• Profiles should not be discernible
WHAT YOU SEE ACTUAL SIGNAL
Krawczyk et al. EG2007
Eurographics 2012, Cagliari, Italy
Construction of Simple Profile (1/2)
Profile from low- pass filtered reference
Size and amplitude adjusted manually
This is unsharp masking
SIGNAL(e.g. TM)
REFERENCE (e.g. HDR)
- +
REFERENCE RESTORED
low-pass filter
Krawczyk et al. EG2007
Eurographics 2012, Cagliari, Italy
Construction of Simple Profile (2/2)
• Well preserved signal is exaggerated by unsharp masking
SIGNAL (texture preserved)
REFERENCE (with texture)
- +
low-pass filter
Krawczyk et al. EG2007
Eurographics 2012, Cagliari, Italy
Multi-resolution Contrast Metric
Reference HDR Image Tone Mapped Image
Measure Lost Contrast
at Several Feature Scales
1 4 7
Contrast ratios
at several scales
Eurographics 2012, Cagliari, Italy
Link: Contrast Metric & Profiles
• Contrast defines the sub-band amplitude
• Contrast for larger scales appears also on smaller scales
• the full profile is always reconstructed (red)
• Scale of contrast defines the profile size
1 4 7
Eurographics 2012, Cagliari, Italy
Adaptive Countershading
final contrast restoration progress of restoration
Eurographics 2012, Cagliari, Italy
Adaptive Countershading
without visual model with visual model
Eurographics 2012, Cagliari, Italy
Restoration of TM Images (1/3)
Eurographics 2012, Cagliari, Italy
Restoration of TM Images (2/3)
reference HDR image (clipped) countershading of tone mapping
countershading profiles tone mapping
Eurographics 2012, Cagliari, Italy
Restoration of TM Images (3/3)
reference HDR image (clipped) countershading of tone mapping
countershading profiles tone mapping
Eurographics 2012, Cagliari, Italy
Countershading Variants
• Traditional countershading:
• Performed in the achromatic channel to enhance perceived luminance contrast
• Cross-modal approach:
• Use depth signal to derive countershading profile
• Countershading over chromatic channels enhances the overall image contrast
• Color2Grey:
• Dimensionality reduction 3->1: may lead to information loss
• Countershading in the achromatic channel used to reproduce lost chromatic contrast
• Disparity: Discussed later
Eurographics 2012, Cagliari, Italy
Purpose: Contrast Restoration
Reference Depth Map Input Image
Measure Missing Contrast
at Several Feature Scales
Enhance Missing Contrast in
The Input Image
Enhanced Image
Luft et al. SIG2008
Eurographics 2012, Cagliari, Italy
Depth Map as Contrast Reference
depth information original image
depth darkening [Luft2006]
adaptive countershading
Luft et al. SIG2008
Eurographics 2012, Cagliari, Italy
Colourfulness Countershading
• “Strasbourg”: Gradient method tone mapping, strong global contrast loss so strong restoration effect.
• Colourfulness contrast at border between sky and buildings
• promotes FG/BG separation
• creates impression of greater dynamic range
• increases impression of depth
Smith et al. EG2006
Eurographics 2012, Cagliari, Italy
Countershading Results (original)
Smith et al. EG2006
Eurographics 2012, Cagliari, Italy
Countershading Results (chroma enhancement)
Smith et al. EG2006
Eurographics 2012, Cagliari, Italy
Color2Grey Application
• Isoluminant color pattern transformed to grey G using Helmholz-Kohlraush effect, which takes into account the contribution of chromatic component into brightness
Smith et al. EG2008
Eurographics 2012, Cagliari, Italy
Color2Grey Application
Smith et al. EG2008
Eurographics 2012, Cagliari, Italy
Color2Grey Application
• G’L* : The effect of adding multi-resolution countershading correction hi(GL*) (upper- left) to the greyscale image GL* (lower-left)
The correction is driven by contrast in chroma channels of the original image I (upper-left)
Smith et al. EG2008
Eurographics 2012, Cagliari, Italy
Color2Grey Application
Smith et al. EG2008
Eurographics 2012, Cagliari, Italy
Countershading in 3D?
Cornsweet in 3D is More plausible → Less of an artefact → Stronger →
Better
D. Purves, A.
Shimpi, R. B. Lotto An empirical explanation of the Cornsweet effect.
J. of Neuroscience 19, 1999
Eurographics 2012, Cagliari, Italy
S. Dalí, Landscape with butterflies
Scene-aligned Countershading
Eurographics 2012, Cagliari, Italy
G. Seurat, Bathers at Asnieres
Scene-aligned Countershading
Eurographics 2012, Cagliari, Italy
Eurographics 2012, Cagliari, Italy
3D Unsharp Masking
• U(S)=S + λ(S – Sσ)
Lit 3D Surface
Contrast Signal
λ x -
Smoothly Lit 3D Surface
σ
Enhanced Image
+
Ritschel et al. SIG2008
Eurographics 2012, Cagliari, Italy
3D unsharp
masking
2D unsharp
masking Original image
3D blurred
signal Enhancement
signal
Mesh
3D Unsharp Masking
Eurographics 2012, Cagliari, Italy
Adjustable Effect
• U(S)=S + λ(S – Sσ)
St re n g th λ
Width σ
Eurographics 2012, Cagliari, Italy
2D vs. 3D Unsharp Masking Comparison
2D 3D
Signal
Image Lit Surface
Smoothing
(Gaussian) Image Blur Laplacian Surface Blur
RepresentationPixels Lit vertices and pixels
Smoothness σImage distance Geodesic world distance
Strength λ
Factor Factor
Eurographics 2012, Cagliari, Italy
3D Unsharp Masking: Scene Coherence
Perspective Perspective Occlusion
Occlusion
Culling Culling
Ritschel et al. SIG2008
Eurographics 2012, Cagliari, Italy
Complex Mesh
Original rendering 3D unsharp masked
rendering
Ritschel et al. SIG2008
Eurographics 2012, Cagliari, Italy
3D unsharp
masking
2D unsharp
masking Original image
3D blurred
signal Enhancement
signal
Mesh
Enhanced Text Contrast in the Shadow
Eurographics 2012, Cagliari, Italy
Results – Legibility
Eurographics 2012, Cagliari, Italy
Normal Enhancement
• Only geometric term
• Shadows ?
• Hightlights ?
• Reflectance ?
• Vertex resolution
• 3D unsharp masking:
Pixel resolution
Cignoni et al. ´05, C & G Vol. 29
Eurographics 2012, Cagliari, Italy
Exaggerated Shading
• Object enhancement
• Illuminate each vertex at grazing angle
• Improves geometry understanding
• Highlights?
• Shadows?
• Scene enhancement
• Change everything
• Both have applications Rusinkiewicz et al., SIGGRAPH‘06
Eurographics 2012, Cagliari, Italy
Specular shading
Eurographics 2012, Cagliari, Italy
Study
• Goals
• Find suitable settings
• See limitations
• Rank preference
• Method of adjustments
• Strength λ: adjustable
• Fixed width σ: low, medium, high
• 4 scenes,15 participants
• Task: Find such λ that:
• Added enhancement is just noticeable
• Added enhancement becomes objectionable
• Image appearance is preferred
Ihrke et al. SPIE2009
Eurographics 2012, Cagliari, Italy
Results
Eurographics 2012, Cagliari, Italy
Ihrke et al. SPIE2009
Results
Eurographics 2012, Cagliari, Italy
Results
• 2 JND
• preferred
• 4 JND
• objectionable
Ihrke et al. SPIE2009
Eurographics 2012, Cagliari, Italy
Countershading parameter effect
Original
Original Enhancement Enhancement Halo Halo
Unsharp masking, countershading and haloes: Enhancements or artifacts?
M.Trentacoste, R. Mantiuk, W. Heidrich, F. Dufrot Eurographics 2012
Unsharp masking, countershading and haloes: Enhancements or artifacts?
M.Trentacoste, R. Mantiuk, W. Heidrich, F. Dufrot Eurographics 2012
Eurographics 2012, Cagliari, Italy
Model of acceptable countershading
Objectionable (halos)
Indistinguishable
Eurographics 2012, Cagliari, Italy
Applications: Image Resizing
Eurographics 2012, Cagliari, Italy
Applications: Viewer-adaptive display
Eurographics 2012, Cagliari, Italy
Applications: Tone-mapping
Eurographics 2012, Cagliari, Italy
Summary
• Better communicate image contents with a minimal change to image appearance
• Application of Cornsweet illusion to image enhancement
• Generalization of unsharp masking
• Automatic enhancement given the reference data:
• HDR image
• depth information
• shading in 3D scene
• Scene consistent 3D unsharp masking leads to even stronger effects
Eurographics 2012, Cagliari, Italy
Glare Illusion [Zavagno and Caputo 2001]
Eurographics 2012, Cagliari, Italy
Glare Illusion
“Alan Wake”© Remedy Entertainment
Eurographics 2012, Cagliari, Italy
Glare Illusion in Different Media
Arts
Photography Computer games
Eurographics 2012, Cagliari, Italy
In Games
• Simple approximation:
convolution with Gaussian
• Already does a good job in conveying brightness Yoshida et al.
(2008)
Eurographics 2012, Cagliari, Italy
In Games
+ +
+ =
+
Kawase: Practical Implementation of High Dynamic Range Rendering. Game Developer’s Conference 2004
Eurographics 2012, Cagliari, Italy
Glare in Realistic Rendering
• Optics-based models for rendering glare illusion
• [Nakamae et al. 1990]
• [Ward Larson et al. 1997]
• [Kakimoto et al. 2004, 2005]
• [Van den Berg et al. 2005]
• [Spencer et al. 1995]
Eurographics 2012, Cagliari, Italy
Dynamic Glare
• Realism
• Movement
• Colors
• Required Model of dynamic human eye to simulate temporal glare
• Study
Can temporal glare boost even further boost brightness?
Ritschel et al. EG2008
Eurographics 2012, Cagliari, Italy
Point Spread Function (PSF)
• Point Spread Function
• Key to glare modeling
• Describes, how a pixel maps to a pattern under an aperture
Ritschel et al. EG2008
Eurographics 2012, Cagliari, Italy
FFT
Our Simplified Model
Tissue
Perceived image Convolution Aperture
plane
Simplification: Fresnel diffraction Real world
Ritschel et al. EG2008
Eurographics 2012, Cagliari, Italy
Diffraction: Single vs. Multi Aperture Planes
Ritschel et al. EG2008
Eurographics 2012, Cagliari, Italy
Diffraction: Fraunhofer vs. Fresnel
Ritschel et al. EG2008
Eurographics 2012, Cagliari, Italy
Temporal Glare Pipeline
Ritschel et al. EG2008
Eurographics 2012, Cagliari, Italy
Aperture: Pupil
Eurographics 2012, Cagliari, Italy
Aperture: Pupil
• Adaptation
• Can convert HDR image into pupil size
• Pupillary hippus
DarkerBrighter
Big pupil Small pupil
Eurographics 2012, Cagliari, Italy
Aperture: Pupil
Eurographics 2012, Cagliari, Italy
Aperture: Lens
Eurographics 2012, Cagliari, Italy
Aperture: Lens
Eurographics 2012, Cagliari, Italy
Aperture: Gratings / Lens fibers
Eurographics 2012, Cagliari, Italy
Aperture: Gratings / Lens fibers
Eurographics 2012, Cagliari, Italy
Aperture: Vitreous Humor
Eurographics 2012, Cagliari, Italy
Aperture: Vitreous Humor
Eurographics 2012, Cagliari, Italy
Aperture: Eyelashes (optional)
Eurographics 2012, Cagliari, Italy
Chromatic Blur
• Compute one wavelength - Get others for free!
380 nm 770 nm
+ + + +
=
575 nm
Eurographics 2012, Cagliari, Italy
Convolution
=
= +
•
Convolution Billboard Bright pixels
HDR image PSF
Eurographics 2012, Cagliari, Italy
Convolution
Convolution Billboard
Eurographics 2012, Cagliari, Italy
Temporal Glare Pipeline
Eurographics 2012, Cagliari, Italy
Eurographics 2012, Cagliari, Italy
Psychophysical Experiment
• Goal:
Measuring the brightness boosts caused by glare illusion
• 2 methods, 6 patterns for each
• Gaussian: blurring kernel Cheap approximation
• Spencer et al.: human eye's PSF (disability glare) Optical correctness
• 10 subjects20 minutes per person
Yoshida et al. APGV2008
Eurographics 2012, Cagliari, Italy
Stimuli
B C D E F
A
a b c d e f
220 cd/m
2150 cd/m
2Yoshida et al. APGV2008 Method 1: Gaussian
Method 2: Spencer et al.
Eurographics 2012, Cagliari, Italy
Perceptual Experiment
Task:
Adjust the target disk luminance as close as possible to that of the Reference, but slightly yet visibly darker/brighter.
Yoshida et al. APGV2008
Eurographics 2012, Cagliari, Italy
Method I (Gaussian)
A
B
C
D
E
F
Yoshida et al. APGV2008
Eurographics 2012, Cagliari, Italy
Method II (Spencer et al.)
a
b c
d
e
f
Yoshida et al. APGV2008
Eurographics 2012, Cagliari, Italy
Trade-offs
Method I (Gaussian)
A
B
C
D
E
F a
b
c
d
e
f
Increasing perceived size!
Method II (Spencer et al.)
Causing Mach-bands!
Measuring brightness boost of glare illusion
–Increasing the perceived luminance by 20–35%–Gaussian blurring is equally effective
Trade-offs Gaussian/human eye's PSF
Yoshida et al. APGV2008
Eurographics 2012, Cagliari, Italy
Summary/Limitations
• Glare illusion might boost apparent brightness up to 30%
• Comprehensible model of light scattering in the eye taking into account dynamic eye elements
• Real-time rendering
• Other temporal low-level eye physics like
• Floaters
• Local adaptation (“After images”)
http://www.mpi-inf.mpg.de/resources/hdr/TemporalGlare/
Eurographics 2012, Cagliari, Italy
Acknowledgements
• I would like to thank Karol Myszkowski,
Grzegorz Krawczyk, Kaleigh Smith, Akiko
Yoshida, and Matthias Ihrke for help in
preparing slides.
Eurographics 2012, Cagliari, Italy
Retargeting Color Content: Color Issues in Tone Mapping
Alessandro Artusi Ramon Cajal Fellow
Eurographics 2012, Cagliari, Italy
Introduction to Color
Eurographics 2012, Cagliari, Italy
What is Color?
Stimulus Object
Source Light Human Visual System
Eurographics 2012, Cagliari, Italy
Quantifying Color
x d
I X ( ) ( ) ( )
0
y d
I Y ( ) ( ) ( )
0
z d
I Z ( ) ( ) ( )
0
I
z y x , ,
SPD of the light Reflectance of the object CIE color matching functions
Eurographics 2012, Cagliari, Italy
How Color is Produced?
Additive
(a) (b)
Subtractive
Eurographics 2012, Cagliari, Italy
Color Space
• Device dependent: the description of color information is related to the characteristics of a particular device
– Set of primaries – Technology
• Device independent: the description of color information is not dependent from the characteristics of a particular device
– CIEXYZ, CIELab, CIELuv etc...
Eurographics 2012, Cagliari, Italy
Chromaticity Diagram and MacAdam’s Ellipses
• MacAdam’s Ellipses
• contains all colors which are indistinguishable to an human observer from the color at the center of the ellipse
• the contour of the ellipse represents the just noticeable differences of chromaticity
Z Y X x X
Z Y X y Y
Eurographics 2012, Cagliari, Italy
Color Attributes by the CIE
• Hue
• Saturation
• Lightness
Perception
• Saturation is the colorfulness of an area judged in proportion to its brightness.
• Lightness Human vision has a nonlinear perceptual response to
luminance: The perceptual response to luminance is called lightness.
16
116
31
*
Y
nL Y
Y
n Y 008856 . 0
• Hue The degree to which a stimulus can be described as similar to
or different from stimuli that are described as red, green, blue, and yellow.
Eurographics 2012, Cagliari, Italy
Color in High Dynamic Range
• Color Ratio (Schlick 1994)
out in
in
out L
L RGB RGB
Mantiuk et al.. “Color Correction for Tone Mapping”, Proceedings Eurographics 2009.
RGB
outRGB
inL
outL
inColor Input
Color Output
Luminance Input
Luminance Output
Eurographics 2012, Cagliari, Italy
Color in High Dynamic Range
• Saturation Control (Thumblin and Turk 1999)
out s
in in
out L
L RGB RGB
s Saturation Parameter
Mantiuk et al.. “Color Correction for Tone Mapping”, Proceedings Eurographics 2009.
c Contrast Compression Under-saturated colors for S=C.
Eurographics 2012, Cagliari, Italy
Color in High Dynamic Range
Mantiuk et al.. “Color Correction for Tone Mapping”, Proceedings Eurographics 2009.
Eurographics 2012, Cagliari, Italy
Color Rendering Pipeline (8 Bit)
Image Acquisition
Device Independent
Displaying Non-Linearity
Device Dependent
Colorimetric Characterization
L
* Gamut’s MismatchC
Eurographics 2012, Cagliari, Italy
Colorimetric Characterisation of a Device
3x3 Matrix
X Y Z R G B
Spectrophotometer Device
3x3 Primaries Matrix X Y Z
R G B
Non-Linearity
linB linG linR blue green red
blue green red
blue green red Z Y X
Z Z Z
Y Y Y
X X X
1
linB deviceB blueTRC
deviceG linG greenTRC
deviceR linR redTRC
1 1 1 Minimisation Process RAW Data are Linear
Eurographics 2012, Cagliari, Italy
Gamma – Curve
...
Gamma Response for RED Gamma Response for BLUE
d b d b
R
i ( 1 )
i
A. Neumann, A. Artusi, L. Neumann, G. Zotti and W. Purgathofer “Accurate Display Gamma Function based on Human Observation“.
Eurographics 2012, Cagliari, Italy
Color Rendering Pipeline in HDR
HDR Image Acquisition
Device Independent
Displaying Non-Linearity
Device Dependent
Colorimetric Characterization Tone Mapping
out s
in in
out L
L RGB RGB
+
Eurographics 2012, Cagliari, Italy
HDR ICC Profile T
1T
n
ICC Profile
X Y Z R G B
Target
Best Exposure Image
Goesele et al. “Color Calibrated High Dynamic Range Imaging with ICC Profiles.”
Eurographics 2012, Cagliari, Italy
HDR Colorimetric Camera Characterization
Min H. Kim et al. “Characterization of High Dynamic Range Imaging.”
Radiance Measurements
HDR Radiance Map
3x3 Matrix
X Y Z
R G B
Eurographics 2012, Cagliari, Italy
Color Gamut
• Device
– Set of colors reproducible by the device
• Image
– Set of colors that compose the image
Eurographics 2012, Cagliari, Italy
Gamut Mapping
a b
Eurographics 2012, Cagliari, Italy
Gamut vs. Tone Mapping
L
wL
d0 . 1
10
5cd / m
210
2Not HDR Content
L
*a 10
510
2Eurographics 2012, Cagliari, Italy
Gamut Mapping Aims (CS)
• Gray axes alignment, mapping white to white and black to black
L
*a L
*0 a
0
Eurographics 2012, Cagliari, Italy
Gamut Mapping Aims (CS)
L
*0 a
• Gray axes alignment, mapping white to white and black to black
L
*0 a
Eurographics 2012, Cagliari, Italy
Gamut Mapping Aims (CS)
• Unchanged the Hue shift, will keep the overall image appearance
L
*C H L
*C Black
White
Eurographics 2012, Cagliari, Italy
Gamut Mapping Aims (CS)
• Limiting out of gamut colours
– Soft clipping can be afterwards adopted to eliminate these extremes
• Increase Image saturation
– Destination gamut has reduced saturation
– Helps maintaining the original chroma differences of the
input Image
Eurographics 2012, Cagliari, Italy
Gamut Mapping Pipeline
Gamut Boundary Descriptor
Line Gamut Boundary
Mapping Color
Space
Hue slice at 1 degree Clipping Compression L
*L
avgC Spatial GMAs
Eurographics 2012, Cagliari, Italy
Color Space Issue
• Gamut Mapping that preserves metric hue angle – No Hue shift after compression or clipping
• CIELab is suffering of non linearity in blue regions, but also in red regions
Braun and Fairchild. “Color Gamut Mapping in Hue-Linearized CIELab Color Space”
Eurographics 2012, Cagliari, Italy
Point-wise Gamut Mapping Techniques
• Clipping
It changes colours which are outside of the destination gamut,
mapping them on the boundaries of the destination gamut – Horizontal (lines of constant lightness)
– Radial to a centre of Gravity
• Centre of lightness axis (Constant)
• Lightness corresponding to the Chroma Cusp (variable) – Distance in CIELab
• To the colour boundary of the destination gamut that
has the smallest distance (HPMin Clipping) E
Eurographics 2012, Cagliari, Italy
Clipping
L
*C L
*C
2
*/ L
Ccusp
L* CcuspCmax
L
*C
L* Preservation Radial to L*/2
Cusp Radial Simultaneous GMAs
Eurographics 2012, Cagliari, Italy
Clipping – Major Drawbacks Erase Local Image variation (Details) L
*C
L
*C Preserve Saturation
Eurographics 2012, Cagliari, Italy
Eurographics 2012, Cagliari, Italy
Eurographics 2012, Cagliari, Italy
• Compression
It makes changes to all the colors of the source gamut to be
accommodated into the destination gamut . – Linear
– Sigmoid – Knee-function
• Parametric
The behaviour change based on the shapes of the two gamut’s
(source and destination) at the hue angle, or it depends from user
parameters. (Clipping and Compression)
Point-wise Gamut Mapping Techniques
Eurographics 2012, Cagliari, Italy
Compression
C
oC
i maxmax o
i
C
C
max
C
oimax
C Linear Knee-Function
Soft
Hard Clipping
Eurographics 2012, Cagliari, Italy
Compression
C
oC
i maxC
omax
C
i
2 2 2
2 1
xbab e
a b
Eurographics 2012, Cagliari, Italy
Parametric
L
*C L
*C L
*C 2 L
*Eurographics 2012, Cagliari, Italy
Preservation of Spatial Details
• Optimization
Making use of Human Visual System Models minimize the perceived differences between the input and output image.
• Multiscale
Re-inserts high-frequency information content in the gamut mapped image (clipped).
• Clipping – loss of details
• General framework has been proposed that includes the different
cases
Eurographics 2012, Cagliari, Italy
Preservation of Spatial Details
L*
C
f
Bonnier et al. “Spatial and Color Adaptive Gamut Mapping: A Mathematical Framework and Two New Algorithms.”
L*
MutliScale
CMerge Manipulation
k
Eurographics 2012, Cagliari, Italy
Automatic estimation of desaturation (s) factor in function of contrast compression (c) (non-linear color correction).
s = f(c) determined based on results of perceptual experiment
Mantiuk et al. “Color Correction for Tone Mapping
Eurographics 2012, Cagliari, Italy
Unchanged luminance value after color correction (luminance preserving solution)