(Investigative Ophthalmology and Visual Science. 2003;44:4693-4697.)
© 2003 by The Association for Research in Vision and Ophthalmology, Inc.
doi:10.1167/iovs.03-0260
Computer-Synthesis of an Interference Color Chart of Human Tear Lipid Layer, by a Colorimetric Approach
Eiki Goto,1,2,3
Murat Dogru,1
Takashi Kojima,1,4 and
Kazuo Tsubota1,2
1From the Department of Ophthalmology, Tokyo Dental College, Chiba, Japan; the
2Department of Ophthalmology, Keio University School of Medicine, Tokyo, Japan; the
3Tokyo Dry Eye Center, Iidabashi Eye Clinic, Tokyo, Japan; and the
4Department of Ophthalmology, Social Insurance Chukyo Hospital, Nagoya, Japan.
 |
Abstract
|
|---|
PURPOSE. To synthesize an interference color chart for the specific tear lipid layer interference camera, DR-1, for the conversion of the tear film lipid layer thickness into color graphic informationthat is, for the quantification of the tear interference imageby a colorimetric approach.
METHODS. Because the color of the tear lipid layer interference image is visualized by a white light source interference phenomena, to produce
-corrected red, green, and blue (RGB) values of a specific interference color at a certain tear lipid film thickness, XYZ tristimulus values of the Commission Internationale de lÉclairage (CIE) were obtained. XYZ tristimulus values were calculated from the light source spectrum of the DR-1 camera, the color-matching function of CIE, and the reflectance of the tear interference image in wavelengths ranging from 380 to 780 nm. These calculated interference colors were synthesized ranging from 0 nm to 1000 nm of lipid film thickness to produce a color chart. The applicability of the new color chart in the analysis of the lipid layer thickness was tested on a healthy control subject with normal tear function and a patient with dry eye who had aqueous tear deficiency and meibomian gland obstruction.
RESULTS. The specific tear interference color chart for the DR-1 camera was obtained with RGB and XYZ tristimulus values. The interference chart ranged from 0 to approximately the 5th interference order. The interference colors from clinical DR-1 images could be converted to lipid thickness data by using the color chart system.
CONCLUSIONS. A new tear interference color chart was developed in this study, which may be of benefit in converting tear interference color information to data describing the thickness of the tear film lipid layer.
Tear lipid layer interferometry is a noninvasive method to visualize the lucent surface lipid layer of the tear film.1 2 3 Tear interference images have been reported in observing the surface phenomena of the tear film,4 5 6 7 8 using the principle first reported by Sir Isaac Newton.9 Existence of the interference phenomena indicates the presence of the optical path difference.10 The interference phenomena inferred from the tear film indicate the presence of a thin film, which is the superficial lipid layer.5 11 The coloration of these interference images has been reported to depend on the lipid films thickness.5 11 This thickness has been considered to affect tear evaporation and lubrication in blinking.12 13 14 Thus, the quantification of the lipid layers thickness is essential in tear interferometry, for objective assessment of tears. However, the translation of the interference color data into film thickness data has not been successfully applied quantitatively, nor has the quantification of the tear interference images itself been achieved.
During examinations for dry eye, interferometry has been used with a color-comparison table to assume lipid layer thickness15 16 17 or with the semiquantitative severity grade scoring system for Sjögrens syndrome18 and dry eye syndrome.19 Thus far, with the color table and the severity grading system, there has been some confusion in interpreting the results of the image. Also, semiquantitative grading was unsuitable for evaluating other than aqueous-deficient dry eye status and might not detect mild cases of dry eye.
Several trials of quantification of interference images have been performed.11 20 21 However, because these trials were not initiated using the interference camera with a normal incident specular angle, the background iris coloration and the image defect on the important central corneal area have always interfered with the quantification itself.
To obtain clear interference images, we used the DR-1 tear interference camera (Kowa Co., Nagoya, Japan), because it is a sophisticated system that clearly presents the image while eliminating background iris color and central image defects, while maintaining a normal incidence specular angle, leading to a clear image.19 22 23 24 25 26 27 28 29 30 Thus far, the correlation of the image from the DR-1 camera to lipid film thickness has been reported only rarely.29 30 Furthermore, the colorimetric approach has not been applied to lipid layer interference colors.
Considering the physics principle of a white light source (broadband) thin film interference with colorimetry and recent advances in computer technology, it should be possible to achieve the quantification of the tear lipid interference image to obtain exact lipid film thickness information.11 21 31 32 Before the conversion of color information to thickness information, a primary requirement would be to convert thickness information to color information.
In this study, we produced a computer-synthesized color chart of a human tear lipid interference image by using the DR-1 camera and the colorimetric approach for the conversion from tear lipid film thickness into interference color information.
 |
