The Frangi filter must perform a three-dimensional Gaussian filter on the volume before applying the three-dimensional hessian. The current code reshapes the 3D volume into a 2D image and then applies a 2D Gaussian.
The segmentation function (psoct_vessel_graphing\vasculature\vesSegment.m) calls the Hessian filter (psoct_vessel_graphing\vasculature\Hessian3D.m), which calls (psoct_vessel_graphing\vasculature\imgaussian.m). This imgaussian function can be replaced with the latest Matlab function imgaussfilt3, which is optimized for applying a Gaussian on a 3D volume.
This function also supports thread-based environments and GPU arrays.
The Frangi filter must perform a three-dimensional Gaussian filter on the volume before applying the three-dimensional hessian. The current code reshapes the 3D volume into a 2D image and then applies a 2D Gaussian.
The segmentation function (psoct_vessel_graphing\vasculature\vesSegment.m) calls the Hessian filter (psoct_vessel_graphing\vasculature\Hessian3D.m), which calls (psoct_vessel_graphing\vasculature\imgaussian.m). This imgaussian function can be replaced with the latest Matlab function imgaussfilt3, which is optimized for applying a Gaussian on a 3D volume.
This function also supports thread-based environments and GPU arrays.