<br clear="all">Hello Every one<br>&nbsp; I get into some problem when using VTK marching cube algorithm to deal with the standard CT HEAD test data, I get a messy visualization result after my marching cuble script Here is my script, and I cant figure out where is the problem, and I cant get any help on how to solve it, can any body help me ?<br>
&nbsp; The script below is my code , any body point me some potential problem? thanks very much<br><br>catch {load <span href="http://www.cadcaecam.com/Forum/tag.php?name=vtk" onclick="tagshow(event)" class="t_tag">vtk</span>tcl}<br>

# get the interactor ui<br>
#source ../../examplesTcl/vtkInt.tcl<br>
#source ../../examplesTcl/colors.tcl<br>
<br>
# Create the RenderWindow, Renderer and both Actors<br>
<br>
vtkRenderer ren1<br>
vtkRenderWindow renWin<br>
&nbsp; &nbsp; renWin AddRenderer ren1<br>
vtkRenderWindowInteractor iren<br>
&nbsp; &nbsp; iren SetRenderWindow renWin<br>
vtkLight lgt<br>
<br>
# create pipeline<br>
#<br>
vtkVolume16Reader v16<br>
&nbsp; &nbsp; v16 SetDataDimensions 256 256 <br>
&nbsp; &nbsp; [v16 GetOutput] SetOrigin 0.0 0.0 0.0<br>
&nbsp; &nbsp; #v16 SetDataByteOrderToLittleEndian<br>
&nbsp; &nbsp; v16 SetFilePrefix &quot;CThead&quot;<br>
&nbsp; &nbsp; v16 SetImageRange 1 113<br>
&nbsp; &nbsp; #v16 SetFilePattern &quot;CThead.%3d&quot;<br>
<br>
&nbsp; &nbsp; v16 SetDataSpacing 1 1 2<br>
vtkMarchingCubes iso<br>
&nbsp; &nbsp; iso SetInput [v16 GetOutput]<br>
&nbsp; &nbsp; iso SetValue 0 100<br>
vtkPolyDataMapper isoMapper<br>
&nbsp; &nbsp; isoMapper SetInput [iso GetOutput]<br>
&nbsp; &nbsp; isoMapper ScalarVisibilityOff<br>
vtkActor isoActor<br>
&nbsp; &nbsp; isoActor SetMapper isoMapper<br>
set isoProp [isoActor GetProperty]<br>
eval $isoProp SetColor 120 120 0<br>
<br>
vtkOutlineFilter outline<br>
&nbsp; &nbsp; outline SetInput [v16 GetOutput]<br>
vtkPolyDataMapper outlineMapper<br>
&nbsp; &nbsp; outlineMapper SetInput [outline GetOutput]<br>
vtkActor outlineActor<br>
&nbsp; &nbsp; outlineActor SetMapper outlineMapper<br>
set outlineProp [outlineActor GetProperty]<br>
#eval $outlineProp SetColor 0 0 0<br>
<br>
# Add the actors to the renderer, set the background and size<br>
#<br>
ren1 AddActor outlineActor<br>
ren1 AddActor isoActor<br>
ren1 SetBackground 1 1 1<br>
ren1 AddLight lgt<br>
renWin SetSize 1000 1000<br>
ren1 SetBackground 0.1 0.2 0.4<br>
<br>
set cam1 [ren1 GetActiveCamera]<br>
$cam1 Elevation 90<br>
$cam1 SetViewUp 0 0 -1<br>
$cam1 Zoom 1.3<br>
eval lgt SetPosition [$cam1 GetPosition]<br>
eval lgt SetFocalPoint [$cam1 GetFocalPoint]<br>
<br>
# render the image<br>
#<br>
#iren SetUserMethod {wm deiconify .vtkInteract}<br>
<br>
renWin Render<br>
#renWin SetFileName &quot;headBone.tcl.ppm&quot;<br>
#renWin SaveImageAsPPM<br>
<br>
# prevent the tk window from showing up then start the event loop<br>
wm withdraw .<br><br><br>your reply will be appreciated.<br>-- <br>~~~~~~~~~~~~~~~<br>no mistakes<br>~~~~~~~~~~~~~~~~~~<br>