if (int(config.KEEP_TMP_FILES)==0):utils.removeCAT(frameCounter) objects = utils.rejectBadObjects(objects) updateMasterObjectList(objects, frameInfo.getWindow(j)) debug.write("Total objects in this window:" + str(len(objects))) averageFrames = [] normalisedImages = [] for j in range(frameInfo.numWindows): averageFrame = numpy.divide(summedFrames[j], frameCounter) averageFrames.append(averageFrame) normalisedImage = utils.percentiles(averageFrame, 25, 99) normalisedImages.append(normalisedImage) # Construct full frame from Windows fullFrame = numpy.zeros((frameInfo.nxmax, frameInfo.nymax)) for j in range(frameInfo.numWindows): xll = frameInfo.getWindow(j).xll yll = frameInfo.getWindow(j).yll xsize = frameInfo.getWindow(j).xsize ysize = frameInfo.getWindow(j).ysize rotatedImage = numpy.flipud(normalisedImages[j]) fullFrame[xll:xll+xsize, yll:yll+ysize] = fullFrame[xll:xll+xsize, yll:yll+ysize] + rotatedImage fullFrame = numpy.fliplr(fullFrame)
FITSHeaders.append(FITSHeader) print "keyword: ", FITSHeader['keyword'], " value: ", FITSHeader['value'], " comment: ", FITSHeader['comment'] headFile.close() FITSInputhdulist = astropy.io.fits.open(fitsFilename) print FITSInputhdulist.info() print FITSInputhdulist[0].data imageData = FITSInputhdulist[0].data if imagePreview: imgplot = plt.imshow(utils.percentiles(imageData, 20, 98), cmap='gray', interpolation='nearest') plt.show() print "Doing some more stuff" FITSInputhdulist.close() debug.write("Writing FITS file: " + outputFilename, level=1) prihdr = astropy.io.fits.Header() prihdr['COMMENT'] = "This file created by makefits.py from the Ultracam pipeline." prihdr['OBJECT'] = "KOI-823" for h in FITSHeaders: print "Adding: ", h["keyword"], ":", h["value"] keyword = str(h["keyword"]) try: prihdr[keyword] = float(h["value"])
# Construct full frame from Windows fullFrame = numpy.zeros((frameInfo.nxmax, frameInfo.nymax)) for j in range(frameInfo.numWindows): frameImage = frameWindows[j + frameInfo.numWindows * i] frameImage = numpy.rot90(frameImage) # All windows need to be rotated by 180 degrees before using the Image library to write to png xll = frameInfo.getWindow(j).xll yll = frameInfo.getWindow(j).yll xsize = frameInfo.getWindow(j).xsize ysize = frameInfo.getWindow(j).ysize rotatedImage = numpy.flipud(frameImage) normalisedImage = utils.percentiles(rotatedImage, 10, 98) fullFrame[xll : xll + xsize, yll : yll + ysize] = ( fullFrame[xll : xll + xsize, yll : yll + ysize] + normalisedImage ) fullFrame = numpy.fliplr(fullFrame) imgData = fullFrame imgSize = numpy.shape(imgData) imgLength = imgSize[0] * imgSize[1] testData = numpy.reshape(imgData, imgLength, order="F") img = Image.new("L", imgSize) img.putdata(testData) if colour == "b": if (frameCounter % rdat.nblue) != 0:
summedFrame = numpy.zeros((width, height)) summedFrames.append(summedFrame) firstFrame = False for j in range(frameInfo.numWindows): frameImage = numpy.copy(frameWindows[j]) utils.createFITS(frameCounter, frameImage) utils.runSex(frameCounter) summedFrame = summedFrames[j] summedFrame = numpy.add(summedFrame, frameImage) summedFrames[j] = summedFrame if mplPreview==True: imgplot = plt.imshow(utils.percentiles(frameImage, 20, 98), cmap='gray', interpolation='nearest') plt.draw() print "Window info:", frameInfo.getWindow(j) if (int(config.KEEP_TMP_FILES)==0):utils.removeFITS(frameCounter); objects = utils.readSexObjects(frameCounter) if (int(config.KEEP_TMP_FILES)==0):utils.removeCAT(frameCounter) objects = utils.rejectBadObjects(objects) updateMasterObjectList(objects, frameInfo.getWindow(j)) debug.write("Total objects in this window:" + str(len(objects)))