1. A method for image enhancement comprising: performing a Gaussian convolution of voxels of an input image; calculating a Hessian matrix of the convolution; performing an Eigen decomposition of the Hessian matrix to derive eigenvectors and eigenvalues of the voxels; sorting the eigenvectors according to the magnitude of the corresponding eigenvalues; and performing a tubularity computation to calculate the probability of a voxel of the image being part of a tubular structure representation in the image whereby voxels determined to be part of tubular structures are used to provide an enhanced depiction of the image.