<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><BR>
<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid">
<DIV id=yiv919008887>
<DIV>Hi, Luis: </DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp; Thanks for your&nbsp; reply.&nbsp; I attach the result for 3D registration for orginal resolution of images and the sub sample version.&nbsp; The sub sample image factor is 2 for x and Y, and 1 for Z. It means that the each slice is half copmared with the original slice. But the number of slices is not changed. </DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp; Regarding your question about what components are you using ?<BR><BR>&nbsp; 1) Image modality ?<BR>&nbsp; 2) Transform ?<BR>&nbsp; 3) Optimizer ?<BR>&nbsp; 4) Metric ?<BR>&nbsp; 5) Interpolator ?<BR></DIV>
<DIV>I attach part of my code in python </DIV>
<DIV>&nbsp;&nbsp;</DIV>
<DIV><BR><FONT color=#0000ff>registration = itk.MultiResolutionImageRegistrationMethod[itk.Image[itk.F,3],itk.Image[itk.F,3]].New()<BR>imageMetric&nbsp; = itk.MattesMutualInformationImageToImageMetric[itk.Image</FONT></DIV>
<DIV><FONT color=#0000ff>[itk.F,3],itk.Image[itk.F,3]].New()<BR></FONT></DIV>
<DIV><BR><FONT color=#0000ff>transform&nbsp;&nbsp;&nbsp; = itk.VersorRigid3DTransform.D.New()<BR>finalTransform= itk.VersorRigid3DTransform.D.New()<BR>optimizer&nbsp;&nbsp;&nbsp; = itk.RegularStepGradientDescentOptimizer.New()<BR>interpolator = itk.LinearInterpolateImageFunction[itk.Image[itk.F,3],itk.D].New() <BR></FONT></DIV>
<DIV><FONT color=#0000ff></FONT>&nbsp;</DIV>
<DIV><FONT color=#0000ff>&nbsp;I also attach the result repsectively using the orginal volume and sub sample volume </FONT></DIV>
<DIV><FONT color=#0000ff></FONT>&nbsp;</DIV>
<DIV><FONT color=#0000ff>&nbsp;for orginal volume I get reasonably registered result. the result is </FONT></DIV>
<DIV><FONT color=#0000ff></FONT>&nbsp;</DIV>
<DIV><FONT color=#0000ff>Final Registration Parameters <BR>versor X&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -0.001937<BR>versor Y&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.003333<BR>versor Z&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -0.002403<BR>Translation X =&nbsp; -0.250085<BR>Translation Y =&nbsp; -4.207276<BR>Translation Z =&nbsp; -782.659578<BR>Iterations&nbsp;&nbsp;&nbsp; =&nbsp; 45.000000<BR>Metric value&nbsp; = -0.124261</FONT></DIV>
<DIV><FONT color=#0000ff></FONT>&nbsp;</DIV>
<DIV><FONT color=#0000ff>for sub sample volume the registered result is empty, which means only defalut pixel values</FONT></DIV>
<DIV><FONT color=#0000ff>Final Registration Parameters <BR>versor X&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.000695<BR>versor Y&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -0.011173<BR>versor Z&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -0.000622<BR>Translation X =&nbsp; -0.034234<BR>Translation Y =&nbsp; -2.751356<BR>Translation Z =&nbsp; -794.621228<BR>Iterations&nbsp;&nbsp;&nbsp; =&nbsp; 45.000000<BR>Metric value&nbsp; = -0.147362</DIV>
<DIV><BR>&nbsp;&nbsp; So I am wondering where goes wrong?&nbsp; Please let me know if there is anything unclear. Thx </DIV>
<DIV>&nbsp;</DIV>
<DIV>Wen </DIV>
<DIV>&nbsp;</DIV></FONT></DIV></BLOCKQUOTE></td></tr></table><br>