def setup_module(module): lsst.utils.tests.init() if os.path.exists(DATAREPO_ROOT): print("Deleting existing repo: %r" % (DATAREPO_ROOT,)) shutil.rmtree(DATAREPO_ROOT) if not haveMeasBase: raise unittest.SkipTest("meas_base could not be imported") # Create a task that creates simulated images and builds a coadd from them mocksTask = lsst.pipe.tasks.mocks.MockCoaddTask() # Create an instance of DetectCoaddSourcesTask to measure on the coadd. # There's no noise in these images, so we set a direct-value threshold, # and the background weighting (when using Approximate) to False detectConfig = DetectCoaddSourcesTask.ConfigClass() # Images have no noise, so we can't use the default DynamicDetectionTask detectConfig.detection.retarget(lsst.meas.algorithms.SourceDetectionTask) detectConfig.detection.thresholdType = "value" detectConfig.detection.thresholdValue = 0.01 detectConfig.detection.background.weighting = False detectTask = DetectCoaddSourcesTask(config=detectConfig) butler = lsst.pipe.tasks.mocks.makeDataRepo(DATAREPO_ROOT) mocksTask.buildAllInputs(butler) addMaskPlanes(butler) mocksTask.buildCoadd(butler) mocksTask.buildMockCoadd(butler) detectTask.writeSchemas(butler) # Now run the seperate multiband tasks on the Coadd to make the reference # catalog for the forced photometry tests. runTaskOnPatches(butler, detectTask, mocksTask) mergeDetConfig = MergeDetectionsTask.ConfigClass() mergeDetConfig.priorityList = ['r', ] mergeDetTask = MergeDetectionsTask(config=mergeDetConfig, butler=butler) mergeDetTask.writeSchemas(butler) runTaskOnPatchList(butler, mergeDetTask, mocksTask) measMergedConfig = MeasureMergedCoaddSourcesTask.ConfigClass() measMergedConfig.measurement.slots.shape = "base_SdssShape" measMergedConfig.measurement.plugins['base_PixelFlags'].masksFpAnywhere = [] measMergedConfig.propagateFlags.flags = {} # Disable flag propagation: no flags to propagate measMergedConfig.doMatchSources = False # We don't have a reference catalog available measMergedTask = MeasureMergedCoaddSourcesTask(config=measMergedConfig, butler=butler) measMergedTask.writeSchemas(butler) runTaskOnPatches(butler, measMergedTask, mocksTask) mergeMeasConfig = MergeMeasurementsTask.ConfigClass() mergeMeasConfig.priorityList = ['r', ] mergeMeasTask = MergeMeasurementsTask(config=mergeMeasConfig, butler=butler) mergeMeasTask.writeSchemas(butler) runTaskOnPatchList(butler, mergeMeasTask, mocksTask) runForcedPhotCoaddTask(butler, mocksTask) runForcedPhotCcdTask(butler)
def setUp(self): # Create a task that creates simulated images and builds a coadd from them self.mocksTask = lsst.pipe.tasks.mocks.MockCoaddTask() # Create an instance of DetectCoaddSourcesTask to measure on the coadd. # There's no noise in these images, so we set a direct-value threshold, # and the background weighting (when using Approximate) to False detectConfig = DetectCoaddSourcesTask.ConfigClass() detectConfig.detection.thresholdType = "value" detectConfig.detection.thresholdValue = 0.01 detectConfig.detection.background.weighting = False self.detectTask = DetectCoaddSourcesTask(config=detectConfig) if REUSE_DATAREPO and os.path.exists(os.path.join(DATAREPO_ROOT, "_mapper")): self.butler = lsst.daf.persistence.Butler(DATAREPO_ROOT) else: self.butler = lsst.pipe.tasks.mocks.makeDataRepo(DATAREPO_ROOT) self.mocksTask.buildAllInputs(self.butler) self.addMaskPlanes() self.checkMaskPlanesExist() self.mocksTask.buildCoadd(self.butler) self.mocksTask.buildMockCoadd(self.butler) self.detectTask.writeSchemas(self.butler) # Now run the seperate multiband tasks on the Coadd to make the reference # catalog for the forced photometry tests. self.runTaskOnPatches(self.detectTask) mergeDetConfig = MergeDetectionsTask.ConfigClass() mergeDetConfig.priorityList = ['r', ] mergeDetTask = MergeDetectionsTask(config=mergeDetConfig, butler=self.butler) mergeDetTask.writeSchemas(self.butler) self.runTaskOnPatchList(mergeDetTask) measMergedConfig = MeasureMergedCoaddSourcesTask.ConfigClass() measMergedConfig.measurement.slots.shape = "base_SdssShape" measMergedConfig.measurement.doApplyApCorr = "no" measMergedConfig.measurement.plugins['base_PixelFlags'].masksFpAnywhere = [] measMergedConfig.propagateFlags.flags = {} # Disable flag propagation: no flags to propagate measMergedConfig.doMatchSources = False # We don't have a reference catalog available measMergedTask = MeasureMergedCoaddSourcesTask(config=measMergedConfig, butler=self.butler) measMergedTask.writeSchemas(self.butler) self.runTaskOnPatches(measMergedTask) mergeMeasConfig = MergeMeasurementsTask.ConfigClass() mergeMeasConfig.priorityList = ['r', ] mergeMeasTask = MergeMeasurementsTask(config=mergeMeasConfig, butler=self.butler) mergeMeasTask.writeSchemas(self.butler) self.runTaskOnPatchList(mergeMeasTask)
def setup_module(module): lsst.utils.tests.init() if os.path.exists(DATAREPO_ROOT): print("Deleting existing repo: %r" % (DATAREPO_ROOT,)) shutil.rmtree(DATAREPO_ROOT) if not haveMeasBase: raise unittest.SkipTest("meas_base could not be imported") # Create a task that creates simulated images and builds a coadd from them mocksTask = lsst.pipe.tasks.mocks.MockCoaddTask() # Create an instance of DetectCoaddSourcesTask to measure on the coadd. # There's no noise in these images, so we set a direct-value threshold, # and the background weighting (when using Approximate) to False detectConfig = DetectCoaddSourcesTask.ConfigClass() # Images have no noise, so we can't use the default DynamicDetectionTask detectConfig.detection.retarget(lsst.meas.algorithms.SourceDetectionTask) detectConfig.detection.thresholdType = "value" detectConfig.detection.thresholdValue = 0.01 detectConfig.detection.background.weighting = False detectTask = DetectCoaddSourcesTask(config=detectConfig) butler = lsst.pipe.tasks.mocks.makeDataRepo(DATAREPO_ROOT) mocksTask.buildAllInputs(butler) addMaskPlanes(butler) mocksTask.buildCoadd(butler) mocksTask.buildMockCoadd(butler) detectTask.writeSchemas(butler) # Now run the seperate multiband tasks on the Coadd to make the reference # catalog for the forced photometry tests. runTaskOnPatches(butler, detectTask, mocksTask) mergeDetConfig = MergeDetectionsTask.ConfigClass() mergeDetConfig.priorityList = ['r', ] mergeDetTask = MergeDetectionsTask(config=mergeDetConfig, butler=butler) mergeDetTask.writeSchemas(butler) runTaskOnPatchList(butler, mergeDetTask, mocksTask) deblendSourcesConfig = DeblendCoaddSourcesTask.ConfigClass() deblendSourcesTask = DeblendCoaddSourcesTask(config=deblendSourcesConfig, butler=butler) deblendSourcesTask.writeSchemas(butler) runTaskOnPatchList(butler, deblendSourcesTask, mocksTask) measMergedConfig = MeasureMergedCoaddSourcesTask.ConfigClass() measMergedConfig.measurement.slots.shape = "base_SdssShape" measMergedConfig.measurement.plugins['base_PixelFlags'].masksFpAnywhere = [] measMergedConfig.propagateFlags.flags = {} # Disable flag propagation: no flags to propagate measMergedConfig.doMatchSources = False # We don't have a reference catalog available measMergedTask = MeasureMergedCoaddSourcesTask(config=measMergedConfig, butler=butler) measMergedTask.writeSchemas(butler) runTaskOnPatches(butler, measMergedTask, mocksTask) mergeMeasConfig = MergeMeasurementsTask.ConfigClass() mergeMeasConfig.priorityList = ['r', ] mergeMeasTask = MergeMeasurementsTask(config=mergeMeasConfig, butler=butler) mergeMeasTask.writeSchemas(butler) runTaskOnPatchList(butler, mergeMeasTask, mocksTask) runForcedPhotCoaddTask(butler, mocksTask) runForcedPhotCcdTask(butler)
#!/usr/bin/env python from lsst.pipe.tasks.multiBand import MergeDetectionsTask MergeDetectionsTask.parseAndRun()
def setUp(self): # Create a task that creates simulated images and builds a coadd from them self.mocksTask = lsst.pipe.tasks.mocks.MockCoaddTask() # Create an instance of DetectCoaddSourcesTask to measure on the coadd. # There's no noise in these images, so we set a direct-value threshold, # and the background weighting (when using Approximate) to False detectConfig = DetectCoaddSourcesTask.ConfigClass() detectConfig.detection.thresholdType = "value" detectConfig.detection.thresholdValue = 0.01 detectConfig.detection.background.weighting = False self.detectTask = DetectCoaddSourcesTask(config=detectConfig) if REUSE_DATAREPO and os.path.exists( os.path.join(DATAREPO_ROOT, "_mapper")): self.butler = lsst.daf.persistence.Butler(DATAREPO_ROOT) else: self.butler = lsst.pipe.tasks.mocks.makeDataRepo(DATAREPO_ROOT) self.mocksTask.buildAllInputs(self.butler) self.addMaskPlanes() self.checkMaskPlanesExist() self.mocksTask.buildCoadd(self.butler) self.mocksTask.buildMockCoadd(self.butler) self.detectTask.writeSchemas(self.butler) # Now run the seperate multiband tasks on the Coadd to make the reference # catalog for the forced photometry tests. self.runTaskOnPatches(self.detectTask) mergeDetConfig = MergeDetectionsTask.ConfigClass() mergeDetConfig.priorityList = [ 'r', ] mergeDetTask = MergeDetectionsTask(config=mergeDetConfig, butler=self.butler) mergeDetTask.writeSchemas(self.butler) self.runTaskOnPatchList(mergeDetTask) measMergedConfig = MeasureMergedCoaddSourcesTask.ConfigClass() measMergedConfig.measurement.slots.shape = "base_SdssShape" measMergedConfig.measurement.doApplyApCorr = "no" measMergedConfig.measurement.plugins[ 'base_PixelFlags'].masksFpAnywhere = [] measMergedConfig.propagateFlags.flags = { } # Disable flag propagation: no flags to propagate measMergedConfig.doMatchSources = False # We don't have a reference catalog available measMergedTask = MeasureMergedCoaddSourcesTask( config=measMergedConfig, butler=self.butler) measMergedTask.writeSchemas(self.butler) self.runTaskOnPatches(measMergedTask) mergeMeasConfig = MergeMeasurementsTask.ConfigClass() mergeMeasConfig.priorityList = [ 'r', ] mergeMeasTask = MergeMeasurementsTask(config=mergeMeasConfig, butler=self.butler) mergeMeasTask.writeSchemas(self.butler) self.runTaskOnPatchList(mergeMeasTask)
def setUpClass(cls): """Create 200MB of test data.""" # Start by assuming nothing failed cls.failed = False if os.path.exists(DATAREPO_ROOT): print(f"Deleting existing repo: {DATAREPO_ROOT}") # Do not ignore errors since failure to clean up is indicative shutil.rmtree(DATAREPO_ROOT) # Create a task that creates simulated images and builds a coadd from them mocksTask = lsst.pipe.tasks.mocks.MockCoaddTask() # Create an instance of DetectCoaddSourcesTask to measure on the coadd. # There's no noise in these images, so we set a direct-value threshold, # and the background weighting (when using Approximate) to False detectConfig = DetectCoaddSourcesTask.ConfigClass() # Images have no noise, so we can't use the default DynamicDetectionTask detectConfig.detection.retarget( lsst.meas.algorithms.SourceDetectionTask) detectConfig.detection.thresholdType = "value" detectConfig.detection.thresholdValue = 0.01 detectConfig.detection.background.weighting = False detectTask = DetectCoaddSourcesTask(config=detectConfig) butler = lsst.pipe.tasks.mocks.makeDataRepo(DATAREPO_ROOT) mocksTask.buildAllInputs(butler) addMaskPlanes(butler) mocksTask.buildCoadd(butler) mocksTask.buildMockCoadd(butler) detectTask.writeSchemas(butler) # Now run the seperate multiband tasks on the Coadd to make the reference # catalog for the forced photometry tests. runTaskOnPatches(butler, detectTask, mocksTask) mergeDetConfig = MergeDetectionsTask.ConfigClass() mergeDetConfig.priorityList = [ 'r', ] mergeDetTask = MergeDetectionsTask(config=mergeDetConfig, butler=butler) mergeDetTask.writeSchemas(butler) runTaskOnPatchList(butler, mergeDetTask, mocksTask) deblendSourcesConfig = DeblendCoaddSourcesTask.ConfigClass() deblendSourcesTask = DeblendCoaddSourcesTask( config=deblendSourcesConfig, butler=butler) deblendSourcesTask.writeSchemas(butler) runTaskOnPatchList(butler, deblendSourcesTask, mocksTask) measMergedConfig = MeasureMergedCoaddSourcesTask.ConfigClass() measMergedConfig.measurement.slots.shape = "base_SdssShape" measMergedConfig.measurement.plugins[ 'base_PixelFlags'].masksFpAnywhere = [] measMergedConfig.propagateFlags.flags = { } # Disable flag propagation: no flags to propagate measMergedConfig.doMatchSources = False # We don't have a reference catalog available measMergedTask = MeasureMergedCoaddSourcesTask(config=measMergedConfig, butler=butler) measMergedTask.writeSchemas(butler) runTaskOnPatches(butler, measMergedTask, mocksTask) mergeMeasConfig = MergeMeasurementsTask.ConfigClass() mergeMeasConfig.priorityList = [ 'r', ] mergeMeasTask = MergeMeasurementsTask(config=mergeMeasConfig, butler=butler) mergeMeasTask.writeSchemas(butler) runTaskOnPatchList(butler, mergeMeasTask, mocksTask) runForcedPhotCoaddTask(butler, mocksTask) runForcedPhotCcdTask(butler)