<html><div style='background-color:'><DIV class=RTE>
<P>I am developing an aplication for endoscopy test of cancer. I am trying to use BSplineTransform, BSplineInterpolate, LBFGSBOptimizer and <FONT size=2>MeanSquaresImageToImageMetric ok.</FONT></P>
<P><FONT size=2>Why I am using that (BSplineTransform, BSplineInterpolate, LBFGSBOptimizer and <FONT size=2>MeanSquaresImageToImageMetric </FONT>)? because I need DEFORMABLEREGISTRATION and the SoftwareGuide says that to do DeformableRegistration they (BSplineTransform, BSplineInterpolate, LBFGSBOptimizer and <FONT size=2>MeanSquaresImageToImageMetric </FONT>) are the way.</P>
<P></FONT>I am working with images *.bmp and *.tif 537X418 pixels and 222KB.</P>
<P>when I saw the DeformableRegistration4,5,6,7,8 (InsightToolkit-2.0.0\Examples\Registration) they say me the way to use BSplineTransform so I change DeformableRegistration7 to 2D images, but the OutPutImage was like MovingImage always. I send my question and you tell me the reason was I did not help the metric (<FONT size=2>MeanSquaresImageToImageMetric</FONT>) so I said you How can I help <FONT size=2>MeanSquaresImageToImageMetric?.</FONT></P>
<P><FONT size=2>Again I send the code but yesterday I find a instruction SetComputeGradient() and the OutPutImage is like FixedImage but is not very well.</FONT></P>
<P><FONT size=2>Ha I do not remeber sometime you say me if I was using MutualInformationImageToImageMetric and I would like to know how to use but I know the answer SoftwareGuide so I saw the Guide and ImageRegistration2.cxx and tell me if I want to use it (MutualInformationImageToImageMetric) alwas do i need <FONT face=Courier size=1>itkGradientDescentOptimizer.h, itkNormalizeImageFilter.h itkDiscreteGaussianImageFilter.h ??</P></FONT></FONT>Thank you for you Help.<BR><BR><BR>>From: Luis Ibanez <luis.ibanez@kitware.com><BR>>To: javier silva bravo <javier_silva_bravo@hotmail.com><BR>>CC: insight-users@itk.org<BR>>Subject: Re: [Insight-users] Where I can find examples or paper about<BR>>Date: Wed, 18 May 2005 21:37:28 -0400<BR>><BR>><BR>><BR>>Hi Javier,<BR>><BR>><BR>><BR>>Please be more specific.<BR>><BR>>Our divinatory skills
are quite limited and<BR>>tend to fade as we are rapidly aging.<BR>><BR>><BR>><BR>><BR>><BR>>When you say:<BR>><BR>><BR>><BR>> "It does not work"<BR>><BR>><BR>><BR>>Do you mean:<BR>><BR>><BR>>1) It doesn't compile ?<BR>>2) It compiles but does not link ?<BR>>3) It links but produces segmentation faults at run time ?<BR>>4) It throws exceptions at run time ?<BR>>5) It runs without crashing but the registration result<BR>> does not aligns the images as expected ?<BR>><BR>><BR>><BR>>We will be more than happy to help you, but we need your<BR>>help in providing enough information about your problem.<BR>><BR>><BR>><BR>><BR>>BTW, note that *ALL* the code shown in the software guide<BR>>is
available in the directory Insight/Examples<BR>><BR>>In particular... the image registration examples are available<BR>>in<BR>><BR>><BR>> Insight/Examples/Registration<BR>><BR>><BR>>You could happily go straight for the files:<BR>><BR>><BR>> ImageRegistration14.cxx<BR>> ImageRegistration15.cxx<BR>> ImageRegistration16.cxx<BR>><BR>><BR>>where that metric is used.<BR>><BR>><BR>><BR>><BR>>These examples are compiled and executed every night in 55<BR>>different platforms ranging from cygwin, to Linux, to
IRIX.<BR>><BR>>http://www.itk.org/Testing/Dashboard/20050517-0500-Nightly/TestDetail/__Examples_Registration_ImageRegistration14Test.html<BR>><BR>><BR>><BR>><BR>><BR>>Regards,<BR>><BR>><BR>><BR>> Luis<BR>><BR>><BR>><BR>><BR>>---------------------------<BR>>javier silva bravo wrote:<BR>><BR>>>Yes Luis but I don undertand why it does not work.<BR>>><BR>>>In the guide says:<BR>>><BR>>>typedef itk::MutualInformationImageToImageMetric<<BR>>><BR>>>InternalImageType,<BR>>><BR>>>InternalImageType > MetricType;<BR>>><BR>>>MetricType::Pointer metric = MetricType::New();<BR>>>registration->SetMetric( metric );<BR>>><BR>>>metric->SetFixedImageStandardDeviation( 0.4 );<BR>>>metric->SetMovingImageStandardDeviation( 0.4
);<BR>>>metric->SetNumberOfSpatialSamples( 50 );<BR>>><BR>>>typedef itk::NormalizeImageFilter<<BR>>>FixedImageType,<BR>>>InternalImageType<BR>>> > FixedNormalizeFilterType;<BR>>>typedef itk::NormalizeImageFilter<<BR>>>MovingImageType,<BR>>>InternalImageType<BR>>> > MovingNormalizeFilterType;<BR>>>FixedNormalizeFilterType::Pointer fixedNormalizer =<BR>>>FixedNormalizeFilterType::New();<BR>>>MovingNormalizeFilterType::Pointer movingNormalizer =<BR>>>MovingNormalizeFilterType::New();<BR>>><BR>>>fixedNormalizer->SetInput( fixedImageReader->GetOutput() );<BR>>><BR>>>movingNormalizer->SetInput( movingImageReader->GetOutput() );<BR>>><BR>>>Now tell me what am i doing wrong?<BR>>><BR>>>typedef
itk::MutualInformationImageToImageMetric<<BR>>> FixedImageType,<BR>>> MovingImageType> MetricType;<BR>>> typedef itk::NormalizeImageFilter<<BR>>>
FixedImageType,<BR>>> InternalImageType> <BR>>>FixedNormalizeFilterType;<BR>>> typedef itk::NormalizeImageFilter<<BR>>> MovingImageType,<BR>>> InternalImageType>
<BR>>>MovingNormalizeFilterType;<BR>>><BR>>>FixedNormalizeFilterType::Pointer fixedNormalizer = <BR>>>FixedNormalizeFilterType::New();<BR>>> MovingNormalizeFilterType::Pointer movingNormalizer = <BR>>>MovingNormalizeFilterType::New();<BR>>><BR>>> MetricType::Pointer metric = MetricType::New();<BR>>><BR>>> fixedNormalizer->SetInput( fixedImage );<BR>>> movingNormalizer->SetInput( movingImageReader->GetOutput() );<BR>>><BR>>> registration->SetFixedImage( fixedNormalizer->GetOutput() );<BR>>> registration->SetMovingImage( movingNormalizer->GetOutput() );<BR>>>
metric->SetFixedImage( fixedNormalizer->GetOutput() );<BR>>> metric->SetMovingImage( movingNormalizer->GetOutput() );<BR>>> interpolator->SetInputImage(movingImageReader->GetOutput());<BR>>><BR>>> FixedImageType::RegionType fixedRegion = <BR>>>fixedImage->GetBufferedRegion();<BR>>> registration->SetFixedImageRegion( fixedRegion );<BR>>> metric->SetFixedImageRegion(fixedRegion);<BR>>><BR>>>metric->SetTransformParameters( transform->GetParameters() );<BR>>><BR>>>metric->SetInterpolator( interpolator
);<BR>>><BR>>>metric->SetFixedImageStandardDeviation(0.49);<BR>>><BR>>>metric->SetMovingImageStandardDeviation(0.49);<BR>>><BR>>>metric->SetNumberOfSpatialSamples( 50 );<BR>>><BR>>>metric->SetComputeGradient(1);<BR>>><BR>>><BR>>><BR>>><BR>>><BR>>><BR>>><BR>>><BR>>> >From: Luis Ibanez <luis.ibanez@kitware.com><BR>>> >To: javier silva bravo <javier_silva_bravo@hotmail.com><BR>>> >CC: insight-users@itk.org<BR>>> >Subject: Re: [Insight-users] Where I can find examples or paper <BR>>>about<BR>>> >Date: Wed, 18 May 2005 17:51:51 -0400<BR>>> ><BR>>> ><BR>>> >Hi
Javier,<BR>>> ><BR>>> >You will find detailed information about both of your<BR>>> >questions in the ITK Software Guide:<BR>>> ><BR>>> ><BR>>> > http://www.itk.org/ItkSoftwareGuide.pdf<BR>>> ><BR>>> ><BR>>> >The guide describe examples on how to use the ImageSeries<BR>>> >and how to do multi-modality registration.<BR>>> ><BR>>> ><BR>>> >The code for this examples is avalable in the directory<BR>>> ><BR>>> ><BR>>> >
Insight/Examples<BR>>> ><BR>>> ><BR>>> ><BR>>> ><BR>>> > Regards,<BR>>> ><BR>>> ><BR>>> > Luis<BR>>> ><BR>>> ><BR>>> >--------------------------<BR>>> >javier silva bravo wrote:<BR>>> ><BR>>> >>itk::NumericSeriesFileNames<BR>>> >>I need read and write a secuence of Images.<BR>>> >>and where i can find too information about<BR>>> >>MuttualInformationImagetoImageMetric? because I dont know how <BR>>>to<BR>>> >>use
them.<BR>>> >> Thank you.<BR>>> >><BR>>> <BR>>> >>------------------------------------------------------------------------<BR>>> >>MSN Premium. Protégete, Comunícate y Diviértete Haz clic aquí<BR>>> >><http://g.msn.com/8HMBESMX/2743??PS=47575><BR>>> >><BR>>> >><BR>>> <BR>>> >>------------------------------------------------------------------------<BR>>> >><BR>>> >>_______________________________________________<BR>>> >>Insight-users mailing
list<BR>>> >>Insight-users@itk.org<BR>>> >>http://www.itk.org/mailman/listinfo/insight-users<BR>>> ><BR>>> ><BR>>> ><BR>>><BR>>>------------------------------------------------------------------------<BR>>>T1msn Hotmail Plus. Más espacio, más funcional Haz clic aquí <BR>>><http://g.msn.com/8HMBESMX/2755??PS=47575><BR>>><BR>>><BR>>>------------------------------------------------------------------------<BR>>><BR>>>_______________________________________________<BR>>>Insight-users mailing list<BR>>>Insight-users@itk.org<BR>>>http://www.itk.org/mailman/listinfo/insight-users<BR>><BR>><BR>><BR></DIV></div><br clear=all><hr>T1msn Search. Todo lo que buscas ahora más rápido <a href="http://g.msn.com/8HMBESMX/2728??PS=47575" target="_top">Haz clic aquí </a> </html>