update deblender nb Fred comments

Author: christinawilliams

DP1/200_Data_Products/206_Deblender_Products/206_2_Deblender_Footprints.ipynb

@@ -697,7 +697,9 @@
    "id": "feb57f86-ccdd-4d84-8546-1e48f644654c",
    "metadata": {},
    "source": [
-    "Load the subset of the image that overlaps with the blend, then use the `blend_data` and the `observation` to create a full `Blend`. All of the model information is the same as the information from the `Blend` that we created earlier. Here, attach real `Observation` data as opposed to an empty `Observation` (as was done in Section 3.2). Notice that the **observed_bands** are passed here, not the **display_bands**, as these bands must match the order of the bands in the `blend`."
+    "Load the subset of the image that overlaps with the blend, then use the `blend_data` and the `observation` to create a full `Blend`. All of the model information is the same as the information from the `Blend` that we created earlier. Here, attach real `Observation` data as opposed to an empty `Observation` (as was done in Section 3.2). Notice that the **observed_bands** are passed here, not the **display_bands**, as these bands must match the order of the bands in the `blend`.\n",
+    "\n",
+    "Model residuals compared to the observations are of interest for exploring how well the blend model matches the data. To do this, first load a multiband exposure as the input observation. Use the patch and tract to reconstruct the multiband deep_coadd for this blend.\n"
    ]
   },
   {
@@ -761,7 +763,7 @@
     "\n",
     "To normalize the image use the Lupton red green blue (RGB) AsinhMapping, which preserves the observed colors.\n",
     "\n",
-    "Note: You'll notice that the colors match the processed color order, not the correct band order. In the future passing `blend[display_bands]` should display everything in the correct order."
+    "Note: You'll notice that the colors match the processed color order, not the correct band order. In the future passing `blend[display_bands]` should display everything in the correct order.\n"
    ]
   },
   {
@@ -821,7 +823,7 @@
    "id": "17e2781f-5f28-4e44-8a6b-550d320c2941",
    "metadata": {},
    "source": [
-    "> Figure 4: Visualization of the individual deblended sources (top to bottom show ID 0 to 8). For each source, left panel shows the PSF-deconvolved model footprint (where the color scale corresponds to the weighted RGB color), middle panel shows the observation of the blend with the source identified by a white box, and right panel shows the flux in each band, from u to y in order of wavelength as defined in Section 3.2."
+    "> Figure 4: Visualization of the individual deblended sources (top to bottom show ID 0 to 8). For each source, left panel shows the PSF-deconvolved model footprint (where the color scale corresponds to the weighted RGB color), middle panel shows the observation of the blend with the source identified by a white box, and right panel shows the flux in each band, from u to y in order of wavelength as defined in Section 3.2. TBD IF THIS WAVELENGTH STATEMENT IS TRUE?"
    ]
   },
   {
@@ -952,140 +954,7 @@
    "id": "cb5c74a8-b1cb-4e89-9ec3-182bbe880900",
    "metadata": {},
    "source": [
-    "> Figure 5:"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "711aaf6f-500d-4215-a6f8-bce00807c878",
-   "metadata": {},
-   "source": [
-    "FRED DELETED EVERYTHING BELOW THIS:\n",
-    "\n",
-    "### 3.3 Blend Model Residuals\n",
-    "\n",
-    "Model residuals compared to the observations are of interest for exploring how well the blend model matches the data. To do this, first load a multiband exposure as the input observation. \n",
-    "\n",
-    "The cell below prepares to build the multiband exposure by first finding the patch and tract that covers the blend of interest."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "8dafe145-8091-422a-9071-263ea0c68e7c",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "skymap = butler.get(\"skyMap\")\n",
-    "radec = geom.SpherePoint(children['coord_ra'][0], children['coord_dec'][0], geom.degrees)\n",
-    "\n",
-    "tractInfo = skymap.findTract(radec)\n",
-    "patchInfo = tractInfo.findPatch(radec)\n",
-    "patch = tractInfo.getSequentialPatchIndex(patchInfo)\n",
-    "tract = tractInfo.getId()\n"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "f2049e77-abd7-40c7-89de-4df0c4a16a0c",
-   "metadata": {},
-   "source": [
-    "Next, use the patch and tract to reconstruct the multiband deep_coadd for this blend, using the bands that were observed."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "72de1035-7f29-469c-9650-2a62e9f82ecf",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "mCoadd = MultibandExposure.fromButler(butler, display_bands,\n",
-    "                                      \"deep_coadd\", tract=tract, patch=patch)\n",
-    "bbox = sl.Box(blend_data.shape, blend_data.origin)\n",
-    "bbox = mes.utils.scarletBoxToBBox(bbox)\n",
-    "\n",
-    "observation = mes.utils.buildObservation(\n",
-    "    modelPsf=model_psf,\n",
-    "    psfCenter=bbox.getCenter(),\n",
-    "    mExposure=mCoadd[:, bbox],\n",
-    ")\n",
-    "\n",
-    "blend = blend_data.to_blend(observation)\n"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "5f16b7ba-1b0b-47e3-a207-ac4d334fa64f",
-   "metadata": {},
-   "source": [
-    "Then, display the scene. The function scarlet `show_scene` will include the model footprints, the model rendered (i.e. convolved with the PSF), the locations on the observed multiband image, and residuals between the real image and the model.\n",
-    "\n",
-    "To normalize the image use the Lupton red green blue (RGB) `AsinhMapping`, which preserves the observed colors. "
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "4ad78ea8-51e4-4956-94c4-7e02a4c5de2d",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "norm = sl.display.AsinhMapping(minimum=0, stretch=1, Q=10)\n",
-    "\n",
-    "sl.display.show_scene(\n",
-    "    blend,\n",
-    "    norm=norm,\n",
-    "    show_model=True,\n",
-    "    show_rendered=True,\n",
-    "    show_observed=True,\n",
-    "    show_residual=True,\n",
-    "    linear=False,\n",
-    "    add_boxes=True,\n",
-    "    add_labels=True,\n",
-    "    figsize=(10, 7),\n",
-    ")\n",
-    "\n",
-    "plt.show()"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "3f5800d3-b018-4685-8555-fcf3e696c16a",
-   "metadata": {},
-   "source": [
-    "> Figure 3: Visualization of the model generated by Scarlet Lite, with x and y axes indicating pixel number from the patch/tract. The panels show the blend model (not PSF-convolved; top left), the blend model rendered (PSF-convolved; top right), the observation (bottom left) and the residual image (observation - model rendered; bottom right). All deblended sources are identified with white boxes and a number indicating its ID within the set of deblended objects."
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "79111d63-105f-408a-ab49-3118c67b2096",
-   "metadata": {},
-   "source": [
-    "Finally, use the `show_sources` function to display each individual deblended object's model (left), its location in the image (middle), and its spectral energy distribution (right)."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "d0ced1f4-3b6c-4bc6-a2dc-80018c3f15f4",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "sl.display.show_sources(blend,\n",
-    "                        norm=norm,\n",
-    "                        show_observed=True,\n",
-    "                        add_boxes=True\n",
-    "                        )\n",
-    "plt.show()"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "895b605f-c964-4dbe-99a3-bdf2bdcec2cc",
-   "metadata": {},
-   "source": [
-    "> Figure 4: Visualization of the individual deblended sources (top to bottom show ID 0 to 8). For each source, left panel shows the PSF-deconvolved model footprint (where the color scale corresponds to the weighted RGB color), middle panel shows the observation of the blend with the source identified by a white box, and right panel shows the flux in each band, from u to y in order of wavelength as defined in Section 3.2."
+    "> Figure 5: Visualization of deblended source 1. Left panel shows the PSF-deconvolved model footprint (where the color scale corresponds to the weighted RGB color), middle panel shows the convolution of the scarlet model with the difference kernel and right panel shows the flux-redistributed model that the science pipelines use for measurements."
    ]
   }
  ],