A computer-implemented method, comprising the steps of: obtaining image data generated by an imaging system; generating data indicative of a degree to which each of a plurality of voxels of the image data corresponds to one or more predetermined geometrical features, wherein each geometrical feature is assigned a different color; generating a signal indicative of a single color value for each of the plurality of voxels based on the degree and the colors; generating a volumetric rending of the image data based on the signal; generating a link between voxels of the volumetric rendering and voxels of the image data; and visually presenting the image data and the volumetric rendering concurrently, wherein the step of generating the data includes: estimating a quadric for each of the plurality of voxels, wherein the quadric includes a matrix of coefficients; determining eigenvalues of the quadrics based on the matrix of coefficients; and determining the degree for the voxels based on the eigenvalues, characterized in that the step of estimating the quadric for the voxel, comprises: casting a plurality of rays from the voxel, wherein each ray is cast in a different direction and traverses neighboring voxels of the image data; terminating a cast ray in response to satisfying predetermined stopping criteria; determining a length of each cast ray; and determining the matrix of coefficients based on the different directions and the determined lengths. obtaining image data generated by an imaging system; generating data indicative of a degree to which each of a plurality of voxels of the image data corresponds to one or more predetermined geometrical features, wherein each geometrical feature is assigned a different color; generating a signal indicative of a single color value for each of the plurality of voxels based on the degree and the colors; generating a volumetric rending of the image data based on the signal; generating a link between voxels of the volumetric rendering and voxels of the image data; and visually presenting the image data and the volumetric rendering concurrently, wherein the step of generating the data includes: estimating a quadric for each of the plurality of voxels, wherein the quadric includes a matrix of coefficients; determining eigenvalues of the quadrics based on the matrix of coefficients; and determining the degree for the voxels based on the eigenvalues, characterized in that the step of estimating the quadric for the voxel, comprises: casting a plurality of rays from the voxel, wherein each ray is cast in a different direction and traverses neighboring voxels of the image data; terminating a cast ray in response to satisfying predetermined stopping criteria; determining a length of each cast ray; and determining the matrix of coefficients based on the different directions and the determined lengths. estimating a quadric for each of the plurality of voxels, wherein the quadric includes a matrix of coefficients; determining eigenvalues of the quadrics based on the matrix of coefficients; and determining the degree for the voxels based on the eigenvalues, characterized in that the step of estimating the quadric for the voxel, comprises: casting a plurality of rays from the voxel, wherein each ray is cast in a different direction and traverses neighboring voxels of the image data; terminating a cast ray in response to satisfying predetermined stopping criteria; determining a length of each cast ray; and determining the matrix of coefficients based on the different directions and the determined lengths. casting a plurality of rays from the voxel, wherein each ray is cast in a different direction and traverses neighboring voxels of the image data; terminating a cast ray in response to satisfying predetermined stopping criteria; determining a length of each cast ray; and determining the matrix of coefficients based on the different directions and the determined lengths.