Hi Camille,<br><br>The problem is that you are creating an empty sample and <br>passing it as input to the filter method:<br><br> filter->SetInputSample(sample);<br><br>The SetInputSample() method is implemented in the base<br>
class: itk::SampleAlgorithmBase as<br><br> void SetInputSample( const TInputSample * sample ) <br> {<br> if ( m_InputSample != sample )<br> {<br> m_InputSample = sample;<br> m_MeasurementVectorSize = m_InputSample->GetMeasurementVectorSize();<br>
this->Modified();<br> }<br> }<br><br><br>In practice you should populate the sample with data <br>before you pass it to a filter.<br><br>See for example:<br><br>Insight/Testing/Code/Numerics/Statistics/<br>
itkSampleMeanShiftClusteringFilterTest.cxx<br><br><br>The class that you may want to use is the<br><br> ImageToListSampleAdaptor.<br><br><br><br><br> Regards,<br><br><br> Luis<br><br><br>
----------------------------------------------------------------------<br><div class="gmail_quote">On Thu, Jul 30, 2009 at 8:57 AM, camille orny <span dir="ltr"><<a href="mailto:camille.orny@gmail.com">camille.orny@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hello,<br><br>I am trying to use the itk::SampleMeanShiftBlurringFilter to get, for each pixel, the value of the mode the mean shift algorithm associates it to.<br>
<br>I have written the following code which is compiled without any error reports. Hower I get a "segmentation error" when I execute it. It seems the problem comes from the line filter->SetInputSample(sample); when I progressively comment the code.<br>
<br>Are there any problems with that filter ?<br><br>Regards,<br><br>typedef itk::Image< PixelType, 2> TInputImage;<br>typedef itk::Image< PixelType, 2> TOutputImage;<br>typedef itk::Image< unsigned char, 2> TOutputImage2;<br>
typedef otb::ImageFileReader<<div>TInputImage> ReaderType;<br>typedef itk::Statistics::ScalarImageToListAdaptor<TInputImage> SampleType;<br><br></div><div>typedef itk::RescaleIntensityImageFilter<TOutputImage,TOutputImage2> RescalerType;<br>
typedef otb::StreamingImageFileWriter<TOutputImage> WriterType;<br>typedef otb::StreamingImageFileWriter<TOutputImage2> WriterType2;<br></div>//typedef std::vector< int > ClusterLabelsType;<div>
<br>typedef itk::Statistics::SampleMeanShiftBlurringFilter<SampleType> FilterType;<br>
<br></div>//typedef itk::FixedArray<int, 1> TVector;<div><br>//typedef itk::Statistics::ListSample < TVector > ListSampleType;<br><br>ReaderType::Pointer reader=ReaderType::New();<br>SampleType::Pointer sample=SampleType::New();<br>
FilterType::Pointer filter=FilterType::New();<br>WriterType::Pointer writer=WriterType::New();<br>WriterType2::Pointer writer2=WriterType2::New();<br>//RescalerType::Pointer rescaler=RescalerType::New();<br><br>reader->SetFileName(argv[1]);<br>
reader->Update();<br>TInputImage::Pointer inimage=reader->GetOutput();<br><br>sample->SetImage(inimage);<br><br>filter->SetInputSample(sample);<br>filter->Update();<br><br>TOutputImage::Pointer limage=TOutputImage::New();<br>
limage->CopyInformation(inimage);<br>limage->SetRegions(limage->GetLargestPossibleRegion());<br>limage->Allocate();<br>limage->FillBuffer(0);<br><br><br>itk::ImageRegionIterator<TOutputImage> lit(limage,limage->GetLargestPossibleRegion());<br>
lit.GoToBegin();<br> <br><br></div>ListSampleType::Iterator it=(filter->GetOutput())-><div>Begin();<br>int i=0;<br><br>std::cout<<"hello";<br><br>while(!lit.IsAtEnd() )<br><br>
{<br> std::cout<<"while";<br>
lit.Set((filter->GetOutput())->GetMeasurementVector(i)[0]);<br> std::cout<<(filter->GetOutput())->GetMeasurementVectorSize();<br> ++lit;<br> ++it;<br> ++i;<br>}<br><br>writer->SetInput(inimage);<div>
<br>
writer->SetFileName(argv[2]);<br>writer->Update();</div></div><br>
<br>_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at: <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br></blockquote></div><br>