updating stamps nb

Author: christinawilliams

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

@@ -225,7 +225,7 @@
    "id": "114a1c30-b0f0-4bb8-a867-1cbec46f63b9",
    "metadata": {},
    "source": [
-    "Within that folder, create a sub-folder named `dp02_13a_temp` in which to save files created by this tutorial.\n",
+    "Within that folder, create a sub-folder named `dp1_103_9_temp` in which to save files created by this tutorial.\n",
     "At the end of the notebook, the last cell will clear the contents and remove this temporary folder."
    ]
   },
@@ -236,7 +236,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "tempdir = userdir + '/dp1_103_5_temp'\n",
+    "tempdir = userdir + '/dp1_103_9_temp'\n",
     "if not os.path.exists(tempdir):\n",
     "    os.makedirs(tempdir)\n",
     "    print('Created ', tempdir)\n",
@@ -369,12 +369,11 @@
     "    #axes = axes.flatten() if n_subplots > 1 else [axes]\n",
     "\n",
     "    for ax in axes:\n",
-    "        ax.axis(\"off\")   # ← removes frame, ticks, labels, everything\n",
+    "        ax.axis(\"off\")\n",
     "\n",
     "    for ax in axes[n_subplots:]:\n",
     "        ax.remove()\n",
     "\n",
-    "    \n",
     "    return fig, axes[:n_subplots]"
    ]
   },
@@ -706,7 +705,8 @@
     "ax2.coords['dec'].set_axislabel('Declination')\n",
     "ax2.coords['dec'].set_axislabel_position('r')\n",
     "ax2.coords['dec'].set_ticklabel_position('r')\n",
-    "ax2.set_title('Template Image')\n"
+    "ax2.set_title('Template Image')\n",
+    "plt.show()"
    ]
   },
   {
@@ -734,28 +734,28 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "num_visits = len(scitab['t_max'])\n",
-    "#print('number of visits:', num_visits)\n",
-    "#for i in range(num_visits):\n",
-    "#    plt.title(f\"MJD {np.round(scitab['t_max'][i], 4)}\")\n",
+    "num_visits = 20 #len(scitab['t_max'])\n",
+    "print('number of visits:', num_visits)\n",
+    "for i in range(num_visits):\n",
+    "    plt.title(f\"MJD {np.round(scitab['t_max'][i], 4)}\")\n",
     "\n",
-    "#    dl_result_sci = DatalinkResults.from_result_url(scitab['access_url'][i],\n",
-    "#                                                    session=get_pyvo_auth())\n",
+    "    dl_result_sci = DatalinkResults.from_result_url(scitab['access_url'][i],\n",
+    "                                                    session=get_pyvo_auth())\n",
     "\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",
-    "#    sciwcs = WCS(scihdul[1].header)\n",
-    "#    reproj_scidata, footprint = reproject_interp((scidata, sciwcs), refwcs)\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",
+    "    sciwcs = WCS(scihdul[1].header)\n",
+    "    reproj_scidata, footprint = reproject_interp((scidata, sciwcs), refwcs)\n",
     "\n",
-    "#    plt.imshow(reproj_scidata, origin='lower', vmin=np.nanpercentile(reproj_scidata, 1),\n",
-    "#               vmax=np.nanpercentile(reproj_scidata, 99))\n",
+    "    plt.imshow(reproj_scidata, origin='lower', vmin=np.nanpercentile(reproj_scidata, 1),\n",
+    "               vmax=np.nanpercentile(reproj_scidata, 99))\n",
     "\n",
-    "#    figname = os.path.join(tempdir, 'cutout_' + str(i) + '.png')\n",
-    "#    if os.path.isfile(figname):\n",
-    "#        os.remove(figname)\n",
-    "#    plt.savefig(figname)\n"
+    "    figname = os.path.join(tempdir, 'cutout_' + str(i) + '.png')\n",
+    "    if os.path.isfile(figname):\n",
+    "        os.remove(figname)\n",
+    "    plt.savefig(figname)\n"
    ]
   },
   {
@@ -773,37 +773,15 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "fig, axes = make_subplot_grid(100)\n",
+    "fig, axes = make_subplot_grid(num_visits)\n",
     "for i, ax in enumerate(axes):\n",
     "\n",
-    "    j = i #+ 5\n",
+    "    j = i\n",
     "    figname = os.path.join(tempdir, 'cutout_' + str(j) + '.png')\n",
     "    img = mpimg.imread(figname)\n",
     "\n",
-    "# 2. Display the image\n",
-    "# imshow interprets the array as an image\n",
     "    ax.imshow(img)\n",
     "\n",
-    "    #ax.set_title(f\"MJD {np.round(scitab['t_max'][j], 4)}\")\n",
-    "\n",
-    "    #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, ra, dec, get_pyvo_auth(), fov)\n",
-    "    #sci_hdu = scihdul[1]\n",
-    "    #sci_header = scihdul[1].header\n",
-    "    #scidata = scihdul[1].data\n",
-    "    #sciwcs = WCS(scihdul[1].header)\n",
-    "    #reproj_scidata, footprint = reproject_interp((scidata, sciwcs), refwcs)\n",
-    "\n",
-    "    #ax.imshow(reproj_scidata, origin='lower', vmin=np.nanpercentile(reproj_scidata, 1),\n",
-    "    #          vmax=np.nanpercentile(reproj_scidata, 99))\n",
-    "\n",
-    "    #figname = os.path.join(tempdir, 'cutout_' + str(i) + '.png')\n",
-    "    #if os.path.isfile(figname):\n",
-    "    #    os.remove(figname)\n",
-    "    #plt.savefig(figname)\n",
-    "\n",
     "plt.tight_layout()\n",
     "plt.show()"
    ]
