This repository contains raw data and analysis routines of the publication “Sub-diffusion flow velocimetry with number fluctuation optical coherence tomography” in Optics Express (doi.org/10.1364/OE.474279). The reader is free to use the scripts and data in this depository if the manuscript is correctly cited in their work. For further questions, feel free to contact the corresponding author. Python 3.7 was used for programming. Keep in mind that running files with larger time series length may take up to 10 minutes and 2D flow profile analysis may take up to one hour. For 1D depth-resolved measurements each dataset includes diffusion, focus (beam shape) calibration, and flow measurements for different discharge rates, Q. For all measurements time series length is 31000 points and the sampling rate is 5.5 kHz. Diffusion measurements are performed on a static sample with a stationary beam. Focus (waist) calibration measurements are performed by moving the OCT beam over the static sample with a known velocity. Flow measurements are performed on the flowing sample with the stationary beam. Each measurement is averaged 6 times. The analysis process is as follows: Firstly, the beam waist (focus) calibration is performed using the script ‘Beam Shape.py’. For improved accuracy it is preferable to perform several measurements and average beam waist values at every depth. Secondly, the Doppler angle is determined using a flow measurement with the largest discharge rate using the script ‘Doppler Angle.py’. Thirdly, the flow profiles are obtained with predetermined calibration parameters using the script ‘Flow Profile.py’. Finally, the particle number density is calculated using the script ‘Number Density.py’. This requires knowledge of particle size for calculating the theoretical number density values. The particle size can be determined using the script ‘Diffusion.py’. All file names are sufficiently descriptive, showing whether it is diffusion, focus (waist) calibration or flow measurement. For 2D depth and laterally resolved measurements each dataset includes diffusion, focus (beam shape) calibration, M-scan and B-scan flow measurements for different discharge rates, Q. Diffusion and focus calibration measurements are same as in 1D. M-scan flow measurements are performed on a flowing sample with a stationary beam. They are same as flow measurements in 1D and are only used for determining the Doppler angle. B-scan flow measurements are performed by moving the OCT beam over the flowing sample with a known velocity. 2D flow profiles can be determined using the script ‘2D Flow Profile.py’. The table below summarizes all datasets and Python scripts uploaded to this repository. Name Usability Description Dataset, 15-03-2022.zip 1D measurements Dataset for Doppler angle of 0.34 deg and alignment angle of 0 deg. Dataset, 16-03-2022.zip 1D measurements Dataset for Doppler angle of 1.74 deg and alignment angle of 2.3 deg. Dataset, 22-03-2022.zip 1D measurements Dataset for Doppler angle of 1.00 deg and alignment angle of 1.15 deg. Dataset, 08-07-2022.zip 2D measurements Dataset for Doppler angle of 1.84 deg and alignment angle of 0 deg. Chirp.data All measurements File containing k-interpolation data ReadOCTFile.py All measurements Written by Jos de Wit, this module reads and imports spectra from raw OCT files. Processing.py All measurements This module contains all analysis and processing routines. Diffusion.py All measurements This script determines particle size from raw OCT spectra. Beam Shape.py All measurements This script determines axial beam shape from raw OCT spectra. Doppler Angle.py All measurements This script determines Doppler angle from raw OCT spectra. Flow Profile.py 1D measurements This script determines flow profiles from raw OCT spectra. Number Density.py 1D measurements This script determines particle number density raw OCT spectra. 2D Flow Profile.py 2D measurements This script determines 2D flow profiles from raw OCT spectra.