54struct DoubleGaussianPsfPersistenceHelper {
59 afw::table::Key<double>
b;
61 static DoubleGaussianPsfPersistenceHelper
const& get() {
62 static DoubleGaussianPsfPersistenceHelper instance;
67 DoubleGaussianPsfPersistenceHelper(
const DoubleGaussianPsfPersistenceHelper&) =
delete;
68 DoubleGaussianPsfPersistenceHelper& operator=(
const DoubleGaussianPsfPersistenceHelper&) =
delete;
71 DoubleGaussianPsfPersistenceHelper(DoubleGaussianPsfPersistenceHelper&&) =
delete;
72 DoubleGaussianPsfPersistenceHelper& operator=(DoubleGaussianPsfPersistenceHelper&&) =
delete;
75 DoubleGaussianPsfPersistenceHelper()
77 dimensions(afw::table::PointKey<
int>::addFields(
schema,
"dimensions",
"width/height of kernel",
79 sigma1(
schema.addField<double>(
"sigma1",
"radius of inner Gaussian",
"pixel")),
80 sigma2(
schema.addField<double>(
"sigma2",
"radius of outer Gaussian",
"pixel")),
81 b(
schema.addField<double>(
"b",
"central amplitude of outer Gaussian (inner amplitude == 1)")) {}
84class DoubleGaussianPsfFactory :
public afw::table::io::PersistableFactory {
87 read(InputArchive
const& archive, CatalogVector
const& catalogs)
const {
88 static DoubleGaussianPsfPersistenceHelper
const&
keys = DoubleGaussianPsfPersistenceHelper::get();
91 afw::table::BaseRecord
const& record = catalogs.front().front();
93 return std::make_shared<DoubleGaussianPsf>(
94 record.get(
keys.dimensions.getX()), record.get(
keys.dimensions.getY()),
95 record.get(
keys.sigma1), record.get(
keys.sigma2), record.get(
keys.b));
98 DoubleGaussianPsfFactory(
std::string const& name) : afw::table::io::PersistableFactory(
name) {}
104makeDoubleGaussianKernel(
int width,
int height,
double sigma1,
double&
sigma2,
double b) {
105 if (
b == 0.0 &&
sigma2 == 0.0) {
110 (boost::format(
"sigma may not be 0: %g, %g") %
sigma1 %
sigma2).
str());
112 afw::math::DoubleGaussianFunction2<double> dg(
sigma1,
sigma2,
b);
117std::string getDoubleGaussianPsfPersistenceName() {
return "DoubleGaussianPsf"; }
119DoubleGaussianPsfFactory registration(getDoubleGaussianPsfPersistenceName());
123DoubleGaussianPsf::DoubleGaussianPsf(
int width,
int height,
double sigma1,
double sigma2,
double b)
130 return std::make_shared<DoubleGaussianPsf>(
getKernel()->getWidth(),
getKernel()->getHeight(), _sigma1,
135 return std::make_shared<DoubleGaussianPsf>(width, height, _sigma1, _sigma2, _b);
141 static DoubleGaussianPsfPersistenceHelper
const&
keys = DoubleGaussianPsfPersistenceHelper::get();
144 (*record).set(
keys.dimensions.getX(),
getKernel()->getWidth());
145 (*record).set(
keys.dimensions.getY(),
getKernel()->getHeight());
149 handle.saveCatalog(catalog);
table::Key< std::string > name
afw::table::PointKey< int > dimensions
afw::table::Schema schema
afw::table::Key< double > sigma2
afw::table::Key< double > sigma1
afw::table::Key< double > b
#define LSST_EXCEPT(type,...)
#define LSST_ARCHIVE_ASSERT(EXPR)
virtual std::shared_ptr< afw::table::io::Persistable > read(afw::table::io::InputArchive const &archive, afw::table::io::CatalogVector const &catalogs) const
std::shared_ptr< RecordT > addNew()
static std::shared_ptr< T > dynamicCast(std::shared_ptr< Persistable > const &ptr)
io::OutputArchiveHandle OutputArchiveHandle
std::string getPersistenceName() const override
void write(OutputArchiveHandle &handle) const override
double getSigma1() const
Return the radius of the inner Gaussian.
double getSigma2() const
Return the radius of the outer Gaussian.
std::shared_ptr< afw::detection::Psf > clone() const override
Polymorphic deep copy. Usually unnecessary, as Psfs are immutable.
double getB() const
Return the ratio of Gaussian peak amplitudes: outer/inner.
std::shared_ptr< afw::detection::Psf > resized(int width, int height) const override
Return a clone with specified kernel dimensions.
A Psf defined by a Kernel.
std::shared_ptr< afw::math::Kernel const > getKernel() const
Return the Kernel used to define this Psf.