Coverage for python/lsst/faro/measurement/example_VisitMeasurementTasks.py: 45%
18 statements
« prev ^ index » next coverage.py v6.5.0, created at 2022-11-05 08:10 +0000
« prev ^ index » next coverage.py v6.5.0, created at 2022-11-05 08:10 +0000
1# This file is part of faro.
2#
3# Developed for the LSST Data Management System.
4# This product includes software developed by the LSST Project
5# (https://www.lsst.org).
6# See the COPYRIGHT file at the top-level directory of this distribution
7# for details of code ownership.
8#
9# This program is free software: you can redistribute it and/or modify
10# it under the terms of the GNU General Public License as published by
11# the Free Software Foundation, either version 3 of the License, or
12# (at your option) any later version.
13#
14# This program is distributed in the hope that it will be useful,
15# but WITHOUT ANY WARRANTY; without even the implied warranty of
16# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17# GNU General Public License for more details.
18#
19# You should have received a copy of the GNU General Public License
20# along with this program. If not, see <https://www.gnu.org/licenses/>.
22import astropy.units as u
23from lsst.pipe.base import Struct, Task
24from lsst.pex.config import Config
25from lsst.verify import Measurement
27__all__ = ("StarFracTask",)
30class StarFracTask(Task):
31 ConfigClass = Config
32 _DefaultName = "starFracTask"
34 def run(self, metricName, catalog, vIds):
35 self.log.info("Measuring %s", metricName)
36 if not catalog.isContiguous():
37 catalog = catalog.copy(deep=True)
38 extended = catalog.get("base_ClassificationExtendedness_value")
39 good_extended = extended[~catalog.get("base_ClassificationExtendedness_flag")]
40 n_gals = sum(good_extended)
41 frac = 100 * (len(good_extended) - n_gals) / len(good_extended)
42 meas = Measurement("starFrac", frac * u.percent)
43 return Struct(measurement=meas)