Measurement of the size and orientation of human masseter and medial pterygoid muscles

To gain a better understanding of biting and chewing performance, the size and orientation of the masseter and medial pterygoid muscles in living humans were studied. Twenty-seven young males having complete dentition, class I dental occlusion and normal muscle and jaw function were examined using magnetic resonance images of the head between the zygomatic arch and hyoid bone. The sections were parallel to the palatal plane, and the thickness was 3 mm without a gap. A computer software program (Medical Dental Image, MDI) was developed to identify and calculate the area of each cross section of the muscle, and the volume of the muscle was then estimated. The axis of the muscle was determined by connecting the centroids of the sections in the lower and upper 1/3 of the whole muscle. The effective muscle cross section area was then calculated by resectioning the muscle perpendicularly to the muscle axis. It was found that the mean masseter muscle volume was around 31 cm3, and that the mean medial pterygoid muscle volume was 11 cm3. Their mean effective cross section areas were around 6.2 cm2 and 3.5 cm2, respectively. The axis of the masseter muscle was more perpendicular to the palatal plane and parallel to the sagittal plane than was the medial pterygoid muscle. The results suggest that the use of magnetic resonance images (MRI) is an effective noninvasive measurement technique for determining the size and orientation of masseter and medial pterygoid muscles. This technique can be employed in future studies on human bite force evaluation and masticatory function.     

Musculoskeletal model of the upper limb based on the visible human male dataset

A mathematical model of the human upper limb was developed based on high-resolution medical images of the muscles and bones obtained from the Visible Human Male (VHM) project. Three-dimensional surfaces of the muscles and bones were reconstructed from Computed Tomography (CT) images and Color Cryosection images obtained from the VHM cadaver. Thirteen degrees of freedom were used to describe the orientations of seven bones in the model: clavicle, scapula, humerus, radius, ulna, carpal bones, and hand. All of the major articulations from the shoulder girdle down to the wrist were included in the model. The model was actuated by 42 muscle bundles, which represented the actions of 26 muscle groups in the upper limb. The paths of the muscles were modeled using a new approach called the Obstacle-set Method [33]. The calculated paths of the muscles were verified by comparing the muscle moment arms computed in the model with the results of anatomical studies reported in the literature. In-vivo measurements of maximum isometric muscle torques developed at the shoulder, elbow, and wrist were also used to estimate the architectural properties of each musculotendon actuator in the model. The entire musculoskeletal model can be reconstructed using the data given in this paper, along with information presented in a companion paper which defines the kinematic structure of the model [26].     

Musculoskeletal Model of the Upper Limb Based on the Visible Human Male Dataset

A mathematical model of Ihe human upper limb was developed based on high-resolution medical images of the muscles and bones obtained from the Visible Human Male ( HM) project. Three-dimensional surfaces of the muscles and bones were reconstructed from Computed Tomography (CT) images and Color Cryosection images obtained from the VHM cadaver. Thirteen degrees of freedom were used to describe the orientations of seven bones in the model: clavicle, scapula, humerus, radius, ulna, carpal bones, and hand. All of the major articulations from the shoulder girdle down to the wrist were included in the model. The model was actuated by 42 muscle bundles, which represented the actions of 26 muscle groups in the upper limb. The paths of the muscles were modeled using a new approach called the Obstacle-set Method (33) The calculated paths of the muscles were verified by comparing the muscle moment arms computed in the model with the results of anatomical studies reported in the literature, In-vivo measurements of maximum isometric muscle torques developed at the shoulder, elbow, and wrist were also used to estimate the architectural properties of each musculotendon actuator in the model. The entire musculoskeletal model can be reconstructed using the data given in this paper, along with information presented in a companion paper which defines the kinematic structure of the model (26)     

Volumic patient-specific reconstruction of muscular system based on a reduced dataset of medical images

