A method to correct for surface scattering in spatial frequency domain imaging (SFDI) is presented. The use of a modified analytical solution of the radiative transfer equation allows calculation of the reflectance and the phase of a rough semi-infinite geometry so that both spatial frequency domain reflectance and phase can be applied for precise retrieval of the bulk optical properties and the surface scattering. For validation of the method, phantoms with different surface roughness were produced. Contrarily, with the modified theory, it was possible to dramatically reduce systematic errors due to surface scattering. The evaluation of these measurements with the state-of-the-art theory and measuring modality, i.e., using crossed linear polarizers, reveals large errors in the determined optical properties, depending on the surface roughness, of up to ≈100 % . These results were confirmed with SFDI measurements on a phantom that has a structured rough surface.

Spatial frequency domain imaging (SFDI) is well established in biology and medicine for non-contact, wide-field imaging of optical properties and 3D topography. Especially for turbid media with displaced, tilted or irregularly shaped surfaces, the reliable quantitative measurement of diffuse reflectance requires efficient calibration and correction methods. In this work, we present the implementation of a generic and hardware independent calibration routine for SFDI setups based on the so-called pinhole camera model for both projection and detection. Using a two-step geometric and intensity calibration, we obtain an imaging model that efficiently and accurately determines 3D topography and diffuse reflectance for subsequently measured samples, taking into account their relative distance and orientation to the camera and projector, as well as the distortions of the optical system. Derived correction procedures for position- and orientation-dependent changes in spatial frequency and intensity allow the determination of the effective scattering coefficient ${\mathsf{\mu }}_s^{\prime }$ and the absorption coefficient ${\mathsf{\mu }}_a$ when measuring a spherical optical phantom at three different measurement positions and at nine wavelengths with an average error of 5% and 12%, respectively. Model-based calibration allows the characterization of the imaging properties of the entire SFDI system without prior knowledge, enabling the future development of a digital twin for synthetic data generation or more robust evaluation methods.