A method to compute a linear medial representation of a complex surface in the 3D discrete space is presented. The method involves voxel classification, surface labeling, anchor point detection, and voxel removal. they border sheets; and internal voxels, otherwise. We call internal sheets the sheets delimited exclusively by junction voxels, and peripheral sheets the sheets at least partially delimited by edge voxels, see Fig. 1. 1