1 snark14.s170710 - 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 PERFECT
177 projection data are noiseless
183 background absorption 0.0000
187 0.438 seconds phantom creation
188 1.219 seconds projection data creation
189 1.656 seconds used for processing command crea
196 EXAMPLE 11 Brain Phantom
200 energy spectrum is monochromatic at energy level 511
204 description of objects
206 numb type x-coord y-coord x-length y-length angle av dens 511
208 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
210 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
212 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
214 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
216 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
218 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
220 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
222 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
224 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
226 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
228 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
230 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
232 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
234 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
236 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
238 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
240 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
242 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
244 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
246 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
248 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
250 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
252 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
254 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
256 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
258 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
260 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
262 scale factor multiplying object densities 0.5100
265 inhomogeneity set to 0.0500
269 density in each pixel is obtained as the average of 7 x 7 points
272 <#> pixe 381 size 0.4000
273 picture size 381 x 381, pixel size 0.4000
276 EXAMPLE 11 Brain Phantom
277 0.062 seconds used for processing command pict
284 EXAMPLE 11 Brain Phantom
288 energy spectrum is monochromatic at energy level 511
292 description of objects
294 numb type x-coord y-coord x-length y-length angle av dens 511
296 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
298 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
300 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
302 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
304 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
306 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
308 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
310 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
312 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
314 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
316 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
318 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
320 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
322 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
324 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
326 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
328 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
330 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
332 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
334 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
336 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
338 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
340 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
342 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
344 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
346 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
348 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
350 scale factor multiplying object densities 0.5100
353 inhomogeneity set to 0.0500
357 projection data are calculated by dividing each ray interval into 1 substrips
359 with aperture (substrip) weights 1
365 <#> dive arc source at 153.0000 det dist 306.0000
366 rays are divergent from point sources
367 source to origin distance 153.0000
368 the detectors lie on an arc with source to detector distance = 306.0000
371 <#> rays user 301 spacing 1.1000
372 number of rays per projection 301
373 snark computed number of rays 437
374 at detector spacing 1.1000
378 total number of projections 500
381 projection angles 0.0000 0.7190 1.4381 2.1571 2.8762 3.5952 4.3142 5.0333 5.7523 6.4713
382 7.1904 7.9094 8.6285 9.3475 10.0665 10.7856 11.5046 12.2236 12.9427 13.6617
383 14.3808 15.0998 15.8188 16.5379 17.2569 17.9760 18.6950 19.4140 20.1331 20.8521
384 21.5711 22.2902 23.0092 23.7283 24.4473 25.1663 25.8854 26.6044 27.3234 28.0425
385 28.7615 29.4806 30.1996 30.9186 31.6377 32.3567 33.0758 33.7948 34.5138 35.2329
386 35.9519 36.6709 37.3900 38.1090 38.8281 39.5471 40.2661 40.9852 41.7042 42.4232
387 43.1423 43.8613 44.5804 45.2994 46.0184 46.7375 47.4565 48.1756 48.8946 49.6136
388 50.3327 51.0517 51.7707 52.4898 53.2088 53.9279 54.6469 55.3659 56.0850 56.8040
389 57.5230 58.2421 58.9611 59.6802 60.3992 61.1182 61.8373 62.5563 63.2754 63.9944
390 64.7134 65.4325 66.1515 66.8705 67.5896 68.3086 69.0277 69.7467 70.4657 71.1848
391 71.9038 72.6228 73.3419 74.0609 74.7800 75.4990 76.2180 76.9371 77.6561 78.3752
392 79.0942 79.8132 80.5323 81.2513 81.9703 82.6894 83.4084 84.1275 84.8465 85.5655
393 86.2846 87.0036 87.7226 88.4417 89.1607 89.8798 90.5988 91.3178 92.0369 92.7559
394 93.4749 94.1940 94.9130 95.6321 96.3511 97.0701 97.7892 98.5082 99.2273 99.9463
395 100.6653 101.3844 102.1034 102.8224 103.5415 104.2605 104.9796 105.6986 106.4176 107.1367
396 107.8557 108.5747 109.2938 110.0128 110.7319 111.4509 112.1699 112.8890 113.6080 114.3271
397 115.0461 115.7651 116.4842 117.2032 117.9222 118.6413 119.3603 120.0794 120.7984 121.5174
398 122.2365 122.9555 123.6745 124.3936 125.1126 125.8317 126.5507 127.2697 127.9888 128.7078
399 129.4269 130.1459 130.8649 131.5840 132.3030 133.0220 133.7411 134.4601 135.1792 135.8982
400 136.6172 137.3363 138.0553 138.7743 139.4934 140.2124 140.9315 141.6505 142.3695 143.0886
401 143.8076 144.5267 145.2457 145.9647 146.6838 147.4028 148.1218 148.8409 149.5599 150.2790
402 150.9980 151.7170 152.4361 153.1551 153.8741 154.5932 155.3122 156.0313 156.7503 157.4693
403 158.1884 158.9074 159.6265 160.3455 161.0645 161.7836 162.5026 163.2216 163.9407 164.6597
404 165.3788 166.0978 166.8168 167.5359 168.2549 168.9739 169.6930 170.4120 171.1311 171.8501
405 172.5691 173.2882 174.0072 174.7263 175.4453 176.1643 176.8834 177.6024 178.3214 179.0405
406 179.7595 180.4786 181.1976 181.9166 182.6357 183.3547 184.0737 184.7928 185.5118 186.2309
407 186.9499 187.6689 188.3880 189.1070 189.8261 190.5451 191.2641 191.9832 192.7022 193.4212
408 194.1403 194.8593 195.5784 196.2974 197.0164 197.7355 198.4545 199.1735 199.8926 200.6116
409 201.3307 202.0497 202.7687 203.4878 204.2068 204.9259 205.6449 206.3639 207.0830 207.8020
410 208.5210 209.2401 209.9591 210.6782 211.3972 212.1162 212.8353 213.5543 214.2733 214.9924
411 215.7114 216.4305 217.1495 217.8685 218.5876 219.3066 220.0257 220.7447 221.4637 222.1828
412 222.9018 223.6208 224.3399 225.0589 225.7780 226.4970 227.2160 227.9351 228.6541 229.3731
413 230.0922 230.8112 231.5303 232.2493 232.9683 233.6874 234.4064 235.1255 235.8445 236.5635
414 237.2826 238.0016 238.7206 239.4397 240.1587 240.8778 241.5968 242.3158 243.0349 243.7539
415 244.4729 245.1920 245.9110 246.6301 247.3491 248.0681 248.7872 249.5062 250.2253 250.9443
416 251.6633 252.3824 253.1014 253.8204 254.5395 255.2585 255.9776 256.6966 257.4156 258.1347
417 258.8537 259.5727 260.2918 261.0108 261.7299 262.4489 263.1679 263.8870 264.6060 265.3251
418 266.0441 266.7631 267.4822 268.2012 268.9202 269.6393 270.3583 271.0774 271.7964 272.5154
419 273.2345 273.9535 274.6725 275.3916 276.1106 276.8297 277.5487 278.2677 278.9868 279.7058
420 280.4248 281.1439 281.8629 282.5820 283.3010 284.0200 284.7391 285.4581 286.1772 286.8962
421 287.6152 288.3343 289.0533 289.7723 290.4914 291.2104 291.9295 292.6485 293.3675 294.0866
422 294.8056 295.5246 296.2437 296.9627 297.6818 298.4008 299.1198 299.8389 300.5579 301.2770
423 301.9960 302.7150 303.4341 304.1531 304.8721 305.5912 306.3102 307.0293 307.7483 308.4673
424 309.1864 309.9054 310.6244 311.3435 312.0625 312.7816 313.5006 314.2196 314.9387 315.6577
425 316.3768 317.0958 317.8148 318.5339 319.2529 319.9719 320.6910 321.4100 322.1291 322.8481
426 323.5671 324.2862 325.0052 325.7242 326.4433 327.1623 327.8814 328.6004 329.3194 330.0385
427 330.7575 331.4766 332.1956 332.9146 333.6337 334.3527 335.0717 335.7908 336.5098 337.2289
428 337.9479 338.6669 339.3860 340.1050 340.8240 341.5431 342.2621 342.9812 343.7002 344.4192
429 345.1383 345.8573 346.5764 347.2954 348.0144 348.7335 349.4525 350.1715 350.8906 351.6096
430 352.3287 353.0477 353.7667 354.4858 355.2048 355.9238 356.6429 357.3619 358.0810 358.8000
434 projection data are noiseless
440 background absorption 0.0000
442 estimate of totlen = 20962469.574466
443 estimate of totden = 9803802.737080
444 estimate of average density = 0.4677
446 EXAMPLE 11 Brain Phantom
447 0.094 seconds used for processing command proj
452 <#> STOP TERMINATION KLDS 100000 RPRT
453 termination test klds
455 reporting file: RPRTklds
456 reporting on every iteration
458 0.000 seconds used for processing command stop
463 <#> SUPERIORIZE 16 0.999 1 SMOO RPRT
464 Superiorization is enabled
468 secondary criterion: smoo
470 reporting file: RPRTsuperiorization
471 reporting on every iteration
472 0.000 seconds used for processing command supe
475 <#> EXECUTE AVERAGE EMAP
477 Example 11 Illustrating the Superiorized MAP EM algorithm for PET
481 -----------------------------------------------------------
483 maximum a-posteriori probability expectation maximization
486 evaluation flag is set
488 -----------------------------------------------------------
490 *** WARNING - There exist negative values in the projection data ***
491 *** Negative values set to 0.0 ***
493 value of phi before algorithm operator: 4.42622e-28
494 value of phi after algorithm operator: 2.39806
495 algorithm executed in iteration 1
496 1.734 seconds for the execution of the algorithm
497 ****Error: Negative raysum in prjfil detected, KLDS is not applicable! Terminating algorithm execution.
498 reconstruction completed after iteration 1
499 1.734 seconds for this iteration
500 1.922 seconds used for processing command exec