<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:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@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 WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
-->
</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=NL link=blue vlink=purple>
<div class=WordSection1>
<p class=MsoNormal><span lang=EN-GB>Dear all,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-GB>I am trying to “translate” the example
<i>DeformableRegistration4.cxx</i> (deformable registration of two 2D images)
into Python. However, I keep getting this error:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:red'>itk::ERROR: Region
ImageRegion (0027F104)</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:red'> Dimension:
2</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:red'> Index: [-1,
-1]</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:red'> Size: [3,
3]</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:red'> is outside
of buffered region ImageRegion (036C44AC)</span><span lang=EN-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:red'> Dimension:
2</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:red'> Index: [0,
0]</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:red'> Size: [8, 8]</span><span lang=EN-GB><o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-GB>What exactly are these two different
regions? The second one seems to be the whole grid and the first one perhaps the
border. But where do the negative indexes come from?<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-GB>Can someone please take a look at my code?
I suspect there may be something wrong in the definition of the grid, but I can’t
figure out what.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><b><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'># Read the
fixed and moving images using filenames</span></b><span lang=EN-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>FixedImageType =
itk.<u>Image</u>[itk.<u>F</u>,</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:maroon'>2</span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>]</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>MovingImageType
= itk.<u>Image</u>[itk.<u>F</u>,</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:maroon'>2</span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>]</span><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>fixedImageReader
= itk.ImageFileReader[FixedImageType].New()</span><span lang=EN-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>movingImageReader
= itk.ImageFileReader[MovingImageType].New()</span><span lang=EN-GB
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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>fixedImageReader.SetFileName(
r</span><i><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:#00AA00'>'C:\Development\Insight\Examples\Data\RatLungSlice1.<u>mha</u>'</span></i><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'> )</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>movingImageReader.SetFileName(
r</span><i><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:#00AA00'>'C:\Development\Insight\Examples\Data\RatLungSlice2.<u>mha</u>'</span></i><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'>
)</span><span lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>fixedImageReader.Update()</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>movingImageReader.Update()</span><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>fixedImage =
fixedImageReader.GetOutput() </span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'>movingImage = movingImageReader.GetOutput()<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'><o:p> </o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'><o:p> </o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><b><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#
Registration framework</span></b><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>registration
= itk.ImageRegistrationMethod[FixedImageType,MovingImageType].New()</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>metric
= itk.MeanSquaresImageToImageMetric[FixedImageType,MovingImageType].New()</span><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>transformType =
itk.<u>BSplineDeformableTransform</u>[itk.<u>D</u>,</span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:maroon'>2</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'>,</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:maroon'>2</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'>]</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>transform
= transformType.New()</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>optimizer
= itk.<u>LBFGSOptimizer</u>.New()</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>interpolator
= itk.LinearInterpolateImageFunction[FixedImageType,itk.<u>D</u>].New()</span><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>registration.SetOptimizer(
optimizer.GetPointer() )</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>registration.SetTransform(
transform.GetPointer() )</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>registration.SetInterpolator(
interpolator.GetPointer() )</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>registration.SetMetric(
metric.GetPointer() )</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>registration.SetFixedImage(
fixedImage )</span><span lang=EN-GB 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'>registration.SetMovingImage( movingImage
)</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'>registration.SetFixedImageRegion(
fixedImage.GetBufferedRegion() )</span><span style='font-size:10.0pt;
font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><b><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'># Transform
parameters </span></b><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'># Define the
grid</span></b><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>bSplineRegion =
itk.<u>ImageRegion</u>[</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:maroon'>2</span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>]()</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>gridSizeOnImage
= itk.<u>Size</u>[</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:maroon'>2</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:black'>]()</span><span lang=EN-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>gridBorderSize =
itk.<u>Size</u>[</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:maroon'>2</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:black'>]()</span><span lang=EN-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>totalGridSize =
itk.<u>Size</u>[</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:maroon'>2</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:black'>]()</span><span lang=EN-GB
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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>gridSizeOnImage.Fill(
</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:maroon'>5</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:black'> )</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>gridBorderSize.Fill(
</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:maroon'>3</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:black'> )</span><span lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:blue'>for</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'> r </span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:blue'>in</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'>
range(</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:maroon'>0</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:black'>,</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:maroon'>2</span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>):</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>
totalGridSize.SetElement(r,gridSizeOnImage.GetElement(r) +
gridBorderSize.GetElement(r))</span><span lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>bSplineRegion.SetSize(totalGridSize)</span><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>spacing =
fixedImage.GetSpacing()</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>origin =
fixedImage.GetOrigin()</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>fixedImageSize =
fixedImage.GetLargestPossibleRegion().GetSize()</span><span lang=EN-GB
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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:blue'>for</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'> r </span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:blue'>in</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'>
range(</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:maroon'>0</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:black'>,</span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:maroon'>2</span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>):</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>
s = mat.floor((fixedImageSize.GetElement(r) - </span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:maroon'>1</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'> )/(
gridSizeOnImage.GetElement(r) - </span><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New";color:maroon'>1</span><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'> ))</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>
</span><span lang=PT style='font-size:10.0pt;font-family:"Courier New";
color:black'>spacing.SetElement(r,spacing.GetElement(r)*s)</span><span lang=PT
style='font-size:10.0pt;font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><span lang=PT style='font-size:
10.0pt;font-family:"Courier New";color:black'>
origin.SetElement(r,origin.GetElement(r)-s)</span><span lang=PT
style='font-size:10.0pt;font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><span lang=PT 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>gridDirection =
fixedImage.GetDirection()</span><span lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>transform.SetGridSpacing(spacing)</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>transform.SetGridOrigin(origin)</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>transform.SetGridRegion(bSplineRegion)</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>transform.SetGridDirection(gridDirection)</span><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>numbOfParameters
= transform.GetNumberOfParameters()</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>parameters =
itk.<u>Array</u>[itk.<u>D</u>]()</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>parameters.SetSize(numbOfParameters)</span><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>parameters.Fill(
</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:maroon'>0.0</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:black'> )</span><span lang=EN-GB 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-GB
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'>transform.SetParameters(parameters)</span><span
style='font-size:10.0pt;font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'>registration.SetInitialTransformParameters(
transform.GetParameters() )<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'><o:p> </o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'><o:p> </o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><b><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'># Define <u>optimizer</u>
parameters</span></b><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>optimizer.SetGradientConvergenceTolerance(
</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:maroon'>0.05</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:black'> );</span><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:black'>optimizer.SetLineSearchAccuracy(
</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
color:maroon'>0.9</span><span lang=EN-GB style='font-size:10.0pt;font-family:
"Courier New";color:black'> );</span><span lang=EN-GB 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'>optimizer.SetDefaultStepLength( </span><span
style='font-size:10.0pt;font-family:"Courier New";color:maroon'>1.5</span><span
style='font-size:10.0pt;font-family:"Courier New";color:black'> );</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'>optimizer.TraceOn();</span><span
style='font-size:10.0pt;font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'>optimizer.SetMaximumNumberOfFunctionEvaluations( </span><span
style='font-size:10.0pt;font-family:"Courier New";color:maroon'>1000</span><span
style='font-size:10.0pt;font-family:"Courier New";color:black'> );<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'><o:p> </o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'><o:p> </o:p></span></p>
<p class=MsoNormal style='text-autospace:none'><b><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'># Start the
registration process</span></b><span lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:#804000'>#</span></b><span
lang=EN-GB 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-GB
style='font-size:10.0pt;font-family:"Courier New";color:blue'>print</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black'> </span><i><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:#00AA00'>"Starting
registration"</span></i><span lang=EN-GB style='font-size:10.0pt;
font-family:"Courier New"'><o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'>registration.StartRegistration()<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:10.0pt;font-family:"Courier New";
color:black'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-GB>Thank you in advance!<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-GB>Rita<o:p></o:p></span></p>
</div>
</body>
</html>