WIP

Author: Leanne Guy

08a_Interactive_Image_Visualization.ipynb

@@ -10,7 +10,7 @@
     "Last verified to run on <b>2021-09-17</b> with LSST Science Pipelines release <b>w_2021_33</b> <br>\n",
     "Contact authors: Leanne Guy <br>\n",
     "Target audience: All DP0 delegates. <br>\n",
-    "Container Size: medium or large <br>\n",
+    "Minimum Container Size: medium <br>\n",
     "Questions welcome at <a href=\"https://community.lsst.org/c/support/dp0\">community.lsst.org/c/support/dp0</a> <br>\n",
     "Find DP0 documentation and resources at <a href=\"https://dp0-1.lsst.io\">dp0-1.lsst.io</a> <br>"
    ]
@@ -110,10 +110,9 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# Prevent some helpful but ancillary warning messages from printing\n",
-    "# during some LSST DM Release calls\n",
-    "warnings.simplefilter(\"ignore\", category=FutureWarning)\n",
-    "warnings.simplefilter(\"ignore\", category=UserWarning)"
+    "# Ignore warnings\n",
+    "import warnings\n",
+    "warnings.simplefilter(\"ignore\", category='UnitsWarning')"
    ]
   },
   {
@@ -166,7 +165,17 @@
     "calexpId = {'visit': 192350, 'detector': 175, 'band': 'i'}\n",
     "calexp = butler.get('calexp', **calexpId)\n",
     "assert calexp is not None\n",
-    "# Source table for this exposure\n",
+    "# 'physical_filter': 'HSC-R'\n",
+    "f\"Visit: {calexp.visitInfo.getExposureId()}, Band:{calexp.getFilterLabel().bandLabel}, Detector: {calexp.detector.getId()}\""
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Source table for this calexp exposure\n",
     "calexpSrc = butler.get('src', **calexpId)"
    ]
   },
@@ -179,7 +188,15 @@
     "# Define a deep coadded image and retrieve it via the Butler\n",
     "coaddId = {'tract': 4226, 'patch': 17, 'band': 'r'}\n",
     "coadd = butler.get('deepCoadd', **coaddId)\n",
-    "assert coadd is not None\n",
+    "assert coadd is not None"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
     "# Source table for this coadd\n",
     "coaddSrc = butler.get('deepCoadd_forced_src', coaddId)"
    ]
@@ -257,7 +274,7 @@
    "outputs": [],
    "source": [
     "# Define some default plot options for the Image\n",
-    "img_opts = dict(height=400, width=500, \n",
+    "img_opts = dict(height=600, width=700, \n",
     "                xaxis=\"bottom\", \n",
     "                padding = 0.01, fontsize={'title': '8pt'},\n",
     "                colorbar=True, toolbar='right', show_grid=True,\n",

08b_Interactive_Catalog_Visualization.ipynb

@@ -10,7 +10,7 @@
     "Last verified to run on <b>2021-10-05</b> with LSST Science Pipelines release <b>w_2021_33</b> <br>\n",
     "Contact authors: Leanne Guy <br>\n",
     "Target audience: All DP0 delegates. <br>\n",
-    "Container Size: medium or large <br>\n",
+    "Minimum Container Size: medium <br>\n",
     "Questions welcome at <a href=\"https://community.lsst.org/c/support/dp0\">community.lsst.org/c/support/dp0</a> <br>\n",
     "Find DP0 documentation and resources at <a href=\"https://dp0-1.lsst.io\">dp0-1.lsst.io</a> <br>"
    ]
@@ -71,7 +71,7 @@
     "from astropy.units.quantity import Quantity\n",
     "\n",
     "# LSST imports\n",
-    "from rubin_jupyter_utils.lab.notebook import get_tap_service\n",
+    "from lsst.rsp import get_tap_service\n",
     "\n",
     "# Bokeh and Holoviews for visualization\n",
     "import bokeh\n",
