<div class="gmail_quote">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></div><br>