adding new cutout sync section 4 to image cutout nb

Author: christinawilliams

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

@@ -683,7 +683,11 @@
    "source": [
     "## 4. Other cutout service data types\n",
     "\n",
-    "As mentioned in Section 3, there are three types of cutout services: cutout-sync-exposure to return an `ExposureF` type image with all LSST image extensions (demonstrated above); cutout-sync which returns just the image extension, and and cutout-sync-mask which contains just the mask extension. The latter 2 will be demonstrated below.\n"
+    "As mentioned in Section 3, there are three types of cutout services: `cutout-sync-exposure` to return an `ExposureF` type image with all LSST image extensions (demonstrated above); `cutout-sync` which returns just the image extension, and and `cutout-sync-mask` which contains just the mask extension. The latter 2 will be demonstrated below.\n",
+    "\n",
+    "### 4.1 Single extension image cutouts\n",
+    "\n",
+    "Below demonstrates how to retrieve single extension images using `cutout-sync`."
    ]
   },
   {
@@ -787,6 +791,99 @@
     "plt.show()"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "918e8779-2939-4978-9ec4-c8a4047f04e9",
+   "metadata": {},
+   "source": [
+    "### 4.2 Single extension mask cutouts\n",
+    "\n",
+    "Below demonstrates how to retrieve the corresponding mask using `cutout-sync-mask`."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "a1213ce5-9049-43ce-8032-f1ea2db9c576",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "\n",
+    "sq5 = SodaQuery.from_resource(dl_result,\n",
+    "                             dl_result.get_adhocservice_by_id(\"cutout-sync-maskedimage\"),\n",
+    "                             session=get_pyvo_auth())\n",
+    "\n",
+    "spherePoint = geom.SpherePoint(target_ra*geom.degrees, target_dec*geom.degrees)\n",
+    "Radius = 0.01 * u.deg\n",
+    "sq5.circle = (spherePoint.getRa().asDegrees() * u.deg,\n",
+    "             spherePoint.getDec().asDegrees() * u.deg,\n",
+    "             Radius)\n",
+    "\n",
+    "cutout_bytes = sq5.execute_stream().read()\n",
+    "sq5.raise_if_error()\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "a2f5ed28-c17d-451b-969c-b7d9186eeb23",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "mem = MemFileManager(len(cutout_bytes))\n",
+    "mem.setData(cutout_bytes, len(cutout_bytes))\n",
+    "\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "770c6bc2-67f5-4d3f-9d3a-98cd27e28b35",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "data = mem.getData()   # returns `bytes`\n",
+    "\n",
+    "# this doesn't work\n",
+    "import lsst.afw.image as afwImage\n",
+    "\n",
+    "#image = afwImage.ImageF(mem)\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "ba884db3-cd3e-45b9-bba0-ecd98cd95358",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "\n",
+    "#data = mem.getData()\n",
+    "hdul = fits.open(io.BytesIO(data))\n",
+    "hdul.info()\n",
+    "\n",
+    "image_hdu = hdul[2]        # IMAGE extension\n",
+    "image_data = image_hdu.data  # NumPy array\n",
+    "#image = afwImage.ImageF(image_data)\n",
+    "image = afwImage.ImageF(image_data)#.astype(\"float32\"))\n",
+    "print(type(image_data))   # numpy.ndarray\n",
+    "print(type(image))        # lsst.afw.image._image.ImageF\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "cceab1fb-7f66-4858-b46f-d71fddb50159",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "\n",
+    "display = afwDisplay.Display()\n",
+    "display.scale('asinh', 'zscale')\n",
+    "display.image(image)\n",
+    "plt.show()"
+   ]
+  },
   {
    "attachments": {
     "923e951d-30c0-4a8f-a1b0-8bb851f1cb5f.png": {