Skip to main content

Patient experience and repeatability of measurements made with the Pentacam HR in patients with keratoconus

Abstract

Background

To investigate whether the repeatability of measurements with the Pentacam HR in patients with keratoconus is improved by patients gaining more experience of the measurement situation. Such an improvement could enhance the accuracy with which progressive keratoconus can be detected.

Methods

Four replicate measurements were performed on Day 0 and on Day 3. Parameters commonly used in the diagnosis of progressive keratoconus were included in the analysis, namely the flattest central keratometry value (K1), the steepest central keratometry value (K2), the maximum keratometry value (Kmax), and the parameters A, B and C from the Belin ABCD Progression Display. In addition, quality parameters used by the Pentacam HR to assess the quality of the measurements were included, namely the analysed area (front + back), 3D (front + back), XY, Z, and eye movements.

Results

Neither the diagnostic parameters nor the quality parameters showed any statistically significant improvement on Day 3 compared to Day 0. The quality parameter “eye movements” deteriorated significantly with increasing Kmax.

Conclusion

Gaining experience of the measurement situation did not increase the accuracy of the measurements. Further investigations should be performed to determine whether the increasing number of eye movements with increasing disease severity has a negative effect on the repeatability of the measurements.

Peer Review reports

Introduction

Corneal crosslinking (CXL) is used worldwide to arrest the progression of keratoconus disease [1]. Pre-clinical investigations have shown that CXL stiffens the cornea [2] while clinical investigations have indicated its clinical efficacy in halting keratoconus disease progression [3], also in the long-term perspective [4]. Furthermore, it has been suggested that CXL reduces the need for corneal transplantation [5]. The most common indication for CXL is progressive keratoconus [6], however, children and adolescents are usually referred to CXL upon the diagnosis of keratoconus as the risk of progression is high in younger patients [3, 7]. Progressive keratoconus is commonly assessed by subjective factors such as medical history and visual acuity [8], in addition to objective parameters measured by corneal tomography. The most commonly used corneal tomographer in scientific investigations is the Pentacam HR (Oculus Optikgeräte, GmbH, Wetzlar, Germany) [1]. Less is known about the equipment used in clinical practice, but a recent survey of northern European countries suggested that the Pentacam HR was the most commonly used tomographer [7].

In terms of objective parameters measured by the Pentacam HR, the most frequently used parameter in scientific investigations is the maximum keratometry value (Kmax) [6, 9,10,11], and sometimes also the steepest central keratometry value (K2) [3]. From the clinical perspective, Kmax has been suggested to be the most important parameter, followed by the parameters A, B and C from the Belin ABCD Progression Display in the Pentacam HR [7]. We therefore included these parameters in the analysis, together with the flattest central keratometry value (K1), as K1 and K2 are commonly used for the objective assessment of astigmatism. All the parameters included reflect the opinions expressed regarding the assessment of progressive keratoconus in the Global Consensus on Keratoconus and Ectatic Diseases [12].

It is of the utmost importance to diagnose progressive keratoconus accurately for timely referral to CXL, and for the reliable recruitment of patients with truly progressive keratoconus in clinical trials [13, 14]. As there is no gold standard for assessing progressive keratoconus, diagnosis must be based on a reliable calculation of the measurement accuracy, i.e., the repeatability of the measurements made with the equipment used for examination [15]. Numerous articles have been published on the general repeatability of measurements using the Pentacam HR in patients with keratoconus and healthy controls [16,17,18]. Furthermore, specific factors that can affect the repeatability have been investigated, such as inter-observer effects [19], the effects of disease severity[20], the number of replicate measurements and inter-day effects [13]. However, an aspect that merits further evaluation is the possible effect of the patient gaining experience of the measurement situation. In other areas of ophthalmology, such as perimetry in glaucoma, it has been demonstrated that the results improve as patient experience of the measurement situation increases [21]. In a previous analysis of the inter-day repeatability of measurements using the Pentacam HR in patients with keratoconus, the keratometric parameters were found to show better repeatability on the second measurement occasion (Day 3 compared to Day 0) [13]. It is thus of interest to determine whether the measurement accuracy can be improved by increasing the patient’s experience of the measurement situation.

