人工晶体度数计算公式推导(1)

Golden Grape

2019-10-13 12:00

按照Thomas Olsen的综述: Calculation of intraocular lens power: a review 进行推导, Thomas Olsen就是Olsen公式的作者, 算是很新的公式了. 这是他在2007年发的一篇综述.

讲人工晶体度数计算的文章, 大概一开始都要八卦一下第一枚人工晶体植入的故事, Ridley爵士在1949年进行了第一枚人工晶体植入, 但是度数算错了, -20D! Ridley用了Gullstrand模型眼的曲率参数, 但忘了PMMA的高折射率. (Sir Harold Ridley and his Fight for Sight )

基础光学¶

四则运算、不超过初中数学, 不必紧张

屈光度¶

$$ F=\frac{n_2-n_1}{r} \tag 1 $$

一个界面的屈光度F, 定义如上单位为D, 也是$M^{-1}$, 单位要千万小心, 如果一直都是字母没什么关系, 但如果代入具体的数字时, 要小心区分毫米、米. 其中:

$n_1$ 是光线入射到界面之前, 介质的折射率.
$n_2$ 是光线从界面出射之后, 介质的折射率.
r 是界面的曲率半径. 如果是向入射光方向凸起的, r的符号为正, 否则为负.

聚散度¶

$$ V=\frac{n}{d} \tag 2 $$

V: 光线的聚散度
d: 从当前面到聚焦点的距离, 注意也是有正负的.
n: 介质的折射率

如果聚散度为$V_1$的光线, 经过屈光度为F的透镜面, 出射光的聚散度$V_2$为 $$ V_2=V_1+F \tag 3 $$

补充:¶

对于平行光, 聚焦于无穷远点, 也就是d=∞, 那么V=0

如果聚散度为$V_0=\frac{n}{d}$的光, 在原来的介质中, 前进了$d_0$距离, 那么还剩$d-d_0$的距离就聚焦了, 所以这时候的聚散度$V=\frac{n}{d-d_0}$, 如果已知$V_0, n$, 那么: $$ d=\frac{n}{V_0} \tag {3.1} $$

将公式(3.1)代入 $$ V=\frac{n}{d-d_0} \tag {3.2} $$

得到: $$ V =\frac{n}{\frac{n}{V_0}-d_0} $$

上下同时除以n化简, 得到 $$ V =\frac{1}{\frac{1}{V_0}-\frac{d_0}{n}} \tag {3.3} $$

作为“薄透镜”的人工晶体¶

如果平行光通过曲率为K的角膜, 眼轴长为Ax, 有效晶体位置为d, 眼内介质的折射率为n, 那么要达成正视眼, IOL的度数应当是多少?

原文中的推导太快, 我稍微慢一点推:

平行光的屈光度=0, 经过曲率为K的角膜, 按照公式(3), $$ V_{离开角膜}=0+K=K \tag {4.1} $$

光线经过距离d, 到达人工晶体的位置, 这是的聚散度$V_1$的计算, 和刚才补充公式(3.3)的过程一致, 将公式(4.1)代入到公式(3.3)中, 得到:

$$ V_{进入人工晶体}=\frac{1}{\frac{1}{V_{离开角膜}}-\frac{d}{n_{房水}}} $$$$ V_{进入人工晶体}=\frac{1}{\frac{1}{K}-\frac{d}{n_{房水}}} \tag 4 $$

此处假设人工晶体的厚度忽略不计, 那么从人工晶体出来, 只要再走(眼轴长Ax - 人工晶体位置d)的距离就到达了视网膜. 为了能够聚焦到视网膜上, 离开人工晶体的光线, 聚散度应当满足如下公式:

$$ V_{离开人工晶体}=\frac{n_{玻璃体}}{Ax-d} \tag 5 $$

根据公式(3), 进入人工晶体和离开人工晶体的屈光度变化, 应当由人工晶体的度数$P_0$提供: $$ V_{离开人工晶体}=V_{进入人工晶体}+P_0 \tag 3 $$

将公式(4)和(5)代入:

$$ \frac{n_{玻璃体}}{Ax-d}=\frac{1}{\frac{1}{K}-\frac{d}{n_{房水}}}+P_0 $$$$ P_0=\frac{n_{玻璃体}}{Ax-d} - \frac{1}{\frac{1}{K}-\frac{d}{n_{房水}}} \tag 6 $$

这就是最简单的薄透镜公式. 里面$n_{玻璃体}, n_{房水}$ 是常数, Ax, K是可以测量到的, d在手术前其实是不知道的, 但可以通过各种方式来预测.

要在薄透镜公式进化,

将角膜和晶体作为有一定厚度的透镜, 前后表面还有着不同的曲率半径.
不仅仅考虑近轴光线, 还要考虑上周边光线造成的像差.

对ELP的预测:¶

有效人工晶体位置, ELP, 上面薄透镜中的d. 在手术之前是不知道的. 也是最关键的误差产生因素, 所以各种公式都在尽量预测ELP, 某种程度上也就是术后的前房深度ACD.

Olsen在该论文里吐槽了各种角膜测量、眼轴测量、超声vs激光测量之后, 列出了一个不错的表格, 说明了各种公式中用了哪些参数来预测ELP(ACD), 并且有各个公式的文献, 我就直接把sci-hub的全文链接附上了

ACD predictor	Formula ⁄ author
Axial length	Binkhorst	Binkhorst 1979
	Hoffer	Hoffer 199380338-0)
	SRK T	Sanders et al. 199080706-7)
	Holladay	Holladay et al. 198880059-2)
	Haigis	Haigis 2004 (原文似乎是到Shammas的书里找)
	Olsen	Olsen 1987a80128-1), 2006;Olsen & Gimbel 199380337-9); Olsen et al. 199580140-x)
Corneal height	Fyodorov	Fyodorov et al. 1975
	SRK T	Sanders et al. 199080706-7)
	Holladay I + II	Holladay et al. 198880059-2)
	Hoffer	Hoffer 199380338-0)
	Olsen	Olsen 1986b80099-2), 2006
Preoperative ACD	Haigis	Haigis 2004
	Holladay II, Olsen	Olsen 1986b80099-2), 2006
Lens thickness	Holladay II, Olsen	Olsen et al. 199580140-x), Olsen 2006
Others
Age	Holladay II	Unpublished data
	Olsen	Olsen 2006
Refraction	Holladay II, Olsen	Olsen 2006

先贤们的文献很多, 慢慢看吧.