Three-dimensional mechanical modelling of muscles is essential for various biomechanical applications and clinical evaluation, but it requires a tedious manual processing of numerous images. A muscle reconstruction method is presented based on a reduced set of images to generate an approximate parametric object from basic dimensions of muscle contours. A regular volumic mesh is constructed based on this parametric object. The approximate object and the corresponding mesh are deformed to fit the exact muscles contours yielding patient-specific geometry. Evaluation was performed by comparison of geometry to that obtained by contouring all computed tomography (CT) slices, and by quantification of the mesh quality criteria. Muscle fatty infiltration was estimated using a threshold between fat and muscle. Volumic fat index (VFI) of a muscle was computed using first all the complete CT scan slices containing the muscle (VFI(ref)) and a second time only the slices used for reconstruction (VFI(recons)). Mean volume error estimation was 2.6% and hexahedron meshes fulfilled quality criteria. VFI(recons) respect the individual variation of fat content.     

Contractile and non-contractile tissue volume and distribution in ankle muscles of young and older adults

Magnetic resonance imaging (MRI) enables accurate in vivo quantification of human muscle volumes, which can be used to estimate subject-specific muscle force capabilities. An important consideration is the amount of contractile and non-contractile tissue in the muscle compartment, which will influence force capability. We quantified age-related differences in the proportion and distribution of contractile and non-contractile tissue in the dorsiflexor and plantar flexor (soleus, and medial and lateral heads of gastrocnemius) muscles, and examined how well these volumes can be estimated from single MRI cross-sections. Axial MRIs of the left leg for 12 young (mean age 27 years) and 12 older (72 years) healthy, active adults were used to compute muscle volumes. Contractile tissue distribution along the leg was characterized by mathematical functions to allow volume prediction from single-slice cross-sectional area (CSA) measurements. Compared to young, older adults had less contractile volume and a greater proportion of non-contractile tissue. In both age groups the proportion of non-contractile tissue increased distally, with the smallest proportion near the maximum compartment CSA. A single CSA measurement predicted contractile volume with 8-11% error, with older adults in the higher end of this range. Using multiple slices improved volume estimates by roughly 50%, with average errors of about 3-4%. These results demonstrate significant age-related differences in non-contractile tissue for the dorsi- and plantar-flexor muscles. Although estimates of contractile volume can be obtained from single CSA measurements, multiple slices are needed for increased accuracy due to inter-individual variations in muscle volume and composition.     

A method for estimating torque-vector directions of shoulder muscles using surface EMGs

In this study, a new method is proposed to estimate the torque-vector directions of each shoulder muscle. The method is based on a multiple regression model that reconstructs shoulder torque, which is calculated from the hand force and posture, from the surface EMG of many muscles recorded simultaneously. The torque-vector directions of eleven shoulder muscles of four subjects were obtained at up to 30 different arm postures with this method. The mean confidence interval (p < 0.05) of the estimated torque-vector direction of each subject was 7.7-10.6 degrees. The correlation coefficient between the measured shoulder torque and reconstructed shoulder torque was between 0.76-0.84. The results for majority of the muscles were in accordance with previous studies, and reasonable from the viewpoint of anatomy. The torque-vector directions of a muscle, which are estimated with this method, have more of a functional meaning than a pure anatomical or mechanical one. These indicate the direction of the shoulder torque accompanying the muscle activation for a normal shoulder action that involves the cooperative contraction of many muscles.     

A method for estimating torque-vector directions of shoulder muscles using surface EMGs

In this study, a new method is proposed to estimate the torque-vector directions of each shoulder muscle. The method is based on a multiple regression model that reconstructs shoulder torque, which is calculated from the hand force and posture, from the surface EMG of many muscles recorded simultaneously. The torque-vector directions of eleven shoulder muscles of four subjects were obtained at up to 30 different arm postures with this method. The mean confidence interval ( p< 0.05) of the estimated torque-vector direction of each subject was 7.7-10.6 degrees. The correlation coefficient between the measured shoulder torque and reconstructed shoulder torque was between 0.76-0.84. The results for majority of the muscles were in accordance with previous studies, and reasonable from the viewpoint of anatomy. The torque-vector directions of a muscle, which are estimated with this method, have more of a functional meaning than a pure anatomical or mechanical one. These indicate the direction of the shoulder torque accompanying the muscle activation for a normal shoulder action that involves the cooperative contraction of many muscles.     

中国数字化可视人体鞍上池的横断面解剖与CT、MRI对照研究

