[Insight-users] Is there a limit for the number of Landmarks -Pairs for Thin Plate Spline Transform ?

[Luis Ibanez]

Hi Siwei,

There is no explicit limit to the number of landmarks in your
KernelTransform class.

The only limitation may arise from the VNL mathematical libraries
that are used for computing the matrix operations of the Transform.





Please do the following:


1) Compile your program in Debug mode,

2) Run it again and interrupt it as soon as you
    start seeing the output in the console.

3) Then look at the callstack.
    (if possible post it to the mailing list).



 From it we should be able to figure out where the
messages are comming from.



    Thanks


      Luis


------------------
Siwei Yang wrote:

> Hi every body,
> when I 'was using itk::ThinPlateSplineKernelTransform(actually include 
> all the classes inherited from itk::KernelTransform), I found ,that 
> there might be  an upper limit for the number of landmark-pairs. This 
> limit is about 150-200 pairs. Sorry , that I can 't say exactly, how much.
> If I exceed this limt, for example 300, then I will get huge amount of 
> standard outputs, which almost floods my screen, just like 
> ............................................................
>        0 52.00000        0 53.00000        0 53.00000        0 
> 54.00000        0 55.00000        0 55.00000        0 56.00000        0 
> 56.00000        0 57.00000        0 57.00000        0 58.00000        0 
> 58.00000        0 59.00000        0 59.00000        0 60.00000        0 
> 60.00000        0 61.00000        0 61.00000        0 61.00000        0 
> 62.00000        0 63.00000        0 64.00000        0 64.00000        0 
> 64.00000        0 65.00000        0 66.00000        0 66.00000        0 
> 67.00000        0 68.00000        0 69.00000        0 70.00000        0 
> 70.00000        0 71.00000        0 72.00000        0 73.00000        0 
> 73.00000        0 74.00000        0 75.00000        0 76.00000        0 
> 77.00000        0 77.00000        0 78.00000        0 79.00000        0 
> 81.00000        0 82.00000        0 83.00000        0 83.00000        0 
> 84.00000        0 85.00000        0 87.00000        0 88.00000        0 
> 89.00000        0 90.00000        0 91.00000        0 92.00000        0 
> 93.00000        0 94.00000        0 95.00000        0 96.00000        0 
> 98.00000        0 99.00000        0 100.00000        0 101.00000        
> 0 102.00000        0 1.00000        0  1.00000        0  1.00000        
> 0  1.00000        0  1.00000        0 
> .................................................................
> ..................................................................................................................................................................................................................................................................................................................................... 
> 
> and I don't know , where they come from. Finally I 'll obtain a wrong 
> interpolation.
> Because my application uses usually more than 3000 pairs landmarks (3D), 
> it would be ideal, if this could be allowed.
> Did anyone have this problem? How can I do ?
> 
> 
> Thank you very much
> best wishes
> 
> Siwei
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
> 
>