Coverage for tests/test_SdssShape.py : 25%

Hot-keys on this page
r m x p toggle line displays
j k next/prev highlighted chunk
0 (zero) top of page
1 (one) first highlighted chunk
# This file is part of meas_base. # # Developed for the LSST Data Management System. # This product includes software developed by the LSST Project # (https://www.lsst.org). # See the COPYRIGHT file at the top-level directory of this distribution # for details of code ownership. # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by # the Free Software Foundation, either version 3 of the License, or # (at your option) any later version. # # This program is distributed in the hope that it will be useful, # but WITHOUT ANY WARRANTY; without even the implied warranty of # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the # GNU General Public License for more details. # # You should have received a copy of the GNU General Public License # along with this program. If not, see <https://www.gnu.org/licenses/>.
self.bbox = lsst.geom.Box2I(lsst.geom.Point2I(-20, -30), lsst.geom.Extent2I(240, 160)) self.dataset = lsst.meas.base.tests.TestDataset(self.bbox) # first source is a point self.dataset.addSource(100000.0, lsst.geom.Point2D(50.1, 49.8)) # second source is extended self.dataset.addSource(100000.0, lsst.geom.Point2D(149.9, 50.3), lsst.afw.geom.Quadrupole(8, 9, 3)) self.config = self.makeSingleFrameMeasurementConfig("base_SdssShape")
del self.bbox del self.dataset del self.config
self.assertTrue(np.isfinite(value), msg="%s is not finite" % (value,))
task = self.makeSingleFrameMeasurementTask("base_SdssShape", config=self.config) exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=0) task.run(catalog, exposure) return exposure, catalog
self.assertFloatsAlmostEqual(result.x, record.get("truth_x"), rtol=1E-2) self.assertFloatsAlmostEqual(result.y, record.get("truth_y"), rtol=1E-2) self.assertFloatsAlmostEqual(result.xx, record.get("truth_xx"), rtol=1E-2) self.assertFloatsAlmostEqual(result.yy, record.get("truth_yy"), rtol=1E-2) self.assertFloatsAlmostEqual(result.xy, record.get("truth_xy"), rtol=1E-1, atol=2E-1) self.assertFinite(result.xxErr) self.assertFinite(result.yyErr) self.assertFinite(result.xyErr) self.assertFinite(result.instFlux_xx_Cov) self.assertFinite(result.instFlux_yy_Cov) self.assertFinite(result.instFlux_xy_Cov) self.assertFinite(result.xx_yy_Cov) self.assertFinite(result.xx_xy_Cov) self.assertFinite(result.yy_xy_Cov) self.assertFalse(result.getFlag(lsst.meas.base.SdssShapeAlgorithm.FAILURE.number)) self.assertFalse(result.getFlag(lsst.meas.base.SdssShapeAlgorithm.UNWEIGHTED_BAD.number)) self.assertFalse(result.getFlag(lsst.meas.base.SdssShapeAlgorithm.UNWEIGHTED.number)) self.assertFalse(result.getFlag(lsst.meas.base.SdssShapeAlgorithm.SHIFT.number)) self.assertFalse(result.getFlag(lsst.meas.base.SdssShapeAlgorithm.MAXITER.number))
self.assertFloatsAlmostEqual(psfResult.getIxx(), psfTruth.getIxx(), rtol=1E-4) self.assertFloatsAlmostEqual(psfResult.getIyy(), psfTruth.getIyy(), rtol=1E-4) self.assertFloatsAlmostEqual(psfResult.getIxy(), psfTruth.getIxy(), rtol=1E-4) self.assertFalse(result.getFlag(lsst.meas.base.SdssShapeAlgorithm.PSF_SHAPE_BAD.number))
"""Test that we measure shapes and record the PSF shape correctly
Note: Given that the PSF model here is constant over the entire image, this test would not catch an error in the potition at which base_SdssShape_psf is computed. Such a test requires a spatially varying PSF model such that different locations can be distinguished by their different PSF model shapes. Such a test exists in meas_algorithms (tests/testSdssShapePsf.py), making use of the PcaPsf algorithm to build the spatially varying PSF. """ exposure, catalog = self._runMeasurementTask() key = lsst.meas.base.SdssShapeResultKey(catalog.schema["base_SdssShape"]) psfTruth = exposure.getPsf().computeShape() for record in catalog: result = record.get(key) self._checkShape(result, record) psfResult = key.getPsfShape(record) self._checkPsfShape(result, psfResult, psfTruth)
"""Test that we measure shapes correctly and do not record the PSF shape when not needed.""" self.config.plugins["base_SdssShape"].doMeasurePsf = False _, catalog = self._runMeasurementTask() key = lsst.meas.base.SdssShapeResultKey(catalog.schema["base_SdssShape"]) for record in catalog: result = record.get(key) self._checkShape(result, record) self.assertNotIn("base_SdssShape_psf_xx", catalog.schema) self.assertNotIn("base_SdssShape_psf_yy", catalog.schema) self.assertNotIn("base_SdssShape_psf_xy", catalog.schema) self.assertNotIn("base_SdssShape_flag_psf", catalog.schema)
"""Test that we measure shapes correctly and set a flag with the PSF is unavailable.""" self.config.plugins["base_SdssShape"].doMeasurePsf = True task = self.makeSingleFrameMeasurementTask("base_SdssShape", config=self.config) exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=1) exposure.setPsf(None) # Set PSF to None to test no PSF case task.run(catalog, exposure) key = lsst.meas.base.SdssShapeResultKey(catalog.schema["base_SdssShape"]) for record in catalog: result = record.get(key) self._checkShape(result, record) self.assertTrue(result.getFlag(lsst.meas.base.SdssShapeAlgorithm.PSF_SHAPE_BAD.number))
lsst.meas.base.tests.CentroidTransformTestCase, lsst.meas.base.tests.SingleFramePluginTransformSetupHelper, lsst.utils.tests.TestCase):
lsst.meas.base.tests.FluxTransformTestCase._setFieldsInRecords(self, records, name) lsst.meas.base.tests.CentroidTransformTestCase._setFieldsInRecords(self, records, name) for record in records: for field in ('xx', 'yy', 'xy', 'xxErr', 'yyErr', 'xyErr', 'psf_xx', 'psf_yy', 'psf_xy'): if record.schema.join(name, field) in record.schema: record[record.schema.join(name, field)] = np.random.random()
lsst.meas.base.tests.FluxTransformTestCase._compareFieldsInRecords(self, inSrc, outSrc, name) lsst.meas.base.tests.CentroidTransformTestCase._compareFieldsInRecords(self, inSrc, outSrc, name)
inShape = lsst.meas.base.ShapeResultKey(inSrc.schema[name]).get(inSrc) outShape = lsst.meas.base.ShapeResultKey(outSrc.schema[name]).get(outSrc)
centroid = lsst.meas.base.CentroidResultKey(inSrc.schema[name]).get(inSrc).getCentroid() xform = self.calexp.getWcs().linearizePixelToSky(centroid, lsst.geom.radians)
trInShape = inShape.getShape().transform(xform.getLinear()) self.assertEqual(trInShape.getIxx(), outShape.getShape().getIxx()) self.assertEqual(trInShape.getIyy(), outShape.getShape().getIyy()) self.assertEqual(trInShape.getIxy(), outShape.getShape().getIxy())
m = lsst.meas.base.makeShapeTransformMatrix(xform.getLinear()) np.testing.assert_array_almost_equal( np.dot(np.dot(m, inShape.getShapeErr()), m.transpose()), outShape.getShapeErr() ) if self.testPsf: inPsfShape = lsst.meas.base.ShapeResultKey( inSrc.schema[inSrc.schema.join(name, "psf")] ).get(inSrc) outPsfShape = lsst.meas.base.ShapeResultKey( outSrc.schema[outSrc.schema.join(name, "psf")] ).get(outSrc) trInPsfShape = inPsfShape.getShape().transform(xform.getLinear())
self.assertEqual(trInPsfShape.getIxx(), outPsfShape.getShape().getIxx()) self.assertEqual(trInPsfShape.getIyy(), outPsfShape.getShape().getIyy()) self.assertEqual(trInPsfShape.getIxy(), outPsfShape.getShape().getIxy()) else: self.assertNotIn(inSrc.schema.join(name, "psf", "xx"), inSrc.schema) self.assertNotIn(inSrc.schema.join(name, "flag", "psf"), inSrc.schema)
ctrl = lsst.meas.base.SdssShapeControl() ctrl.doMeasurePsf = False return ctrl
lsst.utils.tests.init()
lsst.utils.tests.init() unittest.main() |