Inheritance diagram for lsst.pipe.tasks.selectImages.PsfWcsSelectImagesTask:

Public Member Functions
def	runDataRef (self, dataRef, coordList, makeDataRefList=True, selectDataList=[])

def	run (self, coordList)

Static Public Attributes
	ConfigClass = PsfWcsSelectImagesConfig

Detailed Description

Select images using their Wcs and cuts on the PSF properties

Definition at line 262 of file selectImages.py.

Member Function Documentation

◆ run()

def lsst.pipe.tasks.selectImages.BaseSelectImagesTask.run	(	self,
		coordList
	)

inherited

Select images suitable for coaddition in a particular region

@param[in] coordList: list of coordinates defining region of interest; if None then select all images
subclasses may add additional keyword arguments, as required

@return a pipeBase Struct containing:
- exposureInfoList: a list of exposure information objects (subclasses of BaseExposureInfo),
    which have at least the following fields:
    - dataId: data ID dictionary
    - coordList: ICRS coordinates of the corners of the exposure (list of lsst.geom.SpherePoint)

Definition at line 76 of file selectImages.py.

     def run(self, coordList):
         """Select images suitable for coaddition in a particular region
  
         @param[in] coordList: list of coordinates defining region of interest; if None then select all images
         subclasses may add additional keyword arguments, as required
  
         @return a pipeBase Struct containing:
         - exposureInfoList: a list of exposure information objects (subclasses of BaseExposureInfo),
             which have at least the following fields:
             - dataId: data ID dictionary
             - coordList: ICRS coordinates of the corners of the exposure (list of lsst.geom.SpherePoint)
         """
         raise NotImplementedError()
  

◆ runDataRef()

def lsst.pipe.tasks.selectImages.PsfWcsSelectImagesTask.runDataRef	(	self,
		dataRef,
		coordList,
		makeDataRefList = `True`,
		selectDataList = `[]`
	)

Select images in the selectDataList that overlap the patch and satisfy PSF quality critera.

The PSF quality criteria are based on the size and ellipticity residuals from the
adaptive second moments of the star and the PSF.

The criteria are:
  - the median of the ellipticty residuals
  - the robust scatter of the size residuals (using the median absolute deviation)
  - the robust scatter of the size residuals scaled by the square of
    the median size

@param dataRef: Data reference for coadd/tempExp (with tract, patch)
@param coordList: List of ICRS coordinates (lsst.geom.SpherePoint) specifying boundary of patch
@param makeDataRefList: Construct a list of data references?
@param selectDataList: List of SelectStruct, to consider for selection

Reimplemented from lsst.pipe.tasks.selectImages.WcsSelectImagesTask.

Definition at line 268 of file selectImages.py.

     def runDataRef(self, dataRef, coordList, makeDataRefList=True, selectDataList=[]):
         """Select images in the selectDataList that overlap the patch and satisfy PSF quality critera.
  
         The PSF quality criteria are based on the size and ellipticity residuals from the
         adaptive second moments of the star and the PSF.
  
         The criteria are:
           - the median of the ellipticty residuals
           - the robust scatter of the size residuals (using the median absolute deviation)
           - the robust scatter of the size residuals scaled by the square of
             the median size
  
         @param dataRef: Data reference for coadd/tempExp (with tract, patch)
         @param coordList: List of ICRS coordinates (lsst.geom.SpherePoint) specifying boundary of patch
         @param makeDataRefList: Construct a list of data references?
         @param selectDataList: List of SelectStruct, to consider for selection
         """
         result = super(PsfWcsSelectImagesTask, self).runDataRef(dataRef, coordList, makeDataRefList,
                                                                 selectDataList)
  
         dataRefList = []
         exposureInfoList = []
         for dataRef, exposureInfo in zip(result.dataRefList, result.exposureInfoList):
             butler = dataRef.butlerSubset.butler
             srcCatalog = butler.get('src', dataRef.dataId)
             mask = srcCatalog[self.config.starSelection]
  
             starXX = srcCatalog[self.config.starShape+'_xx'][mask]
             starYY = srcCatalog[self.config.starShape+'_yy'][mask]
             starXY = srcCatalog[self.config.starShape+'_xy'][mask]
             psfXX = srcCatalog[self.config.psfShape+'_xx'][mask]
             psfYY = srcCatalog[self.config.psfShape+'_yy'][mask]
             psfXY = srcCatalog[self.config.psfShape+'_xy'][mask]
  
             starSize = np.power(starXX*starYY - starXY**2, 0.25)
             starE1 = (starXX - starYY)/(starXX + starYY)
             starE2 = 2*starXY/(starXX + starYY)
             medianSize = np.median(starSize)
  
             psfSize = np.power(psfXX*psfYY - psfXY**2, 0.25)
             psfE1 = (psfXX - psfYY)/(psfXX + psfYY)
             psfE2 = 2*psfXY/(psfXX + psfYY)
  
             medianE1 = np.abs(np.median(starE1 - psfE1))
             medianE2 = np.abs(np.median(starE2 - psfE2))
             medianE = np.sqrt(medianE1**2 + medianE2**2)
  
             scatterSize = sigmaMad(starSize - psfSize)
             scaledScatterSize = scatterSize/medianSize**2
  
             valid = True
             if self.config.maxEllipResidual and medianE > self.config.maxEllipResidual:
                 self.log.info("Removing visit %s because median e residual too large: %f vs %f" %
                               (dataRef.dataId, medianE, self.config.maxEllipResidual))
                 valid = False
             elif self.config.maxSizeScatter and scatterSize > self.config.maxSizeScatter:
                 self.log.info("Removing visit %s because size scatter is too large: %f vs %f" %
                               (dataRef.dataId, scatterSize, self.config.maxSizeScatter))
                 valid = False
             elif self.config.maxScaledSizeScatter and scaledScatterSize > self.config.maxScaledSizeScatter:
                 self.log.info("Removing visit %s because scaled size scatter is too large: %f vs %f" %
                               (dataRef.dataId, scaledScatterSize, self.config.maxScaledSizeScatter))
                 valid = False
  
             if valid is False:
                 continue
  
             dataRefList.append(dataRef)
             exposureInfoList.append(exposureInfo)
  
         return pipeBase.Struct(
             dataRefList=dataRefList,
             exposureInfoList=exposureInfoList,
         )
  
  

Member Data Documentation

◆ ConfigClass

lsst.pipe.tasks.selectImages.PsfWcsSelectImagesTask.ConfigClass = PsfWcsSelectImagesConfig

static

Definition at line 265 of file selectImages.py.

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

/j/snowflake/release/lsstsw/stack/0.4.1/Linux64/pipe_tasks/21.0.0-43-g9a3faf29+c9d24e8573/python/lsst/pipe/tasks/selectImages.py