The convex wrapping algorithm: A method for identifying muscle paths using the underlying bone mesh

Associating musculoskeletal models to motion analysis data enables the determination of the muscular lengths, lengthening rates and moment arms of the muscles during the studied movement. Therefore, those models must be anatomically personalized and able to identify realistic muscular paths. Different kinds of algorithms exist to achieve this last issue, such as the wired models and the finite elements ones. After having studied the advantages and drawbacks of each one, we present the convex wrapping algorithm. Its purpose is to identify the shortest path from the origin to the insertion of a muscle wrapping over the underlying skeleton mesh while respecting possible non-sliding constraints. After the presentation of the algorithm, the results obtained are compared to a classically used wrapping surface algorithm (obstacle set method) by measuring the length and moment arm of the semitendinosus muscle during an asymptomatic gait. The convex wrapping algorithm gives an efficient and realistic way of identifying the muscular paths with respect to the underlying bones mesh without the need to define simplified geometric forms. It also enables the identification of the centroid path of the muscles if their thickness evolution function is known. All this presents a particular interest when studying populations presenting noticeable bone deformations, such as those observed in cerebral palsy or rheumatic pathologies.