1 snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM
8 <*> CREATION OF PHANTOM AND DIVERGENT PROJECTION DATA
10 <*> RECONSTRUCTION WITH DCON
19 <#> SPECTRUM MONOCHROMATIC 75
20 energy spectrum is monochromatic at energy level 75
24 description of objects
26 numb type x-coord y-coord x-length y-length angle av dens 75
28 1 tria 1.9000 -3.7000 0.5000 4.1000 20.0000 2.5000 2.5000
30 2 segm 4.0000 0.1000 2.9000 4.3000 47.0000 2.2000 2.2000
32 3 segm -2.1000 -2.1000 3.1000 3.0000 -15.0000 2.2000 2.2000
34 4 elip 0.0000 2.5000 1.0000 1.0000 0.0000 1.5000 1.5000
36 5 elip 0.1000 4.3000 1.3000 1.2000 0.0000 1.0000 1.0000
38 6 elip -2.0000 3.5000 1.3000 1.3000 0.0000 1.0000 1.0000
40 7 elip -1.0000 1.5000 1.3000 1.1000 0.0000 1.0000 1.0000
42 8 elip 1.5000 1.6000 1.5000 1.3000 0.0000 1.0000 1.0000
44 9 elip 1.8000 3.3000 1.2000 1.2000 0.0000 1.0000 1.0000
46 scale factor multiplying object densities 1.0000
49 inhomogeneity set to 0.0000
53 this run will generate a phantom
54 density in each pixel is obtained as the average of 3 x 3 points
58 picture size 31 x 31, pixel size 0.4000
63 this run will generate projection data
64 projection data are calculated by dividing each ray interval into 1 substrips
66 with aperture (substrip) weights 1
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
78 <#> RAYS USER 65 150.0
79 number of rays per projection 65
80 at detector spacing 150.0000
83 <#> ANGLES 60 EQUALLY SPACED
84 total number of projections 60
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
94 <#> MEASUREMENT PERFECT
95 projection data are noiseless
101 background absorption 0.0000
105 0.001 seconds phantom creation
106 0.004 seconds projection data creation
107 0.005 seconds used for processing command crea
118 energy spectrum is monochromatic at energy level 75
122 description of objects
124 numb type x-coord y-coord x-length y-length angle av dens 75
126 1 tria 1.9000 -3.7000 0.5000 4.1000 20.0000 2.5000 2.5000
128 2 segm 4.0000 0.1000 2.9000 4.3000 47.0000 2.2000 2.2000
130 3 segm -2.1000 -2.1000 3.1000 3.0000 -15.0000 2.2000 2.2000
132 4 elip 0.0000 2.5000 1.0000 1.0000 0.0000 1.5000 1.5000
134 5 elip 0.1000 4.3000 1.3000 1.2000 0.0000 1.0000 1.0000
136 6 elip -2.0000 3.5000 1.3000 1.3000 0.0000 1.0000 1.0000
138 7 elip -1.0000 1.5000 1.3000 1.1000 0.0000 1.0000 1.0000
140 8 elip 1.5000 1.6000 1.5000 1.3000 0.0000 1.0000 1.0000
142 9 elip 1.8000 3.3000 1.2000 1.2000 0.0000 1.0000 1.0000
144 scale factor multiplying object densities 1.0000
147 inhomogeneity set to 0.0000
151 density in each pixel is obtained as the average of 3 x 3 points
154 <#> pixe 31 size 0.4000
155 picture size 31 x 31, pixel size 0.4000
159 0.000 seconds used for processing command pict
170 energy spectrum is monochromatic at energy level 75
174 description of objects
176 numb type x-coord y-coord x-length y-length angle av dens 75
178 1 tria 1.9000 -3.7000 0.5000 4.1000 20.0000 2.5000 2.5000
180 2 segm 4.0000 0.1000 2.9000 4.3000 47.0000 2.2000 2.2000
182 3 segm -2.1000 -2.1000 3.1000 3.0000 -15.0000 2.2000 2.2000
184 4 elip 0.0000 2.5000 1.0000 1.0000 0.0000 1.5000 1.5000
186 5 elip 0.1000 4.3000 1.3000 1.2000 0.0000 1.0000 1.0000
188 6 elip -2.0000 3.5000 1.3000 1.3000 0.0000 1.0000 1.0000
190 7 elip -1.0000 1.5000 1.3000 1.1000 0.0000 1.0000 1.0000
192 8 elip 1.5000 1.6000 1.5000 1.3000 0.0000 1.0000 1.0000
194 9 elip 1.8000 3.3000 1.2000 1.2000 0.0000 1.0000 1.0000
196 scale factor multiplying object densities 1.0000
199 inhomogeneity set to 0.0000
203 projection data are calculated by dividing each ray interval into 1 substrips
205 with aperture (substrip) weights 1
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
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
224 total number of projections 60
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
236 projection data are noiseless
242 background absorption 0.0000
244 estimate of totlen = 31273.639204
245 estimate of totden = 11129.885164
246 estimate of average density = 0.3559
249 0.000 seconds used for processing command proj
256 EXAMPLE 1 DIVERGENT CONVOLUTION FROM 60 VIEWS
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
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
280 <#> LINES SCALE 100.0 COLUMN 16
284 lines display of reconstruction using DCON iter 1
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
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
327 0.001 seconds used for processing command line
336 reconstruction of EXAMPLE 1 DIVERGENT CONVOLUTIO using DCON iter 1, EXAMPLE 1 DAISY
337 0.000 seconds used for processing command punc