Light Scattering from Human Hair Fibers
Stephen R. Marschner
Cornell University

Henrik Wann Jensen
University of California—San Diego
Steve Worley
Worley Laboratories

Pat Hanrahan
Stanford University
R

Abstract
Light scattering from hair is normally simulated in computer graphics using Kajiya and Kay’s classic phenomenological model. We
have made new measurements of scattering from individual hairfibers that exhibit visually significant effects not predicted by Kajiya and Kay’s model. Our measurements go beyond previous hairmeasurements by examining out-of-plane scattering, and together
with this previous work they show a multiple specular highlight andvariation in scattering with rotation about the fiber axis. We explain
the sources of these effects using a model of a hair fiber as a transparent elliptical cylinder with an absorbing interior and a surface
covered with tilted scales. Based on an analytical scattering function for a circular cylinder, we propose a practical shading modelfor hair that qualitatively matches the scattering behavior shown inthe measurements. In a comparison between a photograph and rendered images, we demonstrate the new model’s ability to match theappearance of real hair.
CR Categories: I.3.7 [Computer Graphics]: Three-Dimensional
Graphics and Realism—Shading

Mike Cammarano
Stanford University

TRT
> 2α

2α

2α

interior:
refrac. index ηabsorption σa< 2αellipticalcross section
axis ratio a:1αtiltedcuticlescalesTTsurface roughness β
root

tip

Figure 1: A schematic of our model for a hair fiber. The dashedlines indicate the scattering angles for a cylinder without tilted surface scales.

Keywords: hair, fibers, optical scattering, rendering

1

Introduction

Realistically rendering hair is essential to portraying people and animals, but achieving realism means confronting many challenges
inherent in the structure and behavior of hair. The geometry of
thousands to millions of strands must be determined; the scatteringof light from the fibers must be simulated; and the resulting thincurves in the image must be sampled without introducing aliasing. Most of the research effort on hair has gone toward modeling and animating the geometry of a collection of fibers [Sourin
et al. 1996; Chen et al. 1999; Hadap and Magnenat-Thalmann
2000; Kim and Neumann 2000; Kim and Neumann 2002; Changet al. 2002; Magnenat-Thalmann et al. 2002], efficient antialiasing and curve-drawing algorithms for scanline rendering [Watanabe and Suenaga 1992; Kong and Nakajima 2000], and approximating shadows within the hair [Lokovic and Veach 2000; Kim and
Neumann 2001]. However, considerably less attention has been
put toward the scattering model that determines the appearance ofa particular assembly of fibers. Sometimes models for reflection
from grooved surfaces [Kajiya 1985; Poulin and Fournier 1990] areused to compute scattering from combed fibers, but most hair renderings have used the classic phenomenological model of Kajiya
and Kay [1989].
Kajiya and Kay’s model was designed to capture the most obvious feature of scattering from a fiber—namely the appearance ofa linear highlight in the image running perpendicular to the fiber
directions. It is based on the observation that the reflection of a parallel beam from the surface of a cylinder will be in a cone centeredon the hair axis. The model places a constant-intensity highlightcentered on that cone. All other scattering is accounted for by a diffuse term that produces radiance proportional to the cosine of theincident angle. Although it has served well for many years, thismodel falls short in several respects. For one thing, it is not energy
conserving, which is important for physically based rendering.
But even when physical correctness is not needed, the KajiyaKay model fails to predict some observed visual effects. Since it
models fibers as opaque cylinders, it does not account for transmission or internal reflection. Hair is a dielectric material, and blond,brown, red, or other light colored hair is very translucent. Goldman has simulated translucency by adding a directional parameterthat controls the relative amount of forward transmission and backward reflection [Goldman 1997]. Kim has extended this model by
proposing a two-term phase function [Kim 2002]. The first term
models surface reflection as a cardioid, based on a ray density argument; our R component is similar but accounts for the Fresnel factor
and handles oblique incidence. The second term models transmission as a forward scattering cosine lobe, supported by Monte Carlocomputations of normal-incidence scattering from a transparent circular cylinder. We extend Kim’s model by accounting for the Fresnel factor and volume absorption in an analytical model that also
includes internal reflection, modeling the separation of highlightsfrom different reflection modes, and approximating the effects ofeccentricity.

Two research groups working in the cosmetics industry havemeasured the scattering properties of hair fibers in order to understand what makes hair look pleasing and healthy. Stamm etal. [1977] measured relative scattering from a rack of individual
fibers as a function of illumination and viewing angles for directions coplanar with the fiber axis (that is, in the incidence plane).
They discovered two deviations from the expected behavior. First,
the primary specular peak occurs at an angle several degrees awayfrom the specular direction. They speculated that the angular shift
is caused by the tilt of the scales forming the cuticle of the hair fiber.
Second, there is a secondary lobe that occurs on the other side of
the specular direction. Bustard and Smith [1991] have observedthat the first peak preserves polarization, whereas the second peakis depolarized. Both groups also observed that the secondary highlight is not present in black hair. Bustard and Smith also report apreliminary observation of azimuthal scattering in the plane normal
to the fiber in which they observed strong peaks, which they conjectured are internal reflection caustics. These experiments support
the hypothesis that the secondary peak is due to internal reflection
off the back side of the fiber.
We have made an experimental and theoretical study of the scattering of light from individual fibers of human hair, and developed
a simple practical model for use in computer graphics. Specifically,• In addition to the incidence-plane measurements of Stamm etal. and Bustard and Smith, we measure scattering in the plane
normal to the fiber. We also report for the first time full 3D
hemispherical scattering measurements. Several new appearance phenomena are visible in the 3D data. First, the primary
specular highlight continues all the way around the hair, while
the secondary highlight is confined to the side of the hair toward the source. Second, a pair of large out-of-plane peaks,
or glints, are present, and as the incidence angle increases thepeaks move closer to the incidence plane, eventually merging
and disappearing.
• The scattering distribution, especially the secondary highlight, depends on the angle of rotation of the hair fiber about
its axis. This is to be expected because hair fibers are notgenerally circular in cross section [Robbins 1994].
• To explain our observations, we propose a simple model basedon the structure and composition of hair fibers (Figure 1). We
approximate the hair as a transparent circular cylinder witha colored interior and a surface composed of rough, tilted
scales. Using geometric optics, we derive predictions for the
far-field scattered light distribution from such a cylinder. Our
theory predicts that three transport modes are significant: surface reflection, transmission, and internal reflection.
• Based on the theory, we propose a shading model for humanhair. We illustrate its use by producing several example renderings of brown hair. The shading model is simple and easily

Figure 2: An electron micrograph of a hair fiber that shows the
structure of the outer cuticle surface, which is composed of thinoverlapping scales [Robbins 1994]. In this image the fiber is oriented with the root at the top and the tip at the bottom.

w
ωi
θi

u

φi

φr

v

θr

ωr

Figure 3: Notation for scattering geometry.
incorporated into standard rendering algorithms. Finally, we
validate the model against experimental data and show thereis good qualitative agreement.

2

FibersWe begin by presenting some background information on humanhair fibers. We also introduce the notation we use in this paper, anddefine the light scattering function for fibers.

