PAROS 1.0.2__tar.gz → 2.0.0__tar.gz

{paros-1.0.2 → paros-2.0.0}/.github/workflows/ci.yml

@@ -8,7 +8,7 @@ jobs:
     runs-on: ubuntu-latest
     strategy:
       matrix:
-        python: &python-versions ["3.9", "3.10", "3.11", "3.12"]
+        python: &python-versions ["3.10", "3.11", "3.12"]
     steps:
       - uses: actions/checkout@v6
       - name: Set up uv

paros-2.0.0/PAROS/data/calibration/Topcon/TRC_50DX.csv

@@ -0,0 +1,10 @@
+experiment_number,lens_1,lens_2,image_size_px,distance_camera_eye_mm,distance_camera_pinhole_mm,distance_pinhole_lens1_mm,distance_lens1_lens2_mm,distance_lens2_grid_mm,axial_length
+0,LB_1471,LB_1945,696.6666667,55,49.61969784,2.880302163,3.027080745,29.69261709,40.71969784
+1,LB_1471,LB_1945,672.6666667,55,49.61969784,2.880302163,3.027080745,31.69261709,42.71969784
+2,LB_1471,LB_1945,660.0,55,49.61969784,2.880302163,3.027080745,33.69261709,44.71969784
+3,LB_1471,LB_1945,637.0,55,49.61969784,2.880302163,3.027080745,35.69261709,46.71969784
+4,LB_1471,LB_1945,625.0,55,49.61969784,2.880302163,3.027080745,37.69261709,48.71969784
+5,LB_1471,LB_1945,608.6666667,55,49.61969784,2.880302163,3.027080745,39.69261709,50.71969784
+6,LB_1471,LB_1945,597.3333333,55,49.61969784,2.880302163,3.027080745,41.69261709,52.71969784
+7,LB_1471,LB_1945,574.6666667,55,49.61969784,2.880302163,3.027080745,43.69261709,54.71969784
+8,LB_1471,LB_1945,561.3333333,55,49.61969784,2.880302163,3.027080745,45.39261709,56.41969784

paros-2.0.0/PAROS/data/calibration/Topcon/TRC_50IX.csv

@@ -0,0 +1,11 @@
+lens_1,lens_2,axial_length,distance_camera_eye_mm,image_size_px,distance_camera_pinhole_mm,distance_pinhole_lens1_mm,distance_lens1_lens2_mm,distance_lens2_grid_mm,experiment_number
+LE_1156,LB_1945,83.9874534,55,331.0,50.2374534,2.262546596,6.344836312,71.54261709,0
+LE_1156,LB_1945,90.7374534,55,313.0,50.2374534,2.262546596,6.344836312,78.29261709,1
+LE_1156,LB_1945,78.0374534,55,346.0,50.2374534,2.262546596,6.344836312,65.59261709,2
+LE_1156,LB_1945,134.6374534,55,224.0,50.2374534,2.262546596,6.344836312,122.1926171,3
+LE_1156,LB_1945,107.9374534,55,269.0,50.2374534,2.262546596,6.344836312,95.49261709,4
+LE_1156,LB_1945,72.0374534,55,371.0,50.2374534,2.262546596,6.344836312,59.59261709,5
+LE_1156,LB_1945,63.6374534,55,408.0,50.2374534,2.262546596,6.344836312,51.19261709,6
+LE_1234,LB_1056,62.51915876,50,428.0,45.31915876,2.180841236,5.205295461,51.1138633,7
+LE_1234,LB_1056,75.91915876,50,366.0,45.31915876,2.180841236,5.205295461,64.5138633,8
+LE_1234,LB_1056,96.51915876,50,298.0,45.31915876,2.180841236,5.205295461,85.1138633,9

{paros-1.0.2 → paros-2.0.0}/PAROS/eyephantom.py

@@ -15,8 +15,6 @@ import sympy as sp
 
 from PAROS.fundscale import spherical_interface, uniform_medium
 
-sp.init_printing(use_unicode=True)
-
 
 class Lens:
     """Ray transfer model of a lens."""

{paros-1.0.2 → paros-2.0.0}/PAROS/fundscale.py

@@ -24,20 +24,32 @@ from __future__ import annotations
 
 import copy
 import math
-from typing import Literal, NamedTuple, Union
+from sys import version_info
+from typing import Literal, NamedTuple, Union, cast
 from warnings import warn
 
 import numpy as np
 import sympy as sp
 
