Пример #1
0
				continue   

			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
			
Пример #2
0
			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