[Insight-users] 3D Vessel Extraction and Measurement

[Luis Ibanez]

Hi Ivan,


You may want to try the application

    http://www.itk.org/HTML/Curve3DExtraction.htm

This application is distributed with the InsightApplications
CVS checkout. (or the .zip and .tgz file in the release).

Note that this application assumes that you are targeting
a particular diameter of vessels.

If what you need is to estimate the diameter of vessels in
the image, then you may want to run the Vesselness filter
for a range of different sigmas (which are associated)
to the vessel diameter) and then choose for each vessel
the sigma where you get the highest response.

Note that this is a typical Multi-Scale analysis.

The drawback is that this may not be a fast processing, since
at all stages you process the entire image.

A more effective approach may involved to just perform the
multi-scale analysis on pixels that you have already spotted
as belonging to a vessel.  In that case you will have a two
pass processing. First, to identify prospective pixels belonging
to vessels, then Second visit to pixels and perform multi-scale
analysis on them, in order to estimate their diameters.


Another method for segmenting vessels is to use itk SpatialObjects
and register/fit  them to your images. For this purpose you may use
the itk Tube SpatialObjects.


  Please let us know if you have further questions.



      Thanks



        Luis



------------------
Iván Macía wrote:
> Hi,
> 
> I'm trying to perform vessel segmentation on MR liver images using some kind
> of medialness function to detect vessel centerlines and estimate radius
> similar to the work of Krissian et al. on "Model Based Detection of Tubular
> Structures on 3D Images". 
> 
> I have found in itk classes for Hessian Matrix computation
> (HessianRecursiveGaussianImageFilter) and measures of vesselness based on
> the eigenvalues of the Hessian Matrix
> (Hessian3DToVesselnessMeasureImageFilter). I have also found a
> VesselTubeSpatialObject suitable to represent the vessel information. But it
> does not seem to be a class suitable for estimation of vessel medialness,
> radius and some other properties.
> 
> I would like to know if there is any algorithm in itk that I could use, if
> someone is working on it or some advice that could be useful to implement
> one (probably in the form of an ImageToImageFilter).
> 
> Thanks in advance for your help
> 
> Iván Macía
>