<div>Hi Lassi,</div>
<div> </div>
<div>I add</div>
<div> </div>
<div> try<br>    {<br>    reader-&gt;Update();<br>    }<br>  catch (itk::ExceptionObject &amp; e)<br>    {<br>    std::cerr &lt;&lt; &quot;exception in file reader &quot; &lt;&lt; std::endl;<br>    std::cerr &lt;&lt; e &lt;&lt; std::endl;<br>
    return EXIT_FAILURE;<br>    }</div>
<div> </div>
<div>after </div>
<div>reader-&gt;SetFileName( argv[1] );</div>
<div> </div>
<div>Now the results look reasonable.</div>
<div> </div>
<div>Thanks for your help.</div>
<div> </div>
<div>Howard<br><br></div>
<div class="gmail_quote">On Tue, Dec 22, 2009 at 6:15 AM, Lassi Paavolainen <span dir="ltr">&lt;<a href="mailto:lassi.paavolainen@jyu.fi" target="_blank">lassi.paavolainen@jyu.fi</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div>Hi Hsiang-Chi,<br><br>Maybe adding: image-&gt;Update; between these lines would work<br></div>
<div>   image = reader-&gt;GetOutput();<br>   IteratorType it( image, image-&gt;GetLargestPossibleRegion() );<br><br></div>
<div>I think you can also change<br></div>
<div>   InputImageType::Pointer image = InputImageType::New();<br></div>
<div>to<br>   InputImageType::Pointer image;<br><br>There is no need to create new image as you get it from reader.<br><br>Lassi<br><br>On Mon, 21 Dec 2009, Hsiang-Chi Kuo wrote:<br><br></div>
<div>
<div></div>
<div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi,<br><br>I am trying to access the maximum and minimum pixel value of a 3D image with<br>the following code, however, the result dose not look right (maximum 0,<br>
minimum 1000). Can anyone give me a hint why this code dose not work?<br><br>Thanks a lot,<br><br>Howard<br><br> <br><br>   typedef unsigned int InputPixelType;<br>  const unsigned int   Dimension = 3;<br> <br>  typedef unsigned int OutputPixelType;<br>
 <br>  typedef itk::Image&lt; OutputPixelType, Dimension &gt;   OutputImageType;<br><br>  typedef itk::Image&lt; InputPixelType, Dimension &gt; InputImageType;<br><br>  typedef itk::ImageFileReader&lt; InputImageType &gt; ReaderType;<br>
<br>  ReaderType::Pointer reader = ReaderType::New();<br>  reader-&gt;SetFileName( argv[1] );<br> <br>  typedef itk::ImageRegionIterator&lt; InputImageType &gt; IteratorType;<br>  InputImageType::Pointer image = InputImageType::New();<br>
  image = reader-&gt;GetOutput();<br>  IteratorType it( image, image-&gt;GetLargestPossibleRegion() );<br>  int max=0;<br>  int min=10000;<br> <br>  for (it.GoToBegin();!it.IsAtEnd();++it)<br>  {<br>  int  temp = it.Get();<br>
   if (temp&gt;max)<br>    max=temp;<br>   if (temp&lt;min)<br>    min=temp;<br>  }<br>  std::cout &lt;&lt; &quot;min=&quot; &lt;&lt; min &lt;&lt;std::endl;<br>  std::cout &lt;&lt; &quot;max=&quot; &lt;&lt; max &lt;&lt;std::endl;<br>
<br> <br><br><br><br></blockquote><br></div></div>
<div>-- <br>Lassi Paavolainen, M.Sc<br>Software Engineer, PhD Student in Computer Science<br>BioImageXD (<a href="http://www.bioimagexd.net/" target="_blank">http://www.bioimagexd.net</a>)<br></div>University of Jyväskylä<br>
<font color="#888888"><a href="mailto:lassi.paavolainen@jyu.fi" target="_blank">lassi.paavolainen@jyu.fi</a></font></blockquote></div><br>