<?xml version="1.0" encoding="UTF-8"?><rss xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:content="http://purl.org/rss/1.0/modules/content/" xmlns:atom="http://www.w3.org/2005/Atom" version="2.0" xmlns:media="http://search.yahoo.com/mrss/"><channel><title><![CDATA[PEARL-DGS Formula]]></title><description><![CDATA[The PEARL project : Resources to understand and design IOL calculation formulas for cataract surgery.]]></description><link>https://news.iolsolver.com/</link><image><url>https://news.iolsolver.com/favicon.png</url><title>PEARL-DGS Formula</title><link>https://news.iolsolver.com/</link></image><generator>Ghost 3.38</generator><lastBuildDate>Wed, 22 Apr 2026 14:47:45 GMT</lastBuildDate><atom:link href="https://news.iolsolver.com/rss/" rel="self" type="application/rss+xml"/><ttl>60</ttl><item><title><![CDATA[How to easily perform the Haigis triple optimization]]></title><description><![CDATA[<p></p><h2 id="the-haigis-formula">The Haigis Formula</h2><p>The Haigis formula is the <em>simplest thin lens formula possible</em>. The ELP (= <em>d</em> value) is predicted using a multiple regression taking ACD and AL as inputs.</p><pre><code class="language-python">def Haigis(AL, R, ACD,a0,a1,a2,power): 
    n = 1.336 
    nc = 1.3315 
    dx = 0.012 
    d = a0 + a1*</code></pre>]]></description><link>https://news.iolsolver.com/how-to-easily-perform-the-haigis-triple-optimization/</link><guid isPermaLink="false">5fcff602cbe72a009e3abd1b</guid><dc:creator><![CDATA[Guillaume Debellemanière]]></dc:creator><pubDate>Thu, 29 Dec 2022 15:54:36 GMT</pubDate><content:encoded><![CDATA[<p></p><h2 id="the-haigis-formula">The Haigis Formula</h2><p>The Haigis formula is the <em>simplest thin lens formula possible</em>. The ELP (= <em>d</em> value) is predicted using a multiple regression taking ACD and AL as inputs.</p><pre><code class="language-python">def Haigis(AL, R, ACD,a0,a1,a2,power): 
    n = 1.336 
    nc = 1.3315 
    dx = 0.012 
    d = a0 + a1*ACD + a2*AL 
    R = R/1000 
    AL = AL/1000 
    d = d/1000 
    Dc = (nc - 1)/(R) 
    qnum = n * (n - power * (AL - d)) 
    qdenum = n * (AL - d) + d * (n - power *(AL-d)) 
    q = qnum / qdenum 
    SEpred = round((q - Dc) / (1 + dx*(q - Dc)),3) 
    return SEpred</code></pre><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Haigis.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Haigis.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Haigis.png 1000w, https://news.iolsolver.com/content/images/2022/09/Haigis.png 1244w" sizes="(min-width: 720px) 720px"><figcaption>Figure from : Debellemanière G, Dubois M, Gauvin M, Wallerstein A, Brenner LF, Rampat R, Saad A, Gatinel D. The PEARL-DGS Formula: The Development of an Open-source Machine Learning-based Thick IOL Calculation Formula. Am J Ophthalmol. 2021 Dec;232:58-69.</figcaption></figure><h2 id="optimizing-the-haigis-formula">Optimizing the Haigis Formula</h2><p>The Haigis formula can be single-optimized (A0 acts like a "potentiometer" <strong>or variable resistor,</strong> like the A constant in SRK/T, until a mean prediction error equal to zero is reached) or triple-optimized.</p><p>The Haigis triple-optimization is in fact a <em>complete retraining</em> of the formula. The perfect d value is back-calculated for each eye of the dataset. This is the d value giving the real postoperative spherical equivalent when using the biometric values and the implanted IOL power in thin lens equations. </p><p>A multiple regression takes ACD and AL as inputs, <strong>then allows fitting to predict </strong>this "perfect" d value. Et voilà : the intercept is a0, the ACD coefficient is a1, and the AL coefficient is a2.</p><h2 id="back-calculation-of-the-d-value">Back-calculation of the d value</h2><p>You can back-calculate the d value for a given eye with the following function (As expected you'll need the postoperative SE, ARC, AL and the implanted IOL power).</p><pre><code class="language-python"># AL, R in meters 
def solve_d_Haigis(AL, R, power, Rx): 
    n = 1.336 
    nc = 1.3315 
    dx = 0.012 
    Dc = (nc - 1) / R 
    z = Dc + Rx/(1-Rx*dx) 
    b = -power * (AL + n/z) 
    c = (n * (AL - n/z) + power * AL * n/z) 
    d = (1 / (2*power)) * (-b - np.sqrt(b**2 - 4 * power * c) ) 
    return d</code></pre><p><em>From :  W.Haigis. "The Haigis formula", in Intraocular Lens Power Calculation, H. John Shammas, Chapter 5, edited by  Slack Incorporated, 2004.</em></p><pre><code class="language-python"># Example : 

val = solve_d_Haigis(0.024, 0.0075, 21, -2) 
print(val) 

#output : 0.005461457048808832</code></pre><p>The output of this function is in meters : multiply by 1000 before fitting your regression.</p><p>Like this : </p><pre><code class="language-python">from sklearn.linear_model import LinearRegression
df['d_Haigis'] = df['d_Haigis'] * 1000
lin = LinearRegression()
lin.fit(df[['ACD','AL']], df['d_Haigis'])
lin.intercept_

# output : 
# -0.07756916259537494</code></pre><p>This is your a0 value ! </p><pre><code class="language-python">lin.coef_

