--- /dev/null
+ snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM
+
+
+ <*> EXAMPLE 1
+
+ <*> DAISY PATTERN
+
+ <*> CREATION OF PHANTOM AND DIVERGENT PROJECTION DATA
+
+ <*> RECONSTRUCTION WITH DCON
+
+ <*> ANALYSIS
+
+ <#> CREATE
+
+ EXAMPLE 1 DAISY
+
+
+ <#> SPECTRUM MONOCHROMATIC 75
+ energy spectrum is monochromatic at energy level 75
+
+
+ <#> OBJECTS
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 75
+
+ 1 tria 1.9000 -3.7000 0.5000 4.1000 20.0000 2.5000 2.5000
+
+ 2 segm 4.0000 0.1000 2.9000 4.3000 47.0000 2.2000 2.2000
+
+ 3 segm -2.1000 -2.1000 3.1000 3.0000 -15.0000 2.2000 2.2000
+
+ 4 elip 0.0000 2.5000 1.0000 1.0000 0.0000 1.5000 1.5000
+
+ 5 elip 0.1000 4.3000 1.3000 1.2000 0.0000 1.0000 1.0000
+
+ 6 elip -2.0000 3.5000 1.3000 1.3000 0.0000 1.0000 1.0000
+
+ 7 elip -1.0000 1.5000 1.3000 1.1000 0.0000 1.0000 1.0000
+
+ 8 elip 1.5000 1.6000 1.5000 1.3000 0.0000 1.0000 1.0000
+
+ 9 elip 1.8000 3.3000 1.2000 1.2000 0.0000 1.0000 1.0000
+
+ scale factor multiplying object densities 1.0000
+
+ seed set to 0
+ inhomogeneity set to 0.0000
+
+ <#> PHANTOM AVERAGE 3
+
+ this run will generate a phantom
+ density in each pixel is obtained as the average of 3 x 3 points
+
+
+ <#> 31 0.4
+ picture size 31 x 31, 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 20.0 10000.0
+ rays are divergent from point sources
+ source to origin distance 20.0000
+ the detectors lie on an arc with source to detector distance = 10000.0000
+
+
+ <#> RAYS USER 65 150.0
+ number of rays per projection 65
+ at detector spacing 150.0000
+
+
+ <#> ANGLES 60 EQUALLY SPACED
+ total number of projections 60
+
+ projection angles 0.0000 5.9746 11.9492 17.9237 23.8983 29.8729 35.8475 41.8220 47.7966 53.7712
+ 59.7458 65.7203 71.6949 77.6695 83.6441 89.6186 95.5932 101.5678 107.5424 113.5169
+ 119.4915 125.4661 131.4407 137.4153 143.3898 149.3644 155.3390 161.3136 167.2881 173.2627
+ 179.2373 185.2119 191.1864 197.1610 203.1356 209.1102 215.0847 221.0593 227.0339 233.0085
+ 238.9831 244.9576 250.9322 256.9068 262.8814 268.8559 274.8305 280.8051 286.7797 292.7542
+ 298.7288 304.7034 310.6780 316.6525 322.6271 328.6017 334.5763 340.5508 346.5254 352.5000
+
+
+ <#> MEASUREMENT PERFECT
+ projection data are noiseless
+
+
+ <#> BACKGROUND 0.0
+ at levels
+ 75
+ background absorption 0.0000
+
+
+ <#> RUN
+ 0.001 seconds phantom creation
+ 0.004 seconds projection data creation
+ 0.005 seconds used for processing command crea
+
+
+ <*>
+
+ <#> PICTURE TEST
+
+ EXAMPLE 1 DAISY
+
+
+ <#> spec mono 75
+ energy spectrum is monochromatic at energy level 75
+
+
+ <#> obje
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 75
+
+ 1 tria 1.9000 -3.7000 0.5000 4.1000 20.0000 2.5000 2.5000
+
+ 2 segm 4.0000 0.1000 2.9000 4.3000 47.0000 2.2000 2.2000
+
+ 3 segm -2.1000 -2.1000 3.1000 3.0000 -15.0000 2.2000 2.2000
+
+ 4 elip 0.0000 2.5000 1.0000 1.0000 0.0000 1.5000 1.5000
+
+ 5 elip 0.1000 4.3000 1.3000 1.2000 0.0000 1.0000 1.0000
+
+ 6 elip -2.0000 3.5000 1.3000 1.3000 0.0000 1.0000 1.0000
+
+ 7 elip -1.0000 1.5000 1.3000 1.1000 0.0000 1.0000 1.0000
+
+ 8 elip 1.5000 1.6000 1.5000 1.3000 0.0000 1.0000 1.0000
+
+ 9 elip 1.8000 3.3000 1.2000 1.2000 0.0000 1.0000 1.0000
+
+ scale factor multiplying object densities 1.0000
+
+ seed set to 0
+ inhomogeneity set to 0.0000
+
+ <#> phan aver 3
+
+ density in each pixel is obtained as the average of 3 x 3 points
+
+
+ <#> pixe 31 size 0.4000
+ picture size 31 x 31, pixel size 0.4000
+
+ test picture read
+ EXAMPLE 1 DAISY
+ 0.000 seconds used for processing command pict
+
+
+ <*>
+
+ <#> PROJECTION REAL
+
+ EXAMPLE 1 DAISY
+
+
+ <#> spec mono 75
+ energy spectrum is monochromatic at energy level 75
+
+
+ <#> obje
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 75
+
+ 1 tria 1.9000 -3.7000 0.5000 4.1000 20.0000 2.5000 2.5000
+
+ 2 segm 4.0000 0.1000 2.9000 4.3000 47.0000 2.2000 2.2000
+
+ 3 segm -2.1000 -2.1000 3.1000 3.0000 -15.0000 2.2000 2.2000
+
+ 4 elip 0.0000 2.5000 1.0000 1.0000 0.0000 1.5000 1.5000
+
+ 5 elip 0.1000 4.3000 1.3000 1.2000 0.0000 1.0000 1.0000
+
+ 6 elip -2.0000 3.5000 1.3000 1.3000 0.0000 1.0000 1.0000
+
+ 7 elip -1.0000 1.5000 1.3000 1.1000 0.0000 1.0000 1.0000
+
+ 8 elip 1.5000 1.6000 1.5000 1.3000 0.0000 1.0000 1.0000
+
+ 9 elip 1.8000 3.3000 1.2000 1.2000 0.0000 1.0000 1.0000
+
+ scale factor multiplying object densities 1.0000
+
+ seed set to 0
+ inhomogeneity set to 0.0000
+
+ <#> 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 20.0000 det dist10000.0000
+ rays are divergent from point sources
+ source to origin distance 20.0000
+ the detectors lie on an arc with source to detector distance = 10000.0000
+
+
+ <#> rays user 65 spacing 150.0000
+ number of rays per projection 65
+ snark computed number of rays 63
+ at detector spacing 150.0000
+
+
+ <#> angl 60
+ total number of projections 60
+
+
+ projection angles 0.0000 5.9746 11.9492 17.9237 23.8983 29.8729 35.8475 41.8220 47.7966 53.7712
+ 59.7458 65.7203 71.6949 77.6695 83.6441 89.6186 95.5932 101.5678 107.5424 113.5169
+ 119.4915 125.4661 131.4407 137.4153 143.3898 149.3644 155.3390 161.3136 167.2881 173.2627
+ 179.2373 185.2119 191.1864 197.1610 203.1356 209.1102 215.0847 221.0593 227.0339 233.0085
+ 238.9831 244.9576 250.9322 256.9068 262.8814 268.8559 274.8305 280.8051 286.7797 292.7542
+ 298.7288 304.7034 310.6780 316.6525 322.6271 328.6017 334.5763 340.5508 346.5254 352.5000
+
+
+ <#> meas perf
+ projection data are noiseless
+
+
+ <#> back 0.0000
+ at levels
+ 75
+ background absorption 0.0000
+
+ estimate of totlen = 31273.639204
+ estimate of totden = 11129.885164
+ estimate of average density = 0.3559
+ projection data read
+ EXAMPLE 1 DAISY
+ 0.000 seconds used for processing command proj
+
+
+ <*>
+
+ <#> EXECUTE DCON
+
+ EXAMPLE 1 DIVERGENT CONVOLUTION FROM 60 VIEWS
+
+ <#> -1 1 -1 0.5 0
+ cubic spline interpolation
+ one out of every 1 projection(s) used for reconstruction
+ 64 points used on each side of convolution
+ pre-smoothing weight on center ray = 0.5000
+ pre-smoothing weight on side rays = 0.2500
+ accurate back-projection
+
+ <#> BANDLIMITING
+ band-limiting filter
+ total time for obtaining smoothed projection data was 0.001
+ total time for convoluting projection data was 0.000
+ total time for back-projecting was 0.001
+ algorithm executed in iteration 1
+ 0.002 seconds for the execution of the algorithm
+ reconstruction completed after iteration 1
+ 0.002 seconds for this iteration
+ 0.002 seconds used for processing command exec
+
+
+ <*>
+
+ <#> LINES SCALE 100.0 COLUMN 16
+
+ last iteration
+
+ lines display of reconstruction using DCON iter 1
+
+ phantom name: EXAMPLE 1 DAISY
+ projection data: EXAMPLE 1 DAISY
+ execution name: EXAMPLE 1 DIVERGENT CONVOLUTION FROM 60 VIEWS
+ scaling factor = 100.00000
+
+ column = 16
+ row original estimate difference
+ 0 0.00000 -7.80005 -7.80005
+ 1 22.22222 28.11589 5.89367
+ 2 100.00000 88.37252 -11.62748
+ 3 100.00000 104.43894 4.43894
+ 4 100.00000 102.69494 2.69494
+ 5 100.00000 107.73658 7.73658
+ 6 100.00000 147.49863 47.49863
+ 7 216.66667 188.04080 -28.62586
+ 8 150.00000 180.51891 30.51891
+ 9 216.66667 211.83020 -4.83647
+ 10 250.00000 245.52265 -4.47735
+ 11 216.66667 199.38366 -17.28301
+ 12 111.11111 120.98082 9.86971
+ 13 66.66667 68.05861 1.39194
+ 14 0.00000 13.57312 13.57312
+ 15 27.77778 12.48964 -15.28813
+ 16 0.00000 27.99191 27.99191
+ 17 0.00000 24.69207 24.69207
+ 18 0.00000 10.47314 10.47314
+ 19 0.00000 -0.41472 -0.41472
+ 20 0.00000 0.60572 0.60572
+ 21 0.00000 10.05671 10.05671
+ 22 122.22222 102.66325 -19.55897
+ 23 97.77778 87.82826 -9.94952
+ 24 0.00000 5.13094 5.13094
+ 25 0.00000 -2.85855 -2.85855
+ 26 0.00000 -3.34366 -3.34366
+ 27 0.00000 1.38876 1.38876
+ 28 0.00000 3.44155 3.44155
+ 29 0.00000 5.01602 5.01602
+ 30 0.00000 -0.04334 -0.04334
+
+
+
+ 0.001 seconds used for processing command line
+
+
+ <*>
+
+ <#> PUNCH PHANTOM
+
+ last iteration
+
+ reconstruction of EXAMPLE 1 DIVERGENT CONVOLUTIO using DCON iter 1, EXAMPLE 1 DAISY
+ 0.000 seconds used for processing command punc
+
+
+ <*>
+
+ <#> END