Hi Marco,<br>The correct way of traversing elements of a image in ITK is through iterators:<br><br>typedef itk::ImageRegionIterator<TITKImgType> IteratorType;<br>IteratorType it(image,image->GetRequestedRegion());<br>
<br>for (it.GoToBegin(); !it.IsAtEnd(); ++it) {<br> it.SetPixel(it.GetPixel()+100);<br>}<br><br>This is an example and I don't know if it will compile or not, but you get the idea on how to use an iterator to traverse all image pixels. Not only the code is much simpler but also faster.<br>
<br>For more info you can look at the itkSoftwareGuide, that has a part dealing with iterators<br><br>Regards<br><br><div class="gmail_quote">On Mon, May 17, 2010 at 4:00 PM, Marco Gheza <span dir="ltr"><<a href="mailto:marcogheza4mailinglists@gmail.com">marcogheza4mailinglists@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi,<div>i'm trying to access every single voxel of a 3d image and change some values. In particular, i want to change that values that are not black; every value has to be increased by a value of 100 for example. If a voxel has value of 10, i want to give it value of 110.</div>
<div><br></div><div>Here is my code:</div><div><br></div><div>imageCT->DisconnectPipeline();</div><div><div>for(i=0;i<512;i++)</div><div> {</div><div><span style="white-space: pre;"> </span>for(j=0;j<512;j++)</div>
<div><span style="white-space: pre;"> </span>{</div><div><span style="white-space: pre;"> </span>for(k=0;k<311;k++)</div><div><span style="white-space: pre;"> </span>{</div>
<div><span style="white-space: pre;"> </span>pixelIndex[0] = i; // x position</div><div><span style="white-space: pre;"> </span>pixelIndex[1] = j; // y position</div><div>
<span style="white-space: pre;"> </span>pixelIndex[2] = k; // z position</div><div><br></div><div><span style="white-space: pre;"> </span>ImageType::PixelType pixelValuePT = imagePT->GetPixel( pixelIndex );</div>
<div><br></div><div><span style="white-space: pre;"> </span>ImageType::PixelType newValue = pixelValueCT+100;</div><div><br></div><div><span style="white-space: pre;"> </span>if(pixelValueCT!=0)</div>
<div><span style="white-space: pre;"> </span>{<span style="white-space: pre;"> </span></div><div><span style="white-space: pre;"> </span>imageCT->SetPixel( pixelIndex, newValue );</div>
<div><span style="white-space: pre;"> </span>}<span style="white-space: pre;"> </span></div><div><span style="white-space: pre;"> </span>}</div><div><span style="white-space: pre;"> </span>}</div>
<div> }</div><div><br></div><div>This code doesn't do the work well. Do you know how to change it?</div><div>Thank you,</div><div>bye,</div><div><br></div><font color="#888888"><div>Marco</div></font></div>
<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>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<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><br clear="all"><br>-- <br>Sergio Vera<br><br> Alma IT Systems<br> C/ Vilana, 4B, 4º 1ª<br> 08022 Barcelona<br> T. (+34) 932 380 592<br> <a href="http://www.alma3d.com">www.alma3d.com</a><br>