# output : 
# array([0.39819972, 0.16012471])</code></pre><p>And the coefficients are a1 and a2, respectively. What a beautifully simple formula! Genius.</p>]]></content:encoded></item><item><title><![CDATA[The (not so) secret mechanisms behind classical thin lens IOL formulas]]></title><description><![CDATA[<p>Third generation formulas are based on thin lens optics. They have the same "skeleton". The value they use to determine the corneal power from the anterior radius of curvature differ - neither do they use the same keratometric index,  nor the same rules to predict the lens position from the</p>]]></description><link>https://news.iolsolver.com/the-not-so-secret-mechanisms-of-cla/</link><guid isPermaLink="false">5fcff34acbe72a009e3abd04</guid><dc:creator><![CDATA[Guillaume Debellemanière]]></dc:creator><pubDate>Sun, 23 Oct 2022 13:09:00 GMT</pubDate><content:encoded><![CDATA[<p>Third generation formulas are based on thin lens optics. They have the same "skeleton". The value they use to determine the corneal power from the anterior radius of curvature differ - neither do they use the same keratometric index,  nor the same rules to predict the lens position from the biometric parameters. However, once those values are calculated, the prediction process is the same for every formula.</p><p>Let's dive into the details :-)</p><h2 id="the-srk-t-formula">The SRK/T Formula</h2><p>Published in 1990 by Sanders, Retzlaff and Kraff, the SRK/T (1,2) has long been the favourite IOL calculation method among ophthalmologists, and is still extremely popular.</p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/SRKT-2.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/SRKT-2.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/SRKT-2.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/SRKT-2.png 1600w, https://news.iolsolver.com/content/images/2022/09/SRKT-2.png 1984w" sizes="(min-width: 720px) 720px"><figcaption>The SRK/T formula is clever and complex. The retinal thickness (RETHICK) value depends on the AL value. AL + RETHICK gives LOPT (Optical AL), the AL value used in the optical equation. Corneal width (Cw) is a function of K and a modified AL value (LCOR). Cw and R are used to calculate the parameter H. ACD const is a transformation of the A constant value : ACD const + H - 3.336 give the lens position value. K is used natively in the optical equation, using a keratometric value of 1.333. Figure from : Debellemanière G, Dubois M, Gauvin M, Wallerstein A, Brenner LF, Rampat R, Saad A, Gatinel D. The PEARL-DGS Formula: The Development of an Open-source Machine Learning-based Thick IOL Calculation Formula. Am J Ophthalmol. 2021 Dec;232:58-69.</figcaption></figure><p>In this formula, the axial length is slightly modified before being used to predict the lens position. It allows to lower its impact on lens position prediction for high AL values.</p><pre><code class="language-python">def setLCOR(LA): 
	if LA &lt;= 24.2: 
    	LCOR = LA 
    else: 
    	LCOR = -3.446 + (1.715 * LA) - (0.0237* (LA**2))
    return LCOR</code></pre><p>To be more precise, the modified AL value (LCOR) is used, with the K value, to predict the Corneal Width (Cw). Cw and R are then used to calculate H, which is added to the transformed A constant to give the ELP.  </p><pre><code class="language-python">Cw = -5.40948 + (0.58412 * LCOR) + (0.098 * K)</code></pre><pre><code class="language-python">def setH(r, Cw):
    x = r**2 - (Cw**2 / 4.0) 
    if x &lt; 0: 
    	x = 0 
    else: 
    	pass 
    H = r - math.sqrt(x) 
    return H</code></pre><p>The influence of the AL value on lens position has a ceiling...</p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-22-a--16.23.14-1.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-22-a--16.23.14-1.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-22-a--16.23.14-1.png 1000w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-22-a--16.23.14-1.png 1360w" sizes="(min-width: 720px) 720px"><figcaption>From : Retzlaff JA, Sanders DR, Kraff MC. Development of the SRK/T intraocular lens implant power calculation formula. J Cataract Refract Surg. 1990 May;16(3):333-40. Using a quadratic function (polynomial of degree 2) to determine the lens position from the AL, avoids the hard cut-off of the Holladay 1 formula.</figcaption></figure><p>... and the predACD value has a strange relationship with the anterior corneal radius value : </p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--09.39.22.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--09.39.22.png 600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--09.39.22.png 784w" sizes="(min-width: 720px) 720px"><figcaption>SRK/T PredACD value as a function of anterior corneal radius (AL = 23.5mm, A const. 119.0)</figcaption></figure><p>This is because the value entered in the square root of the formula allowing the calculation of H is set to zero when negative (3) !</p><p>Let's use a 3D representation : </p><figure class="kg-card kg-image-card"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--08.01.23.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--08.01.23.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-23-a--08.01.23.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-23-a--08.01.23.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--08.01.23.png 1662w" sizes="(min-width: 720px) 720px"></figure><p>Also, note that the AL value used in the optical equation is modified by adding a theoretical Retinal Thickness (RETHICK) value to it. RETHICK is itself a function of AL. </p><pre><code class="language-python">RETHICK = 0.65696 - 0.02029*LA </code></pre><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--06.19.53-1.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--06.19.53-1.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-23-a--06.19.53-1.png 1000w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--06.19.53-1.png 1252w" sizes="(min-width: 720px) 720px"><figcaption>The longer the eye, the thinner the RETHICK value of the SRK/T formula.</figcaption></figure><p>This is the final SRK/T function : </p><pre><code class="language-python">def SRKT_SE(LA, r, aconst, power):
    # a keratometric index value of 1.3375 is used to compute the ACD via the Cw value
    # but : a keratometric index value of 1.333 is used as the real keratometric index in the optical calculation ! (cf ncml value)
    K = 337.5/r 
    V = 12.0 
    na = 1.336 
    nc = 1.333 
    ncml = 0.333
    # the SRK A constant is modified to represent a physical offset
    ACDconst = 0.62467 * aconst - 68.747 
    offset = ACDconst - 3.336
    # The eye-specific lens position is calculated
    LCOR = setLCOR(LA) 
    Cw = -5.40948 + (0.58412 * LCOR) + (0.098 * K) 
    H = setH(r, Cw) 
    # The ELP is added to the modified A constant
    ACDest = H + offset 
    # The optical AL is equal to AL + retinal thickness
    RETHICK = 0.65696 - 0.02029*LA 
    LOPT = LA + RETHICK 
    # the thin lens calculation is finally done
    num = 1000 * na * (na *r - ncml * LOPT)- power * (LOPT-ACDest) * (na * r - ncml * ACDest) 
    denum = na * (V*(na * r - ncml * LOPT) + LOPT * r) - 0.001 * power * (LOPT-ACDest) * (V * (na * r - ncml * ACDest) + ACDest * r) 
    SEpred = num/denum 
    return round(SEpred,3)</code></pre><p></p><h2 id="the-hoffer-q-formula">The Hoffer Q Formula</h2><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-16-a--11.13.53-1.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-16-a--11.13.53-1.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-16-a--11.13.53-1.png 1000w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-16-a--11.13.53-1.png 1570w" sizes="(min-width: 720px) 720px"><figcaption>Figure from : Debellemanière G, Dubois M, Gauvin M, Wallerstein A, Brenner LF, Rampat R, Saad A, Gatinel D. The PEARL-DGS Formula: The Development of an Open-source Machine Learning-based Thick IOL Calculation Formula. Am J Ophthalmol. 2021 Dec;232:58-69.</figcaption></figure><p>The Hoffer Q formula (4-7) is the most difficult, among the third generations formulas, to code correctly. The AL value is used to generate 3 other values, which are used along with the K value to calculate the lens position. pACD + 0.05 is added to the predicted value. K is calculated using a keratometric index of 1.3375. The original K and AL values are used in the optical equations.</p><pre><code class="language-python">def defMG(LA): 
    if LA &lt;= 23.0: 
    	M = 1.0 
        G = 28.0 
    else: 
    	M = -1.0 
        G = 23.5 
        return M, G</code></pre><pre><code class="language-python">def defA(LA): 
    if LA &gt; 31.0: 
    	A = 31.0 
    elif LA &lt; 18.5: 
    	A = 18.5 
    else: 
    	A = LA 
    return A</code></pre><p>The lens position function is complicated : </p><pre><code class="language-python">import math 
