|
22 | 22 | "Data Release: <a href=\"https://dp1.lsst.io/\">Data Preview 1</a>\\\n", |
23 | 23 | "Container Size: medium\\\n", |
24 | 24 | "LSST Science Pipelines version: v29.2.0\\\n", |
25 | | - "Last verified to run: 2026-02-10\\\n", |
| 25 | + "Last verified to run: 2026-03-24\\\n", |
26 | 26 | "Repository: [github.com/lsst/tutorial-notebooks](https://github.com/lsst/tutorial-notebooks)\\\n", |
27 | 27 | "DOI: [10.11578/rubin/dc.20250909.20](https://doi.org/10.11578/rubin/dc.20250909.20)" |
28 | 28 | ] |
|
63 | 63 | "\n", |
64 | 64 | "Two methods for transforming between photometric systems are explored in <a href=\"https://rtn-099.lsst.io/\">RTN-099</a>: polynomial relations and lookup tables. The guiding principle for which is to use is driven by the balance between simplicity and accuracy: in general, simple first- or second-order polynomials are sufficient, but, for higher accuracy, one may instead choose the slightly more difficult-to-employ lookup-table method.\n", |
65 | 65 | "\n", |
66 | | - "This tutorial demonstrates how to apply both methods of photometric transformation. Examples are provided for converting from Rubin DP1 $i$-band photometry to PanSTARRS1 DR2 and DES DR2 $i$ bands, and from Rubin DP1 $g$-band photometry to Gaia DR3 $RP$-band photometry.\n", |
| 66 | + "This tutorial demonstrates how to apply both methods of photometric transformation. Examples are provided for converting from Rubin DP1 $i$-band photometry to PanSTARRS1 DR2 and DES DR2 $i$ bands, and from Rubin DP1 $g$-band photometry to Gaia DR3 $BP$-band photometry.\n", |
67 | 67 | "\n", |
68 | 68 | "**Related tutorials:** The 100-level tutorials demonstrate how to use the butler, the TAP service, and the Firefly image display. The 200-level tutorials introduce the types of image and catalog data.\n" |
69 | 69 | ] |
|
336 | 336 | { |
337 | 337 | "cell_type": "markdown", |
338 | 338 | "id": "686c2789-0b71-4bc4-b15b-0b020a0d965f", |
339 | | - "metadata": {}, |
| 339 | + "metadata": { |
| 340 | + "execution": { |
| 341 | + "iopub.execute_input": "2025-10-17T18:50:22.291274Z", |
| 342 | + "iopub.status.busy": "2025-10-17T18:50:22.291002Z", |
| 343 | + "iopub.status.idle": "2025-10-17T18:50:22.295488Z", |
| 344 | + "shell.execute_reply": "2025-10-17T18:50:22.294633Z", |
| 345 | + "shell.execute_reply.started": "2025-10-17T18:50:22.291257Z" |
| 346 | + } |
| 347 | + }, |
340 | 348 | "source": [ |
341 | 349 | "Calculate color indices using the PSF magnitudes, and plot a color-color diagram. Color-code points by the $r$-band PSF magnitude error." |
342 | 350 | ] |
|
410 | 418 | { |
411 | 419 | "cell_type": "markdown", |
412 | 420 | "id": "e3e1fe82-cea3-4cc2-81d6-5e9310d2ced1", |
413 | | - "metadata": {}, |
| 421 | + "metadata": { |
| 422 | + "execution": { |
| 423 | + "iopub.execute_input": "2025-12-05T20:06:47.537387Z", |
| 424 | + "iopub.status.busy": "2025-12-05T20:06:47.537024Z", |
| 425 | + "iopub.status.idle": "2025-12-05T20:06:47.542498Z", |
| 426 | + "shell.execute_reply": "2025-12-05T20:06:47.541537Z", |
| 427 | + "shell.execute_reply.started": "2025-12-05T20:06:47.537361Z" |
| 428 | + } |
| 429 | + }, |
414 | 430 | "source": [ |
415 | 431 | "\n", |
416 | 432 | "| Conversion | Transformation Equation | RMS | Applicable Color Range | QA Plot |\n", |
|
548 | 564 | { |
549 | 565 | "cell_type": "markdown", |
550 | 566 | "id": "c4f941a7-2394-4dd9-b2cf-ce5fb16f8395", |
551 | | - "metadata": {}, |
| 567 | + "metadata": { |
| 568 | + "execution": { |
| 569 | + "iopub.execute_input": "2025-12-19T23:32:23.571408Z", |
| 570 | + "iopub.status.busy": "2025-12-19T23:32:23.571024Z", |
| 571 | + "iopub.status.idle": "2025-12-19T23:32:23.576885Z", |
| 572 | + "shell.execute_reply": "2025-12-19T23:32:23.575789Z", |
| 573 | + "shell.execute_reply.started": "2025-12-19T23:32:23.571384Z" |
| 574 | + } |
| 575 | + }, |
552 | 576 | "source": [ |
553 | 577 | "Use the `cross_match_catalogs` function defined above to match the two catalogs based on the stars' coordinates." |
554 | 578 | ] |
|
569 | 593 | { |
570 | 594 | "cell_type": "markdown", |
571 | 595 | "id": "8623e190-aa85-4dd0-b186-88e19483743f", |
572 | | - "metadata": {}, |
| 596 | + "metadata": { |
| 597 | + "execution": { |
| 598 | + "iopub.execute_input": "2025-12-16T19:45:53.129870Z", |
| 599 | + "iopub.status.busy": "2025-12-16T19:45:53.129489Z", |
| 600 | + "iopub.status.idle": "2025-12-16T19:45:53.133314Z", |
| 601 | + "shell.execute_reply": "2025-12-16T19:45:53.132627Z", |
| 602 | + "shell.execute_reply.started": "2025-12-16T19:45:53.129846Z" |
| 603 | + } |
| 604 | + }, |
573 | 605 | "source": [ |
574 | 606 | "### 3.3.3 Compare the LSSTComCam-estimated PS1 i-band magnitudes with the PanSTARRS1 DR2 i-band magnitudes" |
575 | 607 | ] |
|
688 | 720 | { |
689 | 721 | "cell_type": "markdown", |
690 | 722 | "id": "0570b618-901f-4c60-bd32-17ea8d69c99f", |
691 | | - "metadata": {}, |
| 723 | + "metadata": { |
| 724 | + "execution": { |
| 725 | + "iopub.execute_input": "2025-12-05T19:57:07.257591Z", |
| 726 | + "iopub.status.busy": "2025-12-05T19:57:07.257182Z", |
| 727 | + "iopub.status.idle": "2025-12-05T19:57:07.261246Z", |
| 728 | + "shell.execute_reply": "2025-12-05T19:57:07.260527Z", |
| 729 | + "shell.execute_reply.started": "2025-12-05T19:57:07.257564Z" |
| 730 | + } |
| 731 | + }, |
692 | 732 | "source": [ |
693 | 733 | "## 4. Piecewise 1st-order polynomial transformation: LSSTComCam $i$ --> DES DR2 $i$\n", |
694 | 734 | "\n", |
|
729 | 769 | { |
730 | 770 | "cell_type": "markdown", |
731 | 771 | "id": "35075b43-4e76-4463-a359-def323fa8520", |
732 | | - "metadata": {}, |
| 772 | + "metadata": { |
| 773 | + "execution": { |
| 774 | + "iopub.execute_input": "2025-12-05T20:10:17.123399Z", |
| 775 | + "iopub.status.busy": "2025-12-05T20:10:17.123015Z", |
| 776 | + "iopub.status.idle": "2025-12-05T20:10:17.128630Z", |
| 777 | + "shell.execute_reply": "2025-12-05T20:10:17.127665Z", |
| 778 | + "shell.execute_reply.started": "2025-12-05T20:10:17.123372Z" |
| 779 | + } |
| 780 | + }, |
733 | 781 | "source": [ |
734 | 782 | "\n", |
735 | 783 | "\n", |
|
1027 | 1075 | "id": "939e52ad-9ce0-45eb-a163-bb9ad9467df8", |
1028 | 1076 | "metadata": {}, |
1029 | 1077 | "source": [ |
1030 | | - "Here are the quality assurance (QA) plots from <a href=\"https://rtn-099.lsst.io/\">RTN-099</a> associated with the second-order polynomial fit of the DP1 $g$-band to Gaia DR3 $RP$-band. Note from the bottom-right plot that a third-order term might benefit the fit, but staying with the simpler second-order polynomial still achieves a 0.027 mag RMS in the fit. For more accuracy, using the lookup table method (described in the next section) is recommended." |
| 1078 | + "Here are the quality assurance (QA) plots from <a href=\"https://rtn-099.lsst.io/\">RTN-099</a> associated with the second-order polynomial fit of the DP1 $g$-band to Gaia DR3 $BP$-band. Note from the bottom-right plot that a third-order term might benefit the fit, but staying with the simpler second-order polynomial still achieves a 0.027 mag RMS in the fit. For more accuracy, using the lookup table method (described in the next section) is recommended." |
1031 | 1079 | ] |
1032 | 1080 | }, |
1033 | 1081 | { |
|
1085 | 1133 | "id": "10f38f41-8bec-418e-a778-18246278f130", |
1086 | 1134 | "metadata": {}, |
1087 | 1135 | "source": [ |
1088 | | - "Test the transformation by comparing the Gaia DR3 $RP$-band magnitudes estimated from the DP1 data in ECDFS against the Gaia DR3 $RP$-band magnitudes from the Gaia DR3 data themselves in ECDFS." |
| 1136 | + "Test the transformation by comparing the Gaia DR3 $BP$-band magnitudes estimated from the DP1 data in ECDFS against the Gaia DR3 $BP$-band magnitudes from the Gaia DR3 data themselves in ECDFS." |
1089 | 1137 | ] |
1090 | 1138 | }, |
1091 | 1139 | { |
|
1180 | 1228 | "id": "febf7500-a25f-4209-b337-86108f4f4d6c", |
1181 | 1229 | "metadata": {}, |
1182 | 1230 | "source": [ |
1183 | | - "Calculate both the differences between the DP1 $g$-band magnitude and the Gaia DR3 $RP$-band magnitude (\"before transformation\") and the differences between DP1-estimated Gaia DR3 $RP$-band magnitude and the actual Gaia DR3 $RP$-band magnitude (\"after transformation\"):" |
| 1231 | + "Calculate both the differences between the DP1 $g$-band magnitude and the Gaia DR3 $BP$-band magnitude (\"before transformation\") and the differences between DP1-estimated Gaia DR3 $BP$-band magnitude and the actual Gaia DR3 $BP$-band magnitude (\"after transformation\"):" |
1184 | 1232 | ] |
1185 | 1233 | }, |
1186 | 1234 | { |
|
1324 | 1372 | "id": "3581e642-4b84-4cc9-b898-be00f0846653", |
1325 | 1373 | "metadata": {}, |
1326 | 1374 | "source": [ |
1327 | | - "Here are the quality‑assurance (QA) plots from <a href=\"https://rtn-099.lsst.io/\">RTN-099</a> that document the lookup‑table transformation from the DP1 $g$‑band to the Gaia DR3 $RP$‑band. The lookup table is constructed by binning matched stars in the ECDFS field by their Rubin DP1 $(g–i)$ color, and then tabulating the median magnitude offset between the Gaia DR3 RP‑band and the Rubin DP1 g‑band for the stars in each color bin." |
| 1375 | + "Here are the quality‑assurance (QA) plots from <a href=\"https://rtn-099.lsst.io/\">RTN-099</a> that document the lookup‑table transformation from the DP1 $g$‑band to the Gaia DR3 $BP$‑band. The lookup table is constructed by binning matched stars in the ECDFS field by their Rubin DP1 $(g–i)$ color, and then tabulating the median magnitude offset between the Gaia DR3 RP‑band and the Rubin DP1 g‑band for the stars in each color bin." |
1328 | 1376 | ] |
1329 | 1377 | }, |
1330 | 1378 | { |
|
1375 | 1423 | "cell_type": "code", |
1376 | 1424 | "execution_count": null, |
1377 | 1425 | "id": "d6f731f9-1084-46fc-a881-80a62b241dd8", |
1378 | | - "metadata": {}, |
| 1426 | + "metadata": { |
| 1427 | + "scrolled": true |
| 1428 | + }, |
1379 | 1429 | "outputs": [], |
1380 | 1430 | "source": [ |
1381 | 1431 | "df_lut = pd.read_csv(lut_url)\n", |
|
1439 | 1489 | "id": "376a20fd-2f9c-4005-9782-6d8607d59d72", |
1440 | 1490 | "metadata": {}, |
1441 | 1491 | "source": [ |
1442 | | - "Test the transformation by comparing the Gaia DR3 $RP$-band magnitudes estimated from the DP1 data in ECDFS against the Gaia DR3 $RP$-band magnitudes from the Gaia DR3 data themselves." |
| 1492 | + "Test the transformation by comparing the Gaia DR3 $BP$-band magnitudes estimated from the DP1 data in ECDFS against the Gaia DR3 $BP$-band magnitudes from the Gaia DR3 data themselves." |
1443 | 1493 | ] |
1444 | 1494 | }, |
1445 | 1495 | { |
|
1553 | 1603 | "plt.grid(True)\n", |
1554 | 1604 | "plt.grid(color='grey')\n", |
1555 | 1605 | "plt.legend(loc='upper left', fontsize=8)\n", |
1556 | | - "plt.title(r\"mag difference between LSSTComCam DP1 $g$-band (and transformed $BP$-band) and Gaia DR3 $RP$-band\")\n", |
| 1606 | + "plt.title(r\"mag difference between LSSTComCam DP1 $g$-band (and transformed $BP$-band) and Gaia DR3 $BP$-band\")\n", |
1557 | 1607 | "plt.show()" |
1558 | 1608 | ] |
1559 | 1609 | }, |
|
0 commit comments