1 snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM
6 <*> SUPERIORIZED MAP EM ALGORITHM FOR EMISSION TOMOGRAPHY. RECONSTRUCTION
8 <*> OF BRAIN PHANTOM. SIMULATING PET GEOMETRY WITH A RING OF 300 DETECTORS
10 <*> WITH EACH DETECTOR IN COINCIDENCE WITH 101 DETECTORS OPPOSITE IT.
12 <*> COEFFICIENT OF PENALTLY TERM IS SET TO 0
18 EXAMPLE 11 Brain Phantom
21 <#> SPECTRUM MONOCHROMATIC 511
22 energy spectrum is monochromatic at energy level 511
26 description of objects
28 numb type x-coord y-coord x-length y-length angle av dens 511
30 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.9500 0.9500
32 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 1.0000 1.0000
34 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 1.0000 1.0000
36 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.9500 0.9500
38 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.9500 0.9500
40 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 1.0000 1.0000
42 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.9500 0.9500
44 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 1.0000 1.0000
46 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 1.0000 1.0000
48 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.9500 0.9500
50 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.9500 0.9500
52 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 1.0000 1.0000
54 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.9500 0.9500
56 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 1.0000 1.0000
58 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 1.0000 1.0000
60 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.9500 0.9500
62 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 1.0000 1.0000
64 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.9500 0.9500
66 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.9500 0.9500
68 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 1.0000 1.0000
70 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 1.0000 1.0000
72 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.9500 0.9500
74 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 1.0000 1.0000
76 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.9500 0.9500
78 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.9500 0.9500
80 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 1.0000 1.0000
82 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 1.0000 1.0000
84 scale factor multiplying object densities 0.5100
87 inhomogeneity set to 0.0500
91 this run will generate a phantom
92 density in each pixel is obtained as the average of 7 x 7 points
95 <#> 381 PIXELS OF SIZE 0.4
96 picture size 381 x 381, pixel size 0.4000
101 this run will generate projection data
102 projection data are calculated by dividing each ray interval into 1 substrips
104 with aperture (substrip) weights 1
110 <#> divergent arc 153 306
111 rays are divergent from point sources
112 source to origin distance 153.0000
113 the detectors lie on an arc with source to detector distance = 306.0000
116 <#> RAYS USER 301 DETECTOR SPACING 1.1
117 number of rays per projection 301
118 at detector spacing 1.1000
121 <#> ANGLES 500 EQUAL SPACING
122 total number of projections 500
124 projection angles 0.0000 0.7190 1.4381 2.1571 2.8762 3.5952 4.3142 5.0333 5.7523 6.4713
125 7.1904 7.9094 8.6285 9.3475 10.0665 10.7856 11.5046 12.2236 12.9427 13.6617
126 14.3808 15.0998 15.8188 16.5379 17.2569 17.9760 18.6950 19.4140 20.1331 20.8521
127 21.5711 22.2902 23.0092 23.7283 24.4473 25.1663 25.8854 26.6044 27.3234 28.0425
128 28.7615 29.4806 30.1996 30.9186 31.6377 32.3567 33.0758 33.7948 34.5138 35.2329
129 35.9519 36.6709 37.3900 38.1090 38.8281 39.5471 40.2661 40.9852 41.7042 42.4232
130 43.1423 43.8613 44.5804 45.2994 46.0184 46.7375 47.4565 48.1756 48.8946 49.6136
131 50.3327 51.0517 51.7707 52.4898 53.2088 53.9279 54.6469 55.3659 56.0850 56.8040
132 57.5230 58.2421 58.9611 59.6802 60.3992 61.1182 61.8373 62.5563 63.2754 63.9944
133 64.7134 65.4325 66.1515 66.8705 67.5896 68.3086 69.0277 69.7467 70.4657 71.1848
134 71.9038 72.6228 73.3419 74.0609 74.7800 75.4990 76.2180 76.9371 77.6561 78.3752
135 79.0942 79.8132 80.5323 81.2513 81.9703 82.6894 83.4084 84.1275 84.8465 85.5655
136 86.2846 87.0036 87.7226 88.4417 89.1607 89.8798 90.5988 91.3178 92.0369 92.7559
137 93.4749 94.1940 94.9130 95.6321 96.3511 97.0701 97.7892 98.5082 99.2273 99.9463
138 100.6653 101.3844 102.1034 102.8224 103.5415 104.2605 104.9796 105.6986 106.4176 107.1367
139 107.8557 108.5747 109.2938 110.0128 110.7319 111.4509 112.1699 112.8890 113.6080 114.3271
140 115.0461 115.7651 116.4842 117.2032 117.9222 118.6413 119.3603 120.0794 120.7984 121.5174
141 122.2365 122.9555 123.6745 124.3936 125.1126 125.8317 126.5507 127.2697 127.9888 128.7078
142 129.4269 130.1459 130.8649 131.5840 132.3030 133.0220 133.7411 134.4601 135.1792 135.8982
143 136.6172 137.3363 138.0553 138.7743 139.4934 140.2124 140.9315 141.6505 142.3695 143.0886
144 143.8076 144.5267 145.2457 145.9647 146.6838 147.4028 148.1218 148.8409 149.5599 150.2790
145 150.9980 151.7170 152.4361 153.1551 153.8741 154.5932 155.3122 156.0313 156.7503 157.4693
146 158.1884 158.9074 159.6265 160.3455 161.0645 161.7836 162.5026 163.2216 163.9407 164.6597
147 165.3788 166.0978 166.8168 167.5359 168.2549 168.9739 169.6930 170.4120 171.1311 171.8501
148 172.5691 173.2882 174.0072 174.7263 175.4453 176.1643 176.8834 177.6024 178.3214 179.0405
149 179.7595 180.4786 181.1976 181.9166 182.6357 183.3547 184.0737 184.7928 185.5118 186.2309
150 186.9499 187.6689 188.3880 189.1070 189.8261 190.5451 191.2641 191.9832 192.7022 193.4212
151 194.1403 194.8593 195.5784 196.2974 197.0164 197.7355 198.4545 199.1735 199.8926 200.6116
152 201.3307 202.0497 202.7687 203.4878 204.2068 204.9259 205.6449 206.3639 207.0830 207.8020
153 208.5210 209.2401 209.9591 210.6782 211.3972 212.1162 212.8353 213.5543 214.2733 214.9924
154 215.7114 216.4305 217.1495 217.8685 218.5876 219.3066 220.0257 220.7447 221.4637 222.1828
155 222.9018 223.6208 224.3399 225.0589 225.7780 226.4970 227.2160 227.9351 228.6541 229.3731
156 230.0922 230.8112 231.5303 232.2493 232.9683 233.6874 234.4064 235.1255 235.8445 236.5635
157 237.2826 238.0016 238.7206 239.4397 240.1587 240.8778 241.5968 242.3158 243.0349 243.7539
158 244.4729 245.1920 245.9110 246.6301 247.3491 248.0681 248.7872 249.5062 250.2253 250.9443
159 251.6633 252.3824 253.1014 253.8204 254.5395 255.2585 255.9776 256.6966 257.4156 258.1347
160 258.8537 259.5727 260.2918 261.0108 261.7299 262.4489 263.1679 263.8870 264.6060 265.3251
161 266.0441 266.7631 267.4822 268.2012 268.9202 269.6393 270.3583 271.0774 271.7964 272.5154
162 273.2345 273.9535 274.6725 275.3916 276.1106 276.8297 277.5487 278.2677 278.9868 279.7058
163 280.4248 281.1439 281.8629 282.5820 283.3010 284.0200 284.7391 285.4581 286.1772 286.8962
164 287.6152 288.3343 289.0533 289.7723 290.4914 291.2104 291.9295 292.6485 293.3675 294.0866
165 294.8056 295.5246 296.2437 296.9627 297.6818 298.4008 299.1198 299.8389 300.5579 301.2770
166 301.9960 302.7150 303.4341 304.1531 304.8721 305.5912 306.3102 307.0293 307.7483 308.4673
167 309.1864 309.9054 310.6244 311.3435 312.0625 312.7816 313.5006 314.2196 314.9387 315.6577
168 316.3768 317.0958 317.8148 318.5339 319.2529 319.9719 320.6910 321.4100 322.1291 322.8481
169 323.5671 324.2862 325.0052 325.7242 326.4433 327.1623 327.8814 328.6004 329.3194 330.0385
170 330.7575 331.4766 332.1956 332.9146 333.6337 334.3527 335.0717 335.7908 336.5098 337.2289
171 337.9479 338.6669 339.3860 340.1050 340.8240 341.5431 342.2621 342.9812 343.7002 344.4192
172 345.1383 345.8573 346.5764 347.2954 348.0144 348.7335 349.4525 350.1715 350.8906 351.6096
173 352.3287 353.0477 353.7667 354.4858 355.2048 355.9238 356.6429 357.3619 358.0810 358.8000
176 <#> MEASUREMENT NOISY
177 noise characteristics of projection data follow
178 nature characteristics
180 <#> QUANTUM 1.0 1.0 CALIBRATION 4
184 seed for random number generator is 0
190 background absorption 0.0000
194 0.358 seconds phantom creation
195 1.289 seconds projection data creation
196 1.647 seconds used for processing command crea
203 EXAMPLE 11 Brain Phantom
207 energy spectrum is monochromatic at energy level 511
211 description of objects
213 numb type x-coord y-coord x-length y-length angle av dens 511
215 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
217 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
219 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
221 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
223 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
225 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
227 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
229 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
231 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
233 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
235 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
237 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
239 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
241 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
243 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
245 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
247 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
249 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
251 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
253 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
255 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
257 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
259 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
261 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
263 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
265 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
267 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
269 scale factor multiplying object densities 0.5100
272 inhomogeneity set to 0.0500
276 density in each pixel is obtained as the average of 7 x 7 points
279 <#> pixe 381 size 0.4000
280 picture size 381 x 381, pixel size 0.4000
283 EXAMPLE 11 Brain Phantom
284 0.042 seconds used for processing command pict
291 EXAMPLE 11 Brain Phantom
295 energy spectrum is monochromatic at energy level 511
299 description of objects
301 numb type x-coord y-coord x-length y-length angle av dens 511
303 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
305 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
307 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
309 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
311 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
313 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
315 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
317 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
319 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
321 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
323 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
325 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
327 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
329 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
331 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
333 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
335 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
337 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
339 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
341 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
343 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
345 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
347 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
349 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
351 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
353 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
355 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
357 scale factor multiplying object densities 0.5100
360 inhomogeneity set to 0.0500
364 projection data are calculated by dividing each ray interval into 1 substrips
366 with aperture (substrip) weights 1
372 <#> dive arc source at 153.0000 det dist 306.0000
373 rays are divergent from point sources
374 source to origin distance 153.0000
375 the detectors lie on an arc with source to detector distance = 306.0000
378 <#> rays user 301 spacing 1.1000
379 number of rays per projection 301
380 snark computed number of rays 437
381 at detector spacing 1.1000
385 total number of projections 500
388 projection angles 0.0000 0.7190 1.4381 2.1571 2.8762 3.5952 4.3142 5.0333 5.7523 6.4713
389 7.1904 7.9094 8.6285 9.3475 10.0665 10.7856 11.5046 12.2236 12.9427 13.6617
390 14.3808 15.0998 15.8188 16.5379 17.2569 17.9760 18.6950 19.4140 20.1331 20.8521
391 21.5711 22.2902 23.0092 23.7283 24.4473 25.1663 25.8854 26.6044 27.3234 28.0425
392 28.7615 29.4806 30.1996 30.9186 31.6377 32.3567 33.0758 33.7948 34.5138 35.2329
393 35.9519 36.6709 37.3900 38.1090 38.8281 39.5471 40.2661 40.9852 41.7042 42.4232
394 43.1423 43.8613 44.5804 45.2994 46.0184 46.7375 47.4565 48.1756 48.8946 49.6136
395 50.3327 51.0517 51.7707 52.4898 53.2088 53.9279 54.6469 55.3659 56.0850 56.8040
396 57.5230 58.2421 58.9611 59.6802 60.3992 61.1182 61.8373 62.5563 63.2754 63.9944
397 64.7134 65.4325 66.1515 66.8705 67.5896 68.3086 69.0277 69.7467 70.4657 71.1848
398 71.9038 72.6228 73.3419 74.0609 74.7800 75.4990 76.2180 76.9371 77.6561 78.3752
399 79.0942 79.8132 80.5323 81.2513 81.9703 82.6894 83.4084 84.1275 84.8465 85.5655
400 86.2846 87.0036 87.7226 88.4417 89.1607 89.8798 90.5988 91.3178 92.0369 92.7559
401 93.4749 94.1940 94.9130 95.6321 96.3511 97.0701 97.7892 98.5082 99.2273 99.9463
402 100.6653 101.3844 102.1034 102.8224 103.5415 104.2605 104.9796 105.6986 106.4176 107.1367
403 107.8557 108.5747 109.2938 110.0128 110.7319 111.4509 112.1699 112.8890 113.6080 114.3271
404 115.0461 115.7651 116.4842 117.2032 117.9222 118.6413 119.3603 120.0794 120.7984 121.5174
405 122.2365 122.9555 123.6745 124.3936 125.1126 125.8317 126.5507 127.2697 127.9888 128.7078
406 129.4269 130.1459 130.8649 131.5840 132.3030 133.0220 133.7411 134.4601 135.1792 135.8982
407 136.6172 137.3363 138.0553 138.7743 139.4934 140.2124 140.9315 141.6505 142.3695 143.0886
408 143.8076 144.5267 145.2457 145.9647 146.6838 147.4028 148.1218 148.8409 149.5599 150.2790
409 150.9980 151.7170 152.4361 153.1551 153.8741 154.5932 155.3122 156.0313 156.7503 157.4693
410 158.1884 158.9074 159.6265 160.3455 161.0645 161.7836 162.5026 163.2216 163.9407 164.6597
411 165.3788 166.0978 166.8168 167.5359 168.2549 168.9739 169.6930 170.4120 171.1311 171.8501
412 172.5691 173.2882 174.0072 174.7263 175.4453 176.1643 176.8834 177.6024 178.3214 179.0405
413 179.7595 180.4786 181.1976 181.9166 182.6357 183.3547 184.0737 184.7928 185.5118 186.2309
414 186.9499 187.6689 188.3880 189.1070 189.8261 190.5451 191.2641 191.9832 192.7022 193.4212
415 194.1403 194.8593 195.5784 196.2974 197.0164 197.7355 198.4545 199.1735 199.8926 200.6116
416 201.3307 202.0497 202.7687 203.4878 204.2068 204.9259 205.6449 206.3639 207.0830 207.8020
417 208.5210 209.2401 209.9591 210.6782 211.3972 212.1162 212.8353 213.5543 214.2733 214.9924
418 215.7114 216.4305 217.1495 217.8685 218.5876 219.3066 220.0257 220.7447 221.4637 222.1828
419 222.9018 223.6208 224.3399 225.0589 225.7780 226.4970 227.2160 227.9351 228.6541 229.3731
420 230.0922 230.8112 231.5303 232.2493 232.9683 233.6874 234.4064 235.1255 235.8445 236.5635
421 237.2826 238.0016 238.7206 239.4397 240.1587 240.8778 241.5968 242.3158 243.0349 243.7539
422 244.4729 245.1920 245.9110 246.6301 247.3491 248.0681 248.7872 249.5062 250.2253 250.9443
423 251.6633 252.3824 253.1014 253.8204 254.5395 255.2585 255.9776 256.6966 257.4156 258.1347
424 258.8537 259.5727 260.2918 261.0108 261.7299 262.4489 263.1679 263.8870 264.6060 265.3251
425 266.0441 266.7631 267.4822 268.2012 268.9202 269.6393 270.3583 271.0774 271.7964 272.5154
426 273.2345 273.9535 274.6725 275.3916 276.1106 276.8297 277.5487 278.2677 278.9868 279.7058
427 280.4248 281.1439 281.8629 282.5820 283.3010 284.0200 284.7391 285.4581 286.1772 286.8962
428 287.6152 288.3343 289.0533 289.7723 290.4914 291.2104 291.9295 292.6485 293.3675 294.0866
429 294.8056 295.5246 296.2437 296.9627 297.6818 298.4008 299.1198 299.8389 300.5579 301.2770
430 301.9960 302.7150 303.4341 304.1531 304.8721 305.5912 306.3102 307.0293 307.7483 308.4673
431 309.1864 309.9054 310.6244 311.3435 312.0625 312.7816 313.5006 314.2196 314.9387 315.6577
432 316.3768 317.0958 317.8148 318.5339 319.2529 319.9719 320.6910 321.4100 322.1291 322.8481
433 323.5671 324.2862 325.0052 325.7242 326.4433 327.1623 327.8814 328.6004 329.3194 330.0385
434 330.7575 331.4766 332.1956 332.9146 333.6337 334.3527 335.0717 335.7908 336.5098 337.2289
435 337.9479 338.6669 339.3860 340.1050 340.8240 341.5431 342.2621 342.9812 343.7002 344.4192
436 345.1383 345.8573 346.5764 347.2954 348.0144 348.7335 349.4525 350.1715 350.8906 351.6096
437 352.3287 353.0477 353.7667 354.4858 355.2048 355.9238 356.6429 357.3619 358.0810 358.8000
441 noise characteristics of projection data follow
442 nature characteristics
444 <#> quan 1.0000 1.0000 cali 4
448 seed for random number generator is 0
454 background absorption 0.0000
456 estimate of totlen = 20962469.574466
457 estimate of totden = 9802000.000000
458 estimate of average density = 0.4676
460 EXAMPLE 11 Brain Phantom
461 0.054 seconds used for processing command proj
466 <#> STOP TERMINATION KLDS 100000 RPRT
467 termination test klds
469 reporting file: RPRTklds
470 reporting on every iteration
472 0.000 seconds used for processing command stop
477 <#> SUPERIORIZE 16 0.999 1 SMOO RPRT
478 Superiorization is enabled
482 secondary criterion: smoo
484 reporting file: RPRTsuperiorization
485 reporting on every iteration
486 0.000 seconds used for processing command supe
489 <#> EXECUTE AVERAGE EMAP
491 Example 11 Illustrating the Superiorized MAP EM algorithm for PET
495 -----------------------------------------------------------
497 maximum a-posteriori probability expectation maximization
500 evaluation flag is set
502 -----------------------------------------------------------
505 value of phi before algorithm operator: 4.42622e-28
506 value of phi after algorithm operator: 3.14467
507 algorithm executed in iteration 1
508 3.847 seconds for the execution of the algorithm
509 current epsilon (KL distance) = 638282
510 iteration 1 completed
511 4.805 seconds for this iteration
513 value of phi before algorithm operator: 0.429607
514 value of phi after algorithm operator: 1.65065
515 algorithm executed in iteration 2
516 2.144 seconds for the execution of the algorithm
517 current epsilon (KL distance) = 446785
518 iteration 2 completed
519 3.149 seconds for this iteration
521 value of phi before algorithm operator: 0.706119
522 value of phi after algorithm operator: 2.49879
523 algorithm executed in iteration 3
524 2.148 seconds for the execution of the algorithm
525 current epsilon (KL distance) = 331430
526 iteration 3 completed
527 3.106 seconds for this iteration
529 value of phi before algorithm operator: 0.564997
530 value of phi after algorithm operator: 2.14542
531 algorithm executed in iteration 4
532 2.143 seconds for the execution of the algorithm
533 current epsilon (KL distance) = 257464
534 iteration 4 completed
535 3.090 seconds for this iteration
537 value of phi before algorithm operator: 0.577605
538 value of phi after algorithm operator: 2.1505
539 algorithm executed in iteration 5
540 2.241 seconds for the execution of the algorithm
541 current epsilon (KL distance) = 207738
542 iteration 5 completed
543 3.451 seconds for this iteration
545 value of phi before algorithm operator: 0.579469
546 value of phi after algorithm operator: 2.11049
547 algorithm executed in iteration 6
548 2.264 seconds for the execution of the algorithm
549 current epsilon (KL distance) = 173224
550 iteration 6 completed
551 3.188 seconds for this iteration
553 value of phi before algorithm operator: 0.581496
554 value of phi after algorithm operator: 2.08666
555 algorithm executed in iteration 7
556 2.374 seconds for the execution of the algorithm
557 current epsilon (KL distance) = 148708
558 iteration 7 completed
559 3.408 seconds for this iteration
561 value of phi before algorithm operator: 0.587932
562 value of phi after algorithm operator: 2.07545
563 algorithm executed in iteration 8
564 2.618 seconds for the execution of the algorithm
565 current epsilon (KL distance) = 130971
566 iteration 8 completed
567 3.687 seconds for this iteration
569 value of phi before algorithm operator: 0.598034
570 value of phi after algorithm operator: 2.07283
571 algorithm executed in iteration 9
572 2.448 seconds for the execution of the algorithm
573 current epsilon (KL distance) = 117930
574 iteration 9 completed
575 3.523 seconds for this iteration
577 value of phi before algorithm operator: 0.610951
578 value of phi after algorithm operator: 2.07635
579 algorithm executed in iteration 10
580 2.199 seconds for the execution of the algorithm
581 current epsilon (KL distance) = 108196
582 iteration 10 completed
583 3.139 seconds for this iteration
585 value of phi before algorithm operator: 0.625866
586 value of phi after algorithm operator: 2.08405
587 algorithm executed in iteration 11
588 2.215 seconds for the execution of the algorithm
589 current epsilon (KL distance) = 100827
590 iteration 11 completed
591 3.599 seconds for this iteration
593 value of phi before algorithm operator: 0.64192
594 value of phi after algorithm operator: 2.09414
595 algorithm executed in iteration 12
596 2.428 seconds for the execution of the algorithm
597 current epsilon (KL distance) = 95173.6
598 reconstruction completed after iteration 12
599 3.422 seconds for this iteration
600 41.565 seconds for all iterations
601 41.706 seconds used for processing command exec