Hi John<div>Remember to reply also to the mailing list so your question may help others in the future.</div><div><br></div><div>When I was referring to the value of the segmented object I meant the label (pixel value) assigned to the region segmented by the region growing filter. (1 or 255 usually)</div>
<div><br></div><div>Doing the segmentation in 3D or in 2D slices may affect the segmentation result, as the same object may be connected in a 3D image but it may appear disconnected in a particular slice. In some cases performing slice by slice segmentation might prevent unwanted leakage but on some other cases, segmentation might be incomplete.</div>
<div>I'd say that 3D region growing is the first thing to try and only use slice by slice on very specific cases.</div><div><br></div><div>Best regards</div><div><br><br><div class="gmail_quote">On Thu, Sep 15, 2011 at 3:05 PM, John Hiller <span dir="ltr"><<a href="mailto:john.hiller53@gmail.com">john.hiller53@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">I've checked paraview, I'm able to see the first mha file (slice), it looks like it segmented that particular slice correctly. I'm reading my dicom files with the image type set to dimension 2 and writing each file into mha's with the same dimension, should I be reading and writing in 3D? I wouldn't think so since I'm processing each slice in turn so a 2D representation is necessary.<div>
<br></div><div>I'm not sure what you mean by the value of my segmented object, what are you referring to as the segmented object (the whole mha series or one slice) and where can I find its value?</div><div><br></div>
<div>I'm doing this slice by slice because I need to extract particular bones from the dicom representation of a skull to get a 3D visualization, but, I want the option to include or exclude fragments of these bones. The major influence for this choice is because of the interaction/HCI I've designed for the user.</div>
<div><br></div><div>Thanks for the help<div><div></div><div class="h5"><br><br><div class="gmail_quote">On Thu, Sep 15, 2011 at 1:45 PM, Sergio Vera <span dir="ltr"><<a href="mailto:sergio.vera@alma3d.com" target="_blank">sergio.vera@alma3d.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<div><br></div><div>First, check if the problem is related with segmentation or visualization.</div><div>Open your mha files with paraview/seg3d2 or whatever 3D viewer you are used to. Check if the segmentation is ok. What value has your segmented object? 1?, 255? create a isosurface to the appropiate value. </div>
<div>Once checked if you are segmenting the slices ok, move to the vtk side and see if the fast marching is working properly or you are using the correct isovalues.</div><div><br></div><div>PS: any special reason to perform the segmentation slice by slice instead of segmenting the whole volume in a single step?</div>
<div>Regards</div><div><br></div><div>Sergio<br><div><br><div class="gmail_quote"><div>On Thu, Sep 15, 2011 at 2:36 PM, John Hiller <span dir="ltr"><<a href="mailto:john.hiller53@gmail.com" target="_blank">john.hiller53@gmail.com</a>></span> wrote:<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div>Hi itk users<div><br></div><div>I'm new to itk and vtk (sorry if this question doesn't apply to itk).</div>
<div>
<br></div><div>Okay, so this is what I'm doing:</div><div>1) read a dicom series, I process each dicom slice individually</div>
<div>2) pick the segmentation points on a particular dicom slice</div><div>3) segment using confidence connected filter and save (or write to a file in other words) the result to a .mha file</div><div>4) Continuously do this for each dicom slice to get series of .mha files</div>
<div>5) Once all that is done, I attempt to use Marching cubes from vtk to render a 3D surface image from the .mha files.</div><div><br></div><div>Everything works with no errors, its just that the vk window comes up black (showing nothing). I'm pretty sure the .mha files contain something (they're round about the same size as the original dicom file). I did the cout << reader->getOutput(), and it outputs information. I'm presuming that I am segmenting these slices incorrectly.</div>
<div><br></div><div>Any advice would be helpful</div><div>Thanks</div><div>John</div>
<br></div></div><div>_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br></div></blockquote></div><font color="#888888"><br><br clear="all"><div><br></div>-- <br>Sergio Vera<br><br> Alma IT Systems<br> C/ Vilana, 4B, 4º 1ª<br> 08022 Barcelona<br> T. <a href="tel:%28%2B34%29%20932%20380%20592" value="+34932380592" target="_blank">(+34) 932 380 592</a><br>
<a href="http://www.alma3d.com" target="_blank">www.alma3d.com</a><br>
</font></div></div>
</blockquote></div><br></div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Sergio Vera<br><br> Alma IT Systems<br> C/ Vilana, 4B, 4º 1ª<br> 08022 Barcelona<br> T. (+34) 932 380 592<br> <a href="http://www.alma3d.com">www.alma3d.com</a><br>
</div>