lsst.ip.diffim.utils.DipoleTestImage Class Reference

Inheritance diagram for lsst.ip.diffim.utils.DipoleTestImage:

Public Member Functions
def	__init__ (self, w=101, h=101, xcenPos=[27.], ycenPos=[25.], xcenNeg=[23.], ycenNeg=[25.], psfSigma=2., flux=[30000.], fluxNeg=None, noise=10., gradientParams=None)
def	fitDipoleSource (self, source, **kwds)
def	detectDipoleSources (self, doMerge=True, diffim=None, detectSigma=5.5, grow=3, minBinSize=32)

Public Attributes
	w
	h
	xcenPos
	ycenPos
	xcenNeg
	ycenNeg
	psfSigma
	flux
	fluxNeg
	noise
	gradientParams

Detailed Description

Utility class for dipole measurement testing.

Generate an image with simulated dipoles and noise; store the original
"pre-subtraction" images and catalogs as well.
Used to generate test data for DMTN-007 (http://dmtn-007.lsst.io).

Definition at line 924 of file utils.py.

Constructor & Destructor Documentation

__init__()

def lsst.ip.diffim.utils.DipoleTestImage.__init__	(		self,
			w = 101,
			h = 101,
			xcenPos = [27.],
			ycenPos = [25.],
			xcenNeg = [23.],
			ycenNeg = [25.],
			psfSigma = 2.,
			flux = [30000.],
			fluxNeg = None,
			noise = 10.,
			gradientParams = None
	)

Definition at line 932 of file utils.py.

                 psfSigma=2., flux=[30000.], fluxNeg=None, noise=10., gradientParams=None):
        self.w = w
        self.h = h
        self.xcenPos = xcenPos
        self.ycenPos = ycenPos
        self.xcenNeg = xcenNeg
        self.ycenNeg = ycenNeg
        self.psfSigma = psfSigma
        self.flux = flux
        self.fluxNeg = fluxNeg
        if fluxNeg is None:
            self.fluxNeg = self.flux
        self.noise = noise
        self.gradientParams = gradientParams
        self._makeDipoleImage()
 

Member Function Documentation

detectDipoleSources()

def lsst.ip.diffim.utils.DipoleTestImage.detectDipoleSources	(		self,
			doMerge = True,
			diffim = None,
			detectSigma = 5.5,
			grow = 3,
			minBinSize = 32
	)

Utility function for detecting dipoles.

Detect pos/neg sources in the diffim, then merge them. A
bigger "grow" parameter leads to a larger footprint which
helps with dipole measurement for faint dipoles.

Parameters
----------
doMerge : `bool`
   Whether to merge the positive and negagive detections into a single
   source table.
diffim : `lsst.afw.image.exposure.exposure.ExposureF`
   Difference image on which to perform detection.
detectSigma : `float`
   Threshold for object detection.
grow : `int`
   Number of pixels to grow the footprints before merging.
minBinSize : `int`
   Minimum bin size for the background (re)estimation (only applies if
   the default leads to min(nBinX, nBinY) < fit order so the default
   config parameter needs to be decreased, but not to a value smaller
   than ``minBinSize``, in which case the fitting algorithm will take
   over and decrease the fit order appropriately.)

Returns
-------
sources : `lsst.afw.table.SourceCatalog`
   If doMerge=True, the merged source catalog is returned OR
detectTask : `lsst.meas.algorithms.SourceDetectionTask`
schema : `lsst.afw.table.Schema`
   If doMerge=False, the source detection task and its schema are
   returned.

Definition at line 1005 of file utils.py.

    def detectDipoleSources(self, doMerge=True, diffim=None, detectSigma=5.5, grow=3, minBinSize=32):
        """Utility function for detecting dipoles.
 
        Detect pos/neg sources in the diffim, then merge them. A
        bigger "grow" parameter leads to a larger footprint which
        helps with dipole measurement for faint dipoles.
 
        Parameters
        ----------
        doMerge : `bool`
           Whether to merge the positive and negagive detections into a single
           source table.
        diffim : `lsst.afw.image.exposure.exposure.ExposureF`
           Difference image on which to perform detection.
        detectSigma : `float`
           Threshold for object detection.
        grow : `int`
           Number of pixels to grow the footprints before merging.
        minBinSize : `int`
           Minimum bin size for the background (re)estimation (only applies if
           the default leads to min(nBinX, nBinY) < fit order so the default
           config parameter needs to be decreased, but not to a value smaller
           than ``minBinSize``, in which case the fitting algorithm will take
           over and decrease the fit order appropriately.)
 
        Returns
        -------
        sources : `lsst.afw.table.SourceCatalog`
           If doMerge=True, the merged source catalog is returned OR
        detectTask : `lsst.meas.algorithms.SourceDetectionTask`
        schema : `lsst.afw.table.Schema`
           If doMerge=False, the source detection task and its schema are
           returned.
        """
        if diffim is None:
            diffim = self.diffim
 
        # Start with a minimal schema - only the fields all SourceCatalogs need
        schema = afwTable.SourceTable.makeMinimalSchema()
 
        # Customize the detection task a bit (optional)
        detectConfig = measAlg.SourceDetectionConfig()
        detectConfig.returnOriginalFootprints = False  # should be the default
 
        psfSigma = diffim.getPsf().computeShape().getDeterminantRadius()
 
        # code from imageDifference.py:
        detectConfig.thresholdPolarity = "both"
        detectConfig.thresholdValue = detectSigma
        # detectConfig.nSigmaToGrow = psfSigma
        detectConfig.reEstimateBackground = True  # if False, will fail often for faint sources on gradients?
        detectConfig.thresholdType = "pixel_stdev"
        # Test images are often quite small, so may need to adjust background binSize
        while ((min(diffim.getWidth(), diffim.getHeight()))//detectConfig.background.binSize
               < detectConfig.background.approxOrderX and detectConfig.background.binSize > minBinSize):
            detectConfig.background.binSize = max(minBinSize, detectConfig.background.binSize//2)
 
        # Create the detection task. We pass the schema so the task can declare a few flag fields
        detectTask = measAlg.SourceDetectionTask(schema, config=detectConfig)
 
        table = afwTable.SourceTable.make(schema)
        catalog = detectTask.run(table, diffim, sigma=psfSigma)
 
        # Now do the merge.
        if doMerge:
            fpSet = catalog.fpSets.positive
            fpSet.merge(catalog.fpSets.negative, grow, grow, False)
            sources = afwTable.SourceCatalog(table)
            fpSet.makeSources(sources)
 
            return sources
 
        else:
            return detectTask, schema

fitDipoleSource()

def lsst.ip.diffim.utils.DipoleTestImage.fitDipoleSource	(		self,
			source,
		**	kwds
	)

Definition at line 1000 of file utils.py.

    def fitDipoleSource(self, source, **kwds):
        alg = DipoleFitAlgorithm(self.diffim, self.posImage, self.negImage)
        fitResult = alg.fitDipole(source, **kwds)
        return fitResult
 

Member Data Documentation

flux

lsst.ip.diffim.utils.DipoleTestImage.flux

Definition at line 941 of file utils.py.

fluxNeg

lsst.ip.diffim.utils.DipoleTestImage.fluxNeg

Definition at line 942 of file utils.py.

gradientParams

lsst.ip.diffim.utils.DipoleTestImage.gradientParams

Definition at line 946 of file utils.py.

h

lsst.ip.diffim.utils.DipoleTestImage.h

Definition at line 935 of file utils.py.

noise

lsst.ip.diffim.utils.DipoleTestImage.noise

Definition at line 945 of file utils.py.

psfSigma

lsst.ip.diffim.utils.DipoleTestImage.psfSigma

Definition at line 940 of file utils.py.

w

lsst.ip.diffim.utils.DipoleTestImage.w

Definition at line 934 of file utils.py.

xcenNeg

lsst.ip.diffim.utils.DipoleTestImage.xcenNeg

Definition at line 938 of file utils.py.

xcenPos

lsst.ip.diffim.utils.DipoleTestImage.xcenPos

Definition at line 936 of file utils.py.

ycenNeg

lsst.ip.diffim.utils.DipoleTestImage.ycenNeg

Definition at line 939 of file utils.py.

ycenPos

lsst.ip.diffim.utils.DipoleTestImage.ycenPos

Definition at line 937 of file utils.py.

---

The documentation for this class was generated from the following file:

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-0.7.0/Linux64/ip_diffim/22.0.1-24-g1ad7a390+a9625a72a8/python/lsst/ip/diffim/utils.py