A clinical utility assessment of the automatic measurement method of the quality of Meibomian glands

Meibomian gland dysfunction (MGD) is one of the most common diseases observed in clinics and is the leading cause of evaporative dry eye. Today, diagnostics of MGD is not fully automatic yet and is based on a qualitative assessment made by an ophthalmologist. Therefore, an automatic analysis method was developed to assess MGD quantiatively. The analysis made use of 228 images of 57 patients recorded by OCULUS Keratograph® 5 M with a resolution of 1024 × 1360 pixels concern 30 eyes of healthy individuals (14 women and 16 men) and 27 eyes of sick patients (10 women and 17 men). The diagnosis of dry eye was made according to the consensus of DED in China (2013). The presented method of analysis is a new, developed method enabling an automatic, reproducible and quantitative assessment of Meibomian glands. The analysis relates to employing the methods of analysis and image processing. The analysis was conducted in the Matlab environment Version 7.11.0.584, R2010b, Java VM Version: Java 1.6.0_17-b04 with Sun Microsystems Inc. with toolboxes: Statistical, Signal Processing and Image Processing. The presented, new method of analysis of Meibomian glands is fully automatic, does not require operator’s intervention, allows obtaining reproducible results and enables a quantitative assessment of Meibomian glands. Compared to the other known methods, particularly with the method described in literature it allows obtaining better sensitivity (98%) and specificity (100%) results by 2%.


BioMedical Engineering OnLine
pathological basis is meibomian gland terminal duct obstruction and (or) qualitative or quantitative change of meibum secretion [2]. It can cause abnormal tear film, eye irritation and ocular surface inflammation. More seriously, MGD can cause the ocular surface damage of visual function. At present, the prevalence of MGD in association with dry eye varies greatly, ranging from 3.5% in the Salisbury Eye Evaluation Study [3] to 60.8% in the Shihpai Eye Study in Taiwan [4], 68.3% in the Beijing Eye Study in China [5], and 56.3% in the Singapore Malay Eye Study [6]. MGD was highly in Asian population and associated with various systemic and ocular conditions, such as postmenopausal women, pinguecula, high diastolic blood pressure, anterior blepharitis, corneal contact lens and Demodex mites [6].
As we all know, the tear film lipid layer (TFLL) forms the outermost layer of the film and is composed of both polar and nonpolar lipids. It has been accepted that the TFLL is an indispensable component of the tear film, providing a smooth optical surface for the cornea and retarding evaporation from the eye [7]. Although the precise aetiology of MGD remain to be undetermined, as reported, MGD is a chronic, diffuse abnormality of the Meibomian glands, commonly characterised by terminal duct obstruction and (or qualitative) quantitative changes in the glandular secretion. It may result in alteration of the tear film, symptoms of eye irritation, clinically apparent inflammation, and ocular surface disease [8]. In a recent review, five separate pathophysiological mechanisms of MGD resulting from tear film instability were reported, eyelid inflammation, conjunctival inflammation, corneal damage, microbiological changes and dry eye disease (DED) [9]. The underlying pathophysiological mechanisms of DED and MGD interact, resulting in a double vicious circle. This is consistent with the high incidence of advanced stage MGD.
Dry eye disease is a common condition that occurs when the eyes don't make enough tears, or the tears evaporate too quickly. It is now generally accepted that there are two major classifications of DED: aqueous deficient and evaporative dry eye [10][11][12]. Rather, DED a term that describes a range of conditions with multiple possible etiologies and comorbidities. Additionally, MGD is one of the most common diseases observed in clinics and is the leading cause of evaporative dry eye [13]. MGD occurs when the Meibomian glands, located in the eyelids, do not sufficiently produce and release the oils needed to protect and maintain a healthy tear film. This causes the watery layer in the tear film to evaporate. Thus, the problem for many dry eye patients is not inadequate tear production, as thought for so many years, but a lack of oil production that ensures the protective integrity of the tear film is maintained on a daily basis. Given the increased recognition of the importance of MGD, a great amount of attention has been paid to diagnosis and therapies targeting this condition.

Subject recruitment
The analysis made use of 228 images of 57 patients recorded by OCULUS Keratograph ® 5 M with a resolution of M × N=1024 × 1360 pixels (where M-the number of rows, N-the number of columns). The obtained images L RGB (m,n,k) (where m-row, ncolumn, k-color) concern 30 eyes of healthy individuals (14 women and 16 men) and 27 eyes of sick patients (10 women and 17 men). The diagnosis of dry eye was made according to the consensus of DED in China (2013): (1) at least 1 of 6 symptoms: dryness, burning, sandiness, tiredness, discomfort, and blurred vision with fluorescein tear breakup time ≤5 s or a non-anesthesia Schirmer I test value ≤5 mm/5 min; (2) at least 1 of 6 symptoms: dryness, burning, sandiness, tiredness, discomfort, and blurred vision with 5 s< fluorescein tear breakup time ≤10 s or 5 mm/5 min< non-anesthesia Schirmer I test ≤10 mm/5 min, accompanied by corneal fluorescein staining score. The diagnostic criteria of MGD were as follows: (1) symptoms of ocular discomfort: itching or photophobia with limitations of activities; (2) clinical signs: at least 1 lid margin feature (plugging, vascularity), meibum quality: grade: ≥4, expressibility: ≥1. Potential subjects were excluded from the study if they had undergone previous corneal or ocular surgery, had any ocular pathology other than DED and DI, had systemic diseases known to affect the eye. Participants were instructed to remove soft contact lenses at least 2 weeks and rigid contact lenses at least 1 month before the examination. Data were collected from August 2015 to May 2016 in Beijing Tongren Hospital, Beijing, China. All participants signed an informed consent form in accordance with the tenets of the Declaration of Helsinki and this study was approved by the institutional review board of Beijing Tongren Hospital, Beijing, China. The analysis was conducted in the Matlab environment Version 7.11.0.584, R2010b, Java VM Version: Java 1.6.0_17-b04 with Sun Microsystems Inc. with toolboxes: Statistical, Signal Processing i Image Processing.

