[Insight-users] FEM warping image error

[verenas at uni-koblenz.de]

Hello,

i'm using DeformableRegistration1.cxx for multimodal registration.
Everything works fine ecceptly one strange error. I't says "Process
aborted. Segmentation fault" while trying to warp the image. The output
image shows the same as the moving image before.

The parameter file:

1	% Number of levels in the multi-res pyramid (1 = single-res)
1	% Highest level to use in the pyramid
1 1 1  	% Scaling at lowest level of pyramid
10 	% Number of pixels per element
1.e4    % Elasticity (E)
1.e4        % Density x capacity (RhoC)
1	% Image energy scaling (gamma) - sets gradient step size
2	% NumberOfIntegrationPoints
1	% WidthOfMetricRegion
20	% MaximumIterations
% -------------------------------
% Parameters for the registration
% -------------------------------
3 0.99  % Similarity metric (0=mean sq (min), 1 = ncc (max), 2=pattern
int, 3=MI (max), 5=demons)
1.0	% Alpha
1	% DescentDirection (1 = max, 0 = min)
2	% DoLineSearch (0=never, 1=always, 2=if needed)2
1.e1	% TimeStep 1.e1
0.5     % Landmark variance
1       % Employ regridding / enforce diffeomorphism ( >= 1 -> true)
% ----------------------------------
% Information about the image inputs
% ----------------------------------
256	% Nx (image x dimension)
256	% Ny (image y dimension)
20	% Nz (image z dimension - not used if 2D)
.../CTAmov.hdr % ReferenceFileName
.../A_MR.hdr % 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
.../CTmov_result                       % ResultsFileName (prefix only)
1       % WriteDisplacementField?
./CTmov__disp                         % DisplacementsFileName (prefix only)
0       % ReadMeshFile?
-                                      % MeshFileName
END


thx
verena