Hi Arnaud,<br><br>thanks for your reply. I used the QuadEdgeMeshNormalFilter to generate the input mesh. Maybe I'm misinterpreting the documentation of the normal filter, but doesn't this filter need a vector type as template parameter? See (in the description, section "todo: <a href="http://www.itk.org/Doxygen/html/classitk_1_1QuadEdgeMeshNormalFilter.html">http://www.itk.org/Doxygen/html/classitk_1_1QuadEdgeMeshNormalFilter.html</a><br>
<br><div>stefan<br><div><br>On Mon, Dec 13, 2010 at 5:04 PM, Arnaud GELAS <<a href="mailto:arnaud_gelas@hms.harvard.edu">arnaud_gelas@hms.harvard.edu</a>> wrote:<br>> Hi Stefan,<br>><br>> The first template argument is used to define the type of the point's<br>
> coordinate (a concept checking may be missing here?).<br>> In your case, I guess it would be double instead of itk::Vector< double, 3<br>>>.<br>><br>> Here, you can find an example of one special<br>
> QuadEdgeMeshToQuadEdgeMeshFilter (for decimation purpose),<br>><br>> <a href="http://www.itk.org/Wiki/ITK/Examples/Meshes/Decimation">http://www.itk.org/Wiki/ITK/Examples/Meshes/Decimation</a><br>><br>> Best,<br>
> Arnaud<br>><br>> On 12/13/2010 10:40 AM, Stefan Dänzer wrote:<br>>><br>>> Hi Luis,<br>>><br>>> thanks for your reply. I have tried to use the<br>>> itkQuadEdgeMeshToQuadEdgeMeshFilter already. The output is a a<br>
>> QuadEdgeMesh with no points and no faces. I'm working on itk 3.20 at<br>>> the moment. This is how I use the filter in my code right now:<br>>><br>>> typedef itk::QuadEdgeMesh<itk::Vector<double,3>, 3><br>
>> meshType;<br>>><br>>> typedef itk::QuadEdgeMeshToQuadEdgeMeshFilter<meshType, meshType><br>>> qeToQeFilterType;<br>>><br>>> qeToQeFilterType::Pointer qeToQeFilter = qeToQeFilterType::New();<br>
>><br>>> qeToQeFilter->SetInput(meanShapeQe);<br>>> qeToQeFilter->Update();<br>>><br>>> meshType::Pointer currentShape = qeToQeFilter->GetOutput();<br>>><br>>> Might the problem be that I use a non scalar PointType<br>
>> (itk::Vector<double,3>) ?<br>>><br>>> stefan<br>>><br>>> On Fri, Dec 10, 2010 at 1:56 PM, Luis Ibanez<<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>><br>
>> wrote:<br>
>>><br>>>> Hi Stefan,<br>>>><br>>>> the itkMeshToMeshFilter is an abstract class<br>>>> that only serves as base class for other filters.<br>>>><br>>>> See the inheritance diagram in:<br>
>>> <a href="http://www.itk.org/Doxygen/html/classitk_1_1MeshToMeshFilter.html">http://www.itk.org/Doxygen/html/classitk_1_1MeshToMeshFilter.html</a><br>>>><br>>>> You may want to try<br>>>><br>
>>> itkQuadEdgeMeshToQuadEdgeMeshFilter.txx<br>>>><br>>>> <a href="http://www.itk.org/Doxygen/html/classitk_1_1QuadEdgeMeshToQuadEdgeMeshFilter.html">http://www.itk.org/Doxygen/html/classitk_1_1QuadEdgeMeshToQuadEdgeMeshFilter.html</a><br>
>>><br>>>> and to use an itkQuadEdgeMesh type for the<br>>>> reader that gets the input mesh from the .meta<br>>>> file.<br>>>><br>>>><br>>>> Regards,<br>
>>><br>>>><br>>>> Luis<br>>>><br>>>><br>>>> ----------------------------------------------------------------<br>>>> On Thu, Dec 2, 2010 at 9:36 PM, Stefan Dänzer<<a href="mailto:stefan.daenzer@gmail.com">stefan.daenzer@gmail.com</a>><br>
>>> wrote:<br>>>>><br>>>>> Hi,<br>>>>> I'm trying to read an itk::Mesh from a .meta file and copy the points<br>>>>> and cells of this itk::Mesh to a itk::QuadEdgeMesh. My first approach<br>
>>>> was to use itk::MeshToMeshFilter to achive this. The problem is that<br>>>>> the output mesh is empty after running the filter. Here is a part of<br>>>>> my code:<br>>>>><br>
>>>> itk::MetaMeshConverter<3,double>* metaConv = new<br>>>>> itk::MetaMeshConverter<3,double>;<br>>>>> itk::MetaMeshConverter<3,double>::SpatialObjectPointer meshSO =<br>
>>>> metaConv->ReadMeta(inFilename.c_str());<br>>>>><br>>>>> // get the itk::mesh<br>>>>> typedef itk::QuadEdgeMesh<double, 3> qeMeshType;<br>>>>> typedef itk::MeshToMeshFilter<MeshType, qeMeshType> qeFilterType;<br>
>>>><br>>>>> qeMeshType::Pointer qeMesh;<br>>>>> qeFilterType::Pointer qeFilter = qeFilterType::New();<br>>>>> qeFilter->SetInput(mesh);<br>>>>> qeFilter->Update();<br>
>>>> qeMesh = qeFilter->GetOutput();<br>>>>><br>>>>> Am I missing something? Thanks for any advice.<br>>>>><br>>>>> Stefan<br>>>>> _____________________________________<br>
>>>> Powered by <a href="http://www.kitware.com">www.kitware.com</a><br>>>>><br>>>>> Visit other Kitware open-source projects at<br>>>>> <a href="http://www.kitware.com/opensource/opensource.html">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">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">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">http://www.itk.org/mailman/listinfo/insight-users</a><br>>>>><br>>><br>>><br>><br>><br><br><br><br>-- <br>--<br>Stefan Daenzer<br>
Körnerplatz 8<br>04107 Leipzig<br><br>Tel.: +49-157-84993879<br><br>"Work like you don't need the money, love like you've never been hurt and dance like no one is watching." - Randall G Leighton<br><br>
</div>
</div>