Analyze_teakreg_shifts.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "bb87181e",
   "metadata": {},
   "source": [
    "### Analyzing tweakreg shifts\n",
    "\n",
    "Author: Karl Gordon\n",
    "\n",
    "5 Jul 2022"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "1503b24f",
   "metadata": {},
   "outputs": [],
   "source": [
    "from tweakwcs import JWSTgWCS\n",
    "from jwst import datamodels\n",
    "from astropy.stats import sigma_clip\n",
    "import numpy as np\n",
    "import glob"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "3ff957e3",
   "metadata": {},
   "outputs": [],
   "source": [
    "n_tiles = 5\n",
    "n_images_per_tile = 10"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "92f030fa",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "./F560W/stage3/jw01040001001_03101_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.1046129  -0.10911081] S1HP090792\n",
      "./F560W/stage3/jw01040001001_03101_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.11033027 -0.11822098] S1HP090792\n",
      "./F560W/stage3/jw01040001001_03101_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.09954632 -0.11549697] S1HP090792\n",
      "./F560W/stage3/jw01040001001_03101_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.10440985 -0.11558269] S1HP090792\n",
      "./F560W/stage3/jw01040001001_03101_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.10531285 -0.11122491] S1HP090792\n",
      "./F560W/stage3/jw01040001001_0310j_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.10934406 -0.11607055] S1HP090792\n",
      "./F560W/stage3/jw01040001001_0310j_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.10699728 -0.12492857] S1HP090792\n",
      "./F560W/stage3/jw01040001001_0310j_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.1052372  -0.11965154] S1HP090792\n",
      "./F560W/stage3/jw01040001001_0310j_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.10873223 -0.12227168] S1HP090792\n",
      "./F560W/stage3/jw01040001001_0310j_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.10924137 -0.12426363] S1HP090792\n",
      "./F560W/stage3/jw01040001002_03101_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.11225891 -0.13882681] S1HP091481\n",
      "./F560W/stage3/jw01040001002_03101_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.11495412 -0.14045409] S1HP091481\n",
      "./F560W/stage3/jw01040001002_03101_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.10680332 -0.13463987] S1HP091481\n",
      "./F560W/stage3/jw01040001002_03101_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.10664813 -0.13479887] S1HP091481\n",
      "./F560W/stage3/jw01040001002_03101_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.11546861 -0.13790822] S1HP091481\n",
      "./F560W/stage3/jw01040001002_0310j_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.10572435 -0.11536618] S1HP091481\n",
      "./F560W/stage3/jw01040001002_0310j_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.10292854 -0.12656802] S1HP091481\n",
      "./F560W/stage3/jw01040001002_0310j_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.10169142 -0.12285754] S1HP091481\n",
      "./F560W/stage3/jw01040001002_0310j_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.10146679 -0.13218628] S1HP091481\n",
      "./F560W/stage3/jw01040001002_0310j_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.10212156 -0.12719029] S1HP091481\n",
      "./F560W/stage3/jw01040001003_03101_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.1060594  -0.10748639] S1HP093172\n",
      "./F560W/stage3/jw01040001003_03101_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.10350294 -0.12131731] S1HP093172\n",
      "./F560W/stage3/jw01040001003_03101_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.10174589 -0.12678009] S1HP093172\n",
      "./F560W/stage3/jw01040001003_03101_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.10026455 -0.12175042] S1HP093172\n",
      "./F560W/stage3/jw01040001003_03101_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.10366742 -0.12346722] S1HP093172\n",
      "./F560W/stage3/jw01040001003_0310j_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.09540509 -0.14431798] S1HP093172\n",
      "./F560W/stage3/jw01040001003_0310j_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.10318174 -0.11834792] S1HP093172\n",
      "./F560W/stage3/jw01040001003_0310j_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.10052989 -0.11568961] S1HP093172\n",
      "./F560W/stage3/jw01040001003_0310j_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.09777084 -0.11814375] S1HP093172\n",
      "./F560W/stage3/jw01040001003_0310j_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.10426461 -0.11791517] S1HP093172\n",
      "./F560W/stage3/jw01040001004_03101_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.08936902 -0.1120277 ] S1HP093172\n",
      "./F560W/stage3/jw01040001004_03101_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.08753498 -0.12052266] S1HP093172\n",
      "./F560W/stage3/jw01040001004_03101_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.08976265 -0.11286233] S1HP093172\n",
      "./F560W/stage3/jw01040001004_03101_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.08844734 -0.11618403] S1HP093172\n",
      "./F560W/stage3/jw01040001004_03101_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.09122827 -0.12278127] S1HP093172\n",
      "./F560W/stage3/jw01040001004_0310j_00001_mirimage_fixed_wcs_tweakreg.fits [ 0.09437473 -0.11548969] S1HP093172\n",
      "./F560W/stage3/jw01040001004_0310j_00002_mirimage_fixed_wcs_tweakreg.fits [ 0.08928152 -0.12696451] S1HP093172\n",
      "./F560W/stage3/jw01040001004_0310j_00003_mirimage_fixed_wcs_tweakreg.fits [ 0.08700749 -0.12716074] S1HP093172\n",
      "./F560W/stage3/jw01040001004_0310j_00004_mirimage_fixed_wcs_tweakreg.fits [ 0.09409717 -0.13180387] S1HP093172\n",
      "./F560W/stage3/jw01040001004_0310j_00005_mirimage_fixed_wcs_tweakreg.fits [ 0.09543207 -0.12869362] S1HP093172\n",
      "./F560W/stage3/jw01040001005_03101_00001_mirimage_fixed_wcs_tweakreg.fits [-0.97318695  1.03764066] S1HP093000\n",
      "./F560W/stage3/jw01040001005_03101_00002_mirimage_fixed_wcs_tweakreg.fits [-0.97524167  1.04695409] S1HP093000\n",
      "./F560W/stage3/jw01040001005_03101_00003_mirimage_fixed_wcs_tweakreg.fits [-0.97705709  1.04906997] S1HP093000\n",
      "./F560W/stage3/jw01040001005_03101_00004_mirimage_fixed_wcs_tweakreg.fits [-0.97934219  1.04616322] S1HP093000\n",
      "./F560W/stage3/jw01040001005_03101_00005_mirimage_fixed_wcs_tweakreg.fits [-0.97394367  1.04842164] S1HP093000\n"
     ]
    }
   ],
   "source": [
    "RAD2ARCSEC = 3600.0 * np.rad2deg(1.0)\n",
    "\n",
    "filter = \"F560W\"\n",
    "twfiles = np.sort(glob.glob(f\"./{filter}/stage3/*wcs_tweakreg.fits\"))\n",
    "\n",
    "shifts = np.zeros((2, len(twfiles)))\n",
    "for k, cfile in enumerate(twfiles):\n",
    "    # extract x,y shifts from the aligned image:\n",
    "    aligned_model = datamodels.open(cfile)\n",
    "    matrix = aligned_model.meta.wcs.forward_transform['tp_affine'].matrix.value\n",
    "    cshift = RAD2ARCSEC * aligned_model.meta.wcs.forward_transform['tp_affine'].translation.value\n",
    "    shifts[:, k] = cshift\n",
    "    print(cfile, cshift, aligned_model.meta.guidestar.gs_id)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "5abffd5d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7fb7bf5cd960>]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAjIElEQVR4nO3dfXAc933f8fcXB+AODwRxICA+gU+SKFGyFVMuLKV246a2ZMlWLEqNE8tpatljj2Y6UZsmzYNUd+zUiTt2M1M7k3ETs7Zi2XEsPySuWYsZVk9pOlVkE3JkPVOk+CDx+UjgABC4A3DAt3/sHnB4xvEOOOD285r5zW93b/fuhwWwn9vdu9/P3B0REYmumko3QEREKktBICIScQoCEZGIUxCIiEScgkBEJOJqK92Ay9He3u7bt2+vdDNERFaVZ5999oK7d0xfviqDYPv27XR3d1e6GSIiq4qZnZhtuS4NiYhEnIJARCTiFAQiIhGnIBARiTgFgYhIxCkIREQiTkEgIhJxq/J7BCIiK1U2l+V4+jjH08e5NHKJbC5LNpdlODcc1GPDjIyNEI/FaapvorGukaa6sK5vIlGbYGRsZGK76eXjuz9OsiFZ1jYrCERELtOBIwd4+s2nOZY+xtHeoxxLH+P0wOklfc33X/1+BYGIyEowMjbCHX91B+M+TmdLJ1cmr+S2q25jR+sOrkxeyfbW7bTEW0jUJiZKvDZOojZBXU0dw2PDDI0OMTgyGNSjQZ0ZzUysV1gaahtI1CZoqm8q+89SliAws4eAXwLOu/tbZ3ncgD8BPgAMAR9z95+Gj90L/Kdw1T9y94fL0SYRkaWUzqYZ8zH+9P1/yv033V/09vkDfFtD2xK0rjjluln8deD2eR5/P7AzLPcBfwZgZm3AZ4CbgZuAz5hZec95RESWQDqbBqA10VrRdpRDWYLA3f8e6JlnlT3ANzzwDNBqZhuB24DH3L3H3XuBx5g/UEREVoS+bB+gICjGZuDNgvmT4bK5ls9gZveZWbeZdadSqSVrqIjIYuiMoALcfa+7d7l7V0fHjO60RUSWVT4I1sbXVrYhZbBcQXAK2FIw3xkum2u5iMiK1jesS0PF2gd81AI/D/S5+xngAPA+M0uGN4nfFy4TEVnRJs4IEqv/jKBcHx/9NvCLQLuZnST4JFAdgLv/ObCf4KOjRwg+Pvrx8LEeM/tD4GD4VJ919/luOouIrAh92T5qrIbm+uZKN6VkZQkCd//IAo878BtzPPYQ8FA52iEislzS2TQt8RZqbNXcap3T6v8JREQqoG+4ryruD4CCQETksqSzaQWBiEiUpbPpqvjoKCgIREQuiy4NiYhEXDqbroqPjoKCQETksvRl+2iNt1a6GWWhIBARKdK4j9M/3K8zAhGRqBoYHsBx3SMQEYmqaup5FBQEIiJFq6aeR0FBICJStGrqeRQUBCIiRaumnkdBQSAiUrRqGqYSFAQiIkXTPQIRkYjL3yPQpaECZna7mR0ysyNm9sAsj3/RzJ4Ly2tmli54bKzgsX3laI+IyFJKZ9M01jVSH6uvdFPKouSBacwsBnwZuBU4CRw0s33u/nJ+HXf/rYL1/y1wY8FTZNx9d6ntEBFZLtXU8yiU54zgJuCIux919xHgEWDPPOt/BPh2GV5XRKQiqqnnUShPEGwG3iyYPxkum8HMtgE7gCcLFifMrNvMnjGzu+Z6ETO7L1yvO5VKlaHZIiKXp5p6HoXlv1l8D/B9dx8rWLbN3buAXwO+ZGZXzbahu+919y537+ro6FiOtoqIzKovqzOC6U4BWwrmO8Nls7mHaZeF3P1UWB8F/o6p9w9ERFacahqmEsoTBAeBnWa2w8zqCQ72Mz79Y2a7gCTwDwXLkmYWD6fbgXcBL0/fVkRkJam2m8Ulf2rI3XNmdj9wAIgBD7n7S2b2WaDb3fOhcA/wiLt7webXAV8xs3GCUPp84aeNRERWGnevupvFJQcBgLvvB/ZPW/bpafN/MMt2TwM3lKMNIiLLIZvLMjI2UlVnBPpmsYhIEaqt51FQEIiIFKXaeh4FBYGISFGqredRUBCIiBSl2oapBAWBiEhRqq0LalAQiIgURTeLRUQiTjeLRUQiri/bR8xiNNU1VbopZaMgEBEpQr7nUTOrdFPKRkEgIlKEauteAhQEIiJFqbaeR0FBICJSlGrreRQUBCIiRdGlIRGRiKu2YSpBQSAiUpS+bB+t8dZKN6OsFAQiIos0Nj7GwMiAzghmY2a3m9khMztiZg/M8vjHzCxlZs+F5ZMFj91rZofDcm852iMishT6h/uB6upeAsowQpmZxYAvA7cCJ4GDZrZvliEnv+Pu90/btg34DNAFOPBsuG1vqe0SESm3aux5FMpzRnATcMTdj7r7CPAIsGeR294GPObuPeHB/zHg9jK0SUSk7Kqx51EoTxBsBt4smD8ZLpvul83seTP7vpltKXJbzOw+M+s2s+5UKlWGZouIFKcaex6F5btZ/L+A7e7+cwTv+h8u9gncfa+7d7l7V0dHR9kbKCKykGrseRTKEwSngC0F853hsgnuftHdh8PZrwL/ZLHbioisFNU4TCWUJwgOAjvNbIeZ1QP3APsKVzCzjQWzdwKvhNMHgPeZWdLMksD7wmUiIitOtd4sLvlTQ+6eM7P7CQ7gMeAhd3/JzD4LdLv7PuDfmdmdQA7oAT4WbttjZn9IECYAn3X3nlLbJCKyFPJB0BJvqWxDyqzkIABw9/3A/mnLPl0w/SDw4BzbPgQ8VI52iIgspb7hPprrm6mtKcuhc8XQN4tFRBapGnseBQWBiMiiVWPPo6AgEBFZtGrseRQUBCIii9aX1RmBiEikVeMwlaAgEBFZNN0sFhGJMHfXzWIRkSgbGh0iN57TGYGISFRVa8+joCAQEVmUau15FBQEIiKLUq09j4KCQERkUaq151FQEIiILEq1DlMJCgIRkUXRzWIRkYjTzeIFmNntZnbIzI6Y2QOzPP7bZvZyOHj9E2a2reCxMTN7Liz7pm8rIrIS9GX7qKupo6G2odJNKbuSR1cwsxjwZeBW4CRw0Mz2ufvLBav9I9Dl7kNm9m+A/wp8OHws4+67S22HiMhSyvczZGaVbkrZleOM4CbgiLsfdfcR4BFgT+EK7v6Uuw+Fs88QDFIvIrJq9A33VeVlIShPEGwG3iyYPxkum8sngL8tmE+YWbeZPWNmd821kZndF67XnUqlSmqwiEixqrXnUSjTmMWLZWa/DnQB/7xg8TZ3P2VmVwJPmtkL7v769G3dfS+wF6Crq8uXpcEiIqFq7XkUynNGcArYUjDfGS6bwsxuAT4F3Onuw/nl7n4qrI8CfwfcWIY2iYiUVbX2PArlCYKDwE4z22Fm9cA9wJRP/5jZjcBXCELgfMHypJnFw+l24F1A4U1mEZEVoZrPCEq+NOTuOTO7HzgAxICH3P0lM/ss0O3u+4A/BpqB74V33N9w9zuB64CvmNk4QSh9ftqnjUREVoRqHaYSynSPwN33A/unLft0wfQtc2z3NHBDOdogIrJURsdGGRwdrNog0DeLRUQWkO9eQh8fFRGJqGrughoUBCIiC6rmnkdBQSAisqBq7nkUFAQiIguq5p5HQUEgIrKgar9HsKxdTKxk4+NgFpRqlsvB4OBkGRoK6kwmmB4amjpdWwtbtgRl61a44gqoKeLtw9gYnD4Nx4/DiRNBnUjAVVdNlqampfppRcqjmoephIgFwfe+B88/D6kUXLgwtb54ERob4ZprgnLttZPlmmugvh56e4PS0zO1Hh4ODrBjY0HJT+dycOkS9PdDX9/UemAAYjFoaAhKIjG1rq+HeDyoC6dra4MD98BA8NyF9dBQ8Lrj4+Ae1PmSD4CRkdL2YV0ddHYGoZBMTg3PwrqnJzjov/FG8Nrz2bBhMhTa24Pfw2wl//PX1c2s8+0obEO+FO6H/P7J75P+/pmlry/Yl/F48LtobJz8PeXLbL+bfF1XF5Tp0+6Tv6/C392lS8E6W7cGZf36ucPWPdi3J08GJZeb2sbCOj8di839XP39cPYsnDsX1BcvBr+Hd7wj+P1KIB8Ea+rXVLYhSyRSQfDNb8KPfgTr1kFHR1B27YJf+IXgADQwAIcOwTPPwHe+E/yjXA6z4J8vFoPmZli7Flpagnrr1qBubg4OStls8A68sE6ngwP2yEgQMoXTuVzwDnrNmuA51qwJyqZNwT9+bW3w+jU1U0ssFmw3W5l+wC08oIyMwJtvBuWNNybrN96A11+f3EeFtTu0tsLNN8OHPwzbtwdl27bg589mg20Ly5Ej8OSTQbAODV3+vi+VWbA/GxqC/Z3JBPVyKgzbrVuDZfkD/8mTQZuKEY/P/P329QUH/vl+tquvhptumiy7dwfbFsM9+DuPxRY+2x4dnTwTzZ+pDg4GQTl92h3a2oL/5ba2yelkMgjdhYyOBuF3+jScORPUly4Fy6eX8XG40NVHS7yFWM0cqbrKmVfqP64EXV1d3t3dXfR2g4PBO+653iEVymSCg9OhQ/Daa8EfQ1tb8Ic2vc4/Z23t4v7gZX7uwQGq8KAwNDT5j5nLTa1HRycDKL99YYnFZoZivm5pmVqam2e+Gx8fD8Irf9kskwlec3pI5+vCg8jIyOQ0TAZ4vuTns9mpIVs4DcGluc7OmSUen3k5L5OZebmvsGQywZuR9euDsmHD5HQyGfzN/+QnQTl4EE6FXUiaBX/r+bOfwjOh2trgZ81mgzI8PDmd/73EYpNnSIVncvl2LnTmuFj5M+vppaEheJ3Tp4OrAPMd+swm2zo4CNf83sfIbHiSN37rjfI0skLM7Fl375q+PFJnBMVci25ogBtuCIosr/wBJ5EIwrbSamom300vpbe9bWmff7E6O+G9752cP3UqCITnngsOioWhl58eHQ0CYbYDcF1dcFYwPcjz77YLz1amn6E2NQVBmT97zU9DcImspye4nFVY9/dPhlC+5AM8mQwue23aBBs3TtYbNwbhmA+pwjeL114Lg7k+2qr0/gBELAhEpHibNwflrrsq3ZKp1q6FHTuW/nWSSUh5mh1V+tFR0MdHRUTmlUzCMNXb8ygoCERE5pVMwkhN9Q5TCQoCEZF5JZMwVle9g9JAmYLAzG43s0NmdsTMHpjl8biZfSd8/Mdmtr3gsQfD5YfM7LZytEdEpFxak47X97E23lrppiyZkoPAzGLAl4H3A9cDHzGz66et9gmg192vBr4IfCHc9nqCoS3fAtwO/Pfw+UREVoTG1ktQM04CnRHM5ybgiLsfdfcR4BFgz7R19gAPh9PfB95rwZiVe4BH3H3Y3Y8BR8LnExFZEeItQT9DtWOtlW3IEipHEGwG3iyYPxkum3Udd88BfcC6RW4LgJndZ2bdZtadSqXK0GwRkYXVrkkDEBttrWg7ltKquVns7nvdvcvduzo6OirdHBGJiFhjcEZAVpeG5nMK2FIw3xkum3UdM6sF1gIXF7mtiEjlJNIAjA+1VrQZS6kcQXAQ2GlmO8ysnuDm775p6+wD7g2nPwQ86UEnR/uAe8JPFe0AdgI/KUObRETKYrw+DUDuUvWeEZTcxYS758zsfuAAEAMecveXzOyzQLe77wO+BnzTzI4APQRhQbjed4GXgRzwG+4+VmqbRETKZaw2uDQ0MtBa2YYsobL0NeTu+4H905Z9umA6C/zKHNt+DvhcOdohIlJuGU8DkE1X7xnBqrlZLCJSCf3DfZCLM9CbqHRTloyCQERkHulsmthoK729lW7J0lEQiIjMIz2cpm5srYJARCSq+rJ9xL26zwgiNTDNH/+/P+b13tfZ0bqDHckdE/W6hnXYPONLjvs44z4+ZZkxuf7g6CDpbHqi9GX7SGfT9A/301DXQDKRpK2hjWRDWCeSNNY1zvuaUZQfNrVS+2VgeICnjj/FgSMHeDH1IrvW7eLtG9/OjRtv5IYrbqChrsgBe1e5kbERUoMpzg+eJzWUIp1N09HYwbbWbXS2dFIfm3twYHenJ9PDyf6TXBi6QENdA831zVNKQ23DZf2ul/vvJJ1Nk6ip7jOCSAXBoYuH+MGrP6An0zNleXN9Mztad9BU30RmNMPQ6NBEyeQyZHPZsrclHovT2dLJ9tbtU8q2tdvY3LKZwZFBLmYucnHo4pS6L9tHvDZOQ20DjXWNNNY10lAXTDfUNhCvjROPxWfUrYlWtrduX3HhMzA8wP7D+/nBqz/g0cOPUldTx672Xexq38V17dcFdcd1bG/dTsxiZHNZBkYG6B/uZ2B4gIGRAQZHBlnfvJ4rk1cW1Wf8uI/z3NnnOHDkAAdeP8DTbz7N6PgoTXVNvPWKt/K9l7/H3p/uBSBmMa7vuJ63b3w7b73iraxvWk97Yzvtje2sa1xHe2M7a+rXYGbkxnP0ZnrpyfTQk+nhYuYiPZkesrks9bH6KSUei1Mfq2fL2i1c3Xb1otqdG8/x+NHH+dYL3+JY7zEyucm/2cK/31hNjIbaBhrqGmbU+TcyjuPuE/W4j9OT6SE1lKJ/uH/ONhjGxjUb2bZ2G9tat7G+aT2poRSn+k9xsv8kpwZOLfh/Yxgt8RbWN69nQ/OGoDRtmJhvqG3gzKUznB44zamBU5weOD1Rxn2ca9ddy/Ud13Nd+3Vc13Ed13dcz9VtV88bUJejb7iP5tjWqg6CSA1en9c/3M/x9HGO9R7jWPrYRJ3NZScPrtMOtLU1k5lZuM8cp6muidZE60RZm1hLa6KVlngLmdEMvdngoJA/OPRme7kwdIE3+9/kePo4J9InOHPpzILtbqhtoCXewuj4KEOjQ0UHVGuila5NXbxj0zsm6s6WzqLDIZvLkhpMMTQ6BEy+M8sfXMyM+lg9TXVNNNU3EY/Fp7xGajDFvkP7+JtX/4bHjz7OyNgIHY0d3HntndTH6nnlwiu8euFVzl46O7FNbU0t7s7YAl8zaU20cmXySna07uDK5JVsadlCJpeZOCgXlpP9J7mYuQjA7g27ue2q27jtqtt455Z3Eq+N4+6c6DvBT8/8dEo5N3hu1teuq6kjUZtgYGSgqP2Zt7NtJ3fsvIM7rrmDd29794wD2s/O/oxv/Owb/NWLf8XZS2dpa2hj94bdU/5eJ+q6BsbGx8jkMmRGM2THsmRGMxPzjk/5fRmGmVFjNSQTSToaO+ho6qCjsYMrmq6go6mD1kQr5y6d40TfCd7oe4MTfSc4kT7Bib4TnB88T0djB50tnXS2dLJ5zeaJ6fbGdobHhrk0cmlGSWfTnBs8x9lLZyfK9ABK1CbYvGYzm9ZsYnPLZjY1bwLglQuv8MqFVziePj6xbm1NLe2N7RM/z/SfMVGbINmQJJlI0ppoJZlIkmwIpg3j0sglBkcHGRwZDOrRQfYf3s+u4V/nhf/yPxgdDcbUXq3mGrw+kkGwEmVzWd7oe4Pj6eOc6j9Fc30z6xrXsa5h3UQ9/dLEuI+TzWUn3gVmchmGc8MMjw0znBsmm8tOTKeGUnSf7qb7dDcvnH+B3HgOgCuaruC69utI1CZmvFutj9UzOj5KajBFaig1UV8auVTUz1ZjNTTXN9NU10RDXQPH08cZ93G2t27n7l13c/euu3nnlncSq5naA3k6m+bVC6/y6oVXee3iawC0xFtYU7+GNfE1E9ONdY2cuXSGY73HONp7lGPpyXpkbAQIDtLrGtfR1tA2Ua5ovIJ3b3s3t151KxuaNyz650ln01wYujCjXBy6SCaXmfIahaWhtoHR8VFGxkYmSv739eL5F3n08KM8dewphseGWVO/hluvupU7dt5Bb6aXbzz/DZ4/9zx1NXXccc0dfPTnPsoHdn6AeG28qN/FajA0OsS5S+fI5DJsWrOJtfG1875ZGRwZ5NDFQ7ySeoWXUy+TGkpNnOUAk2c8+MQbs3Q2TW+md2I6//8A0FjXOPEmpqmuieb6Zt5y4dM89OAH6OuDlpYl3wVLRkEgEzKjGZ4/9zwHTx+k+3Q3r/e+PuXgNDo2ebCqsRo6msJ3hY0dU94pNtU3TTxn/u8of4lhZGyEwdHB4B1W+O4q/25rZ9tO7t51N7s37F7SS1XjPs6FoQs01TWtmnsygyODPHnsSX702o949PCjnBoIut66efPNfPRtH+XDb/kw6xrXVbiV1cXdGRwdBIIQqLGZn6H52tfgk5+E48dh27ZlbmAZzRUEkbpHIIGGugZu7ryZmztvrnRTllSN1XBF0xWVbkZRmuqb+OC1H+SD134Qd+fF8y+SqE2wc93OSjetapkZzfXN866TTAZ1b+/qDoK5KAhEVigz44b1N1S6GcLUIKhG+h6BiMgCFAQiIhGnIBARiTgFgYhIxK1ZA7GYgkBEJLLMoLVVQTArM2szs8fM7HBYJ2dZZ7eZ/YOZvWRmz5vZhwse+7qZHTOz58Kyu5T2iIgslWRSQTCXB4An3H0n8EQ4P90Q8FF3fwtwO/AlM2stePx33X13WJ4rsT0iIktCQTC3PcDD4fTDwF3TV3D319z9cDh9GjgPdJT4uiIiy0pBMLf17p7vLe0ssH6+lc3sJqAeeL1g8efCS0ZfNLM5O04xs/vMrNvMulOpVInNFhEpTqSDwMweN7MXZyl7CtfzoLOZOTsuMrONwDeBj7tPdO7/ILALeAfQBvz+XNu7+15373L3ro4OnVCIyPKq5iBYsIsJd79lrsfM7JyZbXT3M+GB/vwc67UAjwKfcvdnCp47fzYxbGZ/AfxOUa0XEVkm+SBwX91dUc+m1EtD+4B7w+l7gR9OX8HM6oEfAN9w9+9Pe2xjWBvB/YUXS2yPiMiSSCZhbAwuFdcL+6pQahB8HrjVzA4Dt4TzmFmXmX01XOdXgXcDH5vlY6LfMrMXgBeAduCPSmyPiMiSqOZvF5fU+6i7XwTeO8vybuCT4fRfAn85x/bvKeX1RUSWS2EQbN1a2baUm75ZLCKyCNV8RqAgEBFZBAWBiEjEKQhERCJOQSAiEnFr1kBNjYJARCSyamqqtytqBYGIyCJVazcTCgIRkUVSEIiIRJyCQEQk4hQEIiIRpyAQEYm4wq6oq4mCQERkkZJJyOVgcLDSLSkvBYGIyCJV67eLFQQiIoukIJiFmbWZ2WNmdjisk3OsN1YwKM2+guU7zOzHZnbEzL4TjmYmIrIiKQhm9wDwhLvvBJ4I52eTcffdYbmzYPkXgC+6+9VAL/CJEtsjIrJkFASz2wM8HE4/TDDu8KKE4xS/B8iPY1zU9iIiy01BMLv17n4mnD4LrJ9jvYSZdZvZM2Z2V7hsHZB291w4fxLYPNcLmdl94XN0p1KpEpstIlK8ag2CBccsNrPHgQ2zPPSpwhl3dzOb69O129z9lJldCTwZDljfV0xD3X0vsBegq6uryj7FKyKrQUsLmEUwCNz9lrkeM7NzZrbR3c+Y2Ubg/BzPcSqsj5rZ3wE3An8NtJpZbXhW0AmcuoyfQURkWVRrV9SlXhraB9wbTt8L/HD6CmaWNLN4ON0OvAt42d0deAr40Hzbi4isJNXYzUSpQfB54FYzOwzcEs5jZl1m9tVwneuAbjP7GcGB//Pu/nL42O8Dv21mRwjuGXytxPaIiCypagyCBS8NzcfdLwLvnWV5N/DJcPpp4IY5tj8K3FRKG0REllM1BoG+WSwiUgQFgYhIxCkIREQirhq7olYQiIgUIZmE0VEYGqp0S8pHQSAiUoRq/HaxgkBEpAgKAhGRiFMQiIhEnIJARCTiFAQiIhGnIBARibi1a6uvK2oFgYhIEWpqgjBQEIiIRFi1dTOhIBARKZKCQEQk4hQEBcyszcweM7PDYZ2cZZ1/YWbPFZRsfgB7M/u6mR0reGx3Ke0REVkOCoKpHgCecPedwBPh/BTu/pS773b33cB7gCHgfxes8rv5x939uRLbIyKy5BQEU+0BHg6nHwbuWmD9DwF/6+5V1G+fiERNtXVFXWoQrHf3M+H0WWD9AuvfA3x72rLPmdnzZvbF/CD3IiIrWTIJIyOQyVS6JeWxYBCY2eNm9uIsZU/heu7uwJz5aGYbCcYuPlCw+EFgF/AOoI1gMPu5tr/PzLrNrDuVSi3UbBGRJVNt3y5ecPB6d79lrsfM7JyZbXT3M+GB/vw8T/WrwA/cfbTgufNnE8Nm9hfA78zTjr3AXoCurq4qOSETkdWoMAg2b65sW8qh1EtD+4B7w+l7gR/Os+5HmHZZKAwPzMwI7i+8WGJ7RESWXLWdEZQaBJ8HbjWzw8At4Txm1mVmX82vZGbbgS3A/5m2/bfM7AXgBaAd+KMS2yMisuSqLQgWvDQ0H3e/CLx3luXdwCcL5o8DM06g3P09pby+iEgl5IOgp6ey7SgXfbNYRKRI1XZGoCAQESnS2rVBrSAQEYmoWKy6uqJWEIiIXIZq6mZCQSAichkUBCIiEacgEBGJOAWBiEjEKQhERCKumrqiVhCIiFyGauqKWkEgInIZqunbxQoCEZHLoCAQEYk4BYGISMQpCEREIq6agqCk8QhERKIqHwT790MuB/E4JBJT63gc6utnn04kgnmz2Z9/fBzSabh4MSgXLgT1L/8yNDeX92cpKQjM7FeAPwCuA24KB6SZbb3bgT8BYsBX3T0/ktkO4BFgHfAs8K/dfaSUNomILIfWVti0Cb773aBcrkRisjQ0QG1tcJbR0xOEwXRdXfCWt1z+682m1DOCF4F/CXxlrhXMLAZ8GbgVOAkcNLN97v4y8AXgi+7+iJn9OfAJ4M9KbJOIyJKLxeDo0eBd+/BwULLZmdMjI5PLCqfz62SzwXcR8tMjI8HZxrp10N4+tV63DrZuLf/PUupQla8A2FznNoGbgCPufjRc9xFgj5m9ArwH+LVwvYcJzi4UBCKyKsTjsH59pVtRuuW4WbwZeLNg/mS4bB2QdvfctOWzMrP7zKzbzLpTqdSSNVZEJGoWPCMws8eBDbM89Cl3/2H5mzQ7d98L7AXo6uqqgt49RERWhgWDwN1vKfE1TgFbCuY7w2UXgVYzqw3PCvLLRURkGS3HpaGDwE4z22Fm9cA9wD53d+Ap4EPhevcCy3aGISIigZKCwMzuNrOTwD8FHjWzA+HyTWa2HyB8t38/cAB4Bfiuu78UPsXvA79tZkcI7hl8rZT2iIhI8cxXYWfaXV1d3t0961cWRERkDmb2rLt3TV+uLiZERCJOQSAiEnGr8tKQmaWAE5e5eTtwoYzNqQbaJzNpn8xO+2Wm1bRPtrl7x/SFqzIISmFm3bNdI4sy7ZOZtE9mp/0yUzXsE10aEhGJOAWBiEjERTEI9la6ASuQ9slM2iez036ZadXvk8jdIxARkamieEYgIiIFFAQiIhEXqSAws9vN7JCZHTGzByrdnkows4fM7LyZvViwrM3MHjOzw2GdrGQbl5uZbTGzp8zsZTN7ycx+M1we2f1iZgkz+4mZ/SzcJ/85XL7DzH4c/g99J+xIMlLMLGZm/2hmPwrnV/0+iUwQFAyZ+X7geuAjZnZ9ZVtVEV8Hbp+27AHgCXffCTwRzkdJDvgP7n498PPAb4R/G1HeL8PAe9z9bcBu4HYz+3kmh5e9GuglGF42an6ToAPNvFW/TyITBBQMmenuI8AjwJ4Kt2nZufvfAz3TFu8hGCqUsL5rOdtUae5+xt1/Gk4PEPyTbybC+8UDl8LZurA4wfCy3w+XR2qfAJhZJ3AH8NVw3qiCfRKlIJhryEyB9e5+Jpw+C1TBKKyXx8y2AzcCPybi+yW8BPIccB54DHidIoaXrVJfAn4PGA/nixpyd6WKUhDIIoQDBkXyM8Vm1gz8NfDv3b2/8LEo7hd3H3P33QSjB94E7KpsiyrLzH4JOO/uz1a6LeW24FCVVWSuITMFzpnZRnc/Y2YbCd4BRoqZ1RGEwLfc/W/CxZHfLwDunjazpwgGoIry8LLvAu40sw8ACaAF+BOqYJ9E6Yxg1iEzK9ymlWIfwVChEMEhQ8PrvF8DXnH3/1bwUGT3i5l1mFlrON0A3Epw7ySyw8u6+4Pu3unu2wmOH0+6+7+iCvZJpL5ZHCb5l4AY8JC7f66yLVp+ZvZt4BcJus49B3wG+J/Ad4GtBN17/6q7T7+hXLXM7J8B/xd4gclrv/+R4D5BJPeLmf0cwY3PGMEbxu+6+2fN7EqCD1q0Af8I/Lq7D1eupZVhZr8I/I67/1I17JNIBYGIiMwUpUtDIiIyCwWBiEjEKQhERCJOQSAiEnEKAhGRiFMQiIhEnIJARCTi/j8FRJLmdl7fVAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "\n",
    "plt.plot(range(len(twfiles)), shifts[0, :], \"b-\")\n",
    "plt.plot(range(len(twfiles)), shifts[1, :], \"g-\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "15bdab39",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0 0.10637643165419682 -0.11768223374643519\n",
      "1 0.10700657405766205 -0.13107961726782408\n",
      "2 0.10163923831678442 -0.12152158617752193\n",
      "3 0.09065352553736802 -0.12144904023219313\n",
      "4 -0.975754312382031 1.0456499145482434\n"
     ]
    }
   ],
   "source": [
    "# determine the aveage values for each tile\n",
    "#    set the range to be the number of tiles and the \n",
    "for k in range(n_tiles):\n",
    "    k1 = k * n_images_per_tile\n",
    "    k2 = k1 + n_images_per_tile\n",
    "    avex = np.average(sigma_clip(shifts[0, k1:k2]))\n",
    "    avey = np.average(sigma_clip(shifts[1, k1:k2]))\n",
    "    print(k, avex, avey)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "96016508",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python (jwstdev)",
   "language": "python",
   "name": "myenv"
  },
  "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.10.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}