Compute the shortest path tortuosity distribution through the medial axis between boundary faces in an image. This is used to study network flow through a medial axis which connects boundary faces of a 3D image.