<span style="color: rgb(64, 64, 255); font-weight: bold;">Hi all</span><br style="color: rgb(64, 64, 255); font-weight: bold;"> <span style="color: rgb(64, 64, 255); font-weight: bold;">Hi Yan</span><br style="color: rgb(64, 64, 255); font-weight: bold;"> <span style="color: rgb(64, 64, 255); font-weight: bold;">thx for reply</span><br style="color: rgb(64, 64, 255); font-weight: bold;"> <span style="color: rgb(64, 64, 255); font-weight: bold;">but the problem is that i need to browse each pixel of the image extracted, so i have to use iterators, and i did, and there the same problem occur.</span><br style="color: rgb(64, 64, 255); font-weight: bold;"> <span style="color: rgb(64, 64, 255); font-weight: bold;">i make a test after reading the pixel if it's value is equal to 0, if not i process it and so on,</span><br style="color: rgb(64, 64, 255); font-weight: bold;"> <span style="color: rgb(64, 64, 255); font-weight: bold;">but all the test are false, mean he find that all pixel
value are equal to 0( he didn't process any pixel) even many pixel are not equal to 0. if i use the classic method ( using GetPixel() method) , the program work well, but it is very slow!<br> this is the part of code where i use the iterators :</span><br> <span style="font-weight: bold;">......</span><br style="font-weight: bold;"> <span style="font-weight: bold;">ConstIteratorType2 out2( m_DiasstoledImage, m_diasstoledimageregion );</span><br style="font-weight: bold;"> <span style="font-weight: bold;">for ( out2.GoToBegin(); !out2.IsAtEnd(); ++out2)</span><br style="font-weight: bold;"> <span style="font-weight: bold;"> {</span><br style="font-weight: bold;"> <span style="font-weight: bold;"> if(out2.Get()!=0)</span><br style="font-weight: bold;"> <span style="font-weight: bold;"> {</span><br style="font-weight: bold;"> <span style="font-weight:
bold;"> </span><br style="font-weight: bold;"> <span style="font-weight: bold;"> m_diasstoledimageindex=out2.GetIndex();</span><br style="font-weight: bold;"> <span style="font-weight: bold;"> // Coord du pixel consid‚r‚ dans le plan image en millimetres.</span><br style="font-weight: bold;"> <span style="font-weight: bold;"> CoordInRecp[0]=m_diasstoledimageindex[0]*m_diasstoledimagespacing[0];</span><br style="font-weight: bold;"> <span style="font-weight: bold;"> CoordInRecp[1]=m_diasstoledimageindex[1]*m_diasstoledimagespacing[0];</span><br style="font-weight: bold;"><span style="font-weight: bold;">.......<br> <br> <br> <span style="color: rgb(64, 64, 255);">thx in advance</span><br style="color: rgb(64, 64, 255);"> <span
style="color: rgb(64, 64, 255);">any comment is welcome</span><br style="color: rgb(64, 64, 255);"> <span style="color: rgb(64, 64, 255);">amir,</span><br> <br> <br> </span><b><i>Pingkun Yan <pingkun.yan@gmail.com></i></b> a écrit :<blockquote class="replbq" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"> Hi charfeddine,<br><br>I am not sure if the problem is caused by<br>m_thresholdregionsize[2]=0;<br>Maybe you should set it to 1.<br><br>In fact, you have a better choice. You can use the ITK extract slice<br>filter, which will extract one slice from a volume. You can find<br>examples in ITK. Then you just need to use a 2D image iterator to read<br>this 2D slice, which will be much easier.<br><br>HTH,<br>Pingkun<br><br>On 2/9/06, charfeddine amir <charfeddine_amir @yahoo.fr=""> wrote:<br>> Hi all<br>> I'm reading Dicom file containing many slice, from which i extract one<br>> slice.<br>> when i try to show the contenent
of the image extracted (pixel value) using<br>> iterators i get a black screen, and when i try GetPixel method( which is too<br>> slow), i can show the value.<br>> did i make a mistake somewhere!<br>><br>> this is the code<br>><br>> /************************************************************************************/<br>> typedef itk::Image<float ,3=""> InputImageType;<br>> typedef itk::Image< float, 2 > DiasstolImageType;<br>><br>> typedef itk::ImageRegionConstIterator< InputImageType > ConstIteratorType;<br>> typedef itk::ImageRegionIterator< DiasstolImageType> IteratorType;<br>><br>> .....<br>> m_thresholdregion=m_OutputImage->GetLargestPossibleRegion();<br>> ThresholdedImageType::SizeType m_thresholdregionsize=<br>> m_thresholdregion.GetSize();<br>> m_thresholdregionsize[2]=0;<br>><br>> ThresholdedImageType::IndexType m_thresholdregionstart =<br>>
m_thresholdregion.GetIndex();<br>> m_thresholdregionstart[2]=slicenumber;<br>><br>> m_thresholddesiredregion.SetSize( m_thresholdregionsize );<br>> m_thresholddesiredregion.SetIndex( m_thresholdregionstart<br>> );<br>><br>><br>> std::cout<<std::endl =""><<"---------------------------------->incidence<br>> N°"<<card =""><<std::endl ;=""><br>><br>><br>> ConstIteratorType in( m_OutputImage, m_thresholddesiredregion );<br>> IteratorType out( m_DiasstoledImage, m_diasstoledimageregion );<br>> for ( in.GoToBegin(), out.GoToBegin(); !in.IsAtEnd(); ++in, ++out)<br>> {<br>> out.Set( in.Get() );<br>> std::cout<<in.Get (=""><<" - ";<br>> }<br>><br>> /*******************************************************************************/<br>> i get this screen on execution:<br>><br>> La diasstole de l'incidence courante est l'image 0<br>><br>>
---------------------------------->incidence N°0<br>><br>> **********************************Incidence<br>> suivante****************************<br>><br>> La diasstole de l'incidence courante est l'image 0<br>><br>> ---------------------------------->incidence N°1<br>><br>> **********************************Incidence<br>> suivante****************************<br>><br>> /*******************************************************************************************/<br>> But if remplace the use of iterators loop with this classic one :<br>><br>> pixelIndex[2]=slicenumber;<br>> for(j=0; j<m_inputimagesize [1="" j=""><br>> {<br>> pixelIndex[1]=j;<br>> m_diasstoledimageindex[1]=j;<br>> for(i=0; i<m_inputimagesize [0="" i=""><br>> {<br>> m_diasstoledimageindex[0]=i;<br>> pixelIndex[0]=i;<br>> pixelValue = m_OutputImage->GetPixel( pixelIndex
);<br>><br>> m_DiasstoledImage->SetPixel(m_diasstoledimageindex,pixelValue);<br>> std::cout<<pixelValue =""><<" - ";<br>> }<br>> }<br>><br>> I get the next screen:<br>><br>> La diasstole de l'incidence courante est l'image 0<br>> 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 -<br>> 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 -<br>> 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 -<br>> 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 -<br>> 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 -<br>> 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255<br>> - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 -<br>> 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255<br>> - 255 - 255 - 255 - 255 - 255 -
255 - 255 - 255 - 255 - 255 - 255 - 255 -<br>> 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255<br>> - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 - 255 -<br>> 255 - 255 - 255 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0<br>> - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0 - 0<br>> - 0 -.........<br>><br>> any comment is welcome<br>> thx<br>> amir<br>><br>><br>> ________________________________<br>> Nouveau : téléphonez moins cher avec Yahoo! Messenger ! Découvez les tarifs<br>> exceptionnels pour appeler la France et l'international. Téléchargez la<br>> version beta.<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></pixelValue></m_inputimagesize></m_inputimagesize></in.Get></std::endl></card></std::endl></float></charfeddine_amir></blockquote><br><p>
<hr size=1> Nouveau : téléphonez moins cher avec Yahoo! Messenger ! Découvez les tarifs exceptionnels pour appeler la France et l'international.
<a href="http://us.rd.yahoo.com/messenger/mail_taglines/default/*http://fr.beta.messenger.yahoo.com">Téléchargez</a> la version beta.