itkReducedVarianceModelBuilder.h

/*
 * This file is part of the statismo library.
 *
 * Author: Marcel Luethi (marcel.luethi@unibas.ch)
 *
 * Copyright (c) 2011 University of Basel
 * All rights reserved.
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions
 * are met:
 *
 * Redistributions of source code must retain the above copyright notice,
 * this list of conditions and the following disclaimer.
 *
 * Redistributions in binary form must reproduce the above copyright
 * notice, this list of conditions and the following disclaimer in the
 * documentation and/or other materials provided with the distribution.
 *
 * Neither the name of the project's author nor the names of its
 * contributors may be used to endorse or promote products derived from
 * this software without specific prior written permission.
 *
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
 * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
 * HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
 * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED
 * TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 *
 */


#ifndef ITK_PARTIALLY_FIXED_MODELBUILDER_H_
#define ITK_PARTIALLY_FIXED_MODELBUILDER_H_

#include <itkObject.h>

#include "itkDataManager.h"
#include "itkStatisticalModel.h"

#include "ReducedVarianceModelBuilder.h"
#include "statismoITKConfig.h"
#include "StatismoUtils.h"

namespace itk {

template <class Representer>
class ReducedVarianceModelBuilder : public Object {
  public:

    typedef ReducedVarianceModelBuilder            Self;
    typedef Object  Superclass;
    typedef SmartPointer<Self>                Pointer;
    typedef SmartPointer<const Self>          ConstPointer;

    itkNewMacro( Self );
    itkTypeMacro( ReducedVarianceModelBuilder, Object );

    typedef statismo::ReducedVarianceModelBuilder<Representer> ImplType;


    template <class F>
    typename boost::result_of<F()>::type callstatismoImpl(F f) const {
        try {
            return f();
        } catch (statismo::StatisticalModelException& s) {
            itkExceptionMacro(<< s.what());
        }
    }


    ReducedVarianceModelBuilder() : m_impl(ImplType::Create()) {}

    virtual ~ReducedVarianceModelBuilder() {
        if (m_impl) {
            delete m_impl;
            m_impl = 0;
        }
    }



    typename StatisticalModel<Representer>::Pointer BuildNewModelWithLeadingComponents(const StatisticalModel<Representer>* model, unsigned numberOfPrincipalComponents) {
        statismo::StatisticalModel<Representer>* model_statismo = model->GetstatismoImplObj();
        statismo::StatisticalModel<Representer>* new_model_statismo = callstatismoImpl(boost::bind(&ImplType::BuildNewModelWithLeadingComponents, this->m_impl, model_statismo, numberOfPrincipalComponents));
        typename StatisticalModel<Representer>::Pointer model_itk = StatisticalModel<Representer>::New();
        model_itk->SetstatismoImplObj(new_model_statismo);
        return model_itk;
    }

    typename StatisticalModel<Representer>::Pointer BuildNewModelWithVariance(const StatisticalModel<Representer>* model, double totalVariance) {
        statismo::StatisticalModel<Representer>* model_statismo = model->GetstatismoImplObj();
        statismo::StatisticalModel<Representer>* new_model_statismo = callstatismoImpl(boost::bind(&ImplType::BuildNewModelWithVariance, this->m_impl, model_statismo, totalVariance));
        typename StatisticalModel<Representer>::Pointer model_itk = StatisticalModel<Representer>::New();
        model_itk->SetstatismoImplObj(new_model_statismo);
        return model_itk;
    }

    is_deprecated typename StatisticalModel<Representer>::Pointer BuildNewModelFromModel(const StatisticalModel<Representer>* model, double totalVariance) {
        statismo::StatisticalModel<Representer>* model_statismo = model->GetstatismoImplObj();
        statismo::StatisticalModel<Representer>* new_model_statismo = callstatismoImpl(boost::bind(&ImplType::BuildNewModelFromModel, this->m_impl, model_statismo, totalVariance));
        typename StatisticalModel<Representer>::Pointer model_itk = StatisticalModel<Representer>::New();
        model_itk->SetstatismoImplObj(new_model_statismo);
        return model_itk;
    }


  private:
    ReducedVarianceModelBuilder(const ReducedVarianceModelBuilder& orig);
    ReducedVarianceModelBuilder& operator=(const ReducedVarianceModelBuilder& rhs);

    ImplType* m_impl;
};


}

#endif /* ITK_PARTIALLY_FIXED_MODEL_BUILDER */