Utilizator:MagnInd/RFPP

Pupillary light reflex

Pupillary Light Reflex is modeled as a physiologically-based non-linear delay differential equation that describes the changes in the pupil diameter as a function of the environment lighting:

$M(D)=atanh\left({\frac {D-4.9}{3}}\right)$

${\frac {dM}{dD}}{\frac {dD}{dt}}+2.3026\;atanh\left({\frac {D-4.9}{3}}\right)=5.2-0.45\;ln\left[{\frac {\Phi (t-\tau )}{4.8118~\times ~10^{-10}}}\right]\;\;$

where $D$ is the pupil diameter measured in millimeters and $\Phi (t-\tau )$ is the luminous intensity reaching the retina in a time $t$, which can be described as $\Phi =IA$: luminance reaching the eye in lumens/mm² times the pupil area in mm². $\tau$ is the pupillary latency, a time delay between the instant in which the light pulse reaches the retina and the beginning of iridal reaction due nerve transmission, neuro-muscular excitation and activation delays. $dM$, $dD$ and $dt$ are the derivatives for the $M$ function, pupil diameter $D$ and time $t$.

Since the pupil constriction velocity is approximately 3 times faster than (re)dilation velocity , different step sizes in the numerical solver simulation must be used:

$dt_{c}={\frac {T_{c}-T_{p}}{S}}$

$dt_{d}={\frac {T_{c}-T_{p}}{3S}}$

where $dt_{c}$ and $dt_{d}$ are respectively the $dt$ for constriction and dilation measured in milliseconds, $Tc$ and $Tp$ are respectively the current and previous simulation times (times since the simulation started) measured in milliseconds, $S$ is a constant that affects the constriction/dilation velocity and varies among individuals. The higher the $S$ value, the smaller the time step used in the simulation and, consequently, the smaller the pupil constriction/dilation velocity.

In order to improve the realism of the resulting simulations, the hippus effect can be approximated by adding small random variations to the environment light (in the range of 0.05 Hz to 0.3 Hz) as proposed by .