Subjects and methods

The study was conducted at the Department of Ophthalmology at Skåne University Hospital, Lund, Sweden, according to the tenets of the declaration of Helsinki. All participants were given written information on the study, and written consent was obtained The Swedish Ethical Review Authority in Lund, Sweden, approved the studies (No. 2015/373).

The following diagnostic parameters were included as these are frequently used in scientific and clinical practice [7].

  • The flattest central keratometry value in a 15 degree ring around the apex (K1).

  • The steepest central keratometry value in a 15 degree ring around the apex (K2).

  • The maximum keratometry value (Kmax).

  • The anterior curvature of the 3 mm zone over the thinnest point of the cornea (A), the posterior curvature of the 3 mm zone over the thinnest point of the cornea (B) and the thickness at the thinnest point of the cornea (C), from the Belin ABCD Progression Display (22).

The following quality parameters were included as these could be dependent, or partially dependent, on the patient or the operator.

  • Analysed area (front): This parameter ensures that the measured area is sufficient. The measurement of this parameter can be affected if the patient’s eye is not sufficiently open, by blinking, or when the patient has long eye lashes.

  • Analysed area (back): as above, but for the posterior surface.

  • 3D (front and back): This parameter ensures correct 3D modelling of the anterior and posterior surfaces of the cornea. Excessive blinking, loss of fixation or insufficient opening of the eyes affects this parameter negatively.

  • XY and Z: These parameters describe whether the operator of the Pentacam HR has moved the X, Y, Z base slide at the moment the measurement started. These parameters were included to account for effects of the operator on the measurements.

  • Eye movements: this parameter describes excessive eye movements due to loss of fixation.

The following quality parameters were excluded as these were deemed not to be affected by better patient compliance.

  • Valid data: This parameter ensures that a sufficient number of valid data points are found. Factors that can affect this parameter are not blinking before the measurement and the room illumination. As all the patients were instructed to blink before the measurement, and as the illumination in the room was the same on all measurement occasions, this parameter was excluded.

  • Lost segments and lost segments continuous. These parameters are binary, (0 = Approved, 1 = Not approved).

Participant enrolment [13]

Patients with keratoconus fulfilling the inclusion criteria described below were enrolled consecutively. The inclusion criteria were: keratoconus Stage ≤ 2 [23] with no history of, and no current signs of other ocular pathology, including ocular surface disease and external diseases such as dry eyes and atopy. Only subjects who had not previously undergone ocular surgery and who were aged ≥ 18 years were recruited. Pregnant and breastfeeding women were also excluded. Contact lens wear was discontinued at least 2 weeks before the measurements were made. Patients with Stage 3–4 keratoconus were excluded as the purpose was to study those with less advanced disease. Keratoconus was diagnosed clinically and by examination using the Pentacam HR. The sagittal curvature pattern, posterior and anterior elevation maps, and corneal thickness pattern were assessed, in addition to information from the Belin-Ambrosio Enhanced Ectasia Display.

Twenty-five patients were enrolled. If two eyes were eligible for inclusion, both were examined (see Examination below). Computerised randomisation was performed in patients where both eyes met the inclusion criteria to select one eye for inclusion in the study. Twenty-two participants were male and three female, and the mean age of the group was 27 years (range 21–45 years). Twelve right and thirteen left eyes were included.

Equipment [13]

The Pentacam HR is a Scheimpflug-based tomographic system (Pentacam HR, version 1.20r10, Oculus Optikgeräte GmbH, Wetzlar, Germany). The technical features of this system have been described elsewhere [16]. The default setting of 25 pictures per second was used.

Examination [13]

