<HTML>
<HEAD>
<TITLE>Re: [Insight-users] 3d deformable registration</TITLE>
</HEAD>
<BODY>
<FONT FACE="Calibri, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'>The problem is with the parsing of the configuration file. There is a problem in the method ReadConfigFile() in <BR>
itk::FEMRegistrationFilter. I have verified this with some print statements, but I have not yet figured out the problem. I will look into this further over the next couple of days.<BR>
<BR>
Vince<BR>
<BR>
<BR>
On 4/12/10 12:59 PM, "Stephen Yip" <<a href="stephen.fyip1@gmail.com">stephen.fyip1@gmail.com</a>> wrote:<BR>
<BR>
</SPAN></FONT><BLOCKQUOTE><FONT FACE="Calibri, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'>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>
<HR ALIGN=CENTER SIZE="3" WIDTH="95%"></SPAN></FONT><FONT SIZE="2"><FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'>_____________________________________<BR>
Powered by www.kitware.com<BR>
<BR>
Visit other Kitware open-source projects at<BR>
<a href="http://www.kitware.com/opensource/opensource.html">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">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">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">http://www.itk.org/mailman/listinfo/insight-users</a><BR>
</SPAN></FONT></FONT></BLOCKQUOTE><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="vincent-magnotta@uiowa.edu">vincent-magnotta@uiowa.edu</a><BR>
Phone: 319-356-8255 Fax: 319-353-6275<BR>
Website: <a href="http://www.radiology.uiowa.edu">http://www.radiology.uiowa.edu</a><BR>
<BR>
<BR>
</SPAN></FONT>
</BODY>
</HTML>