<p dir="ltr">I'm glad it worked out.</p>
<p dir="ltr">Cheers,<br>
</p>
<div class="gmail_extra"><br><div class="gmail_quote">On Aug 26, 2016 12:12 PM, "amdelaossa" <<a href="mailto:adelaossa@gmail.com">adelaossa@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hello Álvaro,<div><br></div><div>It works! </div><div>Rendering multiple overlapping data volumes with the same grid is possible with <wbr>vtkGPUVolumeRayCastMapper.</div><div>Using numpy arrays was easy to gather the different 3-dimensional data arrays together by adding a 4th index to the combined array. </div><div>The volume properties can be set independently for each one of the components.</div><div>Thanks for the help!</div><div><br></div><div>Best regards,</div><div>Alberto</div><div><br></div><div><br><div><blockquote type="cite"><div>On 26 Aug 2016, at 15:55, amdelaossa <<a href="mailto:adelaossa@gmail.com" target="_blank">adelaossa@gmail.com</a>> wrote:</div><br><div><blockquote type="cite" style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div><br>On 26 Aug 2016, at 15:43, Alvaro Sanchez <<a href="mailto:alvaro.sanchez@kitware.com" target="_blank">alvaro.sanchez@kitware.com</a>> wrote:</div><br><div><div dir="ltr"><div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word"><div><div><div>It is great to read that you are working on this feature for future releases of VTK. </div><div>I think that it will be very useful for many people. </div><div>In the mean time I would like to try your proposed approach of defining multiple components in one volume. </div><div>Let me ask you something:</div><div>If data1, data2, data3 are my three independent three-dimensional data arrays, should I merge them first in a four-dimensional array where the first index is the component index?</div><div><br></div><span></span></div></div></div></blockquote></div><div><br></div><div>You should merge them in a three-dimensional array (each component being dataN).  Checkout <wbr>TestGPURayCastTwoComponentsInd<wbr>ependent, that could serve you as an example.  Going back to your initial question though, as Elvis pointed out, if your volumes are not overlapping you could create instead several vtkVolumes (each with its own mapper), this should also work.</div></div></div></blockquote><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br></div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br></div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">My volumes are overlapping but they are defined in the same (regular) grid. </div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">So I need to use the multi-component trick. </div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">However I still do not understand the format…</div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Lets say that I have two 3D arrays data1 and data2 with identical number of components in each direction, i.e.  data1(Nx,Ny,Nz) and data2(Nx,Ny,Nz).</div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Should merge them into another array which would look like this: </div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">data(2 * Nx, Ny, Nz) </div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">with data(i,j,k) = data1(i,j,k) for i<Nx and </div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">data(i,j,k) = data2(i-Nx,j,k) for i>=Nx  </div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br></div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Is that right?</div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br></div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Thanks!</div><div style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Alberto</div><br style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><blockquote type="cite" style="font-family:Helvetica;font-size:14px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div><div dir="ltr"><div><br></div><div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word"><div><div>Could you give an rough estimation of when multiple volumes rendering with vtkGPUVolumeRayCastMapper will be available in VTK?<br></div></div><div><br></div></div></blockquote></div><div><br></div><div>My best guess would be VTK 8.</div><br><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 26, 2016 at 6:39 AM, amdelaossa<span> </span><span dir="ltr"><<a href="mailto:adelaossa@gmail.com" target="_blank">adelaossa@gmail.<wbr>com</a>></span><span> </span>wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word">Hi Álvaro,<div><br></div><div>Thanks a lot for your prompt answer!</div><div>See my comments below. </div><div><br><div><span><blockquote type="cite"><div>as you have noticed, multi-volume rendering is currently not supported in the vtkGPUVolumeRayCastMapper.  If you are building VTK with OpenGL2 however, it supports independent components so (depending on your use case) you might be able to render up to 4 volumes (one per component), provided you can arrange them in a four component array (this has the limitation that each component will need to have the same dimensions).  That said, we are working on giving support for multiple volumes in the next VTK release.</div></blockquote><div><br></div></span><div>It is great to read that you are working on this feature for future releases of VTK. </div><div>I think that it will be very useful for many people. </div><div>In the mean time I would like to try your proposed approach of defining multiple components in one volume. </div><div>Let me ask you something:</div><div>If data1, data2, data3 are my three independent three-dimensional data arrays, should I merge them first in a four-dimensional array where the first index is the component index?</div><span><div><br></div><blockquote type="cite"><div dir="ltr"><div>With regards to the solution you mentioned (<a href="https://github.com/bozorgi/VTKMultiVolumeRayCaster" target="_blank">https://github.com/bozorgi/VT<wbr>KMultiVolumeRayCaster</a>), you could use it temporarily if it would solve your problem better than fitting the volumes as independent components.  It seems to implement a similar approach as we would do (connecting various inputs/ properties to the mapper) so it should not be extremely difficult to migrate to the official implementation in the future.  I have not tried it myself but please do let us know if you have any remarks once you try it.  As far as I know, if the python bindings are not explicitly disabled for these new classes then they should be also wrapped.</div></div></blockquote><div><br></div></span><div>Alright. However this solution is currently not useful to me since it was implemented for VTK 5 (I think) and it fails to compile against VTK 7. </div><div>So I am temporarily stuck here. I’d like to try your approach with the multi-components…</div><div><br></div><div>Could you give an rough estimation of when multiple volumes rendering with vtkGPUVolumeRayCastMapper will be available in VTK?</div><div><br></div><div>Thanks a lot!</div><div><br></div><div>Regards,</div><div>Alberto</div><div><div><br><blockquote type="cite"><div dir="ltr"><div><br></div><div>Regads,</div><div>Álvaro</div><div><br></div><div><br></div><div><br></div><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-family:verdana,sans-serif;font-size:small;color:rgb(11,83,148)"><br></div><div class="gmail_quote">---------- Forwarded message ----------<br>From:<span> </span><b class="gmail_sendername">amdelaossa</b><span> </span><span dir="ltr"><<a href="mailto:adelaossa@gmail.com" target="_blank">adelaossa@<wbr>gmail.com</a>></span><br>Date: Thu, Aug 25, 2016 at 10:00 AM<br>Subject: [vtkusers] Multiple volumes rendering with vtkGPUVolumeRayCastMapper or similar<br>To:<span> </span><a href="mailto:vtkusers@vtk.org" target="_blank">vtkusers@vtk.org</a><br><br><br><div style="word-wrap:break-word">Dear all, <div>we (a <a href="http://plasma.desy.de/research/simulation_and_theory/index_eng.html" target="_blank">research team</a> in DESY, Germany) are strongly interested in rendering multiple data volumes in the same view with VTK. </div><div>We chose VTK for our 3D visualisation project on beam-plasma interactions due to its extended functionality and features, </div><div>its open source character and because it is written in C++ with bindings for python.<br>However, we will require to show multiple volumes in the same view (representing the data of different particle species in our simulations).<br>We are currently using a vtkGPUVolumeRayCastMapper object for each data set, which is then connected to a vtkVolume. </div><div>The different vtkVolume objects are then added to the vtkRenderer.<br>We have tried to use different vtkRenderer objects for the different vtkVolumes, </div><div>but still only the first volume on the first vtkRenderer is displayed. <br><br>Searching in the web, we have found the following (potential) solution:<br><a href="https://github.com/bozorgi/VTKMultiVolumeRayCaster" target="_blank">https://github.com/bozorgi/VTK<wbr>MultiVolumeRayCaster</a><br>which is published here:<br><a href="http://www.ncbi.nlm.nih.gov/pubmed/24841148" target="_blank">http://www.ncbi.nlm.nih.gov/pu<wbr>bmed/24841148</a><br><br>It consists on new vtk classes to handle multiple volumes rendering.<br>As it might work well for our purposes, we would like to ask your opinion about this new classes.<br>Are you planning to incorporate something like this in the official VTK distribution soon?<br>If not, we will be interested in incorporating this new class into our personal VTK installation.<br>Then the question is: Will be the python bindings automatically generated for these new classes?<br><br>Thanks a lot for your support!</div><div><br></div><div>Best regards,</div><div>Alberto</div><div><br></div></div><br>______________________________<wbr>_________________<br>Powered by<span> </span><a href="http://www.kitware.com/" rel="noreferrer" target="_blank">www.kitware.com</a><br><br>Visit other Kitware open-source projects at<span> </span><a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>opensource/opensource.html</a><br><br>Please keep messages on-topic and check the VTK FAQ at:<span> </span><a href="http://www.vtk.org/Wiki/VTK_FAQ" rel="noreferrer" target="_blank">http://www.vtk.org/Wiki/<wbr>VTK_FAQ</a><br><br>Search the list archives at:<span> </span><a href="http://markmail.org/search/?q=vtkusers" rel="noreferrer" target="_blank">http://markmail.org/<wbr>search/?q=vtkusers</a><br><br>Follow this link to subscribe/unsubscribe:<br><a href="http://public.kitware.com/mailman/listinfo/vtkusers" rel="noreferrer" target="_blank">http://public.kitware.com/mail<wbr>man/listinfo/vtkusers</a><br><br></div><br></div></blockquote></div><br><br clear="all"><div><br></div>--<span> </span><br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Alvaro Sanchez</div><div>Kitware, Inc.</div><div><span style="font-family:Arial;font-size:14.6667px;line-height:20.24px;white-space:pre-wrap">Senior </span>R&D Engineer</div><div>21 Corporate Drive</div><div>Clifton Park, NY 12065-8662</div><div>Phone:<span> </span><a href="tel:518-881-4901" value="+15188814901" target="_blank">518-881-4901</a></div></div></div></div></div></div></div></blockquote></div></div></div><br></div></div></blockquote></div><br><br clear="all"><div><br></div>--<span> </span><br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Alvaro Sanchez</div><div>Kitware, Inc.</div><div><span style="font-family:Arial;font-size:14.6667px;line-height:20.24px;white-space:pre-wrap">Senior </span>R&D Engineer</div><div>21 Corporate Drive</div><div>Clifton Park, NY 12065-8662</div><div>Phone: <a href="tel:518-881-4901" value="+15188814901" target="_blank">518-881-4901</a></div></div></div></div></div></div></div></div></blockquote></div></blockquote></div><br></div></div></blockquote></div></div>