Four replicate measurements were made on two separate occasions (Day 0 and Day 3) by the same examiner (IG). Subjects were instructed to blink between measurements, but not to lean back. Measurements were made during normal working hours. Only examinations deemed “OK” by the Pentacam HR were accepted. The right eye was examined first, then the left, if both eyes were eligible for inclusion. This reflects normal clinical practice, where both the patient’s eyes are usually examined. When recruitment to the study was complete, computerised randomisation was performed to select one eye per subject.

Statistical methods and calculations [13]

IBM SPSS Statistics 22 for Windows (IBM Corporation, Armonk, NY, USA) and SAS Enterprise Guide 6.1 for Windows (SAS Institute Inc., Cary, NC, USA) were used for statistical analyses. A p-value below 0.05 was considered significant. Descriptive statistics are given as subject mean, standard deviation (SD), and minimum and maximum values. Repeatability was assessed by calculating the within-subject standard deviation (Sw), precision, repeatability coefficient, intra-class correlation (ICC) and coefficient of variation (CV%) with associated confidence intervals (CIs) [15, 24, 25]. Kendall’s Tau-b was used to assess the relationship between the mean and SD, and natural-logarithm-transformed data were analysed when appropriate. Differences between coefficients of variation were assessed using a regression test [26]. The Wilcoxon signed ranks test was used for comparisons of the means of the quality parameters on Day 0 and Day 3. A professional medical statistician was consulted and performed the analysis.

Results

Diagnostic parameters

Descriptive statistics and the repeatability of the measurements of the diagnostic parameters on Day 0 and Day 3 are presented in Table 1. The anterior keratometric parameters (K1, K2, Kmax and A) all showed better repeatability on Day 3. Parameter B showed a slightly poorer repeatability, while that of parameter C was unchanged. Kmax appeared to show the greatest improvement in the repeatability of the measurements, with a 0.18 D improvement on Day 3 (Day 0, 0.70 D and Day 3, 0.52 D); the 95% CIs on Day 0 and Day 3 barely overlapped. However, a regression test comparing the coefficients of variation on Day 0 and Day 3 did not result in any statistically significant difference in Kmax, or for any of the other parameters.

Table 1 Descriptive statistics and repeatability of the measurements on Day 0 and Day 3 in subjects with keratoconus

Quality parameters

Descriptive statistics of the quality parameters and the results of the Wilcoxon signed ranks test between Day 0 and Day 3 are presented in Table 2. All the quality parameters had high standard deviations, and the results of the Wilcoxon signed ranks test must thus be interpreted with caution. Apart from the parameter Z, no statistically significant differences were found between Day 0 and Day 3. The value of Z was significantly lower on Day 0 (177) than on Day 3 (220), suggesting that the operator moved the base slide of the instrument more at the moment when the measurement started. However, the high SD (88 on Day 0 and 100 on Day 3) makes the outcome of the Wilcoxon signed ranks test highly unreliable and difficult to interpret.

Table 2 Values of the quality parameters (mean) with associated standard deviations (SD) on Day 0 and Day 3 together with the results of the paired-samples tests for Day 0 and Day 3

The high variation in the quality parameters is illustrated in the figures given in the Supplementary Information: Z (Fig. S1), XY (SI 2), Eye movements (SI 3), 3D front (SI 4), 3D back (SI 5), Area front (SI 6) and Area back (SI 7).

Increasing keratoconus disease severity can lead to increasing difficulties in focusing during the measurements, thus increasing the number of eye movements. A correlation test was thus performed between Kmax and the eye movement parameter, revealing a highly significant association between Kmax and the number of eye movements (Kendall’s Tau–B, 0.440, p = 0.002) (Fig. 1).

Fig. 1
figure 1

The mean value of Kmax for each patient plotted against the mean value of the number of eye movements

Discussion

