This version uses PERFECT data so it can be used to test regression
across different platforms.
+WARNING
+
+This version does not pass regression testing. There is a very small difference in value of phi between platforms:
+
+gcc 5.4 on Windows 10 WSL and Centos original version match
+gcc 7.2 on Ubuntu 2017.10 does not match
--- /dev/null
+ snark14.s170710 - A PICTURE RECONSTRUCTION PROGRAM
+
+
+ <*> *******EXAMPLE 11
+
+ <*> SUPERIORIZED MAP EM ALGORITHM FOR EMISSION TOMOGRAPHY. RECONSTRUCTION
+
+ <*> OF BRAIN PHANTOM. SIMULATING PET GEOMETRY WITH A RING OF 300 DETECTORS
+
+ <*> WITH EACH DETECTOR IN COINCIDENCE WITH 101 DETECTORS OPPOSITE IT.
+
+ <*> COEFFICIENT OF PENALTLY TERM IS SET TO 0
+
+ <*>
+
+ <#> CREATE
+
+ EXAMPLE 11 Brain Phantom
+
+
+ <#> SPECTRUM MONOCHROMATIC 511
+ energy spectrum is monochromatic at energy level 511
+
+
+ <#> OBJECTS
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 511
+
+ 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.9500 0.9500
+
+ 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 1.0000 1.0000
+
+ 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 1.0000 1.0000
+
+ 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.9500 0.9500
+
+ 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.9500 0.9500
+
+ 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 1.0000 1.0000
+
+ 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.9500 0.9500
+
+ 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 1.0000 1.0000
+
+ 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 1.0000 1.0000
+
+ 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.9500 0.9500
+
+ 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.9500 0.9500
+
+ 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 1.0000 1.0000
+
+ 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.9500 0.9500
+
+ 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 1.0000 1.0000
+
+ 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 1.0000 1.0000
+
+ 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.9500 0.9500
+
+ 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 1.0000 1.0000
+
+ 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.9500 0.9500
+
+ 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.9500 0.9500
+
+ 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 1.0000 1.0000
+
+ 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 1.0000 1.0000
+
+ 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.9500 0.9500
+
+ 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 1.0000 1.0000
+
+ 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.9500 0.9500
+
+ 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.9500 0.9500
+
+ 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 1.0000 1.0000
+
+ 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 1.0000 1.0000
+
+ scale factor multiplying object densities 0.5100
+
+ seed set to 1
+ inhomogeneity set to 0.0500
+
+ <#> PHANTOM AVERAGE 7
+
+ this run will generate a phantom
+ density in each pixel is obtained as the average of 7 x 7 points
+
+
+ <#> 381 PIXELS OF SIZE 0.4
+ picture size 381 x 381, pixel size 0.4000
+
+
+ <#> RAYSUM AVERAGE 1
+
+ this run will generate projection data
+ projection data are calculated by dividing each ray interval into 1 substrips
+
+ with aperture (substrip) weights 1
+
+
+ <#> GEOMETRY
+
+
+ <#> divergent arc 153 306
+ rays are divergent from point sources
+ source to origin distance 153.0000
+ the detectors lie on an arc with source to detector distance = 306.0000
+
+
+ <#> RAYS USER 301 DETECTOR SPACING 1.1
+ number of rays per projection 301
+ at detector spacing 1.1000
+
+
+ <#> ANGLES 500 EQUAL SPACING
+ total number of projections 500
+
+ projection angles 0.0000 0.7190 1.4381 2.1571 2.8762 3.5952 4.3142 5.0333 5.7523 6.4713
+ 7.1904 7.9094 8.6285 9.3475 10.0665 10.7856 11.5046 12.2236 12.9427 13.6617
+ 14.3808 15.0998 15.8188 16.5379 17.2569 17.9760 18.6950 19.4140 20.1331 20.8521
+ 21.5711 22.2902 23.0092 23.7283 24.4473 25.1663 25.8854 26.6044 27.3234 28.0425
+ 28.7615 29.4806 30.1996 30.9186 31.6377 32.3567 33.0758 33.7948 34.5138 35.2329
+ 35.9519 36.6709 37.3900 38.1090 38.8281 39.5471 40.2661 40.9852 41.7042 42.4232
+ 43.1423 43.8613 44.5804 45.2994 46.0184 46.7375 47.4565 48.1756 48.8946 49.6136
+ 50.3327 51.0517 51.7707 52.4898 53.2088 53.9279 54.6469 55.3659 56.0850 56.8040
+ 57.5230 58.2421 58.9611 59.6802 60.3992 61.1182 61.8373 62.5563 63.2754 63.9944
+ 64.7134 65.4325 66.1515 66.8705 67.5896 68.3086 69.0277 69.7467 70.4657 71.1848
+ 71.9038 72.6228 73.3419 74.0609 74.7800 75.4990 76.2180 76.9371 77.6561 78.3752
+ 79.0942 79.8132 80.5323 81.2513 81.9703 82.6894 83.4084 84.1275 84.8465 85.5655
+ 86.2846 87.0036 87.7226 88.4417 89.1607 89.8798 90.5988 91.3178 92.0369 92.7559
+ 93.4749 94.1940 94.9130 95.6321 96.3511 97.0701 97.7892 98.5082 99.2273 99.9463
+ 100.6653 101.3844 102.1034 102.8224 103.5415 104.2605 104.9796 105.6986 106.4176 107.1367
+ 107.8557 108.5747 109.2938 110.0128 110.7319 111.4509 112.1699 112.8890 113.6080 114.3271
+ 115.0461 115.7651 116.4842 117.2032 117.9222 118.6413 119.3603 120.0794 120.7984 121.5174
+ 122.2365 122.9555 123.6745 124.3936 125.1126 125.8317 126.5507 127.2697 127.9888 128.7078
+ 129.4269 130.1459 130.8649 131.5840 132.3030 133.0220 133.7411 134.4601 135.1792 135.8982
+ 136.6172 137.3363 138.0553 138.7743 139.4934 140.2124 140.9315 141.6505 142.3695 143.0886
+ 143.8076 144.5267 145.2457 145.9647 146.6838 147.4028 148.1218 148.8409 149.5599 150.2790
+ 150.9980 151.7170 152.4361 153.1551 153.8741 154.5932 155.3122 156.0313 156.7503 157.4693
+ 158.1884 158.9074 159.6265 160.3455 161.0645 161.7836 162.5026 163.2216 163.9407 164.6597
+ 165.3788 166.0978 166.8168 167.5359 168.2549 168.9739 169.6930 170.4120 171.1311 171.8501
+ 172.5691 173.2882 174.0072 174.7263 175.4453 176.1643 176.8834 177.6024 178.3214 179.0405
+ 179.7595 180.4786 181.1976 181.9166 182.6357 183.3547 184.0737 184.7928 185.5118 186.2309
+ 186.9499 187.6689 188.3880 189.1070 189.8261 190.5451 191.2641 191.9832 192.7022 193.4212
+ 194.1403 194.8593 195.5784 196.2974 197.0164 197.7355 198.4545 199.1735 199.8926 200.6116
+ 201.3307 202.0497 202.7687 203.4878 204.2068 204.9259 205.6449 206.3639 207.0830 207.8020
+ 208.5210 209.2401 209.9591 210.6782 211.3972 212.1162 212.8353 213.5543 214.2733 214.9924
+ 215.7114 216.4305 217.1495 217.8685 218.5876 219.3066 220.0257 220.7447 221.4637 222.1828
+ 222.9018 223.6208 224.3399 225.0589 225.7780 226.4970 227.2160 227.9351 228.6541 229.3731
+ 230.0922 230.8112 231.5303 232.2493 232.9683 233.6874 234.4064 235.1255 235.8445 236.5635
+ 237.2826 238.0016 238.7206 239.4397 240.1587 240.8778 241.5968 242.3158 243.0349 243.7539
+ 244.4729 245.1920 245.9110 246.6301 247.3491 248.0681 248.7872 249.5062 250.2253 250.9443
+ 251.6633 252.3824 253.1014 253.8204 254.5395 255.2585 255.9776 256.6966 257.4156 258.1347
+ 258.8537 259.5727 260.2918 261.0108 261.7299 262.4489 263.1679 263.8870 264.6060 265.3251
+ 266.0441 266.7631 267.4822 268.2012 268.9202 269.6393 270.3583 271.0774 271.7964 272.5154
+ 273.2345 273.9535 274.6725 275.3916 276.1106 276.8297 277.5487 278.2677 278.9868 279.7058
+ 280.4248 281.1439 281.8629 282.5820 283.3010 284.0200 284.7391 285.4581 286.1772 286.8962
+ 287.6152 288.3343 289.0533 289.7723 290.4914 291.2104 291.9295 292.6485 293.3675 294.0866
+ 294.8056 295.5246 296.2437 296.9627 297.6818 298.4008 299.1198 299.8389 300.5579 301.2770
+ 301.9960 302.7150 303.4341 304.1531 304.8721 305.5912 306.3102 307.0293 307.7483 308.4673
+ 309.1864 309.9054 310.6244 311.3435 312.0625 312.7816 313.5006 314.2196 314.9387 315.6577
+ 316.3768 317.0958 317.8148 318.5339 319.2529 319.9719 320.6910 321.4100 322.1291 322.8481
+ 323.5671 324.2862 325.0052 325.7242 326.4433 327.1623 327.8814 328.6004 329.3194 330.0385
+ 330.7575 331.4766 332.1956 332.9146 333.6337 334.3527 335.0717 335.7908 336.5098 337.2289
+ 337.9479 338.6669 339.3860 340.1050 340.8240 341.5431 342.2621 342.9812 343.7002 344.4192
+ 345.1383 345.8573 346.5764 347.2954 348.0144 348.7335 349.4525 350.1715 350.8906 351.6096
+ 352.3287 353.0477 353.7667 354.4858 355.2048 355.9238 356.6429 357.3619 358.0810 358.8000
+
+
+ <#> MEASUREMENT PERFECT
+ projection data are noiseless
+
+
+ <#> BACKGROUND 0.0
+ at levels
+ 511
+ background absorption 0.0000
+
+
+ <#> RUN
+ 0.438 seconds phantom creation
+ 1.219 seconds projection data creation
+ 1.656 seconds used for processing command crea
+
+
+ <*>
+
+ <#> PICTURE TEST
+
+ EXAMPLE 11 Brain Phantom
+
+
+ <#> spec mono 511
+ energy spectrum is monochromatic at energy level 511
+
+
+ <#> obje
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 511
+
+ 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
+
+ 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
+
+ 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
+
+ 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
+
+ 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
+
+ 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
+
+ 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
+
+ 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
+
+ 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
+
+ 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
+
+ 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
+
+ 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
+
+ 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
+
+ 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
+
+ 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
+
+ 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
+
+ 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
+
+ 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
+
+ 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
+
+ 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
+
+ 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
+
+ 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
+
+ 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
+
+ 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
+
+ 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
+
+ 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
+
+ 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
+
+ scale factor multiplying object densities 0.5100
+
+ seed set to 1
+ inhomogeneity set to 0.0500
+
+ <#> phan aver 7
+
+ density in each pixel is obtained as the average of 7 x 7 points
+
+
+ <#> pixe 381 size 0.4000
+ picture size 381 x 381, pixel size 0.4000
+
+ test picture read
+ EXAMPLE 11 Brain Phantom
+ 0.062 seconds used for processing command pict
+
+
+ <*>
+
+ <#> PROJECTION REAL
+
+ EXAMPLE 11 Brain Phantom
+
+
+ <#> spec mono 511
+ energy spectrum is monochromatic at energy level 511
+
+
+ <#> obje
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 511
+
+ 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
+
+ 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
+
+ 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
+
+ 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
+
+ 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
+
+ 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
+
+ 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
+
+ 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
+
+ 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
+
+ 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
+
+ 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
+
+ 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
+
+ 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
+
+ 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
+
+ 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
+
+ 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
+
+ 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
+
+ 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
+
+ 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
+
+ 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
+
+ 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
+
+ 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
+
+ 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
+
+ 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
+
+ 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
+
+ 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
+
+ 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
+
+ scale factor multiplying object densities 0.5100
+
+ seed set to 1
+ inhomogeneity set to 0.0500
+
+ <#> rays aver 1
+
+ projection data are calculated by dividing each ray interval into 1 substrips
+
+ with aperture (substrip) weights 1
+
+
+ <#> geom
+
+
+ <#> dive arc source at 153.0000 det dist 306.0000
+ rays are divergent from point sources
+ source to origin distance 153.0000
+ the detectors lie on an arc with source to detector distance = 306.0000
+
+
+ <#> rays user 301 spacing 1.1000
+ number of rays per projection 301
+ snark computed number of rays 437
+ at detector spacing 1.1000
+
+
+ <#> angl 500
+ total number of projections 500
+
+
+ projection angles 0.0000 0.7190 1.4381 2.1571 2.8762 3.5952 4.3142 5.0333 5.7523 6.4713
+ 7.1904 7.9094 8.6285 9.3475 10.0665 10.7856 11.5046 12.2236 12.9427 13.6617
+ 14.3808 15.0998 15.8188 16.5379 17.2569 17.9760 18.6950 19.4140 20.1331 20.8521
+ 21.5711 22.2902 23.0092 23.7283 24.4473 25.1663 25.8854 26.6044 27.3234 28.0425
+ 28.7615 29.4806 30.1996 30.9186 31.6377 32.3567 33.0758 33.7948 34.5138 35.2329
+ 35.9519 36.6709 37.3900 38.1090 38.8281 39.5471 40.2661 40.9852 41.7042 42.4232
+ 43.1423 43.8613 44.5804 45.2994 46.0184 46.7375 47.4565 48.1756 48.8946 49.6136
+ 50.3327 51.0517 51.7707 52.4898 53.2088 53.9279 54.6469 55.3659 56.0850 56.8040
+ 57.5230 58.2421 58.9611 59.6802 60.3992 61.1182 61.8373 62.5563 63.2754 63.9944
+ 64.7134 65.4325 66.1515 66.8705 67.5896 68.3086 69.0277 69.7467 70.4657 71.1848
+ 71.9038 72.6228 73.3419 74.0609 74.7800 75.4990 76.2180 76.9371 77.6561 78.3752
+ 79.0942 79.8132 80.5323 81.2513 81.9703 82.6894 83.4084 84.1275 84.8465 85.5655
+ 86.2846 87.0036 87.7226 88.4417 89.1607 89.8798 90.5988 91.3178 92.0369 92.7559
+ 93.4749 94.1940 94.9130 95.6321 96.3511 97.0701 97.7892 98.5082 99.2273 99.9463
+ 100.6653 101.3844 102.1034 102.8224 103.5415 104.2605 104.9796 105.6986 106.4176 107.1367
+ 107.8557 108.5747 109.2938 110.0128 110.7319 111.4509 112.1699 112.8890 113.6080 114.3271
+ 115.0461 115.7651 116.4842 117.2032 117.9222 118.6413 119.3603 120.0794 120.7984 121.5174
+ 122.2365 122.9555 123.6745 124.3936 125.1126 125.8317 126.5507 127.2697 127.9888 128.7078
+ 129.4269 130.1459 130.8649 131.5840 132.3030 133.0220 133.7411 134.4601 135.1792 135.8982
+ 136.6172 137.3363 138.0553 138.7743 139.4934 140.2124 140.9315 141.6505 142.3695 143.0886
+ 143.8076 144.5267 145.2457 145.9647 146.6838 147.4028 148.1218 148.8409 149.5599 150.2790
+ 150.9980 151.7170 152.4361 153.1551 153.8741 154.5932 155.3122 156.0313 156.7503 157.4693
+ 158.1884 158.9074 159.6265 160.3455 161.0645 161.7836 162.5026 163.2216 163.9407 164.6597
+ 165.3788 166.0978 166.8168 167.5359 168.2549 168.9739 169.6930 170.4120 171.1311 171.8501
+ 172.5691 173.2882 174.0072 174.7263 175.4453 176.1643 176.8834 177.6024 178.3214 179.0405
+ 179.7595 180.4786 181.1976 181.9166 182.6357 183.3547 184.0737 184.7928 185.5118 186.2309
+ 186.9499 187.6689 188.3880 189.1070 189.8261 190.5451 191.2641 191.9832 192.7022 193.4212
+ 194.1403 194.8593 195.5784 196.2974 197.0164 197.7355 198.4545 199.1735 199.8926 200.6116
+ 201.3307 202.0497 202.7687 203.4878 204.2068 204.9259 205.6449 206.3639 207.0830 207.8020
+ 208.5210 209.2401 209.9591 210.6782 211.3972 212.1162 212.8353 213.5543 214.2733 214.9924
+ 215.7114 216.4305 217.1495 217.8685 218.5876 219.3066 220.0257 220.7447 221.4637 222.1828
+ 222.9018 223.6208 224.3399 225.0589 225.7780 226.4970 227.2160 227.9351 228.6541 229.3731
+ 230.0922 230.8112 231.5303 232.2493 232.9683 233.6874 234.4064 235.1255 235.8445 236.5635
+ 237.2826 238.0016 238.7206 239.4397 240.1587 240.8778 241.5968 242.3158 243.0349 243.7539
+ 244.4729 245.1920 245.9110 246.6301 247.3491 248.0681 248.7872 249.5062 250.2253 250.9443
+ 251.6633 252.3824 253.1014 253.8204 254.5395 255.2585 255.9776 256.6966 257.4156 258.1347
+ 258.8537 259.5727 260.2918 261.0108 261.7299 262.4489 263.1679 263.8870 264.6060 265.3251
+ 266.0441 266.7631 267.4822 268.2012 268.9202 269.6393 270.3583 271.0774 271.7964 272.5154
+ 273.2345 273.9535 274.6725 275.3916 276.1106 276.8297 277.5487 278.2677 278.9868 279.7058
+ 280.4248 281.1439 281.8629 282.5820 283.3010 284.0200 284.7391 285.4581 286.1772 286.8962
+ 287.6152 288.3343 289.0533 289.7723 290.4914 291.2104 291.9295 292.6485 293.3675 294.0866
+ 294.8056 295.5246 296.2437 296.9627 297.6818 298.4008 299.1198 299.8389 300.5579 301.2770
+ 301.9960 302.7150 303.4341 304.1531 304.8721 305.5912 306.3102 307.0293 307.7483 308.4673
+ 309.1864 309.9054 310.6244 311.3435 312.0625 312.7816 313.5006 314.2196 314.9387 315.6577
+ 316.3768 317.0958 317.8148 318.5339 319.2529 319.9719 320.6910 321.4100 322.1291 322.8481
+ 323.5671 324.2862 325.0052 325.7242 326.4433 327.1623 327.8814 328.6004 329.3194 330.0385
+ 330.7575 331.4766 332.1956 332.9146 333.6337 334.3527 335.0717 335.7908 336.5098 337.2289
+ 337.9479 338.6669 339.3860 340.1050 340.8240 341.5431 342.2621 342.9812 343.7002 344.4192
+ 345.1383 345.8573 346.5764 347.2954 348.0144 348.7335 349.4525 350.1715 350.8906 351.6096
+ 352.3287 353.0477 353.7667 354.4858 355.2048 355.9238 356.6429 357.3619 358.0810 358.8000
+
+
+ <#> meas perf
+ projection data are noiseless
+
+
+ <#> back 0.0000
+ at levels
+ 511
+ background absorption 0.0000
+
+ estimate of totlen = 20962469.574466
+ estimate of totden = 9803802.737080
+ estimate of average density = 0.4677
+ projection data read
+ EXAMPLE 11 Brain Phantom
+ 0.094 seconds used for processing command proj
+
+
+ <*>
+
+ <#> STOP TERMINATION KLDS 100000 RPRT
+ termination test klds
+ reporting is enabled
+ reporting file: RPRTklds
+ reporting on every iteration
+ epsilon = 100000
+ 0.000 seconds used for processing command stop
+
+
+ <*>
+
+ <#> SUPERIORIZE 16 0.999 1 SMOO RPRT
+ Superiorization is enabled
+ N = 16
+ a = 0.999
+ b = 1
+ secondary criterion: smoo
+ reporting is enabled
+ reporting file: RPRTsuperiorization
+ reporting on every iteration
+ 0.000 seconds used for processing command supe
+
+
+ <#> EXECUTE AVERAGE EMAP
+
+ Example 11 Illustrating the Superiorized MAP EM algorithm for PET
+
+ <#> gamma is 0 EVAL
+
+ -----------------------------------------------------------
+
+ maximum a-posteriori probability expectation maximization
+
+ gamma: 0.000
+ evaluation flag is set
+
+ -----------------------------------------------------------
+
+ *** WARNING - There exist negative values in the projection data ***
+ *** Negative values set to 0.0 ***
+ value of l: 15
+ value of phi before algorithm operator: 4.42622e-28
+ value of phi after algorithm operator: 2.39806
+ algorithm executed in iteration 1
+ 1.734 seconds for the execution of the algorithm
+****Error: Negative raysum in prjfil detected, KLDS is not applicable! Terminating algorithm execution.
+ reconstruction completed after iteration 1
+ 1.734 seconds for this iteration
+ 1.922 seconds used for processing command exec
+
+
+ <*>
+
+ <#> END
--- /dev/null
+ snark14.s170710 - A PICTURE RECONSTRUCTION PROGRAM
+
+
+ <*> *******EXAMPLE 11
+
+ <*> SUPERIORIZED MAP EM ALGORITHM FOR EMISSION TOMOGRAPHY. RECONSTRUCTION
+
+ <*> OF BRAIN PHANTOM. SIMULATING PET GEOMETRY WITH A RING OF 300 DETECTORS
+
+ <*> WITH EACH DETECTOR IN COINCIDENCE WITH 101 DETECTORS OPPOSITE IT.
+
+ <*> COEFFICIENT OF PENALTLY TERM IS SET TO 0
+
+ <*>
+
+ <#> CREATE
+
+ EXAMPLE 11 Brain Phantom
+
+
+ <#> SPECTRUM MONOCHROMATIC 511
+ energy spectrum is monochromatic at energy level 511
+
+
+ <#> OBJECTS
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 511
+
+ 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.9500 0.9500
+
+ 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 1.0000 1.0000
+
+ 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 1.0000 1.0000
+
+ 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.9500 0.9500
+
+ 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.9500 0.9500
+
+ 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 1.0000 1.0000
+
+ 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.9500 0.9500
+
+ 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 1.0000 1.0000
+
+ 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 1.0000 1.0000
+
+ 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.9500 0.9500
+
+ 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.9500 0.9500
+
+ 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 1.0000 1.0000
+
+ 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.9500 0.9500
+
+ 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 1.0000 1.0000
+
+ 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 1.0000 1.0000
+
+ 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.9500 0.9500
+
+ 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 1.0000 1.0000
+
+ 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.9500 0.9500
+
+ 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.9500 0.9500
+
+ 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 1.0000 1.0000
+
+ 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 1.0000 1.0000
+
+ 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.9500 0.9500
+
+ 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 1.0000 1.0000
+
+ 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.9500 0.9500
+
+ 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.9500 0.9500
+
+ 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 1.0000 1.0000
+
+ 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 1.0000 1.0000
+
+ scale factor multiplying object densities 0.5100
+
+ seed set to 1
+ inhomogeneity set to 0.0500
+
+ <#> PHANTOM AVERAGE 7
+
+ this run will generate a phantom
+ density in each pixel is obtained as the average of 7 x 7 points
+
+
+ <#> 381 PIXELS OF SIZE 0.4
+ picture size 381 x 381, pixel size 0.4000
+
+
+ <#> RAYSUM AVERAGE 1
+
+ this run will generate projection data
+ projection data are calculated by dividing each ray interval into 1 substrips
+
+ with aperture (substrip) weights 1
+
+
+ <#> GEOMETRY
+
+
+ <#> divergent arc 153 306
+ rays are divergent from point sources
+ source to origin distance 153.0000
+ the detectors lie on an arc with source to detector distance = 306.0000
+
+
+ <#> RAYS USER 301 DETECTOR SPACING 1.1
+ number of rays per projection 301
+ at detector spacing 1.1000
+
+
+ <#> ANGLES 500 EQUAL SPACING
+ total number of projections 500
+
+ projection angles 0.0000 0.7190 1.4381 2.1571 2.8762 3.5952 4.3142 5.0333 5.7523 6.4713
+ 7.1904 7.9094 8.6285 9.3475 10.0665 10.7856 11.5046 12.2236 12.9427 13.6617
+ 14.3808 15.0998 15.8188 16.5379 17.2569 17.9760 18.6950 19.4140 20.1331 20.8521
+ 21.5711 22.2902 23.0092 23.7283 24.4473 25.1663 25.8854 26.6044 27.3234 28.0425
+ 28.7615 29.4806 30.1996 30.9186 31.6377 32.3567 33.0758 33.7948 34.5138 35.2329
+ 35.9519 36.6709 37.3900 38.1090 38.8281 39.5471 40.2661 40.9852 41.7042 42.4232
+ 43.1423 43.8613 44.5804 45.2994 46.0184 46.7375 47.4565 48.1756 48.8946 49.6136
+ 50.3327 51.0517 51.7707 52.4898 53.2088 53.9279 54.6469 55.3659 56.0850 56.8040
+ 57.5230 58.2421 58.9611 59.6802 60.3992 61.1182 61.8373 62.5563 63.2754 63.9944
+ 64.7134 65.4325 66.1515 66.8705 67.5896 68.3086 69.0277 69.7467 70.4657 71.1848
+ 71.9038 72.6228 73.3419 74.0609 74.7800 75.4990 76.2180 76.9371 77.6561 78.3752
+ 79.0942 79.8132 80.5323 81.2513 81.9703 82.6894 83.4084 84.1275 84.8465 85.5655
+ 86.2846 87.0036 87.7226 88.4417 89.1607 89.8798 90.5988 91.3178 92.0369 92.7559
+ 93.4749 94.1940 94.9130 95.6321 96.3511 97.0701 97.7892 98.5082 99.2273 99.9463
+ 100.6653 101.3844 102.1034 102.8224 103.5415 104.2605 104.9796 105.6986 106.4176 107.1367
+ 107.8557 108.5747 109.2938 110.0128 110.7319 111.4509 112.1699 112.8890 113.6080 114.3271
+ 115.0461 115.7651 116.4842 117.2032 117.9222 118.6413 119.3603 120.0794 120.7984 121.5174
+ 122.2365 122.9555 123.6745 124.3936 125.1126 125.8317 126.5507 127.2697 127.9888 128.7078
+ 129.4269 130.1459 130.8649 131.5840 132.3030 133.0220 133.7411 134.4601 135.1792 135.8982
+ 136.6172 137.3363 138.0553 138.7743 139.4934 140.2124 140.9315 141.6505 142.3695 143.0886
+ 143.8076 144.5267 145.2457 145.9647 146.6838 147.4028 148.1218 148.8409 149.5599 150.2790
+ 150.9980 151.7170 152.4361 153.1551 153.8741 154.5932 155.3122 156.0313 156.7503 157.4693
+ 158.1884 158.9074 159.6265 160.3455 161.0645 161.7836 162.5026 163.2216 163.9407 164.6597
+ 165.3788 166.0978 166.8168 167.5359 168.2549 168.9739 169.6930 170.4120 171.1311 171.8501
+ 172.5691 173.2882 174.0072 174.7263 175.4453 176.1643 176.8834 177.6024 178.3214 179.0405
+ 179.7595 180.4786 181.1976 181.9166 182.6357 183.3547 184.0737 184.7928 185.5118 186.2309
+ 186.9499 187.6689 188.3880 189.1070 189.8261 190.5451 191.2641 191.9832 192.7022 193.4212
+ 194.1403 194.8593 195.5784 196.2974 197.0164 197.7355 198.4545 199.1735 199.8926 200.6116
+ 201.3307 202.0497 202.7687 203.4878 204.2068 204.9259 205.6449 206.3639 207.0830 207.8020
+ 208.5210 209.2401 209.9591 210.6782 211.3972 212.1162 212.8353 213.5543 214.2733 214.9924
+ 215.7114 216.4305 217.1495 217.8685 218.5876 219.3066 220.0257 220.7447 221.4637 222.1828
+ 222.9018 223.6208 224.3399 225.0589 225.7780 226.4970 227.2160 227.9351 228.6541 229.3731
+ 230.0922 230.8112 231.5303 232.2493 232.9683 233.6874 234.4064 235.1255 235.8445 236.5635
+ 237.2826 238.0016 238.7206 239.4397 240.1587 240.8778 241.5968 242.3158 243.0349 243.7539
+ 244.4729 245.1920 245.9110 246.6301 247.3491 248.0681 248.7872 249.5062 250.2253 250.9443
+ 251.6633 252.3824 253.1014 253.8204 254.5395 255.2585 255.9776 256.6966 257.4156 258.1347
+ 258.8537 259.5727 260.2918 261.0108 261.7299 262.4489 263.1679 263.8870 264.6060 265.3251
+ 266.0441 266.7631 267.4822 268.2012 268.9202 269.6393 270.3583 271.0774 271.7964 272.5154
+ 273.2345 273.9535 274.6725 275.3916 276.1106 276.8297 277.5487 278.2677 278.9868 279.7058
+ 280.4248 281.1439 281.8629 282.5820 283.3010 284.0200 284.7391 285.4581 286.1772 286.8962
+ 287.6152 288.3343 289.0533 289.7723 290.4914 291.2104 291.9295 292.6485 293.3675 294.0866
+ 294.8056 295.5246 296.2437 296.9627 297.6818 298.4008 299.1198 299.8389 300.5579 301.2770
+ 301.9960 302.7150 303.4341 304.1531 304.8721 305.5912 306.3102 307.0293 307.7483 308.4673
+ 309.1864 309.9054 310.6244 311.3435 312.0625 312.7816 313.5006 314.2196 314.9387 315.6577
+ 316.3768 317.0958 317.8148 318.5339 319.2529 319.9719 320.6910 321.4100 322.1291 322.8481
+ 323.5671 324.2862 325.0052 325.7242 326.4433 327.1623 327.8814 328.6004 329.3194 330.0385
+ 330.7575 331.4766 332.1956 332.9146 333.6337 334.3527 335.0717 335.7908 336.5098 337.2289
+ 337.9479 338.6669 339.3860 340.1050 340.8240 341.5431 342.2621 342.9812 343.7002 344.4192
+ 345.1383 345.8573 346.5764 347.2954 348.0144 348.7335 349.4525 350.1715 350.8906 351.6096
+ 352.3287 353.0477 353.7667 354.4858 355.2048 355.9238 356.6429 357.3619 358.0810 358.8000
+
+
+ <#> MEASUREMENT PERFECT
+ projection data are noiseless
+
+
+ <#> BACKGROUND 0.0
+ at levels
+ 511
+ background absorption 0.0000
+
+
+ <#> RUN
+ 0.214 seconds phantom creation
+ 0.604 seconds projection data creation
+ 0.820 seconds used for processing command crea
+
+
+ <*>
+
+ <#> PICTURE TEST
+
+ EXAMPLE 11 Brain Phantom
+
+
+ <#> spec mono 511
+ energy spectrum is monochromatic at energy level 511
+
+
+ <#> obje
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 511
+
+ 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
+
+ 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
+
+ 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
+
+ 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
+
+ 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
+
+ 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
+
+ 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
+
+ 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
+
+ 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
+
+ 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
+
+ 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
+
+ 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
+
+ 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
+
+ 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
+
+ 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
+
+ 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
+
+ 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
+
+ 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
+
+ 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
+
+ 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
+
+ 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
+
+ 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
+
+ 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
+
+ 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
+
+ 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
+
+ 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
+
+ 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
+
+ scale factor multiplying object densities 0.5100
+
+ seed set to 1
+ inhomogeneity set to 0.0500
+
+ <#> phan aver 7
+
+ density in each pixel is obtained as the average of 7 x 7 points
+
+
+ <#> pixe 381 size 0.4000
+ picture size 381 x 381, pixel size 0.4000
+
+ test picture read
+ EXAMPLE 11 Brain Phantom
+ 0.043 seconds used for processing command pict
+
+
+ <*>
+
+ <#> PROJECTION REAL
+
+ EXAMPLE 11 Brain Phantom
+
+
+ <#> spec mono 511
+ energy spectrum is monochromatic at energy level 511
+
+
+ <#> obje
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 511
+
+ 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845
+
+ 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100
+
+ 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100
+
+ 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845
+
+ 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845
+
+ 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100
+
+ 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845
+
+ 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100
+
+ 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100
+
+ 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845
+
+ 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845
+
+ 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100
+
+ 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845
+
+ 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100
+
+ 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100
+
+ 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845
+
+ 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100
+
+ 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845
+
+ 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845
+
+ 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100
+
+ 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100
+
+ 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845
+
+ 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100
+
+ 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845
+
+ 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845
+
+ 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100
+
+ 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100
+
+ scale factor multiplying object densities 0.5100
+
+ seed set to 1
+ inhomogeneity set to 0.0500
+
+ <#> rays aver 1
+
+ projection data are calculated by dividing each ray interval into 1 substrips
+
+ with aperture (substrip) weights 1
+
+
+ <#> geom
+
+
+ <#> dive arc source at 153.0000 det dist 306.0000
+ rays are divergent from point sources
+ source to origin distance 153.0000
+ the detectors lie on an arc with source to detector distance = 306.0000
+
+
+ <#> rays user 301 spacing 1.1000
+ number of rays per projection 301
+ snark computed number of rays 437
+ at detector spacing 1.1000
+
+
+ <#> angl 500
+ total number of projections 500
+
+
+ projection angles 0.0000 0.7190 1.4381 2.1571 2.8762 3.5952 4.3142 5.0333 5.7523 6.4713
+ 7.1904 7.9094 8.6285 9.3475 10.0665 10.7856 11.5046 12.2236 12.9427 13.6617
+ 14.3808 15.0998 15.8188 16.5379 17.2569 17.9760 18.6950 19.4140 20.1331 20.8521
+ 21.5711 22.2902 23.0092 23.7283 24.4473 25.1663 25.8854 26.6044 27.3234 28.0425
+ 28.7615 29.4806 30.1996 30.9186 31.6377 32.3567 33.0758 33.7948 34.5138 35.2329
+ 35.9519 36.6709 37.3900 38.1090 38.8281 39.5471 40.2661 40.9852 41.7042 42.4232
+ 43.1423 43.8613 44.5804 45.2994 46.0184 46.7375 47.4565 48.1756 48.8946 49.6136
+ 50.3327 51.0517 51.7707 52.4898 53.2088 53.9279 54.6469 55.3659 56.0850 56.8040
+ 57.5230 58.2421 58.9611 59.6802 60.3992 61.1182 61.8373 62.5563 63.2754 63.9944
+ 64.7134 65.4325 66.1515 66.8705 67.5896 68.3086 69.0277 69.7467 70.4657 71.1848
+ 71.9038 72.6228 73.3419 74.0609 74.7800 75.4990 76.2180 76.9371 77.6561 78.3752
+ 79.0942 79.8132 80.5323 81.2513 81.9703 82.6894 83.4084 84.1275 84.8465 85.5655
+ 86.2846 87.0036 87.7226 88.4417 89.1607 89.8798 90.5988 91.3178 92.0369 92.7559
+ 93.4749 94.1940 94.9130 95.6321 96.3511 97.0701 97.7892 98.5082 99.2273 99.9463
+ 100.6653 101.3844 102.1034 102.8224 103.5415 104.2605 104.9796 105.6986 106.4176 107.1367
+ 107.8557 108.5747 109.2938 110.0128 110.7319 111.4509 112.1699 112.8890 113.6080 114.3271
+ 115.0461 115.7651 116.4842 117.2032 117.9222 118.6413 119.3603 120.0794 120.7984 121.5174
+ 122.2365 122.9555 123.6745 124.3936 125.1126 125.8317 126.5507 127.2697 127.9888 128.7078
+ 129.4269 130.1459 130.8649 131.5840 132.3030 133.0220 133.7411 134.4601 135.1792 135.8982
+ 136.6172 137.3363 138.0553 138.7743 139.4934 140.2124 140.9315 141.6505 142.3695 143.0886
+ 143.8076 144.5267 145.2457 145.9647 146.6838 147.4028 148.1218 148.8409 149.5599 150.2790
+ 150.9980 151.7170 152.4361 153.1551 153.8741 154.5932 155.3122 156.0313 156.7503 157.4693
+ 158.1884 158.9074 159.6265 160.3455 161.0645 161.7836 162.5026 163.2216 163.9407 164.6597
+ 165.3788 166.0978 166.8168 167.5359 168.2549 168.9739 169.6930 170.4120 171.1311 171.8501
+ 172.5691 173.2882 174.0072 174.7263 175.4453 176.1643 176.8834 177.6024 178.3214 179.0405
+ 179.7595 180.4786 181.1976 181.9166 182.6357 183.3547 184.0737 184.7928 185.5118 186.2309
+ 186.9499 187.6689 188.3880 189.1070 189.8261 190.5451 191.2641 191.9832 192.7022 193.4212
+ 194.1403 194.8593 195.5784 196.2974 197.0164 197.7355 198.4545 199.1735 199.8926 200.6116
+ 201.3307 202.0497 202.7687 203.4878 204.2068 204.9259 205.6449 206.3639 207.0830 207.8020
+ 208.5210 209.2401 209.9591 210.6782 211.3972 212.1162 212.8353 213.5543 214.2733 214.9924
+ 215.7114 216.4305 217.1495 217.8685 218.5876 219.3066 220.0257 220.7447 221.4637 222.1828
+ 222.9018 223.6208 224.3399 225.0589 225.7780 226.4970 227.2160 227.9351 228.6541 229.3731
+ 230.0922 230.8112 231.5303 232.2493 232.9683 233.6874 234.4064 235.1255 235.8445 236.5635
+ 237.2826 238.0016 238.7206 239.4397 240.1587 240.8778 241.5968 242.3158 243.0349 243.7539
+ 244.4729 245.1920 245.9110 246.6301 247.3491 248.0681 248.7872 249.5062 250.2253 250.9443
+ 251.6633 252.3824 253.1014 253.8204 254.5395 255.2585 255.9776 256.6966 257.4156 258.1347
+ 258.8537 259.5727 260.2918 261.0108 261.7299 262.4489 263.1679 263.8870 264.6060 265.3251
+ 266.0441 266.7631 267.4822 268.2012 268.9202 269.6393 270.3583 271.0774 271.7964 272.5154
+ 273.2345 273.9535 274.6725 275.3916 276.1106 276.8297 277.5487 278.2677 278.9868 279.7058
+ 280.4248 281.1439 281.8629 282.5820 283.3010 284.0200 284.7391 285.4581 286.1772 286.8962
+ 287.6152 288.3343 289.0533 289.7723 290.4914 291.2104 291.9295 292.6485 293.3675 294.0866
+ 294.8056 295.5246 296.2437 296.9627 297.6818 298.4008 299.1198 299.8389 300.5579 301.2770
+ 301.9960 302.7150 303.4341 304.1531 304.8721 305.5912 306.3102 307.0293 307.7483 308.4673
+ 309.1864 309.9054 310.6244 311.3435 312.0625 312.7816 313.5006 314.2196 314.9387 315.6577
+ 316.3768 317.0958 317.8148 318.5339 319.2529 319.9719 320.6910 321.4100 322.1291 322.8481
+ 323.5671 324.2862 325.0052 325.7242 326.4433 327.1623 327.8814 328.6004 329.3194 330.0385
+ 330.7575 331.4766 332.1956 332.9146 333.6337 334.3527 335.0717 335.7908 336.5098 337.2289
+ 337.9479 338.6669 339.3860 340.1050 340.8240 341.5431 342.2621 342.9812 343.7002 344.4192
+ 345.1383 345.8573 346.5764 347.2954 348.0144 348.7335 349.4525 350.1715 350.8906 351.6096
+ 352.3287 353.0477 353.7667 354.4858 355.2048 355.9238 356.6429 357.3619 358.0810 358.8000
+
+
+ <#> meas perf
+ projection data are noiseless
+
+
+ <#> back 0.0000
+ at levels
+ 511
+ background absorption 0.0000
+
+ estimate of totlen = 20962469.574466
+ estimate of totden = 9804758.593222
+ estimate of average density = 0.4677
+ projection data read
+ EXAMPLE 11 Brain Phantom
+ 0.070 seconds used for processing command proj
+
+
+ <*>
+
+ <#> STOP TERMINATION KLDS 100000 RPRT
+ termination test klds
+ reporting is enabled
+ reporting file: RPRTklds
+ reporting on every iteration
+ epsilon = 100000
+ 0.000 seconds used for processing command stop
+
+
+ <*>
+
+ <#> SUPERIORIZE 16 0.999 1 SMOO RPRT
+ Superiorization is enabled
+ N = 16
+ a = 0.999
+ b = 1
+ secondary criterion: smoo
+ reporting is enabled
+ reporting file: RPRTsuperiorization
+ reporting on every iteration
+ 0.000 seconds used for processing command supe
+
+
+ <#> EXECUTE AVERAGE EMAP
+
+ Example 11 Illustrating the Superiorized MAP EM algorithm for PET
+
+ <#> gamma is 0 EVAL
+
+ -----------------------------------------------------------
+
+ maximum a-posteriori probability expectation maximization
+
+ gamma: 0.000
+ evaluation flag is set
+
+ -----------------------------------------------------------
+
+ *** WARNING - There exist negative values in the projection data ***
+ *** Negative values set to 0.0 ***
+ value of l: 15
+ value of phi before algorithm operator: 0
+ value of phi after algorithm operator: 2.39758
+ algorithm executed in iteration 1
+ 1.266 seconds for the execution of the algorithm
+****Error: Negative raysum in prjfil detected, KLDS is not applicable! Terminating algorithm execution.
+ reconstruction completed after iteration 1
+ 1.268 seconds for this iteration
+ 1.392 seconds used for processing command exec
+
+
+ <*>
+
+ <#> END
<#> RUN
- 0.014 seconds phantom creation
- 0.078 seconds projection data creation
- 0.093 seconds used for processing command crea
+ 0.031 seconds phantom creation
+ 0.109 seconds projection data creation
+ 0.141 seconds used for processing command crea
<*>
test picture read
EXAMPLE 7 Brain Phantom
- 0.003 seconds used for processing command pict
+ 0.000 seconds used for processing command pict
<*>
background absorption 0.0000
estimate of totlen = 4263169.769963
- estimate of totden = 2023083.983824
- estimate of average density = 0.4745
+ estimate of totden = 2020712.203355
+ estimate of average density = 0.4740
projection data read
EXAMPLE 7 Brain Phantom
- 0.017 seconds used for processing command proj
+ 0.016 seconds used for processing command proj
<*>
*** WARNING - There exist negative values in the projection data ***
*** Negative values set to 0.0 ***
algorithm executed in iteration 1
- 0.063 seconds for the execution of the algorithm
+ 0.172 seconds for the execution of the algorithm
iteration 1 completed
- 0.064 seconds for this iteration
+ 0.172 seconds for this iteration
algorithm executed in iteration 2
- 0.040 seconds for the execution of the algorithm
+ 0.062 seconds for the execution of the algorithm
iteration 2 completed
- 0.040 seconds for this iteration
+ 0.062 seconds for this iteration
algorithm executed in iteration 3
- 0.040 seconds for the execution of the algorithm
+ 0.031 seconds for the execution of the algorithm
iteration 3 completed
- 0.040 seconds for this iteration
+ 0.031 seconds for this iteration
algorithm executed in iteration 4
- 0.040 seconds for the execution of the algorithm
+ 0.094 seconds for the execution of the algorithm
iteration 4 completed
- 0.040 seconds for this iteration
+ 0.094 seconds for this iteration
algorithm executed in iteration 5
- 0.041 seconds for the execution of the algorithm
+ 0.094 seconds for the execution of the algorithm
iteration 5 completed
- 0.041 seconds for this iteration
+ 0.094 seconds for this iteration
algorithm executed in iteration 6
- 0.042 seconds for the execution of the algorithm
+ 0.047 seconds for the execution of the algorithm
iteration 6 completed
- 0.042 seconds for this iteration
+ 0.047 seconds for this iteration
algorithm executed in iteration 7
- 0.044 seconds for the execution of the algorithm
+ 0.109 seconds for the execution of the algorithm
iteration 7 completed
- 0.044 seconds for this iteration
+ 0.109 seconds for this iteration
algorithm executed in iteration 8
- 0.041 seconds for the execution of the algorithm
+ 0.047 seconds for the execution of the algorithm
iteration 8 completed
- 0.041 seconds for this iteration
+ 0.047 seconds for this iteration
algorithm executed in iteration 9
- 0.040 seconds for the execution of the algorithm
+ 0.031 seconds for the execution of the algorithm
iteration 9 completed
- 0.040 seconds for this iteration
+ 0.031 seconds for this iteration
algorithm executed in iteration 10
- 0.041 seconds for the execution of the algorithm
+ 0.031 seconds for the execution of the algorithm
iteration 10 completed
- 0.041 seconds for this iteration
+ 0.031 seconds for this iteration
algorithm executed in iteration 11
- 0.040 seconds for the execution of the algorithm
+ 0.062 seconds for the execution of the algorithm
iteration 11 completed
- 0.040 seconds for this iteration
+ 0.062 seconds for this iteration
algorithm executed in iteration 12
- 0.040 seconds for the execution of the algorithm
+ 0.094 seconds for the execution of the algorithm
iteration 12 completed
- 0.040 seconds for this iteration
+ 0.094 seconds for this iteration
algorithm executed in iteration 13
- 0.040 seconds for the execution of the algorithm
+ 0.062 seconds for the execution of the algorithm
iteration 13 completed
- 0.040 seconds for this iteration
+ 0.062 seconds for this iteration
algorithm executed in iteration 14
- 0.040 seconds for the execution of the algorithm
+ 0.031 seconds for the execution of the algorithm
iteration 14 completed
- 0.040 seconds for this iteration
+ 0.031 seconds for this iteration
algorithm executed in iteration 15
- 0.040 seconds for the execution of the algorithm
+ 0.062 seconds for the execution of the algorithm
iteration 15 completed
- 0.040 seconds for this iteration
+ 0.062 seconds for this iteration
algorithm executed in iteration 16
- 0.040 seconds for the execution of the algorithm
+ 0.078 seconds for the execution of the algorithm
iteration 16 completed
- 0.041 seconds for this iteration
+ 0.078 seconds for this iteration
algorithm executed in iteration 17
- 0.040 seconds for the execution of the algorithm
+ 0.078 seconds for the execution of the algorithm
iteration 17 completed
- 0.040 seconds for this iteration
+ 0.078 seconds for this iteration
algorithm executed in iteration 18
- 0.040 seconds for the execution of the algorithm
+ 0.062 seconds for the execution of the algorithm
iteration 18 completed
- 0.040 seconds for this iteration
+ 0.062 seconds for this iteration
algorithm executed in iteration 19
- 0.041 seconds for the execution of the algorithm
+ 0.078 seconds for the execution of the algorithm
iteration 19 completed
- 0.041 seconds for this iteration
+ 0.078 seconds for this iteration
algorithm executed in iteration 20
- 0.041 seconds for the execution of the algorithm
+ 0.062 seconds for the execution of the algorithm
reconstruction completed after iteration 20
- 0.041 seconds for this iteration
- 0.836 seconds for all iterations
- 0.883 seconds used for processing command exec
+ 0.062 seconds for this iteration
+ 1.391 seconds for all iterations
+ 1.453 seconds used for processing command exec
<*>
#!/bin/bash
# Exclude examples 7 and 11 as they use random number generator not compatible across platforms
+# Exclude example 12 as small differences in phi across gcc5.4 vs gcc7.2
-for i in 1 2 3 4 5 6 8 9 10 12 13;
+for i in 1 2 3 4 5 6 8 9 10 13;
do
echo Computing the diff for example b"$i".
cd b"$i"
projects / snark14.git / commitdiff