38 namespace lsst {
namespace meas {
namespace algorithms {
44 struct DoubleGaussianPsfPersistenceHelper {
49 afw::table::Key<double>
b;
51 static DoubleGaussianPsfPersistenceHelper
const &
get() {
52 static DoubleGaussianPsfPersistenceHelper instance;
57 DoubleGaussianPsfPersistenceHelper (
const DoubleGaussianPsfPersistenceHelper&) =
delete;
58 DoubleGaussianPsfPersistenceHelper& operator=(
const DoubleGaussianPsfPersistenceHelper&) =
delete;
61 DoubleGaussianPsfPersistenceHelper (DoubleGaussianPsfPersistenceHelper&&) =
delete;
62 DoubleGaussianPsfPersistenceHelper& operator=(DoubleGaussianPsfPersistenceHelper&&) =
delete;
65 DoubleGaussianPsfPersistenceHelper() :
70 sigma1(schema.addField<
double>(
"sigma1",
"radius of inner Gaussian",
"pixel")),
71 sigma2(schema.addField<
double>(
"sigma2",
"radius of outer Gaussian",
"pixel")),
72 b(schema.addField<
double>(
"b",
"central amplitude of outer Gaussian (inner amplitude == 1)"))
74 schema.getCitizen().markPersistent();
78 class DoubleGaussianPsfFactory :
public afw::table::io::PersistableFactory {
81 virtual PTR(afw::table::io::Persistable)
82 read(InputArchive
const & archive, CatalogVector
const & catalogs)
const {
83 static DoubleGaussianPsfPersistenceHelper
const &
keys = DoubleGaussianPsfPersistenceHelper::get();
86 afw::table::BaseRecord
const & record = catalogs.front().front();
88 return std::make_shared<DoubleGaussianPsf>(
89 record.get(keys.dimensions.getX()),
90 record.get(keys.dimensions.getY()),
91 record.get(keys.sigma1),
92 record.get(keys.sigma2),
97 DoubleGaussianPsfFactory(
std::string const &
name) : afw::table::io::PersistableFactory(name) {}
102 PTR(afw::math::Kernel) makeDoubleGaussianKernel(
103 int width,
int height,
double sigma1,
double &
sigma2,
double b 105 if (b == 0.0 && sigma2 == 0.0) {
108 if (sigma1 <= 0 || sigma2 <= 0) {
110 (boost::format(
"sigma may not be 0: %g, %g") % sigma1 % sigma2).str());
112 afw::math::DoubleGaussianFunction2<double> dg(sigma1, sigma2, b);
113 PTR(afw::math::Kernel) kernel(
new afw::math::AnalyticKernel(width, height, dg));
117 std::string getDoubleGaussianPsfPersistenceName() {
return "DoubleGaussianPsf"; }
119 DoubleGaussianPsfFactory registration(getDoubleGaussianPsfPersistenceName());
124 KernelPsf(makeDoubleGaussianKernel(width, height, sigma1, sigma2, b)),
125 _sigma1(sigma1), _sigma2(sigma2), _b(b)
129 return std::make_shared<DoubleGaussianPsf>(
137 return std::make_shared<DoubleGaussianPsf>(
146 static DoubleGaussianPsfPersistenceHelper
const &
keys = DoubleGaussianPsfPersistenceHelper::get();
149 (*record).set(keys.dimensions.getX(),
getKernel()->getWidth());
150 (*record).set(keys.dimensions.getY(),
getKernel()->getHeight());
153 (*record).set(keys.b,
getB());
154 handle.saveCatalog(catalog);
163 PsfFormatter::doubleGaussianPsfRegistration(
A Psf defined by a Kernel.
virtual void write(OutputArchiveHandle &handle) const
#define LSST_ARCHIVE_ASSERT(EXPR)
static PointKey addFields(Schema &schema, std::string const &name, std::string const &doc, std::string const &unit)
double getSigma1() const
Return the radius of the inner Gaussian.
boost::shared_ptr< afw::math::Kernel const > getKernel() const
Return the Kernel used to define this Psf.
DoubleGaussianPsf(int width, int height, double sigma1, double sigma2=0.0, double b=0.0)
Constructor for a DoubleGaussianPsf.
Represent a Psf as a circularly symmetrical double Gaussian.
virtual std::string getPersistenceName() const
afw::table::PointKey< int > dimensions
virtual boost::shared_ptr< afw::detection::Psf > resized(int width, int height) const
Return a clone with specified kernel dimensions.
double getSigma2() const
Return the radius of the outer Gaussian.
#define LSST_EXCEPT(type,...)
afw::table::Key< double > sigma1
afw::table::Schema schema
double getB() const
Return the ratio of Gaussian peak amplitudes: outer/inner.
io::OutputArchiveHandle OutputArchiveHandle
afw::table::Key< double > b
virtual boost::shared_ptr< afw::detection::Psf > clone() const
Polymorphic deep copy. Usually unnecessary, as Psfs are immutable.
afw::table::Key< double > sigma2
std::shared_ptr< BaseRecord > addNew()