Color image processing Presentation

(Images from Rafael C. Gonzalez and Richard E.
Wood, Digital Image Processing, 2nd Edition.
Spectrum of White Light
1666 Sir Isaac Newton, 24 year old, discovered white light spectrum.

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

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)

Primary and Secondary Colors
Primary
color
Primary
color
Primary
color
Secondary
colors

Primary and Secondary Colors (cont.)
Additive primary colors: RGB
use in the case of light sources
such as color monitors
Subtractive primary colors: CMY
use in the case of pigments in
printing devices
RGB add together to get white
White subtracted by CMY to get
Black

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
Color Characterization
Hue
Saturation
Chromaticity
amount of red (X), green (Y) and blue (Z) to form any particular
color is called tristimulus.

CIE Chromaticity Diagram
Trichromatic coefficients:
ZYX
X
x


ZYX
Y
y


ZYX
Z
z


1 zyx
x
y
Points on the boundary are
fully saturated colors

Color Gamut of Color Monitors and Printing Devices
Color Monitors
Printing devices

RGB Color Model
Purpose of color models: to facilitate the specification of colors in
some standard
RGB color models:
- based on cartesian
coordinate system

RGB Color Cube
R = 8 bits
G = 8 bits
B = 8 bits
Color depth 24 bits
= 16777216 colors
Hidden faces
of the cube

RGB Color Model (cont.)
Red fixed at 127

Safe RGB Colors
Safe RGB colors: a subset of RGB colors.
There are 216 colors common in most operating systems.

RGB 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.

CMY and CMYK Color Models
































B
G
R
Y
M
C
1
1
1
C = Cyan
M = Magenta
Y = Yellow
K = Black

HSI Color Model
RGB, CMY models are not good for human interpreting
HSI Color model:
Hue: Dominant color
Saturation: Relative purity (inversely proportional
to amount of white light added)
Intensity: Brightness
Color carrying
information

Relationship Between RGB and HSI Color Models
RGB HSI

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
Hue Saturation Intensity
RGB Cube

Converting Colors from RGB to HSI






