Parameter estimation for evaporation-driven tear film model in two space dimensions

Technische Samenvatting: Parameterschatting voor een Verdamping-gedreven Traanfilmmodel in Twee Ruimtelijke Dimensies

Probleemstelling
De verstoring van de traanfilm (TF), bekend als traanfilmruptuur (TBU), is een primair mechanisme dat onderliggend is aan droge-oogaandoeningen (DED). Hoewel verdampingsgedreven verdunning wordt erkend als een belangrijke drijvende kracht van TBU, blijven kwantitatieve in vivo schattingen van de fysische parameters die dit proces beheersen beperkt. Eerdere studies hebben gebruik gemaakt van één-dimensionale (1D) modellen of ruimtelijk samengevoegde gewone differentiaalvergelijkingen (ODE-modellen) om fluorescentie (FL) beeldvormingsdata te interpreteren. Echter, 1D-modellen slagen er vaak niet in om de asymmetrische of multi-spot aard van TBU-gebieden vast te leggen, terwijl complexe twee-dimensionale (2D) partiële differentiaalvergelijkingen (PDE-modellen) historisch gezien computationeel onhaalbaar zijn geweest voor fitting aan experimentele data, waardoor de analyse beperkt bleef tot een klein aantal verdunningsgebeurtenissen.

Methodologie
Dit werk presenteert een raamwerk voor het schatten van parameters in een 2D verdampingsgedreven TF-verdunningsmodel met behulp van experimentele FL-beeldvormingsdata van normale proefpersonen. De aanpak integreert de volgende componenten:

1. Wiskundig Model: De studie maakt gebruik van het 2D PDE-systeem ontwikkeld door Braun et al. en uitgebreid door Chen et al., dat de TF modelleert als een Newtonse vloeistof over een vlak hoornvliesoppervlak. Het systeem volgt de TF-dikte ($h$), druk ($p$), osmolariteit ($c$) en fluoresceïneconcentratie ($f$). De besturende vergelijkingen omvatten smeringstheorie, verdamping ($J$) en solutetransport. De verdampingssnelheid wordt gemodelleerd als een ruimtelijk variërende functie, hier geparametriseerd als één of meer gelokaliseerde Gaussische pieken (elliptische vlekken) om inhomogeniteiten in de lipidenlaag weer te geven.
2. Dimensiereductie: Om de computationele kosten van het herhaaldelijk oplossen van de besturende PDE's tijdens optimalisatie te overwinnen, maken de auteurs gebruik van Proper Orthogonal Decomposition (POD). Deze techniek projecteert het hoog-dimensionale dynamische systeem op een laag-dimensionale deelruimte opgespannen door dominante modi die zijn ontleend via Singular Value Decomposition (SVD), waardoor de voorwaartse solver aanzienlijk wordt versneld.
3. Invers Probleem en Optimalisatie: De kerntaak is een invers probleem: het schatten van de parameters van de verdampingsfunctie $J(x,y)$ (inclusief basisverdamping, piekhoogte, locatie, elliptische afmetingen en oriëntatie) die het $L_2$-normverschil minimaliseren tussen de theoretische FL-intensiteit (afgeleid uit het model) en experimentele FL-intensiteitsdata. De optimalisatie maakt gebruik van gradiëntvrije algoritmen (voornamelijk PRAXIS, af en toe Nelder-Mead) uit het NLopt-pakket.
4. Data Preprocessing: Experimentele videodata ondergaat stabilisatie om het bewegende TBU-gebied te volgen, gladmaking via Gaussische filters, en periodieke randvoorwaarden om te voldoen aan de aannames van het model. De data wordt geïnterpoleerd op een $40 \times 40$ ruimtelijk rooster.

Belangrijkste Bijdragen

• 2D Parameterschatting: Het artikel past succesvol een op 2D PDE gebaseerd model toe op in vivo FL-beeldvormingsdata, en gaat voorbij aan de beperkingen van 1D radiale of streepmodellen.
• Computationele Efficiëntie: Door POD te combineren met gradiëntvrije optimalisatie, tonen de auteurs aan dat het fitten van complexe 2D-modellen aan experimentele data computationeel haalbaar is, waarbij typische voorwaartse oplossingen ongeveer 15 seconden duren op een $40 \times 40$ rooster.
• Kwantitatieve Parameterherwinning: Het raamwerk levert schattingen op voor verdampingsgerelateerde en verdunningsparameters binnen TBU-gebieden die eerder niet in vivo zijn bepaald.

Resultaten
De auteurs hebben de methodologie getest op zowel synthetische data als experimentele data van 25 normale proefpersonen (specifiek analyserend Gevallen 1–5 die diverse TBU-morfologieën vertegenwoordigen):

• Synthetische Validatie: Op synthetische data herstelde de 2D-optimalisatie de ware parameters met hoge precisie (tot 8 cijfers). Bij vergelijking met een 1D radiaal model dat naar 2D is opgeheven, toonde de 2D elliptische fit aanzienlijk lagere relatieve fouten (onder de 2% versus groeiend tot meer dan 5% voor het 1D-model).
• Experimentele Gevallen:
  • Geval 1 (Bijna-cirkelvormig): De 2D-fit bereikte relatieve residuen onder de 2%, en presteerde beter dan het 1D-model (~5% fout).
  • Gevallen 2 & 3 (Hoog excentrisch): Voor TBU-vlekken met hoge excentriciteit ($e \approx 0,75-0,77$) hield het 2D elliptische model relatieve fouten onder de 10-12%, terwijl de fouten van het 1D radiale model de 12-14% overschreden.
  • Geval 4 (Ernstige non-uniformiteit): Een enkele elliptische vlek slaagde er niet in de complexe vorm vast te leggen, wat de behoefte aan flexibelere verdampingsparametrisaties benadrukt.
  • Geval 5 (Multi-spot): Het model slaagde erin een gebied met twee distincte vlekken te fitten met een configuratie van 15 parameters, waarbij relatieve fouten voor het grootste deel van de simulatietijd onder de 10% bleven.

Betekenis en Beweringen
De auteurs stellen dat deze resultaten in vivo schattingen opleveren van verdampings- en verdunningsparameters binnen TBU-gebieden, wat het begrip van TF-verdunning en droge-vlekvorming verbetert. Het werk vestigt een kwantitatieve basislijn voor vergelijking met data van patiënten met droge ogen. Het artikel beweert bescheiden dat hoewel de huidige vaste verdampingsparametrisatie de fit voor sterk onregelmatige vormen beperkt, de aangetoonde efficiëntie van het 2D POD-gebaseerde raamwerk de analyse van een groter scala aan traanfilm-dynamica mogelijk maakt dan eerder mogelijk was. Toekomstig werk wordt geïdentificeerd als het uitbreiden van de verdampingsparametrisatie om meer flexibele ruimtelijke structuren mogelijk te maken en het toepassen van het raamwerk op populaties van patiënten met droge ogen om te bepalen of deze parameters dienen als klinisch zinvolle markers.