Contents

Metric Smoothing

Metric smoothing replaces a node's metric (or surface shape) value using a weighted, or unweighted, average of the node's metric value with the metric values of the node's neighbors. Smoothing essentially removes the noise (reduces any local minima or maxima). Smoothing can also be thought of as a low-pass filter.

Parameters Common to All Smoothing Algorithms

Parameters common to all smoothing algorithms:

Strength - Weight applied to average of neighbors (not used in Full Width Half Maximum nor Dilation algorithms)

Iterations - The number of iterations of smoothing

Several algorithms use inter-neighbor distances. It is important that the Fiducial (Anatomical) Surface is used for smoothing.

Average Neighbors

Average neighbors smooths a node by averaging its metric value with the metric values of the node's immediate neighbors.

For node Mx with N neighbors Mi, one iteration of smoothing is average of the node and its neighbors:

The number of iterations provided by the user is the maximum number of iterations of smoothing that are performed using the Full Width Half Maximum Algorithm. Prior to each iteration of smoothing, the Full Width Half Maximum is estimated using equations below. If the calculated FWHM exceeds the FWHM provided by user, smoothing ceases.

dv = Average Euclidean (straight-line) Inter-Neighbor Distance for all nodes.

var(ds) = Variance of node and neighbor metric differences.

var(s) = Variance of node metric values

FWHMsurf =

Gaussian

In gaussian smoothing, the neighboring node is assigned the product of two weights, the normal weight and the tangent weight. For the normal weight, x is the distance from the plane tangent to the node being smoothed to the neighboring node (this distance is parallel to the normal vector). If the node is not within the Normal Below and Normal Above Cutoffs, the normal weight is zero. For the tangent weight, x is the distance from the neighbor to the normal vector. If the the node is not within the Tangent Cutoffs, the tangent weight is zero. In the calculation of both weights, the weight is the gaussian equation where x is the normal or tangent distance.

Sigma Normal - The standard deviation for the gaussian used in the normal weight calculation.

Sigma Tangent - The standard deviation for the gaussian used in the tangent weight calculation.

Normal Below Cutoff - The cutoff distance in a direction opposite of the node's normal vector used durning the normal weight calculation.

Normal Above Cutoff - The cutoff distance in the direction of the node's normal vector used during the normal weight calculation.

Tangent Cutoff - The cutoff distance that is orthogonal to the normal vector.

For node Mx with N neighbors Mi, one iteration of smoothing is:

Wi is the gaussian weight for for neighbor i of the node being smoothed.

(the total of the weights).

Mx = (strength * MW) + ((1.0 − strength) * Mx)

For gaussian smoothing, all nodes within a "neighbor depth" of 5 are considered neighbors of the node being smoothed. If a fiducial surface was used for determining the distance from a node to its neighbor, the Euclidean (straight line) distance may be very short across a sulcus but the distance along the surface (geodesic distance) may be much longer. Unfortunately, the time required to calculate the geodesic distances is too slow. As a result, we have chosen to approximate the geodesic distance by computing the distance between a node and its neighbors using a spherical surface.

Weighted Average Neighbors

Weighted Average Neighbors is similar to Average Neighbors except that neighbors closer to the node being smoothed receive a higher weighting than neighbors further away.

Let:

Di = distance(Mx,Mi)

For node Mx with N neighbors Mi, one iteration of smoothing is :

A Smoothing Analysis

This example uses data from Jason Hill's neonate study. The average neighbors algorithm is used with 10 iterations and strength of 1.0.

Before Smoothing

Iterations=4, Strength=0.5

Iterations=10, Stength=1.0

To further quantify the amount of smoothing this represents, we created a metric map with zero values across the entire inter24 fiducial surface, right hemisphere. Next, we set the values of a handful of isolated nodes in known regions of interest to 10.0. Then, we smoothed the map 10 iterations at 1.0 strength. The resulting map and histogram are shown in the table below: