Space-time inverse-scattering of translation-based motion

This repository contains the MATLAB simulation and reconstruction code for space-time inverse scattering with translation-based motion, as used in the following publication:

Kim, Jeongsoo, and Shwetadwip Chowdhury. "Space-time inverse-scattering of translation-based motion." Optica 12.5 (2025): 643-653.

Abstract: In optical diffraction tomography (ODT), a sample's 3D refractive-index (RI) is often reconstructed after illuminating it from multiple angles, with the assumption that the sample remains static throughout data collection. When the sample undergoes dynamic motion during this data-collection process, significant artifacts and distortions compromise the fidelity of the reconstructed images. In this study, we develop a space-time inverse-scattering technique for ODT that compensates for the translational motion of multiple-scattering samples during data collection. Our approach involves formulating a joint optimization problem to simultaneously estimate a scattering sample's translational position at each measurement and its motion-corrected 3D RI distribution. Experimental results demonstrate the technique's effectiveness, yielding reconstructions with reduced artifacts, enhanced spatial resolution, and improved quantitative accuracy for samples undergoing continuous translational motion during imaging.

Experimental dataset

The experimental dataset for running the code can be downloaded from the link below. This dataset contains angular scattering measurements of a microsphere and scattering phantom undergoing translational motion between frames. It also includes the parameters of the angle-scanning imaging system.

Download Link: https://dataverse.tdl.org/dataset.xhtml?persistentId=doi:10.18738/T8/KYYDHX

Numerical Simulation

To verify the space-time inverse-scattering of translation-based motion through simulation, please follow the steps below.

1. Please download the Recon_corefunction folder, which contains the core functions required to run the reconstruction code.

2. The script Mainsimulation.m is used to run the simulation. Make sure to place this script in the same directory as the Reconcorefunction folder. Once everything is set up, you can run the script to reproduce the simulation results shown below.

3. Since the code contains comments explaining each section, it is recommended to execute the code step by step and read through the related comments when running the code for the first time.

Figure 1: Numerical simulation results

Experiment: Microsphere and Scattering phantom

To verify the space-time inverse-scattering of translation-based motion using a microsphere and scattering phantom, please follow the steps below.

1. Please download the Recon_corefunction folder, which contains the core functions required to run the reconstruction code.

2. Download the data from the link provided in the Experimental Dataset section and save it in the same directory as the Recon_corefunction folder.

3. The scripts MainMicrosphere.m and MainPhantom.m are used to reconstruct the experimental datasets. Place them in the same directory as the Recon_corefunction folder, then run the scripts in MATLAB to reproduce the reconstruction results shown below.

4. Since the code contains comments explaining each section, it is recommended to execute the code step by step and read through the related comments when running the code for the first time.

Figure 2: Experimental results using microsphere

Figure 3: Experimental results using scattering phantom