[vtkusers] Marching Contours Error

[Mathieu Malaterre]

Could not reproduce the problem using mcubes.tcl:

vtkImageData img
   img ShallowCopy [v16 GetOutput]

so please send us the code you were using, thanks.

Mathieu
Script used:

package require vtk
package require vtkinteraction
package require vtktesting


# Create the RenderWindow, Renderer and both Actors
#
vtkRenderer ren1
vtkRenderWindow renWin
     renWin AddRenderer ren1
vtkRenderWindowInteractor iren
     iren SetRenderWindow renWin

# create pipeline
#
vtkVolume16Reader v16
   v16 SetDataDimensions 64 64
   [v16 GetOutput] SetOrigin 0.0 0.0 0.0
   v16 SetDataByteOrderToLittleEndian
   v16 SetFilePrefix "$VTK_DATA_ROOT/Data/headsq/quarter"
   v16 SetImageRange 1 93
   v16 SetDataSpacing 3.2 3.2 1.5
   v16 Update

vtkImageData img
   img ShallowCopy [v16 GetOutput]

vtkMarchingContourFilter iso
   iso SetInput img
   iso SetValue 0 1125

vtkPolyDataMapper isoMapper
   isoMapper SetInput [iso GetOutput]
   isoMapper ScalarVisibilityOff

vtkActor isoActor
   isoActor SetMapper isoMapper
   eval [isoActor GetProperty] SetColor $antique_white

vtkOutlineFilter outline
   outline SetInput [v16 GetOutput]
vtkPolyDataMapper outlineMapper
   outlineMapper SetInput [outline GetOutput]
vtkActor outlineActor
   outlineActor SetMapper outlineMapper
   outlineActor VisibilityOff

# Add the actors to the renderer, set the background and size
#
ren1 AddActor outlineActor
ren1 AddActor isoActor
ren1 SetBackground 0.2 0.3 0.4
renWin SetSize 200 200
[ren1 GetActiveCamera] Elevation 90
[ren1 GetActiveCamera] SetViewUp 0 0 -1
[ren1 GetActiveCamera] Dolly 1.5


ren1 ResetCameraClippingRange

renWin Render

# render the image
#
iren AddObserver UserEvent {wm deiconify .vtkInteract}

# prevent the tk window from showing up then start the event loop
wm withdraw .


gcope at engin.umich.edu wrote:
> Hi,
> 
> I don't seem to have much luck getting responses, but here goes:
> 
> I am getting an "Update Extent does not lie within Whole Extent" error upon
> execution of a Marching Contours filter. This also results in a memory access
> error in vtkImageMarchingCubes.cxx at the line "g[2] = (float)(ptr[inc2]) -
> (float)(ptr[-inc2]);"
> 
> What caused the error is attempting to use a ShallowCopy of a filters output as
> the input to the MarchingContour filter rather than directly connecting them.
> When they are directly connected, everything works fine. Does doing a
> ShallowCopy modify the data extent somehow?
> 
> Any ideas would be appreciated.
> 
> Thanks,
> Jamie
> 
> 
> _______________________________________________
> This is the private VTK discussion list. 
> Please keep messages on-topic. Check the FAQ at: <http://public.kitware.com/cgi-bin/vtkfaq>
> Follow this link to subscribe/unsubscribe:
> http://www.vtk.org/mailman/listinfo/vtkusers
>