mDiffFitExec.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<font style=\"font-size:96px; font-weight:bolder; color:#0040a0\"><img src=\"http://montage.ipac.caltech.edu/docs/M51_logo.png\" alt=\"M\" style=\"float: left; padding: 25px 30px 25px 0px;\" /></font>\n",
    "\n",
    "<i><b>Montage</b> Montage is an astronomical image toolkit with components for reprojection, background matching, coaddition and visualization of FITS files. It can be used as a set of command-line tools (Linux, OS X and Windows), C library calls (Linux and OS X) and as Python binary extension modules.\n",
    "\n",
    "The Montage source is written in ANSI-C and code can be downloaded from GitHub ( https://github.com/Caltech-IPAC/Montage ). The Python package can be installed from PyPI (\"</i>pip install MontagePy<i>\"). The package has no external dependencies. See http://montage.ipac.caltech.edu/ for details on the design and applications of Montage.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# MontagePy.main modules: mDiffFitExec\n",
    "\n",
    "The Montage modules are generally used as steps in a workflow to create a mosaic of a set of input images.  These steps are: determine the geometry of the mosaic on the sky, reproject the images to a common frame and spatial sampling; rectify the backgrounds to a common level, and coadd the images into a mosaic.  This page illustrates the use of one Montage module, mDiffFitExec, which differences and fits a set of image pairs (image overlaps).\n",
    "\n",
    "Visit <a href=\"Mosaic.ipynb\">Building a Mosaic with Montage</a> to see how mDiffFitExec is used as part of a workflow to creage a mosaic (or the <a href=\"Mosaic_oneshot.ipynb\"> one shot </a> version if you just want to see the commands).  See the complete list of Montage Notebooks <a href=\"http://montage.ipac.caltech.edu/MontageNotebooks\">here</a>.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Help on built-in function mDiffFitExec in module MontagePy.main:\n",
      "\n",
      "mDiffFitExec(...)\n",
      "    mDiffFitExec loops over a list of image overlaps, computing then fitting the differences.\n",
      "    \n",
      "    Parameters\n",
      "    ----------\n",
      "    path : str\n",
      "        Path to raw image directory.\n",
      "    tblfile : str\n",
      "        Table file list of input images.\n",
      "    template : str\n",
      "        FITS header file used to define the desired output.\n",
      "    diffdir : str\n",
      "        Directory for temporary difference files.\n",
      "    fitfile : str\n",
      "        Table file for output difference fit info.\n",
      "    keepAll : bool, optional\n",
      "        Flag to keep temporary difference images.\n",
      "    levelOnly : bool, optional\n",
      "        Flag to fit level of differences only, not slopes.\n",
      "    noAreas : bool, optional\n",
      "        Flag indicating there are no area images.\n",
      "    debug : int, optional\n",
      "        Debugging output flag.\n",
      "    \n",
      "    \n",
      "    Returns\n",
      "    -------\n",
      "    status : int\n",
      "        Return status (0: OK, 1:ERROR).\n",
      "    msg : str\n",
      "        Return message (for errors).\n",
      "    count : int\n",
      "        Total number of images.\n",
      "    diff_failed : int\n",
      "        Count where difference failed.\n",
      "    fit_failed : int\n",
      "        Count where fitting failed.\n",
      "    warning : int\n",
      "        Count where fitting produced warnings.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "from MontagePy.main import mDiffFitExec, mViewer\n",
    "\n",
    "help(mDiffFitExec)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## mDiffFitExec Example\n",
    "\n",
    "mDiffFitExec is a general utility for generating and fitting a set of image differences. It's most common use is as part of the process of evaluating overlap area in a set of images as part of background rectification. That process consists of taking the output of the overlaps calculation (mOverlaps) and for each overlap pair doing the difference (mDiff), fitting it to determine offset between those two images (mFitplane).\n",
    "\n",
    "There are two ways to do this in bulk. We can either generate all the differences first, then fit them all or we can loop over the differences, generating and fitting that difference before moving on to the next. This first approach uses more disks space as all the differences have to be stored in between steps. The second can clean up as it goes.\n",
    "\n",
    "Obviously the second approach is preferable unless you want to keep the differences (possibly for quality evaluation).\n",
    "\n",
    "mDiffFitExec is an instantiation of the second approach. It runs mDiff on every pair in the input list and then fits it with a plane.\n",
    "\n",
    "At this point in the processing, we already have a set of \"like\" images (same projection), and have made a list of all the overlaps (normally using mOverlaps).  mDiffFitExec loops over these differences, generating each difference (in a working directory) and fitting a plane to it.\n",
    "\n",
    "mDiffFitExec cleans up the working files as it goes."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'status': '0', 'count': 128, 'diff_failed': 0, 'fit_failed': 0, 'warning': 0}\n"
     ]
    }
   ],
   "source": [
    "rtn = mDiffFitExec('M17/projected', \n",
    "                   'M17/diffs.tbl', \n",
    "                   'M17/M17.hdr',\n",
    "                   'M17/diffs',\n",
    "                   'work/M17/fits.tbl')\n",
    "print(rtn)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The output is a set of planes representing the difference between two original images.  This set will be sent to mBgModel for modelling into a set of corrections to the original images to bring the backgrounds in line."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>plus</th>\n",
       "      <th>minus</th>\n",
       "      <th>a</th>\n",
       "      <th>b</th>\n",
       "      <th>c</th>\n",
       "      <th>crpix1</th>\n",
       "      <th>crpix2</th>\n",
       "      <th>xmin</th>\n",
       "      <th>xmax</th>\n",
       "      <th>ymin</th>\n",
       "      <th>ymax</th>\n",
       "      <th>xcenter</th>\n",
       "      <th>ycenter</th>\n",
       "      <th>npixel</th>\n",
       "      <th>rms</th>\n",
       "      <th>boxx</th>\n",
       "      <th>boxy</th>\n",
       "      <th>boxwidth</th>\n",
       "      <th>boxheight</th>\n",
       "      <th>boxang</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>8</td>\n",
       "      <td>1.268790e-03</td>\n",
       "      <td>4.644410e-04</td>\n",
       "      <td>0.410733</td>\n",
       "      <td>-1711.0</td>\n",
       "      <td>1756.5</td>\n",
       "      <td>1711</td>\n",
       "      <td>1789</td>\n",
       "      <td>-1756</td>\n",
       "      <td>-935</td>\n",
       "      <td>1750.15</td>\n",
       "      <td>-1343.07</td>\n",
       "      <td>52864</td>\n",
       "      <td>7.362550e-01</td>\n",
       "      <td>1749.9</td>\n",
       "      <td>-1346.0</td>\n",
       "      <td>821.2</td>\n",
       "      <td>77.3</td>\n",
       "      <td>-90.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0</td>\n",
       "      <td>9</td>\n",
       "      <td>-2.523660e-03</td>\n",
       "      <td>1.653820e-03</td>\n",
       "      <td>6.969540</td>\n",
       "      <td>-1713.0</td>\n",
       "      <td>988.5</td>\n",
       "      <td>1713</td>\n",
       "      <td>1789</td>\n",
       "      <td>-988</td>\n",
       "      <td>-733</td>\n",
       "      <td>1751.01</td>\n",
       "      <td>-863.07</td>\n",
       "      <td>16232</td>\n",
       "      <td>7.125360e-01</td>\n",
       "      <td>1751.1</td>\n",
       "      <td>-861.0</td>\n",
       "      <td>255.3</td>\n",
       "      <td>76.9</td>\n",
       "      <td>-90.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0</td>\n",
       "      <td>13</td>\n",
       "      <td>2.453210e-03</td>\n",
       "      <td>3.935060e-03</td>\n",
       "      <td>0.277418</td>\n",
       "      <td>-1286.0</td>\n",
       "      <td>986.5</td>\n",
       "      <td>1286</td>\n",
       "      <td>1361</td>\n",
       "      <td>-986</td>\n",
       "      <td>-732</td>\n",
       "      <td>1323.43</td>\n",
       "      <td>-862.04</td>\n",
       "      <td>16218</td>\n",
       "      <td>7.174900e-01</td>\n",
       "      <td>1323.5</td>\n",
       "      <td>-859.5</td>\n",
       "      <td>254.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0</td>\n",
       "      <td>14</td>\n",
       "      <td>-2.490360e-03</td>\n",
       "      <td>8.991680e-04</td>\n",
       "      <td>4.436860</td>\n",
       "      <td>-1283.0</td>\n",
       "      <td>1755.5</td>\n",
       "      <td>1283</td>\n",
       "      <td>1361</td>\n",
       "      <td>-1755</td>\n",
       "      <td>-933</td>\n",
       "      <td>1322.70</td>\n",
       "      <td>-1342.62</td>\n",
       "      <td>52649</td>\n",
       "      <td>7.105830e-01</td>\n",
       "      <td>1322.3</td>\n",
       "      <td>-1344.5</td>\n",
       "      <td>822.2</td>\n",
       "      <td>77.3</td>\n",
       "      <td>89.9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0</td>\n",
       "      <td>21</td>\n",
       "      <td>-3.303240e-10</td>\n",
       "      <td>5.376120e-09</td>\n",
       "      <td>-1.130130</td>\n",
       "      <td>-1286.0</td>\n",
       "      <td>786.5</td>\n",
       "      <td>1286</td>\n",
       "      <td>1790</td>\n",
       "      <td>-786</td>\n",
       "      <td>-732</td>\n",
       "      <td>1537.86</td>\n",
       "      <td>-759.37</td>\n",
       "      <td>27003</td>\n",
       "      <td>7.587320e-06</td>\n",
       "      <td>1538.0</td>\n",
       "      <td>-759.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>504.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>0</td>\n",
       "      <td>24</td>\n",
       "      <td>-1.489270e-10</td>\n",
       "      <td>-5.884170e-09</td>\n",
       "      <td>0.561656</td>\n",
       "      <td>-1284.0</td>\n",
       "      <td>1756.5</td>\n",
       "      <td>1284</td>\n",
       "      <td>1787</td>\n",
       "      <td>-1756</td>\n",
       "      <td>-1702</td>\n",
       "      <td>1535.59</td>\n",
       "      <td>-1729.23</td>\n",
       "      <td>27090</td>\n",
       "      <td>7.108770e-06</td>\n",
       "      <td>1535.5</td>\n",
       "      <td>-1729.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>0</td>\n",
       "      <td>29</td>\n",
       "      <td>-1.230960e-03</td>\n",
       "      <td>5.115690e-04</td>\n",
       "      <td>-79.055300</td>\n",
       "      <td>-1710.0</td>\n",
       "      <td>994.5</td>\n",
       "      <td>1710</td>\n",
       "      <td>1790</td>\n",
       "      <td>-994</td>\n",
       "      <td>-733</td>\n",
       "      <td>1749.78</td>\n",
       "      <td>-865.65</td>\n",
       "      <td>16584</td>\n",
       "      <td>7.530540e-01</td>\n",
       "      <td>1749.9</td>\n",
       "      <td>-864.0</td>\n",
       "      <td>261.3</td>\n",
       "      <td>79.3</td>\n",
       "      <td>-90.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>0</td>\n",
       "      <td>31</td>\n",
       "      <td>2.329150e-03</td>\n",
       "      <td>1.100920e-03</td>\n",
       "      <td>-82.140500</td>\n",
       "      <td>-1708.0</td>\n",
       "      <td>1756.5</td>\n",
       "      <td>1708</td>\n",
       "      <td>1789</td>\n",
       "      <td>-1756</td>\n",
       "      <td>-941</td>\n",
       "      <td>1748.71</td>\n",
       "      <td>-1348.03</td>\n",
       "      <td>53097</td>\n",
       "      <td>7.707710e-01</td>\n",
       "      <td>1748.8</td>\n",
       "      <td>-1349.0</td>\n",
       "      <td>815.2</td>\n",
       "      <td>79.9</td>\n",
       "      <td>89.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>1</td>\n",
       "      <td>6</td>\n",
       "      <td>1.756510e-04</td>\n",
       "      <td>-2.708390e-04</td>\n",
       "      <td>-3.833850</td>\n",
       "      <td>426.0</td>\n",
       "      <td>789.5</td>\n",
       "      <td>-426</td>\n",
       "      <td>-351</td>\n",
       "      <td>-789</td>\n",
       "      <td>32</td>\n",
       "      <td>-387.36</td>\n",
       "      <td>-347.45</td>\n",
       "      <td>42063</td>\n",
       "      <td>9.129390e-01</td>\n",
       "      <td>-388.5</td>\n",
       "      <td>-378.5</td>\n",
       "      <td>822.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>1</td>\n",
       "      <td>11</td>\n",
       "      <td>-2.720230e-03</td>\n",
       "      <td>-5.054520e-04</td>\n",
       "      <td>-0.537130</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>789.5</td>\n",
       "      <td>2</td>\n",
       "      <td>78</td>\n",
       "      <td>-789</td>\n",
       "      <td>31</td>\n",
       "      <td>40.09</td>\n",
       "      <td>-390.72</td>\n",
       "      <td>47921</td>\n",
       "      <td>8.727160e-01</td>\n",
       "      <td>40.0</td>\n",
       "      <td>-379.0</td>\n",
       "      <td>821.0</td>\n",
       "      <td>76.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>1</td>\n",
       "      <td>20</td>\n",
       "      <td>-5.109730e-19</td>\n",
       "      <td>-2.983720e-16</td>\n",
       "      <td>0.763428</td>\n",
       "      <td>425.0</td>\n",
       "      <td>789.5</td>\n",
       "      <td>-425</td>\n",
       "      <td>78</td>\n",
       "      <td>-789</td>\n",
       "      <td>-736</td>\n",
       "      <td>-172.39</td>\n",
       "      <td>-763.09</td>\n",
       "      <td>22413</td>\n",
       "      <td>4.641070e-15</td>\n",
       "      <td>-173.5</td>\n",
       "      <td>-763.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>1</td>\n",
       "      <td>27</td>\n",
       "      <td>1.084200e-19</td>\n",
       "      <td>-5.074070e-17</td>\n",
       "      <td>0.691162</td>\n",
       "      <td>426.0</td>\n",
       "      <td>-180.5</td>\n",
       "      <td>-426</td>\n",
       "      <td>77</td>\n",
       "      <td>180</td>\n",
       "      <td>233</td>\n",
       "      <td>-174.68</td>\n",
       "      <td>207.06</td>\n",
       "      <td>23318</td>\n",
       "      <td>7.921560e-16</td>\n",
       "      <td>-174.5</td>\n",
       "      <td>207.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>1</td>\n",
       "      <td>34</td>\n",
       "      <td>-4.281630e-03</td>\n",
       "      <td>6.650880e-04</td>\n",
       "      <td>-4.570030</td>\n",
       "      <td>426.0</td>\n",
       "      <td>20.5</td>\n",
       "      <td>-426</td>\n",
       "      <td>-351</td>\n",
       "      <td>-20</td>\n",
       "      <td>233</td>\n",
       "      <td>-388.07</td>\n",
       "      <td>106.65</td>\n",
       "      <td>14652</td>\n",
       "      <td>8.112380e-01</td>\n",
       "      <td>-388.5</td>\n",
       "      <td>106.5</td>\n",
       "      <td>254.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>1</td>\n",
       "      <td>41</td>\n",
       "      <td>-9.965350e-03</td>\n",
       "      <td>3.869090e-03</td>\n",
       "      <td>-0.521526</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>21.5</td>\n",
       "      <td>2</td>\n",
       "      <td>77</td>\n",
       "      <td>-21</td>\n",
       "      <td>233</td>\n",
       "      <td>39.46</td>\n",
       "      <td>110.87</td>\n",
       "      <td>15314</td>\n",
       "      <td>9.693480e-01</td>\n",
       "      <td>39.5</td>\n",
       "      <td>106.0</td>\n",
       "      <td>255.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>2</td>\n",
       "      <td>5</td>\n",
       "      <td>5.910760e-03</td>\n",
       "      <td>9.426410e-04</td>\n",
       "      <td>5.666030</td>\n",
       "      <td>854.0</td>\n",
       "      <td>1757.5</td>\n",
       "      <td>-854</td>\n",
       "      <td>-777</td>\n",
       "      <td>-1757</td>\n",
       "      <td>-937</td>\n",
       "      <td>-815.22</td>\n",
       "      <td>-1354.16</td>\n",
       "      <td>49825</td>\n",
       "      <td>7.697600e-01</td>\n",
       "      <td>-815.4</td>\n",
       "      <td>-1347.5</td>\n",
       "      <td>820.1</td>\n",
       "      <td>76.4</td>\n",
       "      <td>90.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>2</td>\n",
       "      <td>6</td>\n",
       "      <td>-7.671010e-04</td>\n",
       "      <td>1.611060e-03</td>\n",
       "      <td>0.543132</td>\n",
       "      <td>854.0</td>\n",
       "      <td>990.5</td>\n",
       "      <td>-854</td>\n",
       "      <td>-778</td>\n",
       "      <td>-990</td>\n",
       "      <td>-734</td>\n",
       "      <td>-815.57</td>\n",
       "      <td>-864.23</td>\n",
       "      <td>14285</td>\n",
       "      <td>9.272300e-01</td>\n",
       "      <td>-816.0</td>\n",
       "      <td>-862.5</td>\n",
       "      <td>256.0</td>\n",
       "      <td>76.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>2</td>\n",
       "      <td>18</td>\n",
       "      <td>4.852270e-10</td>\n",
       "      <td>3.118610e-09</td>\n",
       "      <td>1.600290</td>\n",
       "      <td>1282.0</td>\n",
       "      <td>788.5</td>\n",
       "      <td>-1282</td>\n",
       "      <td>-778</td>\n",
       "      <td>-788</td>\n",
       "      <td>-734</td>\n",
       "      <td>-1029.37</td>\n",
       "      <td>-761.30</td>\n",
       "      <td>26871</td>\n",
       "      <td>7.403790e-06</td>\n",
       "      <td>-1030.0</td>\n",
       "      <td>-761.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>504.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>2</td>\n",
       "      <td>25</td>\n",
       "      <td>-3.896340e-04</td>\n",
       "      <td>9.108010e-04</td>\n",
       "      <td>6.162890</td>\n",
       "      <td>1281.0</td>\n",
       "      <td>1758.5</td>\n",
       "      <td>-1281</td>\n",
       "      <td>-1238</td>\n",
       "      <td>-1758</td>\n",
       "      <td>-933</td>\n",
       "      <td>-1260.25</td>\n",
       "      <td>-1335.75</td>\n",
       "      <td>24848</td>\n",
       "      <td>8.814390e-01</td>\n",
       "      <td>-1260.0</td>\n",
       "      <td>-1346.0</td>\n",
       "      <td>825.1</td>\n",
       "      <td>42.7</td>\n",
       "      <td>-89.9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>2</td>\n",
       "      <td>26</td>\n",
       "      <td>6.081620e-03</td>\n",
       "      <td>1.740740e-03</td>\n",
       "      <td>13.644000</td>\n",
       "      <td>1282.0</td>\n",
       "      <td>986.5</td>\n",
       "      <td>-1282</td>\n",
       "      <td>-1240</td>\n",
       "      <td>-986</td>\n",
       "      <td>-735</td>\n",
       "      <td>-1260.26</td>\n",
       "      <td>-864.96</td>\n",
       "      <td>7125</td>\n",
       "      <td>1.006160e+00</td>\n",
       "      <td>-1261.0</td>\n",
       "      <td>-861.0</td>\n",
       "      <td>251.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>2</td>\n",
       "      <td>28</td>\n",
       "      <td>-3.959880e-20</td>\n",
       "      <td>-1.304840e-16</td>\n",
       "      <td>-0.597885</td>\n",
       "      <td>1280.0</td>\n",
       "      <td>1758.5</td>\n",
       "      <td>-1280</td>\n",
       "      <td>-777</td>\n",
       "      <td>-1758</td>\n",
       "      <td>-1704</td>\n",
       "      <td>-1028.23</td>\n",
       "      <td>-1731.35</td>\n",
       "      <td>24404</td>\n",
       "      <td>2.032750e-15</td>\n",
       "      <td>-1028.5</td>\n",
       "      <td>-1731.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>3</td>\n",
       "      <td>11</td>\n",
       "      <td>2.111740e-03</td>\n",
       "      <td>8.402230e-04</td>\n",
       "      <td>-2.689850</td>\n",
       "      <td>-431.0</td>\n",
       "      <td>784.5</td>\n",
       "      <td>431</td>\n",
       "      <td>505</td>\n",
       "      <td>-784</td>\n",
       "      <td>31</td>\n",
       "      <td>468.00</td>\n",
       "      <td>-378.34</td>\n",
       "      <td>50863</td>\n",
       "      <td>7.444600e-01</td>\n",
       "      <td>468.0</td>\n",
       "      <td>-376.5</td>\n",
       "      <td>816.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>3</td>\n",
       "      <td>13</td>\n",
       "      <td>3.627310e-03</td>\n",
       "      <td>1.314690e-03</td>\n",
       "      <td>-2.519670</td>\n",
       "      <td>-858.0</td>\n",
       "      <td>784.5</td>\n",
       "      <td>858</td>\n",
       "      <td>935</td>\n",
       "      <td>-784</td>\n",
       "      <td>37</td>\n",
       "      <td>896.00</td>\n",
       "      <td>-382.34</td>\n",
       "      <td>50308</td>\n",
       "      <td>7.361090e-01</td>\n",
       "      <td>896.5</td>\n",
       "      <td>-373.5</td>\n",
       "      <td>822.1</td>\n",
       "      <td>76.1</td>\n",
       "      <td>-90.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>3</td>\n",
       "      <td>17</td>\n",
       "      <td>2.483370e-10</td>\n",
       "      <td>-1.152390e-09</td>\n",
       "      <td>-0.965677</td>\n",
       "      <td>-431.0</td>\n",
       "      <td>784.5</td>\n",
       "      <td>431</td>\n",
       "      <td>934</td>\n",
       "      <td>-784</td>\n",
       "      <td>-731</td>\n",
       "      <td>682.61</td>\n",
       "      <td>-757.99</td>\n",
       "      <td>27199</td>\n",
       "      <td>7.623260e-06</td>\n",
       "      <td>682.5</td>\n",
       "      <td>-758.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>3</td>\n",
       "      <td>23</td>\n",
       "      <td>-3.036630e-11</td>\n",
       "      <td>1.690620e-09</td>\n",
       "      <td>-1.258960</td>\n",
       "      <td>-432.0</td>\n",
       "      <td>-185.5</td>\n",
       "      <td>432</td>\n",
       "      <td>935</td>\n",
       "      <td>185</td>\n",
       "      <td>238</td>\n",
       "      <td>683.50</td>\n",
       "      <td>212.00</td>\n",
       "      <td>27216</td>\n",
       "      <td>6.063610e-06</td>\n",
       "      <td>683.5</td>\n",
       "      <td>212.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>3</td>\n",
       "      <td>38</td>\n",
       "      <td>3.505620e-03</td>\n",
       "      <td>1.811340e-04</td>\n",
       "      <td>-3.260770</td>\n",
       "      <td>-860.0</td>\n",
       "      <td>15.5</td>\n",
       "      <td>860</td>\n",
       "      <td>935</td>\n",
       "      <td>-15</td>\n",
       "      <td>238</td>\n",
       "      <td>897.12</td>\n",
       "      <td>111.32</td>\n",
       "      <td>15496</td>\n",
       "      <td>7.787990e-01</td>\n",
       "      <td>897.5</td>\n",
       "      <td>111.5</td>\n",
       "      <td>254.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>3</td>\n",
       "      <td>41</td>\n",
       "      <td>5.486890e-03</td>\n",
       "      <td>1.928290e-03</td>\n",
       "      <td>-4.866100</td>\n",
       "      <td>-431.0</td>\n",
       "      <td>21.5</td>\n",
       "      <td>431</td>\n",
       "      <td>505</td>\n",
       "      <td>-21</td>\n",
       "      <td>238</td>\n",
       "      <td>468.22</td>\n",
       "      <td>108.40</td>\n",
       "      <td>15892</td>\n",
       "      <td>7.351890e-01</td>\n",
       "      <td>468.0</td>\n",
       "      <td>108.5</td>\n",
       "      <td>260.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>4</td>\n",
       "      <td>15</td>\n",
       "      <td>4.008290e-03</td>\n",
       "      <td>-9.678930e-04</td>\n",
       "      <td>6.493410</td>\n",
       "      <td>1747.0</td>\n",
       "      <td>-183.5</td>\n",
       "      <td>-1747</td>\n",
       "      <td>-1672</td>\n",
       "      <td>183</td>\n",
       "      <td>1005</td>\n",
       "      <td>-1708.71</td>\n",
       "      <td>591.93</td>\n",
       "      <td>43190</td>\n",
       "      <td>9.386600e-01</td>\n",
       "      <td>-1709.7</td>\n",
       "      <td>594.6</td>\n",
       "      <td>822.0</td>\n",
       "      <td>73.7</td>\n",
       "      <td>-89.9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>4</td>\n",
       "      <td>18</td>\n",
       "      <td>-3.352750e-04</td>\n",
       "      <td>-2.270560e-03</td>\n",
       "      <td>-3.924300</td>\n",
       "      <td>1283.0</td>\n",
       "      <td>16.5</td>\n",
       "      <td>-1283</td>\n",
       "      <td>-1242</td>\n",
       "      <td>-16</td>\n",
       "      <td>234</td>\n",
       "      <td>-1262.86</td>\n",
       "      <td>106.53</td>\n",
       "      <td>7062</td>\n",
       "      <td>1.012440e+00</td>\n",
       "      <td>-1262.5</td>\n",
       "      <td>109.0</td>\n",
       "      <td>251.0</td>\n",
       "      <td>41.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>4</td>\n",
       "      <td>19</td>\n",
       "      <td>1.435180e-16</td>\n",
       "      <td>-1.116540e-14</td>\n",
       "      <td>-0.190674</td>\n",
       "      <td>1747.0</td>\n",
       "      <td>-952.5</td>\n",
       "      <td>-1747</td>\n",
       "      <td>-1244</td>\n",
       "      <td>974</td>\n",
       "      <td>983</td>\n",
       "      <td>-1493.10</td>\n",
       "      <td>978.90</td>\n",
       "      <td>47</td>\n",
       "      <td>3.698180e-16</td>\n",
       "      <td>-1495.5</td>\n",
       "      <td>979.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>4</td>\n",
       "      <td>26</td>\n",
       "      <td>2.168400e-18</td>\n",
       "      <td>4.445230e-18</td>\n",
       "      <td>1.108510</td>\n",
       "      <td>1745.0</td>\n",
       "      <td>17.5</td>\n",
       "      <td>-1745</td>\n",
       "      <td>-1242</td>\n",
       "      <td>-17</td>\n",
       "      <td>36</td>\n",
       "      <td>-1494.00</td>\n",
       "      <td>9.26</td>\n",
       "      <td>24954</td>\n",
       "      <td>3.432580e-16</td>\n",
       "      <td>-1493.5</td>\n",
       "      <td>9.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98</th>\n",
       "      <td>25</td>\n",
       "      <td>37</td>\n",
       "      <td>2.993400e-03</td>\n",
       "      <td>-9.639940e-04</td>\n",
       "      <td>1.930760</td>\n",
       "      <td>1743.0</td>\n",
       "      <td>1756.5</td>\n",
       "      <td>-1743</td>\n",
       "      <td>-1667</td>\n",
       "      <td>-1756</td>\n",
       "      <td>-934</td>\n",
       "      <td>-1704.44</td>\n",
       "      <td>-1368.87</td>\n",
       "      <td>43558</td>\n",
       "      <td>9.570520e-01</td>\n",
       "      <td>-1705.1</td>\n",
       "      <td>-1345.5</td>\n",
       "      <td>822.1</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99</th>\n",
       "      <td>26</td>\n",
       "      <td>37</td>\n",
       "      <td>2.132510e-03</td>\n",
       "      <td>-7.878530e-04</td>\n",
       "      <td>1.137900</td>\n",
       "      <td>1743.0</td>\n",
       "      <td>987.5</td>\n",
       "      <td>-1743</td>\n",
       "      <td>-1669</td>\n",
       "      <td>-987</td>\n",
       "      <td>-732</td>\n",
       "      <td>-1704.94</td>\n",
       "      <td>-856.26</td>\n",
       "      <td>13175</td>\n",
       "      <td>8.965520e-01</td>\n",
       "      <td>-1706.0</td>\n",
       "      <td>-860.0</td>\n",
       "      <td>255.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>26</td>\n",
       "      <td>42</td>\n",
       "      <td>4.095660e-03</td>\n",
       "      <td>-1.870220e-03</td>\n",
       "      <td>4.653600</td>\n",
       "      <td>1745.0</td>\n",
       "      <td>786.5</td>\n",
       "      <td>-1745</td>\n",
       "      <td>-1670</td>\n",
       "      <td>-786</td>\n",
       "      <td>35</td>\n",
       "      <td>-1705.76</td>\n",
       "      <td>-365.01</td>\n",
       "      <td>44054</td>\n",
       "      <td>8.743800e-01</td>\n",
       "      <td>-1707.4</td>\n",
       "      <td>-375.5</td>\n",
       "      <td>822.0</td>\n",
       "      <td>74.3</td>\n",
       "      <td>90.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>101</th>\n",
       "      <td>27</td>\n",
       "      <td>33</td>\n",
       "      <td>2.527970e-18</td>\n",
       "      <td>-4.167670e-16</td>\n",
       "      <td>1.673900</td>\n",
       "      <td>426.0</td>\n",
       "      <td>-1150.5</td>\n",
       "      <td>-426</td>\n",
       "      <td>77</td>\n",
       "      <td>1150</td>\n",
       "      <td>1203</td>\n",
       "      <td>-173.40</td>\n",
       "      <td>1177.04</td>\n",
       "      <td>21881</td>\n",
       "      <td>6.513750e-15</td>\n",
       "      <td>-174.5</td>\n",
       "      <td>1177.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>102</th>\n",
       "      <td>27</td>\n",
       "      <td>34</td>\n",
       "      <td>3.071140e-03</td>\n",
       "      <td>1.595790e-03</td>\n",
       "      <td>-2.654820</td>\n",
       "      <td>426.0</td>\n",
       "      <td>-180.5</td>\n",
       "      <td>-426</td>\n",
       "      <td>-351</td>\n",
       "      <td>180</td>\n",
       "      <td>1002</td>\n",
       "      <td>-389.66</td>\n",
       "      <td>570.28</td>\n",
       "      <td>44329</td>\n",
       "      <td>9.364340e-01</td>\n",
       "      <td>-388.5</td>\n",
       "      <td>591.5</td>\n",
       "      <td>822.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>103</th>\n",
       "      <td>27</td>\n",
       "      <td>35</td>\n",
       "      <td>1.628030e-03</td>\n",
       "      <td>1.117010e-03</td>\n",
       "      <td>-0.391749</td>\n",
       "      <td>-3.0</td>\n",
       "      <td>-948.5</td>\n",
       "      <td>3</td>\n",
       "      <td>77</td>\n",
       "      <td>948</td>\n",
       "      <td>1203</td>\n",
       "      <td>40.62</td>\n",
       "      <td>1076.31</td>\n",
       "      <td>12992</td>\n",
       "      <td>1.003060e+00</td>\n",
       "      <td>40.0</td>\n",
       "      <td>1076.0</td>\n",
       "      <td>255.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104</th>\n",
       "      <td>27</td>\n",
       "      <td>41</td>\n",
       "      <td>2.676220e-03</td>\n",
       "      <td>1.519470e-03</td>\n",
       "      <td>-1.326430</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>-180.5</td>\n",
       "      <td>2</td>\n",
       "      <td>77</td>\n",
       "      <td>180</td>\n",
       "      <td>1001</td>\n",
       "      <td>40.10</td>\n",
       "      <td>570.35</td>\n",
       "      <td>46199</td>\n",
       "      <td>8.684200e-01</td>\n",
       "      <td>39.5</td>\n",
       "      <td>591.0</td>\n",
       "      <td>821.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>105</th>\n",
       "      <td>27</td>\n",
       "      <td>43</td>\n",
       "      <td>6.693140e-03</td>\n",
       "      <td>3.815910e-03</td>\n",
       "      <td>-3.391830</td>\n",
       "      <td>426.0</td>\n",
       "      <td>-949.5</td>\n",
       "      <td>-426</td>\n",
       "      <td>-352</td>\n",
       "      <td>949</td>\n",
       "      <td>1203</td>\n",
       "      <td>-390.89</td>\n",
       "      <td>1078.05</td>\n",
       "      <td>12730</td>\n",
       "      <td>1.083480e+00</td>\n",
       "      <td>-389.0</td>\n",
       "      <td>1076.5</td>\n",
       "      <td>254.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>106</th>\n",
       "      <td>29</td>\n",
       "      <td>30</td>\n",
       "      <td>-6.003610e-03</td>\n",
       "      <td>7.217340e-03</td>\n",
       "      <td>92.591200</td>\n",
       "      <td>-1716.0</td>\n",
       "      <td>18.5</td>\n",
       "      <td>1716</td>\n",
       "      <td>1798</td>\n",
       "      <td>-18</td>\n",
       "      <td>28</td>\n",
       "      <td>1757.75</td>\n",
       "      <td>4.26</td>\n",
       "      <td>3000</td>\n",
       "      <td>9.457370e-01</td>\n",
       "      <td>1757.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>47.0</td>\n",
       "      <td>82.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>107</th>\n",
       "      <td>29</td>\n",
       "      <td>31</td>\n",
       "      <td>3.228370e-09</td>\n",
       "      <td>3.008890e-09</td>\n",
       "      <td>2.299610</td>\n",
       "      <td>-1710.0</td>\n",
       "      <td>994.5</td>\n",
       "      <td>1710</td>\n",
       "      <td>1798</td>\n",
       "      <td>-994</td>\n",
       "      <td>-941</td>\n",
       "      <td>1754.05</td>\n",
       "      <td>-968.01</td>\n",
       "      <td>4791</td>\n",
       "      <td>7.628790e-06</td>\n",
       "      <td>1754.0</td>\n",
       "      <td>-968.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>88.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>108</th>\n",
       "      <td>30</td>\n",
       "      <td>44</td>\n",
       "      <td>1.178090e-15</td>\n",
       "      <td>1.691360e-17</td>\n",
       "      <td>-3.227810</td>\n",
       "      <td>-1718.0</td>\n",
       "      <td>-951.5</td>\n",
       "      <td>1718</td>\n",
       "      <td>1798</td>\n",
       "      <td>951</td>\n",
       "      <td>1004</td>\n",
       "      <td>1757.89</td>\n",
       "      <td>978.04</td>\n",
       "      <td>3747</td>\n",
       "      <td>2.755180e-14</td>\n",
       "      <td>1758.0</td>\n",
       "      <td>978.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>80.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>109</th>\n",
       "      <td>30</td>\n",
       "      <td>45</td>\n",
       "      <td>5.933540e-03</td>\n",
       "      <td>-1.243340e-03</td>\n",
       "      <td>-13.183500</td>\n",
       "      <td>-1716.0</td>\n",
       "      <td>-183.5</td>\n",
       "      <td>1716</td>\n",
       "      <td>1793</td>\n",
       "      <td>183</td>\n",
       "      <td>1004</td>\n",
       "      <td>1754.93</td>\n",
       "      <td>607.41</td>\n",
       "      <td>46514</td>\n",
       "      <td>8.592740e-01</td>\n",
       "      <td>1754.9</td>\n",
       "      <td>594.0</td>\n",
       "      <td>821.2</td>\n",
       "      <td>76.3</td>\n",
       "      <td>-90.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>110</th>\n",
       "      <td>32</td>\n",
       "      <td>35</td>\n",
       "      <td>1.212370e-03</td>\n",
       "      <td>-2.488770e-03</td>\n",
       "      <td>2.879520</td>\n",
       "      <td>-433.0</td>\n",
       "      <td>-1155.5</td>\n",
       "      <td>433</td>\n",
       "      <td>506</td>\n",
       "      <td>1155</td>\n",
       "      <td>1797</td>\n",
       "      <td>468.94</td>\n",
       "      <td>1573.64</td>\n",
       "      <td>27217</td>\n",
       "      <td>9.886340e-01</td>\n",
       "      <td>469.5</td>\n",
       "      <td>1476.5</td>\n",
       "      <td>642.0</td>\n",
       "      <td>73.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111</th>\n",
       "      <td>32</td>\n",
       "      <td>36</td>\n",
       "      <td>-4.557880e-03</td>\n",
       "      <td>2.532880e-03</td>\n",
       "      <td>0.350670</td>\n",
       "      <td>-862.0</td>\n",
       "      <td>-1155.5</td>\n",
       "      <td>862</td>\n",
       "      <td>936</td>\n",
       "      <td>1155</td>\n",
       "      <td>1797</td>\n",
       "      <td>899.31</td>\n",
       "      <td>1465.39</td>\n",
       "      <td>37138</td>\n",
       "      <td>8.199690e-01</td>\n",
       "      <td>899.0</td>\n",
       "      <td>1476.5</td>\n",
       "      <td>642.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>112</th>\n",
       "      <td>33</td>\n",
       "      <td>35</td>\n",
       "      <td>5.023290e-03</td>\n",
       "      <td>2.110940e-03</td>\n",
       "      <td>-3.150380</td>\n",
       "      <td>-3.0</td>\n",
       "      <td>-1150.5</td>\n",
       "      <td>3</td>\n",
       "      <td>77</td>\n",
       "      <td>1150</td>\n",
       "      <td>1797</td>\n",
       "      <td>40.41</td>\n",
       "      <td>1487.34</td>\n",
       "      <td>35991</td>\n",
       "      <td>8.733230e-01</td>\n",
       "      <td>40.0</td>\n",
       "      <td>1474.0</td>\n",
       "      <td>647.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>113</th>\n",
       "      <td>33</td>\n",
       "      <td>43</td>\n",
       "      <td>3.743590e-03</td>\n",
       "      <td>1.023310e-03</td>\n",
       "      <td>-3.199060</td>\n",
       "      <td>426.0</td>\n",
       "      <td>-1150.5</td>\n",
       "      <td>-426</td>\n",
       "      <td>-352</td>\n",
       "      <td>1150</td>\n",
       "      <td>1797</td>\n",
       "      <td>-389.43</td>\n",
       "      <td>1485.53</td>\n",
       "      <td>33026</td>\n",
       "      <td>9.699490e-01</td>\n",
       "      <td>-389.0</td>\n",
       "      <td>1474.0</td>\n",
       "      <td>647.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>114</th>\n",
       "      <td>34</td>\n",
       "      <td>43</td>\n",
       "      <td>9.529120e-20</td>\n",
       "      <td>-1.970540e-17</td>\n",
       "      <td>-0.217102</td>\n",
       "      <td>855.0</td>\n",
       "      <td>-949.5</td>\n",
       "      <td>-855</td>\n",
       "      <td>-352</td>\n",
       "      <td>949</td>\n",
       "      <td>1002</td>\n",
       "      <td>-603.82</td>\n",
       "      <td>976.05</td>\n",
       "      <td>24455</td>\n",
       "      <td>3.088130e-16</td>\n",
       "      <td>-603.5</td>\n",
       "      <td>976.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>115</th>\n",
       "      <td>34</td>\n",
       "      <td>47</td>\n",
       "      <td>-5.801320e-04</td>\n",
       "      <td>-1.266650e-03</td>\n",
       "      <td>1.172310</td>\n",
       "      <td>855.0</td>\n",
       "      <td>-182.5</td>\n",
       "      <td>-855</td>\n",
       "      <td>-780</td>\n",
       "      <td>182</td>\n",
       "      <td>1002</td>\n",
       "      <td>-817.62</td>\n",
       "      <td>588.56</td>\n",
       "      <td>45873</td>\n",
       "      <td>8.761660e-01</td>\n",
       "      <td>-817.5</td>\n",
       "      <td>592.5</td>\n",
       "      <td>820.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>116</th>\n",
       "      <td>35</td>\n",
       "      <td>41</td>\n",
       "      <td>4.690440e-19</td>\n",
       "      <td>1.317850e-16</td>\n",
       "      <td>-0.366928</td>\n",
       "      <td>-3.0</td>\n",
       "      <td>-948.5</td>\n",
       "      <td>3</td>\n",
       "      <td>506</td>\n",
       "      <td>948</td>\n",
       "      <td>1001</td>\n",
       "      <td>254.91</td>\n",
       "      <td>974.98</td>\n",
       "      <td>26986</td>\n",
       "      <td>2.055550e-15</td>\n",
       "      <td>254.5</td>\n",
       "      <td>975.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>117</th>\n",
       "      <td>36</td>\n",
       "      <td>38</td>\n",
       "      <td>2.612920e-10</td>\n",
       "      <td>-1.019570e-09</td>\n",
       "      <td>1.721730</td>\n",
       "      <td>-861.0</td>\n",
       "      <td>-953.5</td>\n",
       "      <td>861</td>\n",
       "      <td>1364</td>\n",
       "      <td>953</td>\n",
       "      <td>1007</td>\n",
       "      <td>1111.91</td>\n",
       "      <td>980.63</td>\n",
       "      <td>26937</td>\n",
       "      <td>7.228180e-06</td>\n",
       "      <td>1112.5</td>\n",
       "      <td>980.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>118</th>\n",
       "      <td>36</td>\n",
       "      <td>39</td>\n",
       "      <td>-5.956260e-04</td>\n",
       "      <td>-7.872360e-04</td>\n",
       "      <td>-0.036234</td>\n",
       "      <td>-1291.0</td>\n",
       "      <td>-1154.5</td>\n",
       "      <td>1291</td>\n",
       "      <td>1366</td>\n",
       "      <td>1154</td>\n",
       "      <td>1797</td>\n",
       "      <td>1327.21</td>\n",
       "      <td>1483.97</td>\n",
       "      <td>36016</td>\n",
       "      <td>9.384890e-01</td>\n",
       "      <td>1328.3</td>\n",
       "      <td>1476.0</td>\n",
       "      <td>643.1</td>\n",
       "      <td>74.5</td>\n",
       "      <td>-90.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>119</th>\n",
       "      <td>36</td>\n",
       "      <td>45</td>\n",
       "      <td>-5.363510e-03</td>\n",
       "      <td>-4.439150e-03</td>\n",
       "      <td>11.325300</td>\n",
       "      <td>-1290.0</td>\n",
       "      <td>-953.5</td>\n",
       "      <td>1290</td>\n",
       "      <td>1365</td>\n",
       "      <td>953</td>\n",
       "      <td>1207</td>\n",
       "      <td>1327.66</td>\n",
       "      <td>1080.66</td>\n",
       "      <td>13662</td>\n",
       "      <td>9.161630e-01</td>\n",
       "      <td>1327.6</td>\n",
       "      <td>1080.5</td>\n",
       "      <td>254.4</td>\n",
       "      <td>74.4</td>\n",
       "      <td>89.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>120</th>\n",
       "      <td>37</td>\n",
       "      <td>42</td>\n",
       "      <td>1.569380e-17</td>\n",
       "      <td>9.454240e-17</td>\n",
       "      <td>0.765747</td>\n",
       "      <td>1801.0</td>\n",
       "      <td>786.5</td>\n",
       "      <td>-1801</td>\n",
       "      <td>-1670</td>\n",
       "      <td>-786</td>\n",
       "      <td>-732</td>\n",
       "      <td>-1735.80</td>\n",
       "      <td>-759.64</td>\n",
       "      <td>5746</td>\n",
       "      <td>1.590380e-15</td>\n",
       "      <td>-1735.5</td>\n",
       "      <td>-759.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>131.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>121</th>\n",
       "      <td>38</td>\n",
       "      <td>45</td>\n",
       "      <td>-3.689510e-03</td>\n",
       "      <td>-1.696570e-03</td>\n",
       "      <td>4.205390</td>\n",
       "      <td>-1289.0</td>\n",
       "      <td>-184.5</td>\n",
       "      <td>1289</td>\n",
       "      <td>1364</td>\n",
       "      <td>184</td>\n",
       "      <td>1006</td>\n",
       "      <td>1326.44</td>\n",
       "      <td>596.48</td>\n",
       "      <td>48646</td>\n",
       "      <td>8.512290e-01</td>\n",
       "      <td>1326.5</td>\n",
       "      <td>595.5</td>\n",
       "      <td>822.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>122</th>\n",
       "      <td>39</td>\n",
       "      <td>44</td>\n",
       "      <td>2.348420e-03</td>\n",
       "      <td>4.838670e-04</td>\n",
       "      <td>-3.615190</td>\n",
       "      <td>-1718.0</td>\n",
       "      <td>-1153.5</td>\n",
       "      <td>1718</td>\n",
       "      <td>1795</td>\n",
       "      <td>1153</td>\n",
       "      <td>1797</td>\n",
       "      <td>1756.99</td>\n",
       "      <td>1475.15</td>\n",
       "      <td>39354</td>\n",
       "      <td>7.699550e-01</td>\n",
       "      <td>1756.8</td>\n",
       "      <td>1475.5</td>\n",
       "      <td>644.3</td>\n",
       "      <td>76.8</td>\n",
       "      <td>-90.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123</th>\n",
       "      <td>39</td>\n",
       "      <td>45</td>\n",
       "      <td>3.354750e-10</td>\n",
       "      <td>2.823240e-09</td>\n",
       "      <td>0.253280</td>\n",
       "      <td>-1291.0</td>\n",
       "      <td>-1153.5</td>\n",
       "      <td>1291</td>\n",
       "      <td>1794</td>\n",
       "      <td>1153</td>\n",
       "      <td>1207</td>\n",
       "      <td>1542.22</td>\n",
       "      <td>1180.83</td>\n",
       "      <td>27010</td>\n",
       "      <td>6.472080e-06</td>\n",
       "      <td>1542.5</td>\n",
       "      <td>1180.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124</th>\n",
       "      <td>40</td>\n",
       "      <td>43</td>\n",
       "      <td>-1.941730e-03</td>\n",
       "      <td>1.361570e-03</td>\n",
       "      <td>-5.333800</td>\n",
       "      <td>856.0</td>\n",
       "      <td>-1152.5</td>\n",
       "      <td>-856</td>\n",
       "      <td>-781</td>\n",
       "      <td>1152</td>\n",
       "      <td>1797</td>\n",
       "      <td>-817.98</td>\n",
       "      <td>1464.43</td>\n",
       "      <td>34937</td>\n",
       "      <td>9.498680e-01</td>\n",
       "      <td>-818.3</td>\n",
       "      <td>1475.0</td>\n",
       "      <td>645.2</td>\n",
       "      <td>74.7</td>\n",
       "      <td>90.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>125</th>\n",
       "      <td>40</td>\n",
       "      <td>47</td>\n",
       "      <td>-7.749280e-11</td>\n",
       "      <td>2.229940e-09</td>\n",
       "      <td>-2.566560</td>\n",
       "      <td>1284.0</td>\n",
       "      <td>-1151.5</td>\n",
       "      <td>-1284</td>\n",
       "      <td>-781</td>\n",
       "      <td>1151</td>\n",
       "      <td>1205</td>\n",
       "      <td>-1032.04</td>\n",
       "      <td>1178.82</td>\n",
       "      <td>26883</td>\n",
       "      <td>7.354740e-06</td>\n",
       "      <td>-1032.5</td>\n",
       "      <td>1178.5</td>\n",
       "      <td>54.0</td>\n",
       "      <td>503.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>126</th>\n",
       "      <td>43</td>\n",
       "      <td>47</td>\n",
       "      <td>-8.581890e-04</td>\n",
       "      <td>-3.095250e-03</td>\n",
       "      <td>2.637200</td>\n",
       "      <td>855.0</td>\n",
       "      <td>-949.5</td>\n",
       "      <td>-855</td>\n",
       "      <td>-781</td>\n",
       "      <td>949</td>\n",
       "      <td>1205</td>\n",
       "      <td>-818.16</td>\n",
       "      <td>1078.88</td>\n",
       "      <td>14308</td>\n",
       "      <td>8.889450e-01</td>\n",
       "      <td>-818.0</td>\n",
       "      <td>1077.5</td>\n",
       "      <td>256.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>127</th>\n",
       "      <td>44</td>\n",
       "      <td>45</td>\n",
       "      <td>1.303560e-03</td>\n",
       "      <td>-1.999610e-03</td>\n",
       "      <td>-0.969453</td>\n",
       "      <td>-1718.0</td>\n",
       "      <td>-951.5</td>\n",
       "      <td>1718</td>\n",
       "      <td>1794</td>\n",
       "      <td>951</td>\n",
       "      <td>1206</td>\n",
       "      <td>1756.03</td>\n",
       "      <td>1080.04</td>\n",
       "      <td>14963</td>\n",
       "      <td>7.546190e-01</td>\n",
       "      <td>1756.0</td>\n",
       "      <td>1079.0</td>\n",
       "      <td>255.0</td>\n",
       "      <td>76.0</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>128 rows × 20 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     plus  minus             a             b          c  crpix1  crpix2  xmin  \\\n",
       "0       0      8  1.268790e-03  4.644410e-04   0.410733 -1711.0  1756.5  1711   \n",
       "1       0      9 -2.523660e-03  1.653820e-03   6.969540 -1713.0   988.5  1713   \n",
       "2       0     13  2.453210e-03  3.935060e-03   0.277418 -1286.0   986.5  1286   \n",
       "3       0     14 -2.490360e-03  8.991680e-04   4.436860 -1283.0  1755.5  1283   \n",
       "4       0     21 -3.303240e-10  5.376120e-09  -1.130130 -1286.0   786.5  1286   \n",
       "5       0     24 -1.489270e-10 -5.884170e-09   0.561656 -1284.0  1756.5  1284   \n",
       "6       0     29 -1.230960e-03  5.115690e-04 -79.055300 -1710.0   994.5  1710   \n",
       "7       0     31  2.329150e-03  1.100920e-03 -82.140500 -1708.0  1756.5  1708   \n",
       "8       1      6  1.756510e-04 -2.708390e-04  -3.833850   426.0   789.5  -426   \n",
       "9       1     11 -2.720230e-03 -5.054520e-04  -0.537130    -2.0   789.5     2   \n",
       "10      1     20 -5.109730e-19 -2.983720e-16   0.763428   425.0   789.5  -425   \n",
       "11      1     27  1.084200e-19 -5.074070e-17   0.691162   426.0  -180.5  -426   \n",
       "12      1     34 -4.281630e-03  6.650880e-04  -4.570030   426.0    20.5  -426   \n",
       "13      1     41 -9.965350e-03  3.869090e-03  -0.521526    -2.0    21.5     2   \n",
       "14      2      5  5.910760e-03  9.426410e-04   5.666030   854.0  1757.5  -854   \n",
       "15      2      6 -7.671010e-04  1.611060e-03   0.543132   854.0   990.5  -854   \n",
       "16      2     18  4.852270e-10  3.118610e-09   1.600290  1282.0   788.5 -1282   \n",
       "17      2     25 -3.896340e-04  9.108010e-04   6.162890  1281.0  1758.5 -1281   \n",
       "18      2     26  6.081620e-03  1.740740e-03  13.644000  1282.0   986.5 -1282   \n",
       "19      2     28 -3.959880e-20 -1.304840e-16  -0.597885  1280.0  1758.5 -1280   \n",
       "20      3     11  2.111740e-03  8.402230e-04  -2.689850  -431.0   784.5   431   \n",
       "21      3     13  3.627310e-03  1.314690e-03  -2.519670  -858.0   784.5   858   \n",
       "22      3     17  2.483370e-10 -1.152390e-09  -0.965677  -431.0   784.5   431   \n",
       "23      3     23 -3.036630e-11  1.690620e-09  -1.258960  -432.0  -185.5   432   \n",
       "24      3     38  3.505620e-03  1.811340e-04  -3.260770  -860.0    15.5   860   \n",
       "25      3     41  5.486890e-03  1.928290e-03  -4.866100  -431.0    21.5   431   \n",
       "26      4     15  4.008290e-03 -9.678930e-04   6.493410  1747.0  -183.5 -1747   \n",
       "27      4     18 -3.352750e-04 -2.270560e-03  -3.924300  1283.0    16.5 -1283   \n",
       "28      4     19  1.435180e-16 -1.116540e-14  -0.190674  1747.0  -952.5 -1747   \n",
       "29      4     26  2.168400e-18  4.445230e-18   1.108510  1745.0    17.5 -1745   \n",
       "..    ...    ...           ...           ...        ...     ...     ...   ...   \n",
       "98     25     37  2.993400e-03 -9.639940e-04   1.930760  1743.0  1756.5 -1743   \n",
       "99     26     37  2.132510e-03 -7.878530e-04   1.137900  1743.0   987.5 -1743   \n",
       "100    26     42  4.095660e-03 -1.870220e-03   4.653600  1745.0   786.5 -1745   \n",
       "101    27     33  2.527970e-18 -4.167670e-16   1.673900   426.0 -1150.5  -426   \n",
       "102    27     34  3.071140e-03  1.595790e-03  -2.654820   426.0  -180.5  -426   \n",
       "103    27     35  1.628030e-03  1.117010e-03  -0.391749    -3.0  -948.5     3   \n",
       "104    27     41  2.676220e-03  1.519470e-03  -1.326430    -2.0  -180.5     2   \n",
       "105    27     43  6.693140e-03  3.815910e-03  -3.391830   426.0  -949.5  -426   \n",
       "106    29     30 -6.003610e-03  7.217340e-03  92.591200 -1716.0    18.5  1716   \n",
       "107    29     31  3.228370e-09  3.008890e-09   2.299610 -1710.0   994.5  1710   \n",
       "108    30     44  1.178090e-15  1.691360e-17  -3.227810 -1718.0  -951.5  1718   \n",
       "109    30     45  5.933540e-03 -1.243340e-03 -13.183500 -1716.0  -183.5  1716   \n",
       "110    32     35  1.212370e-03 -2.488770e-03   2.879520  -433.0 -1155.5   433   \n",
       "111    32     36 -4.557880e-03  2.532880e-03   0.350670  -862.0 -1155.5   862   \n",
       "112    33     35  5.023290e-03  2.110940e-03  -3.150380    -3.0 -1150.5     3   \n",
       "113    33     43  3.743590e-03  1.023310e-03  -3.199060   426.0 -1150.5  -426   \n",
       "114    34     43  9.529120e-20 -1.970540e-17  -0.217102   855.0  -949.5  -855   \n",
       "115    34     47 -5.801320e-04 -1.266650e-03   1.172310   855.0  -182.5  -855   \n",
       "116    35     41  4.690440e-19  1.317850e-16  -0.366928    -3.0  -948.5     3   \n",
       "117    36     38  2.612920e-10 -1.019570e-09   1.721730  -861.0  -953.5   861   \n",
       "118    36     39 -5.956260e-04 -7.872360e-04  -0.036234 -1291.0 -1154.5  1291   \n",
       "119    36     45 -5.363510e-03 -4.439150e-03  11.325300 -1290.0  -953.5  1290   \n",
       "120    37     42  1.569380e-17  9.454240e-17   0.765747  1801.0   786.5 -1801   \n",
       "121    38     45 -3.689510e-03 -1.696570e-03   4.205390 -1289.0  -184.5  1289   \n",
       "122    39     44  2.348420e-03  4.838670e-04  -3.615190 -1718.0 -1153.5  1718   \n",
       "123    39     45  3.354750e-10  2.823240e-09   0.253280 -1291.0 -1153.5  1291   \n",
       "124    40     43 -1.941730e-03  1.361570e-03  -5.333800   856.0 -1152.5  -856   \n",
       "125    40     47 -7.749280e-11  2.229940e-09  -2.566560  1284.0 -1151.5 -1284   \n",
       "126    43     47 -8.581890e-04 -3.095250e-03   2.637200   855.0  -949.5  -855   \n",
       "127    44     45  1.303560e-03 -1.999610e-03  -0.969453 -1718.0  -951.5  1718   \n",
       "\n",
       "     xmax  ymin  ymax  xcenter  ycenter  npixel           rms    boxx    boxy  \\\n",
       "0    1789 -1756  -935  1750.15 -1343.07   52864  7.362550e-01  1749.9 -1346.0   \n",
       "1    1789  -988  -733  1751.01  -863.07   16232  7.125360e-01  1751.1  -861.0   \n",
       "2    1361  -986  -732  1323.43  -862.04   16218  7.174900e-01  1323.5  -859.5   \n",
       "3    1361 -1755  -933  1322.70 -1342.62   52649  7.105830e-01  1322.3 -1344.5   \n",
       "4    1790  -786  -732  1537.86  -759.37   27003  7.587320e-06  1538.0  -759.5   \n",
       "5    1787 -1756 -1702  1535.59 -1729.23   27090  7.108770e-06  1535.5 -1729.5   \n",
       "6    1790  -994  -733  1749.78  -865.65   16584  7.530540e-01  1749.9  -864.0   \n",
       "7    1789 -1756  -941  1748.71 -1348.03   53097  7.707710e-01  1748.8 -1349.0   \n",
       "8    -351  -789    32  -387.36  -347.45   42063  9.129390e-01  -388.5  -378.5   \n",
       "9      78  -789    31    40.09  -390.72   47921  8.727160e-01    40.0  -379.0   \n",
       "10     78  -789  -736  -172.39  -763.09   22413  4.641070e-15  -173.5  -763.0   \n",
       "11     77   180   233  -174.68   207.06   23318  7.921560e-16  -174.5   207.0   \n",
       "12   -351   -20   233  -388.07   106.65   14652  8.112380e-01  -388.5   106.5   \n",
       "13     77   -21   233    39.46   110.87   15314  9.693480e-01    39.5   106.0   \n",
       "14   -777 -1757  -937  -815.22 -1354.16   49825  7.697600e-01  -815.4 -1347.5   \n",
       "15   -778  -990  -734  -815.57  -864.23   14285  9.272300e-01  -816.0  -862.5   \n",
       "16   -778  -788  -734 -1029.37  -761.30   26871  7.403790e-06 -1030.0  -761.5   \n",
       "17  -1238 -1758  -933 -1260.25 -1335.75   24848  8.814390e-01 -1260.0 -1346.0   \n",
       "18  -1240  -986  -735 -1260.26  -864.96    7125  1.006160e+00 -1261.0  -861.0   \n",
       "19   -777 -1758 -1704 -1028.23 -1731.35   24404  2.032750e-15 -1028.5 -1731.5   \n",
       "20    505  -784    31   468.00  -378.34   50863  7.444600e-01   468.0  -376.5   \n",
       "21    935  -784    37   896.00  -382.34   50308  7.361090e-01   896.5  -373.5   \n",
       "22    934  -784  -731   682.61  -757.99   27199  7.623260e-06   682.5  -758.0   \n",
       "23    935   185   238   683.50   212.00   27216  6.063610e-06   683.5   212.0   \n",
       "24    935   -15   238   897.12   111.32   15496  7.787990e-01   897.5   111.5   \n",
       "25    505   -21   238   468.22   108.40   15892  7.351890e-01   468.0   108.5   \n",
       "26  -1672   183  1005 -1708.71   591.93   43190  9.386600e-01 -1709.7   594.6   \n",
       "27  -1242   -16   234 -1262.86   106.53    7062  1.012440e+00 -1262.5   109.0   \n",
       "28  -1244   974   983 -1493.10   978.90      47  3.698180e-16 -1495.5   979.5   \n",
       "29  -1242   -17    36 -1494.00     9.26   24954  3.432580e-16 -1493.5     9.5   \n",
       "..    ...   ...   ...      ...      ...     ...           ...     ...     ...   \n",
       "98  -1667 -1756  -934 -1704.44 -1368.87   43558  9.570520e-01 -1705.1 -1345.5   \n",
       "99  -1669  -987  -732 -1704.94  -856.26   13175  8.965520e-01 -1706.0  -860.0   \n",
       "100 -1670  -786    35 -1705.76  -365.01   44054  8.743800e-01 -1707.4  -375.5   \n",
       "101    77  1150  1203  -173.40  1177.04   21881  6.513750e-15  -174.5  1177.0   \n",
       "102  -351   180  1002  -389.66   570.28   44329  9.364340e-01  -388.5   591.5   \n",
       "103    77   948  1203    40.62  1076.31   12992  1.003060e+00    40.0  1076.0   \n",
       "104    77   180  1001    40.10   570.35   46199  8.684200e-01    39.5   591.0   \n",
       "105  -352   949  1203  -390.89  1078.05   12730  1.083480e+00  -389.0  1076.5   \n",
       "106  1798   -18    28  1757.75     4.26    3000  9.457370e-01  1757.0     5.0   \n",
       "107  1798  -994  -941  1754.05  -968.01    4791  7.628790e-06  1754.0  -968.0   \n",
       "108  1798   951  1004  1757.89   978.04    3747  2.755180e-14  1758.0   978.0   \n",
       "109  1793   183  1004  1754.93   607.41   46514  8.592740e-01  1754.9   594.0   \n",
       "110   506  1155  1797   468.94  1573.64   27217  9.886340e-01   469.5  1476.5   \n",
       "111   936  1155  1797   899.31  1465.39   37138  8.199690e-01   899.0  1476.5   \n",
       "112    77  1150  1797    40.41  1487.34   35991  8.733230e-01    40.0  1474.0   \n",
       "113  -352  1150  1797  -389.43  1485.53   33026  9.699490e-01  -389.0  1474.0   \n",
       "114  -352   949  1002  -603.82   976.05   24455  3.088130e-16  -603.5   976.0   \n",
       "115  -780   182  1002  -817.62   588.56   45873  8.761660e-01  -817.5   592.5   \n",
       "116   506   948  1001   254.91   974.98   26986  2.055550e-15   254.5   975.0   \n",
       "117  1364   953  1007  1111.91   980.63   26937  7.228180e-06  1112.5   980.5   \n",
       "118  1366  1154  1797  1327.21  1483.97   36016  9.384890e-01  1328.3  1476.0   \n",
       "119  1365   953  1207  1327.66  1080.66   13662  9.161630e-01  1327.6  1080.5   \n",
       "120 -1670  -786  -732 -1735.80  -759.64    5746  1.590380e-15 -1735.5  -759.5   \n",
       "121  1364   184  1006  1326.44   596.48   48646  8.512290e-01  1326.5   595.5   \n",
       "122  1795  1153  1797  1756.99  1475.15   39354  7.699550e-01  1756.8  1475.5   \n",
       "123  1794  1153  1207  1542.22  1180.83   27010  6.472080e-06  1542.5  1180.5   \n",
       "124  -781  1152  1797  -817.98  1464.43   34937  9.498680e-01  -818.3  1475.0   \n",
       "125  -781  1151  1205 -1032.04  1178.82   26883  7.354740e-06 -1032.5  1178.5   \n",
       "126  -781   949  1205  -818.16  1078.88   14308  8.889450e-01  -818.0  1077.5   \n",
       "127  1794   951  1206  1756.03  1080.04   14963  7.546190e-01  1756.0  1079.0   \n",
       "\n",
       "     boxwidth  boxheight  boxang  \n",
       "0       821.2       77.3   -90.1  \n",
       "1       255.3       76.9   -90.2  \n",
       "2       254.0       75.0    90.0  \n",
       "3       822.2       77.3    89.9  \n",
       "4        54.0      504.0    90.0  \n",
       "5        54.0      503.0    90.0  \n",
       "6       261.3       79.3   -90.2  \n",
       "7       815.2       79.9    89.8  \n",
       "8       822.0       75.0    90.0  \n",
       "9       821.0       76.0    90.0  \n",
       "10       53.0      503.0    90.0  \n",
       "11       53.0      503.0    90.0  \n",
       "12      254.0       75.0    90.0  \n",
       "13      255.0       75.0    90.0  \n",
       "14      820.1       76.4    90.1  \n",
       "15      256.0       76.0    90.0  \n",
       "16       54.0      504.0    90.0  \n",
       "17      825.1       42.7   -89.9  \n",
       "18      251.0       42.0    90.0  \n",
       "19       54.0      503.0    90.0  \n",
       "20      816.0       74.0    90.0  \n",
       "21      822.1       76.1   -90.1  \n",
       "22       53.0      503.0    90.0  \n",
       "23       53.0      503.0    90.0  \n",
       "24      254.0       75.0    90.0  \n",
       "25      260.0       74.0    90.0  \n",
       "26      822.0       73.7   -89.9  \n",
       "27      251.0       41.0    90.0  \n",
       "28       54.0      503.0    90.0  \n",
       "29       54.0      503.0    90.0  \n",
       "..        ...        ...     ...  \n",
       "98      822.1       75.0    90.1  \n",
       "99      255.0       74.0    90.0  \n",
       "100     822.0       74.3    90.1  \n",
       "101      53.0      503.0    90.0  \n",
       "102     822.0       75.0    90.0  \n",
       "103     255.0       74.0    90.0  \n",
       "104     821.0       75.0    90.0  \n",
       "105     254.0       74.0    90.0  \n",
       "106      47.0       82.0    90.0  \n",
       "107      53.0       88.0    90.0  \n",
       "108      53.0       80.0    90.0  \n",
       "109     821.2       76.3   -90.2  \n",
       "110     642.0       73.0    90.0  \n",
       "111     642.0       74.0    90.0  \n",
       "112     647.0       74.0    90.0  \n",
       "113     647.0       74.0    90.0  \n",
       "114      53.0      503.0    90.0  \n",
       "115     820.0       75.0    90.0  \n",
       "116      53.0      503.0    90.0  \n",
       "117      54.0      503.0    90.0  \n",
       "118     643.1       74.5   -90.1  \n",
       "119     254.4       74.4    89.7  \n",
       "120      54.0      131.0    90.0  \n",
       "121     822.0       75.0    90.0  \n",
       "122     644.3       76.8   -90.2  \n",
       "123      54.0      503.0    90.0  \n",
       "124     645.2       74.7    90.1  \n",
       "125      54.0      503.0    90.0  \n",
       "126     256.0       74.0    90.0  \n",
       "127     255.0       76.0    90.0  \n",
       "\n",
       "[128 rows x 20 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import pandas as pd\n",
    "from astropy.io import ascii\n",
    "\n",
    "ipactable = ascii.read('work/M17/fits.tbl').to_pandas()\n",
    "\n",
    "ipactable"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "&nbsp;<p/>\n",
    "\n",
    "## Error Handling\n",
    "\n",
    "If mDiffFitExec encounters an error, the return structure will just have two elements: a status of 1 (\"error\") and a message string that tries to diagnose the reason for the error.\n",
    "\n",
    "For instance, if the user specifies a table that doesn't exist: \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'status': '1', 'msg': b'Invalid diffs metadata file: M17/unknown.tbl'}\n"
     ]
    }
   ],
   "source": [
    "rtn = mDiffFitExec('M17/projected', \n",
    "                   'M17/unknown.tbl', \n",
    "                   'M17/M17.hdr',\n",
    "                   'M17/diffs',\n",
    "                   'work/M17/fits.tbl')\n",
    "print(rtn)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "&nbsp;\n",
    "\n",
    "## Classic Montage: mDiffFitExec as a Stand-Alone Program\n",
    "\n",
    "\n",
    "### mDiffFitExec Unix/Windows Command-line Arguments\n",
    "\n",
    "<p>mDiffFitExec can also be run as a command-line tool in Linux, OS X, and Windows:</p>\n",
    "\n",
    "<p><tt>\n",
    "<b>Usage:</b> mDiffFitExec [-d] [-l(evel-only)] [-n(o-areas)] [-p projdir] [-s statusfile] diffs.tbl template.hdr diffdir fits.tbl\n",
    "</tt></p>\n",
    "<p>&nbsp;</p>\n",
    "<p>If you are writing in C/C++, mDiffFitExec can be accessed as a library function:</p>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<pre>\n",
    "/*-*****************************************************************/\n",
    "/*                                                                 */\n",
    "/*  mDiffFitExec                                                   */\n",
    "/*                                                                 */\n",
    "/*  This routine combines the mDiff and mFit functionality and     */\n",
    "/*  optionally discards the difference images as it goes (to       */\n",
    "/*  minimize the file space needed).   It uses the table of        */\n",
    "/*  oerlaps found by mOverlaps, running mDiff, then mFitplane      */\n",
    "/*  on the difference images.  These fits are written to an        */\n",
    "/*  output file which is then used by mBgModel.                    */\n",
    "/*                                                                 */\n",
    "/*   char *path        Path to images to be diffed.                */\n",
    "/*   char *tblfile     Table file list of images to diff.          */\n",
    "/*   char *template    FITS header file used to define the desired */\n",
    "/*                     output.                                     */\n",
    "/*                                                                 */\n",
    "/*   char *diffdir     Directory for temporary output diff files.  */\n",
    "/*   char *fitfile     Table file for output difference fits info. */\n",
    "/*   int keepAll       Flag to keep temporary diff images.         */\n",
    "/*   int levelOnly     Flag to fit level of diff only, not slopes. */\n",
    "/*   int noAreas       Flag indicating there are no area images.   */             \n",
    "/*   int debug         Debug flag.                                 */\n",
    "/*                                                                 */\n",
    "/*******************************************************************/\n",
    "\n",
    "struct mDiffFitExecReturn *mDiffFitExec(char *path, char *tblfile, char *template, char *diffdir,\n",
    "                                        char *fitfile, int keepAll, int levelOnly, int noAreas, int debugin)\n",
    "</pre>\n",
    "<p><b>Return Structure</b></p>\n",
    "<pre>\n",
    "struct mDiffFitExecReturn\n",
    "{\n",
    "   int    status;        // Return status (0: OK, 1:ERROR)\n",
    "   char   msg [1024];    // Return message (for error return)\n",
    "   char   json[4096];    // Return parameters as JSON string\n",
    "   int    count;         // Number of differences                \n",
    "   int    diff_failed;   // Number of differences that failed\n",
    "   int    fit_failed;    // Number of fits to differences that failed\n",
    "   int    warning;       // Number of fits to differences that produced warnings\n",
    "};\n",
    "</pre>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}