In [ ]:

    




    

In [ ]:

    

from lsst.daf.persistence import Butler
import lsst.afw.image as afwImage
import lsst.afw.display.ds9 as ds9
import pickle

DISPLAY=False


    

In [ ]:

    

butler = Butler('/home/afausti/bulge/output/')

f=0

# You may fix the dataId or iterate over all dataIds in the data repository

for dataref in butler.subset(datasetType='calexp', dataId={'ccdnum': 35, 'visit': 205518}):
# for dataref in butler.subset(datasetType='calexp'):

        if dataref.datasetExists(): # processCcd did not fail

            exp=dataref.get()
            src=dataref.get(datasetType='src') # get the src catalog

            if DISPLAY:
                f += 1
  
                ds9.mtv(exp, frame=f, title="CCD")
                   
                for s in src:
                    ds9.dot("+", *s.getCentroid(), size=10, ctype=ds9.RED, frame=f)

                # Remove footprint mask, and plot image

                mask = exp.getMaskedImage().getMask()
                mask &= ~(mask.getPlaneBitMask("DETECTED"))

                f += 1
                ds9.mtv(exp, frame=f, title="Footprint Mask Removed")