Editing STIR-ASIM Howto (section)

==STIR run details==

Use of STIR   FBP2D and FBP3DRP produce the expected images. However OSMAPOSL took a long time to run and produced no output for either the 32 ring or 128 ring configurations.

A willing volunteer is required to pursue use of OSMAPOSL.

'''Example FBP2D.par'''
<pre>
fbp2dparameters :=
input file := asim_nrmd_noisy_sino_5.hs
output filename prefix := img_ias1_asim_FBP2D_5
<nowiki>; </nowiki> output image parameters
<nowiki>; </nowiki> zoom defaults to 1
zoom := 1
<nowiki>; </nowiki> image size defaults to whole FOV
xy output image size (in pixels) := -1
<nowiki>; </nowiki> default means: call SSRB only if no axial compression is already present 
num segments to combine with ssrb := 5
<nowiki>; </nowiki> filter parameters, default to pure ramp
alpha parameter for ramp filter := 0.5
cut-off for ramp filter (in cycles) := 0.5
end :=
</pre>


'''Example FBP3DRP.par'''
<pre>
fbp3drpparameters :=
input file := asim_nrmd_noisy_sino_5.hs
output filename prefix := img_ias1_asim_FBP3DRP_5
<nowiki>; </nowiki> all following parameters are optional
<nowiki>;;;;;;; </nowiki> output image parameters
zoom := 1
<nowiki>; </nowiki> defaults to cover whole FOV
xy output image size (in pixels) := 180
maximum absolute segment number to process := 2
<nowiki>;;;;;;; </nowiki> parameters for initial image you can use an existing image (but be careful about the scale!)
<nowiki>; </nowiki> image to be used for reprojection := some_image
<nowiki>; </nowiki> or you can use FBP on 2D data. following are reconstruction parameters for 2D
<nowiki>; </nowiki> default means: call SSRB only if no axial compression is already present
num segments to combine with ssrb := -1

<nowiki>; </nowiki> filter parameters for initial FBP2D, default to pure ramp
alpha parameter for ramp filter := 1
cut-off for ramp filter (in cycles) := 0.5

<nowiki>; </nowiki> parameters for Colsher filter
alpha parameter for colsher filter in axial direction := 1
cut-off for colsher filter in axial direction (in cycles) := 0.5
alpha parameter for colsher filter in planar direction := 1
cut-off for colsher filter in planar direction (in cycles) := 0.5
<nowiki>; </nowiki> define colsher on finer grid. The higher the number, the better (but slower)
stretch factor for colsher filter definition in axial direction := 2
stretch factor for colsher filter definition in planar direction := 2

<nowiki>;;;;;;;;; </nowiki> other parameters
save intermediate images := 0
display level := 0
end := 
</pre>


===AKNOWLEDGEMENTS===

This work involved assistance from

J.F.,
Software Support,
Centre for Vision, Speech and Signal Processing,
University of Surrey,
U.K.

Brian Elston, 
Senior Computer Specialist, 
University of Washington,
Department of Radiology,
U.S.A.

Please do not contact them directly.
Rather, use the STIR user list and the ASIM website managed by UWMC IRL;
University of Washington.