<html><div style='background-color:'><P><BR><BR></P>
<DIV class=RTE>
<P>I forget send the code here it is (it = code)<BR></P><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. Consigue lo que buscas más rápido <a href="http://g.msn.com/8HMBESMX/2740??PS=47575" target="_top">Haz clic aquí. </a> </html>