LSSTApplications  18.1.0
LSSTDataManagementBasePackage
makeMatchStatistics.cc
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 
3 /*
4  * LSST Data Management System
5  * Copyright 2008, 2009, 2010 LSST Corporation.
6  *
7  * This product includes software developed by the
8  * LSST Project (http://www.lsst.org/).
9  *
10  * This program is free software: you can redistribute it and/or modify
11  * it under the terms of the GNU General Public License as published by
12  * the Free Software Foundation, either version 3 of the License, or
13  * (at your option) any later version.
14  *
15  * This program is distributed in the hope that it will be useful,
16  * but WITHOUT ANY WARRANTY; without even the implied warranty of
17  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18  * GNU General Public License for more details.
19  *
20  * You should have received a copy of the LSST License Statement and
21  * the GNU General Public License along with this program. If not,
22  * see <http://www.lsstcorp.org/LegalNotices/>.
23  */
24 
27 
28 namespace lsst {
29 namespace meas {
30 namespace astrom {
31 
32 template <typename MatchT>
34  afw::math::StatisticsControl const& sctrl) {
35  if (matchList.empty()) {
36  throw LSST_EXCEPT(pexExcept::RuntimeError, "matchList is empty");
37  }
39  val.reserve(matchList.size());
40 
41  for (auto const& match : matchList) {
42  val.push_back(match.distance);
43  }
44  return afw::math::makeStatistics(val, flags, sctrl);
45 }
46 
47 template <typename MatchT>
49  std::vector<MatchT> const& matchList, int const flags,
50  afw::math::StatisticsControl const& sctrl) {
51  if (matchList.empty()) {
52  throw LSST_EXCEPT(pexExcept::RuntimeError, "matchList is empty");
53  }
55  val.reserve(matchList.size());
56 
57  for (auto const& match : matchList) {
58  auto refPtr = match.first;
59  auto srcPtr = match.second;
60  auto srcX = srcPtr->getX();
61  auto srcY = srcPtr->getY();
62  auto refPos = wcs.skyToPixel(refPtr->getCoord());
63  auto refX = refPos[0];
64  auto refY = refPos[1];
65  val.push_back(::hypot(srcX - refX, srcY - refY));
66  }
67  return afw::math::makeStatistics(val, flags, sctrl);
68 }
69 
70 template <typename MatchT>
72  std::vector<MatchT> const& matchList, int const flags,
73  afw::math::StatisticsControl const& sctrl) {
74  if (matchList.empty()) {
75  throw LSST_EXCEPT(pexExcept::RuntimeError, "matchList is empty");
76  }
78  val.reserve(matchList.size());
79 
80  for (auto const& match : matchList) {
81  auto refPtr = match.first;
82  auto srcPtr = match.second;
83  auto refCoord = refPtr->getCoord();
84  auto srcCoord = wcs.pixelToSky(srcPtr->getCentroid());
85  auto angSep = refCoord.separation(srcCoord);
86  val.push_back(angSep.asRadians());
87  }
88  return afw::math::makeStatistics(val, flags, sctrl);
89 }
90 
91 #define INSTANTIATE(MATCH) \
92  template afw::math::Statistics makeMatchStatistics<MATCH>(std::vector<MATCH> const& matchList, \
93  int const flags, \
94  afw::math::StatisticsControl const& sctrl); \
95  template afw::math::Statistics makeMatchStatisticsInPixels<MATCH>( \
96  afw::geom::SkyWcs const& wcs, std::vector<MATCH> const& matchList, int const flags, \
97  afw::math::StatisticsControl const& sctrl); \
98  template afw::math::Statistics makeMatchStatisticsInRadians<MATCH>( \
99  afw::geom::SkyWcs const& wcs, std::vector<MATCH> const& matchList, int const flags, \
100  afw::math::StatisticsControl const& sctrl);
101 
104 
105 } // namespace astrom
106 } // namespace meas
107 } // namespace lsst
afw::math::Statistics makeMatchStatistics(std::vector< MatchT > const &matchList, int const flags, afw::math::StatisticsControl const &sctrl=afw::math::StatisticsControl())
Compute statistics of the distance field of a match list.
T empty(T... args)
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition: SkyWcs.h:117
lsst::geom::Point2D skyToPixel(lsst::geom::SpherePoint const &sky) const
Compute pixel position(s) from sky position(s)
Definition: SkyWcs.h:349
lsst::geom::SpherePoint pixelToSky(lsst::geom::Point2D const &pixel) const
Compute sky position(s) from pixel position(s)
Definition: SkyWcs.h:334
ImageT val
Definition: CR.cc:146
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:71
A class to evaluate image statistics.
Definition: Statistics.h:215
T push_back(T... args)
Pass parameters to a Statistics object.
Definition: Statistics.h:93
A base class for image defects.
afw::math::Statistics makeMatchStatisticsInPixels(afw::geom::SkyWcs const &wcs, std::vector< MatchT > const &matchList, int const flags, afw::math::StatisticsControl const &sctrl=afw::math::StatisticsControl())
Compute statistics of on-detector radial separation for a match list, in pixels.
Statistics makeStatistics(lsst::afw::image::Image< Pixel > const &img, lsst::afw::image::Mask< image::MaskPixel > const &msk, int const flags, StatisticsControl const &sctrl=StatisticsControl())
Handle a watered-down front-end to the constructor (no variance)
Definition: Statistics.h:354
Match< SimpleRecord, SourceRecord > ReferenceMatch
Definition: fwd.h:104
Angle separation(SpherePoint const &other) const noexcept
Angular distance between two points.
Definition: SpherePoint.cc:173
#define INSTANTIATE(MATCH)
T size(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
afw::math::Statistics makeMatchStatisticsInRadians(afw::geom::SkyWcs const &wcs, std::vector< MatchT > const &matchList, int const flags, afw::math::StatisticsControl const &sctrl=afw::math::StatisticsControl())
Compute statistics of on-sky radial separation for a match list, in radians.
Match< SourceRecord, SourceRecord > SourceMatch
Definition: fwd.h:105
T reserve(T... args)
Reports errors that are due to events beyond the control of the program.
Definition: Runtime.h:104