Thanks very much to you all for looking into my "problem". Actually you guys were absolutely correct that there was no "problem" with the sample code. My "problem" was created when I tried to write the output real/imaginary part in a binary file. Sorry for the noise.<div>
<br></div><div>Ming<br><br><div class="gmail_quote">On Mon, Aug 9, 2010 at 3:31 AM, Yu, Yang <span dir="ltr"><<a href="mailto:yang.yu@philips.com">yang.yu@philips.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D">I copy the code you posted and at least it works fine in my program. I think you should pay much attention on the realFilter part. Did you Update the realFilter?</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D"> </span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt">From:</span></b><span style="font-size:10.0pt"> <a href="mailto:insight-users-bounces@itk.org" target="_blank">insight-users-bounces@itk.org</a> [mailto:<a href="mailto:insight-users-bounces@itk.org" target="_blank">insight-users-bounces@itk.org</a>]
<b>On Behalf Of </b>Ming Chao<br>
<b>Sent:</b> 2010 Aug 05 7:21 PM<br>
<b>To:</b> Markus Neuner<br>
<b>Cc:</b> Insight Users<br>
<b>Subject:</b> Re: [Insight-users] FFTImageFilter</span></p>
</div><div><div></div><div class="h5">
<p class="MsoNormal"> </p>
<p class="MsoNormal" style="margin-bottom:12.0pt">If I changed "WriteImageType" to float/double, the code stopped working b/c jpeg only accepts unsigned char/int. </p>
<div>
<p class="MsoNormal">On Thu, Aug 5, 2010 at 3:30 AM, Markus Neuner <<a href="mailto:neuner.markus@gmx.net" target="_blank">neuner.markus@gmx.net</a>> wrote:</p>
<div>
<p class="MsoNormal">I just looked at the code in the example and all seems working (at leasts for me).<br>
What do you expect as output?<br>
Change the "WriteImageType" to float or double and try it again.<br>
Open the image with paraview, imagej or any other suited image viewer.<br>
It may be useful to take the log of the real and complex images for a better contrast.</p>
<div>
<div>
<p class="MsoNormal"><br>
<br>
On 04.08.2010 21:45, Ming Chao wrote: </p>
<p class="MsoNormal" style="margin-bottom:12.0pt">float</p>
<div>
<p class="MsoNormal">On Wed, Aug 4, 2010 at 2:38 PM, Markus Neuner <<a href="mailto:neuner.markus@gmx.net" target="_blank">neuner.markus@gmx.net</a>> wrote:</p>
<div>
<p class="MsoNormal">What is the pixel type of "ImageType"? </p>
<div>
<div>
<p class="MsoNormal"><br>
<br>
On 04.08.2010 21:11, Ming Chao wrote: </p>
</div>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal">I checked the complex image FFT created from the input image was fine. Also I checked the real and imaginary parts in floating point format, respectively. They both are fine. The problem started to happen when rescaling the real/imaginary
image into an interger-type image. Here is part of the rescaling code: </p>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<div>
<p class="MsoNormal"> typedef unsigned char WritePixelType;</p>
</div>
<div>
<p class="MsoNormal"> typedef itk::Image< WritePixelType, Dimension > WriteImageType;</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> typedef itk::RescaleIntensityImageFilter<ImageType, WriteImageType > RescaleFilterType;</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> RescaleFilterType::Pointer intensityRescaler = RescaleFilterType::New();</p>
</div>
<div>
<p class="MsoNormal"> intensityRescaler->SetOutputMinimum( 0 );</p>
</div>
<div>
<p class="MsoNormal"> intensityRescaler->SetOutputMaximum( 255 );</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> intensityRescaler->SetInput( realFilter->GetOutput() );</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> typedef itk::ImageFileWriter< WriteImageType > WriterType;</p>
</div>
<div>
<p class="MsoNormal"> WriterType::Pointer writer = WriterType::New();</p>
</div>
<div>
<p class="MsoNormal"> writer->SetFileName( argv[2] );</p>
</div>
<div>
<p class="MsoNormal"> writer->SetInput( intensityRescaler->GetOutput() );</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> try</p>
</div>
<div>
<p class="MsoNormal"> {</p>
</div>
<div>
<p class="MsoNormal">writer->Update();</p>
</div>
<div>
<p class="MsoNormal"> }</p>
</div>
<div>
<p class="MsoNormal"> catch( itk::ExceptionObject & excp )</p>
</div>
<div>
<p class="MsoNormal"> {</p>
</div>
<div>
<p class="MsoNormal">std::cerr << "Error writing the real image: " << std::endl;</p>
</div>
<div>
<p class="MsoNormal">std::cerr << excp << std::endl;</p>
</div>
<div>
<p class="MsoNormal">return EXIT_FAILURE;</p>
</div>
<div>
<p class="MsoNormal"> }</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> typedef itk::ImageRegionConstIterator< WriteImageType > ConstIteratorType;</p>
</div>
<div>
<p class="MsoNormal"> WriteImageType::Pointer RealImage = intensityRescaler->GetOutput();</p>
</div>
<div>
<p class="MsoNormal"> ConstIteratorType in( RealImage, RealImage->GetRequestedRegion() );</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> for ( in.GoToBegin(); !in.IsAtEnd(); ++in )</p>
</div>
<div>
<p class="MsoNormal"> {</p>
</div>
<div>
<p class="MsoNormal"> std::cout << in.Get() << std::endl;</p>
</div>
<div>
<p class="MsoNormal"> }</p>
</div>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">From the output I got all zeros....Apparently something wrong with the rescaler. Could anybody take a look? Thanks very much,</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">Ming</p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"> </p>
<div>
<p class="MsoNormal">On Wed, Aug 4, 2010 at 9:35 AM, Ming Chao <<a href="mailto:mingchao2005@gmail.com" target="_blank">mingchao2005@gmail.com</a>> wrote:</p>
<p class="MsoNormal">I tried to run FFTImageFilter.cxx to do an FFT of an image. The input image is a 512x512, 24 bit thoracic CT image in jpg format. The output of the example cxx code are the real and imaginary parts of FFT of the input image. However, when
I look at the two output images, they are just uniformly grey/black, depending on
</p>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<div>
<p class="MsoNormal"> intensityRescaler->SetOutputMinimum( 0 );</p>
</div>
<div>
<p class="MsoNormal"> intensityRescaler->SetOutputMaximum( 255 );</p>
</div>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">in the code. What did I do wrong? Many thanks,</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888">Ming</span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"> </span></p>
</div>
</div>
<p class="MsoNormal"> </p>
</div>
</div>
</div>
<pre> </pre>
<pre>_____________________________________</pre>
<pre>Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a></pre>
<pre> </pre>
<pre>Visit other Kitware open-source projects at</pre>
<pre><a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a></pre>
<pre> </pre>
<pre>Kitware offers ITK Training Courses, for more information visit:</pre>
<pre><a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a></pre>
<pre> </pre>
<pre>Please keep messages on-topic and check the ITK FAQ at:</pre>
<pre><a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a></pre>
<pre> </pre>
<pre>Follow this link to subscribe/unsubscribe:</pre>
<pre><a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a></pre>
<pre> </pre>
</blockquote>
</div>
</div>
<p class="MsoNormal"> </p>
</div>
</div>
</div>
</div>
<p class="MsoNormal"> </p>
</div></div></div>
<br>
<hr>
<font face="Arial" color="Gray" size="1">The information contained in this message may be confidential and legally protected under applicable law. The message is intended solely for the addressee(s). If you are not the intended recipient, you are hereby notified
that any use, forwarding, dissemination, or reproduction of this message is strictly prohibited and may be unlawful. If you are not the intended recipient, please contact the sender by return e-mail and destroy all copies of the original message.<br>
</font>
</div>
</blockquote></div><br></div>