/
OSMAPOSL_MRP.par
67 lines (52 loc) · 2.03 KB
/
OSMAPOSL_MRP.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
OSMAPOSLParameters :=
; example file for running MRP (in its One Step Late flavour)
; with subsets
; see OSMAPOSL_osem_with_interfiltering.par for other parameters
objective function type:= PoissonLogLikelihoodWithLinearModelForMeanAndProjData
PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=
input file := test.hs
; if disabled, defaults to maximum segment number in the file
maximum absolute segment number to process := 4
; see User's Guide to see when you need this
zero end planes of segment 0:= 1
; if the next parameter is disabled,
; the sensitivity will be computed
sensitivity filename:= sens.hv
; specify additive projection data to handle randoms or so
; see User's Guide for more info
additive sinogram := 0
; here comes the MRP stuff
prior type := FilterRootPrior
FilterRootPrior Parameters :=
penalisation factor := 1
; you can use any image processor here
; the next parameters specify a 3x3x3 median
Filter type := Median
Median Filter Parameters :=
mask radius x := 1
mask radius y := 1
mask radius z := 1
End Median Filter Parameters:=
END FilterRootPrior Parameters :=
end PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=
initial estimate:= some_image
; enable this when you read an initial estimate with negative data
enforce initial positivity condition:=0
output filename prefix := test_MRP
; Number of subsets should usually be a divisor of num_views/4
number of subsets:= 12
; 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:= 24
save estimates at subiteration intervals:= 12
; following are defaults, so are not necessary
inter-update filter subiteration interval:= 0
inter-update filter type := None
inter-iteration filter subiteration interval:= 0
inter-iteration filter type := None
post-filter type := None
; enable this for multiplicative form of OSMAPOSL-MRP (see User's Guide)
;MAP model := multiplicative
END :=