cortex.polyutils.naive_prism_height

cortex.polyutils.naive_prism_height(flat, wm, pia, polys)

Height of a volume-preserving vertical prism over the flatmap.

The obvious reading of the bumpy flatmap idea: each column keeps its folded volume over a base of the flattened triangle. Kept for comparison – it spikes wherever flattening compressed a triangle hard, and smoothing a field of ratios afterwards does not remove spikes that dominate the mean. See the module docstring for the deeper problem with the denominator.

Parameters:

flat2D ndarray, shape (total_verts, 3)

Location of each vertex in flatmap space.

wm2D ndarray, shape (total_verts, 3)

Location of each vertex on the white matter surface.

pia2D ndarray, shape (total_verts, 3)

Location of each vertex on the pial surface.

polys2D ndarray, shape (total_polys, 3)

Triangle vertex indices, shared by all three surfaces.

Returns:

height1D ndarray, shape (total_verts,)

Height of the pial surface above the flatmap at each vertex, in the units of the input surfaces.