nonrigid.hpp

Go to the documentation of this file.

// cpd - Coherent Point Drift
// Copyright (C) 2017 Pete Gadomski <pete.gadomski@gmail.com>
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License along
// with this program; if not, write to the Free Software Foundation, Inc.,
// 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.


#pragma once

#include <cpd/transform.hpp>

namespace cpd {

const double DEFAULT_BETA = 3.0;
const double DEFAULT_LAMBDA = 3.0;

struct NonrigidResult : public Result {};

class Nonrigid : public Transform<NonrigidResult> {
public:
    Nonrigid()
      : Transform()
      , m_lambda(DEFAULT_LAMBDA)
      , m_beta(DEFAULT_BETA)
      , m_linked(DEFAULT_LINKED) {}

    void init(const Matrix& fixed, const Matrix& moving);

    void modify_probabilities(Probabilities& probabilities) const;

    Nonrigid& beta(double beta) {
        m_beta = beta;
        return *this;
    }

    Nonrigid& lambda(double lambda) {
        m_lambda = lambda;
        return *this;
    }

    NonrigidResult compute_one(const Matrix& fixed, const Matrix& moving,
                               const Probabilities& probabilities,
                               double sigma2) const;

    Nonrigid& linked(bool linked) {
        m_linked = linked;
        return *this;
    }

    virtual bool linked() const { return m_linked; }

private:
    Matrix m_g;
    Matrix m_w;
    double m_lambda;
    double m_beta;
    bool m_linked;
};

NonrigidResult nonrigid(const Matrix& fixed, const Matrix& moving);
} // namespace cpd