Function puma::compute_orientationComparison

Function Documentation

std::pair<double, double> puma::compute_orientationComparison(Workspace *ws, puma::MatVec3<double> *dirs, puma::MatVec3<double> *tangents, puma::Matrix<double> *error, puma::Cutoff cutoff, int numThreads = 0)

Function to compute the error between the correct local orientation (tangents puma matric) of an artificially generated sample and the estimated one (dirs puma matrix). Error is between 0 - 90 degrees. Set error to 100 if NaN inside dirs is found. Set error to 200 if NaN inside tangents is found. Note that these voxels are not considered when computing the Mean angular difference and Standard Deviation (SD)

Parameters

  • ws – the input workspace to serve as a mask according to the cutoff

  • dirs – first orientation field

  • tangents – second orientation field

  • error – matrix containing the angular difference (in degrees)

  • cutoff – the grayscale range to consider, according to the input ws

  • numThreads – number of threads to use

Returns

pair containing the mean (first) and std (second)