def predACD(LA, K, pACD): 
    M, G = defMG(LA) 
    A = defA(LA) 
    p2 = 0.3 * (A - 23.5) 
    p3 = (math.tan(K*(math.pi / 180.0)))**2 
    p4 = 0.1 * M * ((23.5 - A)**2) 
    p5 = math.tan((0.1 * (G - A)**2) * (math.pi / 180.0)) 
    ACD = pACD + p2 + p3 + (p4 * p5) - 0.99166 
    return ACD</code></pre><figure class="kg-card kg-image-card"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--20.18.14.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--20.18.14.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-23-a--20.18.14.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-23-a--20.18.14.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--20.18.14.png 1640w" sizes="(min-width: 720px) 720px"></figure><p>Final Hoffer Q function : </p><pre><code class="language-python">def HofferQ_SE(LA, r, pACD, power): 
    K = 337.5/r 
    ACD = predACD(LA, K, pACD) 
    R = (1.336 / (1.336/(1336 / (LA - ACD - 0.05) - power) + (ACD + 0.05) / 1000.0)) - K 
    SEpred = R / (1.0 + 0.012 * R) 
    return round(SEpred, 3)</code></pre><p></p><h2 id="the-holladay-1-formula">The Holladay 1 Formula</h2><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-16-a--11.14.19-1.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-16-a--11.14.19-1.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-16-a--11.14.19-1.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-16-a--11.14.19-1.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-16-a--11.14.19-1.png 1632w" sizes="(min-width: 720px) 720px"><figcaption>The Holladay 1 formula. The AL and R are transformed, then those transformed values are used to predict the lens position. The SF is added to the predicted value. The "original" R is used in the optical equations, using a keratometric index equal to 4/3 (1.333333.....). Retinal thickness is added to the AL value before being entered in the optical equation. Figure from : Debellemanière G, Dubois M, Gauvin M, Wallerstein A, Brenner LF, Rampat R, Saad A, Gatinel D. The PEARL-DGS Formula: The Development of an Open-source Machine Learning-based Thick IOL Calculation Formula. Am J Ophthalmol. 2021 Dec;232:58-69.</figcaption></figure><p>The Holladay 1 formula (8) uses a keratometric index of 4/3. The lens position (ACD) is inferred from the ARC (with a flooring below 7.0mm) and from the AL (with a ceiling after  25.32mm, approximately).</p><pre><code class="language-python">def setRag(R): 
    if R &lt; 7: 
    	Rag = 7 
    else: 
    	Rag = R 
    return Rag</code></pre><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--06.53.20-1.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--06.53.20-1.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-23-a--06.53.20-1.png 1000w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--06.53.20-1.png 1222w" sizes="(min-width: 720px) 720px"><figcaption>Rag is equal to R, floored at 7 mm.&nbsp;</figcaption></figure><pre><code class="language-python">def setAG(LA): 
    AG = 12.5 * LA / 23.45 
    if AG &gt; 13.5: 
    	AG = 13.5 
    else: 
        AG = AG 
    return AG
    
