--- /dev/null
+ snark14UserDefined.s170426 - A PICTURE RECONSTRUCTION PROGRAM
+
+
+ <*> *** EXAMPLE 2
+
+ <*> Creation of a small test pattern and projection data with execution of
+
+ <*> a user-defined reconstruction algorithm and analysis.
+
+ <*>
+
+ <#> CREATE
+
+ Example 2, A simple test pattern and projection data
+
+
+ <#> SPECTRUM MONOCHROMATIC with energy at 75 KeV
+ energy spectrum is monochromatic at energy level 75
+
+
+ <#> OBJECTS cx cy u v ang den
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 75
+
+ 1 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000
+
+ 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500
+
+ 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000
+
+ scale factor multiplying object densities 1.0000
+
+ seed set to 0
+ inhomogeneity set to 0.0000
+
+ <#> PHANTOM AVERAGE 1 point for estimating pixel density
+
+ this run will generate a phantom
+ density in each pixel is obtained as the average of 1 x 1 points
+
+
+ <#> number of elements is 31 Size of a pixel is 0.07
+ picture size 31 x 31, pixel size 0.0700
+
+
+ <#> RAYSUM AVERAGE 1 line for estimation of ray sum with weight
+
+ this run will generate projection data
+ projection data are calculated by dividing each ray interval into 1 substrips
+
+ with aperture (substrip) weights 1
+
+
+ <#> GEOMETRY
+
+
+ <#> PARALLEL UNIFORM LINE
+ rays are parallel with uniform spacing between rays
+ data collected along lines
+
+
+
+ <#> RAYS USER-defined 63 rays per projection at spacing 0.05
+ number of rays per projection 63
+ at detector spacing 0.0500
+
+
+ <#> ANGLES 30 EQUAL SPACING between the first and last angles which are
+ total number of projections 30
+
+ projection angles 0.0000 6.0000 12.0000 18.0000 24.0000 30.0000 36.0000 42.0000 48.0000 54.0000
+ 60.0000 66.0000 72.0000 78.0000 84.0000 90.0000 96.0000 102.0000 108.0000 114.0000
+ 120.0000 126.0000 132.0000 138.0000 144.0000 150.0000 156.0000 162.0000 168.0000 174.0000
+
+
+ <#> MEASUREMENT PERFECT
+ projection data are noiseless
+
+
+ <#> BACKGROUND absorption is 0.0
+ at levels
+ 75
+ background absorption 0.0000
+
+
+ <#> RUN
+ 0.004 seconds phantom creation
+ 0.000 seconds projection data creation
+ 0.004 seconds used for processing command crea
+
+
+ <*>
+
+ <#> PICTURE TEST
+
+ Example 2, A simple test pattern and projection data
+
+
+ <#> 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 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000
+
+ 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500
+
+ 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000
+
+ scale factor multiplying object densities 1.0000
+
+ seed set to 0
+ inhomogeneity set to 0.0000
+
+ <#> phan aver 1
+
+ density in each pixel is obtained as the average of 1 x 1 points
+
+
+ <#> pixe 31 size 0.0700
+ picture size 31 x 31, pixel size 0.0700
+
+ test picture read
+ Example 2, A simple test pattern and projection data
+ 0.000 seconds used for processing command pict
+
+
+ <*>
+
+ <#> PROJECTION REAL
+
+ Example 2, A simple test pattern and projection data
+
+
+ <#> 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 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000
+
+ 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500
+
+ 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000
+
+ 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
+
+
+ <#> para unif line
+ rays are parallel with uniform spacing between rays
+ data collected along lines
+
+
+
+ <#> rays user 63 spacing 0.0500
+ number of rays per projection 63
+ snark computed number of rays 63
+ at detector spacing 0.0500
+
+
+ <#> angl 30
+ total number of projections 30
+
+
+ projection angles 0.0000 6.0000 12.0000 18.0000 24.0000 30.0000 36.0000 42.0000 48.0000 54.0000
+ 60.0000 66.0000 72.0000 78.0000 84.0000 90.0000 96.0000 102.0000 108.0000 114.0000
+ 120.0000 126.0000 132.0000 138.0000 144.0000 150.0000 156.0000 162.0000 168.0000 174.0000
+
+
+ <#> meas perf
+ projection data are noiseless
+
+
+ <#> back 0.0000
+ at levels
+ 75
+ background absorption 0.0000
+
+ estimate of totlen = 2823.757863
+ estimate of totden = 2677.896607
+ estimate of average density = 0.9483
+ projection data read
+ Example 2, A simple test pattern and projection data
+ 0.001 seconds used for processing command proj
+
+
+ <*>
+
+ <#> MODE LOWER constraint set at 0.0 UPPER constraint set at 2.5
+ lower constraint set to 0.0000
+ upper constraint set to 2.5000
+ 0.000 seconds used for processing command mode
+
+
+ <*>
+
+ <#> SELECT USER RAYSEQ
+
+ <#> STEP projections by 10 and rays by 1
+ sequential ray selection with ray-step 1 and projection-step 10
+ 0.000 seconds used for processing command sele
+
+
+ <*>
+
+ <#> STOP TERMINATION VARIANCE at 0.01
+ termination test vari
+ epsilon = 0.010000
+ 0.000 seconds used for processing command stop
+
+
+ <*>
+
+ <#> EXECUTE AVERAGE ART SMOOTH
+
+ Example 2a underrelaxed art4 with variable tolerance
+
+ <#> smoothing threshold is 0.2 Weights are 2.25 1.5 1.0
+ reconstruction is smoothed after
+
+ iterations 1 2 3 4 5 6 7 8 9 10
+ 11 12 13 14 15 16 17 18 19 20
+ 21 22 23 24 25 26 27 28 29 30
+ 31 32 33 34 35 36 37 38 39 40
+ 41 42 43 44 45 46 47 48 49 50
+ last iteration
+
+ threshold = 0.20000 w1 = 2.25000 w2 = 1.50000 w3 = 1.00000
+
+ <#> ART4 RELAXATION CONSTANT = 0.025 TOLERANCE VARIABLE decreasing
+ art4 method
+ relaxation parameter is 0.0250
+ tolerance is variable
+
+ <#> CONSTRAINT BOUND CONRELAX CONSTANT = 0.5 NOMLIZE at each iteration
+ constraint set to boundaries
+ relax constraints with cr= 0.5000
+ picture is normalized to 0.9483 after each iteration
+ 1890 rays are used for each iteration
+ lower constraint = 0.0000
+ upper constraint = 2.5000
+ algorithm executed in iteration 1
+ 0.001 seconds for the execution of the algorithm
+ iteration 1 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 2
+ 0.001 seconds for the execution of the algorithm
+ iteration 2 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 3
+ 0.002 seconds for the execution of the algorithm
+ iteration 3 completed
+ 0.002 seconds for this iteration
+ algorithm executed in iteration 4
+ 0.001 seconds for the execution of the algorithm
+ iteration 4 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 5
+ 0.003 seconds for the execution of the algorithm
+ iteration 5 completed
+ 0.003 seconds for this iteration
+ algorithm executed in iteration 6
+ 0.002 seconds for the execution of the algorithm
+ iteration 6 completed
+ 0.002 seconds for this iteration
+ algorithm executed in iteration 7
+ 0.001 seconds for the execution of the algorithm
+ iteration 7 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 8
+ 0.001 seconds for the execution of the algorithm
+ iteration 8 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 9
+ 0.001 seconds for the execution of the algorithm
+ iteration 9 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 10
+ 0.002 seconds for the execution of the algorithm
+ iteration 10 completed
+ 0.002 seconds for this iteration
+ algorithm executed in iteration 11
+ 0.001 seconds for the execution of the algorithm
+ iteration 11 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 12
+ 0.002 seconds for the execution of the algorithm
+ iteration 12 completed
+ 0.002 seconds for this iteration
+ algorithm executed in iteration 13
+ 0.001 seconds for the execution of the algorithm
+ iteration 13 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 14
+ 0.001 seconds for the execution of the algorithm
+ iteration 14 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 15
+ 0.001 seconds for the execution of the algorithm
+ iteration 15 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 16
+ 0.000 seconds for the execution of the algorithm
+ iteration 16 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 17
+ 0.002 seconds for the execution of the algorithm
+ iteration 17 completed
+ 0.002 seconds for this iteration
+ algorithm executed in iteration 18
+ 0.001 seconds for the execution of the algorithm
+ iteration 18 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 19
+ 0.001 seconds for the execution of the algorithm
+ iteration 19 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 20
+ 0.001 seconds for the execution of the algorithm
+ iteration 20 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 21
+ 0.002 seconds for the execution of the algorithm
+ iterative process stops at iteration 21
+ the change in variance is less than 0.009852 of the variance
+ reconstruction completed after iteration 21
+ 0.002 seconds for this iteration
+ 0.031 seconds for all iterations
+ 0.031 seconds used for processing command exec
+
+
+ <*>
+
+ <#> STOP ITERATIONS = 30
+ 30 iterations
+ 0.000 seconds used for processing command stop
+
+
+ <*>
+
+ <#> EXECUTE ALP1
+
+ Example 2b underrelaxed art with extrapolation
+ algorithm executed in iteration 1
+ 0.001 seconds for the execution of the algorithm
+ iteration 1 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 2
+ 0.001 seconds for the execution of the algorithm
+ iteration 2 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 5.062718796778e-01
+ algorithm executed in iteration 3
+ 0.001 seconds for the execution of the algorithm
+ iteration 3 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 5.062718796778e-01
+ algorithm executed in iteration 4
+ 0.000 seconds for the execution of the algorithm
+ iteration 4 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 6.560964044719e-01
+ algorithm executed in iteration 5
+ 0.001 seconds for the execution of the algorithm
+ iteration 5 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 6.560964044719e-01
+ algorithm executed in iteration 6
+ 0.000 seconds for the execution of the algorithm
+ iteration 6 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 7.480594958675e-01
+ algorithm executed in iteration 7
+ 0.001 seconds for the execution of the algorithm
+ iteration 7 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 7.480594958675e-01
+ algorithm executed in iteration 8
+ 0.000 seconds for the execution of the algorithm
+ iteration 8 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 7.951514183812e-01
+ algorithm executed in iteration 9
+ 0.001 seconds for the execution of the algorithm
+ iteration 9 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 7.951514183812e-01
+ algorithm executed in iteration 10
+ 0.000 seconds for the execution of the algorithm
+ iteration 10 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.250655739411e-01
+ algorithm executed in iteration 11
+ 0.002 seconds for the execution of the algorithm
+ iteration 11 completed
+ 0.002 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.250655739411e-01
+ algorithm executed in iteration 12
+ 0.000 seconds for the execution of the algorithm
+ iteration 12 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.467169748437e-01
+ algorithm executed in iteration 13
+ 0.001 seconds for the execution of the algorithm
+ iteration 13 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.467169748437e-01
+ algorithm executed in iteration 14
+ 0.000 seconds for the execution of the algorithm
+ iteration 14 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.628058734349e-01
+ algorithm executed in iteration 15
+ 0.002 seconds for the execution of the algorithm
+ iteration 15 completed
+ 0.002 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.628058734349e-01
+ algorithm executed in iteration 16
+ 0.000 seconds for the execution of the algorithm
+ iteration 16 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.749346145383e-01
+ algorithm executed in iteration 17
+ 0.001 seconds for the execution of the algorithm
+ iteration 17 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.749346145383e-01
+ algorithm executed in iteration 18
+ 0.000 seconds for the execution of the algorithm
+ iteration 18 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.843623393369e-01
+ algorithm executed in iteration 19
+ 0.001 seconds for the execution of the algorithm
+ iteration 19 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.843623393369e-01
+ algorithm executed in iteration 20
+ 0.000 seconds for the execution of the algorithm
+ iteration 20 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.920147280181e-01
+ algorithm executed in iteration 21
+ 0.001 seconds for the execution of the algorithm
+ iteration 21 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.920147280181e-01
+ algorithm executed in iteration 22
+ 0.000 seconds for the execution of the algorithm
+ iteration 22 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.985122210336e-01
+ algorithm executed in iteration 23
+ 0.001 seconds for the execution of the algorithm
+ iteration 23 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 8.985122210336e-01
+ algorithm executed in iteration 24
+ 0.000 seconds for the execution of the algorithm
+ iteration 24 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 9.042443038619e-01
+ algorithm executed in iteration 25
+ 0.001 seconds for the execution of the algorithm
+ iteration 25 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 9.042443038619e-01
+ algorithm executed in iteration 26
+ 0.000 seconds for the execution of the algorithm
+ iteration 26 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 9.094448804097e-01
+ algorithm executed in iteration 27
+ 0.001 seconds for the execution of the algorithm
+ iteration 27 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 9.094448804097e-01
+ algorithm executed in iteration 28
+ 0.000 seconds for the execution of the algorithm
+ iteration 28 completed
+ 0.000 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 9.142500229503e-01
+ algorithm executed in iteration 29
+ 0.001 seconds for the execution of the algorithm
+ iteration 29 completed
+ 0.001 seconds for this iteration
+ EXTRAPOLATION PARAMETER IN EXTRAP IS 9.142500229503e-01
+ algorithm executed in iteration 30
+ 0.000 seconds for the execution of the algorithm
+ reconstruction completed after iteration 30
+ 0.000 seconds for this iteration
+ 0.017 seconds for all iterations
+ 0.018 seconds used for processing command exec
+
+
+ <*>
+
+ <*> Reset mode to unconstrained
+
+ <#> MODE
+ 0.000 seconds used for processing command mode
+
+
+ <*>
+
+ <#> EVALUATE RESOLUTION
+
+ Example 2a underrelaxed art4
+
+ <#> WHOLEPIC
+ Region cx cy u v ang t1 t2
+
+ wholepic -1e+20 1e+20
+
+ iterations 1 2 3 4 5 6 7 8 9 10
+ 11 12 13 14 15 16 17 18 19 20
+ 21 22 23 24 25 26 27 28 29 30
+ 31 32 33 34 35 36 37 38 39
+ last iteration
+
+ residual calculation using line integration
+ 0.047 seconds used for processing command eval
+
+
+ <*>
+
+ <#> LINES COLUMNS 15 and 18 are to be printed
+
+ iterations 10 20 30
+ last iteration
+
+ lines display of reconstruction using ART iter 10
+
+ phantom name: Example 2, A simple test pattern and projection data
+ projection data: Example 2, A simple test pattern and projection data
+ execution name: Example 2a underrelaxed art4 with variable tolerance
+ scaling factor = 1.00000
+
+ column = 15 18
+ row original estimate difference original estimate difference
+ 0 0.00000 0.18210 0.18210 0.00000 0.00909 0.00909
+ 1 0.00000 0.20450 0.20450 0.00000 0.24753 0.24753
+ 2 0.00000 1.10539 1.10539 0.00000 0.52569 0.52569
+ 3 2.50000 2.13754 -0.36246 2.50000 1.94767 -0.55233
+ 4 2.50000 2.16437 -0.33563 2.50000 2.18812 -0.31188
+ 5 2.50000 2.16394 -0.33606 2.50000 2.20758 -0.29242
+ 6 2.50000 2.13317 -0.36683 2.50000 2.17610 -0.32390
+ 7 1.25000 1.49667 0.24667 1.25000 1.68287 0.43287
+ 8 1.25000 1.44121 0.19121 1.25000 1.41529 0.16529
+ 9 1.25000 1.37200 0.12200 1.25000 1.38374 0.13374
+ 10 1.25000 1.34328 0.09328 1.25000 1.34276 0.09276
+ 11 1.25000 1.33396 0.08396 1.25000 1.32222 0.07222
+ 12 1.35000 1.34036 -0.00964 1.25000 1.31411 0.06411
+ 13 1.35000 1.34676 -0.00324 1.25000 1.31521 0.06521
+ 14 1.35000 1.35058 0.00058 1.25000 1.31568 0.06568
+ 15 1.35000 1.35053 0.00053 1.35000 1.31710 -0.03290
+ 16 1.35000 1.35058 0.00058 1.25000 1.31607 0.06607
+ 17 1.35000 1.34676 -0.00324 1.25000 1.31587 0.06587
+ 18 1.35000 1.34035 -0.00965 1.25000 1.31417 0.06417
+ 19 1.25000 1.33395 0.08395 1.25000 1.32131 0.07131
+ 20 1.25000 1.34326 0.09326 1.25000 1.34021 0.09021
+ 21 1.25000 1.37196 0.12196 1.25000 1.38008 0.13008
+ 22 1.25000 1.44114 0.19114 1.25000 1.41090 0.16090
+ 23 1.25000 1.49658 0.24658 1.25000 1.68040 0.43040
+ 24 2.50000 2.13333 -0.36667 2.50000 2.17405 -0.32595
+ 25 2.50000 2.16414 -0.33586 2.50000 2.20535 -0.29465
+ 26 2.50000 2.16455 -0.33545 2.50000 2.18655 -0.31345
+ 27 2.50000 2.13774 -0.36226 2.50000 1.94456 -0.55544
+ 28 0.00000 1.10534 1.10534 0.00000 0.44393 0.44393
+ 29 0.00000 0.20439 0.20439 0.00000 0.29112 0.29112
+ 30 0.00000 0.18198 0.18198 0.00000 0.01326 0.01326
+
+
+ lines display of reconstruction using ART iter 20
+
+ phantom name: Example 2, A simple test pattern and projection data
+ projection data: Example 2, A simple test pattern and projection data
+ execution name: Example 2a underrelaxed art4 with variable tolerance
+ scaling factor = 1.00000
+
+ column = 15 18
+ row original estimate difference original estimate difference
+ 0 0.00000 0.00164 0.00164 0.00000 0.02169 0.02169
+ 1 0.00000 0.00132 0.00132 0.00000 0.04480 0.04480
+ 2 0.00000 0.94940 0.94940 0.00000 0.13226 0.13226
+ 3 2.50000 2.43151 -0.06849 2.50000 2.25420 -0.24580
+ 4 2.50000 2.43097 -0.06903 2.50000 2.44794 -0.05206
+ 5 2.50000 2.37936 -0.12064 2.50000 2.44383 -0.05617
+ 6 2.50000 2.31917 -0.18083 2.50000 2.39051 -0.10949
+ 7 1.25000 1.46751 0.21751 1.25000 1.65881 0.40881
+ 8 1.25000 1.31935 0.06935 1.25000 1.33707 0.08707
+ 9 1.25000 1.29770 0.04770 1.25000 1.31051 0.06051
+ 10 1.25000 1.27744 0.02744 1.25000 1.27762 0.02762
+ 11 1.25000 1.28043 0.03043 1.25000 1.26720 0.01720
+ 12 1.35000 1.30399 -0.04601 1.25000 1.26843 0.01843
+ 13 1.35000 1.32095 -0.02905 1.25000 1.27985 0.02985
+ 14 1.35000 1.32945 -0.02055 1.25000 1.28645 0.03645
+ 15 1.35000 1.32847 -0.02153 1.35000 1.29098 -0.05902
+ 16 1.35000 1.32944 -0.02056 1.25000 1.28714 0.03714
+ 17 1.35000 1.32093 -0.02907 1.25000 1.28136 0.03136
+ 18 1.35000 1.30399 -0.04601 1.25000 1.27000 0.02000
+ 19 1.25000 1.28044 0.03044 1.25000 1.26871 0.01871
+ 20 1.25000 1.27744 0.02744 1.25000 1.27815 0.02815
+ 21 1.25000 1.29767 0.04767 1.25000 1.30973 0.05973
+ 22 1.25000 1.31930 0.06930 1.25000 1.33504 0.08504
+ 23 1.25000 1.46737 0.21737 1.25000 1.71927 0.46927
+ 24 2.50000 2.31941 -0.18059 2.50000 2.39181 -0.10819
+ 25 2.50000 2.37969 -0.12031 2.50000 2.44490 -0.05510
+ 26 2.50000 2.43125 -0.06875 2.50000 2.44917 -0.05083
+ 27 2.50000 2.43179 -0.06821 2.50000 2.25444 -0.24556
+ 28 0.00000 0.94928 0.94928 0.00000 0.12809 0.12809
+ 29 0.00000 0.00158 0.00158 0.00000 0.04576 0.04576
+ 30 0.00000 0.00199 0.00199 0.00000 0.02628 0.02628
+
+
+ lines display of reconstruction using ART iter 21
+
+ phantom name: Example 2, A simple test pattern and projection data
+ projection data: Example 2, A simple test pattern and projection data
+ execution name: Example 2a underrelaxed art4 with variable tolerance
+ scaling factor = 1.00000
+
+ column = 15 18
+ row original estimate difference original estimate difference
+ 0 0.00000 0.00100 0.00100 0.00000 0.01742 0.01742
+ 1 0.00000 0.00070 0.00070 0.00000 0.03672 0.03672
+ 2 0.00000 0.93506 0.93506 0.00000 0.11759 0.11759
+ 3 2.50000 2.44949 -0.05051 2.50000 2.27648 -0.22352
+ 4 2.50000 2.44615 -0.05385 2.50000 2.45838 -0.04162
+ 5 2.50000 2.39029 -0.10971 2.50000 2.45205 -0.04795
+ 6 2.50000 2.32845 -0.17155 2.50000 2.39929 -0.10071
+ 7 1.25000 1.46233 0.21233 1.25000 1.65516 0.40516
+ 8 1.25000 1.31367 0.06367 1.25000 1.33085 0.08085
+ 9 1.25000 1.29279 0.04279 1.25000 1.30546 0.05546
+ 10 1.25000 1.27400 0.02400 1.25000 1.27412 0.02412
+ 11 1.25000 1.27844 0.02844 1.25000 1.26496 0.01496
+ 12 1.35000 1.30374 -0.04626 1.25000 1.26714 0.01714
+ 13 1.35000 1.32159 -0.02841 1.25000 1.27948 0.02948
+ 14 1.35000 1.33036 -0.01964 1.25000 1.28660 0.03660
+ 15 1.35000 1.32919 -0.02081 1.35000 1.29140 -0.05860
+ 16 1.35000 1.33034 -0.01966 1.25000 1.28726 0.03726
+ 17 1.35000 1.32157 -0.02843 1.25000 1.28096 0.03096
+ 18 1.35000 1.30372 -0.04628 1.25000 1.26874 0.01874
+ 19 1.25000 1.27844 0.02844 1.25000 1.26657 0.01657
+ 20 1.25000 1.27400 0.02400 1.25000 1.27479 0.02479
+ 21 1.25000 1.29277 0.04277 1.25000 1.30483 0.05483
+ 22 1.25000 1.31362 0.06362 1.25000 1.32898 0.07898
+ 23 1.25000 1.46218 0.21218 1.25000 1.65781 0.40781
+ 24 2.50000 2.32869 -0.17131 2.50000 2.40065 -0.09935
+ 25 2.50000 2.39063 -0.10937 2.50000 2.45318 -0.04682
+ 26 2.50000 2.44644 -0.05356 2.50000 2.45967 -0.04033
+ 27 2.50000 2.44977 -0.05023 2.50000 2.27692 -0.22308
+ 28 0.00000 0.93489 0.93489 0.00000 0.11431 0.11431
+ 29 0.00000 0.00120 0.00120 0.00000 0.03840 0.03840
+ 30 0.00000 0.00161 0.00161 0.00000 0.02249 0.02249
+
+
+ lines display of reconstruction using ALP1 iter 10
+
+ phantom name: Example 2, A simple test pattern and projection data
+ projection data: Example 2, A simple test pattern and projection data
+ execution name: Example 2b underrelaxed art with extrapolation
+ scaling factor = 1.00000
+
+ column = 15 18
+ row original estimate difference original estimate difference
+ 0 0.00000 0.12670 0.12670 0.00000 -0.02567 -0.02567
+ 1 0.00000 0.17187 0.17187 0.00000 0.16958 0.16958
+ 2 0.00000 1.20151 1.20151 0.00000 0.48350 0.48350
+ 3 2.50000 2.23551 -0.26449 2.50000 2.01480 -0.48520
+ 4 2.50000 2.39456 -0.10544 2.50000 2.31508 -0.18492
+ 5 2.50000 2.34966 -0.15034 2.50000 2.46624 -0.03376
+ 6 2.50000 2.31124 -0.18876 2.50000 2.29509 -0.20491
+ 7 1.25000 1.63106 0.38106 1.25000 1.81482 0.56482
+ 8 1.25000 1.35183 0.10183 1.25000 1.44203 0.19203
+ 9 1.25000 1.42885 0.17885 1.25000 1.45015 0.20015
+ 10 1.25000 1.31859 0.06859 1.25000 1.36973 0.11973
+ 11 1.25000 1.42731 0.17731 1.25000 1.34954 0.09954
+ 12 1.35000 1.31939 -0.03061 1.25000 1.31390 0.06390
+ 13 1.35000 1.37323 0.02323 1.25000 1.32415 0.07415
+ 14 1.35000 1.46590 0.11590 1.25000 1.37584 0.12584
+ 15 1.35000 1.13749 -0.21251 1.35000 1.25214 -0.09786
+ 16 1.35000 1.46594 0.11594 1.25000 1.36902 0.11902
+ 17 1.35000 1.37317 0.02317 1.25000 1.30699 0.05699
+ 18 1.35000 1.31944 -0.03056 1.25000 1.34043 0.09043
+ 19 1.25000 1.42727 0.17727 1.25000 1.33982 0.08982
+ 20 1.25000 1.31861 0.06861 1.25000 1.35654 0.10654
+ 21 1.25000 1.42880 0.17880 1.25000 1.44218 0.19218
+ 22 1.25000 1.35172 0.10172 1.25000 1.45375 0.20375
+ 23 1.25000 1.63100 0.38100 1.25000 1.80751 0.55751
+ 24 2.50000 2.31161 -0.18839 2.50000 2.27967 -0.22033
+ 25 2.50000 2.34937 -0.15063 2.50000 2.50118 0.00118
+ 26 2.50000 2.39545 -0.10455 2.50000 2.30942 -0.19058
+ 27 2.50000 2.23526 -0.26474 2.50000 1.99542 -0.50458
+ 28 0.00000 1.20133 1.20133 0.00000 0.49478 0.49478
+ 29 0.00000 0.17173 0.17173 0.00000 0.15734 0.15734
+ 30 0.00000 0.12647 0.12647 0.00000 -0.02832 -0.02832
+
+
+ lines display of reconstruction using ALP1 iter 20
+
+ phantom name: Example 2, A simple test pattern and projection data
+ projection data: Example 2, A simple test pattern and projection data
+ execution name: Example 2b underrelaxed art with extrapolation
+ scaling factor = 1.00000
+
+ column = 15 18
+ row original estimate difference original estimate difference
+ 0 0.00000 -0.09663 -0.09663 0.00000 -0.18986 -0.18986
+ 1 0.00000 -0.12331 -0.12331 0.00000 -0.05402 -0.05402
+ 2 0.00000 1.04579 1.04579 0.00000 0.12373 0.12373
+ 3 2.50000 2.51540 0.01540 2.50000 2.25540 -0.24460
+ 4 2.50000 2.59753 0.09753 2.50000 2.56869 0.06869
+ 5 2.50000 2.55529 0.05529 2.50000 2.67846 0.17846
+ 6 2.50000 2.40841 -0.09159 2.50000 2.45046 -0.04954
+ 7 1.25000 1.52857 0.27857 1.25000 1.73851 0.48851
+ 8 1.25000 1.19882 -0.05118 1.25000 1.29003 0.04003
+ 9 1.25000 1.28220 0.03220 1.25000 1.30029 0.05029
+ 10 1.25000 1.18525 -0.06475 1.25000 1.24281 -0.00719
+ 11 1.25000 1.27819 0.02819 1.25000 1.24380 -0.00620
+ 12 1.35000 1.25011 -0.09989 1.25000 1.22237 -0.02763
+ 13 1.35000 1.30966 -0.04034 1.25000 1.24476 -0.00524
+ 14 1.35000 1.38630 0.03630 1.25000 1.28583 0.03583
+ 15 1.35000 1.11801 -0.23199 1.35000 1.20027 -0.14973
+ 16 1.35000 1.38636 0.03636 1.25000 1.28050 0.03050
+ 17 1.35000 1.30959 -0.04041 1.25000 1.23298 -0.01702
+ 18 1.35000 1.25014 -0.09986 1.25000 1.24848 -0.00152
+ 19 1.25000 1.27817 0.02817 1.25000 1.23597 -0.01403
+ 20 1.25000 1.18526 -0.06474 1.25000 1.23189 -0.01811
+ 21 1.25000 1.28211 0.03211 1.25000 1.29998 0.04998
+ 22 1.25000 1.19865 -0.05135 1.25000 1.30195 0.05195
+ 23 1.25000 1.52839 0.27839 1.25000 1.73195 0.48195
+ 24 2.50000 2.40899 -0.09101 2.50000 2.44156 -0.05844
+ 25 2.50000 2.55477 0.05477 2.50000 2.70939 0.20939
+ 26 2.50000 2.59890 0.09890 2.50000 2.56774 0.06774
+ 27 2.50000 2.51498 0.01498 2.50000 2.24565 -0.25435
+ 28 0.00000 1.04557 1.04557 0.00000 0.14374 0.14374
+ 29 0.00000 -0.12338 -0.12338 0.00000 -0.05331 -0.05331
+ 30 0.00000 -0.09670 -0.09670 0.00000 -0.17428 -0.17428
+
+
+ lines display of reconstruction using ALP1 iter 30
+
+ phantom name: Example 2, A simple test pattern and projection data
+ projection data: Example 2, A simple test pattern and projection data
+ execution name: Example 2b underrelaxed art with extrapolation
+ scaling factor = 1.00000
+
+ column = 15 18
+ row original estimate difference original estimate difference
+ 0 0.00000 -0.10295 -0.10295 0.00000 -0.15511 -0.15511
+ 1 0.00000 -0.21068 -0.21068 0.00000 -0.09236 -0.09236
+ 2 0.00000 0.92855 0.92855 0.00000 -0.06300 -0.06300
+ 3 2.50000 2.64511 0.14511 2.50000 2.37419 -0.12581
+ 4 2.50000 2.61843 0.11843 2.50000 2.63291 0.13291
+ 5 2.50000 2.61277 0.11277 2.50000 2.68002 0.18002
+ 6 2.50000 2.39593 -0.10407 2.50000 2.48350 -0.01650
+ 7 1.25000 1.47574 0.22574 1.25000 1.66721 0.41721
+ 8 1.25000 1.11475 -0.13525 1.25000 1.19787 -0.05213
+ 9 1.25000 1.24880 -0.00120 1.25000 1.25371 0.00371
+ 10 1.25000 1.15550 -0.09450 1.25000 1.22544 -0.02456
+ 11 1.25000 1.25545 0.00545 1.25000 1.24429 -0.00571
+ 12 1.35000 1.26430 -0.08570 1.25000 1.22826 -0.02174
+ 13 1.35000 1.32865 -0.02135 1.25000 1.25347 0.00347
+ 14 1.35000 1.40845 0.05845 1.25000 1.28571 0.03571
+ 15 1.35000 1.14519 -0.20481 1.35000 1.22750 -0.12250
+ 16 1.35000 1.40851 0.05851 1.25000 1.28178 0.03178
+ 17 1.35000 1.32857 -0.02143 1.25000 1.24374 -0.00626
+ 18 1.35000 1.26433 -0.08567 1.25000 1.25314 0.00314
+ 19 1.25000 1.25545 0.00545 1.25000 1.23609 -0.01391
+ 20 1.25000 1.15550 -0.09450 1.25000 1.21469 -0.03531
+ 21 1.25000 1.24869 -0.00131 1.25000 1.25626 0.00626
+ 22 1.25000 1.11456 -0.13544 1.25000 1.20898 -0.04102
+ 23 1.25000 1.47547 0.22547 1.25000 1.65907 0.40907
+ 24 2.50000 2.39664 -0.10336 2.50000 2.47619 -0.02381
+ 25 2.50000 2.61204 0.11204 2.50000 2.70609 0.20609
+ 26 2.50000 2.62009 0.12009 2.50000 2.63048 0.13048
+ 27 2.50000 2.64455 0.14455 2.50000 2.36445 -0.13555
+ 28 0.00000 0.92830 0.92830 0.00000 -0.04433 -0.04433
+ 29 0.00000 -0.21071 -0.21071 0.00000 -0.09436 -0.09436
+ 30 0.00000 -0.10286 -0.10286 0.00000 -0.14676 -0.14676
+
+
+
+ 0.002 seconds used for processing command line
+
+
+ <*>
+
+ <#> END
projects / snark14.git / blobdiff