<!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: itkCoxDeBoorBSplineKernelFunction.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_dd352b1cbb04f9ff09fa5f33693fed9f.html">Review</a> </div> <div class="contents"> <h1>itkCoxDeBoorBSplineKernelFunction.h</h1><a href="itkCoxDeBoorBSplineKernelFunction_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: itkCoxDeBoorBSplineKernelFunction.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-08-04 12:54:50 $</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 __itkCoxDeBoorBSplineKernelFunction_h</span> <a name="l00018"></a>00018 <span class="preprocessor"></span><span class="preprocessor">#define __itkCoxDeBoorBSplineKernelFunction_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span> <a name="l00020"></a>00020 <span class="preprocessor">#include "<a class="code" href="itkKernelFunction_8h.html">itkKernelFunction.h</a>"</span> <a name="l00021"></a>00021 <span class="preprocessor">#include "itkNumericTraits.h"</span> <a name="l00022"></a>00022 <span class="preprocessor">#include "vnl/vnl_math.h"</span> <a name="l00023"></a>00023 <span class="preprocessor">#include "vnl/vnl_real_polynomial.h"</span> <a name="l00024"></a>00024 <span class="preprocessor">#include "vnl/vnl_matrix.h"</span> <a name="l00025"></a>00025 <a name="l00026"></a>00026 <span class="keyword">namespace </span>itk <a name="l00027"></a>00027 { <a name="l00028"></a>00028 <a name="l00056"></a>00056 <span class="keyword">template</span> <<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> VSplineOrder = 3> <a name="l00057"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html">00057</a> <span class="keyword">class </span><a class="code" href="itkWin32Header_8h.html#7d4384322557205d75f0719c012dded6">ITK_EXPORT</a> <a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html" title="BSpline kernel used for density estimation and nonparameteric regression.">CoxDeBoorBSplineKernelFunction</a> <a name="l00058"></a>00058 : <span class="keyword">public</span> <a class="code" href="classitk_1_1KernelFunction.html" title="Kernel used for density estimation and nonparameteric regression.">KernelFunction</a> <a name="l00059"></a>00059 { <a name="l00060"></a>00060 <span class="keyword">public</span>: <a name="l00062"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#e415195412dea6f8a2b8d2e2c1e8e320">00062</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html" title="BSpline kernel used for density estimation and nonparameteric regression.">CoxDeBoorBSplineKernelFunction</a> <a class="code" href="classitk_1_1KernelFunction.html" title="Kernel used for density estimation and nonparameteric regression.">Self</a>; <a name="l00063"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#40e222f663f663fed3f90594d6de160e">00063</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1KernelFunction.html" title="Kernel used for density estimation and nonparameteric regression.">KernelFunction</a> <a class="code" href="classitk_1_1FunctionBase.html">Superclass</a>; <a name="l00064"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#0adb1dea76506da9126b7d7208b70891">00064</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="l00065"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#9726202814fbc9ea46e2b7ab2286cd71">00065</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="l00066"></a>00066 <a name="l00068"></a>00068 <a class="code" href="itkMacro_8h.html#0fdc15685afa271857c57e7cd482f30c">itkNewMacro</a>(<a class="code" href="classitk_1_1KernelFunction.html" title="Kernel used for density estimation and nonparameteric regression.">Self</a>); <a name="l00069"></a>00069 <a name="l00071"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#39cc5b56f4cb0abb7118255f5ce4f821">00071</a> <a class="code" href="itkMacro_8h.html#574c1f86661d409dd32c742307a0e96a">itkTypeMacro</a>( <a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html" title="BSpline kernel used for density estimation and nonparameteric regression.">CoxDeBoorBSplineKernelFunction</a>, <a class="code" href="classitk_1_1KernelFunction.html" title="Kernel used for density estimation and nonparameteric regression.">KernelFunction</a> ); <a name="l00072"></a>00072 <a name="l00073"></a>00073 <span class="keyword">typedef</span> <span class="keywordtype">double</span> RealType; <a name="l00074"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#a1a792fe60184ec0ec17ac6100fb6ba8">00074</a> <span class="keyword">typedef</span> <a class="code" href="classvnl__vector.html">vnl_vector<RealType></a> <a class="code" href="classvnl__vector.html">VectorType</a>; <a name="l00075"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#cecce7e2931967bd3ac2fadbaf089f56">00075</a> <span class="keyword">typedef</span> vnl_real_polynomial PolynomialType; <a name="l00076"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#14ac980a5887a1faf68a51fb2123f3bf">00076</a> <span class="keyword">typedef</span> <a class="code" href="classvnl__matrix.html">vnl_matrix<RealType></a> <a class="code" href="classvnl__matrix.html">MatrixType</a>; <a name="l00077"></a>00077 <a name="l00079"></a>00079 <span class="keywordtype">void</span> SetSplineOrder( <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ); <a name="l00080"></a>00080 <a class="code" href="itkMacro_8h.html#1600305bb01e39e87c472d8bf29bd1cf">itkGetMacro</a>( SplineOrder, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ); <a name="l00082"></a>00082 <a name="l00084"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#f64be7da59a48c6a3fdcce1e84666bc2">00084</a> <span class="keyword">inline</span> RealType Evaluate( <span class="keyword">const</span> RealType & u )<span class="keyword"> const</span> <a name="l00085"></a>00085 <span class="keyword"> </span>{ <a name="l00086"></a>00086 RealType absValue = vnl_math_abs( u ); <a name="l00087"></a>00087 <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> which; <a name="l00088"></a>00088 <span class="keywordflow">if</span>( this->m_SplineOrder % 2 == 0 ) <a name="l00089"></a>00089 { <a name="l00090"></a>00090 which = <span class="keyword">static_cast<</span><span class="keywordtype">unsigned</span> <span class="keywordtype">int</span><span class="keyword">></span>( absValue+0.5 ); <a name="l00091"></a>00091 } <a name="l00092"></a>00092 <span class="keywordflow">else</span> <a name="l00093"></a>00093 { <a name="l00094"></a>00094 which = <span class="keyword">static_cast<</span><span class="keywordtype">unsigned</span> <span class="keywordtype">int</span><span class="keyword">></span>( absValue ); <a name="l00095"></a>00095 } <a name="l00096"></a>00096 <span class="keywordflow">if</span>( which < this->m_BSplineShapeFunctions.rows() ) <a name="l00097"></a>00097 { <a name="l00098"></a>00098 <span class="keywordflow">return</span> PolynomialType( <a name="l00099"></a>00099 this->m_BSplineShapeFunctions.get_row( which ) ).evaluate( absValue ); <a name="l00100"></a>00100 } <a name="l00101"></a>00101 <span class="keywordflow">else</span> <a name="l00102"></a>00102 { <a name="l00103"></a>00103 <span class="keywordflow">return</span> NumericTraits<RealType>::Zero; <a name="l00104"></a>00104 } <a name="l00105"></a>00105 } <a name="l00107"></a>00107 <a name="l00109"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#e67cc271cf18cd7fce48c80f930a23af">00109</a> <span class="keyword">inline</span> RealType EvaluateDerivative( <span class="keyword">const</span> <span class="keywordtype">double</span> & u )<span class="keyword"> const</span> <a name="l00110"></a>00110 <span class="keyword"> </span>{ <a name="l00111"></a>00111 <span class="keywordflow">return</span> this->EvaluateNthDerivative( u, 1 ); <a name="l00112"></a>00112 } <a name="l00113"></a>00113 <a name="l00115"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#182f8aca59b4599eb1649d641e5c414f">00115</a> <span class="keyword">inline</span> RealType EvaluateNthDerivative( <span class="keyword">const</span> <span class="keywordtype">double</span> & u, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> n )<span class="keyword"> const</span> <a name="l00116"></a>00116 <span class="keyword"> </span>{ <a name="l00117"></a>00117 RealType absValue = vnl_math_abs( u ); <a name="l00118"></a>00118 <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> which; <a name="l00119"></a>00119 <span class="keywordflow">if</span>( this->m_SplineOrder % 2 == 0 ) <a name="l00120"></a>00120 { <a name="l00121"></a>00121 which = <span class="keyword">static_cast<</span><span class="keywordtype">unsigned</span> <span class="keywordtype">int</span><span class="keyword">></span>( absValue+0.5 ); <a name="l00122"></a>00122 } <a name="l00123"></a>00123 <span class="keywordflow">else</span> <a name="l00124"></a>00124 { <a name="l00125"></a>00125 which = <span class="keyword">static_cast<</span><span class="keywordtype">unsigned</span> <span class="keywordtype">int</span><span class="keyword">></span>( absValue ); <a name="l00126"></a>00126 } <a name="l00127"></a>00127 <span class="keywordflow">if</span>( which < this->m_BSplineShapeFunctions.rows() ) <a name="l00128"></a>00128 { <a name="l00129"></a>00129 PolynomialType polynomial( <a name="l00130"></a>00130 this->m_BSplineShapeFunctions.get_row( which ) ); <a name="l00131"></a>00131 <span class="keywordflow">for</span>( <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i = 0; i < n; i++ ) <a name="l00132"></a>00132 { <a name="l00133"></a>00133 polynomial = polynomial.derivative(); <a name="l00134"></a>00134 } <a name="l00135"></a>00135 RealType der = polynomial.evaluate( absValue ); <a name="l00136"></a>00136 <span class="keywordflow">if</span>( u < NumericTraits<RealType>::Zero && n % 2 != 0 ) <a name="l00137"></a>00137 { <a name="l00138"></a>00138 <span class="keywordflow">return</span> -der; <a name="l00139"></a>00139 } <a name="l00140"></a>00140 <span class="keywordflow">else</span> <a name="l00141"></a>00141 { <a name="l00142"></a>00142 <span class="keywordflow">return</span> der; <a name="l00143"></a>00143 } <a name="l00144"></a>00144 } <a name="l00145"></a>00145 <span class="keywordflow">else</span> <a name="l00146"></a>00146 { <a name="l00147"></a>00147 <span class="keywordflow">return</span> NumericTraits<RealType>::Zero; <a name="l00148"></a>00148 } <a name="l00149"></a>00149 } <a name="l00151"></a>00151 <a name="l00157"></a><a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html#83e1fcbcede0713138e7def9f4cc1ec2">00157</a> <a class="code" href="classvnl__matrix.html">MatrixType</a> GetShapeFunctions() <a name="l00158"></a>00158 { <a name="l00159"></a>00159 <span class="keywordflow">return</span> this->m_BSplineShapeFunctions; <a name="l00160"></a>00160 } <a name="l00161"></a>00161 <a name="l00166"></a>00166 MatrixType GetShapeFunctionsInZeroToOneInterval(); <a name="l00167"></a>00167 <a name="l00168"></a>00168 <span class="keyword">protected</span>: <a name="l00169"></a>00169 <a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html" title="BSpline kernel used for density estimation and nonparameteric regression.">CoxDeBoorBSplineKernelFunction</a>(); <a name="l00170"></a>00170 ~<a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html" title="BSpline kernel used for density estimation and nonparameteric regression.">CoxDeBoorBSplineKernelFunction</a>(); <a name="l00171"></a>00171 <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="l00172"></a>00172 <a name="l00173"></a>00173 <span class="keyword">private</span>: <a name="l00174"></a>00174 <a class="code" href="classitk_1_1CoxDeBoorBSplineKernelFunction.html" title="BSpline kernel used for density estimation and nonparameteric regression.">CoxDeBoorBSplineKernelFunction</a>( <span class="keyword">const</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#ac7d0120fe5292cdcdddb6163b8cb150">Self</a>& ); <span class="comment">//purposely not implemented</span> <a name="l00175"></a>00175 <span class="keywordtype">void</span> operator=( <span class="keyword">const</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#ac7d0120fe5292cdcdddb6163b8cb150">Self</a>& ); <span class="comment">//purposely not implemented</span> <a name="l00176"></a>00176 <a name="l00181"></a>00181 <span class="keywordtype">void</span> GenerateBSplineShapeFunctions( <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ); <a name="l00182"></a>00182 <a name="l00189"></a>00189 PolynomialType CoxDeBoor( <span class="keywordtype">unsigned</span> <span class="keywordtype">short</span>, <a name="l00190"></a>00190 VectorType, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span>, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ); <a name="l00191"></a>00191 <a name="l00192"></a>00192 MatrixType m_BSplineShapeFunctions; <a name="l00193"></a>00193 <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> m_SplineOrder; <a name="l00194"></a>00194 <a name="l00195"></a>00195 }; <a name="l00196"></a>00196 <a name="l00197"></a>00197 } <span class="comment">// end namespace itk</span> <a name="l00198"></a>00198 <a name="l00199"></a>00199 <span class="preprocessor">#ifndef ITK_MANUAL_INSTANTIATION</span> <a name="l00200"></a>00200 <span class="preprocessor"></span><span class="preprocessor">#include "itkCoxDeBoorBSplineKernelFunction.txx"</span> <a name="l00201"></a>00201 <span class="preprocessor">#endif</span> <a name="l00202"></a>00202 <span class="preprocessor"></span> <a name="l00203"></a>00203 <span class="preprocessor">#endif</span> <a name="l00204"></a>00204 <span class="preprocessor"></span> </pre></div></div> <hr><address><small> Generated at Thu May 7 22:30:21 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>