2.1

Hair fibers and fiber scattering

A fiber of human scalp hair is composed of two main parts: the cuticle and the cortex. The cuticle is a thin protective outside sheaththat surrounds the inner cortex. The cuticle is of particular importance for light scattering, since it forms the interface between thefiber and the air. It is composed of flat cells that overlap like roof
shingles (Figure 2) and make the fiber appear as a nested set ofcones. Because of their overlapped arrangement, the surfaces ofthe scales deviate slightly but systematically from the overall normal of the fiber’s surface, tilting their surfaces toward the root endof the fiber by approximately 3◦ [Bustard and Smith 1991; Robbins
1994].
The cortex forms the bulk of the fiber. At the center is a pigmented core, the medulla. The pigments in the cortex and medulla
determine the hair color. In this paper, we approximate the optical properties of the interior of the fiber with two quantities: anindex of refraction η of approximately 1.55 [Stamm et al. 1977]
and an absorption cross section σa , which is assumed to be uniform
throughout the fiber.
As we will see, the morphology of the hair fiber provides an explanation for the main features of the scattering function. The fiber
is modeled as a dielectric cylinder covered with tilted scales andwith a pigmented interior (Figure 1). Reflection of a directional
beam from a cylinder produces several distinct components of reflected light, all directed into a cone of outgoing directions. Reflection from the surface produces a component that is spread fairly
uniformly around the cone. Transmission through the hair produces
a very bright component that is focused toward directions on the opposite side of the hair. Light that reflects off the inside of the hairsurface produces a back-scattering component. This component ismore complicated and, depending on the shape of the cylinder andthe angle of incidence, may spread the energy over the whole coneor focus it into one or more caustics. Using T and R to stand fortransmission and reflection across a cylinder interface, we denotethese three modes of reflection as R, TT, and TRT.
As illustrated in Figure 1, the tilt of the cuticle scales shifts thecones of the R and TRT components slightly off the ideal specularcone. More importantly, they shift in opposite directions, causing
them to separate into two visually distinguishable highlights. The

10Black hair

Blond hair

510Synthetic

5

10

5

θr

θr

θr

1
.5

1
.5

1
.5

.1

.1

.1

10

10

5

5

10
–30º

30º
5

θr

–60º

θr

1
.5
.1
ROOT

θr

1
.5

1
.5

.1
TIP

.1

ROOT

TIP

ROOT

TIP

Figure 4: Measurements of scattering in the incidence plane: scattering as a function of scattering angle with illumination at 45◦ from the tip
and root ends. Black and blond hair are shown, along with a synthetic fiber from a wig.

R highlight, since it is surface reflection, is white, whereas the TRT
highlight, which is formed by light that passes through the interiorof the fiber, is colored.

2.2

Scattering

Before getting into more detail, let us establish the notation we will
use throughout the paper to describe the scattering geometry (Figure 3). The tangent to the hair is u, pointing in the direction fromthe root toward the tip; the vectors v and w complete a right-handedorthonormal basis, and if the cross section is elliptical v is the major axis and w is the minor axis. We refer to the v–w plane as thenormal plane. The direction of illumination is ωi , and the directionin which scattered light is being computed or measured is ωr ; both
directions point away from the center. We express ωi and ωr in
spherical coordinates. The inclinations with respect to the normal
plane are denoted θi and θr (measured so that 0◦ is perpendicular tothe hair, 90◦ is u, and −90◦ is −u). The azimuths around the hairare denoted φi and φr (measured so that v is 0◦ and w is +90◦).
We also use several derived angles. The difference angle (θr −θi )/2 is denoted θd . The relative azimuth φr − φi is denoted simply
φ . The averages θh = (θi + θr )/2 and φh = (φi + φr )/2 are called
half angles.
The bidirectional scattering function S for a fiber is different
from the bidirectional reflection distribution function fr for a surface, although it shares the same physical units. For the incident
and reflected light we use curve irradiance Ē, or power per unit
length, and curve intensity L̄, or intensity per unit length, respectively. These units are analogous to irradiance (power per unit area)
and radiance (intensity per unit area) on a surface.
S(ωi , ωr ) =d L̄r (ωr ),d Ēi (ωi )where L̄r is the curve intensity scattered from an infinitesimal lengthof fiber, and Ēi is the curve irradiance on that portion of the fiber.
This irradiance is proportional to incoming radiance:d Ēi (ωi ) = DLi (ωi ) cos θi dωi ,where D is the diameter of the fiber (which depends on φi for an
elliptical fiber). Note that the area over which the irradiance is measured is D dl where dl is an infinitesimal arc length along the fiber.

Given this definition, the scattering integral is written asL̄r (ωr ) = D

Z

S(ωi , ωr )Li (ωi ) cos θi dωi(1)Note that, unlike a surface where the integral extends over the upper hemisphere, this integral extends over the entire sphere. The
presence of D in this equation indicates that a thick fiber intercepts
more light, and therefore appears brighter from a distance, than athin fiber.

3

Scattering measurements

The experimental component of our study of hair was intended toprovide a qualitative and quantitative understanding of the phenomena that need to be explained by a scattering model for hair. In thispaper, we briefly outline these measurements; a full description of
our experiments will be provided in a future paper.
In our experiments we illuminated individual hairs with a narrow
beam and measured the scattered light in various directions using a
setup based on a four-axis goniometer that positioned a light source
and a CCD camera at arbitrary directions from the sample. We
used a focused beam to illuminate only a small length of hair, on
the order of 1 to 2 cm, and this illuminated segment defined thelength of hair being measured. This served to reduce the effects ofany variations in properties along the length of the hair.

3.1

Incidence plane

Stamm et al.’s and Bustard and Smith’s measurements were in theincidence plane, meaning that they observed only the 2D slice of thescattering function for which the source and detector are coplanarwith the fiber. To verify these earlier results, we made the samekind of measurements on samples of several different types of hairand on a synthetic fiber from a wig (Figure 4).
For each sample, we set the angle of incidence to 45◦ and measured scattering for varying outgoing angles. We performed thisexperiment twice, once with the illumination from the direction of
the root (with θi = −45◦), and again with the illumination from thedirection of the tip (with θi = 45◦). Each measurement was made inthree color bands across the visible wavelength range. Note that thesmall gap near the incidence angle in each plot is due to the camera
occluding the light source.

θh

15

0º
10

–15

30

15

0º
10

–15

–30 30

15

0º
10

15

–30 30

1

0º
10

15

0º
10

15

0º
40

–30

1

θd = 45º–15–30 301

θd = 30º–15–30 301

θd = 15º–15–30 301

θd = 5º–154θd = 60ºθd = 70º

Figure 5: Measurements of scattering in the incidence plane: scattering as a function of θh for several values of θd . The red and blue curves
correspond to red and blue wavelength bands; for this blond hair the very strong secondary highlight can be seen to be higher in the red. The
solid and dashed curves are two reciprocal measurements combined to extend the instrument’s range.
10θi = 0ºφ60

θi = 15º
30

1

φ

10

60
30

1

0º

φ

10

60

θi = 45º
30

1

0º