@@ -123,10 +123,9 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# Prevent some helpful but ancillary warning messages from printing\n",
-    "# during some LSST DM Release calls\n",
-    "warnings.simplefilter(\"ignore\", category=FutureWarning)\n",
-    "warnings.simplefilter(\"ignore\", category=UserWarning)"
+    "# Ignore warnings\n",
+    "import warnings\n",
+    "warnings.simplefilter(\"ignore\", category='UnitsWarning')"
    ]
   },
   {
@@ -269,7 +268,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "assert len(data) == 23155236"
+    "# assert len(data) == 23155236"
    ]
   },
   {
@@ -416,8 +415,6 @@
     "\n",
     "# Construct a layout using the `+` operator\n",
     "skyplots = skyplot + \\\n",
-    "           skyplot[52.75:60.07, -42.71:-38.95].opts(\n",
-    "    title=\"Skyplot region\", tools=[]) + \\\n",
     "           ra_distribution.options(height=350, width=350)"
    ]
   },
@@ -434,7 +431,7 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
-    "A few things to note about these three plots. The first two are the same plot, with the second one plotting only a subset of the data. When you zooom in on either of these two plots, they will both change in the same manner. The third plot is a different plot object and is not linked to the other two. It will not change in response to actions on either of the first two. Try zooming in on the RA distribution plot, you will notice that the data are not rebinned and that the two skyplots do not change. This is because the plots are not linked. We will see how to link plots in Section 3."
+    "Note that these two plots are not linked, they are two independent plots layed out next to each other. Try zooming in on the skyplot, you will notice that the data are not rebinned in the RA Distribution plot.  We will see how to link plots in Section 3. The tools however do apply to both plots. Try modifying both plots and the used the `reset` tool. Notice that both plots are reset to their original states. "
    ]
   },
   {
@@ -564,7 +561,7 @@
    "outputs": [],
    "source": [
     "%%opts Scatter [tools=['hover'], toolbar='above',height=400, width=400](color='grey')\n",
-    "hv.Scatter(data20K).opts(alpha=0.2, size=2)"
+    "hv.Scatter(data20K).opts(alpha=0.3, size=3)"
    ]
   },
   {
@@ -898,7 +895,6 @@
     "p = dynspread(datashade(points, cmap=\"Viridis\"))\n",
     "p = p.opts(width=800, height=300,\n",
     "    padding=0.05, show_grid=True,\n",
-    "    #xlim=(-2.0, 3.0), ylim=(-2.0, 3.0),\n",
     "    xlim=(-8.0, 18.0), ylim=(-18.0, 8.0),\n",
     "    xlabel=\"(g-r)\", ylabel=\"(r-i)\")"
    ]
@@ -1036,7 +1032,7 @@
     "    return np.log(x) if x > 0 else 0\n",
     "\n",
     "\n",
-    "def update_histogram(bounds=bounds) ->hv.Histogram:\n",
+    "def update_histogram(bounds=bounds) -> hv.Histogram:\n",
     "    selection = (ds.data['ra'] > bounds[0]) & \\\n",
     "                (ds.data['dec'] > bounds[1]) & \\\n",
     "                (ds.data['ra'] < bounds[2]) & \\\n",
@@ -1087,7 +1083,7 @@
     "\n",
     " 1. Holoviews works with a wide range of plotting libraries, Bokeh, matplotlib, plotly, mpld3, pygal to name a few. As an exercise, try changing the Holoviews plotting library to be `matplotlib` instead of `bokeh` in the `Setup` cell at the beginning of the notebook with `hv.extension('matplotlib')`. You will see the holoviews + matplotlib icons displayed when the library is loaded successfully. Run the cells in section 2.1 again and compare the outputs. Try again with some other plotting library. Don't forget to set the plotting library back to whichever you prefer to use for the rest of this tutorial.\n",
     "  \n",
-    " 2. Try writing a different callback function to use in section 6.2. "
+    " 2. Try writing a different callback function to use in section 5.2. "
    ]
   }
  ],