Authors belong to:
- University of LYON http://www.universite-lyon.fr/
- - Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
+ - Léon Bérard cancer center http://www.centreleonberard.fr
- CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
This software is distributed WITHOUT ANY WARRANTY; without even
- BSD See included LICENSE.txt file
- CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
-======================================================================-====*/
+===========================================================================**/
/*=========================================================================
// gets integrated into the main directories.
#include "itkConfigure.h"
-#ifdef ITK_USE_OPTIMIZED_REGISTRATION_METHODS
#include "itkOptMeanSquaresImageToImageMetricFor3DBLUTFFD.txx"
-#else
#include "itkMeanSquaresImageToImageMetricFor3DBLUTFFD.h"
#include "itkImageRegionConstIteratorWithIndex.h"
const unsigned int ParametersDimension = this->GetNumberOfParameters();
derivative = DerivativeType( ParametersDimension );
- derivative.Fill( NumericTraits<ITK_TYPENAME DerivativeType::ValueType>::Zero );
+ derivative.Fill( NumericTraits<typename DerivativeType::ValueType>::Zero );
ti.GoToBegin();
if( this->m_Interpolator->IsInsideBuffer( transformedPoint ) ) {
const RealType movingValue = this->m_Interpolator->Evaluate( transformedPoint );
- const TransformJacobianType & jacobian =
- this->m_Transform->GetJacobian( inputPoint );
-
+ TransformJacobianType jacobian;
+ this->m_Transform->ComputeJacobianWithRespectToParameters( inputPoint, jacobian );
const RealType fixedValue = ti.Value();
this->m_NumberOfPixelsCounted++;
const unsigned int ParametersDimension = this->GetNumberOfParameters();
derivative = DerivativeType( ParametersDimension );
- derivative.Fill( NumericTraits<ITK_TYPENAME DerivativeType::ValueType>::Zero );
+ derivative.Fill( NumericTraits<typename DerivativeType::ValueType>::Zero );
ti.GoToBegin();
if( this->m_Interpolator->IsInsideBuffer( transformedPoint ) ) {
const RealType movingValue = this->m_Interpolator->Evaluate( transformedPoint );
- const TransformJacobianType & jacobian =
- this->m_Transform->GetJacobian( inputPoint );
-
+ TransformJacobianType jacobian;
+ this->m_Transform->ComputeJacobianWithRespectToParameters( inputPoint, jacobian );
const RealType fixedValue = ti.Value();
this->m_NumberOfPixelsCounted++;