Thanks Vince!! It seems to work fine. But I gotta play with the parameters to get a better registration result. <br><br>Does anyone have any suggestion how the parameters should be chosen or they simply vary from case to case?<br>
<br>Thanks,<br>stephen<br><br><div class="gmail_quote">On Wed, Apr 14, 2010 at 7:04 PM, Vincent Magnotta <span dir="ltr"><<a href="mailto:vincent-magnotta@uiowa.edu">vincent-magnotta@uiowa.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div>
<font face="Calibri, Verdana, Helvetica, Arial"><span style="font-size: 11pt;">Stephen,<br>
<br>
The problem was with your parameter file. The attached file worked for me in 3D.<div class="im"><br>
<br>
Vince<br>
<br>
<br>
<br>
On 4/12/10 12:59 PM, "Stephen Yip" <<a href="http://stephen.fyip1@gmail.com" target="_blank">stephen.fyip1@gmail.com</a>> wrote:<br>
<br>
</div></span></font><blockquote><font face="Calibri, Verdana, Helvetica, Arial"><span style="font-size: 11pt;"><div><div></div><div class="h5">Hi Luis, <br>
<br>
Just a quick update. The B-spline 3d registration seemed to work well except it shifted the origin of my image. But I am still having the some issues with the 3d FEM registration. It doesn't run. I am not sure if it is due to the memory or configuration file problem (see below). Thank you very much. Please see below for details. <br>
<br>
By the way, I have a 2 gig ram computer and 32 bit window. Is it a problem for 3d FEM?<br>
<br>
Below are my configuration file (I made sure the file names existed in the directory before running) and the error I got. <br>
<br>
% Parameters for the single- or multi-resolution techniques<br>
% ---------------------------------------------------------<br>
1 % Number of levels in the multi-res pyramid (1 = single-res)<br>
1 % Highest level to use in the pyramid<br>
1 1 1 % Scaling at lowest level of pyramid<br>
3 3 3 % Number of pixels per element<br>
1.e4 1.e4 1.e4 % Elasticity (E)<br>
1.e5 1.e5 1.e5 % Density x capacity (RhoC)<br>
1 % Image energy scaling (gamma) - sets gradient step size 3 3 3 % NumberOfIntegrationPoints<br>
3 % WidthOfMetricRegion<br>
30 % MaximumIterations<br>
<br>
% -------------------------------<br>
% Parameters for the registration<br>
% -------------------------------<br>
0 0.99 % Similarity metric (0=mean sq, 1 = ncc, 2=pattern int, 3=MI, 5=demons)<br>
1.0 % Alpha<br>
0 % DescentDirection (1 = max, 0 = min)<br>
0 % DoLineSearch (0=never, 1=always, 2=if needed)<br>
1.e1 % TimeStep<br>
0.5 % Landmark variance<br>
0 % Employ regridding / enforce diffeomorphism ( >= 1 -> true)<br>
% ----------------------------------<br>
% Information about the image inputs<br>
% ----------------------------------<br>
101 % Nx (image x dimension)<br>
101 % Ny (image y dimension)<br>
30 % Nz (image z dimension - not used if 2D)<br>
G:\ITK_test\deform11\deform11\Debug\A.hdr % ReferenceFileName <br>
G:\ITK_test\deform11\deform11\Debug\B.hdr % TargetFileName<br>
% -------------------------------------------------------------------<br>
% The actions below depend on the values of the flags preceding them.<br>
% For example, to write out the displacement fields, you have to set<br>
% the value of WriteDisplacementField to 1.<br>
% -------------------------------------------------------------------<br>
0 % UseLandmarks? - read the file name below if this is true<br>
- % LandmarkFileName<br>
G:\ITK_test\deform11\deform11\Debug\D % ResultsFileName (prefix only)<br>
1 % WriteDisplacementField?<br>
G:\ITK_test\deform11\deform11\Debug\D % DisplacementsFileName (prefix only)<br>
0 % ReadMeshFile?<br>
- % MeshFileName<br>
END<br>
<br>
Error message I got was:<br>
Reading config file...W.txt <br>
Example configured. E 0 rho 0 <br>
reading moving <br>
reading fixed <br>
Exception caught during reference file reading <br>
<br>
itk::ImageFileReaderException (015BF9CC) <br>
Location: "void __thiscall itk::ImageFileReader<class itk::Image<unsigned char,3>,class itk::DefaultConvertPixelTraits<unsigned char> >::GenerateOutputInformation(void)" <br>
File: c:\itk\src\code\io\itkImageFileReader.txx <br>
Line: 100 <br>
Description: FileName must be specified <br>
<br>
<br>
</div></div><hr width="95%" align="CENTER" size="3"></span></font><div class="im"><font size="2"><font face="Consolas, Courier New, Courier"><span style="font-size: 10pt;">_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
</span></font></font></div></blockquote><div class="im"><font size="2"><font face="Consolas, Courier New, Courier"><span style="font-size: 10pt;"><br>
</span></font></font><font face="Calibri, Verdana, Helvetica, Arial"><span style="font-size: 11pt;">----------------------<br>
Associate Professor<br>
Department of Radiology<br>
0453-D JCP<br>
200 Hawkins Drive<br>
Iowa City, IA 52242<br>
E-mail: <a href="http://vincent-magnotta@uiowa.edu" target="_blank">vincent-magnotta@uiowa.edu</a><br>
Phone: 319-356-8255 Fax: 319-353-6275<br>
Website: <a href="http://www.radiology.uiowa.edu" target="_blank">http://www.radiology.uiowa.edu</a><br>
<br>
<br>
</span></font>
</div></div>
</blockquote></div><br>