<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.E-mailStijl17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="NL" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Dear all,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I have an image mask, from which I need the surface, the volume, the flatness etc.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">For this, I wanted to transform my mask in mesh. In order to calculate the volume, I wanted to transform the normal mesh to a simplex mesh using the TriangleToSimplexFilter.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">This is my code:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">typedef itk::Mesh<float, 3> MeshType;<br>
typedef itk::SimplexMesh<float, 3> TSimplex;<br>
typedef itk::TriangleMeshToSimplexMeshFilter<MeshType, TSimplex> TriangleToSimplexFilter;<br>
typedef itk::SimplexMeshVolumeCalculator<TSimplex> TVolume;<br>
typedef unsigned char PixelType;<br>
typedef itk::BinaryMask3DMeshSource<ImageType, MeshType> MeshSourceType;<br>
typename MeshSourceType::Pointer meshSource = MeshSourceType::New();<br>
const PixelType objectValue = static_cast<PixelType>(1);<br>
meshSource->SetObjectValue(objectValue);<br>
meshSource->SetInput(mask);<br>
try{<br>
meshSource->Update();<br>
}<br>
catch(itk::ExceptionObject &exp){<br>
std::cerr<<"Exception throwing during update()"<<std::endl;<br>
}<br>
// Ensure that all cells of the mesh are triangles.<br>
for (MeshType::CellsContainerIterator it = meshSource->GetOutput()->GetCells()->Begin();<br>
it != meshSource->GetOutput()->GetCells()->End();<br>
++it)<br>
{<br>
MeshType::CellAutoPointer cell;<br>
meshSource->GetOutput()->GetCell(it->Index(), cell);<br>
if (3 != cell->GetNumberOfPoints())<br>
{<br>
std::cerr << "ERROR: All cells must be trianglar." << std::endl;<br>
<br>
}<br>
}<br>
<br>
TriangleToSimplexFilter::Pointer Transform2Simplex = TriangleToSimplexFilter::New();<br>
Transform2Simplex->SetInput(meshSource->GetOutput());<br>
Transform2Simplex->Update();<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">However, when I run the code, it is calculating and never stopping. Would there be another approach? And does someone knows where the problem could be?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Thanks in advance for the help<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Elli<o:p></o:p></span></p>
</div>
<hr>
De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk
gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
<br>
<br>
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held
responsible for incomplete reception or delay of this transferred message.
</body>
</html>