Merge pull request #315 from bhilbert4/imviz-demo-new-data

Imviz demo notebook: Use new data to replace data deleted from MAST

Author: bhilbert4

notebooks/jdaviz_demo/notebooks/Imviz.ipynb

@@ -15,7 +15,7 @@
    "source": [
     "# Jdaviz Imviz Demo\n",
     "**Author: Clare Shanahan, Space Telescope Science Institute**<br>\n",
-    "**Last update: May 27, 2025**\n",
+    "**Last update: Oct 30, 2025**\n",
     "\n",
     "## Tutorial Overview\n",
     "This tutorial will demonstrate an example workflow to show of some key features of [Imviz](https://jdaviz.readthedocs.io/en/latest/imviz/index.html). We will be looking at some basic functionality (starting the app in a notebook, loading data, pan/zoom, subsets/regions of interest) as well as some of the higher level plugin tools including catalog search and aperture photometry.\n",
@@ -73,14 +73,13 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "filenames = ['jw01227-c1002_t005_nircam_clear-f335m_i2d.fits',\n",
-    "             'jw01227-c1002_t005_nircam_clear-f277w_i2d.fits',\n",
-    "             'jw01227-c1002_t005_nircam_clear-f444w_i2d.fits']\n",
-    "\n",
+    "filenames = ['jw01227025001_04101_00001_nrcalong_i2d.fits',  # F277W\n",
+    "             'jw01227025001_02101_00001_nrcalong_i2d.fits',  # F335M\n",
+    "             'jw01227025001_06101_00001_nrcalong_i2d.fits']  # F444W\n",
     "with imviz.batch_load():  # not necessary, but this context manager makes loading multiple files more efficient          \n",
     "    for filename in filenames:\n",
-    "        imviz.load_data(f'mast:JWST/product/{filename}', cache=True)  # to re-download from MAST (or use cached files)\n",
-    "        # imviz.load_data(f'./data/{filename}', cache=True)  # pre-downloaded data"
+    "        imviz.load(f'mast:JWST/product/{filename}', cache=True)  # to re-download from MAST (or use cached files)\n",
+    "        # imviz.load(f'./data/{filename}', cache=True)  # pre-downloaded data"
    ]
   },
   {
@@ -145,7 +144,7 @@
     "plot_options.stretch_preset = '99.5%'\n",
     "\n",
     "# increase vmax to a more suitable value\n",
-    "plot_options.stretch_vmax = 600"
+    "plot_options.stretch_vmax = 6"
    ]
   },
   {
@@ -193,8 +192,8 @@
     "# select Gaia catalog\n",
     "catalogs_plugin.catalog.selected = 'Gaia'\n",
     "\n",
-    "# request only 10 sources\n",
-    "catalogs_plugin.max_sources = 10\n",
+    "# Set an upper limit on the number of sources\n",
+    "catalogs_plugin.max_sources = 500\n",
     "\n",
     "# and run the search\n",
     "catalogs_plugin.search()\n",
@@ -213,13 +212,30 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# now select just the first point\n",
-    "catalogs_plugin.table.selected_rows = catalogs_plugin.table.items[0:1]\n",
+    "# Pick out one source as an example- we will use the source at index 455\n",
+    "idx = 455\n",
+    "\n",
+    "# now select just the matching source\n",
+    "catalogs_plugin.table.selected_rows = catalogs_plugin.table.items[idx: idx+1]\n",
     "\n",
     "# and zoom to that point\n",
     "catalogs_plugin.zoom_to_selected()"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "338e9544-f480-4ab7-8751-69316ae8763d",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Get the RA and Dec of the selected source\n",
+    "ra = float(catalogs_plugin.table.items[idx]['Right Ascension (degrees)'])\n",
+    "dec = float(catalogs_plugin.table.items[idx]['Declination (degrees)'])\n",
+    "print(f\"RA: {ra}\")\n",
+    "print(f\"Dec: {dec}\")"
+   ]
+  },
   {
    "cell_type": "markdown",
    "id": "06025284-99ea-4cd6-bf86-a6e14ffab8f4",
@@ -248,7 +264,8 @@
     "# just for demo sake, shift the coordinates a tiny bit from the catalog position so we\n",
     "# can use the 'recenter' position. This recreates the scenario of freehand drawing a circular subset\n",
     "# rather than placing it at an exact location\n",
-    "circular_region = CircleSkyRegion(center=SkyCoord(14.77039, -72.16949, unit='degree'), radius=0.0001*u.deg)\n",
+    "shift = 3. * 0.062 / 3600.  # shift by 3 pixels\n",
+    "circular_region = CircleSkyRegion(center=SkyCoord(ra + shift, dec, unit='degree'), radius=0.0001*u.deg)\n",
     "subset_plugin.import_region(circular_region)"
    ]
   },
@@ -342,7 +359,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.12.10"
+   "version": "3.11.14"
   }
  },
  "nbformat": 4,