The results of this investigation suggest that there is no need to perform measurements on more than one occasion to increase the patient’s experience of the measurement situation for the purpose of increasing the accuracy with which progressive keratoconus can be detected. Although the repeatability of the measurements of keratometric parameters was better on Day 3, no statistically significant differences were found. Furthermore, the quality parameters given by the Pentacam HR, which are used to assess the quality of the measurements, showed no statistically significant improvement on Day 3 compared to Day 0. The variation of the different quality parameters appears to be random, and not associated with a specific replicate in the chain of measurements, as could be suspected. On the one hand, the initial measurements could be associated with poorer quality scores as the patient is less experienced, and then improve with increasing experience. On the other hand, performing repeated measurements could also tax the patient’s attention, leading to reduced compliance during the measurement procedure.

In previous investigations we have shown the importance of using a mean of replicates when assessing progressive keratoconus, as this increases the accuracy of the diagnosis compared to the use of single measurements [13, 27]. Therefore, it is important to understand that this study demonstrates that up to four replicates can be performed on each occasion without affecting the repeatability of the measurements on each occasion.

In other areas of ophthalmology, such as perimetry in glaucoma, it has been shown that increasing experience in testing leads to better test reliability [21]. The perimetric examination is highly dependent on the individual’s responses, and several factors can affect these responses through learning or experiencing the test procedure. One such factor is the attention span needed to maintain a steady gaze and to give a timely response to each light stimulus during the examination, which often lasts from two to several minutes. The reason why measurements with the Pentacam HR do not improve with increasing patient experience could be the short measurement period, which only lasts about 2 s. However, keeping the gaze fixed on the red spot for 2 s may be challenging for patients with keratoconus. The irregular astigmatism deforms the shape of the red fixation spot, which can provoke eye movements in search of an optimal gaze position. In fact, this study revealed a strong statistically significant association between increasing number of eye movements and increasing values of Kmax. It remains to be elucidated whether the association is a contributing or causative factor to the well-known association between deteriorating repeatability of the measurements and increasing keratoconus disease severity, as we have demonstrated in Pentacam HR measurements in both an intra-day [20] and inter-day setting [13].

The strength of this study is that the measurements were performed by the same examiner, thus avoiding possible inter-examiner effects [19]. In addition, measurements were performed during normal working hours, reducing the likelihood of diurnal effects [28,29,30]. Furthermore, the illumination of the room was the same on all measurement occasions. A limitation is that only subjects with less advanced keratoconus were recruited, therefore, the results of this study are only applicable to this group. A further limitation is that the results of the inter-day Wilcoxon signed ranks test of the XY and Z parameters must be interpreted with caution. The high standard deviation of these parameters makes the results unreliable. A prospective evaluation would, however, require a large number of participants to provide adequate power for the analysis of such data. Future investigations should also be carried out to determine the inter-day repeatability when using optical coherence tomography equipment for the diagnosis of progressive keratoconus.

In summary, the results of this study show that the repeatability and the quality of the measurements with the Pentacam HR do not improve with increasing experience of the measurement situation in patients with less advanced keratoconus. However, it was found that the number of eye movements during the measurement procedure increased with increasing keratoconus disease severity. As it is well-known that the repeatability deteriorates with increasing disease severity, it would be of considerable interest to investigate whether the increasing number of eye movements contributes to the deterioration in the repeatability of the measurements. Accounting for this could lead to increased measurement precision in these subjects, which would be of considerable clinical value.

Data Availability

The datasets with all data are available as supplementary information.

Abbreviations

CI:

Confidence Interval

CV%:

Coefficient of variation

CXL:

Corneal Cross-linking

SD:

Standard Deviation

Sw:

Within-subject standard deviation

