Deep learning-based automated angle measurement for flatfoot diagnosis in weight-bearing lateral radiographs.

This study aimed to develop and evaluate a deep learning-based system for the automatic measurement of angles (specifically, Meary's angle and calcaneal pitch) in weight-bearing lateral radiographs of the foot for flatfoot diagnosis. We utilized 3960 lateral radiographs, either from the left or right foot, sourced from a pool of 4000 patients to construct and evaluate a deep learning-based model. These radiographs were captured between June and November 2021, and patients who had undergone total ankle replacement surgery or ankle arthrodesis surgery were excluded. Various methods, including correlation analysis, Bland-Altman plots, and paired T-tests, were employed to assess the concordance between the angles automatically measured using the system and those assessed by clinical experts. The evaluation dataset comprised 150 weight-bearing radiographs from 150 patients. In all test cases, the angles automatically computed using the deep learning-based system were in good agreement with the reference standards (Meary's angle: Pearson correlation coefficient (PCC) = 0.964, intraclass correlation coefficient (ICC) = 0.963, concordance correlation coefficient (CCC) = 0.963, p-value = 0.632, mean absolute error (MAE) = 1.59°; calcaneal pitch: PCC = 0.988, ICC = 0.987, CCC = 0.987, p-value = 0.055, MAE = 0.63°). The average time required for angle measurement using only the CPU to execute the deep learning-based system was 11 ± 1 s. The deep learning-based automatic angle measurement system, a tool for diagnosing flatfoot, demonstrated comparable accuracy and reliability with the results obtained by medical professionals for patients without internal fixation devices.

A narrative review of the measurement methods for biomechanical properties of plantar soft tissue in patients with diabetic foot.

This article provides an overview of the development history and advantages and disadvantages of measurement methods for soft tissue properties of the plantar foot. The measurement of soft tissue properties is essential for understanding the biomechanical characteristics and function of the foot, as well as for designing and evaluating orthotic devices and footwear. Various methods have been developed to measure the properties of plantar soft tissues, including ultrasound imaging, indentation testing, magnetic resonance elastography, and shear wave elastography. Each method has its own strengths and limitations, and choosing the most appropriate method depends on the specific research or clinical objectives. This review aims to assist researchers and clinicians in selecting the most suitable measurement method for their specific needs.

The Importance of Preconditioning for the Sonographic Assessment of Plantar Fascia Thickness and Shear Wave Velocity.

Plantar fasciopathy is a very common musculoskeletal complaint that leads to reduced physical activity and undermines the quality of life of patients. It is associated with changes in plantar fascia structure and biomechanics which are most often observed between the tissue's middle portion and the calcaneal insertion. Sonographic measurements of thickness and shear wave (SW) elastography are useful tools for detecting such changes and guide clinical decision making. However, their accuracy can be compromised by variability in the tissue's loading history. This study investigates the effect of loading history on plantar fascia measurements to conclude whether mitigation measures are needed for more accurate diagnosis. The plantar fasciae of 29 healthy participants were imaged at baseline and after different clinically relevant loading scenarios. The average (±standard deviation) SW velocity was 6.5 m/s (±1.5 m/s) and it significantly increased with loading. Indicatively, five minutes walking increased SW velocity by 14% (95% CI: -1.192, -0.298, t(27), p = 0.005). Thickness between the calcaneal insertion and the middle of the plantar fascia did not change with the tissues' loading history. These findings suggest that preconditioning protocols are crucial for accurate SW elastography assessments of plantar fasciae and have wider implications for the diagnosis and management of plantar fasciopathy.

The impact of diabetes mellitus on foot perfusion measured by ICG NIR fluorescence imaging.

