examples/b1/b1r.out

   1    snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM
   2
   3
   4      <*>  EXAMPLE 1
   5
   6      <*>  DAISY PATTERN
   7
   8      <*>  CREATION OF PHANTOM AND DIVERGENT PROJECTION DATA
   9
  10      <*>  RECONSTRUCTION WITH DCON
  11
  12      <*>  ANALYSIS
  13
  14      <#> CREATE
  15
  16          EXAMPLE 1 DAISY
  17
  18
  19      <#> SPECTRUM MONOCHROMATIC 75
  20          energy spectrum is monochromatic at energy level    75
  21
  22
  23      <#> OBJECTS
  24          description of objects
  25                                                                 density at levels
  26  numb type  x-coord  y-coord x-length y-length    angle  av dens        75
  27
  28     1 tria   1.9000  -3.7000   0.5000   4.1000  20.0000   2.5000    2.5000
  29
  30     2 segm   4.0000   0.1000   2.9000   4.3000  47.0000   2.2000    2.2000
  31
  32     3 segm  -2.1000  -2.1000   3.1000   3.0000 -15.0000   2.2000    2.2000
  33
  34     4 elip   0.0000   2.5000   1.0000   1.0000   0.0000   1.5000    1.5000
  35
  36     5 elip   0.1000   4.3000   1.3000   1.2000   0.0000   1.0000    1.0000
  37
  38     6 elip  -2.0000   3.5000   1.3000   1.3000   0.0000   1.0000    1.0000
  39
  40     7 elip  -1.0000   1.5000   1.3000   1.1000   0.0000   1.0000    1.0000
  41
  42     8 elip   1.5000   1.6000   1.5000   1.3000   0.0000   1.0000    1.0000
  43
  44     9 elip   1.8000   3.3000   1.2000   1.2000   0.0000   1.0000    1.0000
  45
  46          scale factor multiplying object densities     1.0000
  47
  48          seed set to 0
  49          inhomogeneity set to     0.0000
  50
  51      <#> PHANTOM AVERAGE 3
  52
  53          this run will generate a phantom
  54          density in each pixel is obtained as the average of 3 x 3 points
  55
  56
  57      <#> 31 0.4
  58          picture size 31 x 31,  pixel size     0.4000
  59
  60
  61      <#> RAYSUM AVERAGE 1
  62
  63          this run will generate projection data
  64          projection data are calculated by dividing each ray interval into 1 substrips
  65
  66          with aperture (substrip) weights     1
  67
  68
  69      <#> GEOMETRY
  70
  71
  72      <#> DIVERGENT ARC 20.0 10000.0
  73          rays are divergent from point sources
  74          source to origin distance     20.0000
  75          the detectors lie on an arc with source to detector distance = 10000.0000
  76
  77
  78      <#> RAYS USER 65 150.0
  79          number of rays per projection    65
  80          at detector spacing   150.0000
  81
  82
  83      <#> ANGLES 60 EQUALLY SPACED
  84          total number of projections    60
  85
  86          projection angles    0.0000    5.9746   11.9492   17.9237   23.8983   29.8729   35.8475   41.8220   47.7966   53.7712
  87                              59.7458   65.7203   71.6949   77.6695   83.6441   89.6186   95.5932  101.5678  107.5424  113.5169
  88                             119.4915  125.4661  131.4407  137.4153  143.3898  149.3644  155.3390  161.3136  167.2881  173.2627
  89                             179.2373  185.2119  191.1864  197.1610  203.1356  209.1102  215.0847  221.0593  227.0339  233.0085
  90                             238.9831  244.9576  250.9322  256.9068  262.8814  268.8559  274.8305  280.8051  286.7797  292.7542
  91                             298.7288  304.7034  310.6780  316.6525  322.6271  328.6017  334.5763  340.5508  346.5254  352.5000
  92
  93
  94      <#> MEASUREMENT PERFECT
  95          projection data are noiseless
  96
  97
  98      <#> BACKGROUND 0.0
  99                                at levels
 100                                       75
 101           background absorption   0.0000
 102
 103
 104      <#> RUN
 105               0.001 seconds phantom creation
 106               0.004 seconds projection data creation
 107               0.005 seconds used for processing command crea
 108
 109
 110      <*>
 111
 112      <#> PICTURE TEST
 113
 114          EXAMPLE 1 DAISY
 115
 116
 117      <#> spec    mono   75
 118          energy spectrum is monochromatic at energy level    75
 119
 120
 121      <#> obje
 122          description of objects
 123                                                                 density at levels
 124  numb type  x-coord  y-coord x-length y-length    angle  av dens        75
 125
 126     1 tria   1.9000  -3.7000   0.5000   4.1000  20.0000   2.5000    2.5000
 127
 128     2 segm   4.0000   0.1000   2.9000   4.3000  47.0000   2.2000    2.2000
 129
 130     3 segm  -2.1000  -2.1000   3.1000   3.0000 -15.0000   2.2000    2.2000
 131
 132     4 elip   0.0000   2.5000   1.0000   1.0000   0.0000   1.5000    1.5000
 133
 134     5 elip   0.1000   4.3000   1.3000   1.2000   0.0000   1.0000    1.0000
 135
 136     6 elip  -2.0000   3.5000   1.3000   1.3000   0.0000   1.0000    1.0000
 137
 138     7 elip  -1.0000   1.5000   1.3000   1.1000   0.0000   1.0000    1.0000
 139
 140     8 elip   1.5000   1.6000   1.5000   1.3000   0.0000   1.0000    1.0000
 141
 142     9 elip   1.8000   3.3000   1.2000   1.2000   0.0000   1.0000    1.0000
 143
 144          scale factor multiplying object densities     1.0000
 145
 146          seed set to 0
 147          inhomogeneity set to     0.0000
 148
 149      <#> phan    aver    3
 150
 151          density in each pixel is obtained as the average of 3 x 3 points
 152
 153
 154      <#> pixe       31    size        0.4000
 155          picture size 31 x 31,  pixel size     0.4000
 156
 157          test picture read
 158          EXAMPLE 1 DAISY
 159               0.000 seconds used for processing command pict
 160
 161
 162      <*>
 163
 164      <#> PROJECTION REAL
 165
 166          EXAMPLE 1 DAISY
 167
 168
 169      <#> spec    mono   75
 170          energy spectrum is monochromatic at energy level    75
 171
 172
 173      <#> obje
 174          description of objects
 175                                                                 density at levels
 176  numb type  x-coord  y-coord x-length y-length    angle  av dens        75
 177
 178     1 tria   1.9000  -3.7000   0.5000   4.1000  20.0000   2.5000    2.5000
 179
 180     2 segm   4.0000   0.1000   2.9000   4.3000  47.0000   2.2000    2.2000
 181
 182     3 segm  -2.1000  -2.1000   3.1000   3.0000 -15.0000   2.2000    2.2000
 183
 184     4 elip   0.0000   2.5000   1.0000   1.0000   0.0000   1.5000    1.5000
 185
 186     5 elip   0.1000   4.3000   1.3000   1.2000   0.0000   1.0000    1.0000
 187
 188     6 elip  -2.0000   3.5000   1.3000   1.3000   0.0000   1.0000    1.0000
 189
 190     7 elip  -1.0000   1.5000   1.3000   1.1000   0.0000   1.0000    1.0000
 191
 192     8 elip   1.5000   1.6000   1.5000   1.3000   0.0000   1.0000    1.0000
 193
 194     9 elip   1.8000   3.3000   1.2000   1.2000   0.0000   1.0000    1.0000
 195
 196          scale factor multiplying object densities     1.0000
 197
 198          seed set to 0
 199          inhomogeneity set to     0.0000
 200
 201      <#> rays    aver    1
 202
 203          projection data are calculated by dividing each ray interval into 1 substrips
 204
 205          with aperture (substrip) weights     1
 206
 207
 208      <#> geom
 209
 210
 211      <#> dive    arc     source at   20.0000     det dist10000.0000
 212          rays are divergent from point sources
 213          source to origin distance     20.0000
 214          the detectors lie on an arc with source to detector distance = 10000.0000
 215
 216
 217      <#> rays    user       65    spacing      150.0000
 218          number of rays per projection    65
 219          snark computed number of rays    63
 220          at detector spacing   150.0000
 221
 222
 223      <#> angl       60
 224          total number of projections    60
 225
 226
 227          projection angles    0.0000    5.9746   11.9492   17.9237   23.8983   29.8729   35.8475   41.8220   47.7966   53.7712
 228                              59.7458   65.7203   71.6949   77.6695   83.6441   89.6186   95.5932  101.5678  107.5424  113.5169
 229                             119.4915  125.4661  131.4407  137.4153  143.3898  149.3644  155.3390  161.3136  167.2881  173.2627
 230                             179.2373  185.2119  191.1864  197.1610  203.1356  209.1102  215.0847  221.0593  227.0339  233.0085
 231                             238.9831  244.9576  250.9322  256.9068  262.8814  268.8559  274.8305  280.8051  286.7797  292.7542
 232                             298.7288  304.7034  310.6780  316.6525  322.6271  328.6017  334.5763  340.5508  346.5254  352.5000
 233
 234
 235      <#> meas    perf
 236          projection data are noiseless
 237
 238
 239      <#> back      0.0000
 240                                at levels
 241                                       75
 242           background absorption   0.0000
 243
 244          estimate of totlen =     31273.639204
 245          estimate of totden =     11129.885164
 246          estimate of average density =     0.3559
 247          projection data read
 248          EXAMPLE 1 DAISY
 249               0.000 seconds used for processing command proj
 250
 251
 252      <*>
 253
 254      <#> EXECUTE DCON
 255
 256          EXAMPLE 1 DIVERGENT CONVOLUTION FROM 60 VIEWS
 257
 258      <#> -1 1 -1 0.5 0
 259           cubic spline interpolation
 260           one out of every     1 projection(s) used for reconstruction
 261            64 points used on each side of convolution
 262           pre-smoothing weight on center ray = 0.5000
 263           pre-smoothing weight on side rays  = 0.2500
 264           accurate back-projection
 265
 266      <#> BANDLIMITING
 267           band-limiting filter
 268           total time for obtaining smoothed projection data was     0.001
 269           total time for convoluting projection data was            0.000
 270           total time for back-projecting was                        0.001
 271          algorithm executed in iteration    1
 272               0.002 seconds for the execution of the algorithm
 273          reconstruction completed after iteration    1
 274               0.002 seconds for this iteration
 275               0.002 seconds used for processing command exec
 276
 277
 278      <*>
 279
 280      <#> LINES SCALE 100.0 COLUMN 16
 281
 282          last iteration
 283
 284      lines display of reconstruction using DCON iter 1
 285
 286      phantom name:    EXAMPLE 1 DAISY
 287      projection data: EXAMPLE 1 DAISY
 288      execution name:  EXAMPLE 1 DIVERGENT CONVOLUTION FROM 60 VIEWS
 289      scaling factor =  100.00000
 290
 291  column =              16
 292  row    original   estimate difference
 293    0     0.00000   -7.80005   -7.80005
 294    1    22.22222   28.11589    5.89367
 295    2   100.00000   88.37252  -11.62748
 296    3   100.00000  104.43894    4.43894
 297    4   100.00000  102.69494    2.69494
 298    5   100.00000  107.73658    7.73658
 299    6   100.00000  147.49863   47.49863
 300    7   216.66667  188.04080  -28.62586
 301    8   150.00000  180.51891   30.51891
 302    9   216.66667  211.83020   -4.83647
 303   10   250.00000  245.52265   -4.47735
 304   11   216.66667  199.38366  -17.28301
 305   12   111.11111  120.98082    9.86971
 306   13    66.66667   68.05861    1.39194
 307   14     0.00000   13.57312   13.57312
 308   15    27.77778   12.48964  -15.28813
 309   16     0.00000   27.99191   27.99191
 310   17     0.00000   24.69207   24.69207
 311   18     0.00000   10.47314   10.47314
 312   19     0.00000   -0.41472   -0.41472
 313   20     0.00000    0.60572    0.60572
 314   21     0.00000   10.05671   10.05671
 315   22   122.22222  102.66325  -19.55897
 316   23    97.77778   87.82826   -9.94952
 317   24     0.00000    5.13094    5.13094
 318   25     0.00000   -2.85855   -2.85855
 319   26     0.00000   -3.34366   -3.34366
 320   27     0.00000    1.38876    1.38876
 321   28     0.00000    3.44155    3.44155
 322   29     0.00000    5.01602    5.01602
 323   30     0.00000   -0.04334   -0.04334
 324
 325
 326
 327               0.001 seconds used for processing command line
 328
 329
 330      <*>
 331
 332      <#> PUNCH PHANTOM
 333
 334          last iteration
 335
 336           reconstruction of EXAMPLE 1 DIVERGENT CONVOLUTIO using DCON iter     1,   EXAMPLE 1 DAISY
 337               0.000 seconds used for processing command punc
 338
 339
 340      <*>
 341
 342      <#> END