目的:探讨鞍上池在中国数字化可视人体(Chinese visible human,CVH)与CT、MRI上的横断面解剖形态学表现。方法:选择做64层螺旋CT和MRI头部检查的健康志愿者各60例,获得5mm层厚横断面图像。从第2例中国数字化可视人体数据集中选取与CT、MRI相对应层面的头部薄层连续横断面标本图像,对照观察鞍上池在CVH、MRI与CT图像上的正常解剖形态、毗邻及内部结构。结果:CVH图像上,鞍上池表现为六角形和五角形两种形状。CVH薄层横断面图像能连续、清晰地显示鞍上池的正常形态、毗邻及内部结构。60例CT及MRI图像上,鞍上池全部显示,但解剖结构均不及CVH清晰。鞍上池在CT、MRI横断面图像上形状变化更大,以六角形最多,五角形次之,四角形最少,相应毗邻及内部结构也有所不同。六角形鞍上池在CVH、CT、MRI上有良好的对应关系,五角形鞍上池部分相匹配,CVH图像上无四角形鞍上池。结论:通过与CT、MRI进行对照研究,中国数字化可视人体能为颅脑疾病的影像识别和诊断提供断层解剖学依据。     

Reliability and robustness of muscle architecture measurements obtained using diffusion tensor imaging with anatomically constrained tractography

For detailed analyses of muscle adaptation mechanisms during growth, ageing or disease, reliable measurements of muscle architecture are required. Diffusion tensor imaging (DTI) and DTI tractography have been used to reconstruct the architecture of human muscles in vivo. However, muscle architecture measurements reconstructed with conventional DTI techniques are often anatomically implausible because the reconstructed fascicles do not terminate on aponeuroses, as real muscle fascicles are known to do. In this study, we tested the reliability of an anatomically constrained DTI-based method for measuring three-dimensional muscle architecture. Anatomical magnetic resonance images and diffusion tensor images were obtained from the left legs of eight healthy participants on two occasions one week apart. Muscle volumes, fascicle lengths, pennation angles and fascicle curvatures were measured in the medial and lateral gastrocnemius, soleus and the tibialis anterior muscles. Averaged across muscles, the intraclass correlation coefficient was 0.99 for muscle volume, 0.81 for fascicle length, 0.73 for pennation angle and 0.76 for fascicle curvature. Measurements of muscle architecture obtained using conventional DTI tractography were highly sensitive to variations in the stopping criteria for DTI tractography. The application of anatomical constraints reduced this sensitivity significantly. This study demonstrates that anatomically constrained DTI tractography can provide reliable and robust three-dimensional measurements of whole-muscle architecture. The algorithms used to constrain tractography have been made publicly available.     

Sarcoptes scabiei var. hominis: Three-dimensional structure of a female imago and crusted scabies lesions by X-ray micro-CT

The three-dimensional structure of scabies mites (Sarcoptes scabiei var. hominis) and keratin layers affected by crusted scabies lesions were obtained using X-ray computed tomography at sub-micrometer and micrometer resolution, respectively (X-ray micro-CT). Clear three-dimensional images including internal structure of scabies mites were obtained. Utilizing reconstructed micro-CT data, the sections of the capitulum (head part), digestive organs, and legs are shown. The reconstructed capitulum shows a jaw-like structure capable of penetrating the keratin layer of the skin. The tip of the forelegs of female scabies mites has a flat disk structure that may be used to grasp the skin surface. The keratin layer of a crusted scabies lesion spontaneously exfoliated from a patient was also reconstructed by the X-ray micro-CT technique. Extracted sections from CT data revealed a network structure of tunnels made by scabies mites with numerous larvae and eggs inside the tunnels.     

Quadriceps femoris muscle morphology and function after ACL injury: a differential response in copers versus non-copers

