Skip to content

SP-2689: Update NB that make big deep cutouts #84

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
plazas merged 2 commits into from
Oct 16, 2025
Merged

SP-2689: Update NB that make big deep cutouts #84

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
0603e36
Find bcg tract and patch automatically
plazas Oct 15, 2025
fdf2af5
Update 103_7 notebook
plazas Oct 15, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
33 changes: 19 additions & 14 deletions DP1/100_How_to_Use_RSP_Tools/103_Image_access_and_display/103_6_Color_Composite_Image.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,7 +22,7 @@
"Data Release: <a href=\"https://dp1.lsst.io\">Data Preview 1</a> <br>\n",
"Container Size: large <br>\n",
"LSST Science Pipelines version: r29.2.0 <br>\n",
"Last verified to run: 2025-09-19 <br>\n",
"Last verified to run: 2025-10-15 <br>\n",
"Repository: <a href=\"https://github.com/lsst/tutorial-notebooks\">github.com/lsst/tutorial-notebooks</a> <br>"
]
},
Expand Down Expand Up @@ -816,42 +816,47 @@
"metadata": {},
"source": [
"Use the BCG location to obtain the pixel scale.\n",
"Get the WCS of the patch and tract where the BCG is located.\n",
"From the 300 DP1 tutorials, the patch and tract numbers are 10463 and 62, respectively,"
"First, get the WCS of the patch and tract where the BCG is located."
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "8ccd5259-ffb1-45a9-bfc8-79c110301c48",
"id": "d92642e2-e943-4fb5-a9bb-a8eb138862e6",
"metadata": {},
"outputs": [],
"source": [
"bcg_data_ref = next(\n",
" ref for ref in tract_container\n",
" if ref.dataId[\"tract\"] == 10463 and ref.dataId[\"patch\"] == 62\n",
")"
"bcg_point = SpherePoint(ra_bcg * degrees,\n",
" dec_bcg * degrees)\n",
"bcg_tract = skymap.findTract(bcg_point)\n",
"bcg_patch = bcg_tract.findPatch(bcg_point)\n",
"\n",
"bcg_tract_number = bcg_tract.tract_id\n",
"bcg_patch_number = bcg_patch.getSequentialIndex()\n",
"print(bcg_tract_number, bcg_patch_number)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "85859dbc-faf0-4fa4-a5f6-7a6974fd2b46",
"id": "8ccd5259-ffb1-45a9-bfc8-79c110301c48",
"metadata": {},
"outputs": [],
"source": [
"wcs_bcg = butler.get(bcg_data_ref.makeComponentRef('wcs'))"
"bcg_data_ref = next(\n",
" ref for ref in tract_container\n",
" if ref.dataId[\"tract\"] == bcg_tract_number and ref.dataId[\"patch\"] == bcg_patch_number\n",
")"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "133c1d9d-073b-42e3-8a2f-b782f8a9e691",
"id": "85859dbc-faf0-4fa4-a5f6-7a6974fd2b46",
"metadata": {},
"outputs": [],
"source": [
"bcg_point = SpherePoint(ra_bcg * degrees,\n",
" dec_bcg * degrees)"
"wcs_bcg = butler.get(bcg_data_ref.makeComponentRef('wcs'))"
]
},
{
Expand Down Expand Up @@ -897,7 +902,7 @@
"id": "27dbdb21-0f1a-4c32-94d7-c79f174e7ad5",
"metadata": {},
"source": [
"Set the WCS projection. The WCS of deep_coadd images are TAN projections ([see the DP1 documentation](https://dp1.lsst.io/products/images/deep_coadd.html#wcs))."
"Set the WCS projection. The WCS of `deep_coadd` images are TAN projections ([see the DP1 documentation](https://dp1.lsst.io/products/images/deep_coadd.html#wcs))."
]
},
{
Expand Down
24 changes: 20 additions & 4 deletions DP1/100_How_to_Use_RSP_Tools/103_Image_access_and_display/103_7_Big_deep_coadd_cutout.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,7 +22,7 @@
"Data Release: <a href=\"https://dp1.lsst.io\">Data Preview 1</a> <br>\n",
"Container Size: large <br>\n",
"LSST Science Pipelines version: r29.2.0 <br>\n",
"Last verified to run: 2025-09-22 <br>\n",
"Last verified to run: 2025-10-15 <br>\n",
"Repository: <a href=\"https://github.com/lsst/tutorial-notebooks\">github.com/lsst/tutorial-notebooks</a> <br>"
]
},
Expand Down Expand Up @@ -344,8 +344,24 @@
"metadata": {},
"source": [
"Use the BCG location to obtain the pixel scale there.\n",
"First, get the WCS of the patch and tract where the BCG is located.\n",
"From the 300 DP1 tutorials, the patch and tract numbers are 10463 and 62 respectively."
"First, get the WCS of the patch and tract where the BCG is located."
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "ae116f76-2aca-447f-833f-75bc36a7ebab",
"metadata": {},
"outputs": [],
"source": [
"bcg_point = SpherePoint(ra_bcg * degrees,\n",
" dec_bcg * degrees)\n",
"bcg_tract = skymap.findTract(bcg_point)\n",
"bcg_patch = bcg_tract.findPatch(bcg_point)\n",
"\n",
"bcg_tract_number = bcg_tract.tract_id\n",
"bcg_patch_number = bcg_patch.getSequentialIndex()\n",
"print(bcg_tract_number, bcg_patch_number)"
]
},
{
Expand All @@ -357,7 +373,7 @@
"source": [
"bcg_data_ref = next(\n",
" ref for ref in tract_container\n",
" if ref.dataId[\"tract\"] == 10463 and ref.dataId[\"patch\"] == 62\n",
" if ref.dataId[\"tract\"] == bcg_tract_number and ref.dataId[\"patch\"] == bcg_patch_number\n",
")"
]
},
Expand Down