<div><div dir="auto">Hello,</div></div><div dir="auto"><br></div><div dir="auto">Did you set the dimensions based on your own data?  It looks all dimensions indicate 3.</div><div dir="auto"><br></div><div dir="auto">Thanks</div><div><br><div class="gmail_quote"><div dir="ltr">2018年8月3日(金) 4:50 Lizeth Castellanos <<a href="mailto:castellanoslizan@gmail.com">castellanoslizan@gmail.com</a>>:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>Thanks for help me.</div><div>I have confirmed the out.vti using Paraview's volume rendering. I don't know what is wrong, I get a solid cube again :-(</div><div>I am attaching a screenshot. Could you see, please?</div><div><br></div><div>Thanks</div></div><div dir="ltr"><div>Lizeth<br></div></div><div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 1, 2018 at 10:20 PM, kenichiro yoshimi <span dir="ltr"><<a href="mailto:rccm.kyoshimi@gmail.com" target="_blank">rccm.kyoshimi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
Thanks for trying.<br>
Could you confirm if out.vti is visualized properly using ParaView's<br>
volume rendering?<br>
<br>
Regards<br>
2018年8月2日(木) 6:12 Lizeth Castellanos <<a href="mailto:castellanoslizan@gmail.com" target="_blank">castellanoslizan@gmail.com</a>>:<br>
<div class="m_8285101241440346349HOEnZb"><div class="m_8285101241440346349h5">><br>
> Hi,<br>
><br>
> I got VTK 7.1.1 in my computer, finally! So, I have tested the two solutions that you have suggested:<br>
> 1) using vtkProbeFilter<br>
> 2) using vtkResampleToImage<br>
><br>
> 1) I  have adjusted the image dimensions according to my own data and I am using  vtkSmartVolumeMapper, instead of vtkVolumeRaycastMapper.<br>
> No errors are generated in the terminal. However, I have been waiting for several minutes and nothing happens.<br>
> Do you have any idea what is wrong?<br>
><br>
><br>
> csvfile="directory path to my csv file/sample.csv"<br>
><br>
> reader = vtk.vtkDelimitedTextReader()<br>
> reader.SetFieldDelimiterCharacters(" ")<br>
> reader.DetectNumericColumnsOn()<br>
> reader.SetFileName(csvfile)<br>
> reader.MergeConsecutiveDelimitersOn()<br>
> reader.Update()<br>
><br>
> tableToPoints = vtk.vtkTableToPolyData()<br>
> tableToPoints.SetInputConnection(reader.GetOutputPort())<br>
> tableToPoints.SetXColumn('Field 0')<br>
> tableToPoints.SetYColumn('Field 1')<br>
> tableToPoints.SetZColumn('Field 2')<br>
><br>
> imageData = vtk.vtkImageData()<br>
> imageData.SetDimensions(512, 512, 192)<br>
> imageData.SetOrigin(0.0, 0.0, 0.0)<br>
> imageData.SetSpacing(1.0, 1.0, 1.0)<br>
><br>
> probeFilter = vtk.vtkProbeFilter()<br>
> probeFilter.SetSourceConnection(tableToPoints.GetOutputPort())<br>
> probeFilter.SetInputData(imageData)<br>
> probeFilter.Update()<br>
><br>
> writer = vtk.vtkXMLImageDataWriter()<br>
> writer.SetInputData(probeFilter.GetOutput())<br>
> writer.SetFileName('out.vti')<br>
> writer.Write()<br>
><br>
> probeFilter.GetOutput().GetPointData().SetActiveScalars('Field 3')<br>
> range = probeFilter.GetOutput().GetPointData().GetScalars().GetRange()<br>
><br>
> colors = vtk.vtkNamedColors()<br>
><br>
> # Create the standard renderer, render window<br>
> # and interactor.<br>
> ren = vtk.vtkRenderer()<br>
><br>
> renWin = vtk.vtkRenderWindow()<br>
> renWin.AddRenderer(ren)<br>
><br>
> iren = vtk.vtkRenderWindowInteractor()<br>
> iren.SetRenderWindow(renWin)<br>
><br>
> # Create transfer mapping scalar value to opacity.<br>
> opacityTransferFunction = vtk.vtkPiecewiseFunction()<br>
> opacityTransferFunction.AddPoint(range[0], 0.0)<br>
> opacityTransferFunction.AddPoint(range[1], 1.0)<br>
><br>
> # Create transfer mapping scalar value to color.<br>
> colorTransferFunction = vtk.vtkColorTransferFunction()<br>
> colorTransferFunction.AddRGBPoint(range[0], 0.0, 0.0, 1.0)<br>
> colorTransferFunction.AddRGBPoint((range[0]+range[1])*0.5, 0.0, 1.0, 0.0)<br>
> colorTransferFunction.AddRGBPoint(range[1], 1.0, 0.0, 0.0)<br>
><br>
> # The property describes how the data will look.<br>
> volumeProperty = vtk.vtkVolumeProperty()<br>
> volumeProperty.SetColor(colorTransferFunction)<br>
> volumeProperty.SetScalarOpacity(opacityTransferFunction)<br>
> volumeProperty.SetScalarOpacityUnitDistance(1.73205080757)<br>
> volumeProperty.SetInterpolationTypeToLinear()<br>
><br>
> # The mapper / ray cast function know how to render the data.<br>
> volumeMapper = vtk.vtkSmartVolumeMapper()<br>
> volumeMapper.SetInputConnection(probeFilter.GetOutputPort())<br>
> volumeMapper.SetBlendModeToComposite()<br>
><br>
> # The volume holds the mapper and the property and<br>
> # can be used to position/orient the volume.<br>
> volume = vtk.vtkVolume()<br>
> volume.SetMapper(volumeMapper)<br>
> volume.SetProperty(volumeProperty)<br>
><br>
> ren.AddVolume(volume)<br>
><br>
> ren.SetBackground(colors.GetColor3d("White"))<br>
> ren.GetActiveCamera().ParallelProjectionOn()<br>
> ren.ResetCameraClippingRange()<br>
> ren.ResetCamera()<br>
><br>
> renWin.SetSize(600, 600)<br>
> renWin.Render()<br>
> iren.Start()<br>
><br>
><br>
> 2) I have tested the vtkResampleToImage class, (the solution that you suggest at the beginning of this thread). However, I get a solid a cube as output.<br>
> Do you have any idea why that happens?<br>
><br>
> csvfile="directory path to my csv file/sample.csv"<br>
><br>
> reader = vtk.vtkDelimitedTextReader()<br>
> reader.SetFieldDelimiterCharacters(" ")<br>
> reader.DetectNumericColumnsOn()<br>
> reader.SetFileName(csvfile)<br>
> reader.MergeConsecutiveDelimitersOn()<br>
> reader.Update()<br>
><br>
> tableToPoints = vtk.vtkTableToPolyData()<br>
> tableToPoints.SetInputConnection(reader.GetOutputPort())<br>
> tableToPoints.SetXColumn("Field 0")<br>
> tableToPoints.SetYColumn("Field 1")<br>
> tableToPoints.SetZColumn("Field 2")<br>
> tableToPoints.Update()<br>
><br>
> resample = vtk.vtkResampleToImage()<br>
> resample.SetInputConnection(tableToPoints.GetOutputPort())<br>
> resample.SetSamplingDimensions(3, 3, 3)<br>
> resample.<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div></div></div><div dir="ltr"><div class="gmail_extra">-- <br><div class="m_8285101241440346349gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><br></div><div><br></div></div></div></div></div></div>
</div></div>
</blockquote></div></div>