The morphology (volume and peak cross-sectional area) and voluntary muscle control of 27 athletic people were evaluated with magnetic resonance imaging (MRI) and an established method of testing neuromuscular control in order to explain why some people are able to cope with anterior cruciate ligament (ACL) injury (copers), whereas most cannot (non-copers). Axial spin-echo T1 weighted MRI images were acquired from the level of the ankle mortise to the iliac crest. The subjects' quadriceps, hamstrings, and gastrocnemius muscles were digitally reconstructed from the MRI images. The volume and peak cross-sectional area (CSA) of each muscle were then calculated. Voluntary muscle control was evaluated using an established target-matching protocol that requires subjects to produce and modulate force with control over a range of directions. Electromyographic signals were collected from seven muscles as the subjects performed the experiment. Circular statistics methods were used to calculate a specificity index that describes how focused the activity pattern of each muscle was with respect to its principal direction of action. The results of the non-copers, copers, and uninjured subjects were then compared. The non-copers displayed significantly greater quadriceps atrophy than the copers. The most profound differences were observed in the vastus lateralis muscle. The non-copers also displayed diminished vastus lateralis and lateral gastrocnemius muscle control. Little differences were observed in the results of the copers and uninjured subjects. In general, the copers' results fell between those of the non-copers and uninjured subjects. The results of this study suggest that quadriceps muscle function is a critical factor in the differential response to ACL injury.     

Measuring morphological parameters of the pelvic floor for finite element modelling purposes

The goal of this study was to obtain a complete data set needed for studying the complex biomechanical behaviour of the pelvic floor muscles using a computer model based on the finite element (FE) theory. The model should be able to predict the effect of surgical interventions and give insight into the function of pelvic floor muscles. Because there was a lack of any information concerning morphological parameters of the pelvic floor muscle structures, we performed an experimental measurement to uncover those morphological parameters. Geometric parameters as well as muscle parameters of the pelvic floor muscles were measured on an embalmed female cadaver. A three-dimensional (3D) geometric data set of the pelvic floor including muscle fibre directions was obtained using a palpator device. A 3D surface model based on the experimental data, needed for mathematical modelling of the pelvic floor, was created. For all parts of the diaphragma pelvis, the optimal muscle fibre length was determined by laser diffraction measurements of the sarcomere length. In addition, other muscle parameters such as physiological cross-sectional area and total muscle fibre length were determined. Apart from these measurements we obtained a data set of the pelvic floor structures based on nuclear magnetic resonance imaging (MRI) on the same cadaver specimen. The purpose of this experiment was to discover the relationship between the MRI morphology and geometrical parameters obtained from the previous measurements. The produced data set is not only important for biomechanical modelling of the pelvic floor muscles, but it also describes the geometry of muscle fibres and is useful for functional analysis of the pelvic floor in general. By the use of many reference landmarks all these morphologic data concerning fibre directions and optimal fibre length can be morphed to the geometrical data based on segmentation from MRI scans.These data can be directly used as an input for building a mathematical model based on FE theory.     

IVUS Validation of Patient Coronary Artery Lumen Area Obtained from CT Images

Aims

Accurate computed tomography (CT)-based reconstruction of coronary morphometry (diameters, length, bifurcation angles) is important for construction of patient-specific models to aid diagnosis and therapy. The objective of this study is to validate the accuracy of patient coronary artery lumen area obtained from CT images based on intravascular ultrasound (IVUS).

Methods and Results

Morphometric data of 5 patient CT scans with 11 arteries from IVUS were reconstructed including the lumen cross sectional area (CSA), diameter and length. The volumetric data from CT images were analyzed at sub-pixel accuracy to obtain accurate vessel center lines and CSA. A new center line extraction approach was used where an initial estimated skeleton in discrete value was obtained using a traditional thinning algorithm. The CSA was determined directly without any circular shape assumptions to provide accurate reconstruction of stenosis. The root-mean-square error (RMSE) for CSA and diameter were 16.2% and 9.5% respectively.

Conclusions

The image segmentation and CSA extraction algorithm for reconstruction of coronary arteries proved to be accurate for determination of vessel lumen area. This approach provides fundamental morphometric data for patient-specific models to diagnose and treat coronary artery disease.     

Method and apparatus for soft tissue material parameter estimation using tissue tagged Magnetic Resonance Imaging

