<!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: itkOptLinearInterpolateImageFunction.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>itkOptLinearInterpolateImageFunction.h</h1><a href="itkOptLinearInterpolateImageFunction_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: itkOptLinearInterpolateImageFunction.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-01-10 21:01:46 $</span> <a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.5 $</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 __itkOptLinearInterpolateImageFunction_h</span> <a name="l00018"></a>00018 <span class="preprocessor"></span><span class="preprocessor">#define __itkOptLinearInterpolateImageFunction_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span> <a name="l00020"></a>00020 <span class="preprocessor">#include "<a class="code" href="itkInterpolateImageFunction_8h.html">itkInterpolateImageFunction.h</a>"</span> <a name="l00021"></a>00021 <a name="l00022"></a>00022 <span class="keyword">namespace </span>itk <a name="l00023"></a>00023 { <a name="l00024"></a>00024 <a name="l00042"></a>00042 <span class="keyword">template</span> <<span class="keyword">class</span> TInputImage, <span class="keyword">class</span> TCoordRep = <span class="keywordtype">float</span>> <a name="l00043"></a>00043 <span class="keyword">class </span><a class="code" href="itkWin32Header_8h.html#7d4384322557205d75f0719c012dded6">ITK_EXPORT</a> LinearInterpolateImageFunction : <a name="l00044"></a>00044 <span class="keyword">public</span> InterpolateImageFunction<TInputImage,TCoordRep> <a name="l00045"></a>00045 { <a name="l00046"></a>00046 <span class="keyword">public</span>: <a name="l00048"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#a1ffed024f17ee4e705e5befe1a72698">00048</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1LinearInterpolateImageFunction.html" title="Linearly interpolate an image at specified positions.">LinearInterpolateImageFunction</a> <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a>; <a name="l00049"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#ad5cb7b155ef66378a4e9da383edd5b7">00049</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1InterpolateImageFunction.html">InterpolateImageFunction<TInputImage,TCoordRep></a> <a class="code" href="classitk_1_1InterpolateImageFunction.html">Superclass</a>; <a name="l00050"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#16984144aed0cfcf4682ee43ff32f252">00050</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="l00051"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#4d0c4c9e35a87590e6388257001b6e5b">00051</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="l00052"></a>00052 <a name="l00054"></a>00054 <a class="code" href="itkMacro_8h.html#574c1f86661d409dd32c742307a0e96a">itkTypeMacro</a>(<a class="code" href="classitk_1_1LinearInterpolateImageFunction.html" title="Linearly interpolate an image at specified positions.">LinearInterpolateImageFunction</a>, <a class="code" href="classitk_1_1InterpolateImageFunction.html" title="Base class for all image interpolaters.">InterpolateImageFunction</a>); <a name="l00055"></a>00055 <a name="l00057"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#6f478873c039329091ce8600cd4d0279">00057</a> <a class="code" href="itkMacro_8h.html#0fdc15685afa271857c57e7cd482f30c">itkNewMacro</a>(<a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a>); <a name="l00058"></a>00058 <a name="l00060"></a>00060 <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::OutputType OutputType; <a name="l00061"></a>00061 <a name="l00063"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#49b1e9834e1b34bc62c99f96adf5d1ae">00063</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#3742cbac150bee74daa9aab85ed80623">Superclass::InputImageType</a> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#3742cbac150bee74daa9aab85ed80623">InputImageType</a>; <a name="l00064"></a>00064 <a name="l00066"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#27aed639789aed6096fc65d08ee48c64">00066</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::RealType RealType; <a name="l00067"></a>00067 <a name="l00069"></a>00069 <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="l00070"></a>00070 <a name="l00072"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#9c9fdcbbc9877a69d7074d5256921245">00072</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#be32dc19e71d40322ad9828b9ab40853">Superclass::IndexType</a> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#be32dc19e71d40322ad9828b9ab40853">IndexType</a>; <a name="l00073"></a>00073 <a name="l00075"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#81ba79b4df27cdb2fbe5a91d22d898af">00075</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1ContinuousIndex.html">Superclass::ContinuousIndexType</a> <a class="code" href="classitk_1_1ContinuousIndex.html">ContinuousIndexType</a>; <a name="l00076"></a>00076 <a name="l00085"></a><a class="code" href="classitk_1_1LinearInterpolateImageFunction.html#f113c6c00e7a7af7d30174ecde7b8553">00085</a> <span class="keyword">virtual</span> <span class="keyword">inline</span> OutputType EvaluateAtContinuousIndex( <span class="keyword">const</span> <a name="l00086"></a>00086 <a class="code" href="classitk_1_1ContinuousIndex.html">ContinuousIndexType</a> & <a name="l00087"></a>00087 index )<span class="keyword"> const</span> <a name="l00088"></a>00088 <span class="keyword"> </span>{ <a name="l00089"></a>00089 <span class="keywordflow">return</span> this->EvaluateOptimized( Dispatch< ImageDimension >(), index ); <a name="l00090"></a>00090 } <a name="l00091"></a>00091 <a name="l00092"></a>00092 <span class="keyword">protected</span>: <a name="l00093"></a>00093 <a class="code" href="classitk_1_1LinearInterpolateImageFunction.html" title="Linearly interpolate an image at specified positions.">LinearInterpolateImageFunction</a>(); <a name="l00094"></a>00094 ~<a class="code" href="classitk_1_1LinearInterpolateImageFunction.html" title="Linearly interpolate an image at specified positions.">LinearInterpolateImageFunction</a>(); <a name="l00095"></a>00095 <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="l00096"></a>00096 <a name="l00097"></a>00097 <span class="keyword">private</span>: <a name="l00098"></a>00098 <a class="code" href="classitk_1_1LinearInterpolateImageFunction.html" title="Linearly interpolate an image at specified positions.">LinearInterpolateImageFunction</a>( <span class="keyword">const</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#ac7d0120fe5292cdcdddb6163b8cb150">Self</a>& ); <span class="comment">//purposely not implemented</span> <a name="l00099"></a>00099 <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="l00100"></a>00100 <a name="l00102"></a>00102 <span class="keyword">static</span> <span class="keyword">const</span> <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> m_Neighbors; <a name="l00103"></a>00103 <a name="l00104"></a>00104 <span class="keyword">struct </span>DispatchBase {}; <a name="l00105"></a>00105 <span class="keyword">template</span>< <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> > <span class="keyword">struct </span>Dispatch : DispatchBase {}; <a name="l00106"></a>00106 <a name="l00107"></a>00107 <span class="keyword">inline</span> OutputType EvaluateOptimized( <span class="keyword">const</span> Dispatch<0> &, <a name="l00108"></a>00108 <span class="keyword">const</span> ContinuousIndexType & index)<span class="keyword"> const</span> <a name="l00109"></a>00109 <span class="keyword"> </span>{ <a name="l00110"></a>00110 <span class="keywordflow">return</span> 0; <a name="l00111"></a>00111 } <a name="l00112"></a>00112 <a name="l00113"></a>00113 <span class="keyword">inline</span> OutputType EvaluateOptimized( <span class="keyword">const</span> Dispatch<1>&, <a name="l00114"></a>00114 <span class="keyword">const</span> ContinuousIndexType & index)<span class="keyword"> const</span> <a name="l00115"></a>00115 <span class="keyword"> </span>{ <a name="l00116"></a>00116 <a class="code" href="namespaceHardConnectedComponentImageFilter.html#be32dc19e71d40322ad9828b9ab40853">IndexType</a> basei; <a name="l00117"></a>00117 <a name="l00118"></a>00118 <span class="keywordtype">double</span> i = index[0]; <a name="l00119"></a>00119 basei[0] = (long)i; <a name="l00120"></a>00120 <span class="keywordflow">if</span>( i < 0.0 && <span class="keywordtype">double</span>(basei[0]) != i) <a name="l00121"></a>00121 { <a name="l00122"></a>00122 basei[0]--; <a name="l00123"></a>00123 } <a name="l00124"></a>00124 <a name="l00125"></a>00125 <span class="keywordtype">double</span> distance = i - double(basei[0]); <a name="l00126"></a>00126 <a name="l00127"></a>00127 <span class="keywordtype">double</span> val0 = this->GetInputImage()->GetPixel( basei ); <a name="l00128"></a>00128 <span class="keywordtype">double</span> val1 = val0; <a name="l00129"></a>00129 ++basei[0]; <a name="l00130"></a>00130 val1 = this->GetInputImage()->GetPixel( basei ); <a name="l00131"></a>00131 <a name="l00132"></a>00132 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 + distance * ( val1 - val0 ) ) ); <a name="l00133"></a>00133 } <a name="l00134"></a>00134 <a name="l00135"></a>00135 <span class="keyword">inline</span> OutputType EvaluateOptimized( <span class="keyword">const</span> Dispatch<2>&, <a name="l00136"></a>00136 <span class="keyword">const</span> ContinuousIndexType & index)<span class="keyword"> const</span> <a name="l00137"></a>00137 <span class="keyword"> </span>{ <a name="l00138"></a>00138 <a class="code" href="namespaceHardConnectedComponentImageFilter.html#be32dc19e71d40322ad9828b9ab40853">IndexType</a> basei; <a name="l00139"></a>00139 <a name="l00140"></a>00140 <span class="keywordtype">double</span> i = index[0]; <a name="l00141"></a>00141 basei[0] = (long)i; <a name="l00142"></a>00142 <span class="keywordflow">if</span>( i < 0.0 && <span class="keywordtype">double</span>(basei[0]) != i) <a name="l00143"></a>00143 { <a name="l00144"></a>00144 basei[0]--; <a name="l00145"></a>00145 } <a name="l00146"></a>00146 <span class="keywordtype">double</span> distance0 = i - double(basei[0]); <a name="l00147"></a>00147 <a name="l00148"></a>00148 i = index[1]; <a name="l00149"></a>00149 basei[1] = (long)i; <a name="l00150"></a>00150 <span class="keywordflow">if</span>( i < 0.0 && <span class="keywordtype">double</span>(basei[1]) != i) <a name="l00151"></a>00151 { <a name="l00152"></a>00152 basei[1]--; <a name="l00153"></a>00153 } <a name="l00154"></a>00154 <span class="keywordtype">double</span> distance1 = i - double(basei[1]); <a name="l00155"></a>00155 <a name="l00156"></a>00156 <a name="l00157"></a>00157 <span class="keywordtype">double</span> val00 = this->GetInputImage()->GetPixel( basei ); <a name="l00158"></a>00158 <span class="keywordflow">if</span>(distance0+distance1 == 0) <a name="l00159"></a>00159 { <a name="l00160"></a>00160 <span class="keywordflow">return</span>( static_cast<OutputType>( val00 ) ); <a name="l00161"></a>00161 } <a name="l00162"></a>00162 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(distance1 == 0) <span class="comment">// if they have the same "y"</span> <a name="l00163"></a>00163 { <a name="l00164"></a>00164 ++basei[0]; <span class="comment">// then interpolate across "x"</span> <a name="l00165"></a>00165 <span class="keywordflow">if</span>(basei[0]>this->m_EndIndex[0]) <a name="l00166"></a>00166 { <a name="l00167"></a>00167 <span class="keywordflow">return</span>( static_cast<OutputType>( val00 ) ); <a name="l00168"></a>00168 } <a name="l00169"></a>00169 <span class="keywordtype">double</span> val10 = this->GetInputImage()->GetPixel( basei ); <a name="l00170"></a>00170 <span class="keywordflow">return</span>( static_cast<OutputType>(val00 + distance0 * (val10 - val00)) ); <a name="l00171"></a>00171 } <a name="l00172"></a>00172 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(distance0 == 0) <span class="comment">// if they have the same "x"</span> <a name="l00173"></a>00173 { <a name="l00174"></a>00174 ++basei[1]; <span class="comment">// then interpolate across "y"</span> <a name="l00175"></a>00175 <span class="keywordflow">if</span>(basei[1]>this->m_EndIndex[1]) <a name="l00176"></a>00176 { <a name="l00177"></a>00177 <span class="keywordflow">return</span>( static_cast<OutputType>( val00 ) ); <a name="l00178"></a>00178 } <a name="l00179"></a>00179 <span class="keywordtype">double</span> val01 = this->GetInputImage()->GetPixel( basei ); <a name="l00180"></a>00180 <span class="keywordflow">return</span>( static_cast<OutputType>(val00 + distance1 * (val01 - val00)) ); <a name="l00181"></a>00181 } <a name="l00182"></a>00182 <span class="keywordflow">else</span> <a name="l00183"></a>00183 { <a name="l00184"></a>00184 ++basei[0]; <a name="l00185"></a>00185 <span class="keywordtype">double</span> val10 = this->GetInputImage()->GetPixel( basei ); <a name="l00186"></a>00186 ++basei[1]; <a name="l00187"></a>00187 <span class="keywordtype">double</span> val11 = this->GetInputImage()->GetPixel( basei ); <a name="l00188"></a>00188 --basei[0]; <a name="l00189"></a>00189 <span class="keywordtype">double</span> val01 = this->GetInputImage()->GetPixel( basei ); <a name="l00190"></a>00190 <a name="l00191"></a>00191 <span class="keywordtype">double</span> val0 = val00 + distance1 * ( val01 - val00 ); <span class="comment">// interpolate across "y"</span> <a name="l00192"></a>00192 <span class="keywordtype">double</span> val1 = val10 + distance1 * ( val11 - val10 ); <span class="comment">// interpolate across "y"</span> <a name="l00193"></a>00193 <a name="l00194"></a>00194 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 + distance0 * (val1-val0) ) ); <span class="comment">// interpolate across "X"</span> <a name="l00195"></a>00195 } <a name="l00196"></a>00196 } <a name="l00197"></a>00197 <a name="l00198"></a>00198 <span class="keyword">inline</span> OutputType EvaluateOptimized( <span class="keyword">const</span> Dispatch<3>&, <a name="l00199"></a>00199 <span class="keyword">const</span> ContinuousIndexType & index)<span class="keyword"> const</span> <a name="l00200"></a>00200 <span class="keyword"> </span>{ <a name="l00201"></a>00201 <a class="code" href="namespaceHardConnectedComponentImageFilter.html#be32dc19e71d40322ad9828b9ab40853">IndexType</a> basei; <a name="l00202"></a>00202 <span class="keywordtype">double</span> val[8]; <a name="l00203"></a>00203 <a name="l00204"></a>00204 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> min0; <a name="l00205"></a>00205 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> max0; <a name="l00206"></a>00206 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> min1; <a name="l00207"></a>00207 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> max1; <a name="l00208"></a>00208 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> min2; <a name="l00209"></a>00209 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> max2; <a name="l00210"></a>00210 <a name="l00211"></a>00211 <span class="keywordtype">double</span> i = index[0]; <a name="l00212"></a>00212 basei[0] = (long)i; <a name="l00213"></a>00213 <span class="keywordflow">if</span>( i < 0.0 && <span class="keywordtype">double</span>(basei[0]) != i) <a name="l00214"></a>00214 { <a name="l00215"></a>00215 basei[0]--; <a name="l00216"></a>00216 } <a name="l00217"></a>00217 <span class="keywordtype">double</span> distance0 = i - double(basei[0]); <a name="l00218"></a>00218 <a name="l00219"></a>00219 i = index[1]; <a name="l00220"></a>00220 basei[1] = (long)i; <a name="l00221"></a>00221 <span class="keywordflow">if</span>( i < 0.0 && <span class="keywordtype">double</span>(basei[1]) != i) <a name="l00222"></a>00222 { <a name="l00223"></a>00223 basei[1]--; <a name="l00224"></a>00224 } <a name="l00225"></a>00225 <span class="keywordtype">double</span> distance1 = i - double(basei[1]); <a name="l00226"></a>00226 <a name="l00227"></a>00227 i = index[2]; <a name="l00228"></a>00228 basei[2] = (long)i; <a name="l00229"></a>00229 <span class="keywordflow">if</span>( i < 0.0 && <span class="keywordtype">double</span>(basei[2]) != i) <a name="l00230"></a>00230 { <a name="l00231"></a>00231 basei[2]--; <a name="l00232"></a>00232 } <a name="l00233"></a>00233 <span class="keywordtype">double</span> distance2 = i - double(basei[2]); <a name="l00234"></a>00234 <a name="l00235"></a>00235 val[0] = this->GetInputImage()->GetPixel( basei ); <a name="l00236"></a>00236 <span class="keywordflow">if</span>(distance0+distance1+distance2 == 0) <a name="l00237"></a>00237 { <a name="l00238"></a>00238 <span class="keywordflow">return</span>( static_cast<OutputType>( val[0] ) ); <a name="l00239"></a>00239 } <a name="l00240"></a>00240 <span class="keywordflow">if</span>(distance0 > 0.0) <a name="l00241"></a>00241 { <a name="l00242"></a>00242 min0 = basei[0]; <a name="l00243"></a>00243 max0 = basei[0]+1; <a name="l00244"></a>00244 <span class="keywordflow">if</span>(max0>this->m_EndIndex[0]) <a name="l00245"></a>00245 { <a name="l00246"></a>00246 max0 = this->m_EndIndex[0]; <a name="l00247"></a>00247 } <a name="l00248"></a>00248 } <a name="l00249"></a>00249 <span class="keywordflow">if</span>(distance1 > 0.0) <a name="l00250"></a>00250 { <a name="l00251"></a>00251 min1 = basei[1]; <a name="l00252"></a>00252 max1 = basei[1]+1; <a name="l00253"></a>00253 <span class="keywordflow">if</span>(max1>this->m_EndIndex[1]) <a name="l00254"></a>00254 { <a name="l00255"></a>00255 max1 = this->m_EndIndex[1]; <a name="l00256"></a>00256 } <a name="l00257"></a>00257 } <a name="l00258"></a>00258 <span class="keywordflow">if</span>(distance2 > 0.0) <a name="l00259"></a>00259 { <a name="l00260"></a>00260 min2 = basei[2]; <a name="l00261"></a>00261 max2 = basei[2]+1; <a name="l00262"></a>00262 <span class="keywordflow">if</span>(max2>this->m_EndIndex[2]) <a name="l00263"></a>00263 { <a name="l00264"></a>00264 max2 = this->m_EndIndex[2]; <a name="l00265"></a>00265 } <a name="l00266"></a>00266 } <a name="l00267"></a>00267 <span class="keywordflow">if</span>(distance2 == 0) <a name="l00268"></a>00268 { <a name="l00269"></a>00269 <span class="keywordflow">if</span>(distance1 == 0) <a name="l00270"></a>00270 { <a name="l00271"></a>00271 basei[0] = max0; <a name="l00272"></a>00272 val[1] = this->GetInputImage()->GetPixel( basei ); <a name="l00273"></a>00273 <a name="l00274"></a>00274 <span class="keywordtype">double</span> val0 = val[0] + distance0 * (val[1]-val[0]); <a name="l00275"></a>00275 <a name="l00276"></a>00276 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 ) ); <a name="l00277"></a>00277 } <a name="l00278"></a>00278 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(distance0 == 0) <a name="l00279"></a>00279 { <a name="l00280"></a>00280 basei[1] = max1; <a name="l00281"></a>00281 val[2] = this->GetInputImage()->GetPixel( basei ); <a name="l00282"></a>00282 <a name="l00283"></a>00283 <span class="keywordtype">double</span> val0 = val[0] + distance1 * (val[2]-val[0]); <a name="l00284"></a>00284 <a name="l00285"></a>00285 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 ) ); <a name="l00286"></a>00286 } <a name="l00287"></a>00287 <span class="keywordflow">else</span> <a name="l00288"></a>00288 { <a name="l00289"></a>00289 basei[0] = max0; <a name="l00290"></a>00290 val[1] = this->GetInputImage()->GetPixel( basei ); <a name="l00291"></a>00291 basei[1] = max1; <a name="l00292"></a>00292 val[3] = this->GetInputImage()->GetPixel( basei ); <a name="l00293"></a>00293 basei[0] = min0; <a name="l00294"></a>00294 val[2] = this->GetInputImage()->GetPixel( basei ); <a name="l00295"></a>00295 <a name="l00296"></a>00296 <span class="keywordtype">double</span> val0 = val[0] + distance0 * (val[1]-val[0]); <a name="l00297"></a>00297 <span class="keywordtype">double</span> val1 = val[2] + distance0 * (val[3]-val[2]); <a name="l00298"></a>00298 <a name="l00299"></a>00299 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 + distance1 * (val1 - val0) ) ); <a name="l00300"></a>00300 } <a name="l00301"></a>00301 } <a name="l00302"></a>00302 <span class="keywordflow">else</span> <a name="l00303"></a>00303 { <a name="l00304"></a>00304 basei[2] = max2; <a name="l00305"></a>00305 val[4] = this->GetInputImage()->GetPixel( basei ); <a name="l00306"></a>00306 <span class="keywordflow">if</span>(distance1 == 0) <a name="l00307"></a>00307 { <a name="l00308"></a>00308 <span class="keywordflow">if</span>(distance0 == 0) <a name="l00309"></a>00309 { <a name="l00310"></a>00310 <span class="keywordflow">return</span>( static_cast<OutputType>( val[0] + distance2 <a name="l00311"></a>00311 * (val[4] - val[0]) ) ); <a name="l00312"></a>00312 } <a name="l00313"></a>00313 <span class="keywordflow">else</span> <a name="l00314"></a>00314 { <a name="l00315"></a>00315 basei[0] = max0; <a name="l00316"></a>00316 val[5] = this->GetInputImage()->GetPixel( basei ); <a name="l00317"></a>00317 basei[2] = min2; <a name="l00318"></a>00318 val[1] = this->GetInputImage()->GetPixel( basei ); <a name="l00319"></a>00319 <a name="l00320"></a>00320 <span class="keywordtype">double</span> val0 = val[0] + distance0 * (val[1]-val[0]); <a name="l00321"></a>00321 <span class="keywordtype">double</span> val1 = val[4] + distance0 * (val[5]-val[4]); <a name="l00322"></a>00322 <a name="l00323"></a>00323 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 + distance2 * (val1 - val0) ) ); <a name="l00324"></a>00324 } <a name="l00325"></a>00325 } <a name="l00326"></a>00326 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(distance0 == 0) <a name="l00327"></a>00327 { <a name="l00328"></a>00328 basei[1] = max1; <a name="l00329"></a>00329 val[6] = this->GetInputImage()->GetPixel( basei ); <a name="l00330"></a>00330 basei[2] = min2; <a name="l00331"></a>00331 val[2] = this->GetInputImage()->GetPixel( basei ); <a name="l00332"></a>00332 <a name="l00333"></a>00333 <span class="keywordtype">double</span> val0 = val[0] + distance1 * (val[2]-val[0]); <a name="l00334"></a>00334 <span class="keywordtype">double</span> val1 = val[4] + distance1 * (val[6]-val[4]); <a name="l00335"></a>00335 <a name="l00336"></a>00336 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 + distance2 * (val1 - val0) ) ); <a name="l00337"></a>00337 } <a name="l00338"></a>00338 <span class="keywordflow">else</span> <a name="l00339"></a>00339 { <a name="l00340"></a>00340 basei[0] = max0; <a name="l00341"></a>00341 val[5] = this->GetInputImage()->GetPixel( basei ); <a name="l00342"></a>00342 basei[1] = max1; <a name="l00343"></a>00343 val[7] = this->GetInputImage()->GetPixel( basei ); <a name="l00344"></a>00344 basei[0] = min0; <a name="l00345"></a>00345 val[6] = this->GetInputImage()->GetPixel( basei ); <a name="l00346"></a>00346 basei[2] = min2; <a name="l00347"></a>00347 val[2] = this->GetInputImage()->GetPixel( basei ); <a name="l00348"></a>00348 basei[0] = max0; <a name="l00349"></a>00349 val[3] = this->GetInputImage()->GetPixel( basei ); <a name="l00350"></a>00350 basei[1] = min1; <a name="l00351"></a>00351 val[1] = this->GetInputImage()->GetPixel( basei ); <a name="l00352"></a>00352 <a name="l00353"></a>00353 <span class="keywordtype">double</span> val00 = val[0] + distance0 * (val[1]-val[0]); <a name="l00354"></a>00354 <span class="keywordtype">double</span> val01 = val[2] + distance0 * (val[3]-val[2]); <a name="l00355"></a>00355 <span class="keywordtype">double</span> val0 = val00 + distance1 * (val01-val00); <a name="l00356"></a>00356 <a name="l00357"></a>00357 <span class="keywordtype">double</span> val10 = val[4] + distance0 * (val[5]-val[4]); <a name="l00358"></a>00358 <span class="keywordtype">double</span> val11 = val[6] + distance0 * (val[7]-val[6]); <a name="l00359"></a>00359 <span class="keywordtype">double</span> val1 = val10 + distance1 * (val11-val10); <a name="l00360"></a>00360 <a name="l00361"></a>00361 <span class="keywordflow">return</span>( static_cast<OutputType>( val0 + distance2 * (val1-val0) ) ); <a name="l00362"></a>00362 } <a name="l00363"></a>00363 } <a name="l00364"></a>00364 } <a name="l00365"></a>00365 <a name="l00366"></a>00366 <span class="keyword">inline</span> OutputType EvaluateOptimized( <span class="keyword">const</span> DispatchBase &, <a name="l00367"></a>00367 <span class="keyword">const</span> ContinuousIndexType & index)<span class="keyword"> const</span> <a name="l00368"></a>00368 <span class="keyword"> </span>{ <a name="l00369"></a>00369 <span class="keywordflow">return</span> this->EvaluateUnoptimized( index ); <a name="l00370"></a>00370 } <a name="l00371"></a>00371 <a name="l00372"></a>00372 <span class="keyword">virtual</span> <span class="keyword">inline</span> OutputType EvaluateUnoptimized( <a name="l00373"></a>00373 <span class="keyword">const</span> ContinuousIndexType & index) <span class="keyword">const</span>; <a name="l00374"></a>00374 }; <a name="l00375"></a>00375 <a name="l00376"></a>00376 } <span class="comment">// end namespace itk</span> <a name="l00377"></a>00377 <a name="l00378"></a>00378 <span class="comment">// Define instantiation macro for this template.</span> <a name="l00379"></a><a class="code" href="itkOptLinearInterpolateImageFunction_8h.html#26303dd6e093c190cf3eaac0ac343236">00379</a> <span class="preprocessor">#define ITK_TEMPLATE_LinearInterpolateImageFunction(_, EXPORT, x, y) namespace itk { \</span> <a name="l00380"></a>00380 <span class="preprocessor"> _(2(class EXPORT LinearInterpolateImageFunction< ITK_TEMPLATE_2 x >)) \</span> <a name="l00381"></a>00381 <span class="preprocessor"> namespace Templates { typedef LinearInterpolateImageFunction< ITK_TEMPLATE_2 x > \</span> <a name="l00382"></a>00382 <span class="preprocessor"> LinearInterpolateImageFunction##y; } \</span> <a name="l00383"></a>00383 <span class="preprocessor"> }</span> <a name="l00384"></a>00384 <span class="preprocessor"></span> <a name="l00385"></a>00385 <span class="preprocessor">#if ITK_TEMPLATE_EXPLICIT</span> <a name="l00386"></a>00386 <span class="preprocessor"></span><span class="preprocessor"># include "Templates/itkLinearInterpolateImageFunction+-.h"</span> <a name="l00387"></a>00387 <span class="preprocessor">#endif</span> <a name="l00388"></a>00388 <span class="preprocessor"></span> <a name="l00389"></a>00389 <span class="preprocessor">#if ITK_TEMPLATE_TXX</span> <a name="l00390"></a>00390 <span class="preprocessor"></span><span class="preprocessor"># include "itkOptLinearInterpolateImageFunction.txx"</span> <a name="l00391"></a>00391 <span class="preprocessor">#endif</span> <a name="l00392"></a>00392 <span class="preprocessor"></span> <a name="l00393"></a>00393 <span class="preprocessor">#endif</span> <a name="l00394"></a>00394 <span class="preprocessor"></span> </pre></div></div> <hr><address><small> Generated at Thu May 7 23:44:58 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>