[vtkusers] vtkVolume ( vtkVolumeRayCastMapper ( vtkDICOMImageReader, ... ) ) Seg Faults

[Stewart Dickson]

In Vtk version 4.4.0

vtkVolume (
    vtkVolumeRayCastMapper (  vtkDICOMImageReader, 
vtkVolumeRayCastCompositeFunction),
    vtkVolumeProperty ( vtkPiecewiseFunction, vtkColorTransferFunction, 
vtkPiecewiseFunction)
    )

Produces a Segmentation fault as in the stack trace listed below.

I am using vtkDICOMImageReader SetDirectoryName 
"DICOM/PAT00000/ST000000/SE000002"
where DICOM/PAT00000/ST000000/SE000002 contains  MR0000[00-14]: DICOM 
medical imaging data files.

Attached script is derived from 
VolumeRendering/Testing/Tcl/volVPGradientOpacity.tcl
and Examples/VolumeRendering/Tcl/SimpleRayCast.tcl

At: http://www.vtk.org/doc/nightly/html/classvtkDICOMImageReader.html
it is written:
> ...the reader ... will select the subset corresponding to the first 
> series UID it stumbles across and it will try to build an ordered 
> volume from them based on the slice number. The volume building will 
> be upgraded to something more sophisticated in the future.
I'm suspecting a bug in directory/string parsing in DICOMImageReader
Is there any way to help direct DICOMImageReader on how to get it right?

Thanks,

-Stewart Dickson, Lecturer, Department of Computer Science, UNC-Asheville
Research Associate, National Environmental Modeling and Analysis Center 
http://www.unca.edu/nemac
Technical Director, http://us.imdb.com/Name?Stewart+Dickson
Sculptor, http://emsh.calarts.edu/~mathart/MathArt_siteMap.html


Program received signal EXC_BAD_ACCESS, Could not access memory.
Reason: KERN_INVALID_ADDRESS at address: 0xfffffff4
0x94b55014 in std::basic_stringbuf<char, std::char_traits<char>, 
std::allocator<char> >::_M_sync ()
(gdb) where  # we might ask...
#0  0x94b55014 in std::basic_stringbuf<char, std::char_traits<char>, 
std::allocator<char> >::_M_sync ()
#1  0x02ff8320 in vtkOStrStreamWrapper::vtkOStrStreamWrapper ()
#2  0x029f50d0 in vtkDICOMImageReader::CanReadFile ()
#3  0x029f5658 in vtkDICOMImageReader::ExecuteInformation ()
#4  0x0303a514 in vtkSource::UpdateInformation ()
#5  0x02fac9ec in vtkDataObject::UpdateInformation ()
#6  0x0303a414 in vtkSource::UpdateInformation ()
#7  0x02fac9ec in vtkDataObject::UpdateInformation ()
#8  0x01aa66dc in vtkAbstractVolumeMapper::Update ()
#9  0x01aa6780 in vtkAbstractVolumeMapper::GetBounds ()
#10 0x01b05028 in vtkVolume::GetBounds ()
#11 0x01aee56c in vtkRenderer::ComputeVisiblePropBounds ()
#12 0x01aee7d8 in vtkRenderer::ResetCamera ()
#13 0x01aee230 in vtkRenderer::GetActiveCamera ()
#14 0x01ab68e0 in vtkFrustumCoverageCuller::Cull ()
#15 0x01aedd98 in vtkRenderer::AllocateTime ()
#16 0x01aed6d4 in vtkRenderer::Render ()
#17 0x01af0c2c in vtkRendererCollection::Render ()
#18 0x01aea498 in vtkRenderWindow::DoStereoRender ()
#19 0x01aea408 in vtkRenderWindow::DoFDRender ()
#20 0x01ae9fd0 in vtkRenderWindow::DoAARender ()
#21 0x01ae987c in vtkRenderWindow::Render ()
#22 0x01b31390 in vtkXOpenGLRenderWindow::Render ()
#23 0x01187ad0 in vtkXOpenGLRenderWindowCppCommand ()
#24 0x0a00c818 in TclInvokeStringCommand ()
#25 0x0a00de84 in TclEvalObjvInternal ()
#26 0x0a00f430 in Tcl_EvalEx ()
#27 0x0a04ac18 in Tcl_FSEvalFile ()
#28 0x0a04ad08 in Tcl_EvalFile ()
#29 0x0b01fc28 in Tk_MainEx ()
#30 0x000031d8 in main ()

Machine Model:    PowerMac7,3
  CPU Type:    PowerPC G5  (3.1)
  Number Of CPUs:    2
  CPU Speed:    2.7 GHz
  L2 Cache (per CPU):    512 KB
  Memory:    2 GB
  Bus Speed:    1.35 GHz
 System Version:    Mac OS X 10.4.8 (8L127)
  Kernel Version:    Darwin 8.8.0
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: vtk_SaT1.tcl
URL: <http://www.vtk.org/pipermail/vtkusers/attachments/20070124/f600145a/attachment.txt>