We describe an experimental method and apparatus for the estimation of constitutive parameters of soft tissue using Magnetic Resonance Imaging (MRI), in particular for the estimation of passive myocardial material properties. MRI tissue tagged images were acquired with simultaneous pressure recordings, while the tissue was cyclically deformed using a custom built reciprocating pump actuator A continuous three-dimensional (3D) displacement field was reconstructed from the imaged tag motion. Cavity volume changes and local tissue microstructure were determined from phase contrast velocity and diffusion tensor MR images, respectively. The Finite Element Method (FEM) was used to solve the finite elasticity problem and obtain the displacement field that satisfied the applied boundary conditions and a given set of material parameters. The material parameters which best fit the FEM predicted displacements to the displacements reconstructed from the tagged images were found by nonlinear optimization. The equipment and method were validated using inflation of a deformable silicon gel phantom in the shape of a cylindrical annulus. The silicon gel was well described by a neo-Hookian material law with a single material parameter C1=8.71+/-0.06kPa, estimated independently using a rotational shear apparatus. The MRI derived parameter was allowed to vary regionally and was estimated as C1 =8.80+/-0.86kPa across the model. Preliminary results from the passive inflation of an isolated arrested pig heart are also presented, demonstrating the feasibility of the apparatus and method for isolated heart preparations. FEM based models can therefore estimate constitutive parameters accurately and reliably from MRI tagging data.     

Maximizing quantitative accuracy of lung airway lumen and wall measures obtained from X-ray CT imaging.

To objectively quantify airway geometry from three-dimensional computed tomographic (CT) images, an idealized (circular cross section) airway model is parameterized by airway luminal caliber, wall thickness, and tilt angle. Using a two-dimensional CT slice, an initial guess for the airway center, and the full-width-half-maximum principle, we form an estimate of inner and outer airway wall locations. We then fit ellipses to the inner and outer airway walls via a direct least squares fit and use the major and minor axes of the ellipses to estimate the tilt and in-plane rotation angles. Convolving the airway model, initialized with these estimates, with the three-dimensional scanner point-spread function forms the predicted image. The difference between predicted and actual images is minimized by refining the model parameter estimates via a multidimensional, unconstrained, nonlinear minimization routine. When optimization converges, airway model parameters estimate the airway inner and outer radii and tilt angle. Results using a Plexiglas phantom show that tilt angle is estimated to within +/-4 degrees and both inner and outer radii to within one-half pixel when a "standard" CT reconstruction kernel is used. By opening up the ability to measure airways that are not oriented perpendicular to the scanning plane, this method allows evaluation of a greater sampling of airways in a two-dimensional CT slice than previously possible. In addition, by combining the tilt-angle compensation with the deconvolution method, we provide significant improvement over the previous full-width-half-maximum method for assessing location of the luminal edge but not the outer edge of the airway wall.     

Comparison between line and surface mesh models to represent the rotator cuff muscle geometry in musculoskeletal models

Accurate muscle geometry (muscle length and moment arm) is required to estimate muscle function when using musculoskeletal modelling. In shoulder, muscles are often modelled as a collection of independent line segments, leading to non-physiological muscles trajectory, especially for the rotator cuff muscles. To prevent this, a surface mesh model was developed and validated against 7 MRI positions in one participant. Mean moment arm errors was 11.4% for the line vs. 8.8% for the mesh model. While the model with independent lines led to some non-physiological trajectories, the mesh model gave lower misestimations of muscle lengths and moment arms.     

The validity of ultrasound estimation of muscle volumes

The purpose of this study was to validate ultrasound muscle volume estimation in vivo. To examine validity, vastus lateralis ultrasound images were collected from cadavers before muscle dissection; after dissection, the volumes were determined by hydrostatic weighing. Seven thighs from cadaver specimens were scanned using a 7.5-MHz ultrasound probe (SSD-1000, Aloka, Japan). The perimeter of the vastus lateralis was identified in the ultrasound images and manually digitized. Volumes were then estimated using the Cavalieri principle, by measuring the image areas of sets of parallel two-dimensional slices through the muscles. The muscles were then dissected from the cadavers, and muscle volume was determined via hydrostatic weighing. There was no statistically significant difference between the ultrasound estimation of muscle volume and that estimated using hydrostatic weighing (p > 0.05). The mean percentage error between the two volume estimates was 0.4% +/- 6.9. Three operators all performed four digitizations of all images from one randomly selected muscle; there was no statistical difference between operators or trials and the intraclass correlation was high (>0.8). The results of this study indicate that ultrasound is an accurate method for estimating muscle volumes in vivo.     

