lsst.pipe.tasks.mocks.visualization Namespace Reference

Functions
def	plotObservations (catalog, wcs)
def	plotPatches (tractInfo)
def	plotTruth (catalog, wcs)
def	displayImages (root)
def	makePlots (root)

Function Documentation

displayImages()

def lsst.pipe.tasks.mocks.visualization.displayImages

(

root

)

Display coadd images in different frames, with the bounding boxes of the
observations that went into them overlayed.

Definition at line 73 of file visualization.py.

def displayImages(root):
    """Display coadd images in different frames, with the bounding boxes of the
    observations that went into them overlayed.
    """
    import lsst.afw.display as afwDisplay
    afwDisplay.setDefaultMaskTransparency(75)
 
    butler = lsst.daf.persistence.Butler(root=root)
    skyMap = butler.get("deepCoadd_skyMap")
    tractInfo = skyMap[0]
    task = lsst.pipe.tasks.mocks.MockCoaddTask()
    coadds = [patchRef.get("deepCoadd", immediate=True)
              for patchRef in task.iterPatchRefs(butler, tractInfo)]
    for n, coadd in enumerate(coadds):
        afwDisplay.Display(frame=n + 1).mtv(coadd, title="displayImages: coadd")
    for n, coadd in enumerate(coadds):
        afwDisplay.utils.drawCoaddInputs(coadd, frame=n + 1)
    return butler
 
 

makePlots()

def lsst.pipe.tasks.mocks.visualization.makePlots

(

root

)

Convenience function to make all matplotlib plots.

Definition at line 93 of file visualization.py.

def makePlots(root):
    """Convenience function to make all matplotlib plots.
    """
    import lsst.pipe.tasks.mocks
    import lsst.daf.persistence
    butler = lsst.daf.persistence.Butler(root=root)
    skyMap = butler.get("deepCoadd_skyMap")
    observations = butler.get("observations", tract=0)
    truth = butler.get("truth", tract=0)
    tractInfo = skyMap[0]
    plotPatches(tractInfo)
    plotObservations(observations, tractInfo.getWcs())
    plotTruth(truth, tractInfo.getWcs())
    pyplot.axis("scaled")
    pyplot.show()
    return butler

plotObservations()

def lsst.pipe.tasks.mocks.visualization.plotObservations	(		catalog,
			wcs
	)

Plot the bounding boxes of an observation catalog (see MockCoaddTask.buildObservationCatalog)
using matplotlib, in the coordinates defined by the given Wcs (usually a skymap Wcs).

Definition at line 27 of file visualization.py.

def plotObservations(catalog, wcs):
    """Plot the bounding boxes of an observation catalog (see MockCoaddTask.buildObservationCatalog)
    using matplotlib, in the coordinates defined by the given Wcs (usually a skymap Wcs).
    """
    for record in catalog:
        box = lsst.geom.Box2D(record.getBBox())
        x = []
        y = []
        iWcs = record.getWcs()
        for xi, yi in box.getCorners():
            try:
                coord = iWcs.pixelToSky(xi, yi)
                xo, yo = wcs.skyToPixel(coord)
                x.append(xo)
                y.append(yo)
            except Exception as e:
                print("WARNING: point %d, %d failed: error=%s" % (xi, yi, e))
        pyplot.fill(x, y, facecolor='r', alpha=0.1, edgecolor=None)
 
 

plotPatches()

def lsst.pipe.tasks.mocks.visualization.plotPatches

(

tractInfo

)

Plot the patches in a skymap tract using matplotlib.

Definition at line 47 of file visualization.py.

def plotPatches(tractInfo):
    """Plot the patches in a skymap tract using matplotlib.
    """
    nPatchX, nPatchY = tractInfo.getNumPatches()
    for iPatchX in range(nPatchX):
        for iPatchY in range(nPatchY):
            patchInfo = tractInfo.getPatchInfo((iPatchX, iPatchY))
            xp1, yp1 = list(zip(*patchInfo.getOuterBBox().getCorners()))
            xp2, yp2 = list(zip(*patchInfo.getInnerBBox().getCorners()))
            pyplot.fill(xp1, yp1, fill=False, edgecolor='g', linestyle='dashed')
            pyplot.fill(xp2, yp2, fill=False, edgecolor='g')
 
 

plotTruth()

def lsst.pipe.tasks.mocks.visualization.plotTruth	(		catalog,
			wcs
	)

Plot the objects in a truth catalog as dots using matplotlib, in the coordinate
system defined by the given Wcs.

Definition at line 60 of file visualization.py.

def plotTruth(catalog, wcs):
    """Plot the objects in a truth catalog as dots using matplotlib, in the coordinate
    system defined by the given Wcs.
    """
    xp = []
    yp = []
    for record in catalog:
        x, y = wcs.skyToPixel(record.getCoord())
        xp.append(x)
        yp.append(y)
    pyplot.plot(xp, yp, 'k+')