–30
–60

θi = 30º60

θi = 60º
30

1

0º

–30
–60

φ

10

φ

10

60
30

1

0º

–30
–60

φ

10

60
30

1

0º

–30

–60

θi = 75º

0º

–30
–60

–30
–60

Figure 6: A measurement showing the evolution of glints with incidence angle on a more nearly circular hair. For each θ the source is fixed
while the camera sweeps along the specular cone. A qualitative match to our shading model is plotted with dashed curves.

The most prominent feature in the scattering function is the specular highlight that occurs approximately when θr = −θi . For thesynthetic hair the highlight occurs exactly at the specular direction,
but for the real hairs it occurs 6 to 10 degrees toward the root, regardless from which direction the hair is illuminated. Additionally,
the blond hair, but not the black hair, has a color-dependent asymmetry in the lobe, with a side lobe appearing closer to the specular
direction than the primary highlight. A background of diffuse reflection is present that is larger and more colored for the blond hair.
In other measurements of this type we always observed the secondary highlight (except in black hair), but its strength and widthvaried considerably.
To study the effects of both θi and θr on the scattering function in
the incidence plane, we made a second measurement for one blondfiber. This measurement was parameterized in terms of θh and θd .
Each sweep was made with the source and detector fixed, thereby
holding θd fixed, and the hair rotating to vary θh .
A selection of the results of this full incidence-plane measurement are shown in Figure 5. In these plots, both the primary andsecondary highlight are visible, on either side of zero, which corresponds to the specular direction. This hair in this configuration
has a very strong secondary highlight that is actually stronger thanthe primary highlight for most angles. As the scattering angle increases, the secondary highlight fades out, while the primary highlight maintains more constant amplitude. Both peaks maintain approximately constant width. At high angles the primary highlight
gives way to a sharp peak that emerges very close to the specular
direction. This appears to be the same phenomenon that Stammet al. observed and named the equal-angle peak.

3.2

Normal plane

To begin to understand the full 3D complexities of the scatteringfunction, in particular the secondary highlight, we made anotherplanar 2D measurement, but with the plane of the source and detector oriented approximately perpendicular to the fiber; that is, in thenormal plane. Within this plane, the azimuths φi and φr both varied independently. To be more precise, the incident direction wasin the normal plane but the scattering direction was on the cone at
θr = 10◦. This is because the shift we have observed in the highlights due to the tilt of the cuticular scales would prevent us fromobserving the secondary highlight if we measured exactly in thenormal plane. In this and the following measurements, we show

only the red wavelength band, because that is where the secondary
highlight is most prominent.
The results are plotted as solid lines in Figure 7. We have organized the data in terms of φi and the relative azimuth φ , so eachpolar scattering diagram corresponds to a sweep of the camera withthe source fixed, and moving from diagram to diagram corresponds
to changing the light source angle (or equivalently, rotating the hairaround its axis). The wedges of missing data are caused by themechanical constraints of the goniometer.
This measurement shows that there are two bright out-of-plane
peaks. The strength and φ locations of the peaks vary from hair tohair. We refer to this peak as a “glint” because of the visual effectit produces.
Furthermore, the glints change considerably in brightness andposition as a function of φi , which shows that the hair is not rotationally symmetric. However, two symmetries are evident in the
figure. One is a 180◦ rotational symmetry in φi (diametrically opposite plots are similar); the other is by negating both φi and φ(vertical bilateral symmetry of the whole figure). Together thesesymmetries suggest that the hair has a 180◦ rotational symmetry
and is bilaterally symmetric in cross section.
When the pattern of glints is symmetric about the plane of incidence (φi = 0), it looks a lot like the caustic formation due torefraction in a sphere or cylinder with circular cross section. Moregenerally, the evolution of the peaks as the fiber rotates appears similar to the internal reflection from a transparent elliptical cylinder,
as we discuss in Section 6.1. A microscopic inspection of anotherfiber from the same individual reveals it has an eccentricity of 0.7:1.

3.3

3D hemispherical measurements

3.3.1

Changes in glints with angle of incidence

It is instructive to examine the secondary highlight as a function ofφ while varying the incidence angle. This shows whether and how
this highlight evolves as the inclination of the hair to the light sourcechanges. To make this measurement, we arrange for the source andcamera to be at approximately equal incidence angles θ on opposite
sides of the normal plane, and we measure the φ dependence bysweeping the camera along its cone while keeping the light source
fixed. The results are plotted as solid lines in Figure 6.
As in the previous measurement, the camera needs to be a bit
more toward the tip than the source in order to observe the maxi-

10

φ = 0º
15º
30º
45º
60º

φ

θi = 45º

primary
secondary

5

θr
1
.5
.1

φ

ROOT

TIP

1

φi

5
10

φ

φ = 75º

90º
120º
150º
165º

θr

50

100

Figure 8: A measurement of one hemisphere of the scattered distribution for a fixed incident direction with angle of incidence 45◦.
The primary highlight, the secondary highlight, and the transmitted
component are all visible.
Figure 7: A measurement of scattering in the normal plane from a
blond hair with substantial eccentricity. Illumination comes fromthe right in all plots; the green ellipses indicate the hair orientation.
Bright glints appear whose location and strength depend on the orientation of the hair. The green plots are the results of Monte Carloscattering simulation on the proposed fiber model.

• R: Shift of the primary specular peak toward the root. Thisdeviation is hypothesized to be due to the tilt in the scales on
the hair fiber.
• TT: A strong forward scattering component from light coloredhair. This causes blond, brown, gray, and white hair to lookvery bright when lit from behind.

mum glint intensity, so we maintain the condition that θi + θr = 10◦
as we vary θi .
The phenomenon shown by this measurement is that the azimuthat which the glints occur changes as a function of incidence angle, with the glints moving toward the incidence plane as the incidence moves from normal to grazing. They appear to merge aroundθ = 60◦. This behavior provides another clue to the structure: sucha merging transition is characteristic of the caustics formed by internal reflection in a cylinder [Mount et al. 1998]. The measured
eccentricity of this fiber is approximately 0.8:1.
3.3.2– A colored secondary peak shifted toward the tip fromthe white primary specular peak. In a head of hair,
this leads to the secondary highlight that is visible justabove the primary, sometimes appearing more as a colored fringe on the primary than a separate feature.
– The secondary highlight varies as a function of φ . It
generally contains two peaks, or glints. The locationsof the peaks depend on the angle of incidence, and they
converge to the incidence plane as θ increases. Furthermore, for hair that is not circular in cross section, thestrength and position of the glints depend strongly on
the angle of rotation of the hair around its axis. This
gives the hair fiber a distinctive sparkling appearance.
– Most of the energy in the secondary highlight is contained between the two glints.

Hemispherical scattering