-sp.init_printing(use_unicode=True)
+if version_info <= (3, 11):
+    from typing_extensions import NotRequired, TypedDict, Unpack
+else:
+    from typing import NotRequired, TypedDict, Unpack
 
-NumberOrSymbol = Union[int, float, sp.Symbol]
-EyeModelType = Literal["Navarro", "VughtIOL"]
 
+__all__ = [
+    "Eye",
+    "Camera",
+    "calculate_magnification",
+    "spherical_interface",
+    "uniform_medium",
+    "medium_change",
+    "calculate_piol_curvature",
+    "calculate_piol_matrix",
+]
 
-def _round_expr(expr: sp.Expr, num_digits: int = 3) -> sp.Expr:
-    return expr.xreplace({n: round(n, num_digits) for n in expr.atoms(sp.Number)})
+NumberOrSymbol = Union[int, float, sp.Symbol]
+EyeModelType = Literal["Navarro", "VughtIOL"]
 
 
 def uniform_medium(thickness: float) -> sp.Matrix:
@@ -177,6 +189,67 @@ def calculate_piol_matrix(
     return piol_matrix * distance_matrix
 
 
+class EyeGeometry(TypedDict, total=True):
+    """Eye geometry parameters.
+
+    All curvatures and distances are specified in meters and the spherical equivalent in diopters.
+    PAROS uses an inverted eye model, so for a normal eye the radii of curvature for the cornea and
+    lens front surfaces are negative.
+
+    Attributes
+    ----------
+    R_corF : float
+        Cornea front radius.
+    R_corB : float
+        Cornea back radius.
+    R_lensF : float
+        Lens front radius.
+    R_lensB : float
+        Lens back radius.
+    D_cor : float
+        Cornea thickness.
+    D_ACD : float
+        Anterior chamber depth.
+    D_lens : float
+        Lens thickness.
+    D_vitr : float
+        Vitreous thickness.
+    SE : float
+        Spherical equivalent of refraction of the eye model.
+    """
+
+    R_corF: float
+    R_corB: float
+    R_lensF: float
+    R_lensB: float
+    D_cor: float
+    D_ACD: float
+    D_lens: float
+    D_vitr: float
+    SE: NotRequired[float]
+
+
+class _PartialEyeGeometry(EyeGeometry, total=False):
+    """Eye geometry parameters.
+
+    This class can be used to define a partial eye geometry, in which only a subset
+    of the parameters are specified.
+
+    See Also
+    --------
+    EyeGeometry : Full eye geometry parameters.
+    """
+
+    R_corF: float
+    R_corB: float
+    R_lensF: float
+    R_lensB: float
+    D_cor: float
+    D_ACD: float
+    D_lens: float
+    D_vitr: float
+
+
 class PhakicIOL(NamedTuple):
     """Parameters of a phakic IOL.
 
@@ -198,15 +271,93 @@ class PhakicIOL(NamedTuple):
     lens_distance: NumberOrSymbol
 
 
+class RefractiveIndices(TypedDict):
+    """Refractive indices of the eye media.
+
+    Attributes
+    ----------
+    cor : float
+        Refractive index of the cornea.
+    aq : float
+        Refractive index of the aqueous humor.
+    lens : float
+        Refractive index of the crystalline lens.
+    vit : float
+        Refractive index of the vitreous humor.
+    """
+
+    cor: float
+    aq: float
+    lens: float
+    vit: float
+
+
+class _PartialRefractiveIndices(RefractiveIndices, total=False):
+    """Refractive indices of the eye media.
+
+    This class can be used to define a partial set of refractive indices, in which only a
+    subset of the indices are specified.
+
+    See Also
+    --------
+    RefractiveIndices : Full set of refractive indices.
+    """
+
+    cor: float
+    aq: float
+    lens: float
+    vit: float
+
+
+_DEFAULT_GEOMETRIES: dict[EyeModelType, EyeGeometry] = {
+    "Navarro": EyeGeometry(
+        R_corF=-7.72e-3,
+        R_corB=-6.50e-3,
+        R_lensF=-10.20e-3,
+        R_lensB=+6.00e-3,
+        D_cor=0.55e-3,
+        D_ACD=3.05e-3,
+        D_lens=4.00e-3,
+        D_vitr=16.3203e-3,
+    ),
+    "VughtIOL": EyeGeometry(
+        R_corF=-7.72e-3,
+        R_corB=-6.50e-3,
+        R_lensF=-8.16e-3,
+        R_lensB=+11.18e-3,
+        D_cor=0.55e-3,
+        D_ACD=3.05e-3,
+        D_lens=0.6896e-3,
+        D_vitr=19.3203e-3,
+    ),
+}
+
+_DEFAULT_REFRACTIVE_INDICES: dict[EyeModelType, RefractiveIndices] = {
+    "Navarro": RefractiveIndices(
+        cor=1.3777,
+        aq=1.3391,
+        lens=1.4222,
+        vit=1.3377,
+    ),
+    "VughtIOL": RefractiveIndices(
+        cor=1.3777,
+        aq=1.3391,
+        lens=1.47,
+        vit=1.3377,
+    ),
+}
+
+
 class Eye:
     """Symbolic representation of an eye model."""
 
     def __init__(
         self,
         name: str = "testEye",
-        geometry: dict[str, float] | None = None,
+        geometry: _PartialEyeGeometry | EyeGeometry | None = None,
         model_type: EyeModelType = "Navarro",
         NType: EyeModelType = "Navarro",  # noqa: N803
+        refractive_indices: _PartialRefractiveIndices | None = None,
         refraction: float | None = None,
         pIOL: PhakicIOL | None = None,  # noqa: N803
     ) -> None:
@@ -222,7 +373,7 @@ class Eye:
         ----------
         name : str
             Name of the eye model.
-        geometry : dict[str, float]
+        geometry : _PartialEyeGeometry | None
             Dictionary with numerical geometric parameters of the eye. If `None`, the
             geometry is based on `model_type`. Can be a partial dictionary, in which
             case missing values are based on the eye model specified by `model_type`.
@@ -235,6 +386,10 @@ class Eye:
         refraction : float
             Spherical equivalent of refraction of the eye model. The lens back curvature
             can be adjusted to obtain this refraction with `Eye.adjust_lens_back`.
+        refractive_indices : _PartialRefractiveIndices | None
+            Dictionary with numerical values of the refractive indices of the eye media.
+            If `None`, the refractive indices are based on `NType`. If specified, the
+            refractive indices corresponding to `NType` are updated with the values in this dictionary.
         pIOL : PhakicIOL
             Properties of the pIOL, if present. Tuple of power, thickness, refractive
             index, distance to lens.
@@ -242,28 +397,20 @@ class Eye:
         self.name = name
         self.spherical_equivalent = refraction
         self.is_pseudophakic = model_type == "VughtIOL"
-        self.axial_length = None
         self.pIOL = pIOL
 
         self.geometry = self._build_geometry_dictionary(model_type, geometry)
 
-        if NType == "Navarro":
-            self.refractive_indices = {
-                "cor": 1.3777,
-                "aq": 1.3391,
-                "lens": 1.4222,
-                "vit": 1.3377,
-            }
-        elif NType == "VughtIOL":
-            self.refractive_indices = {
-                "cor": 1.3777,
-                "aq": 1.3391,
-                "lens": 1.47,
-                "vit": 1.3377,
-            }
+        self.refractive_indices: RefractiveIndices
+
+        if NType in _DEFAULT_REFRACTIVE_INDICES:
+            self.refractive_indices = _DEFAULT_REFRACTIVE_INDICES[NType].copy()
         else:
             raise ValueError(f"Model type {NType} is undefined.")
 
+        if refractive_indices is not None:
+            self.refractive_indices.update(refractive_indices)
+
         self.R_corF, self.R_corB, self.R_lensF, self.R_lensB = sp.symbols(
             "R_corF R_corB R_lensF R_lensB"
         )
@@ -271,46 +418,30 @@ class Eye:
             "D_cor D_ACD D_lens D_vitr"
         )
 
-        self.axial_length = (
+        self._update_ray_transfer_matrix()
+
+    @property
+    def axial_length(self) -> float:
+        return (
             self.geometry["D_cor"]
             + self.geometry["D_ACD"]
             + self.geometry["D_lens"]
             + self.geometry["D_vitr"]
         )
-        self.ray_transfer_matrix = self.calculate_ray_transfer_matrix()
 
     def __str__(self) -> str:
         return f"{self.name}"
 
     @staticmethod
     def _build_geometry_dictionary(
-        model_type: EyeModelType, partial_geometry: dict[str, float] | None = None
-    ):
-        if model_type == "Navarro":
-            geometry = {
-                "R_corF": -7.72e-3,
-                "R_corB": -6.50e-3,
-                "R_lensF": -10.20e-3,
-                "R_lensB": +6.00e-3,
-                "D_cor": 0.55e-3,
-                "D_ACD": 3.05e-3,
-                "D_lens": 4.00e-3,
-                "D_vitr": 16.3203e-3,
-            }
-        elif model_type == "VughtIOL":
-            geometry = {
-                "R_corF": -7.72e-3,
-                "R_corB": -6.50e-3,
-                "R_lensF": -8.16e-3,
-                "R_lensB": +11.18e-3,
-                "D_cor": 0.55e-3,
-                "D_ACD": 3.05e-3,
-                "D_lens": 0.6896e-3,
-                "D_vitr": 19.3203e-3,
-            }
-        else:
+        model_type: EyeModelType,
+        partial_geometry: _PartialEyeGeometry | EyeGeometry | None = None,
+    ) -> EyeGeometry:
+        if model_type not in _DEFAULT_GEOMETRIES:
             raise ValueError(f"Model type {model_type} is undefined.")
 
+        geometry = _DEFAULT_GEOMETRIES[model_type].copy()
+
         if partial_geometry is not None:
             # Estimate the cornea back curvature if it is unspecified and the cornea
             # front curvature is specified
@@ -378,6 +509,34 @@ class Eye:
             * uniform_medium(self.D_vitr)
         )
 
+    def update_geometry(self, **geometry: Unpack[_PartialEyeGeometry]) -> None:
+        """Update the eye geometry with new values.
+
+        Parameters
+        ----------
+        geometry : _PartialEyeGeometry
+            Geometry parameters to update as keyword arguments.
+        """
+        self.geometry.update(geometry)
+        self._update_ray_transfer_matrix()
+
+    def update_refractive_indices(
+        self, **refractive_indices: Unpack[_PartialRefractiveIndices]
+    ) -> None:
+        """Update the eye's refractive indices with new values.
+
+        Parameters
+        ----------
+        refractive_indices : _PartialRefractiveIndices
+            Refractive indices to update as keyword arguments.
+        """
+        self.refractive_indices.update(refractive_indices)
+        self._update_ray_transfer_matrix()
+
+    def _update_ray_transfer_matrix(self) -> None:
+        """Update the eye's ray transfer matrix."""
+        self.ray_transfer_matrix = self.calculate_ray_transfer_matrix()
+
     def adjust_lens_back(
         self, target_refraction: float, *, update_model: bool = False
     ) -> tuple[float, float]:
@@ -491,8 +650,8 @@ class Eye:
     ) -> tuple[float, sp.Matrix, float]:
         """Calculate the total magnification of the eye - camera system.
 
-        A structure on the central retina is depicted `magnification` times larger
-        on the camera sensor.
+        The magnification is calculated in units of pixels per millimeter, i.e. a structure
+        of 1 mm on the retina has a size of `magnification` pixels on the camera sensor.
 
         Parameters
         ----------
@@ -505,7 +664,7 @@ class Eye:
         Returns
         -------
         magnification : float
-            Central magnification of the eye - camera system.
+            Central magnification of the eye - camera system, in pixels per millimeter.
         focused_system_matrix : sympy.Matrix
             Ray transfer matrix for the full system.
         focus_lens_radius : float
@@ -516,23 +675,47 @@ class Eye:
             eye_matrix, distance_eye_camera, return_focus_lens_curvature=True
         )
 
-        magnification = focused_system_matrix[0, 0]
+        magnification = focused_system_matrix[0, 0] * camera.pixel_density
 
         return magnification, focused_system_matrix, focus_lens_radius
 
 
+_MISSING = cast(float, object())
+
+
 class Camera:
     def __init__(
         self,
-        F_cond: NumberOrSymbol = None,  # noqa: N803
-        a1: NumberOrSymbol = None,
+        F_cond: NumberOrSymbol | None = None,  # noqa: N803
+        a1: NumberOrSymbol | None = None,
+        pixel_density: float = _MISSING,
         camera_type: Literal["default"] = "default",
     ) -> None:
+        """Create a new camera model.
+
+        Parameters
+        ----------
+        F_cond : NumberOrSymbol
+            Focal length of the condenser lens, in meters.
+        a1 : NumberOrSymbol
+            First order calibration term.
+        camera_type : Literal["default"]
+            Type of camera. Currently only "default" is implemented.
+        pixel_density : float
+            Pixel density of the camera sensor, in pixels per millimeter.
+            Defaults to 100 pixels/mm.
+        """
+        if pixel_density is _MISSING:
+            raise ValueError("Pixel density must be specified.")
+        if pixel_density <= 0:
+            raise ValueError("Pixel density must be a positive float.")
+
         self.camera_type = "lensTaylor"
         self.F_cond = sp.Symbol("F_cond", real=True)
         self.d_CCD = sp.symbols("d_CCD")
         self.R_foc = sp.Symbol("R_foc", real=True)  # in m
         self.a1 = sp.Symbol("a1", real=True)
+        self.pixel_density = pixel_density  # pixels per mm
 
         if camera_type == "default":
             self.d_CCD = self.F_cond
@@ -549,7 +732,10 @@ class Camera:
             self.n_glas, 1.0, -self.R_foc
         ) * spherical_interface(1.0, self.n_glas, self.R_foc)
         self.correction_term = sp.Matrix(
-            [[1 + self.a1 / self.R_foc, 0], [0, 1.0 / (1 + self.a1 / self.R_foc)]]
+            [
+                [1 + self.a1 / self.R_foc, 0],
+                [0, 1.0 / (1 + self.a1 / self.R_foc)],
+            ]
         )
         self.ray_transfer_matrix = (
             self.correction_term
@@ -558,13 +744,15 @@ class Camera:
             * self.condenser_lens
         )
 
-        self.M_camera_alg = self.ray_transfer_matrix
-
         if F_cond and (a1 is not None):
-            self.ray_transfer_matrix = self.M_camera_alg.evalf(
+            self.ray_transfer_matrix = self.ray_transfer_matrix.evalf(
                 subs={self.F_cond: F_cond, self.a1: a1}
             )
 
+    @property
+    def M_camera_alg(self):
+        return self.ray_transfer_matrix
+
     def calculate_focus_lens_radius(
         self, eye_matrix: sp.Matrix, distance_eye_camera: float = 0.05
     ) -> float:
@@ -616,7 +804,7 @@ class Camera:
         )
 
         # For contact cameras such as the Panoret fundus camera, the refractive index (of air) needs to be changed to
-        # that of the medium use dbetween the eye and camera. This can be done by changing the use of uniform_medium()
+        # that of the medium used between the eye and camera. This can be done by changing the use of uniform_medium()
         # to medium_change().
         focus_lens_curvature = self.calculate_focus_lens_radius(
             eye_matrix, distance_eye_camera
@@ -629,6 +817,21 @@ class Camera:
 
         return system_matrix.evalf(subs={self.R_foc: focus_lens_curvature})
 
+    def get_size_on_ccd(self, size_pixels: float) -> float:
+        """Convert an image size in pixels to its physical size on the CCD sensor in meters.
+
+        Parameters
+        ----------
+        size_pixels : float
+            Image size on the CCD sensor in pixels.
+
+        Returns
+        -------
+        float
+            Image size on the CCD sensor in meters.
+        """
+        return size_pixels / (self.pixel_density * 1000)  # convert mm to m
+
 
 # Maximum allowed difference between the calculated and specified refraction of the eye model
 _MAXIMUM_REFRACTION_DEVIATION = 0.05
@@ -647,7 +850,7 @@ def calculate_magnification(
 ) -> tuple[float, float, Eye]:
     """Calculate the total magnification of the eye - camera system.
 
