|
102 | 102 | "from lsst.afw.image import ExposureF, LOCAL, PARENT\n", |
103 | 103 | "import lsst.afw.display as afw_display\n", |
104 | 104 | "\n", |
105 | | - "from astropy.visualization import ZScaleInterval\n" |
| 105 | + "from astropy.visualization import ZScaleInterval" |
106 | 106 | ] |
107 | 107 | }, |
108 | 108 | { |
|
329 | 329 | "metadata": {}, |
330 | 330 | "outputs": [], |
331 | 331 | "source": [ |
332 | | - "subset_box = Box2I(minimum=Point2I(x=subset_xmin, y=subset_ymin), maximum=Point2I(x=subset_xmax, y=subset_ymax))\n" |
| 332 | + "subset_box = Box2I(minimum=Point2I(x=subset_xmin, y=subset_ymin),\n", |
| 333 | + " maximum=Point2I(x=subset_xmax, y=subset_ymax))\n" |
333 | 334 | ] |
334 | 335 | }, |
335 | 336 | { |
|
493 | 494 | "metadata": {}, |
494 | 495 | "outputs": [], |
495 | 496 | "source": [ |
496 | | - "full_exposure = butler.get('deep_coadd', dataId=dataId)\n" |
| 497 | + "full_exposure = butler.get('deep_coadd', dataId=dataId)" |
497 | 498 | ] |
498 | 499 | }, |
499 | 500 | { |
|
680 | 681 | "metadata": {}, |
681 | 682 | "outputs": [], |
682 | 683 | "source": [ |
683 | | - "fig, axes = plt.subplots(nrows=nrows, ncols=ncols, figsize=(16, 4 * nrows))\n", |
| 684 | + "fig, axes = plt.subplots(nrows=nrows, ncols=ncols,\n", |
| 685 | + " figsize=(16, 4 * nrows))\n", |
684 | 686 | "axes = axes.flatten()\n", |
685 | 687 | "\n", |
686 | 688 | "deferred_handles = []\n", |
687 | 689 | "for target in targets:\n", |
688 | 690 | "\n", |
689 | | - " point = geom.SpherePoint(target['ra'], target['dec'], geom.degrees)\n", |
| 691 | + " point = geom.SpherePoint(target['ra'], target['dec'],\n", |
| 692 | + " geom.degrees)\n", |
690 | 693 | "\n", |
691 | 694 | " tract_info = skymap.findTract(point)\n", |
692 | 695 | " patch_info = tract_info.findPatch(point)\n", |
|
702 | 705 | "\n", |
703 | 706 | "for idx, (target, handle) in enumerate(deferred_handles):\n", |
704 | 707 | "\n", |
705 | | - " radec = geom.SpherePoint(target['ra'], target['dec'], geom.degrees)\n", |
| 708 | + " radec = geom.SpherePoint(target['ra'], target['dec'],\n", |
| 709 | + " geom.degrees)\n", |
706 | 710 | " xy = Point2D(wcs.skyToPixel(radec))\n", |
707 | 711 | "\n", |
708 | | - " selection_bbox = Box2I.makeCenteredBox(xy, Extent2I(subset_size, subset_size))\n", |
| 712 | + " selection_bbox = Box2I.makeCenteredBox(xy, Extent2I(subset_size,\n", |
| 713 | + " subset_size))\n", |
709 | 714 | " params = {'bbox': selection_bbox}\n", |
710 | | - " \n", |
| 715 | + "\n", |
711 | 716 | " subset = handle.get(parameters=params)\n", |
712 | | - " \n", |
| 717 | + "\n", |
713 | 718 | " print(f\"Processed ID {target['id']}\")\n", |
714 | | - " \n", |
| 719 | + "\n", |
715 | 720 | " img_array = subset.image.array\n", |
716 | 721 | " vmin, vmax = zscale.get_limits(img_array)\n", |
717 | 722 | " ax = axes[idx]\n", |
|
0 commit comments