Problems with Frangi's filter

4 views (last 30 days)
Alessandro Russo
Alessandro Russo on 14 Mar 2016
Hi all, i have a cardiac CT volume of 512x512x125, and i want to apply Frangi's filter (downloaded from file exchange http://www.mathworks.com/matlabcentral/fileexchange/24409-hessian-based-frangi-vesselness-filter) in order to extract the coronary vesselness map. The problem is, i give as input the approximate vesselness map obtained with global tresholding (a thin band of the external heart wall on which there is coronaries course, with also the ventricles in some slices, all sourrounded by black background), setting the filter parameters as found in literature for this applications (alpha=0.5, beta=10, gamma=10.000, range of sigmas variable between 0 and 2 with step of 0.05), but the result (plotted with vol3d) is the one in figure, totally incorrect.
What am i doing wrong? Of course i tried a lot of different sets of parameters, but the result is always unsatisfying. I attach the code which applies the filter, the workspace with the variables useful for compiling and the tool vol3d for visualization if wanted.
mex eig3volume.c;
options.BlackWhite=false;
options.FrangiScaleRange=[0 3];
options.FrangiScaleRatio=0.01;
options.FrangiC=10000;
options.FrangiAlpha=0.5;
options.FrangiBeta=15;
Vfiltered=FrangiFilter3D(V,options);
figure(1), vol3d('cData',Vfiltered)

Answers (0)

Categories

Find more on Biomedical Imaging in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!