updating image stamps nb

Author: christinawilliams

DP1/100_How_to_Use_RSP_Tools/103_Image_access_and_display/103_5_Image_stamps.ipynb

@@ -718,6 +718,17 @@
     "    return fig, axes[:n_subplots]"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "96462574-05d0-4827-9c6e-95494307f316",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import sys\n",
+    "sys.exit()"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -828,138 +839,6 @@
     "make_gif(tempdir)\n",
     "display(dimg(data=open(tempdir+'/animation.gif', 'rb').read()))"
    ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "8364f14c-5056-413f-b796-ad0dad6e2432",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "import sys\n",
-    "sys.exit()"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "0f78775c-d83d-400d-8d48-1bd9527be028",
-   "metadata": {
-    "scrolled": true
-   },
-   "outputs": [],
-   "source": [
-    "primary_header = scihdul[0].header\n",
-    "wcs_sci = WCS(sci_header)\n",
-    "cd_sci = wcs_sci.wcs.cd\n",
-    "theta_rad = np.arctan2(cd_sci[0,1], cd_sci[0,0])\n",
-    "theta_deg = np.degrees(theta_rad)\n",
-    "\n",
-    "print(\"method 1 for rot angle:\", theta_deg) #cd_sci, theta_rad, \n",
-    "\n",
-    "# Get the pixel scale matrix\n",
-    "cd_sci = wcs_sci.pixel_scale_matrix\n",
-    "# Rotation angle (radians → degrees)\n",
-    "theta_sci = np.degrees(np.arctan2(cd_sci[0,1], cd_sci[0,0]))\n",
-    "print(f\"Rotation angle: {theta_sci:.3f} deg\")\n",
-    "\n",
-    "print('rotation angle in header: ', primary_header['ROTPA'])\n",
-    "#print(sci_header['ROTPA'])\n",
-    "theta_sci = 180 - primary_header['ROTPA']\n",
-    "\n",
-    "\n",
-    "fig = plt.figure(figsize=(6, 6))\n",
-    "ax = fig.add_subplot(111, projection=wcs_sci)\n",
-    "ax.imshow(sci_hdu.data, cmap='gray', norm='asinh', origin='lower')\n",
-    "\n",
-    "wcs_ref = WCS(ref_header)\n",
-    "\n",
-    "# Get the pixel scale matrix\n",
-    "cd_ref = wcs_ref.pixel_scale_matrix\n",
-    "# Rotation angle (radians → degrees)\n",
-    "theta_ref = np.degrees(np.arctan2(cd_ref[0,1], cd_ref[0,0]))\n",
-    "print(f\"Rotation angle: {theta_ref:.3f} deg\")\n",
-    "\n",
-    "# instead of doing the above, grab the header keyword:\n",
-    "# ROTPA which is 180-theta_ref\n",
-    "\n",
-    "fig = plt.figure(figsize=(6, 6))\n",
-    "ax = fig.add_subplot(111, projection=wcs_ref)\n",
-    "ax.imshow(ref_hdu.data, cmap='gray', norm='asinh', origin='lower')\n"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "7313a3a3-6401-4f43-ac8d-9e1f6970da53",
-   "metadata": {},
-   "source": [
-    "### 5.4. Warp the template\n",
-    "\n",
-    "The `template_coadd` needs to be warped to match the orientation of the science and difference cutouts. Perform the image warping of the `template_coadd` cutout using the `Warper` from `lsst.afw.math`. The warped reference image is defined as `warped_ref`."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "7663f1d6-9d07-43eb-83c4-ab34244fdf99",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "warper_config = WarperConfig()\n",
-    "warper = Warper.fromConfig(warper_config)\n",
-    "\n",
-    "sci_wcs = sci.getWcs()\n",
-    "sci_bbox = sci.getBBox()\n",
-    "\n",
-    "warped_ref = warper.warpExposure(sci_wcs, ref, destBBox=sci_bbox)"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "db5865c9-7cd4-4bb2-9bf7-5fb550d4df41",
-   "metadata": {},
-   "source": [
-    "### 5.5. Display the cutout triplet\n",
-    "\n",
-    "Display the science, reference, and difference cutouts for the SNIa candidate associated with the DiaObject 611255759837069401."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "cecaaa36-ec2d-462c-b2ed-cc2f47100b36",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "fig, ax = plt.subplots(1, 3, figsize=(8, 3))\n",
-    "plt.sca(ax[0])\n",
-    "display1 = afwDisplay.Display(frame=fig)\n",
-    "display1.scale('linear', 'zscale')\n",
-    "display1.mtv(sci.image)\n",
-    "plt.title('science')\n",
-    "plt.sca(ax[1])\n",
-    "display2 = afwDisplay.Display(frame=fig)\n",
-    "display2.scale('linear', 'zscale')\n",
-    "display2.mtv(warped_ref.image)\n",
-    "plt.title('template')\n",
-    "plt.sca(ax[2])\n",
-    "display3 = afwDisplay.Display(frame=fig)\n",
-    "display3.scale('linear', 'zscale')\n",
-    "display3.mtv(diff.image)\n",
-    "ax[0].set_axis_off()\n",
-    "ax[1].set_axis_off()\n",
-    "ax[2].set_axis_off()\n",
-    "plt.title('difference')\n",
-    "plt.tight_layout()\n",
-    "fig.suptitle('DiaObject 611255759837069401, MJD='+str(np.round(scitab['t_max'][0], 2)))\n",
-    "plt.show()"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "6f6f513b-cb08-47ba-b8aa-ea18d18fa320",
-   "metadata": {},
-   "source": []
   }
  ],
  "metadata": {