Ocular surface disease index
Before the other clinical examinations, patients were required to complete the Ocular Surface Disease Index (OSDI) Survey. This questionnaire gives a range of 0 (no symptoms) to 100 (severe symptoms).

Corneal fluorescein staining
Corneal fluorescein staining was graded from 0 to 12, a sum of the scores of corneal four quadrants, which were scored individually as 0 (no staining), 1 (mild staining with a few scattered dots of stains), 2 (moderate staining between 1 and 3), and 3 (severe staining with confluent stains or corneal filaments) [13].

Schirmer I test
The SIT was a useful assessment of aqueous tear production. The inferior conjunctival fornix was dried with a cotton stick. One minute later, a standard 5 × 40 mm Schirmer test strip was placed over the junction of the middle and outer third of inferior lid. The patients are instructed to keep their eyes closed during the test. The test lasted 5 min, and the amount wetting was recorded [14].

Fluorescein tear breakup time
BUT is the time interval between the last blink and the appearance of the first random dry spot on the corneal surface after the instillation of fluorescein. The average of three consecutive FBUT values was calculated [15].

Meibomian gland function measurement
The following three parameters were the most commonly used methods to evaluate the morphological characteristics and function of Meibomian glands in clinical practice: abnormalities of lid margins, expression of meibum, and gland dropout degree visualized by meibography [16].

Lid margin
Lid margin abnormality was recorded according to the existence of the following four signs: lid margin irregular, vascular engorgement, glandular orifices obstruction, and anterior or posterior displacement of the mucocutaneous junction, the scored was from 0 to 4 [17].

Meibum expression
The quality of the expression of meibum was assessed semiquantitatively in 8 glands of the central third area of the lower eyelids: Grade 0 clear fluid, Grade 1 cloudy fluid, Grade 2 cloudy, particulate fluid, and Grade 3 inspissated, toothpaste-like fluid [18].

Image acquisition
Image acquisition was conducted in Beijing Tongren Hospital in China [36]. Each patient was sitting still in front of Keratograph 5 M. The ophthalmologist was simultaneously performing the examination and recording the images everting the upper and lower lid of the left and right eye of the patient. In total, 228 images of left and right eye of both upper and lower lid were obtained from 57 patients. The recorded sequence of images EX0000*.zip was automatically decompressed to acquire the images of lower lowerLid. BMP and upper lid upperLid.BMP. The names of the files are fixed for each patient. In the next stage of image processing, the algorithm works independently of the lid (lower or upper) and eye (right or left). The image inserted into the workspace (of lower or upper lid) will be further described as L RGB (m,n,k).

Initial image processing
Initial image processing was linked to three stages: filtering with median filter, normalization and removing uneven illumination [37][38][39]. Filtration with median filter was performed with the use of h MED mask of the size M MED × N MED = 3 × 3 pixels [40]. The output image L MED (m,n,i) was further normalized i.e. the image after normalization L NORM (m,n,i) will be equal to: The last stage of the initial image processing is removal of uneven illumination. Illumination unevenness is clearly visible almost in every image [20]. In practice in order to remove it, the deduction of the result of low-pass filtration L GROUND (m,n,i) from the input image, in this case L NORM (m,n,i), is applied, i.e. the image after removal of unevenness: where: ⊖ erosion, SE structural element of the size of M SE × N SE = 31 × 31 pixels. The output image L NLOW (m,n,i) is further subject to major image processing.