Three-dimensional temporomandibular joint muscle attachment morphometry and its impacts on musculoskeletal modeling

In musculoskeletal models of the human temporomandibular joint (TMJ), muscles are typically represented by force vectors that connect approximate muscle origin and insertion centroids (centroid-to-centroid force vectors). This simplification assumes equivalent moment arms and muscle lengths for all fibers within a muscle even with complex geometry and may result in inaccurate estimations of muscle force and joint loading. The objectives of this study were to quantify the three-dimensional (3D) human TMJ muscle attachment morphometry and examine its impact on TMJ mechanics. 3D muscle attachment surfaces of temporalis, masseter, lateral pterygoid, and medial pterygoid muscles of human cadaveric heads were generated by co-registering measured attachment boundaries with underlying skull models created from cone-beam computerized tomography (CBCT) images. A bounding box technique was used to quantify 3D muscle attachment size, shape, location, and orientation. Musculoskeletal models of the mandible were then developed and validated to assess the impact of 3D muscle attachment morphometry on joint loading during jaw maximal open-close. The 3D morphometry revealed that muscle lengths and moment arms of temporalis and masseter muscles varied substantially among muscle fibers. The values calculated from the centroid-to-centroid model were significantly different from those calculated using the ‘Distributed model’, which considered crucial 3D muscle attachment morphometry. Consequently, joint loading was underestimated by more than 50% in the centroid-to-centroid model. Therefore, it is necessary to consider 3D muscle attachment morphometry, especially for muscles with broad attachments, in TMJ musculoskeletal models to precisely quantify the joint mechanical environment critical for understanding TMJ function and mechanobiology.     

Knee muscle moment arms from MRI and from tendon travel.

We tested magnetic resonance imaging (MRI) as a means to collect geometric data for moment arm estimation. A knee specimen in five successive flexion postures was scanned by MRI, while simultaneously tendon positions of loaded muscles were measured (long head of biceps femoris, lateral and medial gastrocnemius, gracilis, rectus femoris, sartorius, semimembranosus, semitendinosus, and tensor fasciae latae). Discrete rotation centres were derived from MRI pictures. Moment arms were estimated as the distances from these centres to the tendons. The ratio of tendon travel over the increment of joint angulation was the alternative, more reliable estimate of the moment arm. An important principal shortcoming of MRI is the impossibility of accounting for force distribution in taut tissue. As a consequence, for some muscles, considerable inaccuracies in moment arm estimation are found in a relatively small range of joint angulation (up to about 30% for the rectus femoris and semimembranosus). For the tensor fasciae latae, the moment arm cannot be estimated by MRI, while the estimate by tendon travel is unreliable owing to the deformability and attachments of the fascia lata.     

Muscle parameters estimation based on biplanar radiography

The evaluation of muscle and joint forces in vivo is still a challenge. Musculo-Skeletal (musculo-skeletal) models are used to compute forces based on movement analysis. Most of them are built from a scaled-generic model based on cadaver measurements, which provides a low level of personalization, or from Magnetic Resonance Images, which provide a personalized model in lying position. This study proposed an original two steps method to access a subject-specific musculo-skeletal model in 30 min, which is based solely on biplanar X-Rays. First, the subject-specific 3D geometry of bones and skin envelopes were reconstructed from biplanar X-Rays radiography. Then, 2200 corresponding control points were identified between a reference model and the subject-specific X-Rays model. Finally, the shape of 21 lower limb muscles was estimated using a non-linear transformation between the control points in order to fit the muscle shape of the reference model to the X-Rays model. Twelfth musculo-skeletal models were reconstructed and compared to their reference. The muscle volume was not accurately estimated with a standard deviation (SD) ranging from 10 to 68%. However, this method provided an accurate estimation the muscle line of action with a SD of the length difference lower than 2% and a positioning error lower than 20 mm. The moment arm was also well estimated with SD lower than 15% for most muscle, which was significantly better than scaled-generic model for most muscle. This method open the way to a quick modeling method for gait analysis based on biplanar radiography.