updating image stamps nb

Author: christinawilliams

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

@@ -31,14 +31,14 @@
    "id": "ae760c59-bf94-45b3-9f8e-b4551b22e193",
    "metadata": {},
    "source": [
-    "**Learning objective:** An overview of candidate transient (Ia Supernovae) identification in DP1.\n",
+    "**Learning objective:** Generate many image cutouts using the Rubin image cutout service.\n",
     "\n",
     "**LSST data products:** `DiaObject`, `Visit`, `visit_image`, `template_coadd`, `difference_image`\n",
     "\n",
     "**Packages:** `matplotlib`, `numpy`, `lsst.rsp`,`lsst.afw`, `lsst.geom`, `lsst.resources`, `lsst.utils.plotting`, `pyvo.dal.adhoc`, `reproject`\n",
     "\n",
     "**Credit:**\n",
-    "Originally developed by the Rubin Community Science team with feedback from Eric Bellm. This notebook also utilizes a transient detection from LSSTComCam first identified by Dan Taranu.\n",
+    "Originally developed by the Rubin Community Science team with feedback from Gregory Dubois-Feldman. This notebook also utilizes a transient detection from LSSTComCam first identified by Dan Taranu.\n",
     "Please consider acknowledging them if this notebook is used for the preparation of journal articles, software releases, or other notebooks.\n",
     "\n",
     "**Get Support:**\n",
@@ -95,25 +95,20 @@
    "source": [
     "import matplotlib.pyplot as plt\n",
     "import math\n",
-    "\n",
     "import numpy as np\n",
-    "from astropy import units as u\n",
     "\n",
     "from lsst.rsp import get_tap_service\n",
     "from lsst.rsp.service import get_siav2_service\n",
     "from lsst.rsp.utils import get_pyvo_auth\n",
     "\n",
-    "from lsst.utils.plotting import (get_multiband_plot_colors,\n",
-    "                                 get_multiband_plot_symbols)\n",
     "\n",
     "import lsst.afw.display as afwDisplay\n",
     "from lsst.afw.image import ImageF\n",
     "from lsst.afw.fits import MemFileManager\n",
-    "import lsst.geom as geom\n",
-    "import lsst.resources\n",
     "\n",
     "from pyvo.dal.adhoc import DatalinkResults, SodaQuery\n",
     "\n",
+    "from astropy import units as u\n",
     "from astropy.wcs import WCS\n",
     "from astropy.io import fits\n",
     "import io\n",
@@ -173,26 +168,6 @@
     "assert sia_service is not None"
    ]
   },
-  {
-   "cell_type": "markdown",
-   "id": "9b02d9d2-23ce-48c4-8265-460604817467",
-   "metadata": {},
-   "source": [
-    "Define filter names, plot markers, linestyles, and colors for plotting"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "c5f87c79-8798-4306-bb82-495edb1463b2",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "filter_names = ['u', 'g', 'r', 'i', 'z', 'y']\n",
-    "filter_colors = get_multiband_plot_colors()\n",
-    "filter_symbols = get_multiband_plot_symbols()"
-   ]
-  },
   {
    "cell_type": "markdown",
    "id": "cc51bd61-5bea-47d8-b67b-30ff46f6e180",
@@ -334,16 +309,12 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "def get_cutout(dl_result, spherePoint, session, fov):\n",
+    "def get_cutout(dl_result, ra, dec, session, fov):\n",
     "    sq = SodaQuery.from_resource(dl_result,\n",
     "                                 dl_result.get_adhocservice_by_id(\"cutout-sync\"),\n",
     "                                 session=session)\n",
-    "    sphereRadius = fov * u.deg\n",
-    "    sq.circle = (spherePoint.getRa().asDegrees() * u.deg,\n",
-    "                 spherePoint.getDec().asDegrees() * u.deg,\n",
-    "                 sphereRadius)\n",
-    "    # can remove spherepoint stuff using:\n",
-    "    #sq.circle = (cutout_ra, cutout_dec, Radius)\n",
+    "\n",
+    "    sq.circle = (ra * u.deg, dec * u.deg, fov * u.deg)\n",
     "\n",
     "    cutout_bytes = sq.execute_stream().read()\n",
     "    sq.raise_if_error()\n",
@@ -353,6 +324,15 @@
     "    return ImageF(mem), hdul"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "3d998864-d33a-449d-a28a-42e94607f901",
+   "metadata": {},
+   "source": [
+    "#### 1.2.4. subplots\n",
+    "Define a function to organize the display with matplotlib of many cutouts."
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -403,10 +383,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "#SNcandID = 611255759837069401\n",
     "ra = 53.124767650110215\n",
-    "dec = -27.739814591611168\n",
-    "spherePoint = lsst.geom.SpherePoint(ra*geom.degrees, dec*geom.degrees)"
+    "dec = -27.739814591611168"
    ]
   },
   {
@@ -484,8 +462,7 @@
     "tx = np.where((lvl2_table['dataproduct_subtype'] == 'lsst.visit_image')\n",
     "              & (lvl2_table['lsst_band'] == 'r'))[0]\n",
     "scitab = lvl2_table[tx]\n",
-    "scitab.sort('t_max')\n",
-    "print(scitab['t_max'])"
+    "scitab.sort('t_max')"
    ]
   },
   {
@@ -579,9 +556,9 @@
    "outputs": [],
    "source": [
     "fov = 0.002\n",
-    "sci, scihdul = get_cutout(dl_result_sci, spherePoint, get_pyvo_auth(), fov)\n",
-    "ref, refhdul = get_cutout(dl_result_ref, spherePoint, get_pyvo_auth(), fov)\n",
-    "diff, diffhdul = get_cutout(dl_result_diff, spherePoint, get_pyvo_auth(), fov)"
+    "sci, scihdul = get_cutout(dl_result_sci, ra, dec, get_pyvo_auth(), fov)\n",
+    "ref, refhdul = get_cutout(dl_result_ref, ra, dec, get_pyvo_auth(), fov)\n",
+    "diff, diffhdul = get_cutout(dl_result_diff, ra, dec, get_pyvo_auth(), fov)"
    ]
   },
   {
@@ -752,7 +729,7 @@
     "    dl_result_sci = DatalinkResults.from_result_url(scitab['access_url'][j],\n",
     "                                                    session=get_pyvo_auth())\n",
     "\n",
-    "    sci, scihdul = get_cutout(dl_result_sci, spherePoint, get_pyvo_auth(), fov)\n",
+    "    sci, scihdul = get_cutout(dl_result_sci, ra, dec, get_pyvo_auth(), fov)\n",
     "    sci_hdu = scihdul[1]\n",
     "    sci_header = scihdul[1].header\n",
     "    scidata = scihdul[1].data\n",
@@ -804,7 +781,7 @@
     "    dl_result_sci = DatalinkResults.from_result_url(scitab['access_url'][j],\n",
     "                                                    session=get_pyvo_auth())\n",
     "\n",
-    "    sci, scihdul = get_cutout(dl_result_sci, spherePoint, get_pyvo_auth(), fov)\n",
+    "    sci, scihdul = get_cutout(dl_result_sci, ra, dec, get_pyvo_auth(), fov)\n",
     "    sci_hdu = scihdul[1]\n",
     "    sci_header = scihdul[1].header\n",
     "    scidata = scihdul[1].data\n",