<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> <html><head><meta name="robots" content="noindex"> <meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1"> <title>ITK: itkFastSymmetricForcesDemonsRegistrationFunction.h Source File</title> <link href="DoxygenStyle.css" rel="stylesheet" type="text/css"> </head><body bgcolor="#ffffff"> <!-- Section customized for INSIGHT : Tue Jul 17 01:02:45 2001 --> <center> <a href="index.html" class="qindex">Main Page</a> <a href="modules.html" class="qindex">Groups</a> <a href="namespaces.html" class="qindex">Namespace List</a> <a href="hierarchy.html" class="qindex">Class Hierarchy</a> <a href="classes.html" class="qindex">Alphabetical List</a> <a href="annotated.html" class="qindex">Compound List</a> <a href="files.html" class="qindex">File List</a> <a href="namespacemembers.html" class="qindex">Namespace Members</a> <a href="functions.html" class="qindex">Compound Members</a> <a href="globals.html" class="qindex">File Members</a> <a href="pages.html" class="qindex">Concepts</a></center> <!-- Generated by Doxygen 1.5.9 --> <div class="navpath"><a class="el" href="dir_02e23d841c198393a6e961e9a0678fff.html">Code</a> » <a class="el" href="dir_16b0d2d22c345042f35cebae66b36010.html">Algorithms</a> </div> <div class="contents"> <h1>itkFastSymmetricForcesDemonsRegistrationFunction.h</h1><a href="itkFastSymmetricForcesDemonsRegistrationFunction_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*=========================================================================</span> <a name="l00002"></a>00002 <span class="comment"></span> <a name="l00003"></a>00003 <span class="comment"> Program: Insight Segmentation & Registration Toolkit</span> <a name="l00004"></a>00004 <span class="comment"> Module: $RCSfile: itkFastSymmetricForcesDemonsRegistrationFunction.h,v $</span> <a name="l00005"></a>00005 <span class="comment"> Language: C++</span> <a name="l00006"></a>00006 <span class="comment"> Date: $Date: 2008-12-21 19:13:11 $</span> <a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.3 $</span> <a name="l00008"></a>00008 <span class="comment"></span> <a name="l00009"></a>00009 <span class="comment"> Copyright (c) Insight Software Consortium. All rights reserved.</span> <a name="l00010"></a>00010 <span class="comment"> See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.</span> <a name="l00011"></a>00011 <span class="comment"></span> <a name="l00012"></a>00012 <span class="comment"> This software is distributed WITHOUT ANY WARRANTY; without even </span> <a name="l00013"></a>00013 <span class="comment"> the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR </span> <a name="l00014"></a>00014 <span class="comment"> PURPOSE. See the above copyright notices for more information.</span> <a name="l00015"></a>00015 <span class="comment"></span> <a name="l00016"></a>00016 <span class="comment">=========================================================================*/</span> <a name="l00017"></a>00017 <span class="preprocessor">#ifndef __itkFastSymmetricForcesDemonsRegistrationFunction_h</span> <a name="l00018"></a>00018 <span class="preprocessor"></span><span class="preprocessor">#define __itkFastSymmetricForcesDemonsRegistrationFunction_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span> <a name="l00020"></a>00020 <span class="preprocessor">#include "<a class="code" href="itkPDEDeformableRegistrationFunction_8h.html">itkPDEDeformableRegistrationFunction.h</a>"</span> <a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="itkPoint_8h.html">itkPoint.h</a>"</span> <a name="l00022"></a>00022 <span class="preprocessor">#include "<a class="code" href="itkCovariantVector_8h.html">itkCovariantVector.h</a>"</span> <a name="l00023"></a>00023 <span class="preprocessor">#include "<a class="code" href="itkInterpolateImageFunction_8h.html">itkInterpolateImageFunction.h</a>"</span> <a name="l00024"></a>00024 <span class="preprocessor">#include "<a class="code" href="itkLinearInterpolateImageFunction_8h.html">itkLinearInterpolateImageFunction.h</a>"</span> <a name="l00025"></a>00025 <span class="preprocessor">#include "<a class="code" href="itkCentralDifferenceImageFunction_8h.html">itkCentralDifferenceImageFunction.h</a>"</span> <a name="l00026"></a>00026 <span class="preprocessor">#include "<a class="code" href="itkWarpImageFilter_8h.html">itkWarpImageFilter.h</a>"</span> <a name="l00027"></a>00027 <a name="l00028"></a>00028 <span class="keyword">namespace </span>itk { <a name="l00029"></a>00029 <a name="l00045"></a>00045 <span class="keyword">template</span><<span class="keyword">class</span> TFixedImage, <span class="keyword">class</span> TMovingImage, <span class="keyword">class</span> TDeformationField> <a name="l00046"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">00046</a> <span class="keyword">class </span><a class="code" href="itkWin32Header_8h.html#7d4384322557205d75f0719c012dded6">ITK_EXPORT</a> <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">FastSymmetricForcesDemonsRegistrationFunction</a> : <a name="l00047"></a>00047 <span class="keyword">public</span> <a class="code" href="classitk_1_1PDEDeformableRegistrationFunction.html">PDEDeformableRegistrationFunction</a>< TFixedImage, <a name="l00048"></a>00048 TMovingImage, TDeformationField> <a name="l00049"></a>00049 { <a name="l00050"></a>00050 <span class="keyword">public</span>: <a name="l00051"></a>00051 <a name="l00053"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#54718e54a461861b9d9871d5860afefb">00053</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">FastSymmetricForcesDemonsRegistrationFunction</a> <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">Self</a>; <a name="l00054"></a>00054 <span class="keyword">typedef</span> <a class="code" href="classitk_1_1PDEDeformableRegistrationFunction.html">PDEDeformableRegistrationFunction</a>< TFixedImage, <a name="l00055"></a>00055 TMovingImage, TDeformationField > <a name="l00056"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#b3d6bc8f28bdd165601e70968e48b1f8">00056</a> <a class="code" href="classitk_1_1PDEDeformableRegistrationFunction.html">Superclass</a>; <a name="l00057"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#66648ab23d36cfd9a38175e5e8e41a9d">00057</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer<Self></a> <a class="code" href="classitk_1_1SmartPointer.html">Pointer</a>; <a name="l00058"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#036ad73da38ef7ac6009d0359cbc5ebc">00058</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer<const Self></a> <a class="code" href="classitk_1_1SmartPointer.html">ConstPointer</a>; <a name="l00059"></a>00059 <a name="l00061"></a>00061 <a class="code" href="itkMacro_8h.html#0fdc15685afa271857c57e7cd482f30c">itkNewMacro</a>(<a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">Self</a>); <a name="l00062"></a>00062 <a name="l00064"></a>00064 <a class="code" href="itkMacro_8h.html#574c1f86661d409dd32c742307a0e96a">itkTypeMacro</a>( <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">FastSymmetricForcesDemonsRegistrationFunction</a>, <a name="l00065"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#59ceefcc74abe2d1da5c332b51a631a8">00065</a> <a class="code" href="classitk_1_1PDEDeformableRegistrationFunction.html">PDEDeformableRegistrationFunction</a> ); <a name="l00066"></a>00066 <a name="l00068"></a>00068 <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::MovingImageType MovingImageType; <a name="l00069"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#1b1c7c0040f41bcbd219f5fa7e032b31">00069</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::MovingImagePointer MovingImagePointer; <a name="l00070"></a>00070 <a name="l00072"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#1ada585ff044fbcdebc9977d0dc9f8ad">00072</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::FixedImageType FixedImageType; <a name="l00073"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#a21884b56211b11e4f063439e2a0dd3b">00073</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::FixedImagePointer FixedImagePointer; <a name="l00074"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#92d64bd8d7015435b0bfbc8057c71233">00074</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#be32dc19e71d40322ad9828b9ab40853">FixedImageType::IndexType</a> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#be32dc19e71d40322ad9828b9ab40853">IndexType</a>; <a name="l00075"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#6584f3ea4d5a2f7d00ff82ae9b1bf046">00075</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#5bf82236e1e50e653de84c8ead448599">FixedImageType::SizeType</a> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#5bf82236e1e50e653de84c8ead448599">SizeType</a>; <a name="l00076"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#1a256e6a6f01c7e8f29eb0bbb0a60ffd">00076</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> FixedImageType::SpacingType SpacingType; <a name="l00077"></a>00077 <a name="l00079"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#d6f4f85665b031b7e5876a85fbf3f5e4">00079</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::DeformationFieldType DeformationFieldType; <a name="l00080"></a>00080 <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::DeformationFieldTypePointer <a name="l00081"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#aead0b4a004ad10fe2d60a9cd78dd2c4">00081</a> DeformationFieldTypePointer; <a name="l00082"></a>00082 <a name="l00084"></a>00084 <a class="code" href="itkMacro_8h.html#1905cf6b2e88235996fa7ff73b2fb0eb">itkStaticConstMacro</a>(ImageDimension, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span>,Superclass::ImageDimension); <a name="l00085"></a>00085 <a name="l00087"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#31290b41739fcfef3d2fa6fcc6a5fd79">00087</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::PixelType PixelType; <a name="l00088"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#81d8472dd84a5929106924223b2d3288">00088</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::RadiusType RadiusType; <a name="l00089"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#0a4ac39db3ae34ca79387cd6d0f6d6ab">00089</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1ConstNeighborhoodIterator.html" title="Const version of NeighborhoodIterator, defining iteration of a local N-dimensional...">Superclass::NeighborhoodType</a> <a class="code" href="classitk_1_1ConstNeighborhoodIterator.html" title="Const version of NeighborhoodIterator, defining iteration of a local N-dimensional...">NeighborhoodType</a>; <a name="l00090"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#6efea3e73863a22e12ced2b455ec845f">00090</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1Vector.html" title="A templated class holding a n-Dimensional vector.">Superclass::FloatOffsetType</a> <a class="code" href="classitk_1_1Vector.html" title="A templated class holding a n-Dimensional vector.">FloatOffsetType</a>; <a name="l00091"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#9b9ae1948470079ce0b7e49fbeedf02c">00091</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::TimeStepType TimeStepType; <a name="l00092"></a>00092 <a name="l00094"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#87c127775c0dcc92b1414f57dbfbf29d">00094</a> <span class="keyword">typedef</span> <span class="keywordtype">double</span> CoordRepType; <a name="l00095"></a>00095 <span class="keyword">typedef</span> <a class="code" href="classitk_1_1InterpolateImageFunction.html" title="Base class for all image interpolaters.">InterpolateImageFunction<MovingImageType,CoordRepType></a> <a name="l00096"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#b10de693fa5308155df7d804c780b9b7">00096</a> <a class="code" href="classitk_1_1InterpolateImageFunction.html" title="Base class for all image interpolaters.">InterpolatorType</a>; <a name="l00097"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#16839e5b246eb5b8ef6236800525dd81">00097</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1SmartPointer.html">InterpolatorType::Pointer</a> <a class="code" href="classitk_1_1SmartPointer.html">InterpolatorPointer</a>; <a name="l00098"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#b675e09f8665aa0dd4a2c66fc3ad3122">00098</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1Point.html" title="A templated class holding a geometric point in n-Dimensional space.">InterpolatorType::PointType</a> <a class="code" href="classitk_1_1Point.html" title="A templated class holding a geometric point in n-Dimensional space.">PointType</a>; <a name="l00099"></a>00099 <span class="keyword">typedef</span> <a class="code" href="classitk_1_1LinearInterpolateImageFunction.html" title="Linearly interpolate an image at specified positions.">LinearInterpolateImageFunction<MovingImageType,CoordRepType></a> <a name="l00100"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#f068b06f4a2bac1f2d517e0ba8812e82">00100</a> <a class="code" href="classitk_1_1LinearInterpolateImageFunction.html" title="Linearly interpolate an image at specified positions.">DefaultInterpolatorType</a>; <a name="l00101"></a>00101 <a name="l00103"></a>00103 <span class="keyword">typedef</span> <a class="code" href="classitk_1_1WarpImageFilter.html" title="Warps an image using an input deformation field.">WarpImageFilter<MovingImageType,MovingImageType,DeformationFieldType></a> <a name="l00104"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#e10c6f785253e462d540d4223a42d247">00104</a> <a class="code" href="classitk_1_1WarpImageFilter.html" title="Warps an image using an input deformation field.">WarperType</a>; <a name="l00105"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#d4dfe4a011b00b58c0a6caf057462f3a">00105</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1SmartPointer.html">WarperType::Pointer</a> <a class="code" href="classitk_1_1SmartPointer.html">WarperPointer</a>; <a name="l00106"></a>00106 <a name="l00108"></a>00108 <span class="keyword">typedef</span> <a class="code" href="classitk_1_1CovariantVector.html" title="A templated class holding a n-Dimensional covariant vector.">CovariantVector<double,itkGetStaticConstMacro(ImageDimension)></a> <a name="l00109"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#cd2b916a0df0e4d7abefef2d63263eb3">00109</a> <a class="code" href="classitk_1_1CovariantVector.html" title="A templated class holding a n-Dimensional covariant vector.">CovariantVectorType</a>; <a name="l00110"></a>00110 <a name="l00112"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#9706709c0245bac8713dcbd3c46aadb2">00112</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1CentralDifferenceImageFunction.html" title="Calculate the derivative by central differencing.">CentralDifferenceImageFunction<FixedImageType></a> <a class="code" href="classitk_1_1CentralDifferenceImageFunction.html" title="Calculate the derivative by central differencing.">GradientCalculatorType</a>; <a name="l00113"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#e25cf3f1a282e4e4d1438d926b280d21">00113</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1SmartPointer.html">GradientCalculatorType::Pointer</a> <a class="code" href="classitk_1_1SmartPointer.html">GradientCalculatorPointer</a>; <a name="l00114"></a>00114 <a name="l00115"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#80763ea859edef8491baad0803c1d5e9">00115</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1CentralDifferenceImageFunction.html" title="Calculate the derivative by central differencing.">CentralDifferenceImageFunction<MovingImageType></a> <a class="code" href="classitk_1_1CentralDifferenceImageFunction.html" title="Calculate the derivative by central differencing.">MovingGradientCalculatorType</a>; <a name="l00116"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#7960de140568d2bbba7c05a98d9c2ee0">00116</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1SmartPointer.html">MovingGradientCalculatorType::Pointer</a> <a class="code" href="classitk_1_1SmartPointer.html">MovingGradientCalculatorPointer</a>; <a name="l00117"></a>00117 <a name="l00119"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#1bef5207e58ac74c63cb4c3facb1eabe">00119</a> <span class="keywordtype">void</span> SetMovingImageInterpolator( <a class="code" href="classitk_1_1InterpolateImageFunction.html" title="Base class for all image interpolaters.">InterpolatorType</a> * ptr ) <a name="l00120"></a>00120 { <a name="l00121"></a>00121 m_MovingImageInterpolator = ptr; m_MovingImageWarper->SetInterpolator( ptr ); <a name="l00122"></a>00122 } <a name="l00123"></a>00123 <a name="l00125"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#fefb798dd7c0c6239b3c323a4e927eac">00125</a> <a class="code" href="classitk_1_1InterpolateImageFunction.html" title="Base class for all image interpolaters.">InterpolatorType</a> * GetMovingImageInterpolator(<span class="keywordtype">void</span>) <a name="l00126"></a>00126 { <a name="l00127"></a>00127 <span class="keywordflow">return</span> m_MovingImageInterpolator; <a name="l00128"></a>00128 } <a name="l00129"></a>00129 <a name="l00131"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#538010948a672824af4eac5a4b0437ab">00131</a> <span class="keyword">virtual</span> TimeStepType ComputeGlobalTimeStep(<span class="keywordtype">void</span> * <a class="code" href="itkMacro_8h.html#6be3b439282da7a34278c0594589b972">itkNotUsed</a>(GlobalData))<span class="keyword"> const</span> <a name="l00132"></a>00132 <span class="keyword"> </span>{ <a name="l00133"></a>00133 <span class="keywordflow">return</span> m_TimeStep; <a name="l00134"></a>00134 } <a name="l00135"></a>00135 <a name="l00138"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#9b049176c9bdae7128a858c61b80bef3">00138</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> *GetGlobalDataPointer()<span class="keyword"> const</span> <a name="l00139"></a>00139 <span class="keyword"> </span>{ <a name="l00140"></a>00140 <a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html">GlobalDataStruct</a> *global = <span class="keyword">new</span> <a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html">GlobalDataStruct</a>(); <a name="l00141"></a>00141 global-><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html#2628b3da8ac5fc24e152872dfc8a2653">m_SumOfSquaredDifference</a> = 0.0; <a name="l00142"></a>00142 global-><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html#382063026df146a1f98c41d8407f7a8c">m_NumberOfPixelsProcessed</a> = 0L; <a name="l00143"></a>00143 global-><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html#a5ff3f451592dd11beeb4cfb83325335">m_SumOfSquaredChange</a> = 0; <a name="l00144"></a>00144 <span class="keywordflow">return</span> global; <a name="l00145"></a>00145 } <a name="l00147"></a>00147 <a name="l00149"></a>00149 <span class="keyword">virtual</span> <span class="keywordtype">void</span> ReleaseGlobalDataPointer( <span class="keywordtype">void</span> *GlobalData ) <span class="keyword">const</span>; <a name="l00150"></a>00150 <a name="l00152"></a>00152 <span class="keyword">virtual</span> <span class="keywordtype">void</span> InitializeIteration(); <a name="l00153"></a>00153 <a name="l00156"></a>00156 <span class="keyword">virtual</span> PixelType ComputeUpdate(<span class="keyword">const</span> NeighborhoodType &neighborhood, <a name="l00157"></a>00157 <span class="keywordtype">void</span> *globalData, <a name="l00158"></a>00158 <span class="keyword">const</span> FloatOffsetType &offset = FloatOffsetType(0.0)); <a name="l00159"></a>00159 <a name="l00163"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#f7a6986caa68b422c5d84a71fc9b6f70">00163</a> <span class="keyword">virtual</span> <span class="keywordtype">double</span> GetMetric()<span class="keyword"> const</span> <a name="l00164"></a>00164 <span class="keyword"> </span>{ <a name="l00165"></a>00165 <span class="keywordflow">return</span> m_Metric; <a name="l00166"></a>00166 } <a name="l00167"></a>00167 <a name="l00169"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#2d862e782e0a3fc88b23bcf74bf4a0e0">00169</a> <span class="keyword">virtual</span> <span class="keyword">const</span> <span class="keywordtype">double</span> &GetRMSChange()<span class="keyword"> const</span> <a name="l00170"></a>00170 <span class="keyword"> </span>{ <a name="l00171"></a>00171 <span class="keywordflow">return</span> m_RMSChange; <a name="l00172"></a>00172 } <a name="l00173"></a>00173 <a name="l00178"></a>00178 <span class="keyword">virtual</span> <span class="keywordtype">void</span> SetIntensityDifferenceThreshold(<span class="keywordtype">double</span>); <a name="l00179"></a>00179 <span class="keyword">virtual</span> <span class="keywordtype">double</span> GetIntensityDifferenceThreshold() <span class="keyword">const</span>; <a name="l00181"></a>00181 <a name="l00182"></a>00182 <span class="keyword">protected</span>: <a name="l00183"></a>00183 <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">FastSymmetricForcesDemonsRegistrationFunction</a>(); <a name="l00184"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#80f500654a3291a5bff7bb06dc3d4fbf">00184</a> <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#80f500654a3291a5bff7bb06dc3d4fbf">~FastSymmetricForcesDemonsRegistrationFunction</a>() {} <a name="l00185"></a>00185 <span class="keywordtype">void</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#920ae03048a5ad715d803ca2e92b52de">PrintSelf</a>(std::ostream& os, <a class="code" href="classitk_1_1Indent.html" title="Control indentation during Print() invocation.">Indent</a> indent) <span class="keyword">const</span>; <a name="l00186"></a>00186 <a name="l00188"></a><a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html#a7f9f05d1631363a5f800ea1f8abcc6f">00188</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1ConstNeighborhoodIterator.html">ConstNeighborhoodIterator<FixedImageType></a> <a class="code" href="classitk_1_1ConstNeighborhoodIterator.html">FixedImageNeighborhoodIteratorType</a>; <a name="l00189"></a>00189 <a name="l00192"></a><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html">00192</a> <span class="keyword">struct </span><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html">GlobalDataStruct</a> <a name="l00193"></a>00193 { <a name="l00194"></a><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html#2628b3da8ac5fc24e152872dfc8a2653">00194</a> <span class="keywordtype">double</span> m_SumOfSquaredDifference; <a name="l00195"></a><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html#382063026df146a1f98c41d8407f7a8c">00195</a> <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> m_NumberOfPixelsProcessed; <a name="l00196"></a><a class="code" href="structitk_1_1FastSymmetricForcesDemonsRegistrationFunction_1_1GlobalDataStruct.html#a5ff3f451592dd11beeb4cfb83325335">00196</a> <span class="keywordtype">double</span> m_SumOfSquaredChange; <a name="l00197"></a>00197 }; <a name="l00198"></a>00198 <a name="l00199"></a>00199 <span class="keyword">private</span>: <a name="l00200"></a>00200 <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">FastSymmetricForcesDemonsRegistrationFunction</a>(<span class="keyword">const</span> <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">Self</a>&); <span class="comment">//purposely not implemented</span> <a name="l00201"></a>00201 <span class="keywordtype">void</span> operator=(<span class="keyword">const</span> <a class="code" href="classitk_1_1FastSymmetricForcesDemonsRegistrationFunction.html">Self</a>&); <span class="comment">//purposely not implemented</span> <a name="l00202"></a>00202 <a name="l00204"></a>00204 <span class="keywordtype">double</span> m_Normalizer; <a name="l00205"></a>00205 <a name="l00207"></a>00207 <a class="code" href="classitk_1_1SmartPointer.html">GradientCalculatorPointer</a> m_FixedImageGradientCalculator; <a name="l00208"></a>00208 <a name="l00210"></a>00210 <a class="code" href="classitk_1_1SmartPointer.html">MovingGradientCalculatorPointer</a> m_WarpedMovingImageGradientCalculator; <a name="l00211"></a>00211 <a name="l00213"></a>00213 <a class="code" href="classitk_1_1SmartPointer.html">InterpolatorPointer</a> m_MovingImageInterpolator; <a name="l00214"></a>00214 <a name="l00216"></a>00216 <a class="code" href="classitk_1_1SmartPointer.html">WarperPointer</a> m_MovingImageWarper; <a name="l00217"></a>00217 <a name="l00219"></a>00219 TimeStepType m_TimeStep; <a name="l00220"></a>00220 <a name="l00222"></a>00222 <span class="keywordtype">double</span> m_DenominatorThreshold; <a name="l00223"></a>00223 <a name="l00225"></a>00225 <span class="keywordtype">double</span> m_IntensityDifferenceThreshold; <a name="l00226"></a>00226 <a name="l00230"></a>00230 <span class="keyword">mutable</span> <span class="keywordtype">double</span> m_Metric; <a name="l00231"></a>00231 <span class="keyword">mutable</span> <span class="keywordtype">double</span> m_SumOfSquaredDifference; <a name="l00232"></a>00232 <span class="keyword">mutable</span> <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> m_NumberOfPixelsProcessed; <a name="l00233"></a>00233 <span class="keyword">mutable</span> <span class="keywordtype">double</span> m_RMSChange; <a name="l00234"></a>00234 <span class="keyword">mutable</span> <span class="keywordtype">double</span> m_SumOfSquaredChange; <a name="l00235"></a>00235 <a name="l00237"></a>00237 <span class="keyword">mutable</span> <a class="code" href="classitk_1_1SimpleFastMutexLock.html" title="Critical section locking class that can be allocated on the stack.">SimpleFastMutexLock</a> m_MetricCalculationLock; <a name="l00238"></a>00238 <a name="l00239"></a>00239 }; <a name="l00240"></a>00240 <a name="l00241"></a>00241 <a name="l00242"></a>00242 } <span class="comment">// end namespace itk</span> <a name="l00243"></a>00243 <a name="l00244"></a>00244 <span class="preprocessor">#ifndef ITK_MANUAL_INSTANTIATION</span> <a name="l00245"></a>00245 <span class="preprocessor"></span><span class="preprocessor">#include "itkFastSymmetricForcesDemonsRegistrationFunction.txx"</span> <a name="l00246"></a>00246 <span class="preprocessor">#endif</span> <a name="l00247"></a>00247 <span class="preprocessor"></span> <a name="l00248"></a>00248 <span class="preprocessor">#endif</span> <a name="l00249"></a>00249 <span class="preprocessor"></span> </pre></div></div> <hr><address><small> Generated at Thu May 7 22:40:23 2009 for ITK by <a href="http://www.stack.nl/~dimitri/doxygen/index.html"> <img src="http://www.stack.nl/~dimitri/doxygen/doxygen.png" alt="doxygen" align="middle" border=0 width=110 height=53> </a> 1.5.9 written by <a href="mailto:dimitri@stack.nl">Dimitri van Heesch</a>, © 1997-2000</small></address> </body> </html>