Patch Fitting

http://arxiv.org/abs/1305.1256

You find the the example in

/scisoft/ESRF_sw/TEST_PYHST/MYPHANTOM
_images/slice.png

regularised image. 150 projections

_images/patches.png

the overcomplete basis

_images/slice_fbp.png

FBP reconstruction

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 = sino.edf
BACKGROUND_FILE = nienet
FLATFIELD_FILE = niente


NUM_FIRST_IMAGE = 0 
NUM_LAST_IMAGE = 1499
NUMBER_LENGTH_VARIES = NO
LENGTH_OF_NUMERICAL_PART = 4 
FILE_POSTFIX = .edf

NUM_IMAGE_1 = 1024 
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 = 511.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 =   1024+14 ! X-end of reconstruction volume
END_VOXEL_2 =   1024+14 ! 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 = pp_00.edf
MULTIFRAME=1
FBFILTER=1

BETA_TV = 0.001
ITERATIVE_CORRECTIONS_NOPREC=1
ITERATIVE_CORRECTIONS=6000
DENOISING_TYPE=4
PATCHES_FILE = patches_phantom.edf
FISTA=1
STEPFORPATCHES = 3
WEIGHT_OVERLAP=50.0