Hi Xiaoli,
things got busy and this delayed my vmtk list time.
vmtkcenterlines works on any surface obtained in any way, only it needs to be
orientable (i.e. you must be able to define outwards normals, unlike, say, a
Klein bottle, https://en.wikipedia.org/wiki/Klein_bottle - this can also happen
due to small-scale degeneracies) and manifold (i.e. each triangle edge must be
share by at most two triangles).
In your case it was neither due to the segmentation artifacts inside the vessel.
So, yes, you can export an STL file from Amira, read it with vmtk and generate
centerlines.
Luca
On Aug 21, 2013, at 6:39 PM, Xiaoli Zhang wrote:
> Hi Luca,
>
> Thanks for pointing this out. Indeed this is the problem (at least it worked
> on the problematic dataset). It turned out there is some region inside the
> lumen not segmented, and that caused the problem. A quick question through:
> does vmtkcenterlines only work on surface created by
> vmtklevelsetsegmentation? Will it work on a manually segmented surface from
> Amira?
>
> Just FYI, following are my steps for segmentation and centerline creation:
> 1. segment vessel wall from CT image (different from typical cast CT, only
> vessel wall has hyperintensity)
> 2. fill the lumen to include both vessel wall segmented above and its lumen
> 3. use vmtklevelsetsegmentation to generate deformable model from step 2's
> label field (not sure if it is necessary)
> 4. create centerline
>
> Xiaoli
>
>
>
>
> On Fri, Aug 16, 2013 at 7:59 AM, Luca Antiga <luca.ant...@orobix.com> wrote:
> Hi Xiaoli,
> I took a look at your surface. It's not immediately obvious, but there is a
> whole extra portion of surface inside of it (see screenshot taken with some
> translucency). This causes the algorithm to basically not work.
> I'm not sure how you got to this point, but I suggest you try to re-segment
> the surface making sure you obtain a single compact 3D model.
> Best,
>
> Luca
>
> <Screen Shot 2013-08-16 at 1.55.20 PM.png>
>
> On Aug 16, 2013, at 3:55 AM, Xiaoli Zhang wrote:
>
>> Hi Luca,
>>
>> Attached is the surface file I am working on. The levelset segmentation
>> wasn't working very well on the outer curvature of the arch, but that
>> shouldn't be the problem. Even if I am only working on descending aorta, it
>> still won't work.
>>
>> Thanks,
>> Xiaoli
>>
>>
>> On Thu, Aug 15, 2013 at 5:09 PM, Luca Antiga <luca.ant...@gmail.com> wrote:
>> Hi, I see that there is a jagged line on the surface indicating that the
>> surface is apparently not orientable, so normals flip where you see the dark
>> line (since the shader interpolates normals). You can verify this in
>> Paraview by applying the Glyph filter to the normals.
>>
>> This would explain dificulties with the centerline algorithm, since normals
>> are used when extracting the inner Delaunay diagram.
>>
>> Can you send the surface over? I can try to understand what's causing
>> normals to flip.
>>
>> Thanks
>>
>> Luca
>>
>>
>> On 12/ago/2013, at 17:09, Xiaoli Zhang <zhang...@gmail.com> wrote:
>>
>>> Hi Luca,
>>>
>>> Attached is the screenshots of the surface for centerline extraction from
>>> two different views. I have tried not to include aortic arch, but it didn't
>>> work either. Any suggestions will be appreciated.
>>>
>>> Thanks,
>>> Xiaoli
>>>
>>> On Sun, Aug 11, 2013 at 4:24 PM, Luca Antiga <luca.ant...@orobix.com> wrote:
>>> Hi Xiaoli,
>>> there must be a Voronoi diagram degeneracy somewhere. Feel free to post a
>>> screenshot of the model, I can try to point you to the candidate spots.
>>> Best,
>>>
>>> Luca
>>>
>>>
>>> On Aug 8, 2013, at 7:50 PM, Xiaoli Zhang wrote:
>>>
>>>> It seems someone has encountered the same problem before:
>>>> http://www.mail-archive.com/vmtk-users@lists.sourceforge.net/msg00283.html
>>>>
>>>> I followed Luca's suggestions apply smoothing (tried passband of 0.01 and
>>>> 0.1) and add a subdivision step before vmtkcenterlines. Still no luck,
>>>> vmtkcenterlines script will never return. Any suggestions?
>>>>
>>>> On Tue, Aug 6, 2013 at 9:23 PM, Xiaoli Zhang <zhang...@gmail.com> wrote:
>>>> I am new to vmtk and I am trying to extract centerline from CT image of
>>>> the aorta. I prepared the surface in Amira and save as stl (ASCII) format.
>>>> Then I ran vmtkcenterlines script and place a source point at inlet and a
>>>> target point at outlet. vmtk then show "computing centerline" and staying
>>>> for a very long time (>12 hours). The surface I am working on is a small
>>>> file (stl size <1MB).
>>>>
>>>> Anyone can help please? I can send the screenshot of the surface and stl
>>>> file if needed. Thanks.
>>>>
>>>> Xiaoli
>>>>
>>>> ------------------------------------------------------------------------------
>>>> Get 100% visibility into Java/.NET code with AppDynamics Lite!
>>>> It's a free troubleshooting tool designed for production.
>>>> Get down to code-level detail for bottlenecks, with <2% overhead.
>>>> Download for free and get started troubleshooting in minutes.
>>>> http://pubads.g.doubleclick.net/gampad/clk?id=48897031&iu=/4140/ostg.clktrk_______________________________________________
>>>> vmtk-users mailing list
>>>> vmtk-users@lists.sourceforge.net
>>>> https://lists.sourceforge.net/lists/listinfo/vmtk-users
>>>
>>>
>>> <C9N3_surface_centerline.jpg>
>>> ------------------------------------------------------------------------------
>>> Get 100% visibility into Java/.NET code with AppDynamics Lite!
>>> It's a free troubleshooting tool designed for production.
>>> Get down to code-level detail for bottlenecks, with <2% overhead.
>>> Download for free and get started troubleshooting in minutes.
>>> http://pubads.g.doubleclick.net/gampad/clk?id=48897031&iu=/4140/ostg.clktrk
>>> _______________________________________________
>>> vmtk-users mailing list
>>> vmtk-users@lists.sourceforge.net
>>> https://lists.sourceforge.net/lists/listinfo/vmtk-users
>>
>> <C9N3.ls.mc.sm.vtp>
>
>
------------------------------------------------------------------------------
Introducing Performance Central, a new site from SourceForge and
AppDynamics. Performance Central is your source for news, insights,
analysis and resources for efficient Application Performance Management.
Visit us today!
http://pubads.g.doubleclick.net/gampad/clk?id=48897511&iu=/4140/ostg.clktrk
_______________________________________________
vmtk-users mailing list
vmtk-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/vmtk-users