These measurements extended the 1D incidence plane measurements by varying φr as well as θr , so that the scattering direction
varied over a full hemisphere. For a fixed direction of incidence,
this is a complete measurement of the scattering function (or halfof it, which suffices given bilateral symmetry—a reasonable if inexact assumption when we are not allowing φi to vary).
The results are shown in Figure 8. In this plot, several curves
for fixed values of φ are plotted together against θr . A primary
highlight appears at θr = 40◦, and a broad secondary highlight appears around 50◦–80◦. Moving from curve to curve, the primary
highlight maintains a fairly constant amplitude and shift from thespecular cone (marked by a blue line in the plot). The secondary
highlight becomes stronger, peaks around φ = 30◦, and then diesout. The forward-scattering half of the scattering function, shownwith a larger scale on the lower axes, is dominated by a strong
forward-scattered component that is not shifted appreciably fromthe specular cone.

3.4• TRT:

Summary

To summarize, our experiments measuring the scattering of lightfrom hair fibers revealed the following phenomena. Here, we categorize the phenomena in terms of the scattering modes of a dielectric cylinder (Section 2).

4

Theory of scattering from fibers

In this section we present background theory applicable to scattering from smooth cylinders. We show that because of the symmetry of a cylinder, the 4D scattering function can be factored into
a product of two 2D terms. One term, called M, captures the θdependence, and the other, called N, captures the φ dependence.
We also derive N for the case of a circular cylinder. These twoparts form the basis for the complete, practical model developed in
the next section, in which we introduce expressions for M, an approximation to extend N to elliptical cylinders, and procedures for
efficiently approximating N.

4.1

Scattering from cylinders

For fibers formed by translating a fixed cross section along a single
axis we can separate the 4D scattering function to a product of 2D

R (p = 0)

terms. We make use of properties that have been used in previous
work on scattering from fibers [Marcuse 1974; Adler et al. 1998;
Mount et al. 1998]:• A ray that enters a dielectric cylinder at a particular inclination to the axis will always exit at the same inclination, regardless of the sequence of refractions and reflections it undergoes.
This means that a bundle of parallel incident rays coming fromthe direction ωi will produce a collection of scattered rays
whose directions lie on the cone centered on the fiber axis andcontaining −ωi (Appendix A). The directions of the refracted
rays inside the hair also lie on a cone. This effectively reduces
the scattering function from 4D to 3D, because scattering only
occurs when θr = −θi . Of course, this equality holds only
approximately if the surface is rough, and it depends on theassumption that there is no volumetric scattering inside thecylinder.
• The dependence of the scattered distribution on φr can be analyzed by examining only the projection into a plane perpendicular to the hair.
It’s easy to see that the R component can be computed justfrom the projection—the incident and reflected vectors are arranged symmetrically across the surface normal in the projection just as they are in 3D, so a mirror reflection from the 3D
cylinder remains a mirror reflection in the projection.
The same holds for refracted rays, as a consequence of Bravais’s law (derived in Appendix B), which states that if theincident and transmitted vectors at a dielectric interface areprojected into a plane containing the surface normal, the projected vectors still obey Snell’s law, but with the index of refraction η replaced by η 0 (η, θ ) > η. Note that η 0 depends
only on the inclination θ of the incident ray out of the projection plane, and it is independent of the angles in the projection plane. Thus, since we know that all rays originating froma particular incident direction maintain the same inclination
to the normal plane, a 2D analysis that works in the normal
plane for arbitrary refractive index suffices to describe the 3D
scattering function.
These two facts allow us to write the scattering function for asmooth cylinder of any cross section as:
S(φi , θi ; φr , θr ) = δ (θr + θi )N(η 0 (θ ); φi , φr )/ cos2 θ
The first factor says that scattering only occurs in the specular cone;
the second factor N represents the azimuthal scattering function;and cos2 θ accounts for the projected solid angle of the specularcone. We write simply θ to emphasize that θi = θr . Since θ influences N only indirectly through η 0 , we can derive N independently
using a 2D analysis. In the next section, we will generalize thisexpression to
S(φi , θi ; φr , θr ) = M(θi , θr )N(η 0 (θd ); φi , φr )/ cos2 θd ;(2)that is, we will allow for a more complicated dependence on θ . We
will refer to M as the longitudinal scattering function.

4.2

Scattering from a circular cross section

In this section we derive N assuming a circular cross section. Withcircular symmetry, N depends only on φ , the difference between φiand φr , so the 2D azimuthal scattering function is further reduced to
1D, and in this context we will write N with a single angular argument. Scattering from a dielectric circle is well studied, having been

dφR

γi

dh

γi
γt

h

γt
γt

γi
TT (p = 1)

γt

dφTT

dφTRT

γi
TRT (p = 2)

Figure 9: Geometry for scattering from a circular cross section.
first analyzed by Descartes to explain the formation of the rainbow
in a rain shower (for a full discussion, see, e.g., Humphreys [1964]).
Note that the effective refractive index η 0 is a parameter to the entire discussion that follows, so for clarity it is left as an implicit
parameter when it is not being discussed.
We can determine the far-field intensity by tracing rays as theyrefract through a circle. Consider a ray incident on a unit circle
at an offset −1 < h < 1 from its center (see Figure 9). The angle
sin γi = h is the angle of incidence, and η 0 sin γt = h is the angle of
the refracted ray. By following a ray entering a circle as it refracts
and reflects, we can calculate the exit angle φ (h).
From Figure 9 it can be seen that the incident ray deviates by
γt − γi as it enters the circle, by π + 2γt at each internal reflection,and by γt − γi again as it exits, for a total ofφ (p, h) = 2pγt − 2γi + pπ

(3)where p is the number of internal path segments.1 This expression encompasses all three scattering modes we are interested in,for different values of p: surface reflection R (p = 0), refractive
transmission TT (p = 1), and internal reflection TRT (p = 2). For
the rest of this analysis, we will ignore p > 2 terms.
The above ray tracing calculation parameterizes paths by h.
However, for scattering calculations we often need to find all thepaths that contribute to scattering in a given direction φ . The h values for these paths are found by solving for the roots of the function
φ (p, h) − φ = 0. We denote these roots by the function h(p, r, φ ),where different values of r denote different roots. In the p = 0 and
p = 1 cases, there is a single root, and thus a single path. However,for the p = 2 case, there may be one or three roots and hence oneor three paths.
Since the function φ is smooth, the transition from a single exit
ray to three exit rays represents a fold in φ (2, h) and this fold occurs
when dφ
dh = 0. Descartes first showed that this fold occurs when:h2 = (4 − (η 0 )2 )/3.

(4)

This expression is symmetric, and predicts that two symmetric extrema exist. [Humphreys 1964]
Now that we have the expression (3) to relate φ and h, we cancompute the intensity of the scattered light by using the principle ofenergy conservation. When curve irradiance Ē illuminates a fiber,
uniform irradiance E(h) = Ē/2 falls on the cross section across its
1 Note that pπ accounts for p − 1 internal reflections plus the reversal in
direction from ωi , which points outward.

width.2 Setting aside attenuation for the moment, power from asmall interval dh in the incident beam is scattered into an angularinterval dφ in the exitant intensity distribution (Figure 9).
L̄(φ (h)) dφ = E(h) dh = (Ē/2) dhorL̄(φ (h)) = 2

(5)

dφ −1
Ē
dh

Attenuation by absorption and reflection

