the online version of this article ( doi:10.1007/s11517 - 012 - 1016 - 3 ) contains supplementary material , which is available to authorized users .
total hip arthroplasty ( tha ) is one of the most frequently performed reconstructive operations with excellent intermediate to long - term results . however , there are still complications which require the initial procedure to be revised , most commonly due to aseptic loosening and dislocation . both aseptic loosening and dislocation
an overly contained cup may lead to impingement between the neck of the femoral component and the rim of the acetabular cup during terminal motion of the hip .
such contact can create wear particles potentially leading to implant loosening [ 22 , 41 , 54 ] .
further motion beyond the impingement point causes subluxation of the femoral head until the joint dislocates [ 21 , 27 , 31 ] .
in contrast , orientating the prosthetic components to maximise range of motion to prevent impingement would mean only partial containment of the hip joint which risks aseptic loosening and joint dislocation whereby the femoral head  slips out of the acetabular cup [ 21 , 52 ] .
yoshimine and ginbayashi   specified five factors that determine the range of motion which a tha can achieve , four of these are associated with prosthetic component orientation : ( 1 ) acetabular cup anteversion , ( 2 ) acetabular cup inclination , ( 3 ) femoral stem version and ( 4 ) femoral component neck axis away from the transverse plane which is dependent upon femoral stem varus  valgus within the femoral canal and femoral component neck - shaft angle [ 22 , 53 , 56 ] .
a further orientation parameter has been defined by renkawitz et al .   which is significant in femoral components with a non - axis symmetric neck .
femoral tilt where the orientation of the femoral neck in the sagittal plane is controlled by where the femoral stem follows the natural anterior bow of the proximal femur .
all of these factors interact to affect the position of the hip primary arc of movement , as illustrated in fig .  1 .
the final factor termed the oscillation angle (  ) , determines the size of the hip primary arc of movement and is a function of the opening angle of the acetabular liner and the femoral head  neck ratio .
hence , a prosthesis with a large oscillation angle to maximise range of motion and is orientated with good femoral head coverage to achieve a stable joint represents the best balance of these factors .
however , increasing the femoral head diameter to maximise the oscillation angle has been shown to increase the risk of femoral neck fracture in hip resurfacing and has been associated with failures in metal - on - metal implants [ 32 , 45 ] . therefore , achieving the correct prosthetic component orientation to achieve both ideal range of motion and secure containment within the constrained prosthetic impingement limits is vital to operative success.fig . 
1effect of oscillation angle and component orientation on range of motion in the sagittal plane . a poorly orientated acetabular cup  required amount of flexion can not be attained within oscillation angle impingement limits .
b correctly orientated acetabular cup  required amount of flexion can be attained within oscillation angle impingement limits effect of oscillation angle and component orientation on range of motion in the sagittal plane . a poorly orientated acetabular cup  required amount of flexion can not be attained within oscillation angle impingement limits .
b correctly orientated acetabular cup  required amount of flexion can be attained within oscillation angle impingement limits determining the boundary within which an impingement free range of motion is required would not only allow surgeons to determine the optimal prosthetic component orientation in tha but also determine the required bone resection to relieve impingement in a native hip with femoroacetabular impingement ( fai ) [ 1 , 49 ] . at present , specifications for an impingement free range motion outcome have been based on limits of pure joint motion of healthy individuals in the coronal , sagittal and transverse planes [ 11 , 42 , 54 , 56 ] , or from measuring joint rotations for specific activities of daily living [ 16 , 19 , 22 , 31 ] .
other studies have used computer tomography ( ct ) scans of healthy bony anatomy to determine patient range of motion requirement [ 23 , 4649 ] .
these studies have acknowledged that range of motion in the native hip is restricted by osseous impingement , soft tissue impingement as well soft tissue restraint .
however , these studies have not been able to determine the restriction in range of motion due to soft tissue restraint .
consequently , the aim of this study is to address this gap in knowledge by determining a range of motion benchmark which can identify motions that are at risk from impingement as well as motions which are limited by soft tissue restraint .
to identify motions that were at risk from impingement as well as those motions which were limited by soft tissue restraint , data were obtained from literature for 15 activities of daily living of healthy male subjects .
these motions were used as the basis to construct a healthy range of motion benchmark which considered range of motion restriction due to both impingement and soft tissue restraint .
range of motion simulation using three - dimensional models constructed from patient ct scans was then used to identify motion restriction due to osseous impingement . comparing the two motion boundaries we were able to identify those motions at risk of impingement and those which were limited by soft tissue restraint .
the 15 activities were  sitting on the floor cross - legged , kneeling with ankles dorsi - flexed and ankles plantar - flexed , level - walking , standing while turning the upper body away , lying supine and rolling over , squatting both with feet flat and on flexed toes , stand - sit - stand from both a normal and a low seat , ascending and descending stairs , standing then bending to retrieve an object from floor , swinging ones leg back and forth , sitting on a normal seat and bending to tie shoe laces , sitting on a normal seat while crossing legs [ 16 , 19 , 20 , 22 , 31 ] . using the joint coordinate system developed by grood and suntay   and adapted for the hip by wu et al . , an anatomical reference frame was constructed with the following axis definitions : x axis 
anterior / posterior axis : abduction / adduction ( a)y axis  superior / inferior axis : internal / external rotation ( r)z axis 
medial / lateral axis : flexion / extension ( f ) the change in joint angles over the movement cycle has led some researchers to divide a manoeuvre into distinct stages , for example with stand  sit  stand from a normal chair  upright , natural , and leaning forward . for each of the 15 activities ,
key points were selected by identifying the motions of maximum flexion / extension ( f ) , abduction / adduction ( a ) and internal / external rotation ( r ) and then recording the corresponding joint angles at each maximum point in the other two anatomical planes .
( 1 ) , where the initial knee centre position was defined as position vector \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				
\begin{document}$$ p = ( 0 , - 1,0 ) $ $ \end{document}. this position vector accounts for when a person stands in the anatomical neutral posture where they are upright and erect on both legs so that the knee centre position lies directly below the hip centre [ 26 , 39 ] .
this means that the proximal pelvic body segment and the distal femoral body segment are initially aligned so their orientation matrix is a 3    3 identity matrix , i .
consequently , the orientation matrix \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ r = \ { [ ( r_{\text{f } } i)r_{\text{a } } ] r_{\text{r } } \ } $ $ \end{document}can be simplified to define the position of the knee centre using eq . ( 1 ) .
these knee centre positions were then used to calculate individual axes of rotation for each of the data points by solving eq .
( 2 ) , where r defines the orientation matrix \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ r = \ { [ ( r_{\text{f } } i)r_{\text{a } } ] r_{\text{r } } \ } $ $ \end{document } and v defines the fixed axis of rotation [ 6 , 15 , 24 ] .
it was found that 70  % of these manoeuvres occurred about a rotation axis within 15 of the transverse plane indicating the dominant motions for the hip joint were flexion / extension and abduction / adduction coupled with smaller amounts of internal / external rotation .
consequently , rotation axes in the transverse plane were able to be used to simulate range of motion using the ct scans of ten patients for comparison against a healthy range of motion benchmark.1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ p = ( [ \sin_{\text{f } } \cdot \cos_{\text{a } } ] x , [ - \cos_{\text{f } } \cdot \cos_{\text{a } } ] y , [ - \sin_{\text{a } } ] z ) $ $ \end{document}2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ qv = [ r - r^{\text{t } } ] v = \left [ { \begin{array}{*{20}c } 0 & { - q_{z } } & { q_{y } } \\ { q_{z } } & 0 & { - q_{x } } \\ { - q_{y } } & { q_{x } } & 0 \\ \end{array } } \right ]
\cdot \left [ { \begin{array}{*{20}c } { v_{x } } \\ { v_{y } } \\ { v_{z } } \\ \end{array } } \right ] = \left [ { \begin{array}{*{20}c } 0 \\ 0 \\ 0 \\ \end{array } } \right ] $ $ \end{document } to construct a healthy range of motion benchmark , experimental data with regard to pure joint motion of flexion / extension , abduction / adduction and internal / external rotation were used ( table  1 ) .
these pure joint motions have been found to provide a good functional outcome for the patient .
therefore , they correlate positively with being able to perform activities of daily living [ 10 , 19 ] . starting again with an initial knee centre position vector \documentclass[12pt]{minimal } 
				
\usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ p = ( 0 , - 1,0 ) $ $ \end{document } , this vector was rotated about axes in the transverse plane \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( t ) $ $ \end{document } defined in eq .
the angle  was rotated about the transverse plane in 10 increments , producing 36 separate rotation axes within this plane .
the point at which impingement occurred for each of the transverse plane rotation axes was plotted on both a two - dimensional and three - dimensional plot , which permitted range of motion to be represented graphically as a continuum ( figs .  2 , 3 ) .
therefore , the constructed benchmark represented the required range of motion to be achieved to perform daily activities without risk of impingement .
visualisation of how this benchmark was constructed is provided in additional file 1.table  1reference pure joint motions of the hip motionangleflexion120extension30abduction45adduction35internal rotation45external rotation45fig . 
3three - dimensional representation of the healthy range of motion benchmark reference pure joint motions of the hip   two - dimensional representation of the healthy range of motion benchmark three - dimensional representation of the healthy range of motion benchmark in addition to defining the shape of the healthy range of motion benchmark its position was also defined . to do this ,
a technique known as moment of inertia analysis was used , where its position was defined using a directional axis ( fig . 
the directional axis represented the normal vector to the best - fit plane constructed from points taken at the edge of the range of motion benchmark .
the directional axis was determined by calculating the centre of mass of the range of motion benchmark , eq .
the distance of each of the boundary edge points away from the calculated centre of mass was then determined to produce a 3    n matrix , a. the dot product \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ a \cdot a^{\text{t } } $ $ \end{document } shown in eq .
( 5 ) was then solved by finding the eigenvector \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( \upsilon ) $ $ \end{document } which maximised the distance to the boundary edge points . both figs . 
2 and 3 represent range of motion restriction in a healthy hip due to both impingement and soft tissue restraint .
the ct method was then used to highlight the impingement risk within the range of motion benchmark.3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ t = ( \sin \alpha , 0,\cos \alpha ) $ $ \end{document}4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( \bar{x},\bar{y},\bar{z } ) = \frac{{\left ( { \sum\nolimits_{n = 1}^{n } { ( x ) } , \sum\nolimits_{n = 1}^{n } { ( y ) } , \sum\nolimits_{n = 1}^{n } { ( z ) } } \right)}}{n } $ $ \end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \text{if:}}\quad a_{i } = \sum { x_{i } - \bar{x } } , \quad b_{i } = \sum { y_{i } - \bar{y } } , \quad c_{i } = \sum { z_{i } - \bar{z } } $ $ \end{document}5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \text{then}}:\quad a \cdot a^{t } = \left [ { \begin{array}{*{20}c } { \sum { a_{i}^{2 } } } & { \sum { a_{i } b_{i } } } & { \sum { a_{i } c_{i } } } \\ { \sum { b_{i } a_{i } } } & { \sum { b_{i}^{2 } } } & { \sum { b_{i } c_{i } } } \\ { \sum { c_{i } a_{i } } } & { \sum { c_{i } b_{i } } } & { \sum { c_{i}^{2 } } } \\ \end{array } } \right ] $ $ \end{document } for the assessment of range of motion restriction in the hip joint due to osseous impingement , ten ct scans were taken of patients in the supine position .
the ct scans used in this study were available on the university hospitals coventry and warwickshire ( uhcw ) image library .
all patients were scanned for clear clinical reasons and informed consent to use their images for teaching and research purposes was obtained .
the ct scans used in this study were fully anonymised and performed in accordance with the institution s ethical guidelines and with the declaration of helsinki .
the subjects were all male with a mean age of 65.2  years ( 4981 ) and exhibited no evidence of osteoarthritis or abnormal morphology .
the scans were acquired on a general electric lightspeed ct scanner with a slice thickness of 1.25  mm using a soft tissue algorithm , encompassing the complete anatomy of the pelvis and femur .
the ct images in dicom format were imported into the imagej image - processing software ( http://rsbweb.nih.gov/ij/ ) .
the dicom scan slices were converted to binary images and a threshold was applied so that only matter with the same density as bone remained .
each slice was then manually cleaned , removing any non - bone material and filling the gaps in the pelvic or femoral trace .
each cleaned image stack was then imported into the simpleware scanip ( simpleware ltd . ,
pixel spacing was applied to smooth the inconsistencies between slices and a 3d model mesh was then generated for the pelvic and femoral masks .
these three - dimensional models were then imported into the rhino 4.0 nurbs modelling package for use in the motion simulation experiment .
prior to the range of motion simulation both the pelvis and femur had to be orientated correctly in three - dimensional space .
the pelvic coordinate frame was defined using the landmarks of the transverse pelvic plane ( tpp ) .
lateral axis was defined as a line running parallel to the two anterior superior iliac spines ( asis ) running in the positive direction from left to right with the origin at the hip joint centre .
the hip joint centre was defined as the centre of a best - fit sphere of the femoral head .
the transverse plane was defined as a plane containing the two asis and the mid - point of the two posterior superior iliac spines ( psis ) . a line perpendicular to the transverse plane with the origin at the hip joint centre defined the superior - inferior direction .
the coordinate system of the femur was defined according to the standard defined by murphy et al . .
the superior  inferior axis was defined as running in the positive direction from the knee centre to the hip joint centre .
the knee centre was defined by the mid - point of the two femoral epicondyles .
the coronal plane was defined as a plane containing the hip joint centre and a line parallel to the posterior aspect of the femoral condyles located at the knee centre .
the anteroposterior axis was constructed perpendicular to the coronal plane located at the hip joint centre , and the medial 
the femoral 3d model was then aligned so that its axes were coincident with the coordinate frame of the pelvis . the constructed coordinate frame
has previously been found to define the neutral rotation of the femur and when aligned with the pelvic coordinate frame forms an orthogonal joint coordinate frame .
consequently , the subsequent range of motion simulation could be directly compared against the healthy range of motion benchmark . with the pelvis and femur aligned ,
the rhino vbscript language was then used to rotate the femur about axes in the transverse plane defined in eq .
impingement was deemed to have occurred when the femoral triangle mesh intersected with the rim of the pelvic acetabulum .
once impingement had occurred , no further motion was considered possible . the impingement angle for each rotation axis
was plotted on both two - dimensional and three - dimensional plots for analysis , in addition to the ct directional axis .
the average range of motion for the ten simulations were plotted and compared with the healthy range of motion benchmark . to define the type of motion restriction  at any point where
there was impinged motion of less than 5 between healthy range of motion benchmark and the ct benchmark , then this was defined as osseous impingement , in accordance with the findings by tannast et al . . a difference between the healthy and the ct benchmarks of less than 10
was defined as soft tissue impingement and any difference greater than 10 was defined as soft tissue restraint .
to construct a healthy range of motion benchmark , experimental data with regard to pure joint motion of flexion / extension , abduction / adduction and internal / external rotation were used ( table  1 ) .
these pure joint motions have been found to provide a good functional outcome for the patient .
therefore , they correlate positively with being able to perform activities of daily living [ 10 , 19 ] . starting again with an initial knee centre position vector \documentclass[12pt]{minimal } 
				
\usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ p = ( 0 , - 1,0 ) $ $ \end{document } , this vector was rotated about axes in the transverse plane \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( t ) $ $ \end{document } defined in eq .
the angle  was rotated about the transverse plane in 10 increments , producing 36 separate rotation axes within this plane .
the point at which impingement occurred for each of the transverse plane rotation axes was plotted on both a two - dimensional and three - dimensional plot , which permitted range of motion to be represented graphically as a continuum ( figs .  2 , 3 ) .
therefore , the constructed benchmark represented the required range of motion to be achieved to perform daily activities without risk of impingement .
visualisation of how this benchmark was constructed is provided in additional file 1.table  1reference pure joint motions of the hip motionangleflexion120extension30abduction45adduction35internal rotation45external rotation45fig . 
3three - dimensional representation of the healthy range of motion benchmark reference pure joint motions of the hip   two - dimensional representation of the healthy range of motion benchmark three - dimensional representation of the healthy range of motion benchmark in addition to defining the shape of the healthy range of motion benchmark its position was also defined . to do this ,
a technique known as moment of inertia analysis was used , where its position was defined using a directional axis ( fig . 
the directional axis represented the normal vector to the best - fit plane constructed from points taken at the edge of the range of motion benchmark .
the directional axis was determined by calculating the centre of mass of the range of motion benchmark , eq .
the distance of each of the boundary edge points away from the calculated centre of mass was then determined to produce a 3    n matrix , a. the dot product \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ a \cdot a^{\text{t } } $ $ \end{document } shown in eq .
( 5 ) was then solved by finding the eigenvector \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( \upsilon ) $ $ \end{document } which maximised the distance to the boundary edge points . both figs . 
2 and 3 represent range of motion restriction in a healthy hip due to both impingement and soft tissue restraint .
the ct method was then used to highlight the impingement risk within the range of motion benchmark.3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ t = ( \sin \alpha , 0,\cos \alpha ) $ $ \end{document}4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( \bar{x},\bar{y},\bar{z } ) = \frac{{\left ( { \sum\nolimits_{n = 1}^{n } { ( x ) } , \sum\nolimits_{n = 1}^{n } { ( y ) } , \sum\nolimits_{n = 1}^{n } { ( z ) } } \right)}}{n } $ $ \end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \text{if:}}\quad a_{i } = \sum { x_{i } - \bar{x } } , \quad b_{i } = \sum { y_{i } - \bar{y } } , \quad c_{i } = \sum { z_{i } - \bar{z } } $ $ \end{document}5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \text{then}}:\quad a \cdot a^{t } = \left [ { \begin{array}{*{20}c } { \sum { a_{i}^{2 } } } & { \sum { a_{i } b_{i } } } & { \sum { a_{i } c_{i } } } \\ { \sum { b_{i } a_{i } } } & { \sum { b_{i}^{2 } } } & { \sum { b_{i } c_{i } } } \\ { \sum { c_{i } a_{i } } } & { \sum { c_{i } b_{i } } } & { \sum { c_{i}^{2 } } } \\
for the assessment of range of motion restriction in the hip joint due to osseous impingement , ten ct scans were taken of patients in the supine position .
the ct scans used in this study were available on the university hospitals coventry and warwickshire ( uhcw ) image library .
all patients were scanned for clear clinical reasons and informed consent to use their images for teaching and research purposes was obtained .
the ct scans used in this study were fully anonymised and performed in accordance with the institution s ethical guidelines and with the declaration of helsinki .
the subjects were all male with a mean age of 65.2  years ( 4981 ) and exhibited no evidence of osteoarthritis or abnormal morphology .
the scans were acquired on a general electric lightspeed ct scanner with a slice thickness of 1.25  mm using a soft tissue algorithm , encompassing the complete anatomy of the pelvis and femur .
the ct images in dicom format were imported into the imagej image - processing software ( http://rsbweb.nih.gov/ij/ ) .
the dicom scan slices were converted to binary images and a threshold was applied so that only matter with the same density as bone remained .
each slice was then manually cleaned , removing any non - bone material and filling the gaps in the pelvic or femoral trace . each cleaned image stack
pixel spacing was applied to smooth the inconsistencies between slices and a 3d model mesh was then generated for the pelvic and femoral masks .
these three - dimensional models were then imported into the rhino 4.0 nurbs modelling package for use in the motion simulation experiment .
prior to the range of motion simulation both the pelvis and femur had to be orientated correctly in three - dimensional space .
the pelvic coordinate frame was defined using the landmarks of the transverse pelvic plane ( tpp ) .
lateral axis was defined as a line running parallel to the two anterior superior iliac spines ( asis ) running in the positive direction from left to right with the origin at the hip joint centre .
the hip joint centre was defined as the centre of a best - fit sphere of the femoral head .
the transverse plane was defined as a plane containing the two asis and the mid - point of the two posterior superior iliac spines ( psis ) . a line perpendicular to the transverse plane with the origin at the hip joint centre defined the superior - inferior direction .
the coordinate system of the femur was defined according to the standard defined by murphy et al . .
the superior  inferior axis was defined as running in the positive direction from the knee centre to the hip joint centre .
the knee centre was defined by the mid - point of the two femoral epicondyles .
the coronal plane was defined as a plane containing the hip joint centre and a line parallel to the posterior aspect of the femoral condyles located at the knee centre .
the anteroposterior axis was constructed perpendicular to the coronal plane located at the hip joint centre , and the medial 
the femoral 3d model was then aligned so that its axes were coincident with the coordinate frame of the pelvis .
the constructed coordinate frame has previously been found to define the neutral rotation of the femur and when aligned with the pelvic coordinate frame forms an orthogonal joint coordinate frame . consequently , the subsequent range of motion simulation could be directly compared against the healthy range of motion benchmark . with the pelvis and femur aligned ,
the rhino vbscript language was then used to rotate the femur about axes in the transverse plane defined in eq .
impingement was deemed to have occurred when the femoral triangle mesh intersected with the rim of the pelvic acetabulum .
the impingement angle for each rotation axis was plotted on both two - dimensional and three - dimensional plots for analysis , in addition to the ct directional axis .
the average range of motion for the ten simulations were plotted and compared with the healthy range of motion benchmark . to define the type of motion restriction  at any point where
there was impinged motion of less than 5 between healthy range of motion benchmark and the ct benchmark , then this was defined as osseous impingement , in accordance with the findings by tannast et al . . a difference between the healthy and the ct benchmarks of less than 10
was defined as soft tissue impingement and any difference greater than 10 was defined as soft tissue restraint .
a two - dimensional plot of the ct range of motion ( purple ) compared to the healthy range of motion benchmark ( gold ) is shown in fig .  4 .
this plot shows that there is osseous impingement present for the motions of adduction , abduction combined with flexion and pure flexion .
impinged motion has been shown in red where the ct range of motion does not encompass the healthy range of motion benchmark .
there is a significantly larger ct range of motion when compared with the healthy benchmark with regard to motion in pure extension and adduction combined with flexion , which represented motion limitation due to soft tissue restraint . in all other areas ,
the ct range of motion is slightly greater than the healthy range of motion benchmark ( 510 ) representing soft tissue impingement .
figure  5 details the three - dimensional comparison of the two experimental results which is also presented in additional file 1 .
the directional axis of the ct range of motion ( purple ) is less elevated than the healthy range of motion benchmark axis ( red ) and differs by a three - dimensional angle of 15.4. the two - dimensional angle between these axes in the transverse plane is 3.1 , which disregards the elevation difference between these two axes due to motion in extension .
the mean and the standard deviation for the pure joint motions of the 10 ct scans were flexion 120 (   =  10.2 ) , extension 77 (   =  20.1 ) , abduction 55 (   =  9.9 ) and adduction 33 (   =  8.8).fig . 
4two - dimensional representation of ct range of motion comparison with the clinical range of motion benchmarkfig . 
5three - dimensional representation of ct range of motion comparison with the clinical range of motion benchmark two - dimensional representation of ct range of motion comparison with the clinical range of motion benchmark three - dimensional representation of ct range of motion comparison with the clinical range of motion benchmark
studies have attempted to quantify range of motion requirement of the native hip through clinical measurements , gait analysis and ct dynamic simulation .
these studies have acknowledged that range of motion in the constrained hip is restricted by osseous impingement , soft tissue impingement of the labrum and capsule as well soft tissue restraint [ 23 , 25 , 4749 ] .
no study has yet attempted to fully quantify which motions are associated with these separate causes of motion restriction .
this study has provided a comparison of two different methods of analysing range of motion  measurement from previous studies which measured the activities and pure joint motion of healthy individuals which includes all modes of restriction within a person s motion , and secondly ct dynamic simulation which can identify range of motion restriction due to osseous impingement .
consequently , providing a comparison of the results we were able to distinguish which mode of restriction was limiting range of motion for a given manoeuvre .
this was combined with the knowledge that in ct simulations , range of motion restriction slightly over estimates the required range of motion due to the absence of soft tissue by 5 .
firstly , it has been regarded that measurements of healthy individuals between the ages 2070 provide the most stable and realistic sample from which to base a range of motion benchmark [ 18 , 38 ] .
the ct analysis of hip joint range of motion used patients above this 70  year old threshold . however
, this reduction in range of motion is associated with neuro - muscular function rather than morphological changes within the joint and does not affect the result of the ct range of motion experiment [ 28 , 33 ] .
a second limitation was that the two experimental methodologies used data acquired from different subjects .
a cadaver study simulating the range motions presented in combination with ct measurement would make the study findings more objective with regard to the effect that soft tissue restraint has upon range of motion [ 40 , 44 ] .
how much the healthy range of motion benchmark can be generalised needs to be considered .
this is because the range of motion benchmark has used data from male subjects in its construction .
considering ethnicity , although there may be a greater demand from asian and middle eastern cultures to perform high excursion manoeuvres such as kneeling and squatting [ 16 , 29 ] .
these manoeuvres have been incorporated within the range of motion benchmark and do not exceed the pure joint motion values derived from measurements of mainly european or american subjects .
the need for having an alternative range of motion benchmark based on gender is unclear .
data from level - walking studies show that females have in the region of 45 greater motion [ 3 , 4 ] and similar results were found in a limited number of studies measuring higher demand sporting activities [ 12 , 35 ] .
it is unknown whether this difference transfers across all activities to indicate whether females have greater joint mobility or whether , similar to age , joint excursion is dependent upon other factors such as neurological and muscle interactions .
therefore , it is difficult to assess whether a separate range of motion benchmark is required for a female population .
analysis of the results with regard to pure flexion , pure adduction and flexion combined with adduction shows that there was impinged motion in these areas ( fig .  4 ) .
impingement in pure adduction was due to collision between the lesser trochanter of the femur with the pelvis .
in contrast , impingement in pure flexion and flexion combined with adduction was due to collision between the femoral neck and the acetabular rim in the anterosuperior zone .
this is recognised as being the zone in which hip damage occurs , leading to the onset of osteoarthritis [ 8 , 47 ] . consequently ,
motion in this area should be maximised in tha beyond the impingement point to ensure that the femoral neck does not contact with the rim of the acetabular cup .
the impinged motion of 90 flexion combined with adduction was also coupled with internal rotation of 33 which the range of motion comparison does not visualise .
it is regarded that hip pain due to impingement can be replicated by internally rotating the femur at 90 of flexion [ 7 , 34 ] .
if motion in this position is pain free with 30 of internal rotation then this represents an acceptable range of motion for the hip joint .
consequently , the healthy range of motion benchmark incorporates this amount of normal internal rotation in motions with 90 of flexion .
however , knowledge from these previous studies shows that any coupled motion above 90 flexion is a risk area with regard to osseous impingement and should be acknowledged as such [ 7 , 34 ] .
there were a number of motions where the ct range of motion fell within the 510 soft tissue impingement zone defined in sect .
figure  4 shows that abduction is a risk area for soft tissue impingement and is congruent with the findings of kubiak - langer et al .   and tannast et al .  
therefore , in tha it should be ensured that component positioning ensures that the healthy range of motion benchmark is attained , which for pure adduction is 45 , as this signifies the contact point for soft tissue impingement .
there were two areas where the ct patient range of motion was significantly larger than the healthy range of motion benchmark .
it has been found that motion in extension is not limited by osseous impingement ; rather it is limited by soft tissue contracture or limitation in secondary joint motion .
consequently , the patient range of motion in this area can be regarded as clinically non - relevant .
the healthy range of motion directional axis is more elevated than the ct directional axis .
if motion in extension is discounted by measuring the two - dimensional angle between these axes in the transverse plane , then the two directional axes align well with only a difference of 3.1. this demonstrates a close correlation between healthy and ct range of motion , discounting the clinically non - relevant motion in extension and validates the constructed healthy range of motion benchmark .
the extra range of motion found in the area of adduction combined with flexion , when comparing the ct range of motion with the healthy benchmark is a new finding .
it is hypothesised two reasons could cause this extreme deviation between the ct range of motion and the healthy benchmark .
firstly , it is not possible to measure true geometrical adduction as a medial rotation in the coronal plane .
therefore , measurement of hip joint adduction follows a diagonal motion as the adducted leg is moved in front of the stationary leg . therefore , the construction of the healthy range of motion benchmark should have considered the pure joint motion adduction benchmark value of 35 in this diagonal plane rather than the coronal plane .
secondly , the extra motion in adduction combined with flexion as exhibited by the ct range of motion may not be limited by osseous impingement .
this is because in the ct method , the motion of adduction combined with flexion took the femur into the acetabular notch , permitting extra motion .
it is more likely that motion is limited in this area by tension in the adductor muscles .
this study has used two methods for measuring range of motion in the native hip to determine which factors restrict motion .
the results show that motion in pure flexion and flexion combined with adduction are at risk of osseous impingement .
these motions represent where the maximum likely damage will occur in femoroacetabular impingement or are at most risk of prosthetic impingement post - tha .
the study has also shown that motions in extension and adduction combined with flexion are limited by soft tissue restraint , while motions such as pure abduction are a risk for soft tissue impingement .
6 and colour coded to highlight apparent risk  osseous impingement ( red ) , soft tissue impingement ( orange ) and soft tissue restraint ( green ) .
recognising where the impingement risk is within a person s range of motion will provide further guidance with regard to where motion is required to be maximised when performing tha or treating femoroacetabular impingement.fig . 
6illustration of range of motion benchmark with impingement zones illustration of range of motion benchmark with impingement zones
supplementary material 1 ( mp4 12722 kb ) supplementary material 1 ( mp4 12722 kb )
 in total hip arthroplasty , determining the impingement free range of motion requirement is a complex task . this is because in the native hip , motion is restricted by both impingement as well as soft tissue restraint . the aim of this study is to determine a range of motion benchmark which can identify motions which are at risk from impingement and those which are constrained due to soft tissue . 
 two experimental methodologies were used to determine motions which were limited by impingement and those motions which were limited by both impingement and soft tissue restraint . by comparing these two experimental results , 
 motions which were limited by impingement were able to be separated from those motions which were limited by soft tissue restraint . 
 the results show motions in extension as well as flexion combined with adduction are limited by soft tissue restraint . 
 motions in flexion , flexion combined with abduction and adduction are at risk from osseous impingement . 
 consequently , these motions represent where the maximum likely damage will occur in femoroacetabular impingement or at most risk of prosthetic impingement in total hip arthroplasty.electronic supplementary materialthe online version of this article ( doi:10.1007/s11517 - 012 - 1016 - 3 ) contains supplementary material , which is available to authorized users . 
Electronic supplementary material
Introduction
Methods
Healthy range of motion benchmark
CT range of motion simulation
Results
Discussion
Electronic supplementary material