1
|
Huang X, Hargrave A, Bentley J, Dubra A. Biometry study of foveal isoplanatic patch variation for adaptive optics retinal imaging. BIOMEDICAL OPTICS EXPRESS 2024; 15:5674-5690. [PMID: 39421787 PMCID: PMC11482173 DOI: 10.1364/boe.536645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/15/2024] [Revised: 08/15/2024] [Accepted: 08/17/2024] [Indexed: 10/19/2024]
Abstract
The change in ocular wavefront aberrations with visual angle determines the isoplanatic patch, defined as the largest field of view over which diffraction-limited retinal imaging can be achieved. Here, we study how the isoplanatic patch at the foveal center varies across 32 schematic eyes, each individualized with optical biometry estimates of corneal and crystalline lens surface topography, assuming a homogeneous refractive index for the crystalline lens. The foveal isoplanatic patches were calculated using real ray tracing through 2, 4, 6 and 8 mm pupil diameters for wavelengths of 400-1200 nm, simulating five adaptive optics (AO) strategies. Three of these strategies, used in flood illumination, point-scanning, and line-scanning ophthalmoscopes, apply the same wavefront correction across the entire field of view, resulting in almost identical isoplanatic patches. Two time-division multiplexing (TDM) strategies are proposed to increase the isoplanatic patch of AO scanning ophthalmoscopes through field-varying wavefront correction. Results revealed substantial variation in isoplanatic patch size across eyes (40-500%), indicating that the field of view in AO ophthalmoscopes should be adjusted for each eye. The median isoplanatic patch size decreases with increasing pupil diameter, coarsely following a power law. No statistically significant correlations were found between isoplanatic patch size and axial length. The foveal isoplanatic patch increases linearly with wavelength, primarily due to its wavelength-dependent definition (wavefront root-mean-squared, RMS <λ/14), rather than aberration chromatism. Additionally, ray tracing reveals that in strongly ametropic eyes, induced aberrations can result in wavefront RMS errors as large as λ/3 for an 8-mm pupil, with implications for wavefront sensing, open-loop ophthalmic AO, spectacle prescription and refractive surgery.
Collapse
Affiliation(s)
- Xiaojing Huang
- Institute of Optics, University of Rochester, Rochester, New York 14620, USA
- Byers Eye Institute, Stanford University, Palo Alto, California 94303, USA
| | - Aubrey Hargrave
- Byers Eye Institute, Stanford University, Palo Alto, California 94303, USA
| | - Julie Bentley
- Institute of Optics, University of Rochester, Rochester, New York 14620, USA
| | - Alfredo Dubra
- Byers Eye Institute, Stanford University, Palo Alto, California 94303, USA
| |
Collapse
|
2
|
Mrukwa Kominek E, Sarnat-Kucharczyk M, Patel S. The impact of exposure on the magnitude of astigmatism formed within the precorneal tear film over the central optical zone of the cornea in ocular surface disease. Cont Lens Anterior Eye 2020; 43:261-267. [PMID: 31937468 DOI: 10.1016/j.clae.2020.01.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Revised: 12/30/2019] [Accepted: 01/01/2020] [Indexed: 11/30/2022]
Abstract
PURPOSE To estimate the astigmatic power and axis of the tear film over the central optical zone of the cornea by vector analysis of topographic data, in ocular surface disease (OSD) and controls, during blink suppression. METHODS Video-keratoscopic images were captured on opening the eyes after a single blink 5, 10 & 15 s later during blink suppression in OSD patients (mixed aetiology, group 1 age 20 - 50 years, n = 12, group 2 > 50 years, n = 38) and controls (group 3, n = 19). The SimK and axis values were used to calculate the astigmatism (power and axis) that formed in the precorneal tear film during each period. Data were aggregated into 3 periods; T0-T5 (between 0 & 5 s after the blink), T5-T10 (5 & 10 s later, T10-T15 (between 10 & 15 s later). RESULTS Mean (± SD, 95%CI) astigmatic power (DC) formed in the tear film over each period was respectively : Group 1, -0.81 DC (0.99, -1.44 to -0.17), -2.65 DC(1.36, -3.52 to -1.79), -1.37 DC (2.15, -2.73 to -0.01). Group 2, -0.33 DC (0.38, -0.45 to -0.20), -0.57 DC (0.97, -0.91 to -0.24) -0.96 DC (2.10, -1.68 to -0.24), Group 3, -0.57 DC (0.55 -0.76 to -0.38), -0.56 DC (0.57, -0.76 to -0.37), -0.31 DC (0.44, -0.46 to -0.16). Changes were significant in groups 1 (p = 0.013) & 3 (p = 0.033) but not in 2 (p = 0.078). Intergroup differences were significant at all periods (p < 0.05). Significant correlations were detected following vector analysis, e.g Group2 between the astigmatism formed during T5-T10 (y) and ocular surface astigmatism at 5 s was y = 0.281x - 0.834 (r = 0.328, n = 38, p < 0.05). In all three groups apparent changes in the axes of astigmatism were not significant (p > 0.05). CONCLUSIONS Changes in the precorneal tear film after blinking are predominately astigmatic indicating that changes in the central region of the tear film following the natural blink are quasi-orthogonal.
Collapse
Affiliation(s)
- Ewa Mrukwa Kominek
- Department of Ophthalmology, School of Medicine in Katowice, Medical University of Silesia, Katowice, 40-514, Poland
| | - Monika Sarnat-Kucharczyk
- Department of Ophthalmology, School of Medicine in Katowice, Medical University of Silesia, Katowice, 40-514, Poland
| | - Sudi Patel
- NHS National Services Scotland, Edinburgh, EH12 9EB, UK.
| |
Collapse
|
3
|
Aldaba M, Mira-Agudelo A, Barrera Ramírez JF, García-Guerra CE, Ramo JP. Tear film stability assessment by corneal reflex image degradation. JOURNAL OF THE OPTICAL SOCIETY OF AMERICA. A, OPTICS, IMAGE SCIENCE, AND VISION 2019; 36:B110-B115. [PMID: 31044986 DOI: 10.1364/josaa.36.00b110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2018] [Accepted: 02/26/2019] [Indexed: 06/09/2023]
Abstract
Tear film stability assessment is one of the main tests in dry eye diagnosis. However, to date, no test methodology has been adopted as the gold standard due to several reasons, such as the methods being invasive, subjective, or unfeasible for the clinical environment. In this paper, a method that overcomes the above-mentioned limitations for tear film stability measurements is presented, and is based on the degradation of corneal reflex images caused by breakups. The experimental setup, which is based on recording the corneal reflex image or the first Purkinje image, is described, as well as the method used to determine tear film stability by means of the associated breakup time (BUT) using corneal reflex image degradation. Images obtained through simulations of the experimental setup are also shown. Moreover, BUT measurements performed using both the conventional fluorescein method and the proposed method in nine healthy adults are presented. Both the experimental and simulation images show corneal reflex image degradation due to the appearance of breakups in the tear film, highlighting the potential of the method to assess tear film stability. We have shown that the corneal reflex image degrades when the tear film breaks up and, thus, the proposed method can be used to assess tear film stability.
Collapse
|
4
|
Polans J, Keller B, Carrasco-Zevallos OM, LaRocca F, Cole E, Whitson HE, Lad EM, Farsiu S, Izatt JA. Wide-field retinal optical coherence tomography with wavefront sensorless adaptive optics for enhanced imaging of targeted regions. BIOMEDICAL OPTICS EXPRESS 2017; 8:16-37. [PMID: 28101398 PMCID: PMC5231289 DOI: 10.1364/boe.8.000016] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2016] [Revised: 10/28/2016] [Accepted: 11/19/2016] [Indexed: 05/03/2023]
Abstract
The peripheral retina of the human eye offers a unique opportunity for assessment and monitoring of ocular diseases. We have developed a novel wide-field (>70°) optical coherence tomography system (WF-OCT) equipped with wavefront sensorless adaptive optics (WSAO) for enhancing the visualization of smaller (<25°) targeted regions in the peripheral retina. We iterated the WSAO algorithm at the speed of individual OCT B-scans (~20 ms) by using raw spectral interferograms to calculate the optimization metric. Our WSAO approach with a 3 mm beam diameter permitted primarily low- but also high- order peripheral wavefront correction in less than 10 seconds. In preliminary imaging studies in five normal human subjects, we quantified statistically significant changes with WSAO correction, corresponding to a 10.4% improvement in average pixel brightness (signal) and 7.0% improvement in high frequency content (resolution) when visualizing 1 mm (~3.5°) B-scans of the peripheral (>23°) retina. We demonstrated the ability of our WF-OCT system to acquire non wavefront-corrected wide-field images rapidly, which could then be used to locate regions of interest, zoom into targeted features, and visualize the same region at different time points. A pilot clinical study was conducted on seven healthy volunteers and two subjects with prodromal Alzheimer's disease which illustrated the capability to image Drusen-like pathologies as far as 32.5° from the fovea in un-averaged volume scans. This work suggests that the proposed combination of WF-OCT and WSAO may find applications in the diagnosis and treatment of ocular, and potentially neurodegenerative, diseases of the peripheral retina, including diabetes and Alzheimer's disease.
Collapse
Affiliation(s)
- James Polans
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA
| | - Brenton Keller
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA
| | - Oscar M. Carrasco-Zevallos
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA
| | - Francesco LaRocca
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA
| | - Elijah Cole
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA
| | - Heather E. Whitson
- Department of Ophthalmology, Duke University Medical Center, Durham, North Carolina 27710, USA
- Department of Medicine (Geriatrics)/Duke Aging Center, Duke University Medical Center, Durham, North Carolina 27710, USA
- Geriatrics Research Education & Clinical Center, Durham VA Medical Center, Durham, North Carolina 27705, USA
| | - Eleonora M. Lad
- Department of Ophthalmology, Duke University Medical Center, Durham, North Carolina 27710, USA
| | - Sina Farsiu
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA
- Department of Ophthalmology, Duke University Medical Center, Durham, North Carolina 27710, USA
| | - Joseph A. Izatt
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA
- Department of Ophthalmology, Duke University Medical Center, Durham, North Carolina 27710, USA
| |
Collapse
|
5
|
Dynamics and function of the tear film in relation to the blink cycle. Prog Retin Eye Res 2014; 45:132-64. [PMID: 25479602 DOI: 10.1016/j.preteyeres.2014.11.001] [Citation(s) in RCA: 90] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2014] [Revised: 10/24/2014] [Accepted: 11/04/2014] [Indexed: 01/31/2023]
Abstract
Great strides have recently been made in quantitative measurements of tear film thickness and thinning, mathematical modeling thereof and linking these to sensory perception. This paper summarizes recent progress in these areas and reports on new results. The complete blink cycle is used as a framework that attempts to unify the results that are currently available. Understanding of tear film dynamics is aided by combining information from different imaging methods, including fluorescence, retroillumination and a new high-speed stroboscopic imaging system developed for studying the tear film during the blink cycle. During the downstroke of the blink, lipid is compressed as a thick layer just under the upper lid which is often released as a narrow thick band of lipid at the beginning of the upstroke. "Rippling" of the tear film/air interface due to motion of the tear film over the corneal surface, somewhat like the flow of water in a shallow stream over a rocky streambed, was observed during lid motion and treated theoretically here. New mathematical predictions of tear film osmolarity over the exposed ocular surface and in tear breakup are presented; the latter is closely linked to new in vivo observations. Models include the effects of evaporation, osmotic flow through the cornea and conjunctiva, quenching of fluorescence, tangential flow of aqueous tears and diffusion of tear solutes and fluorescein. These and other combinations of experiment and theory increase our understanding of the fluid dynamics of the tear film and its potential impact on the ocular surface.
Collapse
|
6
|
Jayakumar V, Thapa D, Hutchings N, Lakshminarayanan V. Are the fluctuations in dynamic anterior surface aberrations of the human eye chaotic? OPTICS LETTERS 2013; 38:5208-5211. [PMID: 24322219 DOI: 10.1364/ol.38.005208] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
The purpose of the study is to measure chaos in dynamic anterior surface aberrations and examine how it varies between the eyes of an individual. Noninvasive tear breakup time and dynamic corneal surface aberrations were measured for two open-eye intervals of 15 s. The maximal Lyapunov exponent (MLE) was calculated to test the nature of the fluctuations of the dynamic anterior surface aberrations. The average MLE for total higher-order aberration (HOA) was found to be small (+0.0102±0.0072) μm/s. No significant difference in MLE was found between the eyes for HOA (t-test; p=0.131). Data analysis was carried out for individual Zernike coefficients, including vertical prism as it gives a direct measure of the thickness of the tear film over time. The results show that the amount of chaos was small for each Zernike coefficient and not significantly correlated between the eyes.
Collapse
|
7
|
Takeda M. Fourier fringe analysis and its application to metrology of extreme physical phenomena: a review [Invited]. APPLIED OPTICS 2013; 52:20-29. [PMID: 23292372 DOI: 10.1364/ao.52.000020] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2012] [Accepted: 09/14/2012] [Indexed: 06/01/2023]
Abstract
The paper reviews a technique for fringe analysis referred to as Fourier fringe analysis (FFA) or the Fourier transform method, with a particular focus on its application to metrology of extreme physical phenomena. Examples include the measurement of extremely small magnetic fields with subfluxon sensitivity by electron wave interferometry, subnanometer wavefront evaluation of projection optics for extreme UV lithography, the detection of sub-Ångstrom distortion of a crystal lattice, and the measurement of ultrashort optical pulses in the femotsecond to attosecond range, which show how the advantages of FFA are exploited in these cutting edge applications.
Collapse
Affiliation(s)
- Mitsuo Takeda
- Center for Optical Research and Education (CORE), Utsunomiya University, 7-1-2 Yoto, Utsunomiya, Tochigi 321-8585, Japan. ‐u.ac.jp
| |
Collapse
|
8
|
|
9
|
Abstract
PURPOSE To review the methods for dynamic, non-invasive, and objective assessment of tear film surface quality and to outline their current state-of-the-art and their future potential. METHODS Among the methods available, high-speed videokeratoscopy, lateral shearing interferometry, and dynamic wavefront sensing are being considered. RESULTS The principles of operations, their advantages and disadvantages, and limitations of each method are being outlined. The possible future directions of each method are also discussed. CONCLUSIONS To gain a better understanding of tear film, its structure and function, it is essential to combine macroimaging technologies with those focusing on tear film microstructure. In this way, one can envisage a clinical device that could help, in future, early diagnosis of dry eye syndrome and development of better materials for contact lenses and eye lubricants.
Collapse
|
10
|
Dai F, Tang F, Wang X, Feng P, Sasaki O. Use of numerical orthogonal transformation for the Zernike analysis of lateral shearing interferograms. OPTICS EXPRESS 2012; 20:1530-1544. [PMID: 22274496 DOI: 10.1364/oe.20.001530] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
A numerical orthogonal transformation method for reconstructing a wavefront by use of Zernike polynomials in lateral shearing interferometry is proposed. The difference fronts data in two perpendicular directions are fitted to numerical orthonormal polynomials instead of Zernike polynomials, and then the orthonormal coefficients are used to evaluate the Zernike coefficients of the original wavefront by use of a numerical shear matrix. Due to the fact that the dimensions of the shear matrix are finite, the high-order terms of the original wavefront above a certain order have to be neglected. One of advantages of the proposed method is that the impact of the neglected high-order terms on the outcomes of the lower-order terms can be decreased, which leads to a more accurate reconstruction result. Another advantage is that the proposed method can be applied to reconstruct a wavefront on an aperture of arbitrary shape from its difference fronts. Theoretical analysis and numerical simulations shows that the proposed method is correct and its reconstruction error is obviously smaller than that of Rimmer-Wyant method.
Collapse
Affiliation(s)
- Fengzhao Dai
- Laboratory of Information Optics and Opto-Electronic Technology, Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai, China
| | | | | | | | | |
Collapse
|
11
|
Abstract
PURPOSE To date, there have been no measuring techniques available that could clearly identify all phases of tear film surface kinetics in one interblink interval. METHODS Using a series of cases, we show that lateral shearing interferometry equipped with a set of robust parameter estimation techniques is able to characterize up to five different phases of tear film surface kinetics that include: (i) initial fast tear film build-up phase, (ii) further slower tear film build-up phase, (iii) tear film stability, (iv) tear film thinning, and (v), after a detected break-up, subsequent tear film deterioration. RESULTS Several representative examples are given for estimating tear film surface kinetics in measurements in which the subjects were asked to blink and keep their eyes open as long as they could. CONCLUSIONS Lateral shearing interferometry is a noninvasive technique that provides means for temporal characterization of tear film surface kinetics and the opportunity for the analysis of the two-step tear film build-up process.
Collapse
|
12
|
Szczesna DH, Kulas Z, Kasprzak HT, Stenevi U. Examination of tear film smoothness on corneae after refractive surgeries using a noninvasive interferometric method. JOURNAL OF BIOMEDICAL OPTICS 2009; 14:064029. [PMID: 20059267 DOI: 10.1117/1.3275850] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
A lateral shearing interferometer was used to examine the smoothness of the tear film. The information about the distribution and stability of the precorneal tear film is carried out by the wavefront reflected from the surface of tears and coded in interference fringes. Smooth and regular fringes indicate a smooth tear film surface. On corneae after laser in situ keratomileusis (LASIK) or radial keratotomy (RK) surgery, the interference fringes are seldom regular. The fringes are bent on bright lines, which are interpreted as tear film breakups. The high-intensity pattern seems to appear in similar location on the corneal surface after refractive surgery. Our purpose was to extract information about the pattern existing under the interference fringes and calculate its shape reproducibility over time and following eye blinks. A low-pass filter was applied and correlation coefficient was calculated to compare a selected fragment of the template image to each of the following frames in the recorded sequence. High values of the correlation coefficient suggest that irregularities of the corneal epithelium might influence tear film instability and that tear film breakup may be associated with local irregularities of the corneal topography created after the LASIK and RK surgeries.
Collapse
Affiliation(s)
- Dorota H Szczesna
- Wroclaw University of Technology, Institute of Physics, Wybrzeze Wyspianskiego 27, Wroclaw 50-370, Poland.
| | | | | | | |
Collapse
|
13
|
Szczesna DH, Kasprzak HT. Numerical analysis of interferograms for evaluation of tear film build-up time. Ophthalmic Physiol Opt 2009; 29:211-8. [DOI: 10.1111/j.1475-1313.2009.00651.x] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
14
|
Krishnaswami K, Bernacki BE, Hô N, Allen PJ, Anheier NC. Lateral shearing interferometer for measuring photoinduced refractive index change in As(2)S(3). THE REVIEW OF SCIENTIFIC INSTRUMENTS 2008; 79:095101. [PMID: 19044448 DOI: 10.1063/1.2973640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
We have built and demonstrated a lateral shearing interferometer as a process engineering and control tool for the fabrication and characterization of direct-laser-written waveguide structures in chalcogenide glasses. Photoinduced change in refractive index of 0.154+/-0.002 was measured for as-deposited amorphous As(2)S(3) thin films at 633 nm with an estimated measurement uncertainty of 1.3% for this air-gap interferometer configuration. The simple design of this interferometer can easily be adapted to other wavelengths including mid- and long-wave infrared regions to measure changes in refractive index or material inhomogeneities in transmissive materials.
Collapse
Affiliation(s)
- Kannan Krishnaswami
- Pacific Northwest National Laboratory, 902 Battelle Boulevard, Richland, Washington 99352, USA
| | | | | | | | | |
Collapse
|
15
|
Li KY, Yoon G. Changes in aberrations and retinal image quality due to tear film dynamics. OPTICS EXPRESS 2006; 14:12552-9. [PMID: 19529690 DOI: 10.1364/oe.14.012552] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
A reliable and objective method to measure aberration changes due to the tear film is essential in improving clinical assessment of the tear film and in vivo retinal imaging. The tear film of 11 subjects are studied by acquiring continuous wavefront measurements in real-time with a customized Shack-Hartmann wavefront sensor. The device has a high resolution lenslet array (190 mum) and a topographer unit with an infrared pupil illuminator (940 nm). A Fourier transform reconstructor algorithm [1] was used to estimate the eyes' wavefront aberrations from slope measurements. Increasing irregularities in the tear film produced observable wavefront variations. The temporal behavior of tear induced aberrations and retinal image quality was evaluated by the root mean squared (RMS) error of the residual wavefront and the volume modulation transfer function (MTF). Similar trends were observed from both metrics. Our analysis demonstrates the applicability of the SH wavefront sensor to assessing the dynamics of the human tear film.
Collapse
|
16
|
Gruppetta S, Lacombe F, Puget P. Study of the dynamic aberrations of the human tear film. OPTICS EXPRESS 2005; 13:7631-6. [PMID: 19498790 DOI: 10.1364/opex.13.007631] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
The dynamic aberrations introduced by the human tear film are studied by measuring the topography of the tear film surface on 14 subjects using a curvature sensing setup. The RMS wavefront error variation of the data obtained is presented showing the non-negligible contribution of the tear film to overall eye aberrations. The tear film wavefronts are decomposed in their constituent Zernike terms, showing stronger contributions from 4th order terms and terms with vertical symmetry, and the temporal behaviour of these aberrations is analysed.
Collapse
|