{ "cells": [ { "cell_type": "markdown", "id": "0f9a755f", "metadata": {}, "source": [ "# ROEBA\n", "## Row-by-row, Odd/even by amplifier correction\n", "\n", "This is useful for JWST instruments where one want to improve upon the reference pixel step" ] }, { "cell_type": "code", "execution_count": 1, "id": "b2a6afbf", "metadata": {}, "outputs": [], "source": [ "from tshirt.tests import test_phot_algorithms\n", "import numpy as np\n", "from tshirt.pipeline.instrument_specific import rowamp_sub\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "id": "540616c5", "metadata": {}, "source": [ "### Test Simualated Images\n", "Simulate some \"row-like\" noise in an image that looks like 1/f noise.\n", "The simulated data has a noiseless star on top\n", "The following plots are\n", "* \"simdata\" - the simulated data\n", "* \"bkgmask\" - the mask of background pixels\n", "* \"gauss2d\" - the Gaussian star alone" ] }, { "cell_type": "code", "execution_count": 2, "id": "55e7b556", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/everettschlawin/es_programs/tshirt/tshirt/tests/test_phot_algorithms.py:129: UserWarning: Matplotlib is currently using module://matplotlib_inline.backend_inline, which is a non-GUI backend, so cannot show the figure.\n", " fig.show()\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATsAAAEICAYAAAAgMlPEAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAjAElEQVR4nO3deZQkZZnv8e8vsrIqu6uql2qapldBWQRUEFtwXFHUK+pcRq8yLiNuIzojLkfnjuicueocj8sdZXSuR7QZERkHkXFFRRER8aBgowwKNOC0TTe9VO9bVXd1VWXmc/+IqDa7qOrKqs58862K53NOnsqMjIw3Kp6MJ99434g3ZGY459xMl7R6BZxzLgRPds65XPBk55zLBU92zrlc8GTnnMsFT3bOuVxoSbKT9CFJ/9agZZ0oySS1NWJ5bnyS1kt64RjTz5e0qRXrVK/sO3Jyq9fDNZ6kj0j62kTztSTZmdnHzeyvQ5c7HXZK5/JO0ssk3SFpr6Stkv5NUvexLtcPY51zsZkLfAxYApwOLAX++VgX2vRkJ+kDkjZL6pP0sKQLaqudNYehb5a0UdIeSe+Q9HRJv8+y++drlleQ9GlJOyWtA142qrw3S3owK2+dpLdn0zuBHwFLJPVnjyWSzpV0Z1ZOr6TPS2pv9naZxp4uaU0Wp69IKo2eQdK7s3mWZa//Ptu2WyT9de0hpaRrJH1B0o+ymPxS0gmSPpuV8ZCkp9Ys+3JJf8ziu0bSK2reO1nS7ZL2Zd+Pb4z1D0h6dvZdO7/RG2c6k3SOpP/Ktu1/SvqGpI9Jmi/pB5J2ZDH5wUhss88d0bwxav8uSfqapF3ZPna3pEXZe2/K9tE+SY9Iej2AmV1nZj82s4Nmtge4CnhWzfJPyuLcJ+kW4Li6/kEza9oDOA3YCCzJXp8IPAH4CPC1mmkGfBEoAS8GDgHfBY4nzerbgedl878DeAhYDvQAt2Wfb8vef1lWhoDnAQeBc7L3zgc2jVrHpwHPANqydXkQeG8zt8t0fQDrgftrtv0vSX+BD29X4P8A9wALs9cvAbYCZwKzga9l8To5e/8aYGcWhxLwM+AR4BKgkC3/tpp1eDXpL34C/CVwAFicvfd14B+y90rAs2s+Z8DJ2fpsBM5t9faM6QG0AxuA9wBF4JXAULb9FwD/K4tfN/CfwHdHfS9eWPO6dv9+O/D97LOFLM5zgE5gP3BaNt9i4Mxx1u2zwPU1r+8ErgA6gOcCfSPlHfV/bPIGPJk0Ub0QKI6zMU7MvohLa97fBfxlzetvjSSgbGd4R817L6Ym2Y2xDt8F3pM9P7xTHmWd3wt8p9Vfvhgf2Ze6dtu/FPhjtl03Z1/AO4C5NfNcDXxi1HdidLK7qub9dwEP1rx+MrD3KOt0L3BR9vxaYBWwbIz5DPhgtkM/qdXbMrZHljQ2A6qZdgfwsTHmPRvYM+p7MV6yewvwK+Apo5bRCewlTaKzjrJeLwL2AKdmr1cAZaCzZp7r6kl2TT2MNbO1pMnjI8B2SddLWjLO7Ntqng+M8bore76E9Jd5xIbahUi6UNJdknZL2ku6Q45bzZV0alYt3yppP/Dxo83vHrPtR+I5D7iUNLHtq5lndLxqn4+oN/ZIukTSvdkh0V7gSfwpXn9PWqNfLekBSW8ZVc57gRvM7P5x/7v8WgJstix7ZDYCSJot6UuSNmT7yC+AeZIKdSz334GbgeuzZoz/K6loZgdIa+bvAHol/VDSE2s/KOkZpInsVWb2h5r13JN9fsQROWA8TW+zs/T4+9nA40h/XT91jIvsJT2MGrFi5ImkDtJa4KeBRWY2D7iJdAcgK3+0K0kPi08xsznAh2rmd481ettvyZ7vAV4OfEXSs2rm6QWW1byu/fykSHocafvNZcCCLL73k8XLzLaa2dvMbAnp4dMXdOTpJq8G/kLSe6a6DjNYL7BUUu13fyRW7ydtkjov20eem00fmfcA6WHqiBNGnpjZsJl91MzOAJ5J+h25JHvvZjN7Eekh7EOksU0XnLbT3gi8xcxuHbWe87M2+BErqENTk52k0yS9IEtCh0h/pavHuNgbgHdLWiZpPnB5zXvtpMfxO4CypAtJD3NHbAMWSJpbM62btO2gP/tl+ZtjXL+Z7p3Ztu8hbR873AlgZj8HXg98W9K52eQbgDdLOl3SbOAfj6HsTtIfrB2QdkaR1uzIXr+6puF8TzZv7fdtC3AB8B5JHucj3QlUgMsktUm6CBiJYTfpvrs3i/uHR332XuA1koqSVgKvGnlD0vMlPTmrBe4HhoGqpEWSLsqS1iDQTxYrSU8Cfgy8y8y+X1uQmW0AfgN8VFK7pGcDf17PP9jsml0H8EnSBuitpB0OHzzGZV5FWi3+HWlD+LdH3jCzPuDdpDvYHuB1pL8OI+8/RNqIvS47DFoC/F02X1+27DF78Nxh1wE/AdaRttd9rPZNM7uFtJ3m+5LOMbMfAf9K2pG0Frgrm3VwsgWb2RrgM6Q75jbS9rxf1szydODXkvpJ4/4eM1s3ahmPkia8yyUFP9czVmY2RNop8VbStrS/An5AGqfPArNI9+O7SBNRrX8k7RTcA3yU9Dsy4gTgm6SJ7kHgdtJD2wR4H+kP0G7SzsSRH6D3AwuBL+tPZ048ULPM1wHnZZ/7MGlb7YR05CG6c80l6XTSQ88OMyu3en3c+CT9GviimX2l1evSCH5SsWs6Sa+Q1JE1O3wK+L4nuvhIep7ScxzbJL0ReAqPrcVNW57sXAhvJz0F6Y+k7ULeXhan00ibh/aSHkq+ysx6W7pGDdS0ZCfpJUqvmFgr6fKJP+Gmg6nE1cxeYmZzzazHzF4xk3agmULSS0gTXB/puXVPMbMftni1GqopbXZZz8sfSE8I3ATcDbw2a2B205THdWbKS1ybNSzSucDakZ4wSdcDFwFjbrxCV6e1Lehp0qocqbxrN5X+A34e3dRMKq7H9RTsxOXFYCv3298P7jSzhcEKnDkmFVeAdnVYic7x3m6oQxxgyAaPeZ9tVrJbypFnym8i7So+TNKlpGfcU5wzn5Nf874mrcqR1n31iiDlzFCTiuuKpW2svnnK5xBPWmHx2rrOpHePMWFc4cjYlpjNebogyMr9+ohziqeuZQNemtkq0usY6Vix3A6sqAQpt+rjmTRVbVxXnlXy85pmkNrYzlHPtItts5LdZo68LGhZNm1sBYM5w01alTHKclM1ubi66SIXcW1WsrsbOEXSSaQb7TWkZz2PqbhPnHBTmCrXzn3eXHcMJhVXN23kIq5NSXZmVpZ0GellXQXgajN7YLz5yyXYe3KYU/7KtwcpZkaabFzd9JCXuDatzc7MbiIdcWRChSHo3hjm8LIwFKSYGWsycXXTRx7iGsUduUxQ6QhXlnMuf6JIdiRQ6QiUhfwCOedyKYpkV22Dg4vCHMZWo/iPnXOhxbHrGyTDgWp2fuaJc7kURbJLyjBre5gslPjAQs7lUhTJzhKozApTszNvs3Mul6JIdgjKs8KV5ZzLnziSHX546ZxrriiSnSpQ3B+mzU5hxhtwzkUmimRXGDLmbAxTtSsMeXesc3kURbJL+g4x6xcPhSmr/1CQcpxzcYki2Q0dP4tH3/LkMGVdfUuQcpxzcYki2akC7fvCleWcy58okl21CAcXB7pcLNwtEZxzEYki2bUdhOPvqQYpq/dgkGKcc5GJItlZAQbnhrm0wQpBinHORSaOZKdwo5H4eHbO5VMUya4yC3Y/JcxhbOXHQYpxzkUmimSnMnTsCnN8Kb8szblciiLZJRXo2B2uLOdc/sSR7IaMOY+GqXIlfrmYc7kURbIrHBym+97eYGU55/InimRX7mpn97OWhinrR2Fuxu2ci0sUya5ahAOLw5xn51dQOJdPcSS7NhhY6HcXc841TxS7fmEI5jwSriznXP5EkexUgXYfqdg510RRJLvCUJWuLYPBynLO5U8UyU4Vo23PQLCynHP5E0Wyq5QK7Dt9XpiyHvVhT5zLoziSXQfsPzHMqSeVjiDFOOciE0WyKwzCvD+G6TkohGkadM5FJopklwxXmd0bJgslw95B4VweHVOyk7Qe6AMqQNnMVkrqAb4BnAisBy42sz1HW85wZ0Lvn80+llWp2/B/hzlcnu4aFVsXlzzHtRE1u+eb2c6a15cDt5rZJyVdnr3+wNEWUOyrsvT2vgasysQ29nnNbhKOObYuSrmMazMOYy8Czs+efxX4ORNsuGpHQv+KMDW7qtfsjsWkY+umhVzE9ViTnQE/kWTAl8xsFbDIzEbGa9oKLBrrg5IuBS4FaO+cz1BXoBvueK6r15RiWxvXFUujaBJ2R2rIPlsiTOWkkY712/hsM9ss6XjgFkkP1b5pZpZt1MfINvIqgM4Fy60QalBNP6e4XlOKbW1cV55V8q0dn4bss3PUM+1ie0zJzsw2Z3+3S/oOcC6wTdJiM+uVtBjYPtFyCkNVujb55WIxaVRsXVzyHNcpJztJnUBiZn3Z8xcD/wTcCLwR+GT293sTLatcSthzWmmqqzIp5TV+HDuRRsbWxSPvcT2Wmt0i4DuSRpZznZn9WNLdwA2S3gpsAC6ecEkJlGcHuqGr57p6NC62Lia5juuUk52ZrQPOGmP6LuCCyS0MkmFvs4tFQ2PropH3uEbRXaYKdOz18eycc80TR7KzcPdzHbufyTk308WR7KpG20CYXlJVPds5l0dRJDtLRLkUpoPCkkAdIc65qESR7HzUE+dcs0WR7CwRw91hbujqNTvn8imKZFdtE4cWhBkuvdrmyc65PIoi2bUNVJh/375gZTnn8ieKZDfcXaD3ufPDlLXFb7jjXB5FkewKg8a8dcPBynLO5U8Uya7SIfY+PkwHRWW1t9k5l0dRJDsDLNDRpdfrnMunOJJdAYbmhCvLOZc/USS7Yr+x5I4wJxVv6fe6nXN5FEWyq8wSu8/oCFPWfd5m51weRZHsCoPGXO+Ndc41kY/b65zLhShqdpaI4a4wPQd+baxz+RRHsivAUFegIZ68N9a5XIoi2aVtdoFupehtds7lUhTJbqhLbHlOmFspDj3kh7HO5VEUya4wBF0bw9S4CkNBinHORSaKZAf4dVzOuaaKItmVu4wdzyyHKes2z6rO5VEUyY6KKOwPtCoVb7NzLo+iSHZtA7DgvjBlbRsIU45zLi5RJLukbMzaFeYwNin7YaxzeRRFsivPErvODDN4Z/k3fhjrXB5FkewwSMKMA+C9vs7lVBzJDqj6ZVzOuSaKItlZAuXZ4cpyzuVPFMlOFWjfH64s51z+RJHsEFTD9E+A9084l0vxJLtQa+LJzrlcmrAFS9LVkrZLur9mWo+kWyT9d/Z3fjZdkv5V0lpJv5d0Tj0rYVnNLsTDPNkdFiK2LjyP69jqaa6/BnjJqGmXA7ea2SnArdlrgAuBU7LHpcCVda1EBTr2hnkk3mZX6xqaHFvXEtfgcX2MCQ8ezewXkk4cNfki4Pzs+VeBnwMfyKZfa2YG3CVpnqTFZtZ7tDIqJWPfGWGuoKiU/ES7ESFi68LzuI5tqi1li2o2xlZgUfZ8KbCxZr5N2bTHbDhJl5L+klBqn8tpXw5z0eqenZ7sJnBMsa2N64qlcTQJO6DR+yyBzhVroGP+NpqZSZp0BjGzVcAqgO55y+zQ8WFGKraiN9rVayqxrY3ryrO8Gh2jRuyzc9Qz7WI71WS3baSqK2kxsD2bvhlYXjPfsmzaUVWKon9JmFpAxZPdRBoaWxeN3Md1qhnmRuCNwCezv9+rmX6ZpOuB84B99Rz7WxsMLAx0dzE/sppIQ2PropH7uE6460v6OmnD5nGSNgEfJt1gN0h6K7ABuDib/SbgpcBa4CDw5npWQmUoBWpLU5h+kGkhRGxdeB7XsdXTG/vacd66YIx5DXjnZFeiMGx09YY5J6QwPO2aGpomRGxdeB7XsUVxUDc0Dx69qBqmrN8FKcY5F5kokh0mNBhojCe/hMK5XIok2YGGAiUhP4p1LpeiSHZtB2Dhb8KUtf1AmHKcc3GJItmVO2HHykBl3R6mHOdcXKJIdm0HYcG9YQ5jtx8MUoxzLjJRJLtKO/SvCJPsKu1BinHORSaKZCeDwmC4spxz+RNHsquknRShynLxqliVU77zN6x75ZdavSpuhoki2VXbYLAnXFkuLmetfi2zvjn38OvTvnkvz7jrHYdfP/ld93HV8l+2YtXcDBLFrq9qwJpdmAs1XB2u2P14bnn9eSzbtpvy1gcPT68Cc7921+HXW+5YwYVzX8unv3c1Z7bPasGaupkgimTXdtBYeO8kG+1q+zMm0Q634aA32sXgmv3H85OnLcQGH2Si35/y+kcB+LszX8gn7v8ZZ3d0NH8F3YwTRbIrzxbbz6nzC6xRfyFNdnXmsPL9frlYqw1bha+fvhRscj9w1QMH+IfnvJIfrv5hk9bMNYUESlAhuyTUqljVwKpg4Sof9dxwJwzV97AkbXertEO5lP6ttoEV6lyGa7lzVr9hyp+1Awd5/gMXNXBtXFNJqL2dwpwuCgvmU1i0kGRBD4WuTtRWTBNhIFHU7ApD0P3oxI1pVkhHGi7PSmuD1gbJcHoY3DYAybBN2CZXGGrQSrspW3HZXspT/EWv7NlD8onHw9cavFKu8bJEl8ybiy3q4dDxnVRmFWg7UKFjax/J9l1U9/Vh5eEgNbwokt2ElN1btk2UZ8OhBeLQwipWqpAcKFDakVDabbSZSMqWnkvnTXPOtZQKBZLZs7ETFrDnyXPZfYYYnlehY0cHCx4oMvd+SIaGqfRXwJp/TlgUyW64G7Y+9+jZyRKDjirdPQd4zpL1XLxgNacU93Hv4PHcsONc7tmyjIG9JRhOUHX8qvHwrxq99m4yhhvwpZYZw1ahqEDDgrmpUYJmz2JgcRe7niSe9pyHeUHPQ3x369k8Uj2RWTu6aN9VQgcPYgHOkogi2SVDMHvj+F9cy9rbyrOMvsFuHigtZk3nMgDuH1jOw7uP51BvJx27C+mVGDb+lRKJH8a21HP/9zuZ03v3MS2jeOcazrryXaz52y80aK1cw0mQCNoKlGcnlOeXeVHPGl7V/Qg7h7t5eN4KyqUCHUkCCtN1EEWyUxk6dh+lZpd1TLQdFKq2sbm9h28lT+XOzsezoW8+OzfNY9b2Au37IBk6+mGs34Oite789Bd52e0vpbx5y5SXMfTMMz3Rxc4MqgZDwxT3l+nY1s71W57Oo8ct4Cdbnkhpe0KxfxArVwhSrSOWZGfQNtE9sgWFQUMVoUo7G/tOYEPpeJKDCbN2JXTsTTsqkgpHba/za2OdC8MqFWxggFJvPwvum8eWQ8u5bu4yOnaLBQ+VKW7bjx06lJ6GEkAUyS497JygzQ6l19AeNFQVxQMJ1baEZBgKh6Dt0J8S3VGX5cnOuTCsSvXAAMnWHcwtV+jc0kWlVKBwsEzbzj7YuQcbGIBqmAvW40h2gmphgvNtsrdVTWt4yXA2zbKL+y091MXSxDjRclzrPPWHG7n7qW1TOt2gsOh4Xv7/ftqEtXINZ4aVh6nu70eHBmnb2UFbUoBqBTs0iA0OYpVwI3NEkexUgY79dR63Z8mq9r45hw9N69h3fNST1vvwwnv5n13nU+3rm/yH53TxrvkbGr9SrjnMsOGh9Fy6gUMoUUuunoBIkl2lBHtOq/80grFuEFZvW1zFTz1puaIKfOq+W7j8f/wVlYfX1v255Owz+N4P/x3wU06mHTOwSqi+iDFFkexUTYdmD1WWa72ntJd4/Y238bmPX0zP/fux3z4w7rzlFzyNvuXtfOOf/pmiugKupZtJokh2VKFwKFCV1pNdNF7fvYvXf+JKXvfI81n9qz8D0hr6KR9fwx/+4YzD8/3thTfzvp51gCc6N3VRJLukbHRuD9OYlpS9OzY21510G5x02+HXL3ran7P29CtbuEZuJooi2RUGq3Su7w9WlovbLad/v9Wr4GagKJLdcHeB3ufOC1PWZm/cdi6Pokh2h8+VC1SWcy5/okh26X1jw2Qhv2+sc/kURbKjAJWuQG1pfhTrXC5FkeyKfXDCL8Jcx7VzCiftO+emvwmTnaSrgZcD283sSdm0jwBvA3Zks33IzG7K3vsg8FagArzbzG6eqIxyCfaeGmZMq/LtQYqJXoi4utbw2I6tnprdNcDngWtHTf8XM/t07QRJZwCvAc4ElgA/lXSq2dGHp23fX2XZrWFuHLul3mtwZ75raHJcXctcg8f2MSZMdmb2C0kn1rm8i4DrzWwQeETSWuBc4M6jfajantC/PMzNj6tr4rmhWiuFiKtrDY/t2I6lze4ySZcAvwHeb2Z7gKXAXTXzbMqmPYakS4FLAdrmzGf/SWGSUMXvrzyRhsV1xdIomoTdnzQstiVmN3lVG2+qGeZK4AnA2UAv8JnJLsDMVpnZSjNbWSx1UuwnyMMHAjiqhsZ14QLv+o5IY/dZpl+tYUo/vWa2beS5pKuAH2QvNwPLa2Zdlk07qsIhY97a4amsyqQFG3BgGmp0XF08PLZTTHaSFptZb/byFcD92fMbgeskXUHa2HkKsHqi5VWL4uDCMIc81aIPVTyeRsfVxcNjW9+pJ18HzgeOk7QJ+DBwvqSzSS++Wg+8HcDMHpB0A7AGKAPvrKdXR1UoDoQ5vvTD2FSIuLrW8NiOTRZ4aOSxzDphuT3+kvcFKWvdtVcwsHWjV+8CWHlWyVbfvHziGRuksHjtb81sZbACc2yOeuw8XRCkrF/brey33ce8z0bRXVY8UOWEu8IMVbzxgFftnMujKJKdFcTw3GKwspxz+RNHsktEuRTmPDtLPNk5l0dRJLvhbmPzi8K0HQ6vbn0bpXMuvCiSnYbFrM1hVkXDXrNzLo+iSHaFQZi7NkzHQWEwSDHOuchEkewQVNsC1bi8YudcLkWR7FSFtsEwbWl+UrFz+RRFsiscqjDn4X3BynLO5U8UyW5wfoF1r5oXpqwv+EgczuVRFMlOFSgd+9UgdZflnMufKJIdSu9DEaos51z+xJHs8I4D51xzRZHskiHo2hwm2yVDQYpxzkUmjmQ3bHRuCTNScTLsl4s5l0dRJLvybLH9nDBj2pfv90Y75/IoimTXNmAseCBMze7RAa/ZOZdHUSQ7S6DaEabGZX7bWOdyKY5kVxCD3WFO9vXBO53LpyiSHQZJJdDhpR/FOpdLUSS7wmCVOWsPBCvLOZc/USQ7K4ih+WF6Y/0w1rl8iiLZJcNVStvC3F0sGfaanXN5FEWyq5QK7D91Tpiy1vuoJ87lURTJLhmq0rU+TJtdMuQ1O+fyKIpkVykl7H1iV5iy1vqJds7lURTJThWjY3+YgeYU6hQX51xUokh25U6x7elh2tLKq7031rk8iiLZFfuMxXeUg5S1tc9rds7lURTJrtIh+laEWZVKoGtwnXNxiSLZFQaNORvC1OwKgW7Z6JyLSxTJrlISe08uhinrt16zcy6Pokh2qkCxP9BNsv3uYs7l0oTJTtJy4FpgEemYIavM7HOSeoBvACcC64GLzWyPJAGfA14KHATeZGb3HK2MSjv0Lw1T46q0BykmeiHi6sLzuI6vnppdGXi/md0jqRv4raRbgDcBt5rZJyVdDlwOfAC4EDgle5wHXJn9HZcVYGh+mCsbzK8WG9H0uLqW8LiOY8JkZ2a9QG/2vE/Sg8BS4CLg/Gy2rwI/J914FwHXmpkBd0maJ2lxtpyxJUalK9BlXIl3UECguLrgPK7jm1SbnaQTgacCvwYW1WyQraTVZkg37Maaj23Kpo278QoHRc89YapcOw56B8VozYqray2P65HqTnaSuoBvAe81s/3poX7KzEzSpKpMki4FLgVo75xP8UCgDgofB+AIzYzriqVR9H/lUqPjmi3zcGxLzG7UqgZT17dRUpF0w/2HmX07m7xtpLoraTGwPZu+GVhe8/Fl2bQjmNkqYBVAaely23N6oA6KnwUpZlpodlxXnlXyNoMWaEZc4cjYzlHPtIvthEOAZL01XwYeNLMrat66EXhj9vyNwPdqpl+i1DOAfRMd/6sKxX4FeXjNLhUiri48j+v46qnZPQt4A3CfpHuzaR8CPgncIOmtwAbg4uy9m0i7sdeSdmW/eaICkiHo3hAmCyVDQYqZDpoeV9cSHtdx1NMbewcw3jHmBWPMb8A7J7MSScVo7w+U7HyIJyBMXF14HtfxRdGCbAp3/pt5Z6xzuRRFsqsWRf/iMNmuWvRs51weRZHsMEiGw5XlnMufKJKdtcHA8WFqXBbFf+ycCy2OXd8C9pJ6zc65XIoi2SUVKO0Ok4USH+LJuVyKI9kNGd2bwjTaJUNetXMuj6JIdlYQg3PD9MZawXtjncujKJJdtQCD88LcvLrq49k5l0txJLt26FsRriznXP5EkewQdQxJ0MCynHO5E0WyK/YbJ9wZppt0a6Ab+zjn4hJFsrNElGcFOqk48aqdc3kURbIDQJ6EnHPNE0WyqxZgcG6YZOe9sc7lUxTJrjBozF0X5nqxwqC32TmXR1Eku0pJ7D05zDkhld/54bJzeRRFsgOo+pUNzrkmiiLZJWXo3Bbm1JOkHKQY51xkokh2ANU2r9k555onimSXDFXp2jgQrCznXP5EkeyqxYSBE0rBynLO5U8cya4NBnoCjXoSxX/snAstil2/MGR0bwrTc1DwwTudy6Uokl15tthxdjFMWf/lHSHO5VEUyc4SGO4MU+Myb7JzLpeiSHbteyuc9L2+IGVt2+t33HEuj6JIdsPdBTY/f06Ysh71kQCcy6Mokl1ShtlbA91K0a+gcC6Xokh2qkLboTAn+8rPKXYul6JIdpZAeVaYngPvoHAun6JIdhgkw4HOf/PT7JzLpSiSXVI2OgL1kiZlz3bO5VEUya5aFAPHhVmVatFPKnYujybMMJKWA9cCi0gPAleZ2eckfQR4G7Ajm/VDZnZT9pkPAm8FKsC7zezmo5VRKcGe06f8P0xK5dYw5cQuRFxdeB7X8dVTnSoD7zezeyR1A7+VdEv23r+Y2adrZ5Z0BvAa4ExgCfBTSaea2bjHqW0HYNHdYbpJtx4IUsx00PS4upbwuI5jwmRnZr1Ab/a8T9KDwNKjfOQi4HozGwQekbQWOBe4c7wPJENVOjcdnNSKT5WPZ5cKEVcXnsd1fJNqKJN0IvBU4NfAs4DLJF0C/Ib012QP6Ya9q+ZjmxhjY0u6FLgUoNg9n23ndU9l/SetvN7PPRmtWXFdsTSKJuHcamRcs+Udjm2J2c1b8Sape8+X1AV8C3ivme0HrgSeAJxN+kvymckUbGarzGylma0sFjsp7aoGecivoDhCM+O6cIFfmtcqjY4rjNpn6Wjk6gZR10+vpCLphvsPM/s2gJltq3n/KuAH2cvNwPKajy/Lph2lAKi0B+ol9c7Yw5oeV9cSHtex1dMbK+DLwINmdkXN9MVZ+wDAK4D7s+c3AtdJuoK0wfMUYPXRyqi2wWBPmCzkIxWnQsTVhedxHV89u/6zgDcA90m6N5v2IeC1ks4m7d5eD7wdwMwekHQDsIa0Z+idE/bsCMphbkHhNbs/aX5cXSt4XMchs9ZfUSBpB3AA2NngRR83xjIfZ2YLG1yOG0PguILHNhhJfcDDDV5sU+MaRbIDkPQbM1sZ+zLd5HhcZ6bpGFc/D8M5lwue7JxzuRBTsls1TZbpJsfjOjNNu7hG02bnnHPNFFPNzjnnmsaTnXMuF1qe7CS9RNLDktZKunyKy1gu6TZJayQ9IOk92fSPSNos6d7s8dLGrr0bj8d1ZmpEXLPlBI9tS9vsJBWAPwAvIh1t4W7gtWa2ZpLLWQwsrh3DC/gL4GKgf/QYXq65PK4zU6Pimi0reGxbXbM7F1hrZuvMbAi4nnR8rUkxs14zuyd73gdMNIaXay6P68zUkLhCa2Lb6mS3FNhY83rcsbTqNWoML0jH8Pq9pKslzT+WZbu6eVxnpobHFcLFttXJrqGaMYaXaz2P68wVMratTnYNG0trvDG8zKxiZlXgKtJquGs+j+vM1NCx70LHttXJ7m7gFEknSWonvfHHjZNdyNHG8KqZrXYML9dcHteZqSFxhdbEtqVDWZpZWdJlwM1AAbjazB6YwqImNYaXay6P68zUwLhCC2Lrl4s553Kh1YexzjkXhCc751wueLJzzuWCJzvnXC54snPO5YInO+dcLniyc87lwv8HwbTgZPKdoi8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "rab = test_phot_algorithms.rowAmpBacksub()\n", "simDict = rab.sim_data()\n", "rab.show_images()" ] }, { "cell_type": "markdown", "id": "9ebd8555", "metadata": {}, "source": [ "## Now test the simulation with ROEBA applied\n", "The plots are now\n", "* \"simdata\" - the simulated data\n", "* \"bkgmask\" - the mask of background pixels\n", "* \"gauss2d\" - the Gaussian star alone\n", "* \"model\" - the ROEBA model from background pixels\n", "* \"resid\" - the residual of the ROEBA-applied image and the original Gaussian 2D" ] }, { "cell_type": "code", "execution_count": 3, "id": "7cca739a", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWMAAAEICAYAAACK8ZV4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAA1/0lEQVR4nO29eZxkdXnv/37Oqequ3qanZ5gZZgUURMEFENEoEtSYBIzhapS43LgHzRWXG5NIzC+RJCYhv0s03uuVBKLBDZHEDbcoGpefEURFFBhAh2Fg9q1n6bWWc57fH+f0UNPT3dXVdepb36p63q9Xv7rq1KnzPJ/zPeep73m+m6gqhmEYRmsJWu2AYRiGYcHYMAzDCywYG4ZheIAFY8MwDA+wYGwYhuEBFowNwzA8oCXBWETeIyL/ktGxThURFZFcFserw+42Efm1ObZfLCI7XPpSL+n5Or3VfrhARK4WkU+22o+5EJEbReR9i9x3zuvNJ0TkuSLyYKv9cI2IvFpEvrHA598RkTfVOk5LgrGq/q2q1nQua9ohUBoLIyIvEpHvi8hhEdkjIv8iIkOt9qsbmf2jrqr/n6qe2UqfWoGqfkpVf73R41iawmg3hoH3AeuAJwHrgf/VUo+MjsD10/Vsmh6MReTdIrJTRMZE5EEReUH1o2NVmuH1IrJdRA6JyFtE5Bki8vO0BvShquOFInKtiBwQka3Ai2bZe72I3J/a2yoib063DwBfA9aJyHj6t05ELhCR21M7u0XkQyLSs0h5zxCRzanP/yoihTn0vz3dZ0P6/k9SO7tE5E3VtYv0sfXDIvK11L//EpGTReQfUxsPiMi5Vce+SkQeSrVuFpGXVH12uoh8V0SOpOfqM/OUz4Xpeb94IaEicp6I/DS19W8i8hkReZ+IjIjIl0Vkf+rjl2e0pt877vF6VtkXROSTInIwPf8/EpE16WevS8tvTEQeFpFXA6jqTar6H6o6qaqHgBuA51Qd/7RU95iI3AactJCuxZBq+OP0epwQkY+IyJq0nMZE5JsiMpLu+9sicl+q5zsi8qSq45wrInel3/kMUJhl57dE5O70uz8Qkac26vsi9T0p9fVw6vtvp9uPe7xOy+T76evvpZt/ll6rvyuznjzrOW/p/q8RkUfS6+HPZ187TdK+TZIY9XNgIr0ffpCei59V3xfzXZPV5yV9/8L0Xj0iSeySRTmjqk37A84EtgPr0venAo8HrgY+WbVNgX8iuTh/HZgGvgCsJqn57AN+Nd3/LcADwEZgBfDt9Pu59PMXpTYE+FVgEjgv/exiYMcsH58OPAvIpb7cD7xzEdq2AfdW+fFfJDW2YzaAvwDuAlal738T2AOcDfQDn0x9Pz39/EbgQOpTAfhP4GHgNUCYHv/bVT68nKSGGAC/C0wAa9PPPg38WfpZAbiw6nsKnJ76sx24oIbWHuAR4B1AHngpUEr9WQn8TqpnCPg34AuzztOvVb2vLvs3A19KvxumupcBA8BR4Mx0v7XA2fP49o/AzVXvbwfeD/QCFwFjM/YauI63AXcAa3jserwLOLeqnN4LPCEtgxem5+lPgC3p+Zs5h/8z/exlQBl4X2rj3PS4z0zPxWtTu71znccM79F86uN7Uh+fn56zM4HvAG+q2vd1wPdnX0dV7y+m6v5a7HlL9z0LGAcuTP24Nj0/mWueo2zvJrmP1wMHgUtJ7psXpu9XLXRNVp8Xkh//sbR882l5V6rP47y+NFno6WkB/BqQn+eGPDUt1PVVnx8Efrfq/WdJA2RagG+p+uzXqQrGc/jwBeAdc10s8+z/TuDziyzEaj8uBR5KbewkCQjfB4ar9vko8Hezzs/sYHxD1edvA+6vev8U4PACPt0NXJa+/jhwPbBhjv0U+FOS4PDkRWi9KNUkVdu+TxpIZu17DnBo1nmaLxi/AfgB8NRZxxgADpME+b4F/HohcAh4Qvp+U3rhD1TtcxPZBONXz7oer5tVTl8A/hy4pWp7kJ63i9NzuGvWOfwBjwXj64C/nmX3QR6rhBx3HrP6A55LUkEIqrZ9Oi2n79B4MK553tLXfwF8uuqzfpIffBfB+A3p63cDn5j1+ddJfhjnvSY5Phi/Brij6jMBdrCIYNzUNIWqbiEJblcD+0TkZhFZN8/ue6teT83xfjB9vY6kNjfDI9UHEZFLROQOERkVkcMkQXLeR1UReUL6aL1HRI4Cf7vQ/rOY7ceMtuXAFSSB90jVPrN9r349w2LPw8xj3cxj7WHgyVW+/wnJhXBn+uj5hll23kkSOO6dV93xfu/U9Oqq9l1E+kXkn9PHy6PA94DlIhIu4rifILnYb5YkbfP/ikheVSdIavpvAXaLyFdE5InVXxSRZ5EE2pep6i+q/DyUfn+G466PBlhMuayrtqeqMcl5Ws/c57Dat1OAd82UZVqeG3nsmmoW64Dtqa/Vfq3P6PhLuq9VdZKkUuaCGbunAC+fVQYXkjxt1rwmU2brUOa+z0+g6TljTXJ8F5IIVeDvGzzkbpKLdIZNMy9EpJfk1/daYI2qLge+ymM5m7mmqLuOJO1xhqouI3lcW1yO50Q/dqWvDwG/BfyriDynap/dwIaq99XfrwsROYUkX3olsDLVei+p76q6R1V/X1XXkaQDPizHd2d7OfDfROQdizC3G1gvItXnZcb3d5E80j4zPX8XzbiY/p8gqeXMcPLMC1Utq+pfqupZwLNJztlr0s++rqovJHkcfCDVOqP9XOBWkhrNt2b5OSJJ+8AMm3DHLpLrHID0fG0kqR3PdQ6rfdsO/I2qLq/661fVTzvweaOIVMeCTanP85ZdEzju3hCRPpIUmAtm4sJ2kppxdRkMqOo1sPA1WcVx8anqGqhJU4OxiJwpIs9Pg+Q0yS9hXONrtbgFeLuIbEiT/1dVfdZDkivcD1RE5BKSNMYMe4GVIjJctW2IJBc0nv7S/UEdvrw19WMFSX72WCOZqn4HeDXwORG5oMr316cNJv0kj7VLZYDkItoPScMlSc2Y9P3L5bGGtEPpvtXnfhfwAuAdIlJL8+1ABFwpIjkRuQyY0TREUq6H0/Pw3lnfvRt4hYjkReR8klzajI/PE5GnpLXooyQ5wjht5LksDapFklxinH7nycB/AG9T1S9VG1LVR4AfA38pIj0iciHw4hrasuQW4EWSNFLnSX6oiiTpiNtJUihvT8/FS3nsHEJyY79FRJ4pCQOSdONrdre9H5K0q/xJ6tfFJOfsZpKye2n69HM68MZZ390LPC4jP/4deLGIPFuSBvSrWXylKCs+mfrwG5J0FChI0ii5YaFrchZfAc4WkZdK0jvj7SzyR6zZNeNe4BqSRqk9JA1yf9rgMW8gebT9GUljwOdmPlDVMRLxt5AEoFeR1KBmPn+AJB+2NX0MWQf8UbrfWHrsOXsdzMNNwDeArST54uM68KvqbSR50S+JyHmq+jXgf5M0Om4hadyApHDrQlU3A/9AcpPvJckn/1fVLs8Afigi4yTn4B2qunXWMR4lCchXyQKd0lW1RNJo90aSvNl/B76c+v2PQB9JGd9BEiir+XOSBtVDwF+SnLMZTia5CY+SNJx+lyR1EQB/SPKDMUrSEDvzg/EukgaVj8hjvWLuqzrmq0gawUZJfhg+Pp+urFHVB0nOzf8hOR8vBl6sqqWqc/i61Lff5fhr98fA7wMfIjlXW9J9m+1zKfXzktTnDwOvSe+VD5DkbfcCHwM+NevrVwMfS++lyxv04z6SHPLNJLXLcZL2prrvjQZ82A5cRvJ0vJ+kpvzHJNfjQtdk9TEOkDx1XkOSZjmD4+/LeZHjU1iGSyTp9nQvSYt5pdX+1IOI/BD4J1X911b7YnQeIjJI8sN/hqo+3GJ3nGCDPhwjIi8Rkd40xfL3wJfaIRCLyK9K0uc5JyKvBZ7KibVgw1gyIvLiNCUyQNLucw9Jb4euwIKxe95M8vj1EEketp4cdSs5kyQ1dJgkVfAyVd3dUo+MTuMykjTALpLH+1doFz26Ny0Yi8hvSjLibouIXFX7G+3HUjSq6m+q6rCqrlDVl/ge0GY0knSV+4CqDqrqU1X1K632LUvsem09qvqmtAfDsKq+IM3B143vOuelSR2pQ5Ka3+NIejj8DDirGbZa9WcaO+evG3R2g8Z219msiTEuALZo2novIjeTPIJsnmvncHBAcytXNMmV46kcHCUan8iiy0xdGgFOWhHqqRvzGZiuzbbtZQ6MRo3qrFtjj/RqgYH5Ps6UaSYoadF5WQ6O5HXl+hOmIWkaj943fkBVVzV4GCvLOQgHBjS3wlHsGR0lmpg/9jQrGK/n+FEnO0i6Gx1DRK4gGaVGftkIp7/iD5vkyvFs/dj7szpUTY1wvM5N63Pc+fUlj/Ooiwt+Y1GDfmpRt8YC/TxTXpCF7Zr88LjxHg1R1/WaWzVM/1+/MyvbtXnpX2YxitDKMuW4shweYdMf/M+sbC/I9v/7gQU/b9mUcap6PcncCfRu2qgTmyInduPFzseWEdU6z39aoSMbI6o1LpMVHa9x45OH9Y1n3u7M9h87s9R9Zdm7aaNWBtzI1BoTBDQrGO/k+CGAG9JtcxMqLCs3yZU5bGVDfRrbk27QCHXqHAqmeP7A/U13KmOsLD2nWcH4R8AZInIayYl4BcnIqDnJHxFO/qqbKuuBI5mNsKxLY5vSDRqhTp1bD67hlR9z82ibkEkKz8rSc5oSjFW1IiJXkgxbDoGPajLccU4qBTh8upsuz5XvZnOcejW2I92gEerXmRsos+JZe5z5lwVWlv7TtJyxqn6VZMa0moQlGNruJm8TlrI7Vj0a25Vu0Aj16SxP5dl77+ome5Q9VpZ+09I1n2ZQgajXnS3DaIhAifobnXzQMI7Hi2BMAFGvoyhpA8CNRglA+t30/jG6By+CcZyDyTVu0hSxF4qNdqavt8Q5j3vUmb2sliox/MaP0KQQlB3VjDuy56ThkkgDjpT6Wu2G0WF4EYyDCvTtcxMlA+8nqzR8p1TMsfXhNa12w8gCxV0FrYYdL4KxBhD1uakZq+WMjQaRUMkPZtgtx2gZEkN+3E3skRptvl4EYwQqrp76rDeF0SBhGLNieKL2jhnxkDNL3YeGUB7q7OHQdWPpA6NdiGNhouh4khOj4/EiGEsE+aNufp3EeiQZDRJXAsb3u5la0ugevAjGYUlZtt1N1TgsWXcKozGCKWHZZjfzUhvdgxfBOBibpu97D7ixNT7txI7RufSMFln/qV86s3evM0tGK/EiGJdW9/HoG57ixtZHb3Nix+hcKqeG7P/giDuDL3JnymgdXgRjiaDniDtbhtEIlUrIgX3LWu2G0WF4EYzjPEyudTQc2lJ9RoP09ZZ42uMzWdZqUbgbeG20Ei+CcW4SVt/lZhas3ZNOzBgdTPFIL1u+9vhWu2F0GF4EYw2hOOxmaFytjteGUYs4B1OrbQpNI1v8CMbibjY1m8/YaJgAtMe6SBrZ4kUwjvpg9KluahrRfzgxY3Qwq4eOcuXF7nrluFwd2mgdXgRjqUDvQTf5A7Fh10aDHCkX+Nresx1a/IZDW0ar8CIYBxH0jrqzZRiNUCzleWjHqla7YWSAxJCbtFnbjhGUlGWPuqmyBjYc2miQYEoYuKfQajeMDFCByNGcT7Xaq7wIxuFkmaG7dzuzZRiN0HM0YsM3Dzuzt9mZpe5EQz8qaF4E48pgD6PPWe/G1tds6kOjMYorQh66fNidwZ+6M9Vt2OTys4jzMLHWTT9jG4FnNExvTHjGeKu9MDLAJpefRZyDqVW2OrTRHgz1TnPRKVuc2fuFM0tGK/EiNIUlWPawO1uG0QhHJ/r4xl1uZhk0ugcvgrFE0GMrfRhtglSEnv02rr4TkAjyY45yxjVijxfBOCzFDO4qOrNlGI0QFmHYXZbCaCIaQHnQUc64RrOYF8FYIiV3aMqZLcNohKCsDOyxLpIdgeBuxfh26GccFUKOPGm5G1uP2uOl0RjlQWHPrzjsImnzqXQFfgTjXjh6qpuubVGvEzNGByN9ET3nHGq1G0aH4UUwDouw/CE3LWuhm9S00cHoZEj5xw7XwDOahkSQm7AGvGME5Zj+3W6iZFC2BjyjMcIiDG+166gT0ACivg5owBORbcAYEAEVVT1fRFYAnwFOBbYBl6vqgs905YGA3b/S34gri6b8y/rTIVnp9Jlu0AjZ6IxGIo681OEIvE/Vt3s3lGVmGqV2kMwMBw14z1PVA1XvrwK+parXiMhV6ft3L3SA/FjM+u+OZeBKbbaPLblG07DONqAbNEKDOnMHQ1beNNBsHxulG8qyozQ2I01xGXBx+vpjwHeocULi3oDxTW5qxvESasbzULfONqQbNEKdOqMCjD6x7XrldENZtrXGRoOxAt8QEQX+WVWvB9ao6sx8mHuANXN9UUSuAK4A6BkYoTToaEHSpZnJROem9V6k6OcjE40F3PyoNsCSdFZrDEdGKK7wOmfcDWWZicZwxJ+G2Eajw4WqulNEVgO3icgD1R+qqqYn6wTSk3c9wMDKjRq6mvR9aWYy0Xn+0wo+jzjJROMyWeGzRliizmqNhfUbNZz2emXbbijLTDT2btrojcaGgrGq7kz/7xORzwMXAHtFZK2q7haRtcC+WscJSzGDO/wdDp2VTp/pBo2Qjc5wGkYe8OYePoFuKMvMNCpLraDVTw07Sw7GIjIABKo6lr7+deCvgFuB1wLXpP+/WOtYlULAoTPdLGNT2VxfniJLnb7SDRohO52VQdj3LIfB+BOL37UbyrJTNTZSM14DfF5EZo5zk6r+h4j8CLhFRN4IPAJcXvNIAVT6HT321Z8zzk6nv3SDRshKZxgjw97OxdoNZZmdxk6Ym0JVtwJPm2P7QeAF9R0smXzFCXWayVSnp3SDRshQZyzE0342xnZDWWapUWIIi7bs0jEkgt7DNp+x0R5IRejZ58WtYzSIUnvV5ixtLYQXV5QoBI6C5Nztq4ZRBzEEjmpTRpMJIO5xFBTaYj7jWMlNuem3KbFFY6MxJIbcZKu9MDoNL4KxBkKl4KamoYHVaIzGSCaXabUXRqfhRTC2WduMdiIswjKbtc3IGC+CsQZCeSjvzJZhNILmYHqFq6m+jG7Bi2Ac54TplW4mXolzFoyNxojyMLHe2h46gST/b13bjpGbihi554gzW4bRCD1jyvrvVpzZ2+rMUvehAlGvo8nl22FB0vJQyO6L3MyeVN7VdlMfGp5RWgmPvtLhj/pX3ZnqNiSG3ISblFNb1IzDorJ8q5ulz8OiPV4ajRFMCEN3uZlLxWguGkJ5mZvGWK1RD/QiGEe9wuHHuWnAi+60nLHRGPGAMnbedKvdMDoML4KxUvtXI0tbhtEYSpCzrm1GtvgRjEMoLXNnyzAaoScfccrqUWf2HnZmyWglXgTj/Liy7vtuBn3sGre6sdEY0WieIzevb7UbRofhRTCO+oTRs3rd2LrHcsZGY0SDyuiFDuczvsGdKaN1eBGMw6IybL0pjDYhmBSGfuqm8mB0D14EY8NoJwSbitXIHi+CsQZCedBNy5rNTWE0ShxCcXmrvTA6DT+CcQilQUdTaFpvCqNBNAfFk2xYvZEtXgTjJGfspjeF5YyNRslNwkl3uZu17RFnloxW4kUwLg0Ku57rZnhp6QFLUxiNEawo0/+q3e4M/qs7U0br8CIYhyUY3O6mxhp6u8K60S6UpvLsuOfkVrthdBheBGPAxikb7YVdr0bGeBGMK4PK/me7mR+28m27i4zG2LB8lGt++yZn9i5/lzNTRgvxIhgTCeFRR65EljM2GmM86uUHY6c7tHiXQ1tGq/AiGOemYOU9bmztnXJjx+hcDo8N8MX/fKZDi7c4tGW0Ci+CcVBR+g66SVMEFUtTGI0RlGBgpz1hGdniRTCu9AkHz3YzuXzlx3YTGQ0yFBFd5GbNRgA+4M6U0Tq8CMYoBG7mCbJWcKNhVIVS0Y9bx+gcvLmiYhumbLQLouTyNhzayBYvgrEGUOl3Z8swGiEfxqxdftSZvQedWTJaiRfBWCLocXRti1VojAYplUMe3bei1W4YHYYXwRiB2E37XTIZrWE0QBAovQVXjRxGt+BPMHbliQVjo0FElN68BWMjW2pmUEXkoyKyT0Turdq2QkRuE5Ffpv9H0u0iIv9bRLaIyM9F5LzFOKFpzdjFn84TjF3obDXdoBGarzMMlGWForO/Vmj0gW7QWM1i6qM3Ah8CPl617SrgW6p6jYhclb5/N3AJcEb690zguvT/ggQR9B6uy+8lE8yfM76RJuv0gBvpfI3QZJ2lcsije1qeM76Rzi/LG2myRokgN+HmcblWe1XNYKyq3xORU2dtvgy4OH39MeA7JCfkMuDjqqrAHSKyXETWquqCk79GBeXIWW5G4EWFuTsau9DZarpBIzRf57r+I7z3GV/K3O/5eP0c27qhLF1o1ACiPjeDD2r15FpqpnZNlcg9wJr09Xpge9V+O9JtJ5wQEbkCuAKg0DPMmR9xM2nEoQN1nfhMdW5a70eKfhbZliWO+ijWT0M6qzX29C3nhvf8TnO9PY47FrtjN5RlphrDkRF33V1rVMAbjg6qqiL1r5WrqtcD1wMMLd+g06vdrPSh+aU9kmSh8/ynzVMt94QsNC6TFV5rhKXprNbYt2ajHt3k9yilbijLLDT2btrojcalBuO9M48AIrIW2Jdu3wlsrNpvQ7ptQaK8ML7OTa0xqi8YZ6rTU7pBI2SoMy4oR5/gJq1WJ91Qlh2rcakR8FbgtcA16f8vVm2/UkRuJkmeH1lMXkpzMLXK0erQ9SnOVKendINGyFBnrjdizamjzfT1OB5d/K7dUJYdq7FmaBKRT5MkzE8SkR3Ae0lOxC0i8kaSxWsvT3f/KnApsAWYZO62hxNtVKBQXy53ycg8FRoXOltNN2iE5uusFEP2Ptra3hTdUJbdoLGaxfSmeOU8H71gjn0VeGu9ToRlZXC3m3HKYXne3hRN19lqukEjNF9nUBSGftnaxthuKMtu0FiNF837peXw6GWxG1s/c2LG6GA2rjrAtX9wgzN7l1zrzJTRQrwIxqggRUet0/MNwTOMRVLUHFtLqx1avN+hLaNVeBKMQUqOgqQ3HVmMdqWsIbvLy1vthtFheBGMcxOw6sdubO2bcGPH6FxGjwzx6a9e5NDiF2vvYrQ9XgTjygDsP9+Rre+6sWN0LiuGx3j1pe4upL+6ypkpo4V4EYxzk7Dybjdpin2TTswYHczBo0N84psua8a3OrRltAovgnHUA+Ob3ATjqMeJGaOD6ekrs+Epe5zZe9iZJaOVeBGMRSGce9rWptgyjEaoxAEHxgda7YbRYfgRjKOkEc+VrXYi0pgzPv8HbH3pP7faFSMljgImjriZ2MrwBCWZdU1BIjlWq9Nw1ucN4EUwjnNQdDS61NnyTg3wtDtfSd+/Dx97f+a/382z7njLsfdPeds93LDxv1rhmgH05Cucsv6gM3t1zE1hZIzEQCzEPQqiBCVJtomkAViJ8woqhEUhzumSg7IXoUlihzVjNwP9lsT7Rx/Hba9+Jhv2jlLZ81hH/xgY/uRjc9ru+v4mLhl+Jdd+8aOc3dPXAk+7m3IUsmt0WavdMDJAIsiPzRM9VQjKyeRilQhEBYmSGKJp5VjLQlAWFMhNCRKlgXseWwvhRTDOTSqr7q4zaVx9/urIAz8y6WfS+Majq/nG01ehxfup9XtR2ZbUlf7o7F/j7+79T87p7W2+g8YxgqMBy26znHEnoAFUBuaOCXEhRnsjwr4IdvcSTAlBRdBQiXsgnBYkVoIpobgqQjeWiCZzSDEgKJ4Y4Ju10kemVPqFfectMqDIrP+QBONFxtjKvf4Nhy5rxKeftB60vh+keGKCP3vuS/nKnV9pkmdNRgQkQMI08aYxGitoDOrnjyZAPBwxecmYO4P/4s7UkmnTspQYwsk5AmcecuMBmgtA80l6UyDuSVMVkSIKcQhxXsmNBcjhAqGAVIS4V5HyibYWwotgDCw6z6JB+hcm/yXm2KODzymIhTjvzt9j3RLnH9CJSZ5332V8++w2G6UlgvT0EPQVkEIBcjm0XIapaeKpabRS9vYmFmAJC0x0Lm1clhpCZXCWbwqaU+gTtDemsDtHZSBGVCivrBBMhsS9MZWyQAjBZACi9EwGTJ1SIhjPoYEiFTkurmmN6Xe8CMZhCYYerR1JNUxW6qj0JbVpzUFQTtIcuSkIylozIIeljJzOkE1XHqayxIs1OnSI4O8eB5/M2KlmMnPzLh9G16xgevUAUV9IbiKid88Ywb6DxEfGvL2JdSokvme49o7dQJuX5ZykATQsCloJKA8mgVUFwrEQDSCYTnIOUkpzyIFQWqaER3JIJMnCx3U+hHsRjGsiScI8zgmVfpheKUyvitFCRDARUtgfUBhVcioEleTxwSYE8hcJQ4L+fvTklRx6yjCjZwnl5RG9+3tZeV+e4XshKJWJxiNQP/si2uR/Ce1elkm32lmFqSBxgIZKUAnQQNEA4p4T874aJgtWhMWky5to0qMinArR4PiA3BYNeOUh2HPRwtFTA4XemKEVEzx33TYuX3knZ+SPcHdxNbfsv4C7dm1g6nABygESz3+nlH+QtfeNUc7gAhVVyhqRF78XyTyGBEh/H1NrBzn4ZOHpz32Q5694gC/sOYeH41Pp2z9Iz8ECMjmJeph6WjNymHe+3F1a6K1/4cxU/bR5WWoAUd/xsScaqZDbn6dwQJheqVRWVijsyLPsIeXo6ULUq5z5f3fz4JVr6TkU0L9XGTsFihtL9OzO03NYmDo5JhqKCMfC42wthBfBOChB//b5A4kKIFDpU8aKQ9xXWMvmgQ0A3Du1kQdHVzO9e4De0TAZyafzj7QLPEtTXPTHb2XZ7h81dIz87Zt52nVvY/P/+HBGXjUREQgEciGV/oDKSIUXrtjMy4Ye5kB5iAeXb6JSCOkNAhBXa6jXx/7iINdvudChRU9nt+qAskRODJK5A3kQKA1r0hA3HRD1KaNPSYLKGZ8YY8ub1nHardNsvayX4oqkZwXFgKhXKQ0nteTwSC7JPVfZWggvgrFUoHd0gZpxesJyk4LEOXb2rOCzwbncPvA4Hhkb4cCO5fTtC+k5AkFp4TTFfGvgtYrbr/0nXvTdS6ns3LXkY5SefXZ7BGJI8oaxQqlM/miF3r093LzrGTx60kq+seuJFPYF5MeLaCXCy6oUEJVCDu2ynHEnlOW8KEQFRXOKlJPa8Awv/tT3uPa7l/CMf/wJW7/zbCCpXYeTAXEhGQCSDAypz6QfwVghN1VrJwiLmuRloh62j53MI4XVBJMBfQcDeg8nDXlBxIL5YmsEbz0aRejUFIXd46y8Zzm7pjdy0/AGekeFlQ9UyO89ik5PJ12jfCQSckfaJCXUZNq+LOcgzicdAcJpIe5TtFeRohCUhMryiA/+7PlIKeDf7j8P7VHCsSQIJ4E4yR9HfVozLTEbL4JxklaokTMmGf2Sm1QkFvITAXEuIChDOA256ccC8YLHap9ronPRmHhiimDPfoYrEQO7BokKIeFkhdyBMThwCJ2agti/Bp9j2HWU0AllOQdxPulPHEwLxOnIu/SJO46TKBtHkjb2QVASKCUDQjTteyxRfVVjP4KxQBzWcDz9WOKkhhyUqZq4I/mvQfp/oecDD1vBz/3Kdn50bm5JXX/CNav5rf/zzSZ41URU0UqZ+Og4Ml0kd6CXXBBCHKHTRbRYRCOPb14BzVs0Btq/LOcinQwozmsy7FkUzSUDPILpAJ3uTcLIgeR/VHisS62mA0Mkmpm7YvFmvQjGEkHv0UXmlFJx1V2LjqUeFnF/+Dhr23tX3c1vD15MPLaEUV3LBnnbyCPZO9VsVNFyKel/OjWNBNIWI7YAggoU9nvaINUK2rgs50TSPsZBEi80AESRCKLBmNzRkMryCuHRkKg/JpwI0dzMyLzk+1KGuLc+7V4E46gAh85cfA5urj6ei80FR551bQPIS8jf33MbV/3Gfyd6cMuivxeccxZf/MongDbOX6qCRm3VviODFXqe7W7WtrahDctyPqJeRfNK/khAnH/sqVsiQSqQO5Q7bjrNmQa7oAzl5TEyx9wUtfAiGEucLL3kypaPPLWnwKtv/TYf/NvLWXHvUfQn9827b+X5T2dsYw+f+av/RV4GHXppAERRwJEj/a12w2gyUhbKy2KknATWoCzkjkrSc6IoRD1KbiJAIW28U8rDmqx0v4R0qBfBmBjCaUePM54GY4BXDx3k1X93Ha96+Hnc+YNfAZIa/xl/u5lf/NlZx/b7H5d8nT9csRWwQNwSIiEez7faCyMDkjaoBSJn1WdK8lQeTiXBNpxOUxmk/YwVgsr8x2qLiYKCijKwz00yN6j4n8O66bRvw2nfPvb+hU9/MVuedF0LPTKqCUrC4MNe3DpGg8wE2EWRdhggqHovVU1VNY5TK/J4cUWFxZiBbePObLUbtz3pS612wagiN6WsvLdce0fDfyTphubK1kJ4EYzLQyG7L1ruxtbONm7sMrygtELZ/iqHQznbdLrqdkBiyI+56RnTFmmKY32FHdkyjIZQIS5b17ZOQEMoL3PztNwW8xlHPTC+yU2UjHqcmDE6mOG+KV705Hud2bPWgu7Ai2BMmHSmdmXLMBqhLyjx1MHtrXbD6DC8CMb5MTj5e27GKR9wuHSZ0ZnsPbScD3zmvzm0+D2HtoxWUTMYi8hHgd8C9qnqk9NtVwO/D+xPd3uPqn41/exPgTcCEfB2Vf16LRuVAhx+gpscXGWeqWFd6Gw13aARmq8zP1jm5OfsbJL3J/KLObZ1Q1l2g8ZqFlMzvhH4EPDxWds/oKrXVm8QkbOAVwBnA+uAb4rIE1QXXs6i52jMhm9NLNrpRtg1/xwYN9JknR5wI52vEZqsU/fnKf7L2mw9rp8b6fyyvJHO13iMmsFYVb8nIqcu8niXATerahF4WES2ABcAty/0pbgnYHxj3yJNNEa8ee4auAudraYbNELzdVb6YPSs1vam6Iay7AaN1TSSM75SRF4D/Bh4l6oeAtYDd1TtsyPddgIicgVwBUBu2QhHT3NzcUe9dX8lM52b1nuRop+LzDQW8HrOhiXrrNaYX7WMvnNGHbi7JLqhLDPTGI6MNNnVxbPUCHgd8HjgHGA38A/1HkBVr1fV81X1/HxhgPw4Tv7qnCgoU52rVnrZlSPbsqT+XztHNKSzWqMMDHHkcL+zP1caZ+v0tCwz1RgODmTs3tJZUlVNVffOvBaRG4Avp293Ahurdt2QbluQcFpZvsXN8NJ6JiTKWqePdINGyFZnOCEM31nI3MdG6Yay7GSNSwrGIrJWVXenb18CzPSAvxW4SUTeT5JEPwO4s9bx4rwwucrNI3ycX3wXuqx1+kg3aIRsdUa9MH6Kf0M5u6EsO1njYrq2fRq4GDhJRHYA7wUuFpFzSAYXbwPeDKCq94nILcBmoAK8dTGtmRJDfsrNoI/50hQudLaabtAIzdcpEeSPtnb9rm4oy27QWM1ielO8co7NH1lg/78B/qYeJ6IeOHKqm3zqfMOhXehsNd2gEZqvU/oj5NwjS3EtM7qhLLtBYzVeNO/nJ2JOvsPNUh/bJ9pvCk3DL4JDIUOfG2q1G0aH4UUw1lAoD7tZOUFrrUJtGDWI8zC12mZtM7LFj2AcCJWCm4tbAwvGRmNoAGVvu+Aa9SAR5MbdxIRa0wR7EYzLQ8rOF7ppnS7f6V8ruNFerFs5ytWv/ZQze698nzNTXYcGEPW7iQlao77pRTCWstC3040rMyu9GsZS2Vcc4kMPP9+hxZ84tNVlSO0gmaWthfAiGIdFGN7ipmEtLDoxY3QwlYk8e+88udVuGB2GF8EYgTjnqMZqFWOjQVQgthVjjIzxIhhLDLmim7xNnXNTGMYJSAy5SftVN7LFi2AcTkcse9BNJ/pwuq0G5Rgekp9QTv6hm7lUYO7J5Y3Ow4tgXBwJ2fqy5W5sfdjLmdOMNiJcW2LkPdvcGfyaO1NG6/AiGEsEhVE/+voZRi2my3nu3bmu1W4YHYYXwRhJ1sFzZcswGiEMYoaHplrthtFh+BGMsYY1o32IVShVLN1lZIsXwTgoweBON9E4KDkxY3QwOh0y9cDyVrthdBh+BOOyMrDLTet0ULbh0EZjhEUYti4ORsZ4EYwr/cK+89yst1W515LGRmNEwxETl465MzjvDL5GJ+FFMM5NKSvvc1MzfnTKasZGYwRjIT3fXdZqN4wOw4tgrAHEvW5qrM4mBTE6Fg0g6mu1F0YmxBAUHT0t12gW8yMYh0JxyE3rtE0ubzRKnIPpk+wJy8gWL4IxCkHk6OK2e8hoEFGQSqu9MLJAFAJHI9ulRuzxIhiHxZhlWyac2TKMRshNwuqfuLuOHnJmqfuwyeVnoaFQGnHTm8LSFEajxDmYtDXwOgObXP54gnJMYa+b1aGDstWMjcYIS8qyRyxPYWSLF8E4KoQcfYKbrkLRNhvGajRGeUDYe4Gb1cwB+Io7U0br8CIYB6WYwW1ucsZByWrGRmPkJuGkn7ub/s8G+3UHXgTjqBBw+ImDbmxtsVyf0RjlQdjzLIfX0WfdmTJahxfBWCKl96ibmoa46kJndCxBBQoHrCHYyBYvgnFlQNj7DDe53MqddhMZjdG3fJqzL3vAmb3Nf+/MlNFCvAjG+TFl7ffdtE7vGbOasdEY06MFHvz0E1vthtFheBGMo15hbJMbVyJHc2AYnUtUUA4/2bq2GdniRTAOi+76bYZFqxkbjRFOCcvvcXfrPOrMktFKvAjGUUE4fLqbfpvRT6xmbDRGNKAceUax1W4YHYYXwVgiyI+7qbHa6tBGw1QEGXU46MPoCmoGYxHZCHwcWEMy59n1qvpBEVkBfAY4FdgGXK6qh0REgA8ClwKTwOtU9a6FbEQ9ML7eTY016jlxmwuNPtANOl1ozPeVOfmsfc0TMYttc2yzsuwMjdUspmZcAd6lqneJyBDwExG5DXgd8C1VvUZErgKuAt4NXAKckf49E7gu/T8vGkJpxM3IOJ27B13TNXpCN+hsusaBXIlnrdrWPAWzuGPuzVaWnaHxGDWDsaruBnanr8dE5H5gPXAZcHG628eA75CckMuAj6uqAneIyHIRWZseZ24CJRp0NEw5ODEd4kSjB3SDThca+8MS5w0+0jwRi8DKsjM0VlNXzlhETgXOBX4IrKkSuYfkUQKSk7W96ms70m3znpBwUlhxl5tBH/snF06HNEujb3SDzmZp3HV0Of/Pt38nc3/n54cLfmpluXSNEkFuwk2KtFZ71aKDsYgMkoySf6eqHk3SMwmqqiK15rE/4XhXAFcA9AyMkJ9w1IC3QAU8a43pMY/p3LTei/bSppZlgf4sXV0yzdSYWz5C/pAfs/9ZWTamMRwZIepro8nlRSRPcjI+paqfSzfvnXkEEJG1wEyLxk5gY9XXN6TbjndM9XrgeoDC+o166EmOGvD+c+7tzdAIx+s8/2mFlndybnZZLpMVHa9x5Imr9akX/rJp/s9mvpU+rCwb19i7aaO2zeTyaQvlR4D7VfX9VR/dCrwWuCb9/8Wq7VeKyM0kyfMjtXI2EkN+3NGjwhw1YxcafaAbdLrQWIxCth5ambnv9WBl2Rkaq1lMzfg5wO8B94jI3em295CciFtE5I3AI8Dl6WdfJelasoWke8nraxkISjD0iJsGvKA05+ama/SEbtDZdI3xZI6Jn6/I2O26sbLsDI3HWExviu8zfwX7BXPsr8Bb63EiiJSecUfBeI4pNF1o9IFu0OlCo0TQe6i1IzmtLDPSqLhbMb4dVodWmbf/b1NsGUZDSLIoqWFkiReXVJwXxte6icZx3qKx0RhxrzJxxtz5LqO9kBjyY25a8BbqyQWeBGMUgrI7W4bRECpI0Y+ubUZjaADlIUejf7Po2tZsNAdTq93UWNULxUY7k+utsPq0g87stXasX+cTVPx4WvYjNOm8vRyaYsswGiGOhbGp3la7YWSAqLuZHGsNTfEiGAcRFEbdRMnAptA0GiQuhUztdLOaudFcNEhWbnFlayH8CMYlZWiHm6RxULKqsdEY4TQs3+xq2JalKboFL4KxhkJx2E2DiIZ+5IeM9iV22MZhNBnrZ3w8cQjF5W5qGrE1ghsNonkorrF8V6cgkR8/rH4E4x4Y2+TOlmE0Ql9fkac+eZsze5amaB6iEDrqPNAWDXgI4MnMSYZRi0CUQs5Vx3ijmWgAlX5rwDtGflw5+XY3j317HC18anQuU4cKbP7sE1vthtFheBGMNRAqfY4GfQRWNTYaQwOoDLTaCyMLJIL8WJut9NF0xIKk0T44m5DcaD6ezB7mRTCOQygOuzkh1pvCaBTNK9NrKq12w8gADaG8rKUr0x/Di2AcFpXhrW6aNMOi5YyNxggnhZU/dfer/qgzS0Yr8SIYRwXh8Olu+pxFP/PjkcRoX6J+5eDTrZ+xkS1eBGOA2EbGGe2CAHk3j7ZG9+BFMA4qMLDXTU0jsFSf0SBBUeh/yEYPGdniRTAGiHNWMzbaAw0gshk0OwNR4l5HTzk1huB5EYyDUszg9ilntgyjEXJTsGKzXUcdQSzkxh01xsYLVzi9CMZxPmDq5IIzW4bRCFEvHHmcXUedgCiIo9RlW8xNEedgaoWjWdu8UGy0M3EPTG2w3hSdgAYQ9dncFMcIS8rQDjc/T6FNLm80SDgFw5tt9JCRLV4E40q/sP+cvBtbP7WGQqNBhivIb7hbkJQPuTNltA4vgrEGUB7w41HBMGrRE0ZsGj7szN7PnFkyWokXwbjncMRpXxxzYmvvYcv1GY1R2dvD6AdOabUbRofhRTAuD4XsfN4yN7YetVyf0RjRypiDvzfhzuDn3ZkyWocXwTioQP8eN2kKG4FnNEo8HVLa4qbyYHQPXgRjiSE37aYTvVhffaNBJIb8uDUEG9niRTDWACp9blrWrAHPaBQNobTMukga2eJFMEYhKDu6uO0eMhpEYgiLrfbCyAStvRxSlrYWwotgHFSUXke9HIKKRWOjMYIK9O21NEUnIAph0dEaeO0QjOO8MHWSG1fivN1ERmNEeZjYYD/qnYAGUOn3Y4xDzQgoIhuBjwNrSCra16vqB0XkauD3gf3pru9R1a+m3/lT4I1ABLxdVb++kI2oAIeeVMuTbIi+deI2Fxp9oBt0utA4MDTNMy58oEkKTmTrHNusLDtDYzWLqY5WgHep6l0iMgT8RERuSz/7gKpeW72ziJwFvAI4G1gHfFNEnqCq8+YhchOw5kduujnsmbt7aNM1ekI36Gy6xqnDBX7+JUe1h/mxsuwMjceoGYxVdTewO309JiL3A+sX+MplwM2qWgQeFpEtwAXA7fN9ISjFDOyYrMvxpTLXfMYuNPpAN+h0oTGchpFftPb+trKck7bTWE1diVoRORU4F/gh8BzgShF5DfBjkl+wQyQn646qr+1gjhMoIlcAVwDkh0bY+8yhpfhfN5VtCydustSYHu+Yzk3rvUjRA80rywL9zXW8DpqlMVy5nJ2XlZvrfDWfXfhjK8sGynJkpLmO18Gio4OIDJJcFu9U1aMich3w1yS5nL8G/gF4w2KPp6rXA9cDDKzcqIWDjgZ9LDACL2uNcLzO859W8KLVp5lluUxWdLzGwoaNmt/hx7pLVpaNaezdtNELjbDIYCwieZKT8SlV/RyAqu6t+vwG4Mvp253Axqqvb0i3LWAAoh5HvRzmMdN0jZ7QDTqbrVEFIg9+V60sO0PjDIvpTSHAR4D7VfX9VdvXpjkdgJcA96avbwVuEpH3kyTRzwDuXMhGnIPiCjfBeK6VPlxo9IFu0OlCo/TEBOvdtHHM64OVZUdorGYxNePnAL8H3CMid6fb3gO8UkTOIXlU2Aa8GUBV7xORW4DNJK2hb63ZmilQcbME3nw14+Zr9INu0Nl0jWEQMzw43RTn68DKsjM0HkNUvXjc2g9MAAcyPvRJcxzzFFVdlbGdRSEiY8CDGR92Lo3QIp1N0ggelaXj6xU6qyx90+hN7PEiGAOIyI9V9Xzfj9kIptG/4y4VK0t/jtkovui0OcwMwzA8wIKxYRiGB/gUjK9vk2M2gmn077hLxcrSn2M2ihc6vckZG4ZhdDM+1YwNwzC6FgvGhmEYHtDyYCwivykiD4rIFhG5aonH2Cgi3xaRzSJyn4i8I91+tYjsFJG7079Ls/W+Lh87XqdpXPQxvNaY+tLxOr3TqKot+wNC4CHgcUAP8DPgrCUcZy1wXvp6CPgFcBZwNfBHrdTYLTpNY2do7BadPmpsdc34AmCLqm5V1RJwM8mcpHWhqrtV9a709RhQa95T13SDTtO4SDzXCN2h0zuNrQ7G64HtVe/nnRd4scjx855CMu/pz0XkoyLSqslLu0GnaVwCHmqE7tDpncZWB+NMkVnzngLXAY8HziFZMeAfWudddnSDTtPYGRqhO3RmobHVwTiz+UdlnnlPVTVS1Ri4geTRpBV0g07TWAcea4Tu0OmdxlYH4x8BZ4jIaSLSQ7KY4K31HkRk/nlPq3arnvfUNd2g0zQuEs81Qnfo9E9jK1s00xbIS0laIB8C/myJx7iQZG7TnwN3p3+XAp8A7km33wqsNZ2m0TSaTh812nBowzAMD2h1msIwDMPAgrFhGIYXWDA2DMPwAAvGhmEYHmDB2DAMwwMsGBuGYXiABWPDMAwP+P8Biv16Cj9oftgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "outimg, modelimg = rowamp_sub.do_backsub(simDict['simdata'],\n", " amplifiers=1,\n", " backgMask=simDict['bkgmask'])\n", "simDict['model'] = modelimg\n", "simDict['outimg'] = outimg\n", "simDict['resid'] = simDict['gauss2d'] - simDict['outimg']\n", "\n", "rab.show_images(simDict)\n" ] }, { "cell_type": "markdown", "id": "7edd7c9c", "metadata": {}, "source": [ "The residuals don't look perfect. But further examination shows they are within machine errors" ] }, { "cell_type": "code", "execution_count": 4, "id": "b4e2641f", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "4.440892098500626e-16" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.max(np.abs((simDict['resid'])))" ] }, { "cell_type": "markdown", "id": "3a67df7f", "metadata": {}, "source": [ "## How to use it with the JWST pipeline?\n", "\n", "You can run ROEBA after the superbias step like so.\n", "This assumes you have already run the dq init and saturation flagging steps and called the result saturation" ] }, { "cell_type": "markdown", "id": "757ad9ec", "metadata": {}, "source": [ "\n", "``` python\n", " from jwst.superbias import SuperBiasStep\n", "\n", " superbias_step = SuperBiasStep()\n", " # superbias_step.output_dir = output_dir\n", " # superbias_step.save_results = True\n", " \n", " # Call using the the output from the previously-run saturation step\n", " superbias = superbias_step.run(saturation)\n", "\n", " mod_refpix = deepcopy(superbias)\n", " \n", " ngroups = superbias.meta.exposure.ngroups\n", " nints = superbias.data.shape[0] ## could be split into ints per segment\n", " \n", " for oneInt in tqdm.tqdm(np.arange(nints)):\n", " for oneGroup in np.arange(ngroups):\n", " \n", " rowSub, modelImg = rowamp_sub.do_backsub(superbias.data[oneInt,oneGroup,:,:],\n", " backgMask=simDict['bkgmask'],amplifiers=1)\n", " mod_refpix.data[oneInt,oneGroup,:,:] = rowSub\n", "```" ] }, { "cell_type": "code", "execution_count": null, "id": "34350150", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.8.11" } }, "nbformat": 4, "nbformat_minor": 5 }