# HG changeset patch # User Tuomo Valkonen # Date 1714072838 18000 # Node ID b180a4f3b9bd3529cf9c17cbf534198ca4e68492 # Parent bba159cf1636e0db11c9165cd0c2ae340c3bdefd# Parent 75116ad1d2e640452f1d88f1c2ba44194801522d merge diff -r bba159cf1636 -r b180a4f3b9bd src/ImGenerate.jl diff -r bba159cf1636 -r b180a4f3b9bd src/PET/PET.jl --- a/src/PET/PET.jl Thu Apr 25 13:25:48 2024 -0500 +++ b/src/PET/PET.jl Thu Apr 25 14:20:38 2024 -0500 @@ -79,10 +79,10 @@ const p_known₀_pet = default_params ⬿ ( noise_level = 0.5, - shake_noise_level = 0.1, + shake_noise_level = 0.25, shake = 1.0, - rotation_factor = 0.075, - rotation_noise_level = 0.0075, + rotation_factor = 0.15, + rotation_noise_level = 0.035, α = 0.15, ρ̃₀ = 1.0, σ̃₀ = 1.0, @@ -101,27 +101,27 @@ stable_interval = union(Set(1000:2000),Set(3500:4000)), ) -const p_known₀_petb = p_known₀_pet# ⬿ ( seed = 313159, ) +const shepplogan = imgen_shepplogan_radon(p_known₀_pet.sz) +const p_known₀_pets = p_known₀_pet ⬿ ( seed = 314159, ) -const shepplogan = imgen_shepplogan_radon(p_known₀_pet.sz) - +const p_known₀_petb = p_known₀_pet ⬿ ( seed = 9182737465, ) const brainphantom = imgen_brainphantom_radon(p_known₀_pet.sz) const shepplogan_experiments_pdps_known = ( Experiment(AlgorithmNew, DisplacementConstant, shepplogan, - p_known₀_pet ⬿ (predictor=DualScaling(),)), + p_known₀_pets ⬿ (predictor=DualScaling(),)), Experiment(AlgorithmNew, DisplacementConstant, shepplogan, - p_known₀_pet ⬿ (predictor=Greedy(),)), + p_known₀_pets ⬿ (predictor=Greedy(),)), Experiment(AlgorithmNew, DisplacementConstant, shepplogan, - p_known₀_pet ⬿ (predictor=nothing,),), + p_known₀_pets ⬿ (predictor=nothing,),), Experiment(AlgorithmNew, DisplacementConstant, shepplogan, - p_known₀_pet ⬿ (predictor=PrimalOnly(),)), + p_known₀_pets ⬿ (predictor=PrimalOnly(),)), Experiment(AlgorithmProximal, DisplacementConstant, shepplogan, - p_known₀_pet ⬿ (phantom_ρ = 100,)), + p_known₀_pets ⬿ (phantom_ρ = 100,)), Experiment(AlgorithmNew, DisplacementConstant, shepplogan, - p_known₀_pet ⬿ (predictor=Rotation(),)), + p_known₀_pets ⬿ (predictor=Rotation(),)), Experiment(AlgorithmNew, DisplacementConstant, shepplogan, - p_known₀_pet ⬿ (predictor=ZeroDual(),)), + p_known₀_pets ⬿ (predictor=ZeroDual(),)), ) const brainphantom_experiments_pdps_known = (