<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
Hi Christian,<div><br></div><div>I noticed your thread where you were trying to do the same. Similar to your example, I ended up calling generateValues to create N contours and mapped it to colors in a lookup table. Thanks for your example.</div><div><br></div><div>While this method works correctly for a 125x125x50 volume, a larger 312x312x300 volume of short int causes a segmentation fault in vtkDiscreteMarchingCubesComputeGradient(). The volume has approximately 250 disconnected ROIs labeled from 0-N. Is there a limit on how large the volume can be? Or will the labels cause this error since the class expects the scalars to adhere to a specific layout?&nbsp;</div><div><br></div><div>Thanks,</div><div>Thommen<br><br><hr id="stopSpelling">Date: Thu, 24 Feb 2011 14:15:35 +0530<br>Subject: Re: [vtkusers] Visualizing multiple isosurfaces using vtkContourFilter<br>From: karthik.krishnan@kitware.com<br>To: thkorah@hotmail.com<br>CC: vtkusers@vtk.org<br><br><div>&nbsp;Since this is a labelled image with several label (ids), you&nbsp;could use vtkDiscreteMarchingCubes once on each label of interest instead of vtkContourFilter. This will extract just the interface between the cells containing the label of interest and its neighbors, whatever id they may have.</div>

<div>&nbsp;</div>
<div>--</div>
<div>karthik<br><br></div>
<div class="ecxgmail_quote">On Thu, Feb 24, 2011 at 4:39 AM, Thommen Korah <span dir="ltr">&lt;<a href="mailto:thkorah@hotmail.com">thkorah@hotmail.com</a>&gt;</span> wrote:<br>
<blockquote class="ecxgmail_quote" style="padding-left:1ex;border-left:#ccc 1px solid">
<div>Thanks for your help. I decided to use the first method with separate filters. I have multiple segmented regions, and each voxel stores the id that it belongs to. So all cells belonging to object 3 will store the number 3. 
<div><br></div>
<div>But when I create cf3 filter and call cf3-&gt;setValue(0,3), it creates the isosurface using all cells. My intent is to create the surface only using cells storing 3. Any idea why this might be?</div>
<div><br></div>
<div>Thanks.</div>
<div><br>&gt; Date: Wed, 23 Feb 2011 08:23:29 -0500<br>&gt; Subject: Re: [vtkusers] Visualizing multiple isosurfaces using vtkContourFilter<br>&gt; From: <a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a><br>
&gt; To: <a href="mailto:thkorah@hotmail.com">thkorah@hotmail.com</a><br>&gt; CC: <a href="mailto:vtkusers@vtk.org">vtkusers@vtk.org</a> 
<div>
<div></div>
<div class="h5"><br>&gt; <br>&gt; There are two ways to do this:<br>&gt; <br>&gt; 1) Create two ContourFilters. Use cf1-&gt;SetValue(0,val1) and<br>&gt; cf2-&gt;SetValue(0,val2).<br>&gt; You will need 2 mappers and actors.<br>
&gt; 2) Create one ContourFilter. Use cf-&gt;SetValue(0,val1) and cf-&gt;SetValue(1,val2).<br>&gt; You will need 1 mapper and actor.<br>&gt; <br>&gt; Both use share the same memory.<br>&gt; <br>&gt; Coloring is easier in case 1).<br>
&gt; mapper1-&gt;SetInput(cf1-&gt;GetOutput());<br>&gt; mapper1-&gt;ScalarVisibilityOff();<br>&gt; actor1-&gt;SetMapper(mapper1);<br>&gt; actor1-&gt;GetProperty()-&gt;SetDiffuseColor(r,g,b);<br>&gt; mapper2-&gt;SetInput(cf2-&gt;GetOutput());<br>
&gt; mapper2-&gt;ScalarVisibilityOff();<br>&gt; actor2-&gt;SetMapper(mapper2);<br>&gt; actor2-&gt;GetProperty()-&gt;SetDiffuseColor(r,g,b);<br>&gt; <br>&gt; Case 2 is trickier and involves mucking with lookup tables.<br>&gt; <br>
&gt; <br>&gt; <br>&gt; On Wed, Feb 23, 2011 at 1:16 AM, Thommen Korah &lt;<a href="mailto:thkorah@hotmail.com">thkorah@hotmail.com</a>&gt; wrote:<br>&gt; &gt; Hi,<br>&gt; &gt; I have a pointcloud on which I perform segmentation to extract different<br>
&gt; &gt; structures projecting out of a flat surface (a table for example). The<br>&gt; &gt; pointcloud is stored in a voxel data structure. I would now like to<br>&gt; &gt; visualize these segmented regions by creating &nbsp;isosurfaces using<br>
&gt; &gt; vtkContourFilter for each object. How can I generate multiple surfaces (eg.<br>&gt; &gt; book and a cup lying on the table) using the same shared data in memory? I<br>&gt; &gt; would also like to visualize each segmented region in a different color.<br>
&gt; &gt; Thanks.<br>&gt; &gt; _______________________________________________<br>&gt; &gt; Powered by <a href="http://www.kitware.com/" target="_blank">www.kitware.com</a><br>&gt; &gt;<br>&gt; &gt; Visit other Kitware open-source projects at<br>
&gt; &gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>&gt; &gt;<br>&gt; &gt; Please keep messages on-topic and check the VTK FAQ at:<br>
&gt; &gt; <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>&gt; &gt;<br>&gt; &gt; Follow this link to subscribe/unsubscribe:<br>&gt; &gt; <a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
&gt; &gt;<br>&gt; &gt;<br></div></div></div></div><br>_______________________________________________<br>Powered by <a href="http://www.kitware.com/" target="_blank">www.kitware.com</a><br><br>Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br><br>Follow this link to subscribe/unsubscribe:<br><a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
<br></blockquote></div><br></div>                                               </body>
</html>