GB
GB
H
if360
if


 
  











 
2/12
1
))(()(
)()(
2
1
cos
BGBRGR
BRGR

BGR
S


3
1
)(
3
1
BGRI 

Converting Colors from HSI to RGB
)1( SIB 








)60cos(
cos
1
H
HS
IR 
)(1 BRG 
RG sector: 1200  H GB sector: 240120  H
)1( SIR 








)60cos(
cos
1
H
HS
IG 
)(1 GRB 
)1( SIG 








)60cos(
cos
1
H
HS
IB 
)(1 BGR 
BR sector: 360240  H
120 HH
240 HH

Example: HSI Components of RGB Colors
Hue
Saturation Intensity
RGB
Image

Example: Manipulating HSI Components
Hue
RGB
Image Hue Saturation
Intensity RGB
Image

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
Why we need to assign colors to gray scale image?
Answer: Human can distinguish different colors better than different
shades of gray.
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.

Intensity Slicing or Density Slicing






TyxfC
TyxfC
yxg
),(if
),(if
),(
2
1
Formula:
C1 = Color No. 1
C2 = Color No. 2
T
Intensity
Color
C1
C2
T0 L-1
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.

Multi Level Intensity Slicing
kkk lyxflCyxg   ),(for),( 1
Ck = Color No. k
lk = Threshold level k
Intensity
Color
C
1
C2
0
L-1
l1 l2 l3 lklk-1
C3
Ck-1
Ck

Multi Level Intensity Slicing Example
kkk lyxflCyxg   ),(for),( 1
Ck = Color No. k
lk = Threshold level k
An X-ray image of the Picker
Thyroid Phantom.
After density slicing into 8 colors

Color Coding Example
A unique color is assigned to
each intensity value.
Gray-scale image of average
monthly rainfall.
Color coded image
Color
map
South America region

Gray Level to Color Transformation
Assigning colors to gray levels based on specific mapping functions
Red component
Green component
Blue component
Gray scale image

(Images from Rafael C.
Gonzalez and Richard
E. Wood, Digital Image
Processing, 2nd Edition.
Gray Level to Color Transformation Example
An X-ray image of a
garment bag with a
simulated explosive
device
An X-ray image
of a garment bag
Color
coded
images
Transformations

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

Pseudocolor Coding Example
Washington D.C. area
Visible blue
= 0.45-0.52 mm
Max water penetration
Visible green
= 0.52-0.60 mm
Measuring plant
Visible red
= 0.63-0.69 mm
Plant discrimination
Near infrared
= 0.76-0.90 mm
Biomass and shoreline mapping
1 2
3 4 Red =
Green =
Blue =
Color composite images
1
2
3
Red =
Green =
Blue =
1
2
4

Pseudocolor Coding Example
Psuedocolor rendition
of Jupiter moon Io
A close-up
Yellow areas = older sulfur deposits.
Red areas = material ejected from
active volcanoes.

Basics of 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.

Example: Full-Color Image and Variouis Color Space Components
Color image
CMYK components
RGB components
HSI components

Color Transformation
Formulation:
 ),(),( yxfTyxg 
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:
),,,( 21 nii rrrTs  i= 1, 2, …, n
Where ri = color component of f(x,y)
si = color component of g(x,y)
Use to transform colors to colors.
For RGB images, n = 3

Example: Color Transformation
),(),(
),(),(
),(),(
yxkryxs
yxkryxs
yxkryxs
BB
GG
RR



Formula for RGB:
),(),( yxkryxs II 
Formula for CMY:
)1(),(),(
)1(),(),(
)1(),(),(
kyxkryxs
kyxkryxs
kyxkryxs
YY
MM
CC



Formula for HSI:
These 3 transformations give
the same results.
k = 0.7
I H,S

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

Color Complement Transformation Example

Color Slicing Transformation









 
otherwise
2
if5.0
1
i
njany
jj
i
r
W
ar
s
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.
i= 1, 2, …, n
or
 





 
otherwise
if5.0
1
2
0
2
i
n
j
jj
i
r
Rar
s
Set to gray
Keep the original
color
Set to gray
Keep the original
color
i= 1, 2, …, n

Color Slicing Transformation Example
Original image
After color slicing

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.
Power law transformations
Contrast enhancement

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.

Histogram Equalization of a Full-Color Image
 Histogram equalization of a color image can be performed by
adjusting color intensity uniformly while leaving color unchanged.
 The HSI model is suitable for histogram equalization where only
Intensity (I) component is equalized.






k
j
j
k
j
jrkk
N
n
rprTs
0
0
)()(
where r and s are intensity components of input and output color image.

Gonzalez and Richard E.
Wood, Digital Image
Histogram Equalization of a Full-Color Image
Original image
After histogram
equalization
After increasing
saturation component

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
2. Smooth only Intensity component of a HSI image while leaving
H and S unmodified.



























xy
xy
xy
xy
Syx
Syx
Syx
Syx
yxB
K
yxG
K
yxR
K
yx
K
yx
),(
),(
),(
),(
),(
1
),(
1
),(
1
),(
1
),( cc
Note: 2 methods are not equivalent.

Color Image Smoothing Example (cont.)
Color image Red
Green Blue

Color image
HSI Components

Smooth all RGB components Smooth only I component of HSI
(faster)

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

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.

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.

Color Segmentation in HSI Color Space
Hue
Color image
1 2
3 4
Wood, Digital Image

Color Segmentation in HSI Color Space (cont.)
Product of and
5 6
7 8
52Binary thresholding of S component
with T = 10%
Histogram of 6 Segmentation of red color pixels
Red pixels
Wood, Digital Image

Color Segmentation in HSI Color Space (cont.)
Color image Segmented results of red pixels
Wood, Digital Image

Color Segmentation in RGB Vector Space
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






TyxD
TyxD
yxg
T
T
)),,((if0
)),,((if1
),(
cc
cc
cT = color to be segmented.
c(x,y) = RGB vector at pixel (x,y).D(u,v) = distance function

Example: Segmentation in RGB Vector Space
Color image
Results of segmentation in
RGB vector space with Threshold
value
Reference color cT to be segmented
boxtheinpixelofcoloraverageTc
T = 1.25 times the SD of R,G,B values
In the box

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
Edges
We see
4 objects.
We see
2 objects.

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:
  2
1
2sin22cos)()(
2
1
)(






  xyyyxxyyxx gggggF
 








 
yyxx
xy
gg
g2
tan
2
1 1

222
x
B
x
G
x
R
gxx









222
y
B
y
G
y
R
gyy









y
B
x
B
y
G
x
G
y
R
x
R
gxy
















Obtained using
the formula
in the previous
slide
Sum of
gradients of
each color
component
Original
image
Difference
between
2 and 3
2
3
2 3
Gradient of a Color Image Example

Gradients of each color component
Red Green Blue
Gradient of a Color Image Example

Noise in Color Images
Noise can corrupt each color component independently.
Wood, Digital Image
Noise is less
noticeable
in a color
image
AWGN sh
2=800 AWGN sh
2=800
AWGN sh
2=800

Hue
Salt & pepper noise
in Green component

Color Image Compression
JPEG2000 File
Original image
After lossy compression with ratio 230:1

Color image processing Presentation

More Related Content

What's hot (20)

Viewers also liked (10)

Similar to Color image processing Presentation (20)

Recently uploaded (20)

Color image processing Presentation