In this section we add terms for volume absorption in the fiber interior and Fresnel reflection at the interfaces.
For absorption we need to know the path length inside the fiber.
Applying the law of cosines to the triangles in Figure 9 gives thelength of each internal path segment as 2 + 2 cos(2γt ) times thehair’s radius. Let σa be the volume absorption per unit length,with the unit length defined to be the radius of the hair. Theneach of the p segments inside the hair contributes a factor ofT (σa , h) = exp(−2σa (1 + cos(2γt ))) to the absorption.
Now we can introduce an attenuation factor in front of the intensity contributed by a pathL̄(φ ) = A(p, h) 2

dφ −1
Ē
dh

(7)

where

S(φi , θi ; φr ,θr ) =MR (θh )NR (η 0 (η, θd ); φ )/ cos2 θd +MT T (θh )NT T (η 0 (η, θd ); φ )/ cos2 θd +

2A(p, h) = 1 − F(η, γi ) F(1/η, γt ) p−1 T (σa , h) p .
Just as Bravais’s law lets the normal-incidence geometry factor
work for oblique incidence, we can generalize A to oblique incidence simply by changing the parameters. Appendix B shows thatthe correct Fresnel factors may be computed using the usual Fresnel formulas and two virtual indices of refraction, η 0 (θ ) and η 00 (θ ).
Similarly, because all the internal path segments have the same inclination θt to the axis, they lengthen by a factor of 1/ cos θt , so
substituting σa0 (θ ) = σa / cos θt for σa completes the generalization
to 3D:
A(0, h) = F(η 0 , η 00 , γi )
2A(p, h) = 1 − F(η 0 , η 00 , γi ) F



1 1
,
,γ
η 0 η 00 t

 p−1

T (σa0 , h) p

This equation derives from Equation 2 but includes generalizations
that will allow us to introduce several useful approximations in thenext section. The three terms MR , MT T , and MT RT are three separate longitudinal scattering functions. They allow the effects of the
cuticle scales, which are different for the three modes, to be simulated. The three terms NR , NT T , and NT RT are the three modes of
the azimuthal scattering function. The first two of these,
NR (φ ) = N p (0, φ )NT T (φ ) = N p (1, φ ),are defined directly from Equation 8, but the computation of NT RT
will be discussed in more detail in Section 5.2.2. The function η ∗lets the effective index of refraction for NT RT depend on φh , whichallows the approximation for elliptical cross sections presented in
Section 5.2.3.
Note that only the angles θd , φ , θh , and φh appear on the righthand side of Equation 9. Writing the Ms and Ns in terms of sums
and differences in this way emphasizes the natural symmetry of S
and will make it easy to see that the model is reciprocal.

N(φ ) = ∑ N p (p, φ )p

N p (p, φ ) = ∑ A(p, h(p, r, φ )) 2r−1
dφ(p, h(p, r, φ ))dh(8)where the sum is over all the paths of different types (p) including
those with multiple roots (r).
2 For the unit radius circle Ē is spread out across the fiber width of 2. The
fiber diameter does not affect S.

A practical shading model for human hair

In this section we extend the results of the previous section into acomplete, practical model suitable for integration into a renderer.
To do this we introduce empirical functions for the longitudinalscattering functions (MR , etc.), which the theory of smooth cylinders does not predict. We also describe our approach to computing
NT RT in the presence of surface roughness and an approximation toh(p, r, φ ) that is used to compute NR , NT T , and NT RT . Finally, we
introduce convenient parameters for controlling the model.
It is important to recognize that the purpose of the model givenin this section is to capture the phenomena that are important in
rendering. The form of the model and the basic energy distribution
come from the theory presented in the previous section, but there
are many details for which a simple analysis is not available, andwe fill these in empirically.

5.1

The complete normal-plane scattering function is thus

(9)

MT RT (θh )NT RT (η 0 (η ∗ (φh ), θd ); φ )/ cos2 θd .

5
A(0, h) = F(η, γi )Summary

The theory developed in this section predicts that the light from a
directional beam will be scattered into a perfect cone and that the
distribution around the cone is a sum of three distinct scatteringmodes. This leads to the following form for our new scatteringmodel:(6)

This equation implies that intensity goes to infinity at the fold in thefunction φ . This intensity singularity is called a caustic.
In the case of a circular cylinder, the pair of caustics cause theglints. Since in the cylinder the Bravais index η 0 increases with θ ,the two caustics will move closer to the plane as the incident lightbecomes more oblique to the hair. When η 0 reaches 2, the caustic
ceases to form. This caustic merge transition can be seen in themeasurement plotted in Figure 6.

4.3

4.4

The longitudinal scattering function M

The analysis of a smooth cylinder predicts that the reflected lightwill stay exactly in the specular cone. Our hair model includes
two features that cause deviations from this behavior. First, theinterfaces are rough. As rays propagate through the cylinder, theroughness will cause their directions to deviate randomly from thedirections predicted by the model. The overall effect is to blur the
scattered distribution, with the different scattering modes blurred todifferent degrees. Second, the cuticle scales cause a tilt of the surface normals relative to the ideal cylinder, which causes the scattered lobes not to be centered on the specular cone.
Figure 1 illustrates how we can expect the positions of the lobes
to shift: the R mode is displaced by 2α toward the root; the TRT

Parameter

PurposeFiber propertiesη
index of refraction
σaabsorption coefficient (R, G, B)aeccentricity

glint scale factor
azimuthal width of caustic
fade range for caustic merge
caustic intensity limit

30º

0.5 to 5
10◦ to 25◦
0.2 to 0.4
0.5

Table 1: All the parameters of the shading model.
mode is displaced more than that far toward the tip; and the TTmode is displaced toward the tip.3 We approximate the effects ofthe scales by shifting the mean of M differently for each lobe.
MR (θh ) = g(βR ; θh − αR )MT T (θh ) = g(βT T ; θh − αT T )MT RT (θh ) = g(βT RT ; θh − αT RT )where g(β , x) is a unit-integral, zero-mean lobe function with width
β . In our implementation we used a normalized Gaussian function
with standard deviation β .

–30º
5

1.55
0.2 to ∞0.85 to 1

Surface properties
αR
longitudinal shift: R lobe
−10◦ to −5◦
αT T
longitudinal shift: TT lobe
−αR /2αT RTlongitudinal shift: TRT lobe
−3αR /2
βR
longitudinal width (stdev.): R lobe
5◦ to 10◦βT T
longitudinal width (stdev.): TT lobe
βR /2
βT RT
longitudinal width (stdev.): TRT lobe
2βR
Glints
kG
wc
∆η 0
∆hM

10

Typical values

θr

60º
1
.5
.1

Figure 10: A plot of the shading model, under the same geometricconditions used in Figure 4. The parameters have been adjusted tomatch the behavior of the blond hair (top center plot in that figure).
5.2.2

Approximation for TRT

