LSSTApplications  19.0.0-14-gb0260a2+72efe9b372,20.0.0+7927753e06,20.0.0+8829bf0056,20.0.0+995114c5d2,20.0.0+b6f4b2abd1,20.0.0+bddc4f4cbe,20.0.0-1-g253301a+8829bf0056,20.0.0-1-g2b7511a+0d71a2d77f,20.0.0-1-g5b95a8c+7461dd0434,20.0.0-12-g321c96ea+23efe4bbff,20.0.0-16-gfab17e72e+fdf35455f6,20.0.0-2-g0070d88+ba3ffc8f0b,20.0.0-2-g4dae9ad+ee58a624b3,20.0.0-2-g61b8584+5d3db074ba,20.0.0-2-gb780d76+d529cf1a41,20.0.0-2-ged6426c+226a441f5f,20.0.0-2-gf072044+8829bf0056,20.0.0-2-gf1f7952+ee58a624b3,20.0.0-20-geae50cf+e37fec0aee,20.0.0-25-g3dcad98+544a109665,20.0.0-25-g5eafb0f+ee58a624b3,20.0.0-27-g64178ef+f1f297b00a,20.0.0-3-g4cc78c6+e0676b0dc8,20.0.0-3-g8f21e14+4fd2c12c9a,20.0.0-3-gbd60e8c+187b78b4b8,20.0.0-3-gbecbe05+48431fa087,20.0.0-38-ge4adf513+a12e1f8e37,20.0.0-4-g97dc21a+544a109665,20.0.0-4-gb4befbc+087873070b,20.0.0-4-gf910f65+5d3db074ba,20.0.0-5-gdfe0fee+199202a608,20.0.0-5-gfbfe500+d529cf1a41,20.0.0-6-g64f541c+d529cf1a41,20.0.0-6-g9a5b7a1+a1cd37312e,20.0.0-68-ga3f3dda+5fca18c6a4,20.0.0-9-g4aef684+e18322736b,w.2020.45
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: fwd.h:101
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
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