INTRODUCTION: Diabetes Mellitus (DM) is a common chronic disease, affecting 435 million people globally. Impaired vasculature in DM patients leads to complications like lower extremity arterial disease (LEAD) and foot ulcers, often resulting in amputations. DM causes additional peripheral neuropathy leading to multifactorial wound problems. Current diagnostics often deem unreliable, but Near-Infrared Fluorescence with Indocyanine Green (ICG NIR) can be used to assess the foot perfusion. Therefore, this study explores DM's impact on foot perfusion using ICG NIR. METHODS: Baseline ICG NIR fluorescence imaging was performed in LEAD patients with and without DM. Ten perfusion parameters were extracted and analyzed to assess differences in perfusion patterns. RESULTS: Among 109 patients (122 limbs) of the included patients, 32.8 % had DM. Six of ten perfusion parameters, mainly inflow-related, differed significantly between DM and non-DM patients (p-values 0.007-0.039). Fontaine stage 4 DM patients had the highest in- and outflow values, with seven parameters significantly higher (p-values 0.004-0.035). CONCLUSION: DM is associated with increased in- and outflow parameters. Patients with- and without DM should not be compared directly due to different vascular pathophysiology and multifactorial wound problems in DM patients. Quantified ICG NIR fluorescence imaging offers additional insight into the effect of DM on foot perfusion.

Pediatric Foot: Development, Variants, and Related Pathology.

Pediatric foot development throughout childhood and adolescence can present a diagnostic dilemma for radiologists because imaging appearances may be confused with pathology. Understanding pediatric foot development and anatomical variants, such as accessory ossification centers, is essential to interpret musculoskeletal imaging in children correctly, particularly because many of these variants are incidental but others can be symptomatic. We first briefly review foot embryology. After describing common accessory ossification centers of the foot, we explain the different patterns of foot maturation with attention to irregular ossification and bone marrow development. Common pediatric foot variants and pathology are described, such as tarsal coalitions and fifth metatarsal base fractures. We also discuss pediatric foot alignment and various childhood foot alignment deformities.

Ex vivo diffusion tensor imaging of the plantar nerves is feasible in the horse.

OBJECTIVE: The objective of this study was to optimize an MRI-based diffusion tensor imaging (DTI) protocol for imaging the plantar nerves at the level of the tarsus in normal equine limbs. SAMPLE: 12 pelvic cadaver limbs from horses without evidence of proximal suspensory pathology were imaged with a 3T MRI system. METHODS: For diffusion-weighted imaging, b values of 600, 800, and 1,000 s/mm2 were tested. Data were processed with DSI Studio. Cross-sectional areas of the medial and lateral plantar nerve along the plantar tarsus were recorded. The length and number of fiber tracts, signal-to-noise ratio, and DTI variables were recorded. RESULTS: At the level of interest, the mean cross-sectional areas of the plantar nerves ranged from 5.03 to 7.42 mm2. The DTI maps consistently generated tracts in the region of the lateral and medial plantar nerves with DTI values in the range of values reported for peripheral nerves in humans. Our findings demonstrate that DTI of the medial and lateral plantar nerves can be performed successfully and used to generate quantitative parameters including fractional anisotropy and mean, axial, and radial diffusivity. CLINICAL RELEVANCE: Quantitative data generated with this imaging technique can be used to noninvasively characterize the microstructural integrity of neural tissue with possible applications in the evaluation of pathologic changes to the plantar tarsal and metatarsal nerves of horses with proximal suspensory desmopathy.

A finite element model of human hindfoot and its application in supramalleolar osteotomy.

