Digital Image Processing

Chapter 6:
Color Image Processing
Spectrum of White Light

1666 Sir Isaac Newton, 24 year old, discovered white light spectrum.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Electromagnetic Spectrum

Visible light wavelength: from around 400 to 700 nm

1. For an achromatic (monochrome) light source,

there is only 1 attribute to describe the quality: intensity

2. For a chromatic light source, there are 3 attributes to describe

the quality:
Radiance = total amount of energy flow from a light source (Watts)
Luminance = amount of energy received by an observer (lumens)
Brightness = intensity
(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Cross section illustration

The Eye
UMCP ENEE631 Slides (created by M.Wu © 2004)

Figure is from slides at

Gonzalez/ Woods DIP book
website (Chapter 2)
Two Types of Photoreceptors at Retina
• Rods
– Long and thin
– Large quantity (~ 100 million)
– Provide scotopic vision (i.e., dim light vision or at low illumination)
– Only extract luminance information and provide a general overall picture
• Cones
– Short and thick, densely packed in fovea (center of retina)
– Much fewer (~ 6.5 million) and less sensitive to light than rods
– Provide photopic vision (i.e., bright light vision or at high illumination)
– Help resolve fine details as each cone is connected to its own nerve end
– Responsible for color vision
our interest
– Mesopic vision (well-lighted display)
• provided at intermediate illumination by both rod and cones
Sensitivity of Cones in the Human Eye

6-7 millions cones

in a human eye
- 65% sensitive to Red light
- 33% sensitive to Green light
- 2 % sensitive to Blue light

Primary colors:
Defined CIE in 1931
Red = 700 nm
Green = 546.1nm
Blue = 435.8 nm

CIE = Commission Internationale de l’Eclairage

(The International Commission on Illumination) (Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Luminance vs. Brightness

Same lum. Different lum.

Different Similar
brightness brightness

• Luminance (or intensity)

– Independent of the luminance of surroundings

I(x,y,λ) -- spatial light distribution

V(λ) -- relative luminous efficiency func. of visual system ~ bell shape
(different for scotopic vs. photopic vision;
highest for green wavelength, second for red, and least for blue )
• Brightness
– Perceived luminance
– Depends on surrounding luminance
Luminance vs. Brightness (cont’d)

• Example: visible digital watermark

– How to make the watermark
appears the same graylevel
all over the image?

from IBM Watson web page

“Vatican Digital Library”
Look into Simultaneous Contrast Phenomenon

• Human perception more sensitive to luminance

contrast than absolute luminance
• Weber’s Law: | Ls – L0 | / L0 = const
– Luminance of an object (L0) is set to be just noticeable
from luminance of surround (Ls)
– For just-noticeable luminance difference ∆L:
∆L / L ≈ d( log L ) ≈ 0.02 (const)
• equal increments in log luminance are perceived as equally different

• Empirical luminance-to-contrast models

– Assume L ∈ [1, 100], and c ∈ [0, 100]
– c = 50 log10 L (logarithmic law, widely used)
– c = 21.9 L1/3 (cubic root law)
 Mach Bands
UMCP ENEE631 Slides (created by M.Wu © 2004)

Figure is from slides

at Gonzalez/ Woods
DIP book website
(Chapter 2)

• Visual system tends to undershoot or overshoot around the

boundary of regions of different intensities
è Demonstrates the perceived brightness is not a simple
function of light intensity
Color of Light

• Perceived color depends on spectral content (wavelength

UMCP ENEE408G Slides (created by M.Wu & R.Liu © 2002)

composition)
– e.g., 700nm ~ red.
– “spectral color”
• A light with very narrow bandwidth

“Spectrum” from http://www.physics.sfasu.edu/astro/color.html

• A light with equal energy in all visible bands appears

white
Primary and Secondary Colors

Primary
color
Secondary
colors

Primary Primary
color color
Primary and Secondary Colors (cont.)

Additive primary colors: RGB

use in the case of light sources
such as color monitors
RGB add together to get white

Subtractive primary colors: CMY

use in the case of pigments in
printing devices
White subtracted by CMY to get
Black

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Representation by Three Primary Colors

• Any color can be reproduced by mixing an appropriate set of

three primary colors (Thomas Young, 1802)
• Three types of cones in human retina
– Absorption response Si(λ) has peaks around 450nm (blue), 550nm
(green), 620nm (yellow-green)
– Color sensation depends on the spectral response {α1(C), α2(C),
α3(C) } rather than the complete light spectrum C(λ)

∫ S1(λ) C(λ) d λ α1(C)

C(λ)
∫ S2(λ) C(λ) d λ α2(C) Identically
color light perceived colors
α3(C) if αi (C1) = αi (C2)
∫ S3(λ) C(λ) d λ
Example: Seeing Yellow Without Yellow

570nm
520nm 630nm

mix green and red light to obtain perception of

yellow, without shining a single yellow photon

“Seeing Yellow” figure is from B.Liu ELE330 S’01 lecture notes @ Princeton;
R/G/B cone response is from slides at Gonzalez/ Woods DIP book website
Color Matching and Reproduction

• Mixture of three primaries: C = Sum(βk Pk (λ) )

• To match a given color C1
– adjust βk such that αi (C1) = αi (C), i = 1,2,3.

• Tristimulus values Tk (C)

– Tk (C) = βk / wk
wk – the amount of kth primary to match the reference white
• Chromaticity tk = Tk / (T1+T2+T3)
– t1+t2+t3 = 1
– visualize (t1, t2 ) to obtain chromaticity diagram
Color Characterization

Hue: dominant color corresponding to a dominant

wavelength of mixture light wave
Saturation: Relative purity or amount of white light mixed
with a hue (inversely proportional to amount of white
light added)
Brightness: Intensity

Hue
Chromaticity
Saturation

amount of red (X), green (Y) and blue (Z) to form any particular
color is called tristimulus.
 Perceptual Attributes of Color
• Value of Brightness
(perceived luminance)
UMCP ENEE408G Slides (created by M.Wu & R.Liu © 2002)

• Chrominance
– Hue
• specify color tone (redness, greenness, etc.)
• depend on peak wavelength

– Saturation
• describe how pure the color is
• depend on the spread (bandwidth) of light
spectrum
• reflect how much white light is added
HSV circular cone is from online

• RGB ó HSV Conversion ~ nonlinear documentation of Matlab image

processing toolbox
http://www.mathworks.com/access
/helpdesk/help/toolbox/images/col
or10.shtml
CIE Chromaticity Diagram
Trichromatic coefficients:
X
x=
X +Y + Z
Y
y=
X +Y + Z
y Z
z=
X +Y + Z

x + y + z =1

Points on the boundary are

fully saturated colors

x (Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Gamut of Color Monitors and Printing Devices

Color Monitors

Printing devices

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
CIE Color Coordinates (cont’d)
• CIE XYZ system

– hypothetical primary sources to yield all-positive spectral

tristimulus values
– Y ~ luminance
• Color gamut of 3 primaries
– Colors on line C1 and C2 can be
produced by linear mixture of the two
– Colors inside the triangle gamut
can be reproduced by three primaries

From http://www.cs.rit.edu/~ncs/color/t_chroma.html
RGB Color Model
Purpose of color models: to facilitate the specification of colors in
some standard

RGB color models:

- based on cartesian
coordinate system

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
RGB Color Cube

R = 8 bits Color depth 24 bits

G = 8 bits = 16777216 colors
B = 8 bits

Hidden faces
of the cube

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
RGB Color Model (cont.)

Red fixed at 127

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Safe RGB Colors
Safe RGB colors: a subset of RGB colors.
There are 216 colors common in most operating systems.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
RGB Safe-
Safe-color Cube

The RGB Cube is divided into

6 intervals on each axis to achieve
the total 63 = 216 common colors.
However, for 8 bit color
representation, there are the total
256 colors. Therefore, the remaining
40 colors are left to OS.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
CMY and CMYK Color Models
C = Cyan
M = Magenta
• Primary colors for pigment Y = Yellow
– Defined as one that subtracts/absorbs a K = Black
primary color of light & reflects the
other two
• CMY – Cyan, Magenta, Yellow
– Complementary to RGB
– Proper mix of them produces black

 C  1  R 
 M  = 1 − G 
    
 Y  1  B 
HSI Color Model
RGB, CMY models are not good for human interpreting

HSI Color model:

Hue: Dominant color
Color carrying
Saturation: Relative purity (inversely proportional information
to amount of white light added)

Intensity: Brightness

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Relationship Between RGB and HSI Color Models

RGB HSI
(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Hue and Saturation on Color Planes

1. A dot is the plane is an arbitrary color

2. Hue is an angle from a red axis.
3. Saturation is a distance to the point.
HSI Color Model (cont.)

Intensity is given by a position on the vertical axis.

HSI Color Model

Intensity is given by a position on the vertical axis.

Example: HSI Components of RGB Cube

RGB Cube

Hue Saturation Intensity

(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Converting Colors from RGB to HSI

θ if B ≤ G
H =
360 − θ if B > G
 

1
[( R − G ) + ( R − B)] 
θ = cos −1  2
1/ 2 
[ ]
 ( R − G ) + ( R − B )(G − B ) 
2

 
3
S = 1−
R+G+ B
1
I = ( R + G + B)
3
Converting Colors from HSI to RGB
RG sector: 0 ≤ H < 120 GB sector:120 ≤ H < 240
H = H − 120
 S cos H 
R = I 1 +  R = I (1 − S )
 cos( 60 o
− H ) 
B = I (1 − S )  S cos H 
G = I 1 + 
 cos( 60 o
− H ) 
G = 1 − ( R + B)
B = 1 − ( R + G)
BR sector: 240 ≤ H ≤ 360
H = H − 240
 S cos H 
B = I 1 + 
 cos( 60 o
− H ) 
G = I (1 − S )

R = 1 − (G + B )
Example: HSI Components of RGB Colors

RGB
Image Hue

Saturation Intensity

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Example: Manipulating HSI Components
RGB
Image Hue Hue Saturation

Saturation Intensity Intensity RGB

Image
(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Color Coordinates Used in TV Transmission
• Facilitate sending color video via 6MHz mono TV
channel
• YIQ for NTSC (National Television Systems Committee)
transmission system
– Use receiver primary system (RN, GN, BN) as TV receivers
standard
– Transmission system use (Y, I, Q) color coordinate
• Y ~ luminance, I & Q ~ chrominance
• I & Q are transmitted in through orthogonal carriers at the same freq.

• YUV (YCbCr) for PAL and digital video

– Y ~ luminance, Cb and Cr ~ chrominance
Color Coordinates

• RGB of CIE
• XYZ of CIE
• RGB of NTSC
• YIQ of NTSC
• YUV (YCbCr)
• CMY
Examples

RGB

HSV

YUV
Examples

RGB

HSV

YIQ
 Summary
UMCP ENEE631 Slides (created by M.Wu © 2004)

• Monochrome human vision

– visual properties: luminance vs. brightness, etc.
– image fidelity criteria
• Color
– Color representations and three primary colors
– Color coordinates
Color Image Processing

There are 2 types of color image processes

1. Pseudocolor image process: Assigning colors to gray

values based on a specific criterion. Gray scale images to be processed
may be a single image or multiple images such as multispectral images

2. Full color image process: The process to manipulate real

color images such as color photographs.
Pseudocolor Image Processing

Pseudo color = false color : In some case there is no “color” concept

for a gray scale image but we can assign “false” colors to an image.
Why we need to assign colors to gray scale image?

Answer: Human can distinguish different colors better than different

shades of gray.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Intensity Slicing or Density Slicing
Formula:
 C1 if f ( x, y ) ≤ T C1 = Color No. 1
g ( x, y ) =  C2 = Color No. 2
C2 if f ( x, y ) > T

Color
T C2

C1

0 T L-1
Intensity
A gray scale image viewed as a 3D surface.
Intensity Slicing Example

An X-ray image of a weld with cracks

After assigning a yellow color to pixels with

value 255 and a blue color to all other pixels.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Multi Level Intensity Slicing
g ( x , y ) = Ck for lk −1 < f ( x, y ) ≤ lk

Ck = Color No. k
lk = Threshold level k
Color

Ck

Ck-1

C3

C2
C
1
0 l1 l2 l3 lk-1 lk
L-1
Intensity
Multi Level Intensity Slicing Example
Ck = Color No. k
g ( x , y ) = Ck for lk −1 < f ( x, y ) ≤ lk lk = Threshold level k

An X-ray image of the Picker After density slicing into 8 colors

Thyroid Phantom.
(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Color Coding Example
>20

Gray Color
Scale 10 map

Gray-scale image of average

monthly rainfall.
0

Color coded image South America region

(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Gray Level to Color Transformation

Assigning colors to gray levels based on specific mapping functions

Red component

Gray scale image

Green component

Blue component

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Gray Level to Color Transformation Example
An X-ray image An X-ray image of a
of a garment bag garment bag with a
simulated explosive
device

(Images from Rafael C.

Gonzalez and Richard
E. Wood, Digital Image
Processing, 2nd Edition. Transformations

Color
coded
images
Gray Level to Color Transformation Example
An X-ray image An X-ray image of a
of a garment bag garment bag with a
simulated explosive
device

Transformations (Images from Rafael C.

Gonzalez and Richard
E. Wood, Digital Image
Processing, 2nd Edition.

Color
coded
images
Pseudocolor Coding
Used in the case where there are many monochrome images such as multispectral
satellite images.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Pseudocolor Coding Example
Pseudocolor Coding Example
Visible blue Visible green
λ= 0.45-0.52 µm λ= 0.52-0.60 µm
Max water penetration Measuring plant Color composite images
1 2

3 4 Red = 1 Red = 1
Green = 2 Green = 2
Blue = 3 Blue = 4

Better visualization àShow quite

clearly the difference between
biomass (red) and human-made features.

Visible red Near infrared Washington D.C. area

λ= 0.63-0.69 µm λ= 0.76-0.90 µm
(Images from Rafael C. Gonzalez and Richard E.
Plant discrimination Biomass and shoreline mapping Wood, Digital Image Processing, 2nd Edition.
Pseudocolor Coding Example

Psuedocolor rendition
of Jupiter moon Io

Yellow areas = older sulfur deposits.

Red areas = material ejected from
active volcanoes.

A close-up

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Basics of Full-
Full-Color Image Processing
2 Methods:
1. Per-color-component processing: process each component separately.
2. Vector processing: treat each pixel as a vector to be processed.

Example of per-color-component processing: smoothing an image

By smoothing each RGB component separately.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Example: Full
Full--Color Image and Variouis Color Space Components

Color image

CMYK components

RGB components

HSI components

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Transformation
Use to transform colors to colors.

Formulation:
g ( x , y ) = T [ f ( x , y )]

f(x,y) = input color image, g(x,y) = output color image

T = operation on f over a spatial neighborhood of (x,y)

When only data at one pixel is used in the transformation, we

can express the transformation as:

si = Ti ( r1 , r2 ,K, rn ) i= 1, 2, …, n

Where ri = color component of f(x,y) For RGB images, n = 3

si = color component of g(x,y)
Example: Color Transformation
Formula for RGB:
sR ( x, y ) = krR ( x, y )
sG ( x, y ) = krG ( x, y )
sB ( x, y ) = krB ( x, y )

Formula for HSI: k = 0.7

sI ( x, y ) = krI ( x, y )

Formula for CMY: I H,S

sC ( x, y ) = krC ( x, y ) + (1 − k )
sM ( x, y ) = krM ( x, y ) + (1 − k ) These 3 transformations give
sY ( x, y ) = krY ( x, y ) + (1 − k ) the same results.
(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Color Complements
Color complement replaces each color with its opposite color in the
color circle of the Hue component. This operation is analogous to
image negative in a gray scale image.

Color circle

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Complement Transformation Example

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Slicing Transformation

We can perform “slicing” in color space: if the color of each pixel

is far from a desired color more than threshold distance, we set that
color to some specific color such as gray, otherwise we keep the
original color unchanged.

  W Set to gray
0.5 if  rj − a j > 
si =   2  any 1≤ j ≤n
 ri otherwise Keep the original
i= 1, 2, …, n color
or


( )
n
0.5 if ∑ rj − a j > R0
2 2 Set to gray
si =  j =1
 ri otherwise Keep the original
color
i= 1, 2, …, n
Color Slicing Transformation Example

After color slicing

Original image

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Tonal Correction Examples
In these examples, only
brightness and contrast are
adjusted while keeping color
unchanged.
This can be done by
using the same transformation
for all RGB components.

Contrast enhancement

Power law transformations

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Balancing Correction Examples

Color imbalance: primary color components in white area

are not balance. We can measure these components by
using a color spectrometer.

Color balancing can be

performed by adjusting
color components separately
as seen in this slide.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Histogram Equalization of a Full-
Full-Color Image

v Histogram equalization of a color image can be performed by

adjusting color intensity uniformly while leaving color unchanged.

v The HSI model is suitable for histogram equalization where only

Intensity (I) component is equalized.
k
sk = T ( rk ) = ∑ pr ( rj )
j =0
k
nj
=∑
j =0 N

where r and s are intensity components of input and output color image.
Histogram Equalization of a Full-
Full-Color Image
Original image

After histogram After increasing

equalization saturation component

(Images from Rafael C.

Gonzalez and Richard E.
Wood, Digital Image
Processing, 2nd Edition.
Color Image Smoothing

2 Methods:
1. Per-color-plane method: for RGB, CMY color models
Smooth each color plane using moving averaging and
the combine back to RGB

1 
 ∑ R ( x, y ) 
 K ( x , y )∈S xy 
1 1 
c ( x, y ) = ∑ c( x , y ) = 
K ( x , y )∈S xy K
∑ G ( x, y )
 ( x , y )∈S xy 
1 
 K ( x ,∑ B ( x , y ) 
 y )∈S xy 

2. Smooth only Intensity component of a HSI image while leaving

H and S unmodified.
Note: 2 methods are not equivalent.
Color Image Smoothing Example (cont.)

Color image Red

Green Blue

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Image Smoothing Example (cont.)

Color image

HSI Components

Hue Saturation Intensity

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Image Smoothing Example (cont.)

Smooth all RGB components Smooth only I component of HSI

(faster)

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Image Smoothing Example (cont.)

Difference between
smoothed results from 2
methods in the previous
slide.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Image Sharpening
We can do in the same manner as color image smoothing:
1. Per-color-plane method for RGB,CMY images
2. Sharpening only I component of a HSI image

Sharpening all RGB components Sharpening only I component of HSI

Color Image Sharpening Example (cont.)

Difference between
sharpened results from 2
methods in the previous
slide.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Segmentation

2 Methods:
1. Segmented in HSI color space:
A thresholding function based on color information in H and S
Components. We rarely use I component for color image
segmentation.

2. Segmentation in RGB vector space:

A thresholding function based on distance in a color vector space.

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Segmentation in HSI Color Space
Color image Hue

1 2

3 4

Saturation Intensity (Images from Rafael C.

Gonzalez and Richard E.
Wood, Digital Image
Processing, 2nd Edition.
Color Segmentation in HSI Color Space (cont.)
Binary thresholding of S component Product of 2 and 5
with T = 10%

5 6

Red pixels

7 8
(Images from Rafael C.
Gonzalez and Richard E.
Wood, Digital Image Histogram of 6 Segmentation of red color pixels
Processing, 2nd Edition.
Color Segmentation in HSI Color Space (cont.)

Color image Segmented results of red pixels

(Images from Rafael C.

Gonzalez and Richard E.
Wood, Digital Image
Processing, 2nd Edition.
Color Segmentation in RGB Vector Space

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.

1. Each point with (R,G,B) coordinate in the vector space represents

one color.
2. Segmentation is based on distance thresholding in a vector space

1 if D (c( x, y ), cT ) ≤ T
g ( x, y ) = 
0 if D (c( x, y ), cT ) > T
cT = color to be segmented.
D(u,v) = distance function c(x,y) = RGB vector at pixel (x,y).
Example: Segmentation in RGB Vector Space

Color image

Reference color cT to be segmented

cT = average color of pixel in the box

Results of segmentation in
RGB vector space with Threshold
value
T = 1.25 times the SD of R,G,B values
In the box

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Gradient of a Color Image
Since gradient is define only for a scalar image, there is no concept
of gradient for a color image. We can’t compute gradient of each
color component and combine the results to get the gradient of a color
image.
Red Green Blue We see
2 objects.

We see
4 objects.

Edges

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Gradient of a Color Image (cont.)
One way to compute the maximum rate of change of a color image
which is close to the meaning of gradient is to use the following
formula: Gradient computed in RGB color space:
1

F (θ ) =  [( g xx + g yy ) + ( g xx − g yy ) cos 2θ + 2 g xy sin 2θ ]
1  2

2 

1 −1
 2 g xy 
θ = tan  
2  (g xx − g yy )
2 2 2
∂R ∂G ∂B ∂R ∂G ∂B
2 2 2

g xx = + + g yy = + +
∂x ∂x ∂x ∂y ∂y ∂y

∂R ∂R ∂G ∂G ∂B ∂B
g xy = + +
∂x ∂y ∂x ∂y ∂x ∂y
Gradient of a Color Image Example
2

Original Obtained using

image the formula
in the previous
slide

Sum of
gradients of Difference
each color between
component 22 and 33

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Gradient of a Color Image Example

Red Green Blue

Gradients of each color component

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Noise in Color Images
Noise can corrupt each color component independently.

AWGN ση2=800 AWGN ση2=800

AWGN ση2=800
Noise is less
noticeable
in a color
image

(Images from Rafael C.

Gonzalez and Richard E.
Wood, Digital Image
Processing, 2nd Edition.
Noise in Color Images

Hue Saturation Intensity

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Noise in Color Images
Salt & pepper noise Hue
in Green component

Saturation Intensity

(Images from Rafael C. Gonzalez and Richard E.

Wood, Digital Image Processing, 2nd Edition.
Color Image Compression
Original image

JPEG2000 File

After lossy compression with ratio 230:1 (Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.