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
KernelPersistenceHelper.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 
25 #include "lsst/pex/exceptions.h"
26 #include "lsst/afw/math/Kernel.h"
28 
29 namespace lsst {
30 namespace afw {
31 namespace math {
32 
34  : schema(),
35  dimensions(afw::table::PointKey<int>::addFields(schema, "dimensions",
36  "dimensions of a Kernel's images", "pixel")),
37  center(afw::table::PointKey<int>::addFields(schema, "center", "center point in a Kernel image",
38  "pixel")) {
39  if (nSpatialFunctions > 0) {
41  "spatialfunctions", "archive IDs for the Kernel's spatial functions", nSpatialFunctions);
42  }
43 }
44 
46  : schema(schema_), dimensions(schema["dimensions"]), center(schema["center"]) {
47  try {
48  spatialFunctions = schema["spatialfunctions"];
49  } catch (...) {
50  }
51 }
52 
54  afw::table::io::OutputArchiveHandle& handle, Kernel const& kernel) const {
55  afw::table::BaseCatalog catalog = handle.makeCatalog(schema);
57  record->set(dimensions, lsst::geom::Point2I(kernel.getDimensions()));
58  record->set(center, kernel.getCtr());
59  if (spatialFunctions.isValid()) {
60  writeSpatialFunctions(handle, *record, kernel._spatialFunctionList);
61  }
62  handle.saveCatalog(catalog);
63  return record;
64 }
65 
68  std::vector<std::shared_ptr<Kernel::SpatialFunction>> const& spatialFunctionList) const {
69  ndarray::Array<int, 1, 1> array = record[spatialFunctions];
70  for (std::size_t n = 0; n < spatialFunctionList.size(); ++n) {
71  array[n] = handle.put(spatialFunctionList[n]);
72  }
73 }
74 
76  afw::table::io::InputArchive const& archive, afw::table::BaseRecord const& record) const {
77  ndarray::Array<int const, 1, 1> array = record[spatialFunctions];
78  std::vector<std::shared_ptr<Kernel::SpatialFunction>> spatialFunctionList(array.getSize<0>());
79  for (std::size_t n = 0; n < spatialFunctionList.size(); ++n) {
80  spatialFunctionList[n] = archive.get<SpatialFunction>(array[n]);
81  LSST_ARCHIVE_ASSERT(array[n] == 0 || (spatialFunctionList[n]));
82  }
83  return spatialFunctionList;
84 }
85 } // namespace math
86 } // namespace afw
87 } // namespace lsst
afw::table::PointKey< int > dimensions
Definition: GaussianPsf.cc:49
#define LSST_ARCHIVE_ASSERT(EXPR)
An assertion macro used to validate the structure of an InputArchive.
Definition: Persistable.h:48
table::Schema schema
Definition: python.h:134
A Function taking two arguments.
Definition: Function.h:259
Kernels are used for convolution with MaskedImages and (eventually) Images.
Definition: Kernel.h:111
lsst::geom::Extent2I const getDimensions() const
Return the Kernel's dimensions (width, height)
Definition: Kernel.h:213
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:450
lsst::geom::Point2I getCtr() const
Return index of kernel's center.
Definition: Kernel.h:235
Tag types used to declare specialized field types.
Definition: misc.h:32
Base class for all records.
Definition: BaseRecord.h:31
std::shared_ptr< RecordT > addNew()
Create a new record, add it to the end of the catalog, and return a pointer to it.
Definition: Catalog.h:485
Defines the fields and offsets for a table.
Definition: Schema.h:50
Key< T > addField(Field< T > const &field, bool doReplace=false)
Add a new field to the Schema, and return the associated Key.
Definition: Schema.cc:668
A multi-catalog archive object used to load table::io::Persistable objects.
Definition: InputArchive.h:31
std::shared_ptr< Persistable > get(int id) const
Load the Persistable with the given ID and return it.
An object passed to Persistable::write to allow it to persist itself.
void saveCatalog(BaseCatalog const &catalog)
Save a catalog in the archive.
BaseCatalog makeCatalog(Schema const &schema)
Return a new, empty catalog with the given schema.
int put(Persistable const *obj, bool permissive=false)
Save an object to the archive and return a unique ID that can be used to retrieve it from an InputArc...
A base class for image defects.
T size(T... args)
afw::table::Key< afw::table::Array< int > > spatialFunctions
void writeSpatialFunctions(afw::table::io::OutputArchiveHandle &handle, afw::table::BaseRecord &record, std::vector< SpatialFunctionPtr > const &spatialFunctionList) const
std::shared_ptr< afw::table::BaseRecord > write(afw::table::io::OutputArchiveHandle &handle, Kernel const &kernel) const
std::vector< SpatialFunctionPtr > readSpatialFunctions(afw::table::io::InputArchive const &archive, afw::table::BaseRecord const &record) const