LSSTApplications  21.0.0+1b62c9342b,21.0.0+45a059f35e,21.0.0-1-ga51b5d4+ceb9cf20a3,21.0.0-10-g68cce58c5+c7d3cce47e,21.0.0-2-g103fe59+c1ca725317,21.0.0-2-g1367e85+a1c2f7fe71,21.0.0-2-g2909d54+45a059f35e,21.0.0-2-g45278ab+1b62c9342b,21.0.0-2-g4bc9b9f+b2e40a4e47,21.0.0-2-g5242d73+a1c2f7fe71,21.0.0-2-g54e2caa+c00cf99ed0,21.0.0-2-g66bcc37+27b9d7859a,21.0.0-2-g7f82c8f+203cf74700,21.0.0-2-g8dde007+b0df52bfdd,21.0.0-2-g8f08a60+73884b2cf5,21.0.0-2-ga326454+203cf74700,21.0.0-2-ga63a54e+eec04437aa,21.0.0-2-gc738bc1+59028256f4,21.0.0-2-gde069b7+5a8f2956b8,21.0.0-2-ge17e5af+a1c2f7fe71,21.0.0-2-ge712728+9ad031c87e,21.0.0-2-gecfae73+d3766aec80,21.0.0-2-gfc62afb+a1c2f7fe71,21.0.0-20-g4449a12+38dfb87bce,21.0.0-22-gf0532904+afb8e7912b,21.0.0-3-g4c5b185+a403cb96fd,21.0.0-3-g6d51c4a+27b9d7859a,21.0.0-3-g8076721+e873df194c,21.0.0-3-gaa929c8+df5d87f43a,21.0.0-3-gd222c45+afc8332dbe,21.0.0-4-g1383c07+27b9d7859a,21.0.0-4-g3300ddd+1b62c9342b,21.0.0-4-g5873dc9+9a92674037,21.0.0-4-g8a80011+f67daf2f53,21.0.0-5-gcff38f6+bce43c5818,21.0.0-6-g463d161+44134145d4,21.0.0-6-gd3283ba+df5d87f43a,21.0.0-8-g19111d86+d6551531e4,w.2021.04
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
makeMatchStatistics.h
wcs
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:71
std::vector< MatchT >
std::vector::size
T size(T... args)
lsst::afw::geom::SkyWcs
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition: SkyWcs.h:117
val
ImageT val
Definition: CR.cc:146
lsst::meas::astrom::makeMatchStatistics
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.
Definition: makeMatchStatistics.cc:33
lsst::meas::astrom::makeMatchStatisticsInRadians
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.
Definition: makeMatchStatistics.cc:71
lsst::afw::math::makeStatistics
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
lsst::meas::astrom::makeMatchStatisticsInPixels
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.
Definition: makeMatchStatistics.cc:48
lsst::afw::table::Match
Lightweight representation of a geometric match between two records.
Definition: Match.h:67
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
LSST_EXCEPT
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
lsst::afw::math::StatisticsControl
Pass parameters to a Statistics object.
Definition: Statistics.h:93
Runtime.h
INSTANTIATE
#define INSTANTIATE(MATCH)
Definition: makeMatchStatistics.cc:91
std::vector::empty
T empty(T... args)
lsst::afw::math::Statistics
A class to evaluate image statistics.
Definition: Statistics.h:215
lsst.pex::exceptions::RuntimeError
Reports errors that are due to events beyond the control of the program.
Definition: Runtime.h:104