BACKGROUND: The majority of the ankle osteoarthritis cases are posttraumatic and affect younger patients with a longer projected life span. Hence, joint-preserving surgery, such as supramalleolar osteotomy becomes popular among young patients, especially those with asymmetric arthritis due to alignment deformities. However, there is a lack of biomechanical studies on postoperative evaluation of stress at ankle joints. We aimed to construct a verifiable finite element model of the human hindfoot, and to explore the effect of different osteotomy parameters on the treatment of varus ankle arthritis. METHODS: The bones of the hindfoot are reconstructed using normal CT tomography data from healthy volunteers, while the cartilages and ligaments are determined from the literature. The finite element calculation results are compared with the weight-bearing CT (WBCT) data to validate the model. By setting different model parameters, such as the osteotomy height (L) and the osteotomy distraction distance (h), the effects of different surgical parameters on the contact stress of the ankle joint surface are compared. FINDINGS: The alignment and the deformation of hindfoot bones as determined by the finite element analysis aligns closely with the data obtained from WBCT. The maximum contact stress of the ankle joint surface calculated by this model increases with the increase of the varus angle. The maximum contact stresses as a function of the L and h of the ankle joint surface are determined. INTERPRETATION: The relationship between surgical parameters and stress at the ankle joint in our study could further help guiding the planning of the supramalleolar osteotomy according to the varus/valgus alignment of the patients.

Radiological measurements of the feet of normal Straight Egyptian Arabian horses in Qatar.

The Arabian horse has been identified as carrying a risk locus for equine metabolic syndrome, predisposing this breed to development of laminitis. Radigraphy of the equine foot is widely considered the main diagnostic imaging technique for evaluation of the laminitic horse. Knowledge of 'normal' breed values allows assessment of the degree and severity of radiological changes associated with laminitis. The objective of this study was to investigate the normal values for radiological measurements of the feet of the Straight Egyptian Arabian horse in Qatar. The design was a clinical prospective study. Radiographs of the fore and hind feet of 10 clinically normal adult Straight Egyptian Arabian horses were taken. On the lateromedial views, 17 measurements were taken (13 distances and four angles). On the horizontal dorsopalmar/plantar views, two measurements were taken. On the dorsal 45 degree, proximo-palmarodistal oblique projections, four measurements were taken. Normal reference ranges were reported for radiological measurements of the feet of the Straight Egyptian Arabian horse. Several variables showed significant differences between fore and hind feet, including hoof angle, distal wall thickness, and two proximal inner layer measurements (p < 0.05). In addition, the Straight Egyptian Arabian horse was found to have a number of measurements which varied from previously published reports. The results reported within provide a useful reference for normal radiographic measurements of the Straight Egyptian Arabian horse with relevance for laminitis.

Validation of MRI assessment of foot perfusion for improving treatment of patients with peripheral artery disease.