Major image processing
Major image processing was performed on the basis of two images: L GRAY (m,n,i) and L NLOW (m,n,i). The second one was destined to detect the location of Meibomian glands which are usually arranged vertically (although there may also be other arrangements) in relation to the image plane (positioning of the patient). The analysis is based on using Riesz pyramid i.e. h R mask defined as: where σ m ,σ n ,σ α standard deviation of the mean for three dimensions m,n,α, α mask rotation angle, m R , n R , α R normalized values m, n,α up to the range from −0.5 to 0.5, i.e. for example for m R : Riesz pyramid enables creating mask sequences which in turn allow separating Meibomian glands located at various angles relative to image plane and privileged vertical position. In Fig. 1 an exemplary pyramid of h R mask defined for h R σ m = 0.05,σ n = 0.05,σ α = 0.4 and second derivative calculated after variable m is presented. The obtained results, other forms of h R mask, are shown in Fig. 2 for other exemplary values σ m ,σ n and σ α . The obtained results, images L RIESZ (m,n,i), resulting from the h R mask convolution for exemplary eye images L NLOW (m,n,i) are shown in Fig. 3.
The next stage is to designate the area covered by the Meibomian glands. In this case the image L RIESZ (m,n,i) will be subject to averaging based on the structural element SE. This operation will encompass the substitution of areas for which the maximum values of matrix L RIESZ (m,n,i) have exceeded the set threshold p r by the value of '1' . Therefore, the resulting matrix L RES (m,n,i) is equal to:   where ⊕ dilation operator.
The structural element SE2 of the size of M S2 × N S2 = 60 × 60 pixels with the center of coordinate system placed in the middle of the first row.
The last operation is to establish the Bézier curve allowing the separation of the area of upper and lower lid-according to Fig. 4. In Fig. 4a a demonstrative image with characteristic automatically recognized areas is presented: 1,2 area and Bézier curve of the lower eyelid boundary; 3,4-are and Bézier curve of upper eyelid boundary; 5-eyelid edge area. Summing up, the new, proposed automatic method of analysis enables automatic calculation of: • Area occupied by Meibomian glands, • Area occupied by the eyelid edge, • Area occupied by the eyelid.
The conducted analysis, particularly the percentage area occupied by Meibomian glands in relation to the whole surface of the eyelid (upper and lower respectively) constitutes the most important result of the proposed new algorithm. In order to determine if the patient is ill or healthy the‚ cut-off ' threshold-p o was set. Above the p o threshold the patients are considered to be healthy, while below the threshold they are considered ill. The result will be compared with the analysis conducted manually by specialized doctors from Beijing Institute of Ophthalmology, Beijing Tongren Eye Center, Beijing  16:82 Tongren Hospital in China. The results of comparison are shown for exemplary patients in Table 1. In Fig. 4b a demonstrative image L NLOW (m,n,i) was shown with marked digits: 6-area of eyelid edge; 7-the results of manual indication of Meibomian glands area determined by an expert together with the area occupied by them calculated automatically 8. The obtained results in Table 1 concern both healthy and ill patients. According to the expert's decision they were divided into two groups:  Table 2. The ROC curve (receiver operating characteristic) as a relationship between sensitivity and specificity for the changed threshold p o in the range from 0 to 100% is shown in Fig. 5.

Discussion
The developed algorithm was compared with the method proposed in [41]. In the method described in [41] the results of sensitivity were obtained at the level of 99.3% and specificity at 97.5% in the diagnostics of Meibomian glands. The algorithm presented in [41] is not insensitive to the change of parameters, while the speed of image analysis for one patient, according to the Authors' declaration, does not exceed 0.5 s. This attitude however has disadvantages which involve the necessity of analyzing the location and   [41]) but also allows obtaining better sensitivity results equaling 98% and specificity equaling 100%. In Fig. 6 exemplary results comparing two methods are presented: (a) images with automatically marked position of Meibomian glands according to the method described in [41]; (b) area covering Meibomian glands according to the method described in [41] and (c) results of the analysis described in this article (indications complying with indications in Fig. 6b. The key difference visible during comparison of the results is the way of analysis of the areas occupied by Meibomian glands. For the method described in [41] it is also the area occupied by the eyelid edge and therefore included in the area of Meibomian glands. For bigger areas of glands the values are of no significance. However, for small areas occupied by Meibomian glands it is of crucial importance and the resulting consequences are reflected in the values of specificity and sensitivity (97.5 and 99.3% respectively).

Conclusions
In conclusion, this noncontact meibography system is a useful and patient-friendly method to obtain information about the meibomian gland structure. Using this meibography system, the authors found that changes in Meibomian glands increase with age, develop earlier in men than in women, and are more severe in elderly men than in elderly women. The authors hope that this meibography method will contribute to the elucidation of the pathologic mechanisms of meibomian gland diseases and will help to establish better diagnostic criteria for MGD. The presented, new method of analysis of Meibomian glands is fully automatic, does not require operator's intervention, allows obtaining reproducible results and enables a quantitative assessment of Meibomian glands. Compared to the other known methods, particularly with the method described in [41] it allows obtaining better sensitivity and specificity results by 2%. Moreover, it additionally enables removing the area occupied by the eyelid edge. It can be concluded that the analysis of particular Meibomian glands' location is not an indispensable element and feature. In each case the analysis of percentage area occupied by Meibomin glands in relation to the whole surface of the eyelid is meaningful and adequate.
Currently, the method is being tested at the Beijing Institute of Ophthalmology, Beijing Tongren Hospital in China.  6 Exemplary results allowing the qualitative comparison of two methods: a image with automatically marked location of Meibomian glands according to the method described in [41]; b area covering Meibomian glands according to the method described in [41] and c results of analysis described in this article (indications complying with indications in b)