<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.Section1
{page:Section1;}
-->
</style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="DE" link="blue" vlink="purple">
<div class="Section1">
<p class="MsoNormal"><span lang="EN-US">Hello all,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I am trying to compute the first and second order derivatives of a 3D image using some smoothing, so I am using
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black;
background:blue;mso-highlight:blue">itk::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032;background:
silver;mso-highlight:silver">RecursiveGaussianImageFilter</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">.
When I use images of, say, 200x200x200 voxels it works fine, but when I try with, for example 255x255x255 or bigger, it crashes at run time. I follow a simmililar scheme as in the ITK user guide for computing second order derivatives, and apparently the execution
crashes while doing one of the duplicator->Update(); lines (not necessarily the first occurrence of that)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">As the image size gets bigger, it crashes earlier (i.e. with 255x255x255 it crashes while computing Ixz, and with 260x260x260 it crashes at Iyz which
is the previous computation).<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">Anyone know what I may be doing wrong? It looks like memory allocation issues to me, but as I use smart pointers I don’t see where it can be.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">Thank you very much,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">Alberto<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">CODE:<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#7F0055">typedef</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> itk::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032;background:silver;mso-highlight:silver">RecursiveGaussianImageFilter</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"><</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Image3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
> </span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">ShortFilterType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">;</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#7F0055">typedef</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> itk::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032;background:silver;mso-highlight:silver">RecursiveGaussianImageFilter</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"><</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">,</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
> </span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">FilterType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">;</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#7F0055">typedef</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> itk::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">ImageDuplicator</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black"><
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> >
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">DuplicatorType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">;</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#7F0055">typedef</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> itk::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">ImageDuplicator</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black"><
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Image3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> >
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">ShortDuplicatorType</span><span lang="EN-US" style="font-size:
10.0pt;font-family:"Courier New";color:black">;</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DuplicatorType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
duplicator = </span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">DuplicatorType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">::<i>New</i>();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//************** 1st order derivatives ***************************</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">ShortFilterType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
dx = </span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">ShortFilterType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">::<i>New</i>();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->SetDirection( 0 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//<u>dx</u>
works along x</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#7F0055">if</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> (sigmaG != 0)</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> dx->SetSigma( sigmaG );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">// Specify the derivative order in each direction</span><span lang="EN-US" style="font-size:
10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->SetFirstOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">// build the processing pipe</span><span lang="EN-US" style="font-size:10.0pt;font-family:
"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->SetInput( itkImage );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//update and save the output</span><span lang="EN-US" style="font-size:10.0pt;font-family:
"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//*** Ix</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->SetInputImage( dx->GetOutput() );</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Ix = duplicator->GetOutput(); </span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:#3F7F5F">//first order derivative along z</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Iy</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->SetDirection( 1 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//<u>dx</u>
works along y</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Iy = duplicator->GetOutput();</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Iz</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->SetDirection( 2 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//<u>dx</u>
works along z</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">dx->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Iz = duplicator->GetOutput();</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//************** 2nd order derivatives **************</span><span lang="EN-US" style="font-size:
10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//All derivatives are computed
<u>indepentdently</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">ShortFilterType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
ga = </span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">ShortFilterType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">::<i>New</i>();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">FilterType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
gb = </span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">FilterType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">::<i>New</i>();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">FilterType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
gc = </span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#005032">FilterType</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:black">::<i>New</i>();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//DuplicatorType::Pointer duplicator = DuplicatorType::New();</span><span lang="EN-US" style="font-size:
10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">ga->SetDirection( 0 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//<u>ga</u>
works along x</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gb->SetDirection( 1 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//<u>gb</u>
works along y</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->SetDirection( 2 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//<u>gc</u>
works along z</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#7F0055">if</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> (sigma != 0){</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"> ga->SetSigma( sigma );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">gb->SetSigma( sigma );</span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;
font-family:"Courier New";color:black"> gc->SetSigma( sigma );</span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">}</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">// Specify the derivative order in each direction</span><span lang="EN-US" style="font-size:
10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">ga->SetZeroOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//No derivation
in x (only smoothing)</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gb->SetZeroOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">// no derivation
in y (only smoothing)</span><span lang="EN-US" style="font-size:
10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->SetSecondOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">// 2nd order derivation
in z (plus smoothing)</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">// build the processing pipe</span><span lang="EN-US" style="font-size:10.0pt;font-family:
"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">ga->SetInput( itkImage );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gb->SetInput( ga->GetOutput() );</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->SetInput( gb->GetOutput() );</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->SetInputImage( gc->GetOutput() );</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//update and save the output</span><span lang="EN-US" style="font-size:10.0pt;font-family:
"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Izz</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Izz = duplicator->GetOutput(); </span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";color:#3F7F5F">//second order derivative along z</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Iyy</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->SetDirection( 1 );
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//
<u>gc</u> now works along y</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gb->SetDirection( 2 );
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//
<u>gb</u> now works along z</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Iyy = duplicator->GetOutput();</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Ixx</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->SetDirection( 0 );
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//
<u>gc</u> now works along X</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">ga->SetDirection( 1 );
</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
color:#3F7F5F">//
<u>ga</u> now works along Y</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Ixx = duplicator->GetOutput();</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">// Now we compute cross
<u>derivates</u></span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Iyz</u>=<u>Izy</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">ga->SetDirection( 0 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gb->SetDirection( 1 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->SetDirection( 2 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;
font-family:"Courier New";color:black">ga->SetZeroOrder();</span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;
font-family:"Courier New";color:black">gb->SetFirstOrder();</span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;
font-family:"Courier New";color:black">gc->SetFirstOrder();</span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Iyz = duplicator->GetOutput();</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Ixz</u>=<u>Izx</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">ga->SetDirection( 1 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gb->SetDirection( 0 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//<u>gc</u>->SetDirection( 2 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//<u>ga</u>->SetZeroOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">/</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">/<u>gb</u>->SetFirstOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//<u>gc</u>->SetFirstOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Ixz = duplicator->GetOutput();</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//***
<u>Ixy</u> = <u>Iyx</u></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">ga->SetDirection( 2 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//<u>gb</u>->SetDirection( 0 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->SetDirection( 1 );</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//<u>ga</u>->SetZeroOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//<u>gb</u>->SetFirstOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#3F7F5F">//<u>gc</u>->SetFirstOrder();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">gc->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">duplicator->Update();</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">DoubleImage3DType</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">::</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:#005032">Pointer</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">
Ixy = duplicator->GetOutput();</span><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New""><o:p></o:p></span></p>
</div>
<br>
<hr>
<font face="Arial" color="Gray" size="1">The information contained in this message may be confidential and legally protected under applicable law. The message is intended solely for the addressee(s). If you are not the intended recipient, you are hereby notified
that any use, forwarding, dissemination, or reproduction of this message is strictly prohibited and may be unlawful. If you are not the intended recipient, please contact the sender by return e-mail and destroy all copies of the original message.<br>
</font>
</body>
</html>