Hi Aviv, Thank you very much for your help.<br><br>Jian<br><br><br>Date: Mon, 2 Jun 2008 22:02:09 +0300<br>
From: &quot;Aviv Hurvitz&quot; &lt;<a href="mailto:aviv.hurvitz@gmail.com">aviv.hurvitz@gmail.com</a>&gt;<br>
Subject: Re: [Insight-users] Is there way to speed up<br>
 &nbsp; &nbsp; &nbsp; &nbsp;cross-correlation &nbsp; &nbsp; &nbsp; based deformable registration?<br>
To: <a href="mailto:insight-users@itk.org">insight-users@itk.org</a><br>
Message-ID:<br>
 &nbsp; &nbsp; &nbsp; &nbsp;&lt;<a href="mailto:e25acca30806021202l50d68fa3j4a494ac5b4685b15@mail.gmail.com">e25acca30806021202l50d68fa3j4a494ac5b4685b15@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
See if you can use the elastix program (which is based on ITK code).<br>
<a href="http://elastix.isi.uu.nl/index.php" target="_blank">http://elastix.isi.uu.nl/index.php</a><br>
<br>
It has several speed improvements:<br>
1. You can specify samples.<br>
2. It utilizes the sparse Jacobian of the B-spline transform for faster<br>
computation of the metric gradient.<br>
3. You can use a multiple resolution strategy, both for the images and for<br>
the B-spline transform itself.<br>
<br>
- Aviv<br>
<br>
<br>
On Mon, Jun 2, 2008 at 9:20 PM, Jian Yang &lt;<a href="mailto:yaland1977@gmail.com">yaland1977@gmail.com</a>&gt; wrote:<br>
<br>
&gt; Hi all,<br>
&gt;<br>
&gt; I&#39;m currently working on deformable registration of CT images. I&#39;m using<br>
&gt; cross-correlation (DeformableRegistration6.cxx) and mutual information<br>
&gt; (DeformableRegistration15.cxx) metrics with a B-Spline Deformable Transform<br>
&gt; in a 300*208*241 3D image.<br>
&gt;<br>
&gt; It seems that the example DeformableRegistration15.cxx &nbsp;works much faster<br>
&gt; than example DeformableRegistration6.cxx. I guess example<br>
&gt; DeformableRegistration6.cxx &nbsp;is time consuming since it is performed on the<br>
&gt; whole image. And the Mattes Mutual Information metric works faster because<br>
&gt; it allows to take samples (SetNumberOfSpatialSamples) in the image instead<br>
&gt; of taking all voxels. So, I am wondering is there way to do the same thing<br>
&gt; in the cross-correlation metric to make it faster?<br>
&gt;<br>
&gt; Thank you very much for your reply.<br>
&gt;<br>
&gt;<br>
&gt; Best Regards,<br>
&gt;<br>
&gt;<br>
&gt; Jian