lsst.ip.diffim.diaSourceAnalysis Namespace Reference

Classes
class	DiaSourceAnalystConfig
class	DiaSourceAnalyst

Functions
def	parseOptions
def	readSourceSet
def	main

Variables
int	scaling = 5

Function Documentation

def lsst.ip.diffim.diaSourceAnalysis.main

(

)

Main program

Examples:

footprintFunctor.cc, forEachPixel.cc, image1.cc, image2.cc, imageStatistics.cc, maskedImage1.cc, maskedImage2.cc, spatialCellExample.cc, splineInterpolate.cc, and statistics.cc.

Definition at line 193 of file diaSourceAnalysis.py.

def main():
    """Main program"""
    options, args = parseOptions()
    (crDiffSourceFile, crDiffExposureFile) = args

    crDiffSources = readSourceSet(crDiffSourceFile)
    crDiffExposure = afwImage.ExposureF(crDiffExposureFile)
    #import pdb; pdb.set_trace()
    
    analyst = DiaSourceAnalyst()
    
    expX0 = crDiffExposure.getX0()
    expY0 = crDiffExposure.getY0()
    expX1 = expX0 + crDiffExposure.getWidth() - 1
    expY1 = expY0 + crDiffExposure.getHeight() - 1
    
    for i in range(crDiffSources.size()):
        crDiffSource = crDiffSources[i]

        # This segfaults; revisit once the stack stabilizes
        #footprint    = crDiffSource.getFootprint()
        #bbox         = footprint.getBBox()
        
        xAstrom      = crDiffSource.getXAstrom()
        yAstrom      = crDiffSource.getYAstrom()
        Ixx          = max(1.0, crDiffSource.getIxx())
        Iyy          = max(1.0, crDiffSource.getIyy())
        x0           = max(expX0, int(xAstrom - scaling * Ixx))
        x1           = min(expX1, int(xAstrom + scaling * Ixx))
        y0           = max(expY0, int(yAstrom - scaling * Iyy))
        y1           = min(expY1, int(yAstrom + scaling * Iyy))
        bbox         = afwGeom.Box2I(afwGeom.Point2I(x0, y0),
                                     afwGeom.Point2I(x1, y1))
        subExp       = afwImage.ExposureF(crDiffExposure, bbox)
        subMi        = subExp.getMaskedImage()
        imArr, maArr, varArr = subMi.getArrays()

        if analyst.testSource(crDiffSource, subMi):
            ds9.mtv(subExp, frame=1)