|
6 | 6 | "source": [ |
7 | 7 | "<img align=\"left\" src = https://project.lsst.org/sites/default/files/Rubin-O-Logo_0.png width=250 style=\"padding: 10px\"> \n", |
8 | 8 | "<b>Introduction to Source Detection</b> <br>\n", |
9 | | - "Last verified to run on 2021-09-24 with LSST Science Pipelines release w_2021_33 <br>\n", |
| 9 | + "Last verified to run on 2021-10-12 with LSST Science Pipelines release w_2021_40 <br>\n", |
10 | 10 | "Contact author: Alex Drlica-Wagner, Melissa Graham <br>\n", |
11 | 11 | "Target audience: All DP0 delegates. <br>\n", |
12 | 12 | "Container Size: medium <br>\n", |
|
88 | 88 | "from lsst.meas.base import SingleFrameMeasurementTask" |
89 | 89 | ] |
90 | 90 | }, |
| 91 | + { |
| 92 | + "cell_type": "code", |
| 93 | + "execution_count": null, |
| 94 | + "metadata": {}, |
| 95 | + "outputs": [], |
| 96 | + "source": [ |
| 97 | + "# prevent some helpful but ancillary warning messages from printing\n", |
| 98 | + "# during some LSST DM Release calls\n", |
| 99 | + "warnings.simplefilter(\"ignore\", category=UserWarning)\n", |
| 100 | + "warnings.simplefilter(\"ignore\", category=FutureWarning)" |
| 101 | + ] |
| 102 | + }, |
91 | 103 | { |
92 | 104 | "cell_type": "code", |
93 | 105 | "execution_count": null, |
|
302 | 314 | "x_target, y_target = 1700, 2100\n", |
303 | 315 | "width, height = 400, 400\n", |
304 | 316 | "xmin, ymin = x_target-width//2, y_target-height//2\n", |
| 317 | + "point = geom.Point2D(x_target,y_target)\n", |
305 | 318 | "\n", |
306 | 319 | "# Define a small region for a cutout\n", |
307 | 320 | "bbox = geom.Box2I()\n", |
|
493 | 506 | "# Image characterization (this cell may take a few seconds)\n", |
494 | 507 | "result = charImageTask.run(calexp)\n", |
495 | 508 | "\n", |
| 509 | + "# Get the PSF at our point of interest\n", |
496 | 510 | "psf = calexp.getPsf()\n", |
497 | | - "sigma = psf.computeShape().getDeterminantRadius()\n", |
| 511 | + "sigma = psf.computeShape(point).getDeterminantRadius()\n", |
498 | 512 | "pixelScale = calexp.getWcs().getPixelScale().asArcseconds()\n", |
499 | 513 | "\n", |
500 | 514 | "# The factor of 2.355 converts from std to fwhm\n", |
501 | | - "print('psf fwhm = {:.2f} arcsec'.format(sigma*pixelScale*2.355))" |
| 515 | + "#print('psf fwhm = {:.2f} arcsec'.format(sigma*pixelScale*2.355))" |
502 | 516 | ] |
503 | 517 | }, |
504 | 518 | { |
|
0 commit comments