Our theory, which is based on smooth surfaces, predicts that thecaustics in the TRT component produce singularities in S with infinite intensity. Since this is unrealistic, particularly in the presenceof surface roughness, we remove the caustic from NT RT and replace it with a smooth lobe centered at the location of the caustic.
The width of this lobe simulates the blurring of the caustic due toroughness.
The angle φc at which the caustic appears can be computed from
Descartes’s formula (Equation 4) when it is defined, but for incidence angles past the caustic merge it is undefined. To maintaincontinuity we continue to insert the caustic at φc = 0 and fade itout smoothly over a short range of incidence angles past the merge.
It is simplest to precisely describe this component by giving theprocedure used to compute it:
function NT RT (θ , φ ; wc , kG , ∆η 0 , ∆hM )if (η 0 (θ ) < 2)0Compute hc , φc usingq η (θ ) in (4)2

5.2

The azimuthal scattering function N

In this section we introduce three approximations to the azimuthal
scattering function. The first is an efficient method for solving forpaths, the second is a phenomenological model for the secondary
highlights, and the third is a simple model for glints.
5.2.1

Solving for paths

In order to find a path through a circle, we must solve Equation 3 forthe value or values of h(p, r, φ ) that result in the scattering angle φ .
Because the formula for φ involves Snell’s Law and hence arcsines,
it is computationally expensive to solve for h exactly.
To simplify the angular dependence, we approximate Snell’s law
with the cubic polynomial that matches the value and derivative of
the exact expression at ±90◦:
γt =

3c
4c
γ − γ3π i π3 i

where c = sin−1 (1/η 0 ) (in this expression angles γ are measured inradians). The maximum approximation error is less than 0.75◦ forη > 1.5.
With this approximation, φ is a cubic in γi :


6pc
8pc
φˆ (p, γi ) =− 2 γi − 3 γi3 + pπ.
(10)ππThe roots of this equation yield γi from which h is easily found. For
the R and TT cases there will always be exactly one root, whereas
the TRT case may have one or three roots.
3 For the TT and TRT components the shifts depend on the incidence
angle, and for all three components the shift depends on φ , but we ignore
these effects.

∆h = min(∆hM , 2 2wc /| ddhφ2 (hc )|)
t =1
else
φc = 0∆h = ∆hM
t = smoothstep(2, 2 + ∆η 0 , η 0 (θ ))L = N p (2, φ )L = L · (1 − tg(φ − φc , wc )/g(0, wc ))L = L · (1 − tg(φ + φc , wc )/g(0, wc ))L = L + tkG A(2, θ , φ )∆h(g(φ − φc , wc ) + g(φ + φc , wc ))return L

In this procedure we roughly approximate the power that is removed from the intensity distribution by using d 2 φ /dh2 to estimate
the size ∆h of the interval in h that maps within wc of φc . The
limit ∆hM is necessary because d 2 φ /dh2 goes to zero at the caustic
merge. The function smoothstep(a, b, x) is 1 for x < a, 0 for x > b,and smooth in between.
The user-adjustable parameters to this function are wc , the widthof the blur for the caustics, kG , a factor to adjust the strength of theglints, ∆η 0 , the range of η 0 over which to fade out the caustic after
the merge, and ∆hM , a limit on the caustic power.
5.2.3

Approximation for eccentricity

When the fiber’s cross section is elliptical rather than circular, asimple analytical solution for N is not available. However, the appearance of the TRT component is significantly affected by evenmild eccentricity: the angle φc at which the caustics appear changes
by ±100% over the range from a = 0.85 to a = 1/0.85. Because
eccentricities in this range are very common, it is important to approximate the effect of eccentricity on the TRT component. We do
this by using the function η ∗ in Equation 9.

Changing refractive index has effects that are qualitatively similar to changing eccentricity. One case for which it is simple toanalyze this relationship in the TRT component is paraxial paths
that are symmetric across the major axis. For these paths φi = −φrand the reflection occurs on the major axis. If γt0 is the angle ofincidence at the internal reflection it can be shown thatdφr
= 2a2 (η − 1) − η.
dγt0
We can use this first-order approximation to roughly match the behavior of the circular model to the elliptical one by using the index
of refraction η ∗ = 2a2 (η − 1) − η + 2 in the circular model.
This gives us an approximation for how the TRT component of
the scattering function changes with mild eccentricity when φh isaligned with one of the principal axes. To incorporate this approximation into the practical model, we simply interpolate sinusoidally
to define η ∗ for all intermediate values of φh :
η1∗ = 2(η − 1)a2 − η + 2
η2∗ = 2(η − 1)a−2 − η + 2
1(η1∗ + η2∗ ) + cos(2φh )(η1∗ − η2∗ )η ∗ (φh ) =2
As shown in Equation 9, we simulate eccentricity simply bypassing the refractive index η ∗ to the code that computes NT RT
using the formulas derived from the circular case.

5.3

Summary

We have now accounted for all the terms that appear in Equation 9.
The three longitudinal scattering functions MR , MT T , and MT RT
are Gaussians. The azimuthal scattering functions NR and NT T are
evaluated directly from Equation 8, using Equation 10 to computeh. The azimuthal scattering function NT RT is first computed using
Equations 8 and 10 with the refractive index η ∗ , summing over the
possibly multiple values of h, then modified by the procedure in
Section 5.2.2 to smooth out the caustics.
Depending on the type of renderer, it may be helpful in practical use to include a small amount of diffuse reflection. Although
our model does not include a diffuse component, one can easily be
added by introducing a constant term MD into Equation 9.
All the model’s parameters are collected for reference in Table 1.

6

Results

6.1

Validation

In this section we provide evidence that our assumptions about thestructure of hair are valid and that the shading model can qualitatively reproduce the behaviors we have seen in the measurements.
Validating the hair structure model against measurements. The
normal-plane measurement in Figure 7, which best illustrates theglints, behaves entirely consistently with the model of a hair as anelliptical cylinder. To show this, we have used a Monte Carlo procedure to compute the scattering due to the TRT mode in a roughelliptical cylinder with eccentricity matching the measured value
for this individual. The figure shows the results of this numericalexperiment plotted with the results of the physical experiment. The
match in both the angular and the intensity variation in the peaks
is striking, providing clear evidence both for the explanation of theglints as internal reflection caustics and for the model of the hair’s
cross section as an ellipse. The real data of course includes some
contributions from the other scattering modes, which is noticeable
in forward scattering, where the TRT mode does not contribute.

Figure 11: A photograph [Gray 1997] showing brown hair with aprominent secondary highlight caused by a small light source near
the camera. Glints help give the hair its distinctive natural texture.
Validating the shading model against measurements. In order todemonstrate the ability of our shading model to qualitatively matchthe behavior of real hair, we have set the parameters to fit two of
the measurements in Section 3. In Figure 10 we show a plot of themodel adjusted to behave like the blond hair that was measured in
the incidence plane (Figure 4, top center). Note that both the small
side lobe and the broader lobe toward grazing are predicted.
The measurement in Figure 6 shows a caustic merge as inclination changes. The shading model is plotted as a dashed line
along with the data, with a set of parameters that describes thismore mildly eccentric hair. The behavior of two glints merging anddisappearing matches qualitatively.

6.2

Renderings

