Ординатура / Офтальмология / Английские материалы / Handbook of Optical Coherence Tomography_Bouma, Tearney_2002
.pdfPolarization-Sensitive OCT |
261 |
birefringence was completely gone, due to excessive thermal damage. Since Maitland et al. [24] showed that the reduction of birefringence in rat tail tendon is a function of both temperature and time, no quantitative conclusions can be drawn on the reduction of birefringence as a function of temperature alone. However, the presented images demonstrate the capability of PS-OCT to measure the birefringence reduction in collagen due to thermal damage.
9.4IN VIVO BURN DEPTH IMAGING
In this section, PS-OCT cross-sectional images obtained in normal and thermally injured rat skin will be presented alongside histology of the imaged region. A qualitative relationship between the loss of birefrigence and the depth of burn injuries determined by histology will be demonstrated. This research was conducted by MD/ Ph.D. student S. Srinivas and Ph.D. student B. H. Park [29].
9.4.1Introduction
Traditionally, thermal burns have been differentiated into first, second, or third degree injuries. This classification is in many respects retrospective. First degree burns show thermal damage limited to the epidermis and are not considered a clinical problem, because skin will regenerate damaged epithelial cells. A second degree or partial thickness burn shows thermal damage extending into the dermis; however, the skin can still heal by re-epithelialization. A third degree burn shows full thickness damage of the epidermis and dermis and fails to heal by regeneration of epithelium from within the wound margins.
The treatment for full thickness burns is skin grafting, either autologous or transplanted. However, a partial thickness burn, which shows destruction of the epidermis and a portion of the underlying dermis, may require a more complex treatment plan depending on the depth of dermal injury. Currently, if a patient has a suspected superficial partial thickness burn, the surgeon will often wait 2–3 weeks to determine if the wound will heal spontaneously from surviving epithelial appendages. If the burn has not healed within this time period, skin grafting is indicated. Conversely, if a deep partial thickness burn is encountered, then skin grafting should be considered as soon as possible, because there is a lower incidence of infection with early eschar removal [30].
A number of methods have been developed to determine the depth of burn injury, including the use of indocyanine green dye fluorescence [31], vital dyes, fluorescein fluorometry, laser Doppler flowmetry (LDF), thermography, ultrasound, nuclear magnetic resonance imaging, and spectral analysis of light reflectance [32]. Clinical estimation by visual and tactile assessment of the wound remains the gold standard for burn depth determination [32]. However, determining whether a deep burn will heal spontaneously is difficult even for an experienced clinician [31,33].
Skin has an abundance of weekly birefringent collagen molecules in the dermis, and ensuing thermal injury will denature this collagen, resulting in a reduction of birefringence. By comparing changes in normal and burned tissues of SpragueDawley rats, we will show a qualitative link between the depth of thermal damage and the measurable loss of skin birefringence.
262 |
de Boer et al. |
9.4.2Experiment
Figures 14–17 show typical PS-OCT scans and corresponding histology for normal skin and for burned rat skin produced by a preheated 75 C brass rod for 5, 15, or 30 s exposure time, respectively. Details of the experiment are described in Srinivas et al. [29]. The PS-OCT images were formed by gray-scale coding the single-pass bire- fringence-induced phase retardation from 0 (black) to 180 (white). A white contour line indicating a 90 phase retardation level—90 in the figure labels—in the PS-OCT image was plotted and overlaid with the original image. Scans were always 4 mm long by 1 mm deep cross sections. The start and end points of the scan were marked with India ink by tattooing the skin to identify the scanned region in histology. To determine the relationship between the loss of birefringence in skin and actual depth of thermal damage, the PS-OCT scans are shown side by side with the histological sections that were taken from the biopsied tissue.
9.4.3Results
Figure 14a shows a control PS-OCT scan, and Fig. 14b displays the abundance of intact hair follicles found in normal rat skin histology. Figure 15a shows an increased depth of the white 90 phase retardation line, and Fig. 15b shows a corresponding purple region near the surface, indicating damage to the dermis, which would be classified as a superficial partial thickness or second degree burn. Figure 16a has an even deeper white 90 phase retardation line, and the purple region indicating damage to the dermis has increased as well (Fig. 16b). Hair follicle damage can be seen, and the burn in Fig. 16b would be classified as a deep partial thickness or, again, a second degree burn. Figure 17a shows the deepest white 90 phase retardation line, and the regressive H & E stain shows damage over the entire dermis (Fig. 17b). Extensive coagulation can be seen in the dermis with a complete absence of hair follicles, classifying the damage in Fig. 17b as a full thickness or third degree burn.
Figure 14 Typical PS-OCT scan (a) and corresponding histology (b) of normal rat skin.
(a) The PS-OCT scan has been gray-scale coded so that black represents 0 phase retardation (the incident polarization) and white is 180 phase retardation. The white contour line in the scan demarcates the depth at which 90 phase retardation has been reached with respect to the incident polarization. (b) The stain used is regressive H & E (mag 157:5 ). Note the abundance of intact hair follicles. (Reprinted from Ref. 29.)
Polarization-Sensitive OCT |
263 |
Figure 15 PS-OCT scan (a) and corresponding histology (b) of rat skin burned by contact with a 75 C brass rod for 5 s. (a) The PS-OCT scan has been gray-scale coded so that black represents 0 phase retardation (the incident polarization) and white is 180 phase retardation. Also note the increased depth of the white 90 phase retardation line in the scan. (b) The regressive H & E stain (mag 157:5 ) shows a purple region near the surface, indicating damage to the dermis. This would be classified as a superficial partial thickness burn. (Reprinted from Ref. 29.)
As can be observed in Figs. 14–17, all PS-OCT scans consistently show that the incident circular polarization on the sample (corresponding to 0 of phase retardation and coded as black) is preserved through deeper depths as thermal damage increases. The depths of the 90 phase retardation contour line increase as the duration of the thermal injury exposure increases. This correlation demonstrates that PS-OCT can provide quantitative information for burn depth determination.
Figure 16 PS-OCT scan (a) and corresponding histology (b) of rat skin burned by contact with a 75 C brass rod for 15 s. (a) The PS-OCT scan has been gray-scale coded so that black represent 0 phase retardation (the incident polarization) and white is 180 phase retardation. Note the increased depth of the white 90 phase retardation line in the scan. (b) The stain used is regressive H & E (mag 157:5 ). The purple region indicating damage to the dermis has increased in depth. Hair follicle damage can also be seen. This would be classified as a deep partial thickness burn. (Reprinted from Ref. 29.)
264 |
de Boer et al. |
Figure 17 PS-OCT scan (a) and corresponding histology (b) of rat skin burned by contact with a 75 C brass rod for 30 s. (a) The PS-OCT scan has been gray-scale coded so that black represents 0 phase retardation (the incident polarization) and white is 180 phase retardation. Again note the increased depth of the white 90 phase retardation line in the scan. (b) Here the regressive H & E stain (mag 157:5 ) shows damage over the entire dermis. Extensive coagulation can be seen in the dermis with a complete absence of hair follicles. This would be classified as a full thickness burn. (Reprinted from Ref. 29.)
9.5POLARIZATION-SENSITIVE OPTICAL COHERENCE TOMOGRAPHY OF THE RABBIT EYE
In this section, PS-OCT cross-sectional images obtained in the cornea and retina of enucleated rabbit eyes will be presented. Four polarization images corresponding to the four Stokes parameters of light reflected from any position in the tissue will be used to analyze birefringent structures in the eye.
This research was performed by doctoral student Mathieu Ducros from the Biomedical Engineering Department, University of Texas, Austin, during a visit to the Beckman Laser Institute and Medical Clinic [34].
9.5.1Introduction
In the eye, the corneal stroma and the retinal nerve fiber layer (NFL) are birefringent. The crystalline lens is also most likely birefringent but will not be discussed here. The corneal stroma is the central and thickest layer of the cornea. It consists of a large number of stacked lamellae. Each lamella has the form of a thin (a few micrometers thick) ribbon. Principal constituents are collagen fibrils, glycosaminoglycans, and water. The fibrils are oriented parallel to the surfaces of the lamellae. Due to this particular structure, a lamella can be compared to a linear retarder in which the slow and fast axes lie respectively parallel and perpendicular to the fibril direction. Each lamella induces a small phase retardation in incoming polarized light. PS-OCT can be used to investigate the corneal stroma structure. Furthermore, depth-resolved polarization images in the cornea may be useful to observe in vivo thermal damage that could occur during treatments such as photorefractive keratectomy (PRK) and laser in situ keratomileusis (LASIK).
The NFL is the most anterior layer of the retina and consists of glial cells and axons of ganglion cells. The axons converge radially to the optic nerve, where they
Polarization-Sensitive OCT |
265 |
exit the ocular globe. The NFL can be up to 400 m thick at the edge of the optic nerve [35]. The NFL birefringence is due to the dense parallel arrangement of axon membranes and microtubules inside the axons [6,36]. Glaucoma is the second leading cause of blindness worldwide and affects the NFL. In glaucoma, the retinal ganglion cells are destroyed and the axons in the NFL disappear. This process is irreversible. An instrument providing a precise and reliable measurement of the NFL thickness and structural integrity could greatly improve glaucoma detection and aid in followup examinations.
9.5.2Phase Retardation Determination
Gray-scale images of the Stokes parameters (I; Q; U; and V) of the cornea, crystalline lens and retina are presented later in Figs. 19 and 21, respectively. In the I image, pixel brightness is proportional to the local value of I in decibels. In images of Q, U, and V, pixel brightness varies between white for þ1 and black for 1. For each pixel
the phase parameter was calculated as |
ð Þ |
||||
¼ |
" |
|
V |
# |
|
|
|
p |
|
|
|
|
arctan |
Q2 þ U2 |
|
44 |
|
|
|
|
|
|
|
In the Poincare´sphere representation [13], is the angle between the Stokes vectors of light reflected from the surface of the cornea (Q ¼ 0, U ¼ 0, V ¼ 1) and that reflected from a given tissue position. According to the known structure of the corneal stroma and retinal NFL, these tissues can be viewed as linear retarders with fast and slow optic axes perpendicular to the light propagation direction (the illumination beam being perpendicular to the corneal surface). In these conditions, represents the phase retardation between two orthogonal polarization components of light reflected from a given tissue element relative to the phase retardation between two orthogonal components of light incident on the eye. equals 0 for left circularly polarized light, 90 for linearly polarized light, and 180 for right circularly polarized light. When the reflected light is unpolarized, or when the intensity is lower than the system noise level, varies randomly.
9.5.3Methods
Experiments were performed on the enucleated eye of a New Zealand White rabbit. Images were acquired less than 6 h post mortem. The sample holder is diagrammed in Fig. 18. The entry window of the eye holder consisted of a flat glass slide that was gently pressed against the cornea to reduce the refractive power. For imaging the cornea, the beam focus was placed approximately a few hundred micrometers anterior to the apex of the cornea. For imaging the retina, the beam focus was placed anterior to the vitreous/retina interface on the superior and/or inferior region of the optic nerve head where the NFL is the thickest. A red aiming beam, aligned with the broadband infrared SLD beam, was used to determine the imaging position on the retina. The red light was partially transmitted through the choroid and sclera and was visible on the back side of the eye. A blank ink dot was placed on the back of the eye at the lateral edges of the imaged region to identify the corresponding area for histology. After PS-OCT measurements, the eye was fixed in formalin and histology performed.
266 |
de Boer et al. |
Figure 18 Eye holder. The sample was surrounded by moist gauze. A flat glass slide was placed on the cornea to allow good visualization of the eye fundus. For imaging the cornea, the beam focus was placed a few hundred micrometers above the glass slide at the beginning of the scan. For imaging the retina, the eye was oriented to place the optic nerve in the field of view and the beam focus was positioned above the vitreous/retinal interface. (Reprinted from Ref. 34 with permission of the Institute of Electrical and Electronics Engineers, Inc.)
9.5.4 PS-OCT Imaging in the Cornea
Stokes Parameter Images
Figure 19 presents the Stokes parameter images (from left to right I, Q, U, and V) acquired from the cornea of an enucleated rabbit eye. Each image is 1-mm wide by 3.5-mm deep. The letters to the left of the I image refer to the beginning of the scan (a), the air/glass interface (b), the glass/cornea interface (c), the cornea/aqueous interface (d), the aqueous/crystalline lens interface (e), and the end of the scan (f ). As expected, the glass slide (b–c) is transparent (I minimum) and does not induce any change in the light polarization. Q, U, and V appear gray. Their average value is 0, and variations are due only to noise. The reflection from the front and back surfaces of the glass appear bright in the I image, black in the V image, and gray in the Q and U images, meaning that the reflected light polarization is left circular. The helicity of the incident light is reversed as expected for specular reflection. In the cornea (about 0.6 mm thick from c to d), I is relatively low and uniform. No anatomical layers can be clearly differentiated. In contrast, Q, U, and V images of the cornea are nonuniform: the polarization state of reflected light varies as a function of depth and lateral position. Variations of Q, U, and V values in depth in the corneal stroma are due to the cumulative polarization effect of all lamellae, with each lamella acting as a linear retarder. A prominent feature indicating an increase in the V value (from black to white) is observed in the right half of the cornea on the V image, meaning that polarization of reflected light changes from left circular to approximately right circular. Therefore, in this region, the whole corneal stroma acts as a quarter-wave retarder in single pass (half-wave retarder in double pass). A different feature is observed in the left half of the V image. There, backscattered light polarization state changes from left circular at the glass/cornea interface (point c) to approximately linear at
Polarization-Sensitive OCT |
267 |
Figure 19 From left to right: Stokes parameters I and normalized parameters Q, U, and V in the rabbit cornea. Each image is 1 mm wide by 3.5 mm deep. The pixel size is 10 m 10 m. The letters a and f correspond to the beginning and end of the scan, and b, c, d, and e indicate the air/glass, glass/cornea, cornea/aqueous, and aqueous/crystalline lens interfaces. Five zones can be seen from top to bottom: air (a–b), glass slide (b–c), cornea (c–d), anterior chamber (d–e), and part of the crystalline lens (e–f ). In the cornea, the I image is relatively uniform, but the Q, U, and V images show regions where the reflected light has different polarization states indicating structural differences in the corneal stroma. (Reprinted from Ref. 34 with permission of the Institute of Electrical and Electronics Engineers, Inc.)
45 (U 1) at the cornea/aqueous interface (point d). Therefore, the corneal birefringence varies between the left and right sides of the imaged region. In the cornea, we also observe that the Q parameter value decreases (from gray to dark pixels) and increases (from gray to bright pixels) at different lateral positions. This can be explained only by variations in the local optic axis direction as a function of lateral position. Indeed, if the optic axis direction was constant over the imaged region, the Q parameter would vary monotonically as a function of depth but could not increase and decrease.
Our results are consistent with the observations of Chang et al. [23], who used small-angle light scattering. They found that the local birefringence and optic axis of the rabbit cornea depend on lateral position. The glass slide modified the anatomical structure of the cornea and may have induced stress birefringence. However, we expect that the effect was small, inasmuch as we observed similar local variations in the polarization state of backscattered light from corneas that were unstressed. In the Q, U, and V images we see that the polarization state of light is not modified by the anterior chamber fluid and remains approximately constant in the anterior lens. Therefore the aqueous fluid and lens do not appear birefringent.
268 |
de Boer et al. |
A-Scan Averages of the Cornea
The amplitude of the Stokes parameters as a function of depth, i.e., in A-scan mode, has a significant noise component due to speckle effects, detector shot noise, and electronic amplification noise. To decrease this noise, we averaged laterally adjacent pixel values (after reregistering A-scan depth positions to align laterally on the glass slide surface). Figure 20 shows an example of an average of 25 adjacent A-scans for the intensity I and phase retardation calculated according to Eq. (44). The left y axis corresponds to I in decibels (solid line) and the right y axis refers to in degrees (dashed line). The letters above the x axis refer to the depth positions described in Fig. 19. In the air (a–b), glass slide (b–c), and aqueous (d–e), the intensity I is minimum. Peak values of I are observed at the air/glass (b), glass/cornea (c), cornea/aqueous (d), and aqueous/crystalline lens (e) interfaces. In the cornea (c–d), I varies between 70 dB and 80 dB, and in the anterior part of the crystalline lens it decreases slowly down to the noise level ( 92 dB). In air (a–b), glass slide (b–c), and aqueous (d–e), the reflected light intensity is below the noise level of the PS-OCT system. In these regions, the value of is not related to the polarization state of the reflected light but varies randomly between 0 and 180 on individual A-scans. averages 65 in Fig. 20, instead of 90 (all values have equal probability), probably because of an imbalance in the detector noise figures. At the air/glass (b), glass/ cornea (c), and aqueous/crystalline lens (e) interfaces, we observed that when the
Figure 20 Plots of I and parameters versus depth. The letters a, b, c, d, e, and f indicate the depth positions introduced on Fig. 19. Plots were calculated by averaging 25 adjacent A- scans of the I and parameter images. was calculated according to Eq. (44). At all interfaces a peak of intensity (I) can be observed. In the air, glass, and aqueous, no light is reflected; and I is about 92 dB and ¼ 65 . In the cornea, I varies between 70 dB and 80 dB, and in the crystalline lens it slowly drops down to the noise level. At the front and back glass slide surfaces, drops to almost 0 ; in the cornea (c–d) it increases from 0 to 120 . It is still 120 at the aqueous/lens interface, showing that the aqueous did not affect the parameter value. In the anterior part of the crystalline lens it decreases from 120 to 100 . (Reprinted from Ref. 34 with permission of the Institute of Electrical and Electronics Engineers, Inc.)
Polarization-Sensitive OCT |
269 |
intensity increases by 6 dB above the noise level, varied by more than 70% of its initial value. Therefore, we decided to consider that was representative of the polarization state of the reflected light when the intensity I is 6 dB or more above the noise level ( 86 dB in Fig. 20).
In the cornea, increases from 0 to about 120 , indicating that the polarization state of the reflected light varies from left circular to some right elliptical state between the anterior and posterior surfaces of the cornea. The aqueous does not modify the phase retardation, because equals 120 at the posterior cornea and anterior crystalline lens. Within the crystalline lens decreases. The value of computed from pixels 250–350 is unreliable because the intensity is lower than86 dB.
9.5.5PS-OCT Imaging in the Retina
Stokes Parameter Images
The Stokes parameter images I, Q, U, and V acquired in the retina of an enucleated rabbit eye, as well as the phase retardation calculated with Eq. (44), are presented in Fig. 21. To reduce speckle noise, each pixel value in the Q, U, and V images was multiplied by the I image pixel value, a low-pass filter was applied to the resulting image, and each pixel was divided by the local average intensity. Each image represents a 3 mm wide by 2 mm deep section. The left edge corresponds approximately to the center of the optic never. Above the optic nerve the posterior hyaloid phase scatters some light and shows up as bright and dark ‘‘clouds’’ in the I and U images, respectively. The dark and light bands following the contour of the retina in the U and V images, respectively, mean that the light reflected from this region is in an elliptical polarization state between right circular and linear at 45 . In the image we notice a region on top of the retina (right half) and in the optic nerve (left third) that changes approximately from gray to white as a function of depth, corresponding to an increase approximately from gray to white as a function of depth, corresponding to an increase in phase retardation. In the region where the retina curves into the optic nerve head, this increase is not as pronounced. At these positions, the axons are oriented downward and the optic axis of the tissue is nearly parallel to the incident light propagation direction.
A-Scan Averages of the Retina
Plots of I and as a function of depth in the retina were obtained by averaging 50 laterally adjacent A-scans. An example is presented in Fig. 22. The intensity increases from 104 dB to 81 dB at the vitreous/retina interface, is maintained at 81 dB over approximately 100 m, then decreases exponentially to 102 dB in 1 mm. Between pixels 31 and 90 the intensity is greater than 6 dB above the noise level ( 102 dB). In this region, describes reliably the polarization state of the reflected light. At the vitreous/retina interface the intensity I increase sharply anddrops to about 60 . Although noise is observed in the plot, the increase after the vitreous/retina interface is significantly higher than the noise amplitude. The distance between the minimum and maximum of the phase retardation was measured on a smoothed plot (a low-pass or median filter can be used). For example, in Fig. 22, increased over a depth of 90 m (pixels 35–44). Since the NFL is
Figure 21 Stokes parameters I, normalized parameters Q, U, and V, and phase retardationimages in the rabbit retina. The left edge of each image is at approximately the center of the optic nerve. A significant change in phase retardation can be observed in the image at the position of the NFL. (Reprinted from Ref. 34 with permission of the Institute of Electrical and Electronics Engineers, Inc.)
Figure 22 The values of I in decibels (solid line) and in degrees (dashed line) are plotted as a function of depth in the rabbit retina. Each plot is an average of 50 depth scans, corresponding to a lateral scanning distance of 500 m. The intensity I, i.e., the conventional OCT signal, increases from 104 dB to 81 dB at the vitreous/retina interface. The retardation , calculated with Eq. (44), increases from about 60 to 115 over a depth of 90 m. (Reprinted from Ref. 34 with permission of the Institute of Electrical and Electronics Engineers, Inc.)
270
