def getImage(self, ccd, amp=None, imageFactory=None):
"""Return an image of dlnI/dr"""
ccdNum = ccd.getId().getSerial()
try:
if self.kwargs.get("ccd") is not None and not ccdNum in self.kwargs.get("ccd"):
raise RuntimeError
dataId = self.kwargs
if dataId.has_key("ccd"):
dataId = self.kwargs.copy()
del dataId["ccd"]
raw = cgUtils.trimExposure(self.butler.get("raw", ccd=ccdNum, immediate=True, **dataId).convertF(),
subtractBias=True, rotate=True)
if False:
msk = raw.getMaskedImage().getMask()
BAD=msk.getPlaneBitMask("BAD")
msk |= BAD
msk[15:-15, 0:-20] &= ~BAD
except Exception, e:
if self.verbose and str(e):
print e
ccdImage = afwImage.ImageF(*ccd.getAllPixels(self.isTrimmed).getDimensions())
if self.bin:
ccdImage = afwMath.binImage(ccdImage, self.bin)
return ccdImage
def getImage(self, ccd, amp=None, imageFactory=None):
"""Return an image of dlnI/dr"""
ccdNum = ccd.getId().getSerial()
try:
if self.kwargs.get("ccd") is not None and ccdNum not in self.kwargs.get("ccd"):
raise RuntimeError
dataId = self.kwargs
if "ccd" in dataId:
dataId = self.kwargs.copy()
del dataId["ccd"]
raw = cgUtils.trimExposure(self.butler.get("raw", ccd=ccdNum, immediate=True,
**dataId).convertF(),
subtractBias=True, rotate=True)
if False:
msk = raw.getMaskedImage().getMask()
BAD = msk.getPlaneBitMask("BAD")
msk |= BAD
msk[15:-15, 0:-20] &= ~BAD
except Exception as e:
if self.verbose and str(e):
print(e)
ccdImage = afwImage.ImageF(*ccd.getAllPixels(self.isTrimmed).getDimensions())
if self.bin:
ccdImage = afwMath.binImage(ccdImage, self.bin)
return ccdImage
return fitPatches(raw, bin=self.bin, returnResidualImage=False)[2]
def __call__(self, dataId):
if self.verbose:
print "Visit %(visit)d CCD%(ccd)03d\r" % dataId
raw = self.butler.get("raw", immediate=True, **dataId)
try:
raw = cgUtils.trimExposure(raw, subtractBias=True, rotate=True)
except:
return None
if False:
msk=raw.getMaskedImage().getMask()
BAD=msk.getPlaneBitMask("BAD"); msk |= BAD; msk[15:-15, 20:-20] &= ~BAD
r, lnGrad, theta = [], [], []
fitPatches(raw, bin=self.bin, r=r, lnGrad=lnGrad, theta=theta)
return r, theta, lnGrad
def __call__(self, dataId):
if self.verbose:
print "Visit %(visit)d CCD%(ccd)03d\r" % dataId
raw = self.butler.get("raw", immediate=True, **dataId)
try:
raw = cgUtils.trimExposure(raw, subtractBias=True, rotate=True)
except:
return None
if False:
msk = raw.getMaskedImage().getMask()
BAD = msk.getPlaneBitMask("BAD")
msk |= BAD
msk[15:-15, 20:-20] &= ~BAD
r, lnGrad, theta = [], [], []
fitPatches(raw, bin=self.bin, r=r, lnGrad=lnGrad, theta=theta)
return r, theta, lnGrad