@@ -813,7 +791,7 @@
    "id": "2b54d1d3-9564-4836-8265-4c8eedc4502a",
    "metadata": {},
    "source": [
-    "> **Figure 3**: Many cutouts of the science visit images of the DIA object, in chronological order, after being reprojected onto the same WCS as the template image cutout (right). Every Nth visit of the 117 existing visits is plotted for brevity.\n"
+    "> **Figure 3**: Many cutouts of the science visit images of the DIA object, in chronological order, after being reprojected onto the same WCS as the template image cutout (right). "
    ]
   },
   {
@@ -829,43 +807,38 @@
    ]
   },
   {
-   "cell_type": "code",
-   "execution_count": null,
+   "cell_type": "markdown",
+   "id": "1c07a8b1-f24a-4c03-a56b-515a381cfa42",
+   "metadata": {},
+   "source": [
+    "> **Figure 4**: A gif animation generated from the reproject visit images in Figure 3."
+   ]
+  },
+  {
+   "cell_type": "markdown",
    "id": "172dc7bb-67ee-49f9-bee1-a6f50ad96c89",
    "metadata": {},
-   "outputs": [],
-   "source": []
+   "source": [
+    "## 4. Future functionality\n",
+    "\n",
+    "when rotate, don't do grid in final version, then get to a releasable version if you explain that the wcs rotate and just produce a rotated image (png, no wcs info) rotated visit image stamp that you can make a png image out of. Say a few things about how we're working on how to overplot grid lines and wcs functionality. this gets us to a publishable thing "
+   ]
   },
   {
    "cell_type": "markdown",
    "id": "a8eb043b-5afe-4990-9050-bd059aaf5abf",
    "metadata": {},
    "source": [
     "\n",
-    "when rotate, don't do grid in final version, then get to a releasable version if you explain that the wcs rotate and just produce a rotated image (png, no wcs info) rotated visit image stamp that you can make a png image out of. Say a few things about how we're working on how to overplot grid lines and wcs functionality. this gets us to a publishable thing \n",
+    "\n",
     "\n",
     "only for epochs where it was detected in visit images. users can extrapolate to get all images. \n",
     "\n",
     "show how to do smaller more constrained part (i.e timeline over which visits are retrieved; where a dia ). \n",
     "\n",
     "\n",
-    " We want to show how to store a bunch of png files. Once you have those on file then you can just use a simple\n",
-    " function to make the gif\n",
-    "\n",
-    "\n",
-    " note: store png files to scratch directory\n",
-    " cutouts could be even smaller than this.\n",
-    "\n",
     "start with an ra/dec and date range where you already know. don't even need to refer to the other notebook\n"
    ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "1c07a8b1-f24a-4c03-a56b-515a381cfa42",
-   "metadata": {},
-   "source": [
-    "> **Figure 4**: A gif animation generated from the reproject visit images in Figure 3."
-   ]
   }
  ],
  "metadata": {