X-Git-Url: http://git.kpe.io/?p=snark14.git;a=blobdiff_plain;f=examples%2Fb7%2Fb7r.out;fp=examples%2Fb7%2Fb7r.out;h=0000000000000000000000000000000000000000;hp=fec918204308bcad2ffc74ecf5c871a9d952553d;hb=8fbca36cd1abc71c89e6d49611dbd58424881a96;hpb=5de7494cf1d3672c9419ac0a3d6f26ae2b9855af diff --git a/examples/b7/b7r.out b/examples/b7/b7r.out deleted file mode 100644 index fec9182..0000000 --- a/examples/b7/b7r.out +++ /dev/null @@ -1,534 +0,0 @@ - snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM - - - <*> *******EXAMPLE 7 - - <*> MAP EM ALGORITHM FOR EMISSION TOMOGRAPHY. RECONSTRUCTION OF BRAIN - - <*> PHANTOM. SIMULATING PET GEOMETRY WITH A RING OF 300 DETECTORS WITH - - <*> EACH DETECTOR IN COINCIDENCE WITH 101 DETECTORS OPPOSITE IT. - - <*> COEFFICIENT OF PENALTLY TERM IS SET TO 10.0 - - <*> - - <#> CREATE - - EXAMPLE 7 Brain Phantom - - - <#> SPECTRUM MONOCHROMATIC 511 - energy spectrum is monochromatic at energy level 511 - - - <#> OBJECTS - description of objects - density at levels - numb type x-coord y-coord x-length y-length angle av dens 511 - - 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.9500 0.9500 - - 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 1.0000 1.0000 - - 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 1.0000 1.0000 - - 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.9500 0.9500 - - 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.9500 0.9500 - - 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 1.0000 1.0000 - - 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.9500 0.9500 - - 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 1.0000 1.0000 - - 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 1.0000 1.0000 - - 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.9500 0.9500 - - 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.9500 0.9500 - - 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 1.0000 1.0000 - - 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.9500 0.9500 - - 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 1.0000 1.0000 - - 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 1.0000 1.0000 - - 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.9500 0.9500 - - 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 1.0000 1.0000 - - 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.9500 0.9500 - - 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.9500 0.9500 - - 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 1.0000 1.0000 - - 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 1.0000 1.0000 - - 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.9500 0.9500 - - 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 1.0000 1.0000 - - 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.9500 0.9500 - - 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.9500 0.9500 - - 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 1.0000 1.0000 - - 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 1.0000 1.0000 - - scale factor multiplying object densities 0.5100 - - seed set to 1 - inhomogeneity set to 0.0500 - - <#> PHANTOM AVERAGE 7 - - this run will generate a phantom - density in each pixel is obtained as the average of 7 x 7 points - - - <#> 95 PIXELS OF SIZE 1.6 - picture size 95 x 95, pixel size 1.6000 - - - <#> 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 153 306 - rays are divergent from point sources - source to origin distance 153.0000 - the detectors lie on an arc with source to detector distance = 306.0000 - - - <#> RAYS USER 101 DETECTOR SPACING 3.2 - number of rays per projection 101 - at detector spacing 3.2000 - - - <#> ANGLES 300 EQUAL SPACING - total number of projections 300 - - projection angles 0.0000 1.2000 2.4000 3.6000 4.8000 6.0000 7.2000 8.4000 9.6000 10.8000 - 12.0000 13.2000 14.4000 15.6000 16.8000 18.0000 19.2000 20.4000 21.6000 22.8000 - 24.0000 25.2000 26.4000 27.6000 28.8000 30.0000 31.2000 32.4000 33.6000 34.8000 - 36.0000 37.2000 38.4000 39.6000 40.8000 42.0000 43.2000 44.4000 45.6000 46.8000 - 48.0000 49.2000 50.4000 51.6000 52.8000 54.0000 55.2000 56.4000 57.6000 58.8000 - 60.0000 61.2000 62.4000 63.6000 64.8000 66.0000 67.2000 68.4000 69.6000 70.8000 - 72.0000 73.2000 74.4000 75.6000 76.8000 78.0000 79.2000 80.4000 81.6000 82.8000 - 84.0000 85.2000 86.4000 87.6000 88.8000 90.0000 91.2000 92.4000 93.6000 94.8000 - 96.0000 97.2000 98.4000 99.6000 100.8000 102.0000 103.2000 104.4000 105.6000 106.8000 - 108.0000 109.2000 110.4000 111.6000 112.8000 114.0000 115.2000 116.4000 117.6000 118.8000 - 120.0000 121.2000 122.4000 123.6000 124.8000 126.0000 127.2000 128.4000 129.6000 130.8000 - 132.0000 133.2000 134.4000 135.6000 136.8000 138.0000 139.2000 140.4000 141.6000 142.8000 - 144.0000 145.2000 146.4000 147.6000 148.8000 150.0000 151.2000 152.4000 153.6000 154.8000 - 156.0000 157.2000 158.4000 159.6000 160.8000 162.0000 163.2000 164.4000 165.6000 166.8000 - 168.0000 169.2000 170.4000 171.6000 172.8000 174.0000 175.2000 176.4000 177.6000 178.8000 - 180.0000 181.2000 182.4000 183.6000 184.8000 186.0000 187.2000 188.4000 189.6000 190.8000 - 192.0000 193.2000 194.4000 195.6000 196.8000 198.0000 199.2000 200.4000 201.6000 202.8000 - 204.0000 205.2000 206.4000 207.6000 208.8000 210.0000 211.2000 212.4000 213.6000 214.8000 - 216.0000 217.2000 218.4000 219.6000 220.8000 222.0000 223.2000 224.4000 225.6000 226.8000 - 228.0000 229.2000 230.4000 231.6000 232.8000 234.0000 235.2000 236.4000 237.6000 238.8000 - 240.0000 241.2000 242.4000 243.6000 244.8000 246.0000 247.2000 248.4000 249.6000 250.8000 - 252.0000 253.2000 254.4000 255.6000 256.8000 258.0000 259.2000 260.4000 261.6000 262.8000 - 264.0000 265.2000 266.4000 267.6000 268.8000 270.0000 271.2000 272.4000 273.6000 274.8000 - 276.0000 277.2000 278.4000 279.6000 280.8000 282.0000 283.2000 284.4000 285.6000 286.8000 - 288.0000 289.2000 290.4000 291.6000 292.8000 294.0000 295.2000 296.4000 297.6000 298.8000 - 300.0000 301.2000 302.4000 303.6000 304.8000 306.0000 307.2000 308.4000 309.6000 310.8000 - 312.0000 313.2000 314.4000 315.6000 316.8000 318.0000 319.2000 320.4000 321.6000 322.8000 - 324.0000 325.2000 326.4000 327.6000 328.8000 330.0000 331.2000 332.4000 333.6000 334.8000 - 336.0000 337.2000 338.4000 339.6000 340.8000 342.0000 343.2000 344.4000 345.6000 346.8000 - 348.0000 349.2000 350.4000 351.6000 352.8000 354.0000 355.2000 356.4000 357.6000 358.8000 - - - <#> MEASUREMENT NOISY - noise characteristics of projection data follow - nature characteristics - - <#> QUANTUM 1.0 1.0 CALIBRATION 4 - Emission tomography - - <#> SEED 0 - seed for random number generator is 0 - - - <#> BACKGROUND 0.0 - at levels - 511 - background absorption 0.0000 - - - <#> RUN - 0.025 seconds phantom creation - 0.100 seconds projection data creation - 0.125 seconds used for processing command crea - - - <*> - - <#> PICTURE TEST - - EXAMPLE 7 Brain Phantom - - - <#> spec mono 511 - energy spectrum is monochromatic at energy level 511 - - - <#> obje - description of objects - density at levels - numb type x-coord y-coord x-length y-length angle av dens 511 - - 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845 - - 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100 - - 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100 - - 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845 - - 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845 - - 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100 - - 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845 - - 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100 - - 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100 - - 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845 - - 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845 - - 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100 - - 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845 - - 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100 - - 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100 - - 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845 - - 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100 - - 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845 - - 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845 - - 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100 - - 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100 - - 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845 - - 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100 - - 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845 - - 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845 - - 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100 - - 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100 - - scale factor multiplying object densities 0.5100 - - seed set to 1 - inhomogeneity set to 0.0500 - - <#> phan aver 7 - - density in each pixel is obtained as the average of 7 x 7 points - - - <#> pixe 95 size 1.6000 - picture size 95 x 95, pixel size 1.6000 - - test picture read - EXAMPLE 7 Brain Phantom - 0.000 seconds used for processing command pict - - - <*> - - <#> PROJECTION REAL - - EXAMPLE 7 Brain Phantom - - - <#> spec mono 511 - energy spectrum is monochromatic at energy level 511 - - - <#> obje - description of objects - density at levels - numb type x-coord y-coord x-length y-length angle av dens 511 - - 1 elip -7.0000 46.0000 3.0000 6.0000 17.0000 0.4845 0.4845 - - 2 elip 7.0000 46.0000 3.0000 6.0000 -17.0000 0.5100 0.5100 - - 3 rect -12.0000 64.0000 7.5000 4.5000 5.0000 0.5100 0.5100 - - 4 rect 12.0000 64.0000 7.5000 4.5000 -5.0000 0.4845 0.4845 - - 5 rect -38.0000 51.0000 3.5000 13.0000 -39.0000 0.4845 0.4845 - - 6 rect 38.0000 51.0000 3.5000 13.0000 39.0000 0.5100 0.5100 - - 7 rect -46.0000 24.0000 6.5000 6.0000 -18.0000 0.4845 0.4845 - - 8 rect 46.0000 24.0000 6.5000 6.0000 18.0000 0.5100 0.5100 - - 9 rect -49.0000 6.0000 2.5000 10.0000 63.0000 0.5100 0.5100 - - 10 rect 49.0000 6.0000 2.5000 10.0000 -63.0000 0.4845 0.4845 - - 11 rect -52.0000 -14.0000 9.0000 7.0000 -14.0000 0.4845 0.4845 - - 12 rect 52.0000 -14.0000 9.0000 7.0000 14.0000 0.5100 0.5100 - - 13 rect -10.0000 -56.0000 5.5000 10.0000 -1.0000 0.4845 0.4845 - - 14 rect 10.0000 -56.0000 5.5000 10.0000 1.0000 0.5100 0.5100 - - 15 elip -40.0000 -47.0000 9.0000 22.5000 48.0000 0.5100 0.5100 - - 16 elip 40.0000 -47.0000 9.0000 22.5000 -48.0000 0.4845 0.4845 - - 17 elip -8.0000 -22.0000 3.5000 15.5000 -9.0000 0.5100 0.5100 - - 18 elip 8.0000 -22.0000 3.5000 15.5000 9.0000 0.4845 0.4845 - - 19 elip -27.0000 -6.0000 5.5000 23.5000 -5.0000 0.4845 0.4845 - - 20 elip 27.0000 -6.0000 5.5000 23.5000 5.0000 0.5100 0.5100 - - 21 elip -25.0000 38.0000 6.5000 10.5000 -14.0000 0.5100 0.5100 - - 22 elip 25.0000 38.0000 6.5000 10.5000 14.0000 0.4845 0.4845 - - 23 rect -8.0000 32.0000 1.5000 6.5000 38.0000 0.5100 0.5100 - - 24 rect 8.0000 32.0000 1.5000 6.5000 -38.0000 0.4845 0.4845 - - 25 rect -8.0000 3.0000 1.0000 9.0000 -33.0000 0.4845 0.4845 - - 26 rect 8.0000 3.0000 1.0000 9.0000 33.0000 0.5100 0.5100 - - 27 elip 0.0000 0.0000 66.5000 74.0000 0.0000 0.5100 0.5100 - - scale factor multiplying object densities 0.5100 - - seed set to 1 - inhomogeneity set to 0.0500 - - <#> 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 153.0000 det dist 306.0000 - rays are divergent from point sources - source to origin distance 153.0000 - the detectors lie on an arc with source to detector distance = 306.0000 - - - <#> rays user 101 spacing 3.2000 - number of rays per projection 101 - snark computed number of rays 151 - at detector spacing 3.2000 - - - <#> angl 300 - total number of projections 300 - - - projection angles 0.0000 1.2000 2.4000 3.6000 4.8000 6.0000 7.2000 8.4000 9.6000 10.8000 - 12.0000 13.2000 14.4000 15.6000 16.8000 18.0000 19.2000 20.4000 21.6000 22.8000 - 24.0000 25.2000 26.4000 27.6000 28.8000 30.0000 31.2000 32.4000 33.6000 34.8000 - 36.0000 37.2000 38.4000 39.6000 40.8000 42.0000 43.2000 44.4000 45.6000 46.8000 - 48.0000 49.2000 50.4000 51.6000 52.8000 54.0000 55.2000 56.4000 57.6000 58.8000 - 60.0000 61.2000 62.4000 63.6000 64.8000 66.0000 67.2000 68.4000 69.6000 70.8000 - 72.0000 73.2000 74.4000 75.6000 76.8000 78.0000 79.2000 80.4000 81.6000 82.8000 - 84.0000 85.2000 86.4000 87.6000 88.8000 90.0000 91.2000 92.4000 93.6000 94.8000 - 96.0000 97.2000 98.4000 99.6000 100.8000 102.0000 103.2000 104.4000 105.6000 106.8000 - 108.0000 109.2000 110.4000 111.6000 112.8000 114.0000 115.2000 116.4000 117.6000 118.8000 - 120.0000 121.2000 122.4000 123.6000 124.8000 126.0000 127.2000 128.4000 129.6000 130.8000 - 132.0000 133.2000 134.4000 135.6000 136.8000 138.0000 139.2000 140.4000 141.6000 142.8000 - 144.0000 145.2000 146.4000 147.6000 148.8000 150.0000 151.2000 152.4000 153.6000 154.8000 - 156.0000 157.2000 158.4000 159.6000 160.8000 162.0000 163.2000 164.4000 165.6000 166.8000 - 168.0000 169.2000 170.4000 171.6000 172.8000 174.0000 175.2000 176.4000 177.6000 178.8000 - 180.0000 181.2000 182.4000 183.6000 184.8000 186.0000 187.2000 188.4000 189.6000 190.8000 - 192.0000 193.2000 194.4000 195.6000 196.8000 198.0000 199.2000 200.4000 201.6000 202.8000 - 204.0000 205.2000 206.4000 207.6000 208.8000 210.0000 211.2000 212.4000 213.6000 214.8000 - 216.0000 217.2000 218.4000 219.6000 220.8000 222.0000 223.2000 224.4000 225.6000 226.8000 - 228.0000 229.2000 230.4000 231.6000 232.8000 234.0000 235.2000 236.4000 237.6000 238.8000 - 240.0000 241.2000 242.4000 243.6000 244.8000 246.0000 247.2000 248.4000 249.6000 250.8000 - 252.0000 253.2000 254.4000 255.6000 256.8000 258.0000 259.2000 260.4000 261.6000 262.8000 - 264.0000 265.2000 266.4000 267.6000 268.8000 270.0000 271.2000 272.4000 273.6000 274.8000 - 276.0000 277.2000 278.4000 279.6000 280.8000 282.0000 283.2000 284.4000 285.6000 286.8000 - 288.0000 289.2000 290.4000 291.6000 292.8000 294.0000 295.2000 296.4000 297.6000 298.8000 - 300.0000 301.2000 302.4000 303.6000 304.8000 306.0000 307.2000 308.4000 309.6000 310.8000 - 312.0000 313.2000 314.4000 315.6000 316.8000 318.0000 319.2000 320.4000 321.6000 322.8000 - 324.0000 325.2000 326.4000 327.6000 328.8000 330.0000 331.2000 332.4000 333.6000 334.8000 - 336.0000 337.2000 338.4000 339.6000 340.8000 342.0000 343.2000 344.4000 345.6000 346.8000 - 348.0000 349.2000 350.4000 351.6000 352.8000 354.0000 355.2000 356.4000 357.6000 358.8000 - - - <#> meas nois - noise characteristics of projection data follow - nature characteristics - - <#> quan 1.0000 1.0000 cali 4 - Emission tomography - - <#> seed 0 - seed for random number generator is 0 - - - <#> back 0.0000 - at levels - 511 - background absorption 0.0000 - - estimate of totlen = 4263169.769963 - estimate of totden = 2020362.000000 - estimate of average density = 0.4739 - projection data read - EXAMPLE 7 Brain Phantom - 0.012 seconds used for processing command proj - - - <*> - - <#> STOP ITERATION 20 - 20 iterations - 0.000 seconds used for processing command stop - - - <*> - - <#> EXECUTE AVERAGE EMAP - - Example 7 Illustrating the MAP EM algorithm for PET - - <#> gamma is 10.0 EVAL - - ----------------------------------------------------------- - - maximum a-posteriori probability expectation maximization - - gamma: 10.000 - evaluation flag is set - - ----------------------------------------------------------- - - algorithm executed in iteration 1 - 0.167 seconds for the execution of the algorithm - iteration 1 completed - 0.167 seconds for this iteration - algorithm executed in iteration 2 - 0.096 seconds for the execution of the algorithm - iteration 2 completed - 0.096 seconds for this iteration - algorithm executed in iteration 3 - 0.106 seconds for the execution of the algorithm - iteration 3 completed - 0.106 seconds for this iteration - algorithm executed in iteration 4 - 0.097 seconds for the execution of the algorithm - iteration 4 completed - 0.097 seconds for this iteration - algorithm executed in iteration 5 - 0.096 seconds for the execution of the algorithm - iteration 5 completed - 0.096 seconds for this iteration - algorithm executed in iteration 6 - 0.097 seconds for the execution of the algorithm - iteration 6 completed - 0.097 seconds for this iteration - algorithm executed in iteration 7 - 0.096 seconds for the execution of the algorithm - iteration 7 completed - 0.097 seconds for this iteration - algorithm executed in iteration 8 - 0.097 seconds for the execution of the algorithm - iteration 8 completed - 0.097 seconds for this iteration - algorithm executed in iteration 9 - 0.112 seconds for the execution of the algorithm - iteration 9 completed - 0.112 seconds for this iteration - algorithm executed in iteration 10 - 0.116 seconds for the execution of the algorithm - iteration 10 completed - 0.116 seconds for this iteration - algorithm executed in iteration 11 - 0.100 seconds for the execution of the algorithm - iteration 11 completed - 0.100 seconds for this iteration - algorithm executed in iteration 12 - 0.106 seconds for the execution of the algorithm - iteration 12 completed - 0.106 seconds for this iteration - algorithm executed in iteration 13 - 0.117 seconds for the execution of the algorithm - iteration 13 completed - 0.117 seconds for this iteration - algorithm executed in iteration 14 - 0.115 seconds for the execution of the algorithm - iteration 14 completed - 0.115 seconds for this iteration - algorithm executed in iteration 15 - 0.115 seconds for the execution of the algorithm - iteration 15 completed - 0.115 seconds for this iteration - algorithm executed in iteration 16 - 0.109 seconds for the execution of the algorithm - iteration 16 completed - 0.109 seconds for this iteration - algorithm executed in iteration 17 - 0.102 seconds for the execution of the algorithm - iteration 17 completed - 0.102 seconds for this iteration - algorithm executed in iteration 18 - 0.089 seconds for the execution of the algorithm - iteration 18 completed - 0.089 seconds for this iteration - algorithm executed in iteration 19 - 0.101 seconds for the execution of the algorithm - iteration 19 completed - 0.101 seconds for this iteration - algorithm executed in iteration 20 - 0.099 seconds for the execution of the algorithm - reconstruction completed after iteration 20 - 0.099 seconds for this iteration - 2.135 seconds for all iterations - 2.183 seconds used for processing command exec - - - <*> - - <#> END