snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM <*> *******EXAMPLE 10 <*> ROTATED RECTANGLE WITH PSEUDO PROJECTION DATA <*> RECONSTRUCTION BY ART WITH BLOB BASIS <*> <*> TRACE 1 <#> CREATE EXAMPLE 10 ROTATED RECTANGLE <#> SPECTRUM MONOCHROMATIC 10 energy spectrum is monochromatic at energy level 10 <#> OBJECTS description of objects density at levels numb type x-coord y-coord x-length y-length angle av dens 10 1 rect 0.0000 0.0000 2.0000 10.0000 -30.0000 1.0000 1.0000 scale factor multiplying object densities 1.0000 seed set to 0 inhomogeneity set to 0.0000 <#> PHANTOM AVERAGE 1 this run will generate a phantom density in each pixel is obtained as the average of 1 x 1 points <#> 25 PIXELS OF SIZE 1.0 picture size 25 x 25, pixel size 1.0000 <#> RAYSUM 0.000 seconds phantom creation 0.000 seconds used for processing command crea <*> RUN <*> <#> PICTURE TEST EXAMPLE 10 ROTATED RECTANGLE <#> spec mono 10 energy spectrum is monochromatic at energy level 10 <#> obje description of objects density at levels numb type x-coord y-coord x-length y-length angle av dens 10 1 rect 0.0000 0.0000 2.0000 10.0000 -30.0000 1.0000 1.0000 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 25 size 1.0000 picture size 25 x 25, pixel size 1.0000 test picture read EXAMPLE 10 ROTATED RECTANGLE 0.000 seconds used for processing command pict <*> <#> PROJECTION PSEUDO EXAMPLE 10 PSEUDO PROJECTION DATA FOR TWELVE VIEWS <#> GEOMETRY <#> PARALLEL UNIFORM LINE rays are parallel with uniform spacing between rays data collected along lines <#> RAYS PROGRAM 25 2.0 DETECTOR SPACING AT 2.0 number of rays per projection 37 snark computed number of rays 19 at detector spacing 2.0000 <#> ANGLES 12 EQUALLY SPACED total number of projections 12 projection angles 0.0000 15.0000 30.0000 45.0000 60.0000 75.0000 90.0000 105.0000 120.0000 135.0000 150.0000 165.0000 <#> MEASUREMENT PERFECT projection data are noiseless <#> BACKGROUND 50.0 at levels 10 background absorption 50.0000 estimate of totlen = 3776.330760 estimate of totden = 498.725717 estimate of average density = 0.1321 projection data read EXAMPLE 10 PSEUDO PROJECTION DATA FOR TWELVE VIEWS 0.000 seconds used for processing command proj <*> <#> BASIS BLOB hexagonal grid: there are 25 rows, the middle row has 21 grid points and the row next to the middle one has 22 grid points blob basis with parameters: support = 2.0629 shape = 11.2829 delta = 1.1547 0.083 seconds used for processing command basi <*> <#> SELECT SNARK RAYSEQ <#> STEP 3 1 sequential ray selection with ray-step 1 and projection-step 3 0.000 seconds used for processing command sele <*> <#> STOP ITERATION 10 10 iterations 0.000 seconds used for processing command stop <*> <#> EXECUTE AVERAGE ART ART WITH BLOB BASIS <#> ART3 art3 method relaxation parameter is 1.0 norm is 2 tolerance is 0.0 <#> CONSTRAINT ART2 constraint art2 way relaxation constraints with cr = 1.0 picture is not normalized 228 rays are used for each iteration 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.000 seconds for the execution of the algorithm iteration 2 completed 0.001 seconds for this iteration algorithm executed in iteration 3 0.000 seconds for the execution of the algorithm iteration 3 completed 0.001 seconds for this iteration algorithm executed in iteration 4 0.000 seconds for the execution of the algorithm iteration 4 completed 0.001 seconds for this iteration algorithm executed in iteration 5 0.000 seconds for the execution of the algorithm iteration 5 completed 0.001 seconds for this iteration algorithm executed in iteration 6 0.000 seconds for the execution of the algorithm iteration 6 completed 0.001 seconds for this iteration algorithm executed in iteration 7 0.000 seconds for the execution of the algorithm iteration 7 completed 0.001 seconds for this iteration algorithm executed in iteration 8 0.000 seconds for the execution of the algorithm iteration 8 completed 0.001 seconds for this iteration algorithm executed in iteration 9 0.000 seconds for the execution of the algorithm iteration 9 completed 0.001 seconds for this iteration algorithm executed in iteration 10 0.000 seconds for the execution of the algorithm reconstruction completed after iteration 10 0.001 seconds for this iteration 0.007 seconds for all iterations 0.008 seconds used for processing command exec <*> <#> BASIS PIXEL pixel basis 0.000 seconds used for processing command basi <#> EXECUTE AVERAGE ART ART WITH PIXEL BASIS <#> ART3 art3 method relaxation parameter is 1.0 norm is 2 tolerance is 0.0 <#> CONSTRAINT ART2 constraint art2 way relaxation constraints with cr = 1.0 picture is not normalized 228 rays are used for each iteration algorithm executed in iteration 1 0.000 seconds for the execution of the algorithm iteration 1 completed 0.000 seconds for this iteration algorithm executed in iteration 2 0.000 seconds for the execution of the algorithm iteration 2 completed 0.000 seconds for this iteration algorithm executed in iteration 3 0.000 seconds for the execution of the algorithm iteration 3 completed 0.000 seconds for this iteration algorithm executed in iteration 4 0.000 seconds for the execution of the algorithm iteration 4 completed 0.000 seconds for this iteration algorithm executed in iteration 5 0.000 seconds for the execution of the algorithm iteration 5 completed 0.000 seconds for this iteration algorithm executed in iteration 6 0.000 seconds for the execution of the algorithm iteration 6 completed 0.000 seconds for this iteration algorithm executed in iteration 7 0.000 seconds for the execution of the algorithm iteration 7 completed 0.000 seconds for this iteration algorithm executed in iteration 8 0.000 seconds for the execution of the algorithm iteration 8 completed 0.000 seconds for this iteration algorithm executed in iteration 9 0.000 seconds for the execution of the algorithm iteration 9 completed 0.000 seconds for this iteration algorithm executed in iteration 10 0.000 seconds for the execution of the algorithm reconstruction completed after iteration 10 0.000 seconds for this iteration 0.001 seconds for all iterations 0.002 seconds used for processing command exec <*> <#> EVALUATE BOTH COMPARING BLOB AND PIXEL BASIS <#> WHOLEPIC Region cx cy u v ang t1 t2 wholepic -1e+20 1e+20 iterations 1 2 3 4 5 6 7 8 9 last iteration residual calculation using line integration 0.004 seconds used for processing command eval <#> END