INTRODUCTION: Information on tissue perfusion in the foot is important when treating patients with chronic limb-threatening ischemia. This study aims to test the reliability of different magnetic resonance sequences when measuring perfusion in the foot. METHODS: Sixteen healthy volunteers had their right foot scanned in a test/retest study with six different magnetic resonance sequences (BOLD, multi-echo gradient echo (mGRE), 2D and 3D pCASL, PASL FAIR, and DWI with intravoxel incoherent motion (IVIM) with quantitative measurements of perfusion. For five sequences, cuff-induced ischemia followed by a hyperactive response was measured. Images of the feet were segmented into angiosomes and perfusion data were extracted from the five angiosomes. RESULTS: BOLD, PASL FAIR, mGRE, and DWI with IVIM had low mean differences between the first and second scans, while the results of 2D and 3D pCASL had the highest differences. Based on a paired t-test, BOLD, and FAIR were able to distinguish between perfusion and no perfusion in all angiosomes with p-values below 0.01. This was not the case with 2D and 3D pCASL with p-values above 0.05 in all angiosomes. The mGRE could not distinguish between perfusion and no perfusion in the lateral side of the foot. CONCLUSION: BOLD, mGRE, pASL FAIR, and DWI with IVIM seem to give more robust results compared to 2D and 3D pCASL. Further studies on patients with peripheral artery disease should explore if the sequences can have clinical relevance when assessing tissue ischemia and results of revascularization. IMPLICATIONS FOR PRACTICE: This study provides knowledge that could be used to improve the diagnosis of patient with chronic limb-threatening ischemia to explore tissue perfusion.

Hindfoot kinematics and kinetics - A combined in vivo and in silico analysis approach.

BACKGROUND: The complex anatomical structure of the foot-ankle imposes challenges to accurately quantify detailed hindfoot kinematics and estimate musculoskeletal loading parameters. Most systems used to capture or estimate dynamic joint function oversimplify the anatomical structure by reducing its complexity. RESEARCH QUESTION: Can four dimensional computed tomography (4D CT) imaging in combination with an innovative foot manipulator capture in vivo hindfoot kinematics during a simulated stance phase of walking and can talocrural and subtalar articular joint mechanics be estimated based on a detailed in silico musculoskeletal foot-ankle model. METHODS: A foot manipulator imposed plantar/dorsiflexion and inversion/eversion representing a healthy stance phase of gait in 12 healthy participants while simultaneously acquiring 4D CT images. Participant-specific 3D hindfoot rotations and translations were calculated based on bone-specific anatomical coordinate systems. Articular cartilage contact area and contact pressure of the talocrural and subtalar joints were estimated using an extended foot-ankle model updated with an elastic foundation contact model upon prescribing the participant-specific rotations measured in the 4D CT measurement. RESULTS: Plantar/dorsiflexion predominantly occurred at the talocrural joint (RoM 15.9±3.9°), while inversion/eversion (RoM 5.9±3.9°) occurred mostly at the subtalar joint, with the contact area being larger at the subtalar than at the talocrural joint. Contact pressure was evenly distributed between the talocrural and subtalar joint at the beginning of the simulated stance phase but was then redistributed from the talocrural to the subtalar joint with increasing dorsiflexion. SIGNIFICANCE: In a clinical case study, the healthy participants were compared with four patients after surgically treaded intra-articular calcaneal fracture. The proposed workflow was able to detect small but meaningful differences in hindfoot kinematics and kinetics, indicative of remaining hindfoot pathomechanics that may influence the onset and progression of degenerative joint diseases.

Recommendation of minimal distal tibial length for long axis coordinate system definitions.

Accurate anatomical coordinate systems for the foot and ankle are critical for interpreting their complex biomechanics. The tibial superior-inferior axis is crucial for analyzing joint kinematics, influencing bone motion analysis during gait using CT imaging and biplane fluoroscopy. However, the lack of consensus on how to define the tibial axis has led to variability in research, hindering generalizability. Even as advanced imaging techniques evolve, including biplane fluoroscopy and weightbearing CT, there exist limitations to imaging the entire foot together with the full length of the tibia. These limitations highlight the need to refine axis definitions. This study investigated various superior-inferior axes using multiple distal tibia lengths to determine the minimal field of view for representing the full tibia long-axis. Twenty human cadaver tibias were imaged and segmented to generate 3D bone models. Axes were calculated based on coordinate definitions that required user manual input, and a gold standard mean superior-inferior axis was calculated based on the population's principal component analysis axis. Four manually calculated superior-inferior tibial axes groups were established based on landmarks and geometric fittings. Statistical analysis revealed that geometrically fitting a cylinder 1.5 times the mediolateral tibial width, starting 5 cm above the tibial plafond, yielded the smallest angular deviation from the gold standard. From these findings, we recommend a minimum field of view that includes 1.5 times the mediolateral tibial width, starting 5 cm above the tibial plafond for tibial long-axis definitions. Implementing these findings will help improve foot and ankle research generalizability and impact clinical decisions.

Use of computed tomography to identify muscle quality subgroups, spatial mapping, and preliminary relationships to function in those with diabetic peripheral neuropathy.

BACKGROUND: Decreased muscle volume and increased muscle-associated adipose tissue (MAAT, sum of intra and inter-muscular adipose tissue) of the foot intrinsic muscle compartment are associated with deformity, decreased function, and increased risk of ulceration and amputation in those with diabetic peripheral neuropathy (DPN). RESEARCH QUESTION: What is the muscle quality (normal, abnormal muscle, and adipose volumes) of the DPN foot intrinsic compartment, how does it change over time, and is muscle quality related to gait and foot function? METHODS: Computed tomography was performed on the intrinsic foot muscle compartment of 45 subjects with DPN (mean age: 67.2 ± 6.4 years) at baseline and 3.6 years. Images were processed to obtain volumes of MAAT, highly abnormal, mildly abnormal, and normal muscle. For each category, annual rates of change were calculated. Paired t-tests compared baseline and follow-up. Foot function during gait was assessed using 3D motion analysis and the Foot and Ankle Ability Measure. Correlations between muscle compartment and foot function during gait were analyzed using Pearson's correlations. RESULTS: Total muscle volume decreased, driven by a loss of normal muscle and mildly abnormal muscle (p<0.05). MAAT and the adipose-muscle ratio increased. At baseline, 51.5% of the compartment was abnormal muscle or MAAT, increasing to 55.0% at follow-up. Decreased total muscle volume correlated with greater midfoot collapse during gait (r = -0.40, p = 0.02). Greater volumes of highly abnormal muscle correlated with a lower FAAM score (r = -0.33, p = 0.03). SIGNIFICANCE: Muscle volume loss may progress in parallel with MAAT accumulation, impacting contractile performance in individuals with DPN. Only 48.5% of the DPN intrinsic foot muscle compartment consists of normal muscle and greater abnormal muscle is associated with worse foot function. These changes identify an important target for rehabilitative intervention to slow or prevent muscle deterioration and poor foot outcomes.

Impairments in peroneal muscle size and activation in individuals with patellofemoral pain in weight-bearing position.

BACKGROUND: Patellofemoral pain (PFP) is characterized by chronic pain in the anterior aspect of the knee during loading activities. Many studies investigating muscle morphology changes for individuals with PFP focus on the proximal joints, however, few studies have investigated muscles of the foot and ankle complex. This study aimed to explore the differences in peroneal muscle size and activation between individuals with PFP and healthy controls using ultrasound imaging in weight-bearing. METHODS: A case-control study in a university lab setting was conducted. Thirty individuals with PFP (age: 20.23 ± 3.30 years, mass: 74.70 ± 27.63 kgs, height: 161.32 ± 11.72 cm) and 30 healthy individuals (age: 20.33 ± 3.37 years, mass: 64.02 ± 11.00 kgs, height: 169.31 ± 9.30 cm) participated. Cross-sectional area (CSA) images of the peroneal muscles were taken in non-weight bearing and weight-bearing positions. The functional activation ratio from lying to single-leg standing (SLS) was calculated. RESULTS: There was a statistically significant (p = 0.041) group (PFP, healthy) by position (non-weight-bearing, weight-bearing) interaction for the peroneal muscle CSA with a Cohen's d effect size of 0.2 in non-weight-bearing position and 0.7 in weight-bearing position. The functional activation ratio for the healthy group was significantly more (p = 0.01) than the PFP group. CONCLUSION: Peroneal muscles were found to be smaller in size in those with PFP compared to the healthy subjects in the weight-bearing SLS position. This study found that those with PFP have lower activation of peroneal muscles in functional position.

Differences in ankle and knee muscle architecture and plantar pressure distribution among women with knee osteoarthritis.

BACKGROUND: The aim of this study was to compare the plantar pressure distribution and knee and ankle muscle architecture in women with and without knee osteoarthritis (OA). METHODS: Fifty women with knee OA (mean age = 52.11 ± 4.96 years, mean Body mass index (BMI) = 30.94 ± 4.23 kg/m2) and 50 healthy women as a control group (mean age = 50.93 ± 3.78 years, mean BMI = 29.06 ± 4.82 kg/m2) were included in the study. Ultrasonography was used to evaluate knee and ankle muscles architecture and femoral cartilage thickness. The plantar pressure distribution was evaluated using the Digital Biometry Scanning System and Milleri software (DIASU, Italy). Static foot posture was evaluated using the Foot Posture Index (FPI), and pain severity was assessed using the Visual Analog Scale. RESULTS: The OA group exhibited lower muscle thickness in Rectus Femoris (RF) (p = 0.003), Vastus Medialis (VM) (p = 0.004), Vastus Lateralis (p = 0.023), and Peroneus Longus (p = 0.002), as well as lower Medial Gastrocnemius pennation angle (p = 0.049) and higher Fat thickness (FT) in RF (p = 0.033) and VM (p = 0.037) compared to the control group. The OA group showed thinner femoral cartilage thickness (p = 0.001) and higher pain severity (p = 0.001) than the control groups. FPI scores were higher (p = 0.001) in OA group compared to the control group. The plantar pressure distribution results indicated an increase in total surface (p = 0.027), total load (p = 0.002), medial load (p = 0.005), and lateral load (p = 0.002) on dominant side in OA group compared to the control group. CONCLUSIONS: Knee and ankle muscle architecture, knee extensor muscle FT, and plantar pressure distribution in the dominant foot differed in individuals with knee OA compared to the control group.

Ankle and hindfoot motion of healthy adults during running revealed by dynamic biplane radiography: Side-to-side symmetry, sex-specific differences, and comparison with walking.

This study aimed to characterize ankle and hindfoot kinematics of healthy men and women during overground running using biplane radiography, and to compare these data to those previously obtained in the same cohort during overground walking. Participants ran across an elevated platform at a self-selected pace while synchronized biplane radiographs of their ankle and hindfoot were acquired. Motion of the tibia, talus, and calcaneus was tracked using a validated volumetric model-based tracking process. Tibiotalar and subtalar 6DOF kinematics were obtained. Absolute side-to-side differences in ROM and kinematics waveforms were calculated. Side-to-side and sex-specific differences were evaluated at 10 % increments of stance phase with mixed model analysis. Pearson correlation coefficients were used to assess the relationship between stance-phase running and walking kinematics. 20 participants comprised the study cohort (10 men, mean age 30.8 ± 6.3 years, mean BMI 24.1 ± 3.1). Average absolute side-to-side differences in running kinematics waveforms were 5.6°/2.0 mm or less at the tibiotalar joint and 5.2°/3.2 mm or less at the subtalar joint. No differences in running kinematics waveforms between sides or between men and women were detected. Correlations were stronger at the tibiotalar joint (42/66 [64 %] of correlations were p < 0.05), than at the tibiotalar joint (38/66 [58 %] of correlations were p < 0.05). These results provide a normative reference for evaluating native ankle and hindfoot kinematics which may be informative in surgical or rehabilitation contexts. Sex-specific differences in ankle kinematics during overground running are likely not clinically or etiologically significant. Associations seen between walking and running kinematics suggest one could be used to predict the other.

A method for automated masking and plantar pressure analysis using segmented computed tomography scans.

BACKGROUND: Plantar pressure, a common gait and foot biomechanics measurement, is typically analyzed using proprietary commercial software packages. Regional plantar pressure analysis is often reported in terms of underlying bony geometry, and recent advances in image processing and accessibility have made computed tomography, radiographs, magnetic resonance imaging, or other imaging methods more popular for incorporating bone analyses in biomechanics. RESEARCH QUESTION: Can a computed tomography-based regional mask provide comparable regional analysis to commercial plantar pressure software and can the increased flexibility of an in-house method obtain additional insight from common measurements? METHODS: A plantar pressure analysis method was developed based on bony geometry from computed tomography scans to calculate peak pressure, pressure time integral incorporating sub-peak values, force time integral, pressure gradient, and pressure gradient angle. Static and dynamic plantar pressure were acquired for 4 subjects (male, 65 ± 2.4 years). Plantar pressure variables were calculated using commercial and computed tomography-based systems. RESULTS: Dynamic peak pressure, pressure time integral, and force-time integral computed using the bone-based software was 5 % (9kPa), 7 % (0.3kPa-s) and 13 % (0.3 N-s) different than the commercial software on average. Region masks of the metatarsals and toes differed between commercial and computed tomography-based software due to subject-specific bone geometry and toe shape. Pressure time integral values incorporating sub-peak pressure were higher and demonstrated higher relative hindfoot values compared to those without. Removing step-on frames to static pressure analysis decreased forefoot pressures. Regional maps of peak pressure and maximum pressure gradient demonstrate different peak locations. SIGNIFICANCE: Computed tomography-based regional masks are comparable to commercial masks. Inclusion of static step-on frames and sub-peak pressures may change regional plantar pressure patterns. Differences in location of maximum pressure gradient and peak pressure may be useful for assessing subject specific injury risk.

Skin marker-based versus bone morphology-based coordinate systems of the hindfoot and forefoot.

Segment coordinate systems (CSs) of marker-based multi-segment foot models are used to measure foot kinematics, however their relationship to the underlying bony anatomy is barely studied. The aim of this study was to compare marker-based CSs (MCSs) with bone morphology-based CSs (BCSs) for the hindfoot and forefoot. Markers were placed on the right foot of fifteen healthy adults according to the Oxford, Rizzoli and Amsterdam Foot Model (OFM, RFM and AFM, respectively). A CT scan was made while the foot was loaded in a simulated weight-bearing device. BCSs were based on axes of inertia. The orientation difference between BCSs and MCSs was quantified in helical and 3D Euler angles. To determine whether the marker models were able to capture inter-subject variability in bone poses, linear regressions were performed. Compared to the hindfoot BCS, all MCSs were more toward plantar flexion and internal rotation, and RFM was also oriented toward more inversion. Compared to the forefoot BCS, OFM and RFM were oriented more toward dorsal and plantar flexion, respectively, and internal rotation, while AFM was not statistically different in the sagittal and transverse plane. In the frontal plane, OFM was more toward eversion and RFM and AFM more toward inversion compared to BCS. Inter-subject bone pose variability was captured with RFM and AFM in most planes of the hindfoot and forefoot, while this variability was not captured by OFM. When interpreting multi-segment foot model data it is important to realize that MCSs and BCSs do not always align.

Non-contrast MRI of micro-vascularity of the feet and toes.

PURPOSE: This study aimed to develop novel non-contrast MR perfusion techniques for assessing micro-vascularity of the foot in human subjects. METHODS: All experiments were performed on a clinical 3 T scanner using arterial spin labeling (ASL). Seven healthy subjects (30-72 years old, 5 males and 2 females) were enrolled and bilateral feet were imaged with tag-on and tag-off alternating inversion recovery spin labeling for determining micro-vascularity. We compared an ASL technique with 1-tag against 4-tag pulses. For perfusion, we determined signal increase ratio (SIR) at varying inversion times (TI) from 0.5 to 2 s. SIR versus TI data were fit to determine perfusion metrics of peak height (PH), time to peak (TTP), full width at half maximum (FWHM), area under the curve (AUC), and apparent blood flow (aBF) in the distal foot and individual toes. Using analysis of variance (ANOVA), effects of tag pulse and region of interest (ROI) on the mean perfusion metrics were assessed. In addition, a 4-tag pulse perfusion experiment was performed on patients with peripheral artery disease (PAD) and Raynaud's disease. RESULTS: Using our MR perfusion techniques, SIR versus TI data showed well-defined leading and trailing edges, with a peak near TI of 0.75-1.0 s and subsiding quickly to near zero by TI of 2 s, particularly when 4-tag pulses were used. When imaged with 4-tag pulse, we found significantly greater values in perfusion metrics, as compared to 1-tag pulse. The patients with PAD and Raynaud's disease showed a reduced or scattered perfusion curves compared to the healthy control. CONCLUSION: MR perfusion imaging of the distal foot shows greater SIR and perfusion metrics with the 4-tag pulse compared to the 1-tag pulse technique. This will likely benefit those with low perfusion due to aging, PAD, diabetic foot, and other vascular diseases.