]> Import two-point file

Import two-point file

Read the *.2pt file generated by either tMCimg or tFDimg. The two-point file records the photon density throughout the medium (possibly including its phase or time course, depending on the details of the simulation).

Function Summary

Syntax:	[data,seg] = loadMC2pt(filename);
Inputs:	filename	Filename to read (without the .2pt extension)
Outputs:	data	Photon density data
Outputs:	seg	Contents of the segmentation file

Detailed Descriptions

To generate a forward forward matrix from Monte Carlo two-point files, multiply the source photon density (or their gradients) by the detector sensitivity for each measurement. Since the reciprocity theorem allows us to treat detectors as time-reversed sources, the "detector" sensitivities are equivalent to source photon densities for a source located at the same position.

A (r_{src}, r_{vox}, r_{det}) \sim G (r_{vox} - r_{src}) G (r_{det} - r_{vox})

where the Green's function $G$ and the photon density returned by loadMC2pt() are the same to within a constant scale factor.

loadMC2pt needs an accompanying .inp file in order to get the system size and other simulation parameters. If this file is not found, the function will return an error.