OSMAPOSL_test_PMFromFile_QPweights.par

OSMAPOSLParameters :=
; test file for OSMAPOSL with a quadratic prior (and ray tracing projection matrix)
objective function type:= PoissonLogLikelihoodWithLinearModelForMeanAndProjData
PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

input file := Utahscat600k_ca_seg4.hs
zero end planes of segment 0:= 1
; if disabled, defaults to maximum segment number in the file
maximum absolute segment number to process := 3

; change to STIR 2.x default for compatibility 
use subset sensitivities:=0
sensitivity filename:= RPTsens_seg3_PM.hv

projector pair type := Matrix
  Projector Pair Using Matrix Parameters :=
  Matrix type := From File
Projection Matrix By Bin From File Parameters:=
Version := 1.0
symmetries type := PET_CartesianGrid
 PET_CartesianGrid symmetries parameters:=
  do_symmetry_90degrees_min_phi:= 1
  do_symmetry_180degrees_min_phi:= 1
  do_symmetry_swap_segment:= 1
  do_symmetry_swap_s:= 1
  do_symmetry_shift_z:= 1
 End PET_CartesianGrid symmetries parameters:=
template proj data filename:=my_PMRT_template_proj_data.hs
template density filename:=my_PMRT_template_density.hv
data_filename:=my_PMRT.pm
End Projection Matrix By Bin From File Parameters:=
  End Projector Pair Using Matrix Parameters :=

prior type := quadratic
  Quadratic Prior Parameters:=
  penalisation factor := .9
  weights:={{{0,1,0},{1,0,1},{0,1,0}}}
  END Quadratic Prior Parameters:=


end PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

output filename prefix := my_test_image_PMFromFile_QPweights

number of subsets:= 12
start at subset:= 0
number of subiterations:= 6
save estimates at subiteration intervals:= 6

initial estimate := my_uniform_image_circular.hv

map model := multiplicative

; for compatibility with STIR 1.x
inter-iteration filter subiteration interval:= 1
inter-iteration filter type := Truncate To Cylindrical FOV
Truncate To Cylindrical FOV Parameters:=
End Truncate To Cylindrical FOV Parameters:=

END :=