Inheritance diagram for lsst.ip.diffim.dipoleMeasurement.DipoleAnalysis:

Public Member Functions
def	__init__ (self)

def	__call__ (self, source)

def	getSn (self, source)

def	getCentroid (self, source)

def	getOrientation (self, source)

def	displayDipoles (self, exposure, sources)

Detailed Description

Functor class that provides (S/N, position, orientation) of measured dipoles

Definition at line 377 of file dipoleMeasurement.py.

Constructor & Destructor Documentation

◆ init()

def lsst.ip.diffim.dipoleMeasurement.DipoleAnalysis.__init__ ( self )

Definition at line 380 of file dipoleMeasurement.py.

     def __init__(self):
         pass
  

Member Function Documentation

◆ call()

def lsst.ip.diffim.dipoleMeasurement.DipoleAnalysis.__call__	(	self,
		source
	)

Parse information returned from dipole measurement

Parameters
----------
source : `lsst.afw.table.SourceRecord`
    The source that will be examined

Definition at line 383 of file dipoleMeasurement.py.

     def __call__(self, source):
         """Parse information returned from dipole measurement
  
         Parameters
         ----------
         source : `lsst.afw.table.SourceRecord`
             The source that will be examined"""
         return self.getSn(source), self.getCentroid(source), self.getOrientation(source)
  

◆ displayDipoles()

def lsst.ip.diffim.dipoleMeasurement.DipoleAnalysis.displayDipoles	(	self,
		exposure,
		sources
	)

Display debugging information on the detected dipoles

Parameters
----------
exposure : `lsst.afw.image.Exposure`
    Image the dipoles were measured on
sources : `lsst.afw.table.SourceCatalog`
    The set of diaSources that were measured

Definition at line 449 of file dipoleMeasurement.py.

     def displayDipoles(self, exposure, sources):
         """Display debugging information on the detected dipoles
  
         Parameters
         ----------
         exposure : `lsst.afw.image.Exposure`
             Image the dipoles were measured on
         sources : `lsst.afw.table.SourceCatalog`
             The set of diaSources that were measured"""
  
         import lsstDebug
         display = lsstDebug.Info(__name__).display
         displayDiaSources = lsstDebug.Info(__name__).displayDiaSources
         maskTransparency = lsstDebug.Info(__name__).maskTransparency
         if not maskTransparency:
             maskTransparency = 90
         disp = afwDisplay.Display(frame=lsstDebug.frame)
         disp.setMaskTransparency(maskTransparency)
         disp.mtv(exposure)
  
         if display and displayDiaSources:
             with disp.Buffering():
                 for source in sources:
                     cenX, cenY = source.get("ipdiffim_DipolePsfFlux_centroid")
                     if np.isinf(cenX) or np.isinf(cenY):
                         cenX, cenY = source.getCentroid()
  
                     isdipole = source.get("ip_diffim_ClassificationDipole_value")
                     if isdipole and np.isfinite(isdipole):
                         # Dipole
                         ctype = afwDisplay.GREEN
                     else:
                         # Not dipole
                         ctype = afwDisplay.RED
  
                     disp.dot("o", cenX, cenY, size=2, ctype=ctype)
  
                     negCenX = source.get("ip_diffim_PsfDipoleFlux_neg_centroid_x")
                     negCenY = source.get("ip_diffim_PsfDipoleFlux_neg_centroid_y")
                     posCenX = source.get("ip_diffim_PsfDipoleFlux_pos_centroid_x")
                     posCenY = source.get("ip_diffim_PsfDipoleFlux_pos_centroid_y")
                     if (np.isinf(negCenX) or np.isinf(negCenY) or np.isinf(posCenX) or np.isinf(posCenY)):
                         continue
  
                     disp.line([(negCenX, negCenY), (posCenX, posCenY)], ctype=afwDisplay.YELLOW)
  
             lsstDebug.frame += 1
  
  

◆ getCentroid()

def lsst.ip.diffim.dipoleMeasurement.DipoleAnalysis.getCentroid	(	self,
		source
	)

Get the centroid of the dipole; average of positive and negative lobe

Parameters
----------
source : `lsst.afw.table.SourceRecord`
    The source that will be examined

Definition at line 411 of file dipoleMeasurement.py.

     def getCentroid(self, source):
         """Get the centroid of the dipole; average of positive and negative lobe
  
         Parameters
         ----------
         source : `lsst.afw.table.SourceRecord`
             The source that will be examined"""
  
         negCenX = source.get("ip_diffim_PsfDipoleFlux_neg_centroid_x")
         negCenY = source.get("ip_diffim_PsfDipoleFlux_neg_centroid_y")
         posCenX = source.get("ip_diffim_PsfDipoleFlux_pos_centroid_x")
         posCenY = source.get("ip_diffim_PsfDipoleFlux_pos_centroid_y")
         if (np.isinf(negCenX) or np.isinf(negCenY) or np.isinf(posCenX) or np.isinf(posCenY)):
             return None
  
         center = geom.Point2D(0.5*(negCenX+posCenX),
                               0.5*(negCenY+posCenY))
         return center
  

◆ getOrientation()

def lsst.ip.diffim.dipoleMeasurement.DipoleAnalysis.getOrientation	(	self,
		source
	)

Calculate the orientation of dipole; vector from negative to positive lobe

Parameters
----------
source : `lsst.afw.table.SourceRecord`
    The source that will be examined

Definition at line 430 of file dipoleMeasurement.py.

     def getOrientation(self, source):
         """Calculate the orientation of dipole; vector from negative to positive lobe
  
         Parameters
         ----------
         source : `lsst.afw.table.SourceRecord`
             The source that will be examined"""
  
         negCenX = source.get("ip_diffim_PsfDipoleFlux_neg_centroid_x")
         negCenY = source.get("ip_diffim_PsfDipoleFlux_neg_centroid_y")
         posCenX = source.get("ip_diffim_PsfDipoleFlux_pos_centroid_x")
         posCenY = source.get("ip_diffim_PsfDipoleFlux_pos_centroid_y")
         if (np.isinf(negCenX) or np.isinf(negCenY) or np.isinf(posCenX) or np.isinf(posCenY)):
             return None
  
         dx, dy = posCenX-negCenX, posCenY-negCenY
         angle = geom.Angle(np.arctan2(dx, dy), geom.radians)
         return angle
  

◆ getSn()

def lsst.ip.diffim.dipoleMeasurement.DipoleAnalysis.getSn	(	self,
		source
	)

Get the total signal-to-noise of the dipole; total S/N is from positive and negative lobe

Parameters
----------
source : `lsst.afw.table.SourceRecord`
    The source that will be examined

Definition at line 392 of file dipoleMeasurement.py.

     def getSn(self, source):
         """Get the total signal-to-noise of the dipole; total S/N is from positive and negative lobe
  
         Parameters
         ----------
         source : `lsst.afw.table.SourceRecord`
             The source that will be examined"""
  
         posflux = source.get("ip_diffim_PsfDipoleFlux_pos_instFlux")
         posfluxErr = source.get("ip_diffim_PsfDipoleFlux_pos_instFluxErr")
         negflux = source.get("ip_diffim_PsfDipoleFlux_neg_instFlux")
         negfluxErr = source.get("ip_diffim_PsfDipoleFlux_neg_instFluxErr")
  
         # Not a dipole!
         if (posflux < 0) is (negflux < 0):
             return 0
  
         return np.sqrt((posflux/posfluxErr)**2 + (negflux/negfluxErr)**2)
  

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/dipoleMeasurement.py

Public Member Functions