[Insight-users] Re: Need help with itkPointSetToPointSetRegistration

[Luis Ibanez]

Hi Eve,

There seems to be an actual bug on the way the MultivaluedVnlOptimizer
is reporting progress updates.

We recently add a similar Command observer the example:


          Examples/Patented/
              IterativeClosestPoint1.cxx


and it is not reporting the values as expected either.

Please add this problem as a bug entry to the bug tracker:

     http://public.kitware.com/Bug/index.php


let us know if you find any problem creating the bug entry.



    Thanks


      Luis



-----------------
Eve Heyes wrote:
> Hi Luis,
> 
> I am trying to observe the evolution of the optimization. I have used an observer to update the position values (see below). However, nothing was observed; as if the optimization ended in the first iteration. I have used very few sets of points; could this be why?
> 
> Thanks,
> Eve
> 
> 
> Set1
> 
> -0.268955 5.36033 30.5068 
> -0.268955 5.08468 27.3331 
> -0.816196 8.54178 24.9906 
> -10.7404 2.91095 25.7686 
> 10.8408 2.91095 25.7686 
> -4.81812 11.6438 25.7686 
> 4.81812 11.7442 25.7686 
> 
> Set2
> 
> -0.0752423 12.2746 -82.7269 
> 24.9043 11.6886 -85.8659 
> 0.196281 14.329 -88.1486 
> -10.8841 9.06179 -87.67 
> 11.1226 8.97906 -87.67 
> -4.59649 18.6587 -87.67 
> 4.58675 18.6587 -87.67 
> 
> 
> class CommandIterationUpdate : public itk::Command 
> {
> public:
>   typedef  CommandIterationUpdate   Self;
>   typedef  itk::Command             Superclass;
>   typedef itk::SmartPointer<Self>  Pointer;
>   
>   itkNewMacro( Self );
> protected:
>   CommandIterationUpdate() {};
> public:	
> 
>   typedef itk::LevenbergMarquardtOptimizer OptimizerType;
>   typedef  const OptimizerType   *    OptimizerPointer;
>   
>   void Execute(itk::Object *caller, const itk::EventObject & event)
>     {
>       Execute( (const itk::Object *)caller, event);
>     }
> 
>   void Execute(const itk::Object * object, const itk::EventObject & event)
>   {      
>         OptimizerPointer optimizer = dynamic_cast< OptimizerPointer >( object );
>       
> 	    if( typeid( event ) != typeid( itk::IterationEvent ) )
>         {
> 			return;
>         }
> 		
> 	//	std::cout << optimizer->GetValue();
> 		std::cout << ") ";
> 		std::cout << optimizer->GetCurrentPosition ();		
>   }
> 		
> };
> 
> 
> 
> Could this be due to the simplicity of the problem?
> 
> --- luis.ibanez at kitware.com wrote:
> 
> From: Luis Ibanez <luis.ibanez at kitware.com>
> To:  evelyn_heyes at homeworking.org
> CC: Insight Users <insight-users at itk.org>
> Subject: Re: Need help with itkPointSetToPointSetRegistration
> Date: Wed, 11 Apr 2007 15:18:03 -0400
> 
> 
> Hi Evelyn,
> 
> What version of ITK are you using ?
> 
> What event are you observing ?
> 
> 
> Please let us know,
> 
> 
>    Thanks
> 
> 
>       Luis
> 
> 
> --------------------------------------
> evelyn_heyes at homeworking.org wrote:
> 
>>Hello Luis,
>>
>>Re: IterativeClosestPoint2 Example
>>
>>I have added an event observer to the LM optimizer but failed to observe any events. What is the cause of this? I noted that the same problem was posted at http://public.kitware.com/pipermail/insight-users/2006-May/018044.html. Has a solution been found? 
>>
>>Thanks,
>>Eve
>>
>>
>>
>>
>>
>>
> 
> 
> 
> 
> 
> _____________________________________________________________
> Visit http://www.homeworking.com ~ To help you work at home
>