LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
MaskDict.cc
Go to the documentation of this file.
1 /*
2  * Developed for the LSST Data Management System.
3  * This product includes software developed by the LSST Project
4  * (https://www.lsst.org).
5  * See the COPYRIGHT file at the top-level directory of this distribution
6  * for details of code ownership.
7  *
8  * This program is free software: you can redistribute it and/or modify
9  * it under the terms of the GNU General Public License as published by
10  * the Free Software Foundation, either version 3 of the License, or
11  * (at your option) any later version.
12  *
13  * This program is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16  * GNU General Public License for more details.
17  *
18  * You should have received a copy of the GNU General Public License
19  * along with this program. If not, see <https://www.gnu.org/licenses/>.
20  */
21 
22 #include <mutex>
23 #include <set>
24 
25 #include "boost/functional/hash.hpp"
27 
28 namespace lsst { namespace afw { namespace image { namespace detail {
29 
30 // A thread-safe singleton that manages MaskDict's global state.
32 public:
33 
34  static GlobalState & get() {
35  static GlobalState instance;
36  return instance;
37  }
38 
41  if (!mpd.empty()) {
42  std::shared_ptr<MaskDict> dict(new MaskDict(mpd));
43  _allMaskDicts.insert(dict);
44  return dict;
45  }
46  return _defaultMaskDict;
47  }
48 
52  _allMaskDicts.insert(result);
53  return result;
54  }
55 
57  return _defaultMaskDict;
58  }
59 
61  _defaultMaskDict = std::move(dict);
62  }
63 
66  _defaultMaskDict = copy(*_defaultMaskDict);
67  return _defaultMaskDict;
68  }
69 
70  template <typename Functor>
71  void forEachMaskDict(Functor functor) {
73  _prune(); // guarantees dereference below is safe
74  for (auto const & ptr : _allMaskDicts) {
75  functor(*ptr.lock());
76  }
77  }
78 
79 private:
80 
81  GlobalState() : _defaultMaskDict(new MaskDict()) {
82  _allMaskDicts.insert(_defaultMaskDict);
83  _defaultMaskDict->_addInitialMaskPlanes();
84  }
85 
86  GlobalState(GlobalState const &) = delete;
87  GlobalState(GlobalState &&) = delete;
88 
89  GlobalState & operator=(GlobalState const &) = delete;
90  GlobalState & operator=(GlobalState &&) = delete;
91 
92  ~GlobalState() = default;
93 
94  // Prune expired weak_ptrs from _allMaskDicts. Not thread safe; should
95  // only be called by routines that have already locked the mutex.
96  void _prune() {
97  for (auto iter = _allMaskDicts.begin(); iter != _allMaskDicts.end(); ) {
98  if (iter->expired()) {
99  iter = _allMaskDicts.erase(iter);
100  } else {
101  ++iter;
102  }
103  }
104  }
105 
106  std::recursive_mutex _mutex; // guards _allMaskDicts and synchronizes updates to it and _defaultMaskDict
108  std::shared_ptr<MaskDict> _defaultMaskDict;
109 };
110 
112  return GlobalState::get().copyOrGetDefault(mpd);
113 }
114 
116  return GlobalState::get().getDefault();
117 }
118 
121 }
122 
124  return GlobalState::get().detachDefault();
125 }
126 
127 void MaskDict::addAllMasksPlane(std::string const & name, int bitId) {
129  [&name, bitId](MaskDict & dict) {
130  auto const found = std::find_if(
131  dict.begin(), dict.end(),
132  [bitId](auto const & item) { return item.second == bitId; }
133  );
134  if (found == dict.end()) {
135  // is name already in use?
136  if (dict.find(name) == dict.end()) {
137  dict.add(name, bitId);
138  }
139  }
140  }
141  );
142 }
143 
144 MaskDict::~MaskDict() noexcept = default;
145 
146 std::shared_ptr<MaskDict> MaskDict::clone() const {
147  return GlobalState::get().copy(*this);
148 }
149 
151  if (empty()) {
152  return 0;
153  }
154 
155  auto const maxIter = std::max_element(
156  begin(), end(),
157  [](auto const & a, auto const & b) { return a.second < b.second; }
158  );
159  assert(maxIter != end());
160  int id = maxIter->second + 1; // The maskPlane to use if there are no gaps
161 
162  for (int i = 0; i < id; ++i) {
163  // is i already used in this Mask?
164  auto const sameIter = std::find_if(
165  begin(), end(),
166  [i](auto const & item) { return item.second == i; }
167  );
168  if (sameIter == end()) { // Not used; so we'll use it
169  return i;
170  }
171  }
172 
173  return id;
174 }
175 
177  auto iter = find(name);
178  return (iter == end()) ? -1 : iter->second;
179 }
180 
181 void MaskDict::print() const {
182  for (auto const & item : *this) {
183  std::cout << "Plane " << item.second << " -> " << item.first << std::endl;
184  }
185 }
186 
187 bool MaskDict::operator==(MaskDict const& rhs) const {
188  return this == &rhs || (_hash == rhs._hash && _dict == rhs._dict);
189 }
190 
191 void MaskDict::add(std::string const & name, int bitId) {
192  _dict[name] = bitId;
193  _hash = boost::hash<MaskPlaneDict>()(_dict);
194 }
195 
197  _dict.erase(name);
198  _hash = boost::hash<MaskPlaneDict>()(_dict);
199 }
200 
202  _dict.clear();
203  _hash = 0x0;
204 }
205 
206 void MaskDict::_addInitialMaskPlanes() {
207  int i = -1;
208  _dict["BAD"] = ++i;
209  _dict["SAT"] = ++i;
210  _dict["INTRP"] = ++i;
211  _dict["CR"] = ++i;
212  _dict["EDGE"] = ++i;
213  _dict["DETECTED"] = ++i;
214  _dict["DETECTED_NEGATIVE"] = ++i;
215  _dict["SUSPECT"] = ++i;
216  _dict["NO_DATA"] = ++i;
217  _hash = boost::hash<MaskPlaneDict>()(_dict);
218 }
219 
220 MaskDict::MaskDict() : _dict(), _hash(0x0) {}
221 
222 MaskDict::MaskDict(MaskPlaneDict const & dict) : _dict(dict), _hash(boost::hash<MaskPlaneDict>()(_dict)) {}
223 
224 }}}} // lsst::afw::image::detail
py::object result
Definition: _schema.cc:430
table::Key< std::string > name
Definition: Amplifier.cc:116
table::Key< int > id
Definition: Detector.cc:162
uint64_t * ptr
Definition: RangeSet.cc:88
table::Key< int > b
table::Key< int > a
std::shared_ptr< MaskDict > detachDefault()
Definition: MaskDict.cc:64
std::shared_ptr< MaskDict > getDefault() const noexcept
Definition: MaskDict.cc:56
std::shared_ptr< MaskDict > copy(MaskDict const &dict)
Definition: MaskDict.cc:49
std::shared_ptr< MaskDict > copyOrGetDefault(MaskPlaneDict const &mpd)
Definition: MaskDict.cc:39
void setDefault(std::shared_ptr< MaskDict > dict)
Definition: MaskDict.cc:60
int getMaskPlane(std::string const &name) const
Definition: MaskDict.cc:176
const_iterator begin() const noexcept
Definition: MaskDict.h:131
static void addAllMasksPlane(std::string const &name, int bitId)
Definition: MaskDict.cc:127
static std::shared_ptr< MaskDict > detachDefault()
Definition: MaskDict.cc:123
const_iterator end() const noexcept
Definition: MaskDict.h:132
static std::shared_ptr< MaskDict > getDefault()
Definition: MaskDict.cc:115
static void setDefault(std::shared_ptr< MaskDict > dict)
Definition: MaskDict.cc:119
bool operator==(MaskDict const &rhs) const
Definition: MaskDict.cc:187
const_iterator find(std::string const &name) const
Definition: MaskDict.h:135
void add(std::string const &name, int bitId)
Definition: MaskDict.cc:191
static std::shared_ptr< MaskDict > copyOrGetDefault(MaskPlaneDict const &dict)
Definition: MaskDict.cc:111
bool empty() const noexcept
Definition: MaskDict.h:141
void erase(std::string const &name)
Definition: MaskDict.cc:196
T clear(T... args)
T empty(T... args)
T endl(T... args)
T erase(T... args)
T find_if(T... args)
T max_element(T... args)
T move(T... args)
Definition: Polygon.cc:25
std::map< std::string, int > MaskPlaneDict
Definition: Mask.h:58
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
class[[deprecated("Removed with no replacement (but see lsst::afw::image::TransmissionCurve). Will be " "removed after v22.")]] FilterProperty final
Describe the properties of a Filter (e.g.
Definition: Filter.h:53
A base class for image defects.
STL namespace.