prevCatalog = channelTempCatalogs[channel] masterObjectList = allObjects[channel] singleChannelFrame = wholeFrame[channel] assembledChannelFrame = numpy.zeros((frameInfo.nxmax, frameInfo.nymax)) #stackedFrame = numpy.zeros((frameInfo.nxmax, frameInfo.nymax)) stackedImageObject = stackedImages[channel] newObjects = [] for j in range(frameInfo.numWindows): windowImage = singleChannelFrame[j] tmpFilename = ultracamutils.createFITS(trueFrameNumber, j, channel, windowImage, arg.runname) catFilename = ultracamutils.runSex(tmpFilename, arg.apertures) newObjectsinWindow = ultracamutils.readSexObjects(catFilename, config.SEX_MAGNITUDE) newObjectsinWindow = ultracamutils.rejectBadObjects(newObjectsinWindow) if config.KEEP_TMP_FILES!="1": ultracamutils.removeTMPFile(tmpFilename) ultracamutils.removeTMPFile(catFilename) xll = frameInfo.getWindow(j).xll yll = frameInfo.getWindow(j).yll xsize = frameInfo.getWindow(j).xsize ysize = frameInfo.getWindow(j).ysize #print "windowImage.shape:", windowImage.shape #print "assembledChannelFrame.shape:", assembledChannelFrame.shape leftPixel = xll-1
prevCatalog = channelTempCatalogs[channel] masterObjectList = allObjects[channel] singleChannelFrame = wholeFrame[channel] assembledChannelFrame = numpy.zeros((frameInfo.nxmax, frameInfo.nymax)) assembledChannelFrameUnequalised = numpy.zeros((frameInfo.nxmax, frameInfo.nymax)) #stackedFrame = numpy.zeros((frameInfo.nxmax, frameInfo.nymax)) stackedImageObject = stackedImages[channel] newObjects = [] for j in range(frameInfo.numWindows): windowImage = singleChannelFrame[j] tmpFilename = ultracamutils.createFITS(trueFrameNumber, j, channel, windowImage) catFilename = ultracamutils.runSex(tmpFilename) newObjectsinWindow = ultracamutils.readSexObjects(catFilename) newObjectsinWindow = ultracamutils.rejectBadObjects(newObjectsinWindow) if config.KEEP_TMP_FILES!="1": ultracamutils.removeTMPFile(tmpFilename) ultracamutils.removeTMPFile(catFilename) xll = frameInfo.getWindow(j).xll yll = frameInfo.getWindow(j).yll xsize = frameInfo.getWindow(j).xsize ysize = frameInfo.getWindow(j).ysize assembledChannelFrame[xll:xll+xsize, yll:yll+ysize] = assembledChannelFrame[xll:xll+xsize, yll:yll+ysize] + ultracamutils.percentiles(windowImage, 20, 98) assembledChannelFrameUnequalised[xll:xll+xsize, yll:yll+ysize] = assembledChannelFrameUnequalised[xll:xll+xsize, yll:yll+ysize] + windowImage #stackedFrame[xll:xll+xsize, yll:yll+ysize] = stackedFrame[xll:xll+xsize, yll:yll+ysize] + windowImage