PROJECTION = p ('Modified Plate Carree')
latitude and longitude of origin = 43.000 -87.000

x correction = 1.000975       (xcor)
    rotation = 0.000000 deg.  (theta)

Input File ('MTHICK.xyc') is of type from Scanned Image:
     ( Output from program 'imvis.c' )

scale =  0.66853 km/pixel   or     3.9480 km/mm   or 3948012:1
image scale = 150.000 Pixels/Inch
image root mean error =      0.000 mm
image root mean error =      0.000 pixels
earth root mean error =      0.095 m =      0.312 ft
X0 =   -110.614  Y0 =   -249.453

------------------------------------------------------------------------
         input*                output**                 errors (m) 
-- ------------------   -------------------   --------------------------
      Lat      Lon         Lat       Lon         X        Y       Total
-- -------- ---------   --------- ---------   -------  -------   -------
 1  44.5006  -87.9107 ->  44.5006  -87.9107:     0.05     0.07      0.09
 2  44.5006  -86.2065 ->  44.5006  -86.2065:    -0.07     0.07      0.10
 3  41.6132  -86.2065 ->  41.6132  -86.2065:    -0.07    -0.07      0.10
 4  41.6132  -87.9107 ->  41.6132  -87.9107:     0.05    -0.07      0.09

------------------------------------------------------------------------
 *  control point latitude and longitude coordinates
**  computed from the corresponding x and y coordinates

X and Y values from the '.xyc' file are corrected before
converting to lat/lon for rotation (theta), directional
stretching (xcor), translation (X0 and Y0), scale and
for the fact that Y is down int the '.xyc' file.

   x = (X+X0)*xcor*cos(-theta)-(Y+Y0)*sin(-theta)
   y =-(X+X0)*xcor*sin(-theta)-(Y+Y0)*cos(-theta)

   X = x*cos(theta)/xcor-y*sin(theta)/xcor-X0
   Y =-x*sin(theta)-y*cos(theta)-Y0

   where:
      X = column of image or distance right on tablet
      Y = row of image or distance left on tablet
     X0 = right offset to projected map origin
     Y0 = down offset to projected map origin
      x = meters left on projected map
      y = meters up on projected map

   functions used to translate systems:

int xy2ll(double x,double y,double *lat,double *lon,struct control cc)
int ll2xy(double lat,double lon,double *x,double *y,struct control cc)