LSSTApplications: lsst.meas.astrom.catalogStarSelector.CatalogStarSelector Class Reference

Inheritance diagram for lsst.meas.astrom.catalogStarSelector.CatalogStarSelector:

Public Member Functions
def	__init__

def	selectStars

Static Public Attributes
	ConfigClass = CatalogStarSelectorConfig

Private Attributes
	_kernelSize

	_borderWidth

	_fluxLim

	_fluxMax

	_badStarPixelFlags

Detailed Description

Definition at line 78 of file catalogStarSelector.py.

Constructor & Destructor Documentation

def lsst.meas.astrom.catalogStarSelector.CatalogStarSelector.__init__	(	self,
		config = `None`
	)

Construct a star selector that uses second moments

This is a naive algorithm and should be used with caution.

@param[in] config: An instance of CatalogStarSelectorConfig

Definition at line 81 of file catalogStarSelector.py.

 
     def __init__(self, config=None):
         """Construct a star selector that uses second moments
         
         This is a naive algorithm and should be used with caution.
         
         @param[in] config: An instance of CatalogStarSelectorConfig
         """
         if not config:
             config = CatalogStarSelector.ConfigClass()
 
         self._kernelSize  = config.kernelSize
         self._borderWidth = config.borderWidth
         self._fluxLim  = config.fluxLim
         self._fluxMax  = config.fluxMax
         self._badStarPixelFlags = config.badStarPixelFlags

lsst::meas.astrom.catalogStarSelector.CatalogStarSelector.__init__

def __init__

Definition: catalogStarSelector.py:81

lsst::meas.astrom.catalogStarSelector.CatalogStarSelector._fluxLim

_fluxLim

Definition: catalogStarSelector.py:93

lsst::meas.astrom.catalogStarSelector.CatalogStarSelector._kernelSize

_kernelSize

Definition: catalogStarSelector.py:91

lsst::meas.astrom.catalogStarSelector.CatalogStarSelector._badStarPixelFlags

_badStarPixelFlags

Definition: catalogStarSelector.py:95

lsst::meas.astrom.catalogStarSelector.CatalogStarSelector._fluxMax

_fluxMax

Definition: catalogStarSelector.py:94

lsst::meas.astrom.catalogStarSelector.CatalogStarSelector._borderWidth

_borderWidth

Definition: catalogStarSelector.py:92

def lsst.meas.astrom.catalogStarSelector.CatalogStarSelector.selectStars	(	self,
		exposure,
		sources,
		matches = `None`
	)

Return a list of PSF candidates that represent likely stars

A list of PSF candidates may be used by a PSF fitter to construct a PSF.

@param[in] exposure: the exposure containing the sources
@param[in] sources: a source list containing sources that may be stars
@param[in] matches: a match vector as produced by meas_astrom; not actually optional
            (passing None just allows us to handle the exception better here
            than in calling code)

@return psfCandidateList: a list of PSF candidates.

Definition at line 97 of file catalogStarSelector.py.

 
     def selectStars(self, exposure, sources, matches=None):
         """Return a list of PSF candidates that represent likely stars
         
         A list of PSF candidates may be used by a PSF fitter to construct a PSF.
         
         @param[in] exposure: the exposure containing the sources
         @param[in] sources: a source list containing sources that may be stars
         @param[in] matches: a match vector as produced by meas_astrom; not actually optional
                             (passing None just allows us to handle the exception better here
                             than in calling code)
         
         @return psfCandidateList: a list of PSF candidates.
         """
         import lsstDebug
         display = lsstDebug.Info(__name__).display
         displayExposure = lsstDebug.Info(__name__).displayExposure     # display the Exposure + spatialCells
         pauseAtEnd = lsstDebug.Info(__name__).pauseAtEnd               # pause when done
 
         if matches is None:
             raise RuntimeError(
                 "Cannot use catalog star selector without running astrometry."
                 )
 
         mi = exposure.getMaskedImage()
         
         if display:
             frames = {}
             if displayExposure:
                 frames["displayExposure"] = 1
                 ds9.mtv(mi, frame=frames["displayExposure"], title="PSF candidates")
         #
         # Read the reference catalogue
         #
         isGoodSource = CheckSource(sources, self._fluxLim, self._fluxMax, self._badStarPixelFlags)
 
         #
         # Go through and find all the PSFs in the catalogue
         #
         # We'll split the image into a number of cells, each of which contributes only
         # one PSF candidate star
         #
         psfCandidateList = []
 
         with ds9.Buffering():
             for ref, source, d in matches:
                 if not ref.get("resolved"):
                     if not isGoodSource(source):
                         symb, ctype = "+", ds9.RED
                     else:
                         try:
                             psfCandidate = measAlg.makePsfCandidate(source, exposure)
 
                             # The setXXX methods are class static, but it's convenient to call them on
                             # an instance as we don't know Exposure's pixel type
                             # (and hence psfCandidate's exact type)
                             if psfCandidate.getWidth() == 0:
                                 psfCandidate.setBorderWidth(self._borderWidth)
                                 psfCandidate.setWidth(self._kernelSize + 2*self._borderWidth)
                                 psfCandidate.setHeight(self._kernelSize + 2*self._borderWidth)
 
                             im = psfCandidate.getMaskedImage().getImage()
                             max = afwMath.makeStatistics(im, afwMath.MAX).getValue()
                             if not numpy.isfinite(max):
                                 continue
                             psfCandidateList.append(psfCandidate)
 
                             symb, ctype = "+", ds9.GREEN
                         except Exception as err:
                             symb, ctype = "o", ds9.RED
                             print "RHL", err
                             pass # FIXME: should log this!
                 else:
                     symb, ctype = "o", ds9.BLUE
 
                 if display and displayExposure:
                     ds9.dot(symb, source.getX() - mi.getX0(), source.getY() - mi.getY0(),
                             size=4, frame=frames["displayExposure"], ctype=ctype)
 
         if display and pauseAtEnd:
             raw_input("Continue? y[es] p[db] ")
 
         return psfCandidateList
 
 measAlg.starSelectorRegistry.register("catalog", CatalogStarSelector)
 

Public Member Functions

Static Public Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation