[Insight-users] FEM registration plus landmarks crash

[Markus Weigert]

Dear insight-users,

whenever I use the FEMRegistrationFilter with
user supplied landmarks (3D), the program crashes
if there was at least one landmark correspondence in the file.

The file with landmarks had only one entry:

<LoadLandmark>
3 259 260 32    3 261 260 28   0.1


I used valgrind to track the error and it seems to raise
at "itk::fem::SolverCrankNicolson::AssembleKandM()".

Did somebody experience the same problem???

I use itk3.2.0 on a suse linux 10.2 machine.

Valgrind produced the following output:

Reading config file.../data/Weigert/RIP/FusionJ/FEMConfig.txt
Example configured. E 10000 1000 500 rho 10000 1000 500
 beginning level 0
 scaling 8
 scaling 8
 scaling 8
 ElementsPerDim 16 16 1.5
 generating regular mesh
 generating regular mesh done
 DO NOT init interpolation grid : im sz 65 65 7 MeshSize 64 64 6
 applying loads
Landmarks.txt
Try loading landmarks...
Try loading landmarks...
 num of LM loads 1
 prescale Pt 259 260 32
 postscale Pt 32.375 32.5 4 scale 8
 landmarks done
 which node 1
 edge coord 4 0 0
 which node 1
 edge coord 4 0 0
 which node 1
 edge coord 4 0 0
 which node 3
 edge coord 0 4 0
 which node 3
 edge coord 0 4 0
 which node 3
 edge coord 0 4 0
 which node 4
 edge coord 0 0 6
 which node 4
 edge coord 0 0 6
 which node 4
 edge coord 0 0 6
 node 1
 allocating deformation field
 load sizes [64, 64, 6]  image [64, 64, 6]
==29084==
==29084== Invalid read of size 4
==29084==    at 0x8927ED2: itk::fem::SolverCrankNicolson::AssembleKandM() 
(in /data/Weigert/RIP/FusionJ/RegistrationImagingPlattform)
==29084==    by 0x859BFB0: 
FEMRegFrameworkImp::FEMRegistrationFilter<itk::Image<float, 3>, 
itk::Image<float, 3> >::MultiResSolve() 
(in /data/Weigert/RIP/FusionJ/RegistrationImagingPlattform)
==29084==    by 0x859C98A: 
FEMRegFrameworkImp::FEMRegistrationFilter<itk::Image<float, 3>, 
itk::Image<float, 3> >::RunRegistration() 
(in /data/Weigert/RIP/FusionJ/RegistrationImagingPlattform)
Segmentation fault


Best regards,
Markus


-----------------------------------------------------------------------------------------
-----------------------------------------------------------------------------------------
Forschungszentrum Juelich GmbH
52425 Juelich

Sitz der Gesellschaft: Juelich
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
Vorsitzende des Aufsichtsrats: MinDirig'in Baerbel Brumme-Bothe
Vorstand: Prof. Dr. Achim Bachem (Vorsitzender), Dr. Ulrich Krafft (stellv. 
Vorsitzender)
-----------------------------------------------------------------------------------------
-----------------------------------------------------------------------------------------