frms,
flux,
int(noise),
filMode,
recoMode.lower(),
"-{}tv".format(tvGrad) if tvGrad else ""
)
fpath = os.path.join(VOLUME_OUTPUT_DIR, fname)
tmp_fpath = os.path.join(TMP_VOLUME_OUTPUT_DIR, fname)
# Create the workspace file.
algorithms = ''.join([
STRING_CBSIMULATION.format(
frames=frms,
i0=float(flux)/(65535.0*frms),
addPoissonNoise=int(noise),
sdd=GEOMETRIES[GEOMETRY].sdd,
spd=GEOMETRIES[GEOMETRY].spd
),
#STRING_PREFILTER if prefilter else "",
STRING_CBRECONSTRUCTION.format(
recoMode=int(recoMode == "FDK"),
fourierFilter=STRING_FOURIERFILTER if recoMode == "FDK" else "",
setSize=numRuns*"20 ",
numRuns=numRuns,
mult=1 if recoMode == "MLEM" else 0,
compute = 3 if prefilter else 2,
tvGrad=tvGrad
),
STRING_EXPORTRAW.format(
fname=tmp_fpath,
# Output file name.
fname = "{}frms_{:.0e}flux_{}noise_{}filter_{}{}".format(
frms,
flux,
int(noise),
filMode,
recoMode.lower(),
"-tv" if tvGrad else ""
)
fpath = os.path.join(OUTPUT_DIR, fname)
# Create the workspace file.
algorithms = ''.join([
STRING_CBSIMULATION.format(
frames=frms,
i0=float(flux)/(65535.0*frms),
addPoissonNoise=int(noise)
),
STRING_PREFILTER if prefilter else "",
STRING_CBRECONSTRUCTION.format(
recoMode=int(recoMode == "FDK"),
fourierFilter=STRING_FOURIERFILTER if recoMode == "FDK" else "",
setSize=numRuns*"20 ",
numRuns=numRuns,
mult=1 if recoMode == "MLEM" else 0,
compute = 3 if prefilter else 2,
tvGrad=tvGrad
),
#STRING_LINEEXTRACTION.format(
#fname=fpath,
#compute = 4 if prefilter else 3),
