<div dir="ltr"><div class="gmail_extra">>  Do you have to use paraview gui?</div><div class="gmail_extra"><br></div><div class="gmail_extra">No, the task on hand is already accomplished so no rush. I merely used the occasion to further teach myself a bit about scientific visualization software and the alternatives.</div><div class="gmail_extra"><br></div><div class="gmail_extra">I wonder, is it possible to have the volume data as HSV? That is the "value" is natural the intensity and equals the alpha value and a second "label" array the gives Hue and Saturation? That would feel more natural to me than having a Red/Green/Blue volume data and having it to code back RGB using 1:1 colorTransfers.</div><div class="gmail_extra"><br></div><div class="gmail_extra">- Axel</div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 16, 2015 at 4:18 PM, Aashish Chaudhary <span dir="ltr"><<a href="mailto:aashish.chaudhary@kitware.com" target="_blank">aashish.chaudhary@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Axel<br>
<br>
We were just discussing this few days ago. Do you have to use paraview gui? Or paraview python will work? There might be a way in pv python route.<br>
<br>
Thanks<br>
<br>
Sent from my iPhone<br>
<div><div class="h5"><br>
> On Dec 16, 2015, at 5:01 AM, Axel Kittenberger <<a href="mailto:axkibe@gmail.com">axkibe@gmail.com</a>> wrote:<br>
><br>
> Hello,<br>
><br>
> the Python VTK script below gives results, like a expect. Now I'm trying to reproduce it in Paraview and I fail.<br>
><br>
> I get rather confused as soon I give it a source with more than one scalar component.<br>
> Paraview creates one color map; I would expect 3 color maps, 1 for every component.<br>
> What ever I do effective opacity is always 1. That is I see a black cube with it's surface the images<br>
> of the data on its surface (so far correct), however no effective volume rendering -- I selected Volume, "Ray Cast Only"<br>
> which should be vtkFixedPointVolumeRayCastMapper as far I understood.<br>
><br>
> I tried creating 3 readers with the different scalar components extracted first, but doesn't give<br>
> correct rendering -- neither in VTK nor in paraview, since multiple renderer seem to compete which gets to<br>
> draw a pixel instead of composing them together<br>
> like a red transparent + a green transparent should give a yellow transparent impression.<br>
><br>
> How can I achieve the same result with Paraview?<br>
><br>
> Kind regards, Axel<br>
><br>
> --------------------------<br>
> #!/usr/bin/env python<br>
><br>
> import vtk<br>
><br>
> # Create the standard renderer, render window and interactor<br>
> ren = vtk.vtkRenderer( )<br>
> renWin = vtk.vtkRenderWindow( )<br>
> renWin.AddRenderer( ren )<br>
> iren = vtk.vtkRenderWindowInteractor( )<br>
> iren.SetRenderWindow( renWin )<br>
><br>
> reader = vtk.vtkImageReader2( )<br>
> reader.SetFileName( "./data.raw" );<br>
> reader.SetDataScalarTypeToUnsignedChar( );<br>
> reader.SetNumberOfScalarComponents( 3 );<br>
> reader.SetFileDimensionality( 3 );<br>
> reader.SetDataExtent( 0, 261, 0, 511, 0, 21 );<br>
> reader.SetDataSpacing( 1, 1, 3.5 );<br>
><br>
> # Create transfer mapping scalar value to opacity<br>
> opacityTransferFunction = vtk.vtkPiecewiseFunction( )<br>
> opacityTransferFunction.AddPoint(   0, 0.0 )<br>
> opacityTransferFunction.AddPoint(  20, 0.0 )<br>
> opacityTransferFunction.AddPoint(  40, 0.1 )<br>
> opacityTransferFunction.AddPoint( 255, 0.1 )<br>
><br>
> # Create transfer mapping scalar value to color<br>
><br>
> colorRedTransferFunction = vtk.vtkColorTransferFunction( )<br>
> colorRedTransferFunction.AddRGBPoint(   0.0, 0.2, 0.0, 0.0 )<br>
> colorRedTransferFunction.AddRGBPoint( 255.0, 1.0, 0.0, 0.0 )<br>
><br>
> colorGreenTransferFunction = vtk.vtkColorTransferFunction( )<br>
> colorGreenTransferFunction.AddRGBPoint(   0.0, 0.0, 0.2, 0.0 )<br>
> colorGreenTransferFunction.AddRGBPoint( 255.0, 0.0, 1.0, 0.0 )<br>
><br>
> colorBlueTransferFunction = vtk.vtkColorTransferFunction( )<br>
> colorBlueTransferFunction.AddRGBPoint(   0.0, 0.0, 0.0, 0.0 )<br>
> colorBlueTransferFunction.AddRGBPoint( 255.0, 0.0, 0.0, 1.0 )<br>
><br>
> volumeProperty = vtk.vtkVolumeProperty( )<br>
> volumeProperty.SetColor( 0, colorRedTransferFunction )<br>
> volumeProperty.SetColor( 1, colorGreenTransferFunction )<br>
> volumeProperty.SetColor( 2, colorBlueTransferFunction )<br>
><br>
> volumeProperty.SetScalarOpacity( 0, opacityTransferFunction )<br>
> volumeProperty.SetScalarOpacity( 1, opacityTransferFunction )<br>
> volumeProperty.SetScalarOpacity( 2, opacityTransferFunction )<br>
> volumeProperty.SetInterpolationTypeToLinear( )<br>
><br>
> volumeMapper = vtk.vtkFixedPointVolumeRayCastMapper( )<br>
> volumeMapper.SetInputConnection( reader.GetOutputPort( ) )<br>
> volumeMapper.SetNumberOfThreads( 1 )<br>
><br>
> volume = vtk.vtkVolume( )<br>
> volume.SetMapper( volumeMapper )<br>
> volume.SetProperty( volumeProperty )<br>
><br>
> ren.AddVolume( volume )<br>
> ren.SetBackground( 1, 1, 1 )<br>
> renWin.SetSize( 600, 600 )<br>
><br>
> def CheckAbort( obj, event ):<br>
>     if obj.GetEventPending( ) != 0:<br>
>         obj.SetAbortRender( 1 )<br>
><br>
> renWin.AddObserver( "AbortCheckEvent", CheckAbort )<br>
><br>
> iren.Initialize( )<br>
> renWin.Render( )<br>
> iren.Start( )<br>
> --------------------------<br>
</div></div>> _______________________________________________<br>
> Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" rel="noreferrer" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
><br>
> Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" rel="noreferrer" target="_blank">http://markmail.org/search/?q=ParaView</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://public.kitware.com/mailman/listinfo/paraview" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><br>
</blockquote></div><br></div></div>