lsst.afw  22.0.1-32-g5ddfab5d3+23bd69c089
ApCorrMap.cc
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008-2014 LSST Corporation.
5  *
6  * This product includes software developed by the
7  * LSST Project (http://www.lsst.org/).
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 LSST License Statement and
20  * the GNU General Public License along with this program. If not,
21  * see <http://www.lsstcorp.org/LegalNotices/>.
22  */
23 
24 #include <memory>
25 
31 
32 namespace lsst {
33 namespace afw {
34 
37 
38 namespace image {
39 
40 // Even though this static const member is set in the header, it needs to be declared here if we need
41 // to take its address (and Swig might).
42 std::size_t const ApCorrMap::MAX_NAME_LENGTH;
43 
44 std::shared_ptr<math::BoundedField> const ApCorrMap::operator[](std::string const& name) const {
45  Iterator i = _internal.find(name);
46  if (i == _internal.end()) {
48  (boost::format("Aperture correction with name '%s' not found") % name).str());
49  }
50  return i->second;
51 }
52 
53 std::shared_ptr<math::BoundedField> const ApCorrMap::get(std::string const& name) const {
54  Iterator i = _internal.find(name);
55  if (i == _internal.end()) {
57  }
58  return i->second;
59 }
60 
62  if (name.size() > MAX_NAME_LENGTH) {
63  throw LSST_EXCEPT(
65  (boost::format("Aperture correction name '%s' exceeds size limit of %d characters") % name %
66  MAX_NAME_LENGTH)
67  .str());
68  }
69  _internal.insert(std::make_pair(name, field));
70 }
71 
72 namespace {
73 
74 struct PersistenceHelper {
75  table::Schema schema;
76  table::Key<std::string> name;
77  table::Key<int> field;
78 
79  static PersistenceHelper const& get() {
80  static PersistenceHelper const instance;
81  return instance;
82  }
83 
84 private:
85  PersistenceHelper()
86  : schema(),
87  name(schema.addField<std::string>("name", "name of the aperture correction",
88  ApCorrMap::MAX_NAME_LENGTH)),
89  field(schema.addField<int>("field", "archive ID of the BoundedField object")) {}
90 };
91 
92 class ApCorrMapFactory : public table::io::PersistableFactory {
93 public:
94  std::shared_ptr<table::io::Persistable> read(InputArchive const& archive,
95  CatalogVector const& catalogs) const override {
96  PersistenceHelper const& keys = PersistenceHelper::get();
97  LSST_ARCHIVE_ASSERT(catalogs.size() == 1u);
98  LSST_ARCHIVE_ASSERT(catalogs.front().getSchema() == keys.schema);
99  std::shared_ptr<ApCorrMap> result = std::make_shared<ApCorrMap>();
100  for (table::BaseCatalog::const_iterator i = catalogs.front().begin(); i != catalogs.front().end();
101  ++i) {
102  result->set(i->get(keys.name), archive.get<math::BoundedField>(i->get(keys.field)));
103  }
104  return result;
105  }
106 
107  ApCorrMapFactory(std::string const& name) : afw::table::io::PersistableFactory(name) {}
108 };
109 
110 std::string getApCorrMapPersistenceName() { return "ApCorrMap"; }
111 
112 ApCorrMapFactory registration(getApCorrMapPersistenceName());
113 
114 } // namespace
115 
116 bool ApCorrMap::isPersistable() const noexcept {
117  for (auto const &i : *this) {
118  if (!i.second->isPersistable()) return false;
119  }
120  return true;
121 }
122 
123 std::string ApCorrMap::getPersistenceName() const { return getApCorrMapPersistenceName(); }
124 
125 std::string ApCorrMap::getPythonModule() const { return "lsst.afw.image"; }
126 
127 void ApCorrMap::write(OutputArchiveHandle& handle) const {
128  PersistenceHelper const& keys = PersistenceHelper::get();
129  table::BaseCatalog catalog = handle.makeCatalog(keys.schema);
130  for (auto const &i : *this) {
131  std::shared_ptr<table::BaseRecord> record = catalog.addNew();
132  record->set(keys.name, i.first);
133  record->set(keys.field, handle.put(i.second));
134  }
135  handle.saveCatalog(catalog);
136 }
137 
138 ApCorrMap& ApCorrMap::operator*=(double const scale) {
139  Internal replacement;
140  for (auto const &i : *this) {
141  replacement[i.first] = (*i.second) * scale;
142  }
143  _internal = replacement;
144  return *this;
145 }
146 
147 std::shared_ptr<typehandling::Storable> ApCorrMap::cloneStorable() const {
148  return std::make_unique<ApCorrMap>(*this);
149 }
150 
151 } // namespace image
152 } // namespace afw
153 } // namespace lsst
table::Key< std::string > name
Definition: ApCorrMap.cc:76
table::Schema schema
Definition: ApCorrMap.cc:75
table::Key< int > field
Definition: ApCorrMap.cc:77
#define LSST_EXCEPT(type,...)
#define LSST_ARCHIVE_ASSERT(EXPR)
An assertion macro used to validate the structure of an InputArchive.
Definition: Persistable.h:48
A thin wrapper around std::map to allow aperture corrections to be attached to Exposures.
Definition: ApCorrMap.h:45
Internal::const_iterator Iterator
Iterator type returned by begin() and end().
Definition: ApCorrMap.h:54
CatalogIterator< typename Internal::const_iterator > const_iterator
Definition: Catalog.h:112
static std::shared_ptr< T > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18
daf::base::PropertySet * set
Definition: fits.cc:912
T make_pair(T... args)
def keys(self)
def scale(algorithm, min, max=None, frame=None)
Definition: ds9.py:108
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
std::string getPythonModule() const override
std::string getPersistenceName() const override
void write(OutputArchiveHandle &handle) const override
CatalogT< BaseRecord > BaseCatalog
Definition: fwd.h:72
A base class for image defects.
STL namespace.
std::shared_ptr< table::io::Persistable > read(table::io::InputArchive const &archive, table::io::CatalogVector const &catalogs) const override
Definition: warpExposure.cc:0