hessian3D {vesselr}

3D Volume Hessian

Description

This function returns the eigenvalues of the hessian matrices for a 3D array or NIfTI volume.

Usage

hessian3D(image, mask, radius = 1, parallel = FALSE, cores = 2)

Arguments

image

a 3D array or image of class nifti

mask

an array or nifti mask of voxels for which vesselness will be calculated, with more selective masking improving speed significantly. Note that mask should be in the same space as the image volume

radius

an integer specifying radius of the neighborhood (in voxels) for which the hessian should be calculated

parallel

is a logical value that indicates whether the user's computer is Linux or Unix (i.e. macOS), and should run the code in parallel

cores

if parallel = TRUE, cores is an integer value that indicates how many cores the function should be run on

Value

A list of three eigenvalue volumes.

Examples

## Not run: 
library(neurobase)
epi <- readnii('path/to/epi')
mask <- epi!=0
hesseigs <- hessian3D(image = epi, mask = mask) 
## End(Not run)