examples/b11/b11r.out

   1    snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM
   2
   3
   4      <*> *******EXAMPLE 11
   5
   6      <*>  SUPERIORIZED MAP EM ALGORITHM FOR EMISSION TOMOGRAPHY. RECONSTRUCTION
   7
   8      <*>  OF BRAIN PHANTOM. SIMULATING PET GEOMETRY WITH A RING OF 300 DETECTORS
   9
  10      <*>  WITH EACH DETECTOR IN COINCIDENCE WITH 101 DETECTORS OPPOSITE IT.
  11
  12      <*>  COEFFICIENT OF PENALTLY TERM IS SET TO 0
  13
  14      <*>
  15
  16      <#> CREATE
  17
  18          EXAMPLE 11 Brain Phantom
  19
  20
  21      <#> SPECTRUM MONOCHROMATIC 511
  22          energy spectrum is monochromatic at energy level   511
  23
  24
  25      <#> OBJECTS
  26          description of objects
  27                                                                 density at levels
  28  numb type  x-coord  y-coord x-length y-length    angle  av dens       511
  29
  30     1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.9500    0.9500
  31
  32     2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   1.0000    1.0000
  33
  34     3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   1.0000    1.0000
  35
  36     4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.9500    0.9500
  37
  38     5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.9500    0.9500
  39
  40     6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   1.0000    1.0000
  41
  42     7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.9500    0.9500
  43
  44     8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   1.0000    1.0000
  45
  46     9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   1.0000    1.0000
  47
  48    10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.9500    0.9500
  49
  50    11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.9500    0.9500
  51
  52    12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   1.0000    1.0000
  53
  54    13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.9500    0.9500
  55
  56    14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   1.0000    1.0000
  57
  58    15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   1.0000    1.0000
  59
  60    16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.9500    0.9500
  61
  62    17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   1.0000    1.0000
  63
  64    18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.9500    0.9500
  65
  66    19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.9500    0.9500
  67
  68    20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   1.0000    1.0000
  69
  70    21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   1.0000    1.0000
  71
  72    22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.9500    0.9500
  73
  74    23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   1.0000    1.0000
  75
  76    24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.9500    0.9500
  77
  78    25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.9500    0.9500
  79
  80    26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   1.0000    1.0000
  81
  82    27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   1.0000    1.0000
  83
  84          scale factor multiplying object densities     0.5100
  85
  86          seed set to 1
  87          inhomogeneity set to     0.0500
  88
  89      <#> PHANTOM AVERAGE 7
  90
  91          this run will generate a phantom
  92          density in each pixel is obtained as the average of 7 x 7 points
  93
  94
  95      <#> 381 PIXELS OF SIZE 0.4
  96          picture size 381 x 381,  pixel size     0.4000
  97
  98
  99      <#> RAYSUM AVERAGE 1
 100
 101          this run will generate projection data
 102          projection data are calculated by dividing each ray interval into 1 substrips
 103
 104          with aperture (substrip) weights     1
 105
 106
 107      <#> GEOMETRY
 108
 109
 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
 114
 115
 116      <#> RAYS USER 301 DETECTOR SPACING 1.1
 117          number of rays per projection   301
 118          at detector spacing     1.1000
 119
 120
 121      <#> ANGLES 500 EQUAL SPACING
 122          total number of projections   500
 123
 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
 174
 175
 176      <#> MEASUREMENT NOISY
 177          noise characteristics of projection data follow
 178                nature          characteristics
 179
 180      <#> QUANTUM 1.0 1.0 CALIBRATION 4
 181               Emission tomography
 182
 183      <#> SEED 0
 184           seed for random number generator is              0
 185
 186
 187      <#> BACKGROUND 0.0
 188                                at levels
 189                                      511
 190           background absorption   0.0000
 191
 192
 193      <#> RUN
 194               0.358 seconds phantom creation
 195               1.289 seconds projection data creation
 196               1.647 seconds used for processing command crea
 197
 198
 199      <*>
 200
 201      <#> PICTURE TEST
 202
 203          EXAMPLE 11 Brain Phantom
 204
 205
 206      <#> spec    mono  511
 207          energy spectrum is monochromatic at energy level   511
 208
 209
 210      <#> obje
 211          description of objects
 212                                                                 density at levels
 213  numb type  x-coord  y-coord x-length y-length    angle  av dens       511
 214
 215     1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.4845    0.4845
 216
 217     2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   0.5100    0.5100
 218
 219     3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   0.5100    0.5100
 220
 221     4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.4845    0.4845
 222
 223     5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.4845    0.4845
 224
 225     6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   0.5100    0.5100
 226
 227     7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.4845    0.4845
 228
 229     8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   0.5100    0.5100
 230
 231     9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   0.5100    0.5100
 232
 233    10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.4845    0.4845
 234
 235    11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.4845    0.4845
 236
 237    12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   0.5100    0.5100
 238
 239    13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.4845    0.4845
 240
 241    14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   0.5100    0.5100
 242
 243    15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   0.5100    0.5100
 244
 245    16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.4845    0.4845
 246
 247    17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   0.5100    0.5100
 248
 249    18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.4845    0.4845
 250
 251    19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.4845    0.4845
 252
 253    20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   0.5100    0.5100
 254
 255    21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   0.5100    0.5100
 256
 257    22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.4845    0.4845
 258
 259    23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   0.5100    0.5100
 260
 261    24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.4845    0.4845
 262
 263    25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.4845    0.4845
 264
 265    26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   0.5100    0.5100
 266
 267    27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   0.5100    0.5100
 268
 269          scale factor multiplying object densities     0.5100
 270
 271          seed set to 1
 272          inhomogeneity set to     0.0500
 273
 274      <#> phan    aver    7
 275
 276          density in each pixel is obtained as the average of 7 x 7 points
 277
 278
 279      <#> pixe      381    size        0.4000
 280          picture size 381 x 381,  pixel size     0.4000
 281
 282          test picture read
 283          EXAMPLE 11 Brain Phantom
 284               0.042 seconds used for processing command pict
 285
 286
 287      <*>
 288
 289      <#> PROJECTION REAL
 290
 291          EXAMPLE 11 Brain Phantom
 292
 293
 294      <#> spec    mono  511
 295          energy spectrum is monochromatic at energy level   511
 296
 297
 298      <#> obje
 299          description of objects
 300                                                                 density at levels
 301  numb type  x-coord  y-coord x-length y-length    angle  av dens       511
 302
 303     1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.4845    0.4845
 304
 305     2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   0.5100    0.5100
 306
 307     3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   0.5100    0.5100
 308
 309     4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.4845    0.4845
 310
 311     5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.4845    0.4845
 312
 313     6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   0.5100    0.5100
 314
 315     7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.4845    0.4845
 316
 317     8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   0.5100    0.5100
 318
 319     9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   0.5100    0.5100
 320
 321    10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.4845    0.4845
 322
 323    11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.4845    0.4845
 324
 325    12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   0.5100    0.5100
 326
 327    13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.4845    0.4845
 328
 329    14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   0.5100    0.5100
 330
 331    15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   0.5100    0.5100
 332
 333    16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.4845    0.4845
 334
 335    17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   0.5100    0.5100
 336
 337    18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.4845    0.4845
 338
 339    19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.4845    0.4845
 340
 341    20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   0.5100    0.5100
 342
 343    21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   0.5100    0.5100
 344
 345    22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.4845    0.4845
 346
 347    23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   0.5100    0.5100
 348
 349    24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.4845    0.4845
 350
 351    25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.4845    0.4845
 352
 353    26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   0.5100    0.5100
 354
 355    27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   0.5100    0.5100
 356
 357          scale factor multiplying object densities     0.5100
 358
 359          seed set to 1
 360          inhomogeneity set to     0.0500
 361
 362      <#> rays    aver    1
 363
 364          projection data are calculated by dividing each ray interval into 1 substrips
 365
 366          with aperture (substrip) weights     1
 367
 368
 369      <#> geom
 370
 371
 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
 376
 377
 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
 382
 383
 384      <#> angl      500
 385          total number of projections   500
 386
 387
 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
 438
 439
 440      <#> meas    nois
 441          noise characteristics of projection data follow
 442                nature          characteristics
 443
 444      <#> quan            1.0000            1.0000    cali  4
 445               Emission tomography
 446
 447      <#> seed       0
 448           seed for random number generator is              0
 449
 450
 451      <#> back      0.0000
 452                                at levels
 453                                      511
 454           background absorption   0.0000
 455
 456          estimate of totlen =  20962469.574466
 457          estimate of totden =   9802000.000000
 458          estimate of average density =     0.4676
 459          projection data read
 460          EXAMPLE 11 Brain Phantom
 461               0.054 seconds used for processing command proj
 462
 463
 464      <*>
 465
 466      <#> STOP TERMINATION KLDS 100000 RPRT
 467          termination test klds
 468          reporting is enabled
 469          reporting file: RPRTklds
 470          reporting on every iteration
 471          epsilon = 100000
 472               0.000 seconds used for processing command stop
 473
 474
 475      <*>
 476
 477      <#> SUPERIORIZE 16 0.999 1 SMOO RPRT
 478          Superiorization is enabled
 479          N = 16
 480          a = 0.999
 481          b = 1
 482          secondary criterion: smoo
 483          reporting is enabled
 484          reporting file: RPRTsuperiorization
 485          reporting on every iteration
 486               0.000 seconds used for processing command supe
 487
 488
 489      <#> EXECUTE AVERAGE EMAP
 490
 491          Example 11 Illustrating the Superiorized MAP EM algorithm for PET
 492
 493      <#> gamma is 0 EVAL
 494
 495   -----------------------------------------------------------
 496
 497    maximum a-posteriori probability expectation maximization
 498
 499             gamma:   0.000
 500             evaluation flag is set
 501
 502   -----------------------------------------------------------
 503
 504               value of l: 15
 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
 512               value of l: 31
 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
 520               value of l: 47
 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
 528               value of l: 63
 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
 536               value of l: 79
 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
 544               value of l: 95
 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
 552               value of l: 111
 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
 560               value of l: 127
 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
 568               value of l: 143
 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
 576               value of l: 159
 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
 584               value of l: 175
 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
 592               value of l: 191
 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
 602
 603
 604      <*>
 605
 606      <#> END