<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
<br>I would like to make a 3D Affine image registration code, if there is none. I found a similar code in ImageRegistration9.cxx and would like to alter to make it applicable for 3D images. Are these the only changes needed or is there anything else that I need to make. <br>After Line#529, I dont know if there is anything to change.<br>Thanks,<br>
Ganesh<br><br>****<br>Line 146: &nbsp; const&nbsp;&nbsp;&nbsp; unsigned int&nbsp;&nbsp;&nbsp; Dimension = 3;<br><br>Line 285 to 290: Need to remove lines 289, 290 and add 6 more lines to make a (4x3) matrix for 3D images.<br>&nbsp; optimizerScales[4] =&nbsp; 1.0;&nbsp; optimizerScales[5] =&nbsp; 1.0;<br>&nbsp; optimizerScales[6] =&nbsp; 1.0;&nbsp; optimizerScales[7] =&nbsp; 1.0;<br>&nbsp; optimizerScales[8] =&nbsp; 1.0;<br>&nbsp; optimizerScales[9] =&nbsp; translationScale;<br>&nbsp; optimizerScales[10] =&nbsp; translationScale;<br>&nbsp; optimizerScales[11] =&nbsp; translationScale;<br><br>After line 395, add <br>&nbsp; const double finalRotationCenterZ=transform-&gt;GetCenter()[2];<br><br>After line 397, add <br>&nbsp; const double finalTranslationZ =finalParameters[6];<br><br>After line 408, add <br>&nbsp; std::cout &lt;&lt; " Center Z &nbsp;&nbsp;&nbsp;&nbsp; = " &lt;&lt; finalRotationCenterZ&nbsp; &lt;&lt; std::endl;<br><br>After line 410, add <br>&nbsp; std::cout &lt;&lt; " Translation Z = " &lt;&lt; finalTranslationZ&nbsp; &lt;&lt; std::endl;&nbsp; <br><br>Change line 418 into <br>&nbsp; vnl_matrix&lt;double&gt; p(3, 3);<br><br>After line 420 add<br>&nbsp; p[0][2] = (double) finalParameters[2];&nbsp; p[1][0] = (double) finalParameters[3];&nbsp; p[1][1] = (double) finalParameters[4];<br>&nbsp; p[1][2] = (double) finalParameters[5];&nbsp; p[2][0] = (double) finalParameters[6];&nbsp; p[2][1] = (double) finalParameters[7];<br>
&nbsp; p[2][2] = (double) finalParameters[8];<br><br><br>Change line 425 into <br>&nbsp; vnl_matrix&lt;double&gt; r(3, 3);<br><br>Add after line 432 <br>&nbsp; std::cout &lt;&lt; " Scale 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = " &lt;&lt; svd.W(2)&lt;&lt; std::endl;<br><br><br>***<br><br><br><br>                                               <br /><hr />Hotmail: Trusted email with Microsoft’s powerful SPAM protection. <a href='http://clk.atdmt.com/GBL/go/201469226/direct/01/' target='_new'>Sign up now.</a></body>
</html>