-    A structure on the central retina is depicted `magnification` times larger
+    A structure of 1 mm  on the central retina has a size of `magnification` pixels
     on the camera sensor.
 
     Parameters
@@ -669,7 +872,7 @@ def calculate_magnification(
     Returns
     -------
     magnification : float
-        Central magnification of the eye - camera system.
+        Central magnification of the eye - camera system, in pixels per millimeter.
     glasses_power : float
         Power of the corrective lens required to correct for the eye's refraction, in
         diopters.
@@ -681,7 +884,7 @@ def calculate_magnification(
     When the calculated glasses power differs significantly from
     `eye.spherical_equivalent`, a warning is displayed.
     """
-    glasses_curvature, glasses_power = eye.calculate_refraction(display=False)
+    _, glasses_power = eye.calculate_refraction(display=False)
 
     if (
         not suppress_warnings
@@ -713,7 +916,7 @@ def calculate_magnification(
         subs=({camera.R_foc: focus_lens_radius})
     )
 
-    magnification: float = focused_system_matrix[0, 0]
+    magnification: float = focused_system_matrix[0, 0] * camera.pixel_density
 
     if abs(focused_system_matrix[0, 1]) > _MAXIMUM_FOCUS_DEVIATION:
         warn(f"focused_system_matrix not in focus for patient {eye.name}")

{paros-1.0.2 → paros-2.0.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PAROS
-Version: 1.0.2
+Version: 2.0.0
 Summary: PARaxial Optical fundus Scaling (PAROS)
 Author-email: "Lennart J. Pors" <l.j.pors@lumc.nl>, "Jan-Willem M. Beenakker" <j.w.m.beenakker@lumc.nl>, Corné Haasjes <c.haasjes@lumc.nl>
 License: MIT License
@@ -30,12 +30,13 @@ Classifier: Development Status :: 5 - Production/Stable
 Classifier: Programming Language :: Python :: 3
 Classifier: Topic :: Scientific/Engineering :: Physics
 Classifier: Typing :: Typed
-Requires-Python: >=3.9
+Requires-Python: >=3.10
 Requires-Dist: matplotlib
 Requires-Dist: numpy
 Requires-Dist: pandas
 Requires-Dist: scipy
 Requires-Dist: sympy
+Requires-Dist: typing-extensions>=4.12.2; python_full_version < '3.11'
 Description-Content-Type: text/markdown
 
 # PARaxial Optical fundus Scaling (PAROS)
@@ -49,6 +50,7 @@ The code is provided as is, without any warranty, under the [MIT license](LICENS
 ## Basic usage
 
 ### Try it online
+
 A basic version of PAROS can most easily be run online here: [basic version](https://demo.mreye.nl/paros)
 
 [![Try PAROS online](.github/assets/screenshot.png)](https://demo.mreye.nl/PAROS)
@@ -73,10 +75,15 @@ The implementation of PAROS in this repository is functional for the camera in o
 ## Camera constants
 
 Please find here a table of known camera calibration constants; this can be added upon by other contributors. 
-| Camera type     | CCD type | Condenser lens power (D) | First order calibration term |
-| --------------- | -------- | ------------------------ | ---------------------------- |
-| Topcon TRC-50DX |          | 37.6                     | 0.03481                      |
-| Topcon TRC-50IX |          | 37.7                     | 0.03226                      |
+
+| Camera type     | CCD type | Condenser lens power [D] | First order calibration term | Pixel density [px/mm] |
+| :-------------- | -------- | -----------------------: | ---------------------------: | --------------------: |
+| Topcon TRC-50DX |          | 37.6                     | 0.03481                      | 100                   |
+| Topcon TRC-50IX |          | 37.7                     | 0.03226                      | 100                   |
+
+> [!NOTE]
+> The pixel density for the Topcon system is estimated.
+> Since this value is only used to fit the other parameters, it does not need to be exact, but the same value should be used for both calibration and magnification calculation.
 
 ## Referencing
 

{paros-1.0.2 → paros-2.0.0}/README.md

@@ -9,6 +9,7 @@ The code is provided as is, without any warranty, under the [MIT license](LICENS
 ## Basic usage
 
 ### Try it online
+
 A basic version of PAROS can most easily be run online here: [basic version](https://demo.mreye.nl/paros)
 
 [![Try PAROS online](.github/assets/screenshot.png)](https://demo.mreye.nl/PAROS)
@@ -33,10 +34,15 @@ The implementation of PAROS in this repository is functional for the camera in o
 ## Camera constants
 
 Please find here a table of known camera calibration constants; this can be added upon by other contributors. 
-| Camera type     | CCD type | Condenser lens power (D) | First order calibration term |
-| --------------- | -------- | ------------------------ | ---------------------------- |
-| Topcon TRC-50DX |          | 37.6                     | 0.03481                      |
-| Topcon TRC-50IX |          | 37.7                     | 0.03226                      |
+
+| Camera type     | CCD type | Condenser lens power [D] | First order calibration term | Pixel density [px/mm] |
+| :-------------- | -------- | -----------------------: | ---------------------------: | --------------------: |
+| Topcon TRC-50DX |          | 37.6                     | 0.03481                      | 100                   |
+| Topcon TRC-50IX |          | 37.7                     | 0.03226                      | 100                   |
+
+> [!NOTE]
+> The pixel density for the Topcon system is estimated.
+> Since this value is only used to fit the other parameters, it does not need to be exact, but the same value should be used for both calibration and magnification calculation.
 
 ## Referencing