<div dir="ltr">Romain,<div><br></div><div>Let's pop up a level and discuss what you are trying to do. It sounds like you are trying to voxelize a geometry into a 3D image. Is that right?</div><div><br></div><div>You don't have to use the vtkCutter at all but rather set your geometry as the input for <span style="font-size:12.8px">vtkPolyDataToImageStencil. I have used this class for voxelization in the past and it works very well.</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">HTH,</span></div><div><span style="font-size:12.8px">Cory</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 6, 2015 at 9:32 AM, Romain LEGUAY <span dir="ltr"><<a href="mailto:romain.leguay@gmail.com" target="_blank">romain.leguay@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thanks David for your answer!<div><br></div><div>I just saw on the vtkPolyDataToImageStencil documentation that the contour must be align throw the Z planes.<br></div><div><br></div><div>I use vtkCutter to cut my mesh along Z-axis but all the pictures saved are black. I attached my code and the CMakeLists.txt.</div><div><br></div><div>Is there another way to "project" my filled contour on an image?</div><div><br></div><div>Thanks again.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Romain</div><div><br></div><div><br></div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">2015-11-06 14:15 GMT+01:00 David Gobbi <span dir="ltr"><<a href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Romain,<div><br></div><div>The vtkPolyDataToImageStencil class needs a "hollow" shape as input,</div><div>either a 3D surface or a 2D contour. If you fill the contour, then it can't</div><div>be used by vtkPolyDataToImageStencil.</div><span><font color="#888888"><div><br></div><div> - David</div></font></span><div><div><div><br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 6, 2015 at 2:22 AM, Romain LEGUAY <span dir="ltr"><<a href="mailto:romain.leguay@gmail.com" target="_blank">romain.leguay@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello everyone,<div><br></div><div>In purpose to create a series of pictures (with the same size 512*512 for example), I cut a mesh along an axis (X, Y or Z) with vtkCutter, then I triangulate those contours with vtkContourTriangulator. I obtain some "flat" vtkPolydata<i>s</i> (one per cut). Now, I want to save all those "flat" meshes to 2D images with fixed size.</div><div>I saw on the wiki (<a href="http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/PolyDataToImageData" target="_blank">http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/PolyDataToImageData</a>) how to convert a vtkPolydata in vtkImageData, but in this case, my vtkPolydata is in a plane (the output of vtkContourTriangulator) so the dimension of the white image is 0 or negative because of this part of code:</div><br><div><div>vtkSmartPointer<vtkImageData> whiteImage = vtkSmartPointer<vtkImageData>::New(); </div><div>double bounds[6];</div><div>polydata->GetBounds(bounds);</div><div>double spacing[3]; // desired volume spacing</div><div>spacing[0] = 0.5;</div><div>spacing[1] = 0.5;</div><div>spacing[2] = 0.5;</div><div>whiteImage->SetSpacing(spacing);</div><div><br></div><div>// compute dimensions</div><div>int dim[3];</div><div>for (int i = 0; i < 3; i++)</div><div>{</div><div> dim[i] = static_cast<int>(ceil((bounds[i * 2 + 1] - bounds[i * 2]) / spacing[i])); // One of the dimension is 0</div><div>}</div><div>whiteImage->SetDimensions(dim);</div><div>whiteImage->SetExtent(0, dim[0] - 1, 0, dim[1] - 1, 0, dim[2] - 1); // Here </div></div><div><br></div><div><br></div><div><div>I don't know if there is other choice to convert a polydata to an image but this code always give some black pictures.</div></div><div>I'm not sure if I use the good filter to convert my polydata to a 2D image but I don't find anything.</div><div><br></div><div>Do you have any clues?</div><div><br></div><div>Thanks,</div><div><br></div><div>Romain</div></div></blockquote></div></div></div></div></div></div>
</blockquote></div><br></div>
</div></div><br>_______________________________________________<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 VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" rel="noreferrer" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=vtkusers" rel="noreferrer" target="_blank">http://markmail.org/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/mailman/listinfo/vtkusers</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">Cory Quammen<br>R&D Engineer<br>Kitware, Inc.</div>
</div>