Hi Luis,<br>I convert the Histogram into a 1DImage in this way:<br><br>typedef unsigned short PixelType1D;<br>typedef itk::GDCMImageIO ImageIOType;<br>typedef itk::Image< PixelType1D, 1 > ImageType1D;<br>typedef itk::Image< PixelType1D, 1 > ImageType1D_2;<br>
typedef itk::ImageFileReader< ImageType1D > ReaderType1D;<br>typedef itk::ImageRegionIterator< ImageType1D> IteratorType;<br><br>ReaderType::Pointer readerHis = ReaderType::New();<br>ImageIOType::Pointer gdcmImageIOHis= ImageIOType::New();<br>
HistogramGeneratorType::Pointer histogramGenerator = HistogramGeneratorType::New();<br>ReaderType1D::Pointer reader1D = ReaderType1D::New();<br>ImageType1D::Pointer line = ImageType1D::New();<br>ImageType1D::SizeType size1D;<br>
ImageType1D::IndexType pixelValue;<br>ImageType1D::IndexType start1D;<br>ImageType1D::RegionType region1D;<br><br><br><br>void createImage(int *valoresHisto){<br><br> size1D.Fill( 2286 );<br> start1D.Fill( 0 );<br> region1D.SetSize( size1D );<br>
region1D.SetIndex( start1D );<br> line->SetRegions( region1D );<br> line->Allocate();<br><br> IteratorType lineIt(line, line->GetRequestedRegion());<br><br> int i = 0;<br> for (lineIt.GoToBegin(); !lineIt.IsAtEnd();<br>
++lineIt)<br> {<br> lineIt.Set( valoresHisto[i]);<br><br> i=i+1;<br> }<br><br><br> typedef itk::ImageFileWriter< ImageType1D > WriterType1D;<br> WriterType1D::Pointer writer1D = WriterType1D::New();<br>
writer1D->SetFileName( "Image1D.png" );<br><br> writer1D->SetInput( line );<br> try<br> {<br> writer1D->Update();<br> }<br> catch (itk::ExceptionObject &e)<br> {<br> std::cerr << e << std::endl;<br>
<br> }<br>}<br><br>Now I want to apply SmoothingRecursiveGaussianImageFilter2.cxx, but I have a problem because if need a 2d Image. How can I resolve it?<br><br>Thanks<br><br>Raquel <br><br><br><div class="gmail_quote">
2009/7/4 Luis Ibanez <span dir="ltr"><<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Hi Raquel,<br>
<br>
You can convert the Histogram into a 1D Image,<br>
and use any of the ITK Smoothing filters on<br>
this 1D image, then convert the image back<br>
to a Histogram.<br>
<br>
See the class<br>
Insight/Code/Numerics/Statistics/itkHistogramToImageFilter.h<br>
<br>
and its related classes<br>
<br>
itkHistogramToEntropyImageFilter.h<br>
itkHistogramToIntensityImageFilter.h<br>
itkHistogramToLogProbabilityImageFilter.h<br>
itkHistogramToProbabilityImageFilter.h<br>
<br>
Then look at the section of Smoothing image filters<br>
in the ITK Software Guide<br>
<br>
<br>
<a href="http://www.itk.org/ItkSoftwareGuide.pdf" target="_blank">http://www.itk.org/ItkSoftwareGuide.pdf</a><br>
<br>
<br>
<br>
Regards,<br>
<br>
<br>
Luis<br>
<br>
<br>
<br>
--------------------<br>
Raquel Itk wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="h5">
Hi,<br>
<br>
I have computed the histogram of an image with Examples/Statistics/ImageHistogram2.cxx, now I want to smooth the array , where I have saved the values of an Histogram ( valoresHisto[i]) because I have a lot of peaks. I would like to use a Gaussian method, I have been read about Gaussian but only for Images.<br>
<br>
while( itr != end )<br>
<br>
{<br>
<br>
valoresHisto[i]=itr.GetFrequency();<br>
<br>
++itr;<br>
<br>
++binNumber;<br>
<br>
++i;<br>
<br>
}<br>
<br>
return(valoresHisto);<br>
<br>
<br>
Thanks<br>
<br>
<br>
Raquel<br>
<br>
<br></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>
</blockquote>
</blockquote></div><br>