Methods
|
|---|
Synthesis of the Color Chart
Color consists of three primary hues: red (R), green (G), and blue (B). To synthesize interference color at a certain film thickness (Fig. 1 , d, in nanometers),
-corrected RGB values (R'G'B', equation 1 ) of the colorthat is, corrected for human retinal response to colorswere required, because the DR-1 interference images were output into general National Television Standards Committee (NTSC) signals.19 34 35
 | (1) |

View larger version (12K):
[in this window]
[in a new window]
|
FIGURE 1. Interference caused by reflection from the two surfaces, lipid and aqueous layers: r1 and r2, the reflections from the lipid and aqueous layers; , the amplitude reflectance of interference from these two reflected waves, d, lipid film thickness; , angle of refraction in the layer. In the DR-1 assembly, approaches 0. n, refractive indices of the air, lipid layer, and aqueous layer. (This is an illustration of a simplified model. The actual tear lipid layer causes interference by multiple reflections.10 33 ) Assumptions: DR-1 was originally assembled so that all light from the DR-1 light source through the convex lens would be reflected at the surface and at the back of the tear lipid layer with the normal incidence specular angle within a range of the 8-mm diameter of the cornea.19 Because individual corneal curvatures are variable, it would not always be possible to achieve a normal incidence. However, under the practical usage of the DR-1, we think that these problems would be negligible.
|
|
To obtain RGB values, the following matrix was used to transform X, Y, and Z tristimulus values, which are sets of three linear light components that conform to the Commission Internationale de lÉclairage (CIE) color-matching functions32 34 36 37 38 39 :
 | (2) |
CIE X, Y, and Z tristimulus values of the color of the reflection were obtained as follows32 36 :
 | (3) |
where S(
) (in microwatts/square centimeter) is the spectrum from the DR-1 light source at the wavelength
(in nanometers). The DR-1 camera was assembled with a halogen bulb of 3000 K of color temperature through a heat-absorbing filter.19 Each
,
, and
(no unit) is the color matching function of XYZ standard colorimetric system of CIE 1931.32 36 37 R(
) (no unit) is the energy light reflectance at a light wavelength
.33 36 R(
) was obtained as follows (Fig. 1) .
The amplitude reflectance of interference,
(no unit), in multiple reflection at a thin film monolayer is expressed as (i: an imaginary number)
 | (4) |
where r1 and r2 are Fresnels indices of reflection at each lipid layer surface and aqueous layer interface (Fig. 1) and are given as follows by Fresnels equation, when the angle of specular reflection is normal incidence (DR-1 assembly)11 19 33 40 :
 | (5) |
Also, n0, n1, and n2 are the refractive indices of air, lipid layer, and aqueous layer and are reported as n0 = 1, n1 = 1.48,41 42 and n2 = 1.33.12
Noting that energy is proportional to the square of amplitude, here R(
) is expressed as11 33 43
 | (6) |
where
* is the conjugate complex numbers of
.
Thus, using Eulers equation,
 | (7) |
A phase difference of two waves, r1 and r2that is, 2
1is expressed as:
 | (8) |
where
1 is the angle of refraction (Fig. 1) and is assumed to be 0 in the DR-1 optical system. Thus, cos
1 = 1.11 19 33 40
Thus
 | (9) |
At this point, R(
) can be obtained only from the actual numbers, n0, n1, and n2, d, and
, and could be calculated to obtain the X, Y, and Z tristimulus values in the broadband light source for RGB and R'G'B' values, which were transformed into the synthesized interference color chart. Computer programs (Excel X; Microsoft, Redmond, WA, and ImageJ; image processing and analysis software; http://rsb.info.nih.gov/nih-image, developed by Wayne Rasband of the National Institutes of Health, Bethesda, MD) were used to obtain RGB profiles of the interference colors and to synthesize the interference color chart.
Actual DR-1 images were also obtained from a healthy control subject with normal tear function and a patient with dry eye who had aqueous tear deficiency and meibomian gland obstruction. Interference image data from random pixel points were then converted to lipid layer thicknesses data using the principles of the new color chart.
 |
Results
|
|---|
A synthesized interference color chart for the DR-1 camera (Fig. 2C) was obtained along with
-corrected RGB (R'G'B') profiles (Fig. 2A) and X, Y, and Z tristimulus values (Fig. 2B) , for the range of tear lipid film thickness from 0 to 1000 nm (0 to approximately the 5th interference order).

View larger version (28K):
[in this window]
[in a new window]
|
FIGURE 2. Synthesized interference color chart with -corrected RGB (R'G'B') profile and X, Y, and Z tristimulus values. The interference chart was obtained as shown in (C). From 0 to 1000 nm film thickness, the color forms a line of dark gray-brown, bright gray-brown, almost white, and brown in the 1st interference order. After that, a cyclical pattern of blue, green-yellow, and red is seen, with the change of color intensity according to the increasing lipid layer thickness until approximately the 5th interference order. The -corrected RGB (R'G'B') values (A, red line, R'; green line, G'; and blue line, B'); and X, Y, and Z tristimulus values (B, cyan line X; magenta, Y; and yellow, Z) are also shown. Each RGB value synthesizes the interference color (C). (Please note that K was overridden for photogenic demography, because the calculated color chart was dark. Please see RGB values for the actual results. Please also note that the color chart presented at was printed using CMYK graphics and may differ from the original color using the RGB system. The intensity of the DR-1 light source was relatively low (maximum level with minimal saturated white-out for image analysis), and the chart in general was colored redder than the chart by King-Smith et al.,11 who used a simulated equal-energy spectrum.)
|
|
For tear films thinner than 92.5 nm (interference order lower than 0.5), there was a gray-brown monochromatic increase in intensity. Then, for the film thickness of approximately 185 nm (interference order 0.51), a brown color appeared. After that, a cyclical pattern of blue, green-yellow, and red was observed with the change of color intensity according to increasing lipid layer thickness.
Figure 3 and Table 1 show the result of the conversion of interference colors on the DR-1 image to lipid layer thickness. The central, top, right, bottom, and bottom left points on the cornea were selected, and the corresponding lipid layer thicknesses were displayed as shown in the figure.

View larger version (56K):
[in this window]
[in a new window]
|
FIGURE 3. Clinical application of our method to an actually acquired DR-1 image mapped to show the distribution of the tear lipid layer. The interference color of the DR-1 image was converted to lipid layer thickness data as described in Table 1 . The numbers indicate the tear lipid layer thickness in nanometers. (A) DR-1 image from a 31-year-old normal male subject. Schirmer test, 15 mm; fluorescein score, 0; rose bengal score, 0; tear break up time (BUT), 7 seconds; and meibomian gland orifice obstruction score, 0 points. (B) DR-1 image from a 73-year-old woman with non-Sjögren dry eye. Schirmer test 3 mm. Fluorescein score: 1 point; rose bengal score, 2 points; BUT, 1 second; and meibomian gland orifice obstruction score: 3 points. DR-1 image consists of 640 x 480 pixels. Numbers on the horizontal and vertical axes of (A) and (B) indicate pixel position.
|
|
 |
Discussion
|
|---|
In this study, the tear lipid interference color chart for the specific interference camera, DR-1, which visualizes clear interference images, was synthesized in a colorimetric approach. The colors in the chart were apparently similar to the reported interference colors using the DR-1.19 30 This chart is expected to contribute to the exact quantification of the interference images.
The principles of tear interference colors including the meaning and specifics of each color have not been discussed satisfactorily until now. King-Smith et al.11 first showed the principle of tear interference in detail using a simulated equal-energy spectrum or the spectrum of the Tear Scope (Keeler, Windsor, UK), with fundamental retinal spectral sensitivity. However, usage of the Tear Scope was impractical, because the image always overlapped with the background iris color, producing an unclear image, because it could not achieve the normal incident specular reflection.
Goto and Tseng29 30 reported the trial quantification of the lipid layer thickness with the color comparison method using the DR-1 and the color chart of King-Smith et al.11 based on a simulated equal-energy spectrum, which was not based on the spectrum from the DR-1. However, progress and changes in physical optics and colorimetry stimulated us to develop the current color chart system, which would be more suitable for DR-1 camera optics. The DR-1 camera, having superior specular angle capabilities, was the only camera system appropriate for the quantification of the interference images through its sophisticated optics.19
We believe that our efforts in devising the current chart are the first steps in accurate quantification of the interference images. Our chart differs from that of King-Smith et al.11 in color intensity and interference orders, which makes it more suitable for the DR-1 system. Quantification of RGB color intensities at specific spots on the precorneal tear film may pave the way for the development of the quantification of tear film lipid layer thickness in the near future.
Although tear interferometry has been designed principally for the study of the tear lipid layer, up until now it has been used to investigate the changes in the aqueous layer status, lipidaqueous layer interaction,19 22 and combined aqueous tear and lipid deficiency in certain patients with dry eye.30 Although real-time topographical interference displays in the DR-1 using our color chart have not been achieved at present, conversion of interference colors to lipid layer thickness data based on our logical color chart system has been realized. It is our belief that conducting further clinical studies to determine repeatability of results obtained with our color chart in patients with dry eye and ocular surface disorders would be highly interesting. Our method can also have application in the evaluation of objective and quantitative parameters with different therapeutic modalities.
In conclusion, we developed a new color chart that we believe will be of benefit in converting tear interference color information to tear lipid layer film thickness data.
 |
Acknowledgements
|
|---|
The authors thank Fabrice Manns, PhD, Ocular Biophysics Center, Bascom Palmer Eye Institute, University of Miami (Miami, FL), and Naoshi Shinozaki, Executive Director, Cornea Center and Eye Bank, Tokyo Dental College (Chiba, Japan), for instruction on the principles of the interference phenomena and Hirayuki Sato and Koji Endo, PhD, Analytical Research Center, KAO Corp. (Tochigi, Japan), for providing expertise in linking colorimetry and interferometry.
 |
Footnotes
|
|---|
Supported by grants from the Japanese Ministry of Education, Culture, Sports, Science, and Technology and the Hightech Research Center at Tokyo Dental College.
Submitted for publication March 13, 2003; revised July 3, 2003; accepted July 13, 2003.
Disclosure: E. Goto, None; M. Dogru, None; T. Kojima, None; K. Tsubota, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Eiki Goto, Department of Ophthalmology, Tokyo Dental College, 5-11-13 Sugano, Ichikawa-shi, Chiba 72-8513, Japan; eikigoto{at}nifty.com.
 |
References
|
|---|
- Pflugfelder, SC, Tseng, SCG, Sanabria, O, et al (1998) Evaluation of subjective assessments and objective diagnostic tests for diagnosing tear-film disorders known to cause ocular irritation Cornea 17,38-56[CrossRef][Medline][Order article via Infotrieve]
- Chen, H-B, Yamabayashi, S, Tanaka, Y, Ohno, S, Tsukahara, S. (1997) Structure and composition of rat precorneal tear film. A study by an in vitro cryofixation Invest Ophthalmol Vis Sci 38,381-387[Abstract/Free Full Text]
- Wolff, E. (1946) The muco-cutaneous junction of the lid margin and the distribution of the tear fluid Trans Ophthalmol Soc UK 66,291-308
- McDonald, JE. (1968) Surface phenomena of tear films Trans Am Ophthalmol Soc 66,905-939[Medline][Order article via Infotrieve]
- Norn, MS. (1979) Semiquantitative interference study of fatty layer of precorneal film Acta Ophthalmol (Copenh) 57,766-774
- Guillon, JP. (1982) Tear film photography and contact lens wear J Br Contact Lens Assoc 5,84-87[CrossRef]
- Olsen, T. (1985) Reflectometry of the precorneal film Acta Ophthalmol (Copenh) 63,432-438
- Bron, AJ, Tiffany, JM. (1998) The meibomian glands and tear film lipids: structure, function, and control Adv Exp Med Biol 438,281-295[Medline][Order article via Infotrieve]
- Newton, I. (1718) Optiks: A treatise of the Reflections, Refractions, Inflections and Colors of Light 2nd ed. ,168-206 London.
- Guenther, R. (1990) Interference Modern Optics ,87-128 John Wiley & sons New York.
- King-Smith, PE, Fink, BA, Fogt, N. (1999) Three interferometric methods for measuring the thickness of layers of the tear film Optom Vis Sci 76,19-32[CrossRef][Medline][Order article via Infotrieve]
- Tiffany, JM. (1987) The lipid secretion of the meibomian glands Adv Lipid Res 22,1-62[Medline][Order article via Infotrieve]
- Lemp, MA. (1995) Report of the National Eye Institute/Industry workshop on Clinical Trials in Dry Eyes CLAO J 21,221-232[Medline][Order article via Infotrieve]
- Driver, PJ, Lemp, MA. (1996) Meibomian gland dysfunction Surv Ophthalmol 40,343-367[Medline][Order article via Infotrieve]
- Korb, DR, Baron, DF, Herman, JP, et al (1994) Tear film lipid layer thickness as a function of blinking Cornea 13,354-359[CrossRef][Medline][Order article via Infotrieve]
- Korb, DR, Greiner, JV. (1994) Increase in tear film lipid layer thickness following treatment of meibomian gland dysfunction Adv Exp Med Biol 350,293-298[Medline][Order article via Infotrieve]
- Isreb, MA, Greiner, JV, Korb, DR, et al (2003) Correlation of lipid layer thickness measurements with fluorescein tear film break-up time and Schirmers test Eye 17,79-83[CrossRef][Medline][Order article via Infotrieve]
- Danjo, Y, Hamano, T. (1995) Observation of precorneal tear film in patients with Sjögrens syndrome Acta Ophthalmol Scand 73,501-505[Medline][Order article via Infotrieve]
- Yokoi, N, Takehisa, Y, Kinoshita, S. (1996) Correlation of tear lipid layer interference patterns with the diagnosis and severity of dry eye Am J Ophthalmol 122,818-824[Medline][Order article via Infotrieve]
- Doane, MG, Lee, ME. (1998) Tear film interferometry as a diagnostic tool for evaluating normal and dry-eye tear film Adv Exp Med Biol 438,297-303[Medline][Order article via Infotrieve]
- Khamene, A, Negahdaripour, S, Tseng, SC. (2000) A spectral-discrimination method for tear-film lipid-layer thickness estimation from fringe pattern images IEEE Trans Biomed Eng 47,249-258[CrossRef][Medline][Order article via Infotrieve]
- Yokoi, N, Mossa, F, Tiffany, JM, Bron, AJ. (1999) Assessment of meibomian gland function in dry eye using meibometry Arch Ophthalmol 117,723-729[Abstract/Free Full Text]
- Komuro, A, Yokoi, N, Takehisa, Y, Kinoshita, S. (1998) Association of tear lipid layer interference patterns with superficial punctate keratopathy Adv Exp Med Biol 438,315-317[Medline][Order article via Infotrieve]
- Dogru, M, Katakami, C, Miyashita, M, et al (2000) Ocular surface changes after excimer laser phototherapeutic keratectomy Ophthalmology 107,1144-1152[CrossRef][Medline][Order article via Infotrieve]
- Inoue, K, Kato, S, Ohara, C, et al (2001) Ocular and systemic factors relevant to diabetic keratoepitheliopathy Cornea 20,798-801[CrossRef][Medline][Order article via Infotrieve]
- Goto, E, Shimazaki, J, Monden, Y, et al (2002) Low-concentration homogenized castor oil eye drops for noninflamed obstructive meibomian gland dysfunction Ophthalmology 109,2030-2035[CrossRef][Medline][Order article via Infotrieve]
- Kubo, M, Sakuraba, T, Arai, Y, Nakazawa, M. (2001) Tear lipid layer interference changes after dacryocystorhinostomy Jpn J Ophthalmol 45,653-656[CrossRef][Medline][Order article via Infotrieve]
- Tsubota, K, Dogru, M. (2002) Changing perspectives for the treatment of dry eye Contemp Ophthalmol 1,1-8
- Goto, E, Tseng, SC. (2003) Differentiation of lipid tear deficiency dry eye by kinetic analysis of tear interference images Arch Ophthalmol 121,173-180[Abstract/Free Full Text]
- Goto, E, Tseng, SC. (2003) Kinetic analysis of tear interference images in aqueous tear deficiency dry eye before and after punctal occlusion Invest Ophthalmol Vis Sci 44,1897-1905[Abstract/Free Full Text]
- Toyooka, S, Hayasaka, N, Kataguchi, M, Kobayashi, F. (1998) Spectral based estimation of optical path difference of interference colors Opt Rev 5,207-211
- Kato, T. (2001) A method to synthesize interference color chart with personal computer J Geol Soc Jpn 107,64-67
- Yoshida, S, Yajima, H. (1994) Optics of Thin Film. Optical Thin Films and devices (in Japanese) ,14-17 University of Tokyo Press Tokyo.
- Kobayashi, M. (2002) Color reproduction/color management/color appearance (in Japanese) J Color Sci Assoc Jpn 26,18-29
- Poynton, C. (1996) Gamma Collection on the Apple Macintosh: A Technical Introduction to Digital Video ,1-6 Wiley New York.
- Specification of Colours According to the CIE 1931 Standard Colorimetric System and the CIE 1964 Supplementary Standard Colorimetric System (in Japanese). Japanese Industrial Standards Committee 1999,1-28 Japanese Standards Association Tokyo. Publication JIS Z 8701.
- . CIE (1932) Commission Internationale de lEclairage, Proceedings, 1931. 8th Session ,19-29 Cambridge University Press Cambridge.
- . Commission Internationale de lÉclairage (CIE) (1968) Colorimetry 2nd ed. ,74 Publication CIE 15.2-1986. Vienna, Austria.
- . The Color Science Association of Japan (1998) Handbook of Color Science: New Edition (in Japanese) ,1141 University of Tokyo Press Tokyo.
- Jenkins, FA. (1976) Fundamentals of Optics 4th ed. ,286-289 McGraw-Hill New York.
- Mishima, S, Maurice, DM. (1961) The oily layer of the tear film and evaporation from the corneal surface Exp Eye Res 1,39-45[Medline][Order article via Infotrieve]
- Tiffany, JM. (1986) Refractive index of meibomian and other lipids Curr Eye Res 5,887-889[Medline][Order article via Infotrieve]
- Danjo, Y, Nakamura, M, Hamano, T. (1994) Measurement of the precorneal tear film thickness with a non-contact optical interferometry film thickness measurement system Jpn J Ophthalmol 38,260-266
This article has been cited by other articles:

|
 |

|
 |
 
J. J. Nichols, G. L. Mitchell, and P. E. King-Smith
Thinning Rate of the Precorneal and Prelens Tear Films
Invest. Ophthalmol. Vis. Sci.,
July 1, 2005;
46(7):
2353 - 2361.
[Abstract]
[Full Text]
[PDF]
|
 |
|