Hello,<br><br>I'm using itk 2.8.1 and I have a problem reading some
DICOM images with itkGDCMImageIO. For most of images I opened since
now, it works ok, but I have images from several studies that is not able to open
because it crashes with a SIGFPE (floating point exception). When
debugging, this is the message after crashing<br>
<br><span style="font-family: courier new,monospace;">Program received signal SIGFPE, Arithmetic exception.</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">[Switching to Thread -1223084336 (LWP 6551)]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">0xb7de14e6 in gdcm::PixelReadConvert::ConvertReArrangeBits () from /usr/lib/InsightToolkit/libitkgdcm.so.2.8</span><br><br><br>and the backtrace is the following:<br><br>
<br>
<span style="font-family: courier new,monospace;">#0 0xb7de14e6 in gdcm::PixelReadConvert::ConvertReArrangeBits () from /usr/lib/InsightToolkit/libitkgdcm.so.2.8</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#1 0xb7de9c31 in gdcm::PixelReadConvert::ReadAndDecompressPixelData () from /usr/lib/InsightToolkit/libitkgdcm.so.2.8</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#2 0xb7dbd69c in gdcm::FileHelper::GetRaw () from /usr/lib/InsightToolkit/libitkgdcm.so.2.8</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#3 0xb7dbebf2 in gdcm::FileHelper::GetImageData () from /usr/lib/InsightToolkit/libitkgdcm.so.2.8</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#4 0xb7ee1892 in itk::GDCMImageIO::Read () from /usr/lib/InsightToolkit/libITKIO.so.2.8</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#5 0x08061964 in itk::ImageFileReader<itk::Image<int, 3u>, itk::DefaultConvertPixelTraits<int> >::GenerateData ()</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#6 0xb7a62380 in itk::ProcessObject::UpdateOutputData () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#7 0xb7a1f969 in itk::DataObject::UpdateOutputData () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#8 0xb7a1f770 in itk::DataObject::Update () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#9 0xb7a5ffca in itk::ProcessObject::UpdateLargestPossibleRegion () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#10 0x0805d6c0 in itk::ImageSeriesReader<itk::Image<int, 3u> >::GenerateData ()</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#11 0xb7a62380 in itk::ProcessObject::UpdateOutputData () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#12 0xb7a1f969 in itk::DataObject::UpdateOutputData () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#13 0xb7a1f770 in itk::DataObject::Update () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#14 0xb7a60062 in itk::ProcessObject::Update () from /usr/lib/InsightToolkit/libITKCommon.so.2.8</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#15 0x0804e0f4 in main ()</span><br style="font-family: courier new,monospace;">
<br>To test this problem I've writen a little program I've attached
(gdcmBug.cpp) reproducing the problem. You can use the following anonymized study ( <a href="http://hades.udg.edu/%7Echus/study1.tar.gz" target="_blank">http://hades.udg.edu/~chus/study1.tar.gz</a> ) with
some series to test this bug. It crashes with any of the images of the study.<br>
<br>I also tried to read the same images with other libraries such as
vtk or dcmtk and they're able to open the same images successfully.<br><br>Assuming ITKINCLUDEDIR is <span style="font-family: courier new,monospace;">/usr/include/InsightToolkit</span> and ITKLIBDIR is <font style="font-family: arial,sans-serif;" size="2"><font size="4"><font size="4"><tt>/usr/lib/InsightToolkit</tt></font>, </font>the compile command for g++ is the following</font><br>
<br><tt><font size="4"><b> g++ -o gdcmBug gdcmBug.o -I. -I/usr/include/InsightToolkit
-I/usr/include/InsightToolkit/IO -I/usr/include/InsightToolkit/gdcm/src
-L/usr/lib/InsightToolkit -lITKIO -litkgdcm</b></font><br></tt><br>To run this test you have to pass by argument the files you want to read. <br><br>For example, if you want to read 3 images (im1, im2, im3) from the same series:<br>
<br><span style="font-family: courier new,monospace;">$>gdcmBug im1 im2 im3</span><br><br>If you want to test with only one single image<br><br><span style="font-family: courier new,monospace;">$>gdcmbug im1 </span><br>
<br>and so on...<br><br>I tested this under Linux Mandriva 2008.0 and Kubuntu Feisty 7.04.<br><br>Maybe
this bug is resolved in a more recent version of itk, but I don't know
it. If so, let me know which version has solved this problem.<br>
<br>I hope this test and the images I provided can help to find the
problem. I have more images to reproduce the problem, if necessary.<br><br><br>Thanks in advance!<br>