********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
OBJECTS
   1 elip   -7.0   46.0   3.0   6.0   17.0  0.95  1.0
   2 elip    7.0   46.0   3.0   6.0  -17.0  1.0   0.95
   3 rect  -12.0   64.0   7.5   4.5    5.0  1.0   1.0
   4 rect   12.0   64.0   7.5   4.5   -5.0  0.95  0.95
   5 rect  -38.0   51.0   3.5  13.0  -39.0  0.95  1.0
   6 rect   38.0   51.0   3.5  13.0   39.0  1.0   0.95
   7 rect  -46.0   24.0   6.5   6.0  -18.0  0.95  1.0
   8 rect   46.0   24.0   6.5   6.0   18.0  1.0   0.95
   9 rect  -49.0    6.0   2.5  10.0   63.0  1.0   1.0
  10 rect   49.0    6.0   2.5  10.0  -63.0  0.95  0.95
  11 rect  -52.0  -14.0   9.0   7.0  -14.0  0.95  1.0
  12 rect   52.0  -14.0   9.0   7.0   14.0  1.0   0.95
  13 rect  -10.0  -56.0   5.5  10.0   -1.0  0.95  1.0
  14 rect   10.0  -56.0   5.5  10.0    1.0  1.0   0.95
  15 elip  -40.0  -47.0   9.0  22.5   48.0  1.0   1.0
  16 elip   40.0  -47.0   9.0  22.5  -48.0  0.95  0.95
  17 elip   -8.0  -22.0   3.5  15.5   -9.0  1.0   1.0
  18 elip    8.0  -22.0   3.5  15.5    9.0  0.95  0.95
  19 elip  -27.0   -6.0   5.5  23.5   -5.0  0.95  1.0
  20 elip   27.0   -6.0   5.5  23.5    5.0  1.0   0.95
  21 elip  -25.0   38.0   6.5  10.5  -14.0  1.0   1.0
  22 elip   25.0   38.0   6.5  10.5   14.0  0.95  0.95
  23 rect   -8.0   32.0   1.5   6.5   38.0  1.0   1.0
  24 rect    8.0   32.0   1.5   6.5  -38.0  0.95  0.95
  25 rect   -8.0    3.0   1.0   9.0  -33.0  0.95  1.0
  26 rect    8.0    3.0   1.0   9.0   33.0  1.0   0.95
  27 elip    0.0    0.0  66.5  74.0    0.0  1.0   1.0
LAST .51    1   0.05  
PHANTOM AVERAGE 7
95 PIXELS OF SIZE 1.6
RAYSUM AVERAGE 1
1
GEOMETRY
divergent arc 153 306
RAYS USER 101 DETECTOR SPACING 3.2
ANGLES 300 EQUAL SPACING
0.0 358.8
MEASUREMENT PERFECT
BACKGROUND 0.0
RUN
*
PICTURE TEST
*
PROJECTION REAL
*
STOP ITERATION 20
*
EXECUTE AVERAGE EMAP
Example 7 Illustrating the MAP EM algorithm for PET
gamma is 10.0 EVAL
*
END