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><div class="gmail_quote"><div>Below are my configuration file (I made sure the file names existed in the directory before running) and the error I got. <br>
<br><span style="color: rgb(0, 0, 153);">% Parameters for the single- or multi-resolution techniques</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% ---------------------------------------------------------</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">1 % Number of levels in the multi-res pyramid (1 = single-res)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">1 % Highest level to use in the pyramid</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">1 1 1 % Scaling at lowest level of pyramid</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">3 3 3 % Number of pixels per element</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">1.e4 1.e4 1.e4 % Elasticity (E)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">1.e5 1.e5 1.e5 % Density x capacity (RhoC)</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">1 % Image energy scaling (gamma) - sets gradient step size 3 3 3 % NumberOfIntegrationPoints</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">3 % WidthOfMetricRegion</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">30 % MaximumIterations</span><br style="color: rgb(0, 0, 153);"><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% -------------------------------</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">% Parameters for the registration</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% -------------------------------</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">0 0.99 % Similarity metric (0=mean sq, 1 = ncc, 2=pattern int, 3=MI, 5=demons)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">1.0 % Alpha</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">0 % DescentDirection (1 = max, 0 = min)</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">0 % DoLineSearch (0=never, 1=always, 2=if needed)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">1.e1 % TimeStep</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">0.5 % Landmark variance</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">0 % Employ regridding / enforce diffeomorphism ( >= 1 -> true)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% ----------------------------------</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">% Information about the image inputs</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% ----------------------------------</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">101 % Nx (image x dimension)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">101 % Ny (image y dimension)</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">30 % Nz (image z dimension - not used if 2D)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">G:\ITK_test\deform11\deform11\Debug\A.hdr % ReferenceFileName </span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">G:\ITK_test\deform11\deform11\Debug\B.hdr % TargetFileName</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% -------------------------------------------------------------------</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">% The actions below depend on the values of the flags preceding them.</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% For example, to write out the displacement fields, you have to set</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">% the value of WriteDisplacementField to 1.</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">% -------------------------------------------------------------------</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">0 % UseLandmarks? - read the file name below if this is true</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">- % LandmarkFileName</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">G:\ITK_test\deform11\deform11\Debug\D % ResultsFileName (prefix only)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">1 % WriteDisplacementField?</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">G:\ITK_test\deform11\deform11\Debug\D % DisplacementsFileName (prefix only)</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">0 % ReadMeshFile?</span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">- % MeshFileName</span><br style="color: rgb(0, 0, 153);">
<span style="color: rgb(0, 0, 153);">END</span><br><br>Error message I got was:<br><span style="color: rgb(255, 0, 0);">Reading config file...W.txt </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">Example configured. E 0 rho 0 </span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> reading moving </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);"> reading fixed </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">Exception caught during reference file reading </span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">itk::ImageFileReaderException (015BF9CC) </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">Location: "void __thiscall itk::ImageFileReader<class itk::Image<unsigned char,3>,class itk::DefaultConvertPixelTraits<unsigned char> >::GenerateOutputInformation(void)" </span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);">File: c:\itk\src\code\io\itkImageFileReader.txx </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">Line: 100 </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">Description: FileName must be specified </span><br>
<br></div></div>