We have implemented the shading model described in the previous section and used it to simulate light reflection from procedurally generated assemblies of hair fibers. We have used thecommercial “Sasquatch” rendering software (Worley Laboratories,
http://www.worley.com) to model the hair, and each model consisted of 50,000–100,000 spline curves, corresponding to the average number of hairs on a human head. The shading parameters wereconstant per fiber, but in some cases (noted below) we used a noisefunction on the scalp to assign a random absorption to each fiber.
The azimuthal orientation of the elliptical fibers was random fromfiber to fiber but smooth along each fiber. All renderings have η =1.55, and ∆hM = 0.5. The parameters αT T and αT RT are definedfrom αR by the ratios recommended in Table 1.
We used ray tracing to render the images, accounting for occlusion and shadowing but not multiple scattering between the hairs.
We represented the hairs as flat ribbons that always face the ray
direction, with a width equal to the hair diameter. This is more
efficient than using cylinders, and since our scattering function abstracts away all interactions below the level of strands it is appropriate.
In the following results, computation was dominated by ray intersection so that even though our model is substantially more expensive to compute than Kajiya’s very simple model, the increase
in rendering time was modest. For instance, switching to the newmodel increased the rendering time for Figure 12 from 6 minutes to
8 minutes on a two-processor 1GHz Pentium III system.
The first set of renderings is inspired by the photograph shownin Figure 11. Because the hair is combed into a fairly flat surface,Figure 12: A comparison of Kajiya and Kay’s model (left) under a single point source, our proposed model (center) with the same lighting,and the hair from the photograph in Figure 11 (removed from context to simplify the comparison). The Kajiya model’s diffuse term results ina flat appearance, while the secondary highlight in our model correctly captures the colored shading of the real hair.

the secondary highlight is quite prominent; there is very little diffuse scattering, as evidenced by the dark appearance at the top andbottom of the hair.
We constructed an idealized geometric model to roughly matchthe overall shape of the hair, without attempting to reproduce theplacement of the fibers exactly. In Figure 12 we show a renderingof this model using settings appropriate for mildly elliptical brown
hair. The result of rendering the same model in the same lighting
using Kajiya’s model is shown for comparison. The diffuse component required to produce an overall brown color does not adequately
predict the variation in shading that occurs in the real hair, where
the brown color comes from the secondary highlight. The sameprocedural color variation is present in both images; the additional
texture from the proposed model is due to eccentricity. The parameters for this images are αR = 3◦, βR = 14◦, βT T = 8◦, βT RT = 22◦,
σa = (0.432, 0.612, 0.98) – (0.36, 0.72, 0.94), kG = 0.4, wc = 1.5◦,
and a = 0.9.
To illustrate the components of the model and the effect of eccentricity, Figure 14 shows the two highlights separately and a comparison between circular and elliptical hair with all other parameters
held constant. The parameters for these images are α = 3◦, βR = 8◦,
βT T = 6◦, βT RT = 15◦, σa = (0.44, 0.64, 0.9), kG = 0.4, wc = 1.5◦,
and a = 1 and 0.9.
To further demonstrate the capabilities of the model we include
additional results for varying light direction in Figure 15 and forcomplex illumination (using the method of Agarwal et al. [2003])
in Figure 13.

7

comes primarily from diffuse scattering is not supported by the evidence in the case of fibers, and our images show that this assumption limits the quality of results from current hair shading models.
The shading model we have introduced could be made more accurate in future work. Two avenues for improvement are eccentricity and surface roughness. It would be useful to handle more extreme deviations from circularity, in order to allow us to accurately
model more types of hair. Also, the current model approximates theeffects of eccentricity on the ray density factor but not the Fresnel
factor, which limits the quantitative accuracy of the energy distribution. Our model accounts for the effects of surface roughness
and the scale angle in the simplest way that is consistent with theappearance phenomena. The true effects of these surface features
are complicated and dependent on the scattering geometry. A moreexact model of these effects could lead to a scattering model that ismore accurate for grazing angles.
Our analysis applies to many kinds of transparent fibers, so it has
implications for rendering animal hair, synthetic hair or fur, and all
kinds of cloth.

Conclusions

In this paper we have looked at light scattering from human hairin reality, in theory, and in rendering practice. Our measurements
show that the scattering behavior of individual hair fibers differs
significantly from what has previously been assumed in shading
models; the theory of scattering from dielectric cylinders explains
the sources of these differences; and we present a shading modelthat captures the relevant visual effects, including those due to eccentricity, in a practical, adjustable model for use in rendering.
The most important implication of this work is that realistic rendering of hair needs to account for the multiple modes of scattering that occur in cylinders. The standard assumption that the color

Figure 13: A hair model (62K fibers) illuminated by a complex
environment map.

(a)

(b)

(c)

(d)

(e)

Figure 14: (a–c) The components of our model for circular hair (a = 1): (a) the R component (primary highlight) alone; (b) the TRT
component (secondary highlight) alone; (c) the full model. (d) The full model for a = 0.9, showing texture due to glints; (e) an enlargementof (d).

Figure 15: A hair model under different illumination angles. As the relative strengths of the primary and secondary highlights change, theyalter the color, brightness, and sheen of the hair.

Acknowledgements
The authors would like to thank Heidi Marschner, Joyce Ma, andJames Davis for donating the hairs we measured, François Guimbretière for the presentation idea for Figure 7, and Dr. John Gray forthe use of the photograph in Figure 11. This research was supported
by an NSF ITR grant (IIS-0085864), and the third author was supported by a 3Com Corporation Stanford Graduate Fellowship andan NSF Graduate Fellowship. The scattering measurements were
performed on the Stanford Spherical Gantry, which was designed
by Marc Levoy, Brian Curless, and Szymon Rusinkiewicz and constructed using funds from Interval Research.

References
A DLER , C. L., L OCK , J. A., AND S TONE , B. R. 1998. Rainbow scattering
by a cylinder with a nearly elliptical cross section. Applied Optics 37, 9,1540–1550.
AGARWAL , S., R AMAMOORTHI , R., B ELONGIE , S., AND J ENSEN , H. W.
2003. Structured importance sampling of environment maps. ACM
Transactions on Graphics 22, 3.
B USTARD , H., AND S MITH , R. 1991. Investigation into the scattering oflight by human hair. Applied Optics 24, 30, 3485–3491.
C HANG , J., J IN , J., AND Y U , Y. 2002. A practical model for hair mutual interactions. In Proceedings of the ACM SIGGRAPH Symposium on
Computer Animation, ACM Press, 73–80.
C HEN , L.-H., S AEYOR , S., D OHI , H., AND I SHIZUKA , M. 1999. A
system of 3D hair style synthesis based on the wisp model. The VisualComputer 15, 4, 159–170.
G OLDMAN , D. B. 1997. Fake fur rendering. In Proceedings of SIGGRAPH
1997, ACM Press, 127–134.
G RAY, J. 1997. World of Hair. Hairdressing Training Board.

H ADAP, S., AND M AGNENAT-T HALMANN , N. 2000. Interactive hair styler
based on fluid flow. In Computer Animation and Simulation 2000, Eurographics, 87–99.
H UMPHREYS , W. J. 1964. Physics of the Air. Dover, London.
K AJIYA , J. T., AND K AY, T. L. 1989. Rendering fur with three dimensional
textures. In Proceedings of SIGGRAPH 1989, 271–280.
K AJIYA , J. T. 1985. Anisotropic reflection models. In Proceedings of
SIGGRAPH 85, ACM Press, 15–21.
K IM , T.-Y., AND N EUMANN , U. 2000. A thin shell volume for modelinghuman hair. In Computer Animation 2000, IEEE CS Press, 104–111.
K IM , T.-Y., AND N EUMANN , U. 2001. Opacity shadow maps. In Rendering Techniques 2001: Proceedings of Eurographics Rendering Workshop, Springer Verlag Wien, 177–182.
K IM , T.-Y., AND N EUMANN , U. 2002. Interactive multiresolution hairmodeling and editing. ACM Transactions on Graphics 21, 3, 620–629.
K IM , T.-Y. 2002. Modeling, Rendering, and Animating Human Hair. PhD
thesis, University of Southern California.
KONG , W., AND NAKAJIMA , M. 2000. Hair rendering by jittering andpseudo shadow. In Computer Graphics International 2000, 287–294.
L OKOVIC , T., AND V EACH , E. 2000. Deep shadow maps. In Proceedingsof SIGGRAPH 2000, ACM Press, 385–392.
M AGNENAT-T HALMANN , N., H ADAP, S., AND K ALRA , P. 2002. State of
the art in hair simulation. In International Workshop on Human Modeling and Animation, Korea Computer Graphics Society, 3–9.
M ARCUSE , D. 1974. Light scattering from elliptical fibers. Applied Optics
13, 1903–1905.
M OUNT, C. M., T HIESSEN , D. B., AND M ARSTON , P. L. 1998. Scatteringobservations for tilted transparent fibers. Applied Optics 37, 9, 1534–
1539.
P OULIN , P., AND F OURNIER , A. 1990. A model for anisotropic reflection.
In Proceedings of SIGGRAPH 90, 273–282.

(a)

(b)
vi
n′
n

vi
vi′

vi

n′

n

(c)

hi

si

ht

n

st′
st

si′

vt
vt′

vr′

γ

st′= si′/ η
s t = si / η

δ

vr

vt

vt′

Figure 16: Reflected (a) and refracted (b) vectors from a cylinder lie in a cone; derivation of Bravais’s law.

ROBBINS , C. R. 1994. Chemical and Physical Behavior of Human Hair,
third ed. Springer-Verlag, New York.
S OURIN , A., PASKO , A., AND S AVCHENKO , V. 1996. Using real functions
with application to hair modelling. Computers & Graphics 20, 1, 11–19.
S TAMM , R. F., G ARCIA , M. L., AND F UCHS , J. J. 1977. The optical properties of human hair I. fundamental considerations and goniophotometer
curves. J. Soc. Cosmet. Chem., 28, 571–600.
T RICKER , R. A. R. 1970. Introduction to Meteorological Optics. Mills &Boon, London.
WATANABE , Y., AND S UENAGA , Y. 1992. A trigonal prism-based method
for hair image generation. IEEE Computer Graphics & Applications 12,
1, 47–53.

A

Loci of reflections

In this appendix, the geometry of reflection and refraction from acylindrical fiber is analyzed. The set of normal vectors to a cylindrical fiber lie in a plane perpendicular to the axis of the fiber; thisplane is called the normal plane. The geometry is most easily visualized using a unit sphere; each point on the unit sphere represents
a different direction. On this sphere, the axis of the fiber is at thenorth pole and the set of normal vectors is contained in the equator(see Figure 16a).
To analyze reflection for a given incident direction vi , consider
a particular surface normal n and the corresponding reflection vector vr . The law of reflection states that vi , n, and vr are coplanarand n · vi = n · vr . By finding similar triangles, the reflection rulealso forces the distances of vi and vr from the equator to be equal.
This equidistance condition must hold for any normal vector. Thus,
the set of reflection vectors formed for a given incident direction
must lie on a horizontal circle at constant height that is at the same
distance from the equator as vi . This is a formal proof of the wellknown result that, in 3D, the set of reflections from a cylinder form
a cone.
Figure 16b is analogous to the previous figure, but shows therefracted directions instead. The rule for refraction (Snell’s law)is similar to reflection. First, vi , n, and vt are coplanar. Second,
Snell’s Law says that η sin θt = sin θi , and using similar triangles,
this law implies that the perpendicular distances to the equator, hi
and ht , also have the ratio η. This means the transmitted vectors alllie on a circle that is a factor of η closer to the equator than vi .
By the same argument, the rays that refract again on the way outof the cylinder will obey the equidistance condition. Furthermore,internal reflections are no different than external reflections. The
conclusion is that all rays that exit a cylindrical fiber must lie in thesame cone as the reflected rays.

B

The Bravais index

A cylindrical fiber is formed by sweeping a cross section along an
axis. Bravais has shown that the optics of a 3D cylindrical fiber may
be reduced to the 2D analysis of the optics of its cross section. The
incident direction is first projected into the cross-section plane andthen reflected and refracted in 2D. However, to model the correct
physics in 2D, the index of refraction must be changed as a functionof the angle of incidence. Bravais’s law is often used to analyze
refraction through crystals (as in ice haloes [Tricker 1970]).
Figure 16c shows the geometry of the refraction across a horizontal surface. The direction vi is refracted in the direction vt . To
reduce this to a 2D problem, the vectors vi and vt are projected onto
the vertical plane containing the normal. The projected vectors arelabeled v0i and vt0 . The goal is to find an effective index of refraction, η 0 , such that if v0i is refracted at the boundary it will yield vt0 .
First note that Snell’s Law states that sin θi = η sin θt . That is, thelength of the projection of vi onto the horizontal plane is η times
the length of the projection of vt . Since vi and vt lie in a plane,
the two triangles shown in the horizontal plane are similar. The
length of each edge of the left (incident-side) triangle is η times
the length of the corresponding edge in the right (transmitted-side)
triangle. Second, an equivalent Snell’s Law is obeyed in the vertical normal plane. Here, sin θi0 = η 0 sin θt0 . But sin θi0 = sin θi cos γ,
where γ is the angle between the incident vector and its projection
on the normal plane. Similarly, sin θt0 = sin θt cos δ , where δ is theangle between the transmitted direction and its projection on thenormal plane. Substituting, yields sin θi cos γ = η 0 sin θt cos δ , andsubstituting Snell’s Law into this equation yields η cos γ = η 0 cos δ ,Equivalently,

q
0
2
2
η (γ) = η − sin γ cos γ.
The implication of this is that for an incident ray that makes anangle γ with the projection plane, we can compute the projection of
the refracted direction from the projection of the incident direction
using the usual Snell’s law but substituting the effective index ofrefraction η 0 .
The Fresnel attenuation at the interface may be computed in asimilar way. From Fresnel’s formulas, it is straightforward to verify
that the factor for oblique incidence may be computed from the
projected angles by using η 0 for the perpendicular component of
the reflectance and a second index,.qη 00 (γ) = η 2 cos γ
η 2 − sin2 γ
for the parallel component. We will use the notation F(η 0 , η 00 , γ)for this slightly generalized Fresnel function.