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