|
| 1 | +.. _dia: |
| 2 | + |
| 3 | +############################### |
| 4 | +Difference image analysis (DIA) |
| 5 | +############################### |
| 6 | + |
| 7 | +.. important:: |
| 8 | + |
| 9 | + **Placeholder.** Prompt data processing is under construction. |
| 10 | + |
| 11 | + |
| 12 | +DIA uses the decorrelated `Alard and Lupton (1998) <https://ui.adsabs.harvard.edu/abs/1998ApJ...503..325A/abstract>`_ |
| 13 | +image differencing algorithm, |
| 14 | +as described in "Implementation of Image Difference Decorrelation" (`dmtn-021.lsst.io <https://dmtn-021.lsst.io/>`_) |
| 15 | +and in the LSST Science Pipelines documentation for |
| 16 | +the `Alard-Lupton implementation <https://pipelines.lsst.io/modules/lsst.ip.diffim/AL_implementation.html>`_. |
| 17 | + |
| 18 | +Key terminology: |
| 19 | + |
| 20 | +* ``DiaSource``: a detection, positive or negative, in a single difference image |
| 21 | +* ``DiaObject``: the astrophysical transient or variable object at a static sky coordinate |
| 22 | + |
| 23 | + |
| 24 | +**Light curves:** the forced photometry fluxes should be used for light curves. |
| 25 | + |
| 26 | +All measurements are stored in the catalog :doc:`/products/index`. |
| 27 | + |
| 28 | + |
| 29 | +Difference images |
| 30 | +================= |
| 31 | + |
| 32 | +Template coadd images are warped, scaled, and PSF-matched to a visit image, |
| 33 | +and then subtracted from the visit image to create a difference image. |
| 34 | + |
| 35 | +In the context of DIA, the visit image is sometimes referred to as the |
| 36 | +"direct" or "science" image. |
| 37 | + |
| 38 | + |
| 39 | +DIA source detection |
| 40 | +==================== |
| 41 | + |
| 42 | +DIA source detection refers to the process of finding regions with above-threshold |
| 43 | +absolute flux in difference images, where the threshold is a signal-to-noise ratio of 5. |
| 44 | + |
| 45 | +Difference images contain sources of residual flux from the subtraction, both |
| 46 | +positive and negative. |
| 47 | +Real astrophysical transients, variables, and moving objects appear as sources with |
| 48 | +a flux that is equal to the difference between their flux in the template and visit images. |
| 49 | + |
| 50 | +Spurious sources are also present due to unavoidably imperfect template matching. |
| 51 | + |
| 52 | +Deblending is not run on difference image sources. |
| 53 | + |
| 54 | +The results of DIA source detection are stored in the ``DiaSource`` catalog. |
| 55 | + |
| 56 | + |
| 57 | +Reliability score |
| 58 | +----------------- |
| 59 | + |
| 60 | +Both real and spurious DIA sources can be detected with a signal-to-noise ratio of 5. |
| 61 | +An estimate of how likely a source is to be real is provided to enable filtering of DIA sources by users. |
| 62 | + |
| 63 | +Dipoles |
| 64 | +------- |
| 65 | + |
| 66 | +Sources with footprints containing both positive and negative peaks are fit with a |
| 67 | +dipole centroid code. |
| 68 | +Dipoles represent sources that have moved only slightly (e.g., less than the PSF) |
| 69 | +between template and visit image. |
| 70 | + |
| 71 | +Pixel flags |
| 72 | +----------- |
| 73 | + |
| 74 | +If the footprint contains one or more flagged pixels, |
| 75 | +e.g., for cosmic rays, detector edge, bad pixels, |
| 76 | +known artifacts, non-astrophysical trails, |
| 77 | +the source is also flagged. |
| 78 | + |
| 79 | + |
| 80 | + |
| 81 | +DIA object association |
| 82 | +====================== |
| 83 | + |
| 84 | +DIA sources are associated into DIA objects using a radius of 1 arcsecond. |
| 85 | + |
| 86 | +The DIA objects are stored in the ``DiaObject`` catalog. |
| 87 | + |
| 88 | + |
| 89 | +Variability characterization |
| 90 | +---------------------------- |
| 91 | + |
| 92 | +Light curve characterization parameters are calculated for all DIA objects, |
| 93 | +such as the minimum, mean, maximum, standard deviation, and skew in the |
| 94 | +difference image fluxes, the light curve's slope, percentiles, and ``StetsonJ`` |
| 95 | +parameter. |
| 96 | + |
| 97 | + |
| 98 | +DIA forced photometry |
| 99 | +===================== |
| 100 | + |
| 101 | +In general, "forced" photometry means a measurement made at a fixed coordinate in an image, |
| 102 | +regardless of whether an above-threshold region was detected in that particular image. |
| 103 | + |
| 104 | +Forced PSF photometry measurements are made on all visit and difference images |
| 105 | +at the locations of all DIA objects. |
| 106 | +The results are stored in the ``ForcedSourceOnDiaObject`` catalog. |
| 107 | + |
0 commit comments