/
POSMAPOSL.par
90 lines (74 loc) · 3.15 KB
/
POSMAPOSL.par
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
OSMAPOSLParameters :=
; example file for running OSEM with interfiltering
; filter specific things are at the end
objective function type:=PoissonLogLikelihoodWithLinearKineticModelAndDynamicProjectionData
PoissonLogLikelihoodWithLinearKineticModelAndDynamicProjectionData Parameters:=
input file := SOME_DYNAMIC_PROJDATA.S
; if disabled, defaults to maximum segment number in the file
maximum absolute segment number to process := -1
; see User's Guide to see when you need this
zero end planes of segment 0:= 0
projector pair type := Matrix
Projector Pair Using Matrix Parameters :=
Matrix type := Ray Tracing
Ray tracing matrix parameters :=
End Ray tracing matrix parameters :=
End Projector Pair Using Matrix Parameters :=
; normalisation (and attenuation correction)
Bin Normalisation type := Chained
Chained Bin Normalisation Parameters:=
Bin Normalisation to apply first:= From ECAT7
Bin Normalisation From ECAT7:=
normalisation_ECAT7_filename:= NORMFILE.n
singles rates := Singles From ECAT7
Singles Rates From ECAT7:=
ECAT7_filename := SOME_DYNAMIC_PROJDATA.S
End Singles Rates From ECAT7:=
End Bin Normalisation From ECAT7:=
Bin Normalisation to apply second:= From ProjData
Bin Normalisation From ProjData :=
normalisation projdata filename:= AN_ATTENUATION_SINOGRAM_tr1_meas_zoom1_a_f1g1d0b0.hs
End Bin Normalisation From ProjData:=
END Chained Bin Normalisation Parameters:=
; if the next parameter is disabled,
; the sensitivity will be computed using the normalisation object
sensitivity filename:= sens.img
; if next is set to 1, sensitivity will be recomputed
; and written to file (if "sensitivity filename" is set)
recompute sensitivity := 1
; specify additive projection data to handle randoms or so
; see User's Guide for more info
additive sinograms := SCATTER_AND_RANDOMS_ADDED_ECAT7.S
; patlak related files
Kinetic Model type := Patlak Plot
Patlak Plot Parameters :=
time frame definition filename := TIMEFRAMES.fdef
starting frame := FRAME_NUM
calibration factor := 1
blood data filename := TEXTFILE-TIME-PLASMA-BLOOD.if
; Time Shift in seconds
Time Shift := 0
; boolean - direct reconstructions and STIR, 0: CTI scale
In total counts := 0
; 0 (default), if 1 it automatically scale with xy_voxel_size/bin_size
;In correct scale := 1
end Patlak Plot Parameters :=
end PoissonLogLikelihoodWithLinearKineticModelAndDynamicProjectionData:=
; Number of subsets should be a divisor of num_views/4
number of subsets:=NUM_SUBSETS
; Use for starting the numbering from something else than 1
start at subiteration number:=1
; Use if you want to start from another subset than 0 (but why?)
start at subset:= 0
number of subiterations:= NUM_SUBITERATIONS
save estimates at subiteration intervals:= NUM_SUBITERATIONS_INTERVALS
initial estimate := PARAMETRIC_IMAGE_OF_TWO_PARAMETERS.img
; enable this when you read an initial estimate with negative data
enforce initial positivity condition:=1
write update image := 0
inter-update filter subiteration interval:= 0
inter-update filter type := None
inter-iteration filter subiteration interval:= 0
post-filter type := None
output filename prefix:= OUTPUT_PARAMETRIC_IMAGE.img
END :=