References

  1. Ng SM, Ren M, Lindsley KB, Hawkins BS, Kuo IC. Transepithelial versus epithelium-off corneal crosslinking for progressive keratoconus. Cochrane Database Syst Rev. 2021;3(3):Cd013512.

    PubMed  Google Scholar 

  2. Spörl E, Huhle M, Kasper M, Seiler T. Increased rigidity of the cornea caused by intrastromal cross-linking. Ophthalmologe. 1997;94(12):902–6.

    Article  PubMed  Google Scholar 

  3. Larkin DFP, Chowdhury K, Burr JM, Raynor M, Edwards M, Tuft SJ, et al. Effect of corneal cross-linking versus Standard Care on Keratoconus Progression in Young Patients: the KERALINK Randomized Controlled Trial. Ophthalmology. 2021;128(11):1516–26.

    Article  PubMed  Google Scholar 

  4. Raiskup F, Herber R, Lenk J, Ramm L, Wittig D, Pillunat LE, et al. Corneal crosslinking with Riboflavin and UVA Light in Progressive Keratoconus: fifteen-year results. Am J Ophthalmol. 2023;250:95–102.

    Article  CAS  PubMed  Google Scholar 

  5. Godefrooij DA, Gans R, Imhof SM, Wisse RP. Nationwide reduction in the number of corneal transplantations for keratoconus following the implementation of cross-linking. Acta Ophthalmol. 2016;94(7):675–8.

    Article  PubMed  Google Scholar 

  6. Ng SM, Hawkins BS, Kuo IC. Transepithelial Versus epithelium-off corneal crosslinking for Progressive Keratoconus: findings from a Cochrane systematic review. Am J Ophthalmol. 2021;229:274–87.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Gustafsson I, Vicente A, Bergström A, Stenevi U, Ivarsen A, Hjortdal J. Current clinical practice in corneal crosslinking for treatment of progressive keratoconus in four nordic countries. Acta Ophthalmol. 2023;101(1):109–16.

    Article  CAS  PubMed  Google Scholar 

  8. Wisse RPL, Simons RWP, van der Vossen MJB, Muijzer MB, Soeters N, Nuijts R, et al. Clinical evaluation and validation of the dutch crosslinking for Keratoconus score. JAMA Ophthalmol. 2019;137(6):610–6.

    Article  PubMed  Google Scholar 

  9. Kobashi H, Tsubota K. Accelerated Versus Standard corneal cross-linking for Progressive Keratoconus: a Meta-analysis of Randomized controlled trials. Cornea. 2020;39(2):172–80.

    Article  PubMed  Google Scholar 

  10. Achiron A, El-Hadad O, Leadbetter D, Hecht I, Hamiel U, Avadhanam V, et al. Progression of Pediatric Keratoconus after corneal cross-linking: a systematic review and pooled analysis. Cornea. 2022;41(7):874–8.

    Article  PubMed  Google Scholar 

  11. Wen D, Li Q, Song B, Tu R, Wang Q, O’Brart DPS, et al. Comparison of Standard Versus accelerated corneal collagen cross-linking for Keratoconus: a Meta-analysis. Invest Ophthalmol Vis Sci. 2018;59(10):3920–31.

    Article  CAS  PubMed  Google Scholar 

  12. Gomes JA, Tan D, Rapuano CJ, Belin MW, Ambrósio R Jr, Guell JL, et al. Global consensus on keratoconus and ectatic diseases. Cornea. 2015;34(4):359–69.

    Article  PubMed  Google Scholar 

  13. Gustafsson I, Bergstrom A, Cardiakides A, Ivarsen A, Hjortdal JO. The inter-day repeatability of parameters for the assessment of progressive disease in subjects with less advanced keratoconus. Am J Ophthalmol. 2021;225:38–46.

    Article  PubMed  Google Scholar 

  14. Wagner H, Barr JT, Zadnik K. Collaborative longitudinal evaluation of Keratoconus (CLEK) study: methods and findings to date. Cont Lens Anterior Eye. 2007;30(4):223–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Bland JM, Altman DG. Measurement error. BMJ. 1996;313(7059):744.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. de Luis Eguileor B, Escudero Argaluza J, Pijoán Zubizarreta JI, Santamaria Carro A, Etxebarria Ecenarro J. Evaluation of the reliability and repeatability of Scheimpflug System Measurement in Keratoconus. Cornea. 2018;37(2):177–81.

    Article  PubMed  Google Scholar 

  17. Guber I, McAlinden C, Majo F, Bergin C. Identifying more reliable parameters for the detection of change during the follow-up of mild to moderate keratoconus patients. Eye Vis (Lond). 2017;4:24.

    Article  PubMed  Google Scholar 

  18. McAlinden C, Khadka J, Pesudovs K. A comprehensive evaluation of the precision (repeatability and reproducibility) of the Oculus Pentacam HR. Invest Ophthalmol Vis Sci. 2011;52(10):7731–7.

    Article  PubMed  Google Scholar 

  19. Brunner M, Czanner G, Vinciguerra R, Romano V, Ahmad S, Batterbury M, et al. Improving precision for detecting change in the shape of the cornea in patients with keratoconus. Sci Rep. 2018;8(1):12345.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Gustafsson I, Bergstrom A, Myers AC, Ivarsen A, Hjortdal J. Association between keratoconus disease severity and repeatability in measurements of parameters for the assessment of progressive disease. PLoS ONE. 2020;15(2):e0228992.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Heijl A, Bengtsson B. The effect of perimetric experience in patients with glaucoma. Arch Ophthalmol. 1996;114(1):19–22.

    Article  CAS  PubMed  Google Scholar 

  22. Belin MW, Duncan JK. Keratoconus: the ABCD Grading System. Klin Monbl Augenheilkd. 2016;233(6):701–7.

    Article  CAS  PubMed  Google Scholar 

  23. Amsler M. [Classic keratocene and crude keratocene; Unitary arguments]. Ophthalmologica. 1946;111(2–3):96–101.

    Article  CAS  PubMed  Google Scholar 

  24. Bland JM, Altman DG. Measurement error proportional to the mean. BMJ. 1996;313(7049):106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. McGraw KO, Wong SP. Forming inferences about some intraclass correlation coefficients. Psychol Methods. 1996;1(1):30–46.

    Article  Google Scholar 

  26. Shoukri MM, Colak D, Kaya N, Donner A. Comparison of two dependent within subject coefficients of variation to evaluate the reproducibility of measurement devices. BMC Med Res Methodol. 2008;8:24.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Gustafsson I, Faxén T, Vicente A, Bergström A, Ivarsen A, Hjortdal J. An inter-day assessment of the ABC parameters in the evaluation of progressive keratoconus. Sci Rep. 2021;11(1):16037.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. du Toit R, Vega JA, Fonn D, Simpson T. Diurnal variation of corneal sensitivity and thickness. Cornea. 2003;22(3):205–9.

    Article  PubMed  Google Scholar 

  29. Read SA, Collins MJ. Diurnal variation of corneal shape and thickness. Optom Vis Sci. 2009;86(3):170–80.

    Article  PubMed  Google Scholar 

  30. Cicek A, Demirtas AA, Ozsaygili C, Duru Z, Ulusoy DM, Duru N, et al. Diurnal variation of anterior segment parameters handled with Scheimpflug imaging in keratoconus patients. Int Ophthalmol. 2020;40(6):1481–5.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

N/A.

Funding

Open access funding provided by Lund University.

Author information

Authors and Affiliations

Authors

Contributions

I.G., J.H. and A.I. proposed the idea. I.G. made the measurements with the Pentacam HR. I.G. and D.B. analysed data and drafted the manuscript. I.G., D.B., J.H. and A.I. analysed data and revised the manuscript. J.H. and A.I. supervised the project. All authors have read, reviewed and approved the final manuscript.

Corresponding author

Correspondence to Ingemar Gustafsson.

Ethics declarations

Ethics approval and consent to participate

All participants were given written information on the study, and an informed written consent was obtained. The study was conducted at the Department of Ophthalmology at Skåne University Hospital, Lund, Sweden, according to the declaration of Helsinki. The Swedish Ethical Review Authority in Lund, Sweden, approved the studies (No. 2015/373).

Consent for publication

N/A.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gustafsson, I., Bizios, D., Ivarsen, A. et al. Patient experience and repeatability of measurements made with the Pentacam HR in patients with keratoconus. BMC Ophthalmol 23, 201 (2023). https://doi.org/10.1186/s12886-023-02930-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12886-023-02930-4

Keywords