[Insight-users] segfault in multires FEMRegistration

[Felix Bollenbeck]

Hello everybody,

Im struggeling along with the FEM registration method and i cant proceed 
to multires registration sheme.
I used the default parameter file so far, which uses only one resolution 
level. I changed the parameter file according to entries in this mailing 
list

after allocating the deformtion field, the program exits with a 
segmentation fault.

can anybody give me a hint what might be my mistake here?

Cheers

-Felix.

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %
 % Parameters for the single- or multi-resolution techniques
 %

 3        % Number of levels in the multi-res pyramid (1 = single-res)
 3        % Highest level to use in the pyramid
 8    % Scaling at lowest level of pyramid
 4 4 4    % Number of pixels per element
 1.e5 1.e4 1.e4 % Elasticity (E)
 1.e4 1.e4 1.e4 % Density x capacity (RhoC)
 1 1 1    % Image energy scaling (gamma) - sets gradient step size
 2 2 2    % NumberOfIntegrationPoints
 4 4 4    % WidthOfMetricRegion
 10 1 1   % MaximumIterations
% -------------------------------
% Parameters for the registration
% -------------------------------
0 0.99  % Similarity metric (0=mean sq, 1 = ncc, 2=pattern int, 3=MI, 
5=demons)
1.0    % Alpha
0    % DescentDirection (1 = max, 0 = min)
0    % DoLineSearch (0=never, 1=always, 2=if needed)
1.e1    % TimeStep
0.5     % Landmark variance
0       % Employ regridding / enforce diffeomorphism ( = 1 - true)
% ----------------------------------
% Information about the image inputs
% ----------------------------------
1800    % Nx (image x dimension)
1500    % Ny (image y dimension)
0      % Nz (image z dimension - not used if 2D)
../../Insight/Examples/Data/RatLungSlice1.mha  % ReferenceFileName
../../Insight/Examples/Data/RatLungSlice2.mha  % TargetFileName
% -------------------------------------------------------------------
% The actions below depend on the values of the flags preceding them.
% For example, to write out the displacement fields, you have to set
% the value of WriteDisplacementField to 1.
% -------------------------------------------------------------------
0    % UseLandmarks? - read the file name below if this is true
-    % LandmarkFileName
./result                       % ResultsFileName (prefix only)
1       % WriteDisplacementField?
./disp                         % DisplacementsFileName (prefix only)
0       % ReadMeshFile?
-                                      % MeshFileName
END



-- 
----------------------------------------------------------------
contact:
bollen at ipk-gatersleben.de
Leibniz-Institute of Plant Genetics and Crop Plant Research (IPK)
Pattern Recognition Group (PRG) / Department of Cytogenetics
Corrensstr. 3, D-06466 Gatersleben, Germany.   (Room A-0.06)
Tel./Fax: ++49 (0) 39482  - 5 344 / - 5 137.