# the AG has a ceiling when greater than 13.5 : this corresponds to an AL value of 25.32mm</code></pre><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--06.59.55.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--06.59.55.png 600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--06.59.55.png 816w" sizes="(min-width: 720px) 720px"><figcaption>AG is a function of AL, with a ceiling at 13.5</figcaption></figure><p>The lens position, ACD, is a function of Rag and AG : </p><pre><code class="language-python">def setACD(Rag, AG): 
    ACD = 0.56 + Rag - (math.sqrt( (Rag**2) -((AG**2) /4) )) 
    return ACD</code></pre><figure class="kg-card kg-image-card"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--07.27.05.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--07.27.05.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-23-a--07.27.05.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-23-a--07.27.05.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--07.27.05.png 1736w" sizes="(min-width: 720px) 720px"></figure><p>A retinal thickness value of 0.2mm is added to the AL.</p><p>Final Holladay 1 function : </p><pre><code class="language-python">def Holladay_SE(LA, R, SF, power): 
    Rag = setRag(R) 
    AG = setAG(LA) 
    ACD = setACD(Rag, AG) 
    nc = 4/3 
    na = 1.336 
    RT = 0.2
    V = 12 
    ALm = LA + RT 
    num = 1000*na*(na*R - (nc - 1)*ALm) - power*(ALm - ACD - SF)*(na*R - (nc - 1)*(ACD + SF)) denum = na*(V*(na*R-(nc-1)*ALm)+ALm*R)-0.001*power*(ALm-ACD-SF)*(V*(na*R-(nc-1)*(ACD+SF))+(ACD+SF)*R) result = round(num / denum, 3) 
    return result</code></pre><p></p><h2 id="the-haigis-formula">The Haigis Formula</h2><p>The Haigis formula (9) is the most simple third generation formula... and also consistently has the best outcomes when using its triple-optimized version. The lens position d is obtained using a multiple regression taking ACD and AL as inputs :</p><p><strong><em>d = a0 + a1 x ACD + a2 x AL</em></strong></p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Haigis-3.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Haigis-3.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Haigis-3.png 1000w, https://news.iolsolver.com/content/images/2022/09/Haigis-3.png 1244w" sizes="(min-width: 720px) 720px"><figcaption>The Haigis formula is beautifully simple. The lens position value (d) is a function of ACD and AL, using simple coefficients (a1 and a2, respectively) and a constant (a0). K is calculated from R using a keratometric index of 1.3315. Neither AL or K are transformed before being used in the optical equation.</figcaption></figure><figure class="kg-card kg-image-card"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--08.08.28.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-23-a--08.08.28.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-23-a--08.08.28.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-23-a--08.08.28.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-23-a--08.08.28.png 1672w" sizes="(min-width: 720px) 720px"></figure><p>What's your favourite approach ?</p><h2 id="references">References</h2><p>1. Retzlaff JA, Sanders DR, Kraff MC. Development of the SRK/T intraocular lens implant power calculation formula. J Cataract Refract Surg. 1990 May;16(3):333-40.</p><p>2. Retzlaff JA, Sanders DR, Kraff MC. Development of the SRK/T intraocular lens implant power calculation formula: Erratum. J Cataract Refract Surg. 1990;16(4):528.</p><p>3. Haigis W. Occurrence of erroneous anterior chamber depth in the SRK/T formula. J Cataract Refract Surg. 1993 May;19(3):442-6.</p><p>4. Hoffer, Kenneth J. M.D.<sup>a</sup>. The Hoffer Q formula: A comparison of theoretic and regression formulas. Journal of Cataract &amp; Refractive Surgery: November 1993 - Volume 19 - Issue 6 - p 700-712</p><p>5. Hoffer, Kenneth J. MD. Reply: Errata in printed Hoffer Q formula. Journal of Cataract &amp; Refractive Surgery: January 2007 - Volume 33 - Issue 1 - p 2-3<br>doi: 10.1016/j.jcrs.2006.08.056</p><p>6. Zuberbuhler B, Morrell AJ. Errata in printed Hoffer Q formula. J Cataract Refract Surg. 2007 Jan;33(1):2; author reply 2-3.</p><p>7. Hoffer, K. Errors in self-programming the Hoffer Q formula. <em>Eye</em><strong>21</strong>, 429 (2007).</p><p>8. Holladay JT, Musgrove KH, Prager TC, Lewis JW, Chandler TY, Ruiz RS. A three-part system for refining intraocular lens power calculations. J Cataract Refract Surg. 1988;14(1):17–24.</p><p>9. Haigis W, Lege B, Miller N, Schneider B. Comparison of immersion ultrasound biometry and partial coherence interferometry for intraocular lens calculation according to Haigis. Graefes Arch Clin Exp Ophthalmol. 2000;238(9): 765–773.</p>]]></content:encoded></item><item><title><![CDATA[IOL calculation before / after DMEK surgery]]></title><description><![CDATA[<h2></h2><p>Descemet membrane endothelial keratoplasty (DMEK) was introduced in 2006 by Melles et al (1). This procedure revolutionized the management of endothelial pathologies by allowing a near perfect anatomical replacement of the diseased endothelium.</p><p>In most of the cases, the diseased (pre-DMEK) posterior corneal surface is (fortunately !) not taken into account</p>]]></description><link>https://news.iolsolver.com/iol-calculation-before-after-dmek-surgery/</link><guid isPermaLink="false">6028c78ecbe72a009e3abde8</guid><dc:creator><![CDATA[Guillaume Debellemanière]]></dc:creator><pubDate>Mon, 17 Oct 2022 07:58:39 GMT</pubDate><content:encoded><![CDATA[<h2></h2><p>Descemet membrane endothelial keratoplasty (DMEK) was introduced in 2006 by Melles et al (1). This procedure revolutionized the management of endothelial pathologies by allowing a near perfect anatomical replacement of the diseased endothelium.</p><p>In most of the cases, the diseased (pre-DMEK) posterior corneal surface is (fortunately !) not taken into account in IOL formula calculations. As this posterior surface is restored back to its physiological shape by DMEK surgery, precise IOL calculation in triple procedure (i.e. prediction of the adequate IOL power for a post-DMEK eye using measurements made on a pre-DMEK eye) should be possible without any surprise. </p><p>However, a so-called "hyperopic shift" is observed both in standalone DMEK and in DMEK triple procedure. Why ? What is the difference between those phenomena ?</p><h2 id="hyperopic-shift-in-dmek-standalone-procedure">Hyperopic shift in DMEK standalone procedure</h2><p>The refractive shift observed in DMEK single procedure is logical and easy to understand. The pathological, flattened posterior corneal surface is restored to its physiological value by the surgery. The posterior corneal radius decreases, thus decreasing the corneal power*. The eye is more hyperopic after the surgery (2-4).</p><p><em>*An increase in anterior corneal radius decreases the corneal power ; an increase in posterior corneal radius <strong>increases</strong> the corneal power.</em></p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.26.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.26.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.26.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.26.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.26.png 1778w" sizes="(min-width: 720px) 720px"><figcaption>A healthy cornea. The green line represents the physiological posterior curvature&nbsp;</figcaption></figure><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.34.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.34.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.34.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.34.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.34.png 1862w" sizes="(min-width: 720px) 720px"><figcaption>A decompensated cornea. The posterior radius is flatter, the corneal power is higher.</figcaption></figure><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.42.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.42.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.42.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.42.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.42.png 1786w" sizes="(min-width: 720px) 720px"><figcaption>After DMEK surgery : the posterior corneal radius decreases and returns to its physiological state, thus decreasing the corneal power. An hyperopic shift happens (but no IOL calculation error should happen in this case !)</figcaption></figure><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.49.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.49.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.49.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.49.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.49.png 1798w" sizes="(min-width: 720px) 720px"><figcaption>Theoretically, if the post-DMEK posterior radius is smaller than the physiological radius, a greater hyperopic shift should be observed - or a small hyperopic IOL calculation error.</figcaption></figure><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.56.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.56.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.56.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.56.png 1600w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.11.56.png 1844w" sizes="(min-width: 720px) 720px"><figcaption>And if the post-DMEK posterior radius is higher than the physiological radius, a small hyperopic shift and a small myopic IOL calculation error should be observed.</figcaption></figure><p></p><p>Does this phenomenon have an impact in IOL calculation? Theoretically, no : let's see why.</p><h2 id="restoring-the-posterior-corneal-surface-to-its-physiological-value-does-not-impact-iol-calculation">Restoring the posterior corneal surface to its physiological value does not impact IOL calculation</h2><p>As explained in 2011 by Ham et al. (5), biometric measurements do not usually take into account the posterior corneal radius. The corneal power is extrapolated from the anterior corneal radius, and (it was initially thought that) no modification of the anterior corneal radius happens after DMEK surgery. If the anterior radius is not modified and the posterior radius goes back to its physiological value, a <em>refractive shift</em> (in a single procedure) should be observed, but <em>no IOL calculation error</em> (in DMEK triple procedure) should occur.</p><p>Note : Of course this is <strong>not</strong> true for ray-tracing methods or total corneal power calculations using a measured posterior corneal radius, available in the most recent biometers. <strong>Those methods should by avoided at all costs in DMEK-triple procedure</strong>, because the posterior corneal radius <em>will</em> change after DMEK surgery.</p><p>Since 2011, various authors have in fact described a modification of the anterior corneal radius, and steeper-than-normal posterior corneal radii (6-11); IOL calculation error in triple-DMEK has also been well described by different teams (12-15). If the anterior radius is indeed different at the time of the biometry and after DMEK, it is logical to observe an error. Similarly, if the posterior corneal radius is steeper-than-normal after DMEK, the usual relationship between the anterior radius and the posterior radius is no longer valid and standard keratometric indices will be flawed. Both phenomena can explain an IOL calculation error.</p><p>As explained in a recent paper by our team, it is then necessary to distinguish two different phenomena (17) : </p><ul><li><strong>DMEK-induced refractive shift (DIRS)</strong> : the (usually hyperopic) refractive shift observed in standalone procedure, which is explained by the difference between the pre-DMEK corneal power and the post-DMEK corneal power. <em>It is a pure modification of corneal power.</em></li><li><strong>DMEK-induced IOL calculation error (DICE)</strong> : the IOL calculation error happening because of anterior corneal radius modifications before and after DMEK, and non-physiological post-DMEK posterior radii. <em>It is the difference between the corneal power estimated by IOL formulas using pre-DMEK measurements and the real post-DMEK corneal power. </em>IOL formulas are flawed because the anterior radius changes (or is incorrectly measured pre-DMEK), or because the anterior/posterior relationship is broken, or both.</li></ul><p><strong><em>Surprisingly, both phenomena are usually referred to as "hyperopic shift" in the literature. They should be differentiated !</em></strong></p><p>NB : it could also be hypothesized that the refractive index of the cornea is also non-physiological after DMEK. This phenomenon is very difficult to isolate in clinical practice, but could also partly explain both DIRS and DICE.</p><h2 id="relative-impact-of-anterior-and-posterior-corneal-radii-variations-on-corneal-power">Relative impact of anterior and posterior corneal radii variations on corneal power</h2><p>Anterior and posterior corneal radii variations both influence the corneal power. But what is the relative importance of both phenomena ?</p><p>It should be remembered that the difference between the refractive index of the air (n_air = 1) and of the cornea (n_cor = 1.376) is much greater than the difference between the refractive index of the cornea and of the aqueous (n_aq = 1.336) : 0.376 and 0.04, respectively. </p><p>For a theoretical cornea (ARC = 7.5mm, PRC = 6.5mm, CCT = 0.520mm), a variation of 0.1mm in ARC leads to a difference in corneal power of 0.68D, while a the same variation of PRC modifies it by 0.1D. This is almost a sevenfold difference! </p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.22.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.22.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.22.png 1000w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.22.png 1410w" sizes="(min-width: 720px) 720px"><figcaption>The variation of the posterior corneal radius alone is not enough in itself to explain the refractive shift or the IOL calculation error observed in DMEK.</figcaption></figure><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.31.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.31.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.31.png 1000w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.31.png 1412w" sizes="(min-width: 720px) 720px"><figcaption>If the anterior radius decreases after DMEK, the cornea is more powerful and a myopic shift occurs.</figcaption></figure><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.38.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.38.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.38.png 1000w, https://news.iolsolver.com/content/images/2022/09/Capture-d-e-cran-2022-09-11-a--12.08.38.png 1406w" sizes="(min-width: 720px) 720px"><figcaption>If the anterior radius increases after DMEK, the cornea is less powerful and an hyperopic shift occurs.‌</figcaption></figure><p></p><h2 id="is-it-possible-to-predict-the-variation-in-anterior-corneal-radius-and-or-the-post-dmek-anterior-posterior-ratio">Is it possible to predict the variation in anterior corneal radius and/or the post-DMEK anterior/posterior ratio ?</h2><p>DIRS can't be prevented : the corneal power will always be modified by DMEK surgery. By contrast, DICE could theoretically be avoided (by correcting the IOL calculation) if its underlying mechanisms were clearly understood.</p><p>Flatter and more oblate posterior corneal surfaces, as well as a higher posterior to anterior corneal curvature radii ratio have been cited as DIRS or DICE risk factors (8,15,16). </p><p>In a recent study (17), we did not find any strong predictor for DIRS or DICE. However, we found, like Schoenberg et al. (14), a lower variation of the ARC measured before and after DMEK when using a Scheimpflug camera, compared to the measurement obtained with the IOLMaster. Could the "variation" in ARC be partly explained by an ARC measurement error made by our biometers, which could be lowered by using the ARC measured by Scheimpflug camera in our IOL formulas ?</p><h2 id="conclusion">Conclusion</h2><ul><li>DIRS is the DMEK-induced variation in corneal power, observed in DMEK single procedure, which can't be avoided. Any anterior or posterior corneal radius modification induces a DIRS. It is not preventable.</li><li>DICE is the IOL calculation error induced by modifications/measurement errors of the ARC, and/or non-physiological post-DMEK posterior corneal radii. Any anterior corneal radius modification induces a DICE, but, unlike DIRS, posterior corneal radius modifications are involved in DIRS only when this radius is non-physiological after DMEK. It is theoretically preventable.</li><li>It could be wise to aim for a more myopic target when the posterior corneal surface is very flat and oblate. The ARC value measured using the Pentacam device should be preferred to the ARC value obtained using a biometer. However, no strong DICE predictor has been found to date.</li></ul><h2 id="references-">References : </h2><ol><li>Melles GR. Posterior lamellar keratoplasty: DLEK to DSEK to DMEK. Cornea. 2006 Sep;25(8):879-81. </li><li>Dunker SL, Veldman MHJ, Winkens B, van den Biggelaar FJHM, Nuijts RMMA, Kruit PJ, Dickman MM; Dutch Cornea Consortium. Real-World Outcomes of DMEK: A Prospective Dutch registry study. Am J Ophthalmol. 2021 Feb;222:218-225.</li><li>Price MO, Giebel AW, Fairchild KM, Price FW Jr. Descemet's membrane endothelial keratoplasty: prospective multicenter study of visual and refractive outcomes and endothelial survival. Ophthalmology. 2009</li><li>Agha B, Ahmad N, Dawson DG, Kohnen T, Schmack I. Refractive outcome and tomographic changes after Descemet membrane endothelial keratoplasty in pseudophakic eyes with Fuchs' endothelial dystrophy. Int Ophthalmol. 2021 Aug;41(8):2897-2904.</li><li>Ham L, Dapena I, Moutsouris K, Balachandran C, Frank LE, van Dijk K, Melles GR. Refractive change and stability after Descemet membrane endothelial keratoplasty. Effect of corneal dehydration-induced hyperopic shift on intraocular lens power calculation. J Cataract Refract Surg. 2011 Aug;37(8):1455-64.</li><li>Alnawaiseh M, Rosentreter A, Eter N, Zumhagen L. Changes in Corneal Refractive Power for Patients With Fuchs Endothelial Dystrophy After DMEK. Cornea. 2016 Aug;35(8):1073-7.</li><li>Arnalich-Montiel F, Mingo-Botín D, Diaz-Montealegre A. Keratometric, Pachymetric, and Surface Elevation Characterization of Corneas With Fuchs Endothelial Corneal Dystrophy Treated With DMEK. Cornea. 2019;38:535-541.</li><li>Cheung AY, Chachare DY, Eslani M, et al. Tomographic changes in eyes with hyperopic shift after triple Descemet membrane endothelial keratoplasty. J Cataract Refract Surg. 2018;44:738-744.</li><li>Parker J, Dirisamer M, Naveiras M, et al. Outcomes of Descemet membrane endothelial keratoplasty in phakic eyes. J Cataract Refract Surg. 2012;38:871-877.</li><li>van Dijk K, Ham L, Tse WHW, et al. Near complete visual recovery and refractive stability in modern corneal transplantation: Descemet membrane endothelial keratoplasty (DMEK). Cont Lens Anterior Eye. 2013;36:13-21.</li><li>van Dijk K, Rodriguez-Calvo-de-Mora M, van Esch H, et al. Two-Year Refractive Outcomes After Descemet Membrane Endothelial Keratoplasty. Cornea. 2016;35:1548-1555.</li><li>Augustin VA, Weller JM, Kruse FE, Tourtas T. Refractive Outcomes After Descemet Membrane Endothelial Keratoplasty + Cataract/Intraocular Lens Triple Procedure: A Fellow Eye Comparison. Cornea. 2021 Jul 1;40(7):883-887</li><li>Bae SS, Ching G, Holland S, McCarthy M, Ritenour R, Iovieno A, Yeung SN. Refractive Outcomes of Descemet Membrane Endothelial Keratoplasty Combined With Cataract Surgery in Fuchs Endothelial Dystrophy. J Refract Surg. 2020 Oct 1;36(10):661-666.</li><li>Schoenberg ED, Price FW Jr, Miller J, et al. Refractive outcomes of Descemet membrane endothelial keratoplasty triple procedures (combined with cataract surgery). J Cataract Refract Surg. 2015;41:1182-1189.</li><li>Fritz M, Grewing V, Böhringer D, et al. Avoiding Hyperopic Surprises After Descemet Membrane Endothelial Keratoplasty in Fuchs Dystrophy Eyes by Assessing Corneal Shape. Am J Ophthalmol. 2019;197:1-6.</li><li>Diener R, Eter N, Alnawaiseh M. Using the posterior to anterior corneal curvature radii ratio to minimize the risk of a postoperative hyperopic shift after Descemet membrane endothelial keratoplasty. Graefes Arch Clin Exp Ophthalmol. January 2020. </li><li>Debellemanière G, Ghazal W, Dubois M, et al. Descemet Membrane Endothelial Keratoplasty–Induced Refractive Shift and Descemet Membrane Endothelial Keratoplasty–Induced Intraocular Lens Calculation Error. Cornea 2022.</li></ol>]]></content:encoded></item><item><title><![CDATA[How to optimize the SRK/T A constant]]></title><description><![CDATA[<p></p><h2 id="the-srk-t-formula">The SRK/T formula</h2><p>The SRK/T formula (1,2) has long been one of the most popular IOL calculation methods : it is still widely used and its optimized constant is often used to tune other formulas.</p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/SRKT-1.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/SRKT-1.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/SRKT-1.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/SRKT-1.png 1600w, https://news.iolsolver.com/content/images/2022/09/SRKT-1.png 1984w" sizes="(min-width: 720px) 720px"><figcaption>Inner working of the SRK/T formula. Figure from : Debellemanière G, Dubois M,</figcaption></figure>]]></description><link>https://news.iolsolver.com/how-to-optimize-the-srk-t-a-constant/</link><guid isPermaLink="false">5fcff646cbe72a009e3abd1f</guid><dc:creator><![CDATA[Guillaume Debellemanière]]></dc:creator><pubDate>Thu, 06 Oct 2022 07:17:54 GMT</pubDate><content:encoded><![CDATA[<p></p><h2 id="the-srk-t-formula">The SRK/T formula</h2><p>The SRK/T formula (1,2) has long been one of the most popular IOL calculation methods : it is still widely used and its optimized constant is often used to tune other formulas.</p><figure class="kg-card kg-image-card kg-card-hascaption"><img src="https://news.iolsolver.com/content/images/2022/09/SRKT-1.png" class="kg-image" alt srcset="https://news.iolsolver.com/content/images/size/w600/2022/09/SRKT-1.png 600w, https://news.iolsolver.com/content/images/size/w1000/2022/09/SRKT-1.png 1000w, https://news.iolsolver.com/content/images/size/w1600/2022/09/SRKT-1.png 1600w, https://news.iolsolver.com/content/images/2022/09/SRKT-1.png 1984w" sizes="(min-width: 720px) 720px"><figcaption>Inner working of the SRK/T formula. Figure from : Debellemanière G, Dubois M, Gauvin M, Wallerstein A, Brenner LF, Rampat R, Saad A, Gatinel D. The PEARL-DGS Formula: The Development of an Open-source Machine Learning-based Thick IOL Calculation Formula. Am J Ophthalmol. 2021 Dec;232:58-69.</figcaption></figure><h2 id="what-is-an-optimization">What is an optimization ?</h2><p>Optimizing a formula involves shifting its constant up or down, to obtain a mean prediction error equal to zero. </p><h2 id="srk-t-optimization-process">SRK/T Optimization process</h2><p>The following code allows you to obtain the optimized A constant for a given dataset. </p><p>The OptimizeSRKT() function takes various variables as inputs : dataframe is the pd.DataFrame() object containing the data; col_AL , col_R, col_power, col_SE are the names of the related dataframe columns (string);  aconst_start (float) should be the closest value you can guess for the adjusted A constant. The SRKTES() function calculates the SRK/T SE prediction for a given eye / implanted IOL power. </p><p>The function will calculate the predicted SRK/T spherical equivalent, calculate the prediction error for each eye, then its mean for the whole dataset ; by adding the mean error to the A constant (it's a trick !), we'll quickly find the A constant leading to a prediction error equal to zero, which is the adjusted A constant.</p><pre><code class="language-python">def OptimizeSRKT(dataframe, col_AL, col_R, col_power, col_SE, aconst_start):
    dataframe['aconst'] = aconst_start
    dataframe['SRKTES'] = dataframe[[col_AL, col_R, 'aconst', col_power]].apply(lambda x : SRKTES(*x), axis = 1)
    dataframe['error'] = dataframe[col_SE] - dataframe['SRKTES']
    mean_error = dataframe.error.mean()
    while mean_error &gt; 0.00001:
        dataframe['aconst'] += mean_error
        dataframe['SRKTES'] = dataframe[[col_AL, col_R, 'aconst', col_power]].apply(lambda x : SRKTES(*x), axis = 1)
        dataframe['error'] = dataframe[col_SE] - df['SRKTES']
        mean_error = dataframe.error.mean()
        print('Mean error : ', mean_error)
    print('Optimized SRK/T A constant : ', dataframe.aconst.unique()[0])</code></pre><pre><code class="language-python"># Example : 

