<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Hello, <br><br>Thank you for your reply again. <br>I am using ITK-3.20.0 version. <br>The source code is as follows, in it, the "ABC Dicom Loader" is a VTK based data loader. Hope it could be helpful to figure out its error.: <br>******************************************************************************************<br>#if defined(_MSC_VER)<br>#pragma warning ( disable : 4786 )<br>#endif<br><br>#ifdef __BORLANDC__<br>#define ITK_LEAN_AND_MEAN<br>#endif<br><br>//for DCMTK<br>#define HAVE_CONFIG_H 1 <br>#include <osconfig.h><br><br>#include <iostream><br>#include <fstream><br>#include <sstream><br>#include <cstring><br>#include <string><br>#include <stdio.h><br>#include <stdlib.h><br><br>#include "BaseData.h"<br><br>//for casting to float and converting to ITK image<br>#include "vtkImageCast.h"<br>#include
"./Headerfiles/itkVTKImageToImageFilter.h"<br>#include "itkCastImageFilter.h"<br><br><br><br>//for geodesic active contor segmentation<br>#include "itkGradientAnisotropicDiffusionImageFilter.h"<br>#include "itkScalarImageToHistogramGenerator.h"<br><br>//for itk images<br>#include "itkImage.h"<br>#include "itkImageRegionIterator.h"<br>#include "itkImageRegionConstIterator.h"<br>#include "itkCovariantVector.h"<br>#include "itkImageRegionConstIterator.h"<br>#include "itkImageRegionIterator.h"<br><br>//for itk vector<br>#include "itkVector.h"<br><br>//modified ITK classes<br>#include "./Headerfiles/VavGradientVectorFlowImageFilter.h"<br>#include "./Headerfiles/VavEdgeEnhancingDiffusionImageFilter.h"<br><br>//header files for reading/writing dicom series<br>#include "itkImage.h"<br>#include "itkGDCMImageIO.h"<br>#include "itkOrientedImage.h"<br>#include "itkImageFileReader.h"<br>#include "itkImageFileWriter.h"<br>#include "itkImageSeriesReader.h"<br>#include
"itkGDCMSeriesFileNames.h"<br><br>//header files for Demons algorithm registration<br>#include "itkCastImageFilter.h"<br>#include "itkWarpImageFilter.h"<br>#include "itkHistogramMatchingImageFilter.h"<br>#include "itkLinearInterpolateImageFunction.h"<br>#include "itkSymmetricForcesDemonsRegistrationFilter.h"<br>#include "itkDemonsRegistrationFilter.h" <br><br>//for reading images<br>#include <QtGui><br>#include "DICOMLoader.h"<br>//#include "dicomviewer2d.h"<br>#include "ContentTree.h"<br>#include "ContentItem.h"<br>#include <QtUiTools><br>#include "/usr/local/include/vtk-5.2/QVTKWidget.h"<br>#include <QMainWindow><br>#include <QVBoxLayout><br>#include <QSlider><br>#include <QtGui><br><br>//for converting itk result back to vtk<br>#include "itkImageToVTKImageFilter.h"<br><br>//for interpolator <br>#include "itkContinuousIndex.h"<br>#include "itkLinearInterpolateImageFunction.h"<br><br>// image types<br>const unsigned
int Dimension = 3;<br>typedef unsigned short PixelType;<br>typedef float InternalPixelType;<br>typedef itk::Image< InternalPixelType, Dimension > InternalImageType;<br>typedef itk::Image< float, 3 > ImageType;<br>typedef itk::Image< float, 3 > ShapeImageType;<br>typedef itk::Image< int, 3 > IntegerImageType;<br>typedef itk::VTKImageToImageFilter<ImageType> VTKImageToImageType;<br>typedef itk::Image< unsigned char, 3 > OutputImageType;// 2D output image for .png writers<br>typedef itk::CovariantVector<float, 3> GradientPixelType;<br>typedef itk::Image<GradientPixelType,3> GradientImageType;<br>typedef itk::Image< PixelType, Dimension > FixedImageType;<br>typedef itk::Image< PixelType, Dimension > MovingImageType;<br>// region iterators<br>typedef itk::ImageRegionConstIterator<ImageType> ImageConstIterator;<br>typedef
itk::ImageRegionConstIterator<OutputImageType> ImageConstIteratorChar;<br>typedef itk::ImageRegionIterator<ImageType> ImageIterator;<br>typedef itk::ImageRegionIterator<GradientImageType> DerivativeImageIterator;<br><br><br>//define filters<br>typedef itk::VTKImageToImageFilter < InternalImageType > VTK2ITKConnectorFilterType;<br>//matchingfilter type<br>typedef itk::HistogramMatchingImageFilter<InternalImageType, InternalImageType > MatchingFilterType;<br><br>//deformation field type<br>typedef itk::Vector< float, Dimension > VectorPixelType;<br>typedef itk::Image< VectorPixelType, Dimension > DeformationFieldType;<br><br>// writers<br>typedef itk::ImageFileWriter< OutputImageType > WriterType;<br>typedef itk::ImageFileWriter< ImageType > WriterTypeDouble;<br>typedef itk::ImageFileWriter< GradientImageType >
WriterTypeGradient;<br>typedef itk::ImageFileWriter< DeformationFieldType > FieldWriterType;<br><br>// reader<br>typedef itk::ImageFileReader< ImageType > ReaderType;<br>typedef itk::ImageFileReader< GradientImageType > ReaderTypeGradient;<br><br><br>//registration filter type<br> typedef itk::SymmetricForcesDemonsRegistrationFilter<br> <InternalImageType,InternalImageType,DeformationFieldType> RegistrationFilterType;<br><br>//warper type<br> typedef itk::WarpImageFilter<<br> ImageType, <br>
ImageType,<br> DeformationFieldType> WarperType;<br><br>//interpolator type<br> typedef itk::LinearInterpolateImageFunction<<br> ImageType, double> InterpolatorType;<br><br>//output image type<br> typedef unsigned char OutputPixelType;<br> typedef itk::Image< OutputPixelType, Dimension > OutputImageType;<br><br>//castfilter type<br> typedef itk::CastImageFilter< <br>
ImageType,<br> OutputImageType > CastFilterType;<br><br>//itk to vtk connector<br>typedef itk::ImageToVTKImageFilter< ImageType > ConnectorType;<br><br><br><br><br>using namespace std;<br><br>//implements a command observer<br>class CommandIterationUpdate : public itk::Command <br>{<br> public:<br> typedef CommandIterationUpdate Self;<br> typedef itk::Command Superclass;<br> typedef itk::SmartPointer<CommandIterationUpdate> Pointer;<br> itkNewMacro( CommandIterationUpdate );<br> protected:<br> CommandIterationUpdate() {};<br><br>// typedef itk::Image< float, 3 > InternalImageType;<br><br> typedef itk::Image< InternalPixelType, Dimension >
InternalImageType;<br> typedef itk::Vector< float, 3 > VectorPixelType;<br> typedef itk::Image< VectorPixelType, 3 > DeformationFieldType;<br><br><br> public:<br><br> void Execute(itk::Object *caller, const itk::EventObject & event)<br> {<br> Execute( (const itk::Object *)caller, event);<br> }<br><br>void Execute(const itk::Object * object, const itk::EventObject & event)<br> {<br> const RegistrationFilterType * filter = <br> dynamic_cast< const RegistrationFilterType * >( object );<br> if( !(itk::IterationEvent().CheckEvent( &event )) )<br> {<br> return;<br>
}<br> cout << filter->GetMetric() <<endl;<br> }<br> };<br><br><br>int main( int argc, char * argv[] )<br>{<br><br><br> //////////////////////////////////////////////<br> //////////// READING IMAGES /////////////////<br> ////////////////////////////////////////////<br><br><br>/////////////////////////////////////<br>/// Reading the fixed image ////////<br>///////////////////////////////////<br><br> ABC::Loader::DICOMLoader* fixedImageLoader = new ABC::Loader::DICOMLoader();<br> ABC::Data::BaseData* fixedImageData= new ABC::Data::BaseData();<br> ContentTree *fixedImageTree;<br> ContentItem* fixedImageItem;<br><br><br><br> fixedImageTree=NULL;<br><br> fixedImageTree =
fixedImageLoader->LoadDataContent(argv[1]);<br> <br><br><br> fixedImageItem = fixedImageTree->GetFistItem()->GetFirstItem(); // first study first series<br> cout << " 1.2" << fixedImageItem->GetDescription() << endl;<br> cout << " 1.2" << fixedImageItem->GetID() << endl;<br> cout << " 1.2" << fixedImageItem->GetFSPath() << endl;<br><br> /// LOAD DATA !!!<br><br> fixedImageData = fixedImageLoader->LoadData(fixedImageItem);<br><br> vector<int> dimensionsFixed = fixedImageData->GetScalar()->GetDimensions();<br> printf("Dimensions %d %d %d \n", dimensionsFixed[0], dimensionsFixed[1], dimensionsFixed[2]);<br><br> vector<double> SpacingFixed =
fixedImageData->GetScalar()->GetSpacing();<br> printf("Spacing %f %f %f \n", SpacingFixed[0], SpacingFixed[1], SpacingFixed[2]);<br>cout<<"The end! "<<endl<<endl<<endl;<br> double* bndFixed = fixedImageData->GetScalar()->GetImageData()->GetBounds();<br> printf("Bounds1 %lf %lf %lf %lf %lf %lf \n", bndFixed[0], bndFixed[1], bndFixed[2], bndFixed[3], bndFixed[4], bndFixed[5]);<br><br><br><br> vtkImageCast* VtkCasterFixed = vtkImageCast::New(); // Casting to float in VTK is needed<br> VtkCasterFixed->SetInput(fixedImageData -> GetScalar() -> GetImageData());<br> VtkCasterFixed->SetOutputScalarTypeToFloat();<br><br> VTK2ITKConnectorFilterType::Pointer VTK2ITKconnectorFixed = VTKImageToImageType::New();<br> VTK2ITKconnectorFixed->SetInput( VtkCasterFixed->GetOutput());<br>
VTK2ITKconnectorFixed->GetImporter()->Update();<br><br><br><br>/////////////////////////////////////<br>/// Reading the moving image ////////<br>///////////////////////////////////<br><br> ABC::Loader::DICOMLoader* movingImageLoader = new ABC::Loader::DICOMLoader();<br> ABC::Data::BaseData* movingImageData= new ABC::Data::BaseData();<br> ContentTree *movingImageTree;<br> ContentItem* movingImageItem;<br><br> movingImageTree = movingImageLoader->LoadDataContent(argv[2]);<br><br> movingImageItem = movingImageTree->GetFistItem()->GetFirstItem(); // first study first series<br> cout << " 1.2" << movingImageItem->GetDescription() << endl;<br> cout << " 1.2" << movingImageItem->GetID() << endl;<br> cout << "
1.2" << movingImageItem->GetFSPath() << endl;<br><br> /// LOAD DATA !!!<br><br> movingImageData = movingImageLoader->LoadData(fixedImageItem);<br><br> vector<int> dimensionsMoving = movingImageData->GetScalar()->GetDimensions();<br> printf("Dimensions %d %d %d \n", dimensionsMoving[0], dimensionsMoving[1], dimensionsMoving[2]);<br><br> vector<double> SpacingMoving = movingImageData->GetScalar()->GetSpacing();<br> printf("Spacing %f %f %f \n", SpacingMoving[0], SpacingMoving[1], SpacingMoving[2]);<br><br> double* bndMoving = movingImageData->GetScalar()->GetImageData()->GetBounds();<br> printf("Bounds1 %lf %lf %lf %lf %lf %lf \n", bndMoving[0], bndMoving[1], bndMoving[2], bndMoving[3], bndMoving[4], bndMoving[5]);<br><br>//<br><br> vtkImageCast*
VtkCasterMoving = vtkImageCast::New(); // Casting to float in VTK is needed<br> VtkCasterMoving->SetInput(movingImageData -> GetScalar() -> GetImageData());<br> VtkCasterMoving->SetOutputScalarTypeToFloat();<br> VTK2ITKConnectorFilterType::Pointer VTK2ITKconnectorMoving = VTKImageToImageType::New();<br> VTK2ITKconnectorMoving->SetInput( VtkCasterMoving->GetOutput());<br> VTK2ITKconnectorMoving->GetImporter()->Update();<br><br><br>/////////////////////////////////////////////<br>// Demons registration starts from here ////<br>///////////////////////////////////////////<br><br>// The basic idea is to match the histograms of the two images at a user-specified number of quantile values. <br> MatchingFilterType::Pointer matcher = MatchingFilterType::New();<br> matcher->SetInput( VTK2ITKconnectorMoving->GetOutput() );<br>
matcher->SetReferenceImage( VTK2ITKconnectorFixed->GetOutput() );<br><br><br><br><br><br>//select the number of bins to represent the histograms and the<br>// number of points or quantile values where the histogram is to be<br>// matched.<br> matcher->SetNumberOfHistogramLevels( 1024 );<br> matcher->SetNumberOfMatchPoints( 7 );<br><br>// Simple background extraction is done by thresholding at the mean intensity.<br> matcher->ThresholdAtMeanIntensityOn();<br><br>//the deformation field is represented as an image whose pixels are floating point vectors.<br> RegistrationFilterType::Pointer filter = RegistrationFilterType::New();<br><br>//Create the Command observer and register it with the registration filter.<br> CommandIterationUpdate::Pointer observer = CommandIterationUpdate::New();<br> filter->AddObserver( itk::IterationEvent(), observer );<br> filter->SetFixedImage(
VTK2ITKconnectorFixed->GetOutput() );<br> filter->SetMovingImage( matcher->GetOutput() );<br><br><br> filter->SetNumberOfIterations( 1.0 );<br> filter->SetStandardDeviations( 1.0 );<br><br><br> filter->Update(); //the filter output is the computed deformation field <br> <br><br>//warp the moving image with the output deformation field <br> WarperType::Pointer warper = WarperType::New();<br> InterpolatorType::Pointer interpolator = InterpolatorType::New();<br><br> typedef itk::CastImageFilter< ImageType, InternalImageType > FixedImageCasterType;<br> typedef itk::CastImageFilter< ImageType, InternalImageType > MovingImageCasterType;<br><br><br> FixedImageCasterType::Pointer fixedImageCaster = FixedImageCasterType::New();<br> MovingImageCasterType::Pointer movingImageCaster = MovingImageCasterType::New();<br><br>
fixedImageCaster->SetInput( VTK2ITKconnectorFixed->GetOutput() );<br> movingImageCaster->SetInput( VTK2ITKconnectorMoving->GetOutput() ); <br><br> ImageType::Pointer fixedImage = fixedImageCaster->GetOutput();<br> ImageType::Pointer movingImage = movingImageCaster->GetOutput();<br><br> warper->SetInput( movingImage );<br> warper->SetInterpolator( interpolator );<br> warper->SetOutputSpacing( fixedImage->GetSpacing() );<br> warper->SetOutputOrigin( fixedImage->GetOrigin() );<br> warper->SetOutputDirection( fixedImage->GetDirection() );<br> warper->SetDeformationField( filter->GetOutput() );<br><br>// Write warped image out to file<br> WriterType::Pointer writer = WriterType::New();<br> CastFilterType::Pointer caster = CastFilterType::New();<br><br> writer->SetFileName( argv[3]
);<br> <br> caster->SetInput( warper->GetOutput() );<br> writer->SetInput( caster->GetOutput() );<br> writer->Update();<br><br><br><br>if( argc > 4 ) <br> {<br><br> FieldWriterType::Pointer fieldWriter = FieldWriterType::New();<br> fieldWriter->SetFileName( argv[6] );<br><br> fieldWriter->SetInput( filter->GetOutput() );<br> fieldWriter->Update();<br> }<br><br><br>return EXIT_SUCCESS;<br><br> }<br><br><br><br><br><br>--- On <b>Fri, 3/25/11, Ariel Hernán Curiale <i><curiale@gmail.com></i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Ariel Hernán Curiale <curiale@gmail.com><br>Subject: Re: [Insight-users] Itk Segmentation Fault<br>To: "akbarjan bayiz" <abayiz@yahoo.com><br>Cc: insight-users@itk.org<br>Date: Friday, March 25, 2011,
7:21 AM<br><br><div id="yiv707525447">What ITK version are running ?<div><br></div><div>Another things is that "DeformationPixelType" need to be a float or double.</div><div><br></div><div>Please post the source code.</div><div><br></div><div><br><div>
<span class="yiv707525447Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; font-size: medium;"><div>Saludos</div><div>Ariel Hernán Curiale</div><div><a rel="nofollow" target="_blank" href="http://www.curiale.com.ar">www.curiale.com.ar</a></div><div><br></div></span><br class="yiv707525447Apple-interchange-newline">
</div>
<br><div><div>El 25/03/2011, a las 12:54, akbarjan bayiz escribió:</div><br class="yiv707525447Apple-interchange-newline"><blockquote type="cite"><table style="" border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font: inherit;" valign="top">Hello!<br><br>Thank you for your reply. <br>I've done it as you said, changed the data types, but unfortunately it gives out a full screen of error now. I am trying to read two dicom series from two different directories, and perform Demons registration, then save the output as a new dicom file. Maybe this is different from what you are doing, and that is why your data type does not work for my code I think. <br><br>Do you have any other good suggestions? <br><br>Best Regards, <br><br>--- On <b>Fri, 3/25/11, Ariel Hernán Curiale <i><<a rel="nofollow" ymailto="mailto:curiale@gmail.com" target="_blank" href="/mc/compose?to=curiale@gmail.com">curiale@gmail.com</a>></i></b> wrote:<br><blockquote
style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Ariel Hernán Curiale <<a rel="nofollow" ymailto="mailto:curiale@gmail.com" target="_blank" href="/mc/compose?to=curiale@gmail.com">curiale@gmail.com</a>><br>Subject: Re: [Insight-users] Itk Segmentation Fault<br>To: "akbarjan bayiz" <<a rel="nofollow" ymailto="mailto:abayiz@yahoo.com" target="_blank" href="/mc/compose?to=abayiz@yahoo.com">abayiz@yahoo.com</a>><br>Cc:
<a rel="nofollow" ymailto="mailto:insight-users@itk.org" target="_blank" href="/mc/compose?to=insight-users@itk.org">insight-users@itk.org</a><br>Date: Friday, March 25, 2011, 5:34 AM<br><br><div id="yiv707525447">Hi, Im recently used this class and I had no errors. This is my code:<div><br></div><div><span class="yiv707525447Apple-tab-span" style="white-space: pre;"> </span> <font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">typedef</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size:
11px;"><font class="yiv707525447Apple-style-span" color="#010101">itk</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">::</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">Image</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101"><</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">PixelType</font></span></font><font
class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">,</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">3</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">> </font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">ImageType</font></span></font><font class="yiv707525447Apple-style-span" face="Menlo" size="3"><span
class="yiv707525447Apple-style-span" style="font-size: 11px;"><font class="yiv707525447Apple-style-span" color="#010101">;</font></span></font></div><div><div style="margin: 0px; font: 11px Menlo; color: rgb(0, 179, 249);"><span style="color: rgb(227, 57, 161);"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"> </span><font class="yiv707525447Apple-style-span" color="#010101">typedef</font></span><font class="yiv707525447Apple-style-span" color="#010101"> itk::Image< itk::Vector<DeformationPixelType, 3>, 3> DeformationFieldType</font><span style="color: rgb(255, 255, 255);">;</span></div></div><div><div style="margin: 0px; font: 11px Menlo; color: rgb(0, 179, 249);"><font class="yiv707525447Apple-style-span" color="#ffffff"><div style="margin: 0px; font: 11px Menlo;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font
class="yiv707525447Apple-style-span" color="#010101">typedef
itk::DiffeomorphicDemonsRegistrationFilter <ImageType, ImageType, DeformationFieldType> FilterType;</font></div><div style="margin: 0px; font: 11px Menlo; min-height: 13px;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font class="yiv707525447Apple-style-span" color="#010101"><br class="yiv707525447webkit-block-placeholder"></font></div><div style="margin: 0px; font: 11px Menlo;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font class="yiv707525447Apple-style-span" color="#010101">typename FilterType::Pointer filter = FilterType::New();</font></div><div style="margin: 0px; font: 11px Menlo; min-height: 13px;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font
class="yiv707525447Apple-style-span" color="#010101"><br class="yiv707525447webkit-block-placeholder"></font></div><div style="margin: 0px; font: 11px Menlo; min-height: 13px;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font class="yiv707525447Apple-style-span" color="#010101"><br class="yiv707525447webkit-block-placeholder"></font></div><div style="margin: 0px; font: 11px Menlo;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font class="yiv707525447Apple-style-span" color="#010101">filter->SetFixedImage(img_diastole->GetOutput());</font></div><div style="margin: 0px; font: 11px Menlo;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font
class="yiv707525447Apple-style-span" color="#010101">filter->SetMovingImage(img_sistole->GetOutput());</font></div><div style="margin: 0px; font: 11px Menlo;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font class="yiv707525447Apple-style-span" color="#010101">filter->Update();</font></div><div style="margin: 0px; font: 11px Menlo; min-height: 13px;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font class="yiv707525447Apple-style-span" color="#010101"><br class="yiv707525447webkit-block-placeholder"></font></div><div style="margin: 0px; font: 11px Menlo;"><span class="yiv707525447Apple-tab-span" style="white-space: pre;"><font class="yiv707525447Apple-style-span" color="#010101"> </font></span><font class="yiv707525447Apple-style-span"
color="#010101">typename DeformationFieldType::Pointer deformacion = filter->GetDeformationField();</font></div><div><span style="color: rgb(255, 255, 255);"><br></span></div></font></div></div><div><br><div>
<div>Saludos</div><div>Ariel Hernán Curiale</div><div><a rel="nofollow" target="_blank" href="http://www.curiale.com.ar/">www.curiale.com.ar</a></div><div><br></div><br class="yiv707525447Apple-interchange-newline">
</div>
<br><div><div>El 25/03/2011, a las 10:00, akbarjan bayiz escribió:</div><br class="yiv707525447Apple-interchange-newline"><blockquote type="cite"><div>Hello all!<br><br>I am using Demons algorithm to perform registration of two DICOM series. My code compiles without any error. And when I run it, at the final step, when it updates the filter, there is a "0" on the screen, and gives error message saying "Segmentation Fault". <br><br>What should I do for it? <br><br>Thank you in advance! <br><br><br><br>_____________________________________<br>Powered by <a rel="nofollow" target="_blank" href="http://www.kitware.com/">www.kitware.com</a><br><br>Visit other Kitware open-source projects at<br><a rel="nofollow" target="_blank" href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a><br><br>Kitware offers ITK Training Courses, for more information
visit:<br><a rel="nofollow" target="_blank" href="http://www.kitware.com/products/protraining.html">http://www.kitware.com/products/protraining.html</a><br><br>Please keep messages on-topic and check the ITK FAQ at:<br>http://www.itk.org/Wiki/ITK_FAQ<br><br>Follow this link to subscribe/unsubscribe:<br>http://www.itk.org/mailman/listinfo/insight-users<br></div></blockquote></div><br></div></div></blockquote></td></tr></tbody></table><br>
</blockquote></div><br></div></div></blockquote></td></tr></table><br>