Total VariationΒΆ
You find the the example in
/scisoft/ESRF_sw/TEST_PYHST/TV_PHANTOM
Input File:
! the number of used projections is reduced increasing NRED
NRED=10
FILE_INTERVAL = 1*NRED
ANGLE_BETWEEN_PROJECTIONS = 0.120000*NRED
FILE_PREFIX = testsdata/proj
NUM_FIRST_IMAGE = 0
NUM_LAST_IMAGE = 1499
NUMBER_LENGTH_VARIES = NO
LENGTH_OF_NUMERICAL_PART = 4
FILE_POSTFIX = .edf
NUM_IMAGE_1 = 2048
NUM_IMAGE_2 = 1
IMAGE_PIXEL_SIZE_1 = 10000.000000
IMAGE_PIXEL_SIZE_2 = 10000.000000
SUBTRACT_BACKGROUND = NO ! No background subtraction
BACKGROUND_FILE = N.A.
CORRECT_FLATFIELD = NO
FLATFIELD_CHANGING = N.A.
FLATFIELD_FILE = N.A.
FF_PREFIX = N.A.
FF_NUM_FIRST_IMAGE = N.A.
FF_NUM_LAST_IMAGE = N.A.
FF_NUMBER_LENGTH_VARIES = N.A.
FF_LENGTH_OF_NUMERICAL_PART = N.A.
FF_POSTFIX = N.A.
FF_FILE_INTERVAL = N.A.
TAKE_LOGARITHM = NO
ROTATION_VERTICAL = YES
ROTATION_AXIS_POSITION = 1023.500000
OUTPUT_SINOGRAMS = NO ! Output sinograms to files or not
OUTPUT_RECONSTRUCTION = YES ! Reconstruct and save or not
START_VOXEL_1 = 1 ! X-start of reconstruction volume
START_VOXEL_2 = 1 ! Y-start of reconstruction volume
START_VOXEL_3 = 1 ! Z-start of reconstruction volume
END_VOXEL_1 = 2048 ! X-end of reconstruction volume
END_VOXEL_2 = 2048 ! Y-end of reconstruction volume
END_VOXEL_3 = 1 ! Z-end of reconstruction volume
OVERSAMPLING_FACTOR = 4 ! 0 = Linear, 1 = Nearest pixel
ANGLE_OFFSET = 0.000000 ! Reconstruction rotation offset angle in degrees
SINOGRAM_MEGABYTES = 400 ! Maximum size of sinogram storage (megabytes)
DO_CCD_FILTER = NO ! CCD filter (spikes)
CCD_FILTER = "CCD_Filter"
CCD_FILTER_PARA = {"threshold": 100.000000 }
DO_SINO_FILTER = NO ! Sinogram filter (rings)
SINO_FILTER = "SINO_Filter"
ar = Numeric.ones(1748,'f')
ar[0]=0.0
ar[2:18]=0.0
SINO_FILTER_PARA = {"FILTER": ar }
DO_AXIS_CORRECTION = NO ! Axis correction
AXIS_CORRECTION_FILE = 'filename'
OPTIONS= { 'padding':'E' , 'axis_to_the_center':'Y' } ! Padding and position axis
OUTPUT_FILE = slice_bp.vol
BETA_TV = 10
FISTA=1
ITERATIVE_CORRECTIONS=1000
ITERATIVE_CORRECTIONS_NOPREC=1
DENOISING_TYPE=1
Input File:
! the number of used projections is reduced increasing NRED
NRED=10
FILE_INTERVAL = 1*NRED
ANGLE_BETWEEN_PROJECTIONS = 0.120000*NRED
FILE_PREFIX = testsdata/proj
NUM_FIRST_IMAGE = 0
NUM_LAST_IMAGE = 1499
NUMBER_LENGTH_VARIES = NO
LENGTH_OF_NUMERICAL_PART = 4
FILE_POSTFIX = .edf
NUM_IMAGE_1 = 2048
NUM_IMAGE_2 = 1
IMAGE_PIXEL_SIZE_1 = 10000.000000
IMAGE_PIXEL_SIZE_2 = 10000.000000
SUBTRACT_BACKGROUND = NO ! No background subtraction
BACKGROUND_FILE = N.A.
CORRECT_FLATFIELD = NO
FLATFIELD_CHANGING = N.A.
FLATFIELD_FILE = N.A.
FF_PREFIX = N.A.
FF_NUM_FIRST_IMAGE = N.A.
FF_NUM_LAST_IMAGE = N.A.
FF_NUMBER_LENGTH_VARIES = N.A.
FF_LENGTH_OF_NUMERICAL_PART = N.A.
FF_POSTFIX = N.A.
FF_FILE_INTERVAL = N.A.
TAKE_LOGARITHM = NO
ROTATION_VERTICAL = YES
ROTATION_AXIS_POSITION = 1023.500000
OUTPUT_SINOGRAMS = NO ! Output sinograms to files or not
OUTPUT_RECONSTRUCTION = YES ! Reconstruct and save or not
START_VOXEL_1 = 1 ! X-start of reconstruction volume
START_VOXEL_2 = 1 ! Y-start of reconstruction volume
START_VOXEL_3 = 1 ! Z-start of reconstruction volume
END_VOXEL_1 = 2048 ! X-end of reconstruction volume
END_VOXEL_2 = 2048 ! Y-end of reconstruction volume
END_VOXEL_3 = 1 ! Z-end of reconstruction volume
OVERSAMPLING_FACTOR = 4 ! 0 = Linear, 1 = Nearest pixel
ANGLE_OFFSET = 0.000000 ! Reconstruction rotation offset angle in degrees
SINOGRAM_MEGABYTES = 400 ! Maximum size of sinogram storage (megabytes)
DO_CCD_FILTER = NO ! CCD filter (spikes)
CCD_FILTER = "CCD_Filter"
CCD_FILTER_PARA = {"threshold": 100.000000 }
DO_SINO_FILTER = NO ! Sinogram filter (rings)
SINO_FILTER = "SINO_Filter"
ar = Numeric.ones(1748,'f')
ar[0]=0.0
ar[2:18]=0.0
SINO_FILTER_PARA = {"FILTER": ar }
DO_AXIS_CORRECTION = NO ! Axis correction
AXIS_CORRECTION_FILE = 'filename'
OPTIONS= { 'padding':'E' , 'axis_to_the_center':'Y' } ! Padding and position axis
OUTPUT_FILE = slice_fbp.vol
BETA_TV = 0.01
FISTA=1
ITERATIVE_CORRECTIONS=100
ITERATIVE_CORRECTIONS_NOPREC=0
DENOISING_TYPE=1