OptimizeSRKT(df, 'AL', 'R', 'Power', 'ES_postop', 118.0)

# Output : 
# Mean error :  0.0773036659614707
# Mean error :  0.007425076192446311
# Mean error :  0.0007240370609196819
# Mean error :  7.070440730831753e-05
# Mean error :  6.90546869455154e-06
# Optimized SRK/T A constant :  119.08180447676536</code></pre><h2 id="references">References</h2><p>1- Retzlaff JA, Sanders DR, Kraff MC. Development of the SRK/T intraocular lens implant power calculation formula. J Cataract Refract Surg. 1990 May;16(3):333-40.</p><p>2- Retzlaff JA, Sanders DR, Kraff MC. Development of the SRK/T intraocular lens implant power calculation formula: Erratum. J Cataract Refract Surg. 1990;16(4):528.</p>]]></content:encoded></item><item><title><![CDATA[PEARL Online Calculator User Manual]]></title><description><![CDATA[<p>In this post, we will review the PEARL online calculator features, discuss the mechanisms behind them, and their applications.</p><h3 id="regular-eyes">Regular Eyes</h3><p>This page allows the user to perform a calculation for 'standard eyes', i.e. eyes which have not undergone refractive surgery, have no corneal scar, and for which information</p>]]></description><link>https://news.iolsolver.com/pearl-online-calculator-user-manual/</link><guid isPermaLink="false">63194c6dcbe72a009e3ac2d2</guid><dc:creator><![CDATA[Guillaume Debellemanière]]></dc:creator><pubDate>Wed, 28 Sep 2022 09:53:13 GMT</pubDate><content:encoded><![CDATA[<p>In this post, we will review the PEARL online calculator features, discuss the mechanisms behind them, and their applications.</p><h3 id="regular-eyes">Regular Eyes</h3><p>This page allows the user to perform a calculation for 'standard eyes', i.e. eyes which have not undergone refractive surgery, have no corneal scar, and for which information about the first operated eye is not available.</p><p>The <strong>"keratometric index"</strong> field allows selection of the index used by the user's biometer. In most cases, it is 1.3375 ; however, some biometers use a different value. The PEARL methodology is based on thick lens equations and does not use the K value delivered by the biometer in its equations, but rather the anterior corneal radius of curvature measurement. It is then fundamental to select the right keratometric index value, because it is used by the calculator to determine this anterior corneal radius of curvature. </p><p>The <strong>"Biometer"</strong> field is used to apply slight changes to the AL measurement. If "Argos" or "US (immersion)" are selected, the approximation of the Sum-of-Segments AL utilised in the PEARL formula is not used. If US (applanation) is selected, a correction factor of 0.15 mm (1-4) is added to the AL measurement to compensate for the underestimation due to inadvertent compression.</p><p>The <strong>"IOL Model"</strong> field allows the calculator to determine which correction factor should be applied to modify the AL value used in the Theoretical Internal Lens Position (TILP) prediction algorithm. It also allows to perform a calculation for a meniscus IOL (e.g. Alcon MA60). In this case, the A constant can't be modified, which is normal : the geometrical particularities of this meniscus lens is integrated into the formula and customized for each lens power along the -5 to +5 range.</p><p>You may notice that the "<strong>LT", "CCT" and "WTW" fields</strong> are optional, this is to allow for calculation even with biometers which do not perform those measurements. However, it is strongly recommended to fill them in if they are available, especially the LT and CCT.</p><h3 id="complex-eyes">Complex Eyes</h3><p>The "Complex Eyes" platform allows the user to perform a calculation for eyes which have undergone previous corneal refractive surgery, have corneal scarring, or an Implantable Collamer Lens (ICL) in-situ.</p><p>If the "<strong>Corneal Refractive Surgery"</strong> field is activated, it is mandatory to determine which LVC procedure (myopic or hyperopic) was performed, as the corneal power adjustment will not be the same. Selecting this option also automatically removes the ARC and CCT from the lens position predictors.</p><p>If available, the <strong>"Pre-LVC anterior corneal radius"</strong> or the <strong>"Corrected ametropia"</strong> fields will allow to respectively obtain or approximate the pre-LVC anterior corneal radius, which will then be used to predict the posterior corneal radius more accurately.</p><p>If the <strong>"Radial Keratotomy"</strong> option is activated, a different method is used to predict the corneal power. ARC and CCT are also removed from the lens position predictors in this scenario. This option can be combined with the Corneal Refractive Surgery fields in order to perform a calculation for eyes with a history of RK and subsequent PRK. </p><p>The <strong>"Scarred cornea"</strong> option allows removal of the anterior corneal radius from the lens position predictors. It is wise to use it in those cases, to avoid adding a lens position prediction error to the corneal power prediction error which will probably occur.</p><p>The<strong> "ICL"</strong> option is used to apply a small correction to the AL measurement, due to the AL measurement error induced by those lenses, as described in the literature (5,6).</p><h3 id="second-eyes">Second Eyes</h3><p>In this section, it is possible for the user to fill in the information allowing to back-calculate the lens position of the first operated eye, which will then be used as a predictor to enhance the lens position prediction of the second eye.</p><p>The options here are the same as the options available for the 'Complex Eyes' section. </p><p>Note that the Second Eye function is not available for eyes with a discordant corneal surgery history, neither for MA60 IOLs.</p><h2 id="references-">References : </h2><p>1. Rose LT, Moshegov CN. Comparison of the Zeiss IOLMaster and applanation A-scan ultrasound: biometry for intraocular lens calculation. Clin Exp Ophthalmol. 2003 Apr;31(2):121-4.</p><p>2. Bjeloš Rončević M, Bušić M, Cima I, Kuzmanović Elabjer B, Bosnar D, Miletić D. Comparison of optical low-coherence reflectometry and applanation ultrasound biometry on intraocular lens power calculation. Graefes Arch Clin Exp Ophthalmol. 2011 Jan;249(1):69-75.</p><p>3. Rajan MS, Keilhorn I, Bell JA. Partial coherence laser interferometry vs conventional ultrasound biometry in intraocular lens power calculations. Eye (Lond). 2002 Sep;16(5):552-6.</p><p>4. Pereira A, Popovic M, Lloyd JC, El-Defrawy S, Schlenker MB. Preoperative measurements for cataract surgery: a comparison of ultrasound and optical biometric devices. Int Ophthalmol. 2021 Apr;41(4):1521-1530.</p><p>5. Chen X, Zhang D, Liu Z, Liu Y, Cai H, Wu Q, Zhang Y. Effect of Implantable Collamer Lens on Anterior Segment Measurement and Intraocular Lens Power Calculation Based on IOLMaster 700 and Sirius. J Ophthalmol. 2021 Dec</p><p>6. Pitault G, Leboeuf C, Leroux les Jardins S, Auclin F, Chong-Sit D, Baudouin C. Biométrie optique des yeux avec implants phaques [Optical biometry of eyes corrected by phakic intraocular lenses]. J Fr Ophtalmol. 2005 Dec;28(10):1052-7.</p>]]></content:encoded></item></channel></rss>