WO2007067987A2 - Systems and methods for elastography imaging - Google Patents
Systems and methods for elastography imaging Download PDFInfo
- Publication number
- WO2007067987A2 WO2007067987A2 PCT/US2006/061809 US2006061809W WO2007067987A2 WO 2007067987 A2 WO2007067987 A2 WO 2007067987A2 US 2006061809 W US2006061809 W US 2006061809W WO 2007067987 A2 WO2007067987 A2 WO 2007067987A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- strain
- displacement
- ultrasound
- information
- interest
- Prior art date
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4514—Cartilage
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4523—Tendons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4533—Ligaments
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56358—Elastography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/40—Animals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F2/00—Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
- A61F2/02—Prostheses implantable into the body
- A61F2/08—Muscles; Tendons; Ligaments
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F2/00—Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
- A61F2/02—Prostheses implantable into the body
- A61F2/30—Joints
- A61F2/30756—Cartilage endoprostheses
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F2/00—Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
- A61F2/02—Prostheses implantable into the body
- A61F2/30—Joints
- A61F2/3094—Designing or manufacturing processes
- A61F2/30942—Designing or manufacturing processes for designing or making customized prostheses, e.g. using templates, CT or NMR scans, finite-element analysis or CAD-CAM techniques
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F2/00—Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
- A61F2/02—Prostheses implantable into the body
- A61F2/30—Joints
- A61F2/46—Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor
- A61F2/468—Testing instruments for artificial joints
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Definitions
- the present invention relates to the use of cystography in the mechanical characterization of anatomical structures associated with mammalian joints and tendons.
- ACL human anterior cruciate ligament
- BPTB bone-patellar tendon-bone
- HT hamstring tendon
- HT grafts which are the most commonly used grafts for ACL reconstruction. While harvesting of HT grafts results in significantly less donor site morbidity, these grafts arc mechanically anchored, and their clinical success is limited by the lack of biological graft integration with the subchondral bone.
- Native ACL inserts into bone through a direct insertion consisting of a linear transition from ligament to fibrocartilage to bone.
- the f ⁇ brocartilage zone is further divided into non-mineralized and mineralized fibrocartilage regions.
- the ACL-bone interface Due to the presence of several types of tissue, the ACL-bone interface is expected to vary in cellular, chemical, and mechanical properties. It is believed this controlled heterogeneity permits the transition of mechanical load between bone and soft tissue and minimizes the formation of stress concentrations. This interface, however, is not re-established after tendon graft-based ACL reconstruction. Without a stable interface, the fixation site of the HT grafts to bone becomes the weak link in the reconstructed graft, a leading cause of graft failure and resulting in revision surgery.
- Osteoarthritis is another condition that affects many millions of people worldwide. Osteoarthritis is a disease process involving articular cartilage. Articular cartilage is poroelastic and bears load in articular joints. The use of radiography and physical examination to examine nascent osteoarthritis is quite limited, however.
- the present invention includes the use of ultrasound elastography to determine strain distribution of joint structures and tendons.
- Joints may include but are not limited to those of the foot, ankle, hip, temporomandibular joint (TMJ), shoulder, elbow, hand and wrist and corresponding anatomical structures in non-human mammals.
- Joint structures may include, for example, ligaments such as the ACL, cartilage (especially articular cartilage), and the medial and lateral menisci of, for example, the tibiofemoral joint.
- Tendons may include, for example, the achilles tendon and flexor and extensor tendons of mammalian extremities.
- Exemplary embodiments of the present invention provide methods for obtaining information about the mechanical behaviour of structures associated with. mammalian joints where such methods include, creating deformation in a joint structure of interest, using an ultrasound scanner and a linear array to acquire sequences of ultrasound data of the joint structure, and estimating the axial displacement between a reference frame of the data and successive frames of the data.
- Deformation may be either active or passive and can include, for example, tension, compression, relaxation and combination thereof.
- Another exemplary embodiment of the present invention provides methods for obtaining information about the mechanical behaviour of tendons wherein such methods include, creating deformation in a tendon of interest, using an ultrasound scanner and a linear array to acquire sequences of ultrasound data of the tendon, and estimating the axial displacement between a reference frame of the data and successive frames of the data.
- deformation may be either active or passive and can include, for example, tension,
- the present invention provides methods which include estimating axial displacement and strain using a ID cross-correlation algorithm.
- the present invention provides methods which include estimating 2D and/or 3D axial displacement and strain.
- the present invention includes the use of cross-correlation algorithms to determine time-shifts between two backscattered signals by cross-correlating sliding windows over a 2D ultrasound image.
- a further exemplary embodiment of the present invention includes using information concerning the mechanical characterization of structures associated with joints and tendons to inform the design of tissue grafts.
- Exemplary embodiments of the present invention also include tissue grafts produced using information obtained by the mechanical characterization of structures associated with, joints and tendons.
- Exemplary embodiments of the present invention also allow for imaging of displacement and strain as well as estimation of displacement and strain.
- Figure 1 illustrates a flow chart for a representative recorrelation technique.
- Figures 2A and 2B provide a representative image produced by tracking an RF segment in 2D for axial ( ⁇ s) and lateral ( ⁇ l) displacement estimation.
- A's and B's are RF lines corresponding to consecutive frames in time.
- FIG. 3A shows an exemplary equipment configuration for ultrasound data acquisition during tensile testing.
- Neonatal bovine patellofemoral joint is loaded on a mechanical testing system (MTS) modified with a cylindrical polycarbonate tank.
- the tank is filled with physiologic saline and the ultrasound transducer (arrow) is mounted inside the tank such that the ACL and insertions can be scanned posteriorly.
- MTS mechanical testing system
- Figure 3B shows the neonatal bovine patellofemoral joint shown in Figure 3 A mounted in the MTS with a tibial orientation and 0° flexion.
- Figure 4A shows a posterior view of scanned ACL and insertions.
- Figure 4B shows an ultrasound image of tibial insertion.
- Figure 4C shows a corresponding y-displacement map (mm), with blue to red indicating small to large displacements, respectively.
- Figure 4D shows a corresponding elastogram with compressive strain (not in
- Figures 5A and B show, respectively, displacement (5A) and strain (5B) at the tibial insertion as collected from point-wise temporal analysis. Point of data analysis is indicated by the arrow shown in Figure 4.
- Figure 7A shows an ultrasound image of an ACL and insertions with the transducer rotated so that the face of the transducer was aligned along the principal axis of the
- Figures 7B and 7C show, respectively, axial displacement (mm) (7B) and strain (not in %) (7C) during tensile stretching in the lateral direction versus the axial beam propagation.
- a compressive (blue) strain can be seen in 7B corresponding to the ACL-bone interface (orange region in Figure 7A).
- Figure 8 shows a representative design of a compression apparatus and image acquisition arrangement for elastographic imaging of cartilage.
- Figure 9 graphically illustrates a representative load versus time for an imposed pre-strain of 10% followed by an additional 2% strain 90 seconds later
- the load axis is on the right of the graph and the strain axis is on the right.
- Figures IQA and IQB show, respectively, a representative grey-scale RF signal of a cartilage sample within a compression apparatus with a 1/8" inch opening (10A) and a
- Figure 1 1 A provides a representative displacement image of a femoral condyle cartilage sample showing uniform displacement in the region of interest, using a 1/8" opening with a scale ranging from -0. lmm to 0.1 mm.
- Figure 1 IA where the sample shows essentially zero local strain except at the surface and interface of zone 1 and 2 with a scale ranging from -0.9% to 0.9%.
- Figure 12A shows a representative displacement image of a femoral condyle cartilage sample showing a slight displacement gradient in the region of interest, using a
- the loading plate is indicated by an arrow.
- the scale ranges from -0. lmm to
- Figure 12B shows a representative elastogram of the same sample shown in
- Figure 12A where the sample shows essentially zero local strain except at the surface and interface of zone 1 and 2 with a scale ranging from -0.2% to 0.2%.
- ACL grafts used for ACL reconstruction is limited by healing of the graft with bone, which results in non-anatomical fibrovascular scar tissue at the interface between the tendon graft and. bone.
- the success of ACL grafts depends on reforming the native anatomical tendon-bone interface.
- Methods according to the present invention allow for the determination of the structure- function relationship of joint structures, such as the ACL. This information then may be used to design scaffold systems, such as tendon grafts to bone, which mimic the native tissue in morphology, chemical composition, cellular distribution, and mechanical properties.
- Methods according to the present invention are particularly useful for characterizing strain applied to structures associated with joints such as ligaments and cartilage, as well as tendons.
- Much of the description below and in the Examples section is directed to methods involving characterization of the anterior cruciate ligament and the characterization of articular cartilage of the femoral condyles. It should be understood, however, that these methods are readily applicable to other anatomical structures, for example, other ligaments, cartilage and tendons, for which detailed characterization of strain responses is desired.
- Minimally invasive conditions include, for example, endoscopic procedures.
- Methods according to the present invention also can be used for veterinary applications in, for example, equine, canine and feline species.
- Methods according to the present invention are ideal for examining, for example, the ACL-bone interface as such methods permit the characterization of relatively small areas (on the order of about 0.1-2 mm, depending on the ultrasound frequency used) with complex stress distributions.
- an ultrasound transducer scans a region of interest while an external load is applied to induce strain. Speckle tracking techniques may be employed to analyze the collected radio-frequency ultrasonic data before and after incremental loading and to estimate the resulting strain and strain distributions. Standard ultrasound scanners (e.g. Terason 2000, Teratech, Framingham MA) or similar devices may be used.
- Ultrasound frequencies may range from, for example, 2-40 MHz, although higher or lower frequencies may be appropriate in certain circumstances, as will be appreciated by those of ordinary skill in the art. Sequences of RF data may be acquired during loading of the joint. Axial displacement between a reference and successive frames can be estimated using cross-correlation and recorrelation techniques, exemplary
- Axial displacements, or displacements occurring in the direction orthogonal to the face of the transducer and parallel to the direction of ultrasound propagation can be estimated for each RF frame with respect to a reference frame by using a ID cross-correlation algorithm. Strain distribution can be computed by differentiating the displacement map along the axial direction. For numerical differentiation, a least-squares regression method may be used. Displacement and strain then may be estimated relative to a reference frame to obtain a temporal profile and map of the cumulative deformation at the ligament and the insertion. These methods may be used to generate maps of cumulative deformation and strain in, for example, the ACL and tibial insertion during tensile loading. These techniques can be used to generate detailed information about displacements and strains associated with the tibiofemoral joint and other joints and tendons.
- Cross-correlation and recorrelation algorithms may be used to obtain detailed information as described below.
- Figure 1 provides a representative sequence of cross- correlation and recorrelation algorithms. These techniques may be applied between consecutive echo (RF, envelope-detected or B-mode) axial (ID) segments in 2D (in-plane) or 3D (in-plane and out-of-plane).
- each frame may be considered as a set of echo segments.
- each echo segment in a first frame may be cross-correlated with echo segments in a second frame to find the best match.
- the 2D or 3D resultant that denotes the path of motion between segments 1 and 2 may be broken into its individual components in the axial (in-plane, along the propagation axis; ⁇ s in Figure 2), lateral (in-plane, orthogonal to the propagation axis; ⁇ l in Figure 2) and elevational (out-of-planc).
- rccorrclation i.e., after correction, or removal, for the motion occurring in a direction perpendicular to that of estimation.
- the lateral component may be removed first from original Frame 2 (i.e., frame 2_1 may be generated) for the estimation of the new axial component (axial estimate 2; Figure 1) and the second lateral component may be estimated (lateral estimate 2; Figure 1) after removal of the new axial component (i.e., frame 2_2 may be generated) and so on until reaching high enough correlation for the iterations to no longer be useful ( Figure 1).
- cross-correlation technique allows for effective decoupling of the 2D or 3D components and higher quality displacement and strain estimates and images of the tissue under deformation. It should be noted that the cross-correlation technique also may be employed without the use of recorrelation methods in the case where estimation of one or all components is considered to be of sufficient quality (e.g., high signal-to-noise ratio or high contrast-to-noise ratio).
- Information generated by methods according to the present invention may be used to inform the design and material selection for the production of tissue grafts.
- biodegradable scaffolds can be found in the literature, including Lu ct al., Biomatcrials 26 (2005) 4805-4816, the contents of which arc incorporated herein by reference in its entirety.
- Methods according to the present invention allow for graft designs and material selection to take into consideration the detailed strain response for a particular anatomical structure, such as an ACL. Both graft design and material selection are very important for long term clinical success and involve a balance between scaffold
- Exemplary materials include those comprising poly-alpha-hydroxyesters such as polyglycoKc acid (PGA), poly-L-lactic acid (PLLA), and polylactic-co-glycolic copolymer (PLAGA), all of which have been approved by the FDA for these purposes. These types of degradable polymers do not elicit a permanent foreign body reaction and are gradually reabsorbed and replaced by natural tissue.
- poly-alpha-hydroxyesters such as polyglycoKc acid (PGA), poly-L-lactic acid (PLLA), and polylactic-co-glycolic copolymer (PLAGA)
- Neonatal bovine calf up to one week old tibiofemoral joints obtained from an abattoir (Fresh Farm Beef, Vermont) were used for this Example. After removal of surrounding muscle and adipose tissue, the joint capsule was opened. Fascia lata and connective tissue were removed from the joint capsule with the ACL and posterior cruciate ligament (PCL) undisturbed. The PCL was maintained intact until immediately prior to testing in order to maintain joint stability and prevent premature damage to the ACL. During all joint preparation procedures, the ACL and surrounding tissues were kept hydrated with physiologic saline.
- the femur and tibia were cut to approximately 12 cm from the joint with a hacksaw, the periosteum removed, and bone marrow extracted from the intramedullary cavity to improve cement fixation of the joint. Subsequently, the tibia and femur were secured with custom anchors and cement to prevent slippage during testing.
- the joint was then mounted on a uniaxial material testing system (MTS 858 Bionix Testing System; MTS, Eden Prairie, MN) fitted with a custom cylindrical polycarbonate tank, the PCL was severed, and the medial femoral condyle was removed with a hacksaw to improve line-of-sight access to the ACL and insertions for the ultrasound transducer ( Figures 3 A and 3B). These procedures may be applied in vivo for veterinary medicine purposes or in humans for assessing injury or age-related diseases such as osteoporosis. Materials and Methods Tensile Testing
- FTC Fibre Channel tibial alignment
- the femur and tibia were aligned along the tensile axis with 0° of knee flexion, and the sample was submerged in degassed physiologic saline.
- the saline provided a medium for ultrasound propagation.
- a preload of 2 N was applied for one minute, and the joint was preconditioned by cyclic sawtooth loading from 0— 0.75 mm for 10 cycles at 20 rnm/min followed by a rest of 1 min.
- the joint was cyclically loaded from 0 - 2 mm at 20 mm/min, with 0 mm being the displacement during the preload. Following a 30 minute rest, the joint was cyclically reloaded from 0 - 3 mm, with additional displacement applied during this testing regimen to ensure a detectable amount of deformation occurred across the insertions. Finally, after an additional 30 minute rest, the joint was loaded to failure at 10 mm/min.
- RF Radio Frequency
- the ultrasound transducer was mounted inside the saline tank and positioned to image the ACL and insertions. Sequences of ultrasound RF data were acquired continuously during the applied loading repeatedly for periods of 3 seconds at 54 frames/s (128 RF lines, sampling frequency: 10 MHz). The axial displacement between a reference and successive frames was estimated offline and imaged using cross-correlation and recorrelation techniques with a window size of 3 mm and a window overlap of 80%.
- time-shifts between two backscattered signals are determined by the cross-correlation of small sliding windows over the entire 2D ultrasound image.
- recorrelation techniques were employed.
- strain distribution was computed by differentiating the displacement map along the axial direction. For the numerical differentiation, a least-squares regression method was used. Displacement and strain were estimated relative to a reference frame, which was captured at the beginning of the application of tensile load, in order to obtain a temporal profile and map of the cumulative deformation at the ligament and the insertion.
- Ligament and bone are less dense and therefore less echogenic, enabling the structure of these tissues to be discernable on the B-mode images.
- ACL a narrow band of high strain in the middle and along the length of ACL was noted that also corresponded to a highly echogenic area on the B-scan images. This may reflect the parallel bundle- organization of the ACL. Distinctions between soft and hard tissue signatures were used to identify the ACL insertions into bone.
- Figure 4C shows the distribution of deformation throughout the FATC, with magnitudes of deformation represented according to a colormap, with small deformations blue and large deformations red.
- the magnitude of displacement was found to be the highest (red in Figure 4C) within the ACL proper and decreased in value in a gradual transition (orange, yellow, and green) from ligament (red) to bone (blue).
- elastographic analysis revealed through strain maps that the strain profile at the tibial insertion was highly complex as the FATC was loaded in tension (Figure 4D). Both compressive and tensile strains were visualized at the tibial insertion site, indicated by the green-blue and yellow-red regions, respectively, on the elastogram in Figure 4D.
- fibrocartilaginous transitional tissue between ligament and bone demonstrates that a compressive strain component exists in that region during physiological loading.
- collagen fibers extending from ligament into bone at the insertions when loaded in tension, transmit shear and compressive stresses through the fibrocartilage zones of the insertions.
- Example 1 Techniques applied in Example 1 to anterior cruciate ligaments in vitro, arc applied in situ to characterize ACL ligaments. Highly detailed data characterizing strain responses of ACL ligaments arc obtained.
- Example 1 Techniques applied in Example 1 to anterior cruciate ligaments are applied to other structures of the tibiofemoral joint, in vitro and in situ, including the posterior cruciate ligament, cartilage, and medial and lateral menisci. Highly detailed data characterizing strain responses of these structures are obtained.
- Example 1 Techniques applied in Example 1 are applied to other joints of the upper and lower extremities including the foot, ankle, hip, temporomandibular joint (TMJ), shoulder, elbow, hand and wrist, in vitro and in situ. Highly detailed data characterizing strain responses of structures associated with these joints are obtained.
- TMJ temporomandibular joint
- Example 1 Techniques applied in Example 1 are applied to tendons, such as the achilles tendon and flexor and extensor tendons of mammalian extremities, in vitro and in situ.
- the specimens were oriented such that the deep portion of the cartilage contacted an aluminum loading plate and the articular surface rested upon another rigid, impermeable surface containing a 3 mm opening to serve as the acoustic window for the high-resolution ultrasound transducer (f/2, 8mm focus, 55 MHz, 46 Hz frame rate, Vevo 770, Visualsonics, Toronto, Canada).
- the ultrasound probe was separated by 3-mm from the surface of the articular cartilage.
- a tare strain of 0.1% based on the measurement of the undeformed cartilage plugs was sustained for 30 seconds, followed by a ramp to strains ranging from 0.5 to 4.0% strain at 0.1 mm/sec for two femoral condyle samples and one femoral head sample, the results of which are graphically described in Figure 9.
- B-mode ultrasound scans were acquired immediately after the tare strain and immediately after ramped compression of the cartilage.
- a pre-strain of 10% based on the measurement of the undeformed cartilage plugs was applied for 30 seconds, followed by a ramp to strains ranging from 2% at 0.1 mm/sec for two femoral head samples.
- Radio frequency (RF) ultrasound signals of these samples were acquired once equilibrium was attained.
- Displacement images and elastograms were generated using ID crosscorrelation techniques and gradient operators on the RF signals (window size 0.3 mm, 85% overlap), respectively. Median filtering of the displacement data was also implemented. This arrangement simulates a device in which an ultrasound transducer is incorporated into an arthroscopic indentation device.
- This Example confirms the usefulness of high resolution ultrasound elastographic imaging of articular cartilage, for example, for the early diagnosis and monitoring of treatment for articular cartilage pathologies, such as osteoarthritis.
- Example 6 Techniques applied in Example 6 to articular cartilage in vitro, arc applied in situ. Highly detailed data images characterizing local strains to articular cartilage are obtained.
- Example 6 Techniques applied in Example 6 are applied to articular cartilage of the joints of the upper and lower extremities including the ankle, hip, shoulder, elbow and wrist, in vitro and in situ. Highly detailed data images characterizing local strains of articular cartilage associated with these joints are obtained.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Rheumatology (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Rehabilitation Therapy (AREA)
- Quality & Reliability (AREA)
- Vascular Medicine (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Methods for obtaining information about the mechanical behaviour of structures associated with mammalian joints and tendons are provided. Embodiments of such methods include creating deformation in a joint structure (such as ligaments and articular cartilage) or tendon of interest, using an ultrasound scanner and a single element or array of elements to acquire sequences of ultrasound data of the joint structure or tendon, estimating one, two or three components of the resulting displacement and strain between a reference frame of ultrasound data and successive frames of ultrasound data, and using a cross-correlation algorithm to estimate the displacement and strain components. This information may be used to inform the design of tissue grafts. Tissue grafts produced using this information are also provided. The same method can be used in situ together with noninvasive or invasive procedures.
Description
SYSTEMS AND METHODS POR ELASTOGRAPHY IMAGING
RELATED APPTJCATIONS
[00001] This application claims priority to and benefit of U.S. provisional application no. 60/749,245 filed December 9, 2005 and U.S. provisional application no. 60/751,122 filed December 16, 2005, the contents of each of which is incorporated herein in its entirety.
FIELD OF THE INVENTION
[00002] The present invention relates to the use of cystography in the mechanical characterization of anatomical structures associated with mammalian joints and tendons.
BACKGROUND
[00003] Elasticity imaging, or methods for the mapping of mechanical responses or properties using ultrasound or MRI images acquired before and after a mechanical excitation, was initially developed in the early 1990s as an alternative tool for early tumor diagnosis based on the principle of palpation. It has since been utilized in intravascular and
cardiovascular applications in vivo, as well as in the guidance of thermal therapy procedures and robotic surgery. Elastography has also been used to determine mechanical properties of tissues by measuring both elastic and viscoelastic properties.
[00004] Despite these recent advances, a need still exists for obtaining detailed mechanical characterization of joint structures and tendons in both human and veterinary medicine. For example, detailed characterization of the manner in which a human anterior cruciate ligament (hereinafter "ACL") responds to strain would be beneficial in designing appropriate replacement grafts. In the United States alone approximately 75,000-100,000 ACL reconstructions are performed each year.
[00005] Due to the poor healing potential of the ACL, surgical intervention is often required following ligament injury. Surgery most often constitutes complete reconstruction of the ACL using an autologous tendon graft. Two autografts commonly utilized to replace
the ACL are the bone-patellar tendon-bone (BPTB) graft and the semitendinosus or hamstring tendon (HT) graft. While the BPTB graft, which consists of the central third of the patellar tendon, has the advantage of possessing bony ends that facilitate integration with bone in the femoral and tibial bone tunnels, harvesting of these autografts often results in significant donor site morbidity. Consequently, there has been a shift toward the usage of HT grafts, which are the most commonly used grafts for ACL reconstruction. While harvesting of HT grafts results in significantly less donor site morbidity, these grafts arc mechanically anchored, and their clinical success is limited by the lack of biological graft integration with the subchondral bone.
[00006] Native ACL inserts into bone through a direct insertion consisting of a linear transition from ligament to fibrocartilage to bone. The fϊbrocartilage zone is further divided into non-mineralized and mineralized fibrocartilage regions. Due to the presence of several types of tissue, the ACL-bone interface is expected to vary in cellular, chemical, and mechanical properties. It is believed this controlled heterogeneity permits the transition of mechanical load between bone and soft tissue and minimizes the formation of stress concentrations. This interface, however, is not re-established after tendon graft-based ACL reconstruction. Without a stable interface, the fixation site of the HT grafts to bone becomes the weak link in the reconstructed graft, a leading cause of graft failure and resulting in revision surgery.
[00007] Promoting clinically successful graft-bone integration via the regeneration of the native ligament-bone interface of the ACL (or other joint structures) requires a detailed understanding of the mechanical properties of both the ligament and the ligament-bone interface. This would allow for the proper consideration of material selection and interface regeneration via tissue engineering methods. Elucidating the structure-function relationship at the insertions would be critical in the design of the next generation of graft fixation
devices, enabling biological fixation of tendon grafts to bone through the reestablishment of the native tendon to bone interface.
100008] Osteoarthritis is another condition that affects many millions of people worldwide. Osteoarthritis is a disease process involving articular cartilage. Articular cartilage is poroelastic and bears load in articular joints. The use of radiography and physical examination to examine nascent osteoarthritis is quite limited, however.
[00009] Consequently, there have been attempts to develop arthroscopic indentation devices that enable measurement of the biomcchanical properties of degenerate cartilage, with the aim of early diagnosis. In particular, some efforts have been made to combine arthroscopic indenters with ultrasound probes, mainly to obtain more accurate estimates of cartilage thickness. But in order for arthroscopy to become the gold standard for diagnosing cartilage pathology, data acquisition and presentation must be improved upon greatly.
[00010] Detailed mechanical characterization, data acquisition and presentation of structures associated with joints (such as ligaments, cartilage and menisci) and tendons would be very beneficial.
SUMMARY
[00011] The present invention includes the use of ultrasound elastography to determine strain distribution of joint structures and tendons. Joints may include but are not limited to those of the foot, ankle, hip, temporomandibular joint (TMJ), shoulder, elbow, hand and wrist and corresponding anatomical structures in non-human mammals. Joint structures may include, for example, ligaments such as the ACL, cartilage (especially articular cartilage), and the medial and lateral menisci of, for example, the tibiofemoral joint. Tendons may include, for example, the achilles tendon and flexor and extensor tendons of mammalian extremities.
[00012] Exemplary embodiments of the present invention provide methods for obtaining information about the mechanical behaviour of structures associated with.
mammalian joints where such methods include, creating deformation in a joint structure of interest, using an ultrasound scanner and a linear array to acquire sequences of ultrasound data of the joint structure, and estimating the axial displacement between a reference frame of the data and successive frames of the data. Deformation may be either active or passive and can include, for example, tension, compression, relaxation and combination thereof.
[00013] Another exemplary embodiment of the present invention provides methods for obtaining information about the mechanical behaviour of tendons wherein such methods include, creating deformation in a tendon of interest, using an ultrasound scanner and a linear array to acquire sequences of ultrasound data of the tendon, and estimating the axial displacement between a reference frame of the data and successive frames of the data. Again, deformation may be either active or passive and can include, for example, tension,
compression, relaxation and combination thereof.
[00014] In yet another exemplary embodiment, the present invention provides methods which include estimating axial displacement and strain using a ID cross-correlation algorithm.
[00015] In still another exemplary embodiment, the present invention provides methods which include estimating 2D and/or 3D axial displacement and strain.
[00016] In another exemplary embodiment, the present invention includes the use of cross-correlation algorithms to determine time-shifts between two backscattered signals by cross-correlating sliding windows over a 2D ultrasound image.
[00017] A further exemplary embodiment of the present invention includes using information concerning the mechanical characterization of structures associated with joints and tendons to inform the design of tissue grafts.
[00018] Exemplary embodiments of the present invention also include tissue grafts produced using information obtained by the mechanical characterization of structures associated with, joints and tendons.
[00019J Exemplary embodiments of the present invention also allow for imaging of displacement and strain as well as estimation of displacement and strain.
BRIEF DESCRIPTION OF THE DRAWINGS
[00020] Figure 1 illustrates a flow chart for a representative recorrelation technique.
[00021] Figures 2A and 2B provide a representative image produced by tracking an RF segment in 2D for axial (Δs) and lateral (Δl) displacement estimation. A's and B's are RF lines corresponding to consecutive frames in time.
[00022] Figure 3A shows an exemplary equipment configuration for ultrasound data acquisition during tensile testing. Neonatal bovine patellofemoral joint is loaded on a mechanical testing system (MTS) modified with a cylindrical polycarbonate tank. The tank is filled with physiologic saline and the ultrasound transducer (arrow) is mounted inside the tank such that the ACL and insertions can be scanned posteriorly.
[00023] Figure 3B shows the neonatal bovine patellofemoral joint shown in Figure 3 A mounted in the MTS with a tibial orientation and 0° flexion.
[00024] Figure 4A shows a posterior view of scanned ACL and insertions.
[00025] Figure 4B shows an ultrasound image of tibial insertion.
[00026] Figure 4C shows a corresponding y-displacement map (mm), with blue to red indicating small to large displacements, respectively.
[00027] Figure 4D shows a corresponding elastogram with compressive strain (not in
%) indicated by blue regions and tensile strain by red regions. Negative values denote compressive strain.
[00028] Figures 5A and B show, respectively, displacement (5A) and strain (5B) at the tibial insertion as collected from point-wise temporal analysis. Point of data analysis is indicated by the arrow shown in Figure 4.
[00029 J Figure 6 shows sequences (1-16) of a) displacement images and b)
elastograms at 15 frames apart during tensile loading. Blue denotes high compressive strain and red denotes high tensile strain on the elastograms. Development of high displacements in the ligament compared to bone during tensile loading arc shown as indicated by the change from blue to red in the ligament proper.
[00030] Figure 7A shows an ultrasound image of an ACL and insertions with the transducer rotated so that the face of the transducer was aligned along the principal axis of the
ACL.
[00031] Figures 7B and 7C show, respectively, axial displacement (mm) (7B) and strain (not in %) (7C) during tensile stretching in the lateral direction versus the axial beam propagation. A compressive (blue) strain can be seen in 7B corresponding to the ACL-bone interface (orange region in Figure 7A).
[00032] Figure 8 shows a representative design of a compression apparatus and image acquisition arrangement for elastographic imaging of cartilage.
[00033] Figure 9 graphically illustrates a representative load versus time for an imposed pre-strain of 10% followed by an additional 2% strain 90 seconds later
demonstrating stress-relaxation of articular cartilage. The load axis is on the right of the graph and the strain axis is on the right.
[00034] Figures IQA and IQB show, respectively, a representative grey-scale RF signal of a cartilage sample within a compression apparatus with a 1/8" inch opening (10A) and a
1/32" opening (10B). In Figures 1OA and 1OB the bracket represents the opening in the loading plate, while he arrow show the highly reflective metal loading plate.
[00035] Figure 1 1 A provides a representative displacement image of a femoral condyle cartilage sample showing uniform displacement in the region of interest, using a 1/8" opening with a scale ranging from -0. lmm to 0.1 mm.
[00036J Figure HB provides a representative elastogram of the same sample shown in
Figure 1 IA where the sample shows essentially zero local strain except at the surface and interface of zone 1 and 2 with a scale ranging from -0.9% to 0.9%.
[00037] Figure 12A shows a representative displacement image of a femoral condyle cartilage sample showing a slight displacement gradient in the region of interest, using a
1/32" opening. The loading plate is indicated by an arrow. The scale ranges from -0. lmm to
0.1 mm.
[00038] Figure 12B shows a representative elastogram of the same sample shown in
Figure 12A where the sample shows essentially zero local strain except at the surface and interface of zone 1 and 2 with a scale ranging from -0.2% to 0.2%.
DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS
[00039] Understanding the mechanical behavior exhibited by joint structures and tendons is crucial in the diagnostic evaluation of conditions created by pathology and/or injury In addition, such understanding is necessary for the logical implementation of replacement grafts as well as in the design of biomimetic scaffold systems.
[00040] For example, understanding the mechanical behavior of native ACL and insertions provides valuable information concerning the native functional ligament interface at the junctions between the graft and bone. Currently, the long term integrity of
semitendinosus and hamstring tendon grafts used for ACL reconstruction is limited by healing of the graft with bone, which results in non-anatomical fibrovascular scar tissue at the interface between the tendon graft and. bone. The success of ACL grafts depends on reforming the native anatomical tendon-bone interface. Methods according to the present
invention allow for the determination of the structure- function relationship of joint structures, such as the ACL. This information then may be used to design scaffold systems, such as tendon grafts to bone, which mimic the native tissue in morphology, chemical composition, cellular distribution, and mechanical properties.
[00041] Methods according to the present invention are particularly useful for characterizing strain applied to structures associated with joints such as ligaments and cartilage, as well as tendons. Much of the description below and in the Examples section is directed to methods involving characterization of the anterior cruciate ligament and the characterization of articular cartilage of the femoral condyles. It should be understood, however, that these methods are readily applicable to other anatomical structures, for example, other ligaments, cartilage and tendons, for which detailed characterization of strain responses is desired.
[00042] These methods are especially appropriate where such structures are subject to non-uniform strain and may be used in vitro, in situ and under minimally invasive conditions. Minimally invasive conditions include, for example, endoscopic procedures. Methods according to the present invention also can be used for veterinary applications in, for example, equine, canine and feline species.
[00043] Methods according to the present invention are ideal for examining, for example, the ACL-bone interface as such methods permit the characterization of relatively small areas (on the order of about 0.1-2 mm, depending on the ultrasound frequency used) with complex stress distributions. In an exemplary embodiment, an ultrasound transducer scans a region of interest while an external load is applied to induce strain. Speckle tracking techniques may be employed to analyze the collected radio-frequency ultrasonic data before and after incremental loading and to estimate the resulting strain and strain distributions. Standard ultrasound scanners (e.g. Terason 2000, Teratech, Framingham MA) or similar
devices may be used. Ultrasound frequencies may range from, for example, 2-40 MHz, although higher or lower frequencies may be appropriate in certain circumstances, as will be appreciated by those of ordinary skill in the art. Sequences of RF data may be acquired during loading of the joint. Axial displacement between a reference and successive frames can be estimated using cross-correlation and recorrelation techniques, exemplary
embodiments of which are described in more detail below.
[00044] Axial displacements, or displacements occurring in the direction orthogonal to the face of the transducer and parallel to the direction of ultrasound propagation can be estimated for each RF frame with respect to a reference frame by using a ID cross-correlation algorithm. Strain distribution can be computed by differentiating the displacement map along the axial direction. For numerical differentiation, a least-squares regression method may be used. Displacement and strain then may be estimated relative to a reference frame to obtain a temporal profile and map of the cumulative deformation at the ligament and the insertion. These methods may be used to generate maps of cumulative deformation and strain in, for example, the ACL and tibial insertion during tensile loading. These techniques can be used to generate detailed information about displacements and strains associated with the tibiofemoral joint and other joints and tendons.
[00045] Cross-correlation and recorrelation algorithms may be used to obtain detailed information as described below. Figure 1 provides a representative sequence of cross- correlation and recorrelation algorithms. These techniques may be applied between consecutive echo (RF, envelope-detected or B-mode) axial (ID) segments in 2D (in-plane) or 3D (in-plane and out-of-plane). In a representative embodiment, each frame may be considered as a set of echo segments. As shown in Figures 2A and 2B each echo segment in a first frame may be cross-correlated with echo segments in a second frame to find the best match. After the best match is found between axial segments An (Figure 2A) and Bn (Figure
2B) the 2D or 3D resultant that denotes the path of motion between segments 1 and 2 may be broken into its individual components in the axial (in-plane, along the propagation axis; Δs in Figure 2), lateral (in-plane, orthogonal to the propagation axis; Δl in Figure 2) and elevational (out-of-planc). In its first iteration, the same technique is applied after rccorrclation, i.e., after correction, or removal, for the motion occurring in a direction perpendicular to that of estimation. This provides an increase of the correlation coefficient and, thus, increases the signal-to-noise ratio associated with the resulting motion estimate. As a result, the largest source of decorrelation is removed, which is typically in the direction of deformation. For example, if the direction of mechanical tension is within the ultrasound image plane but in the lateral direction, the lateral component may be removed first from original Frame 2 (i.e., frame 2_1 may be generated) for the estimation of the new axial component (axial estimate 2; Figure 1) and the second lateral component may be estimated (lateral estimate 2; Figure 1) after removal of the new axial component (i.e., frame 2_2 may be generated) and so on until reaching high enough correlation for the iterations to no longer be useful (Figure 1). This approach allows for effective decoupling of the 2D or 3D components and higher quality displacement and strain estimates and images of the tissue under deformation. It should be noted that the cross-correlation technique also may be employed without the use of recorrelation methods in the case where estimation of one or all components is considered to be of sufficient quality (e.g., high signal-to-noise ratio or high contrast-to-noise ratio).
[00046] The algorithms described immediately above can be modified in some circumstances for improved implementation of the invention, as will be appreciated by those of ordinary skill in the art. In addition, cross-correlation and recorrelation algorithms also may include those described in U.S. Patent No. 6,270,459, the contents of which are incorporated by reference in their entirety.
[00047] Tt also should be noted that non-RF data such as B-mode and envelope- detected data may be obtained and applied to methods according to the present invention. 2D and 3D displacement and strain estimation also may be obtained by utilizing methods described herein.
[00048] Information generated by methods according to the present invention may be used to inform the design and material selection for the production of tissue grafts.
Information on biodegradable scaffolds can be found in the literature, including Lu ct al., Biomatcrials 26 (2005) 4805-4816, the contents of which arc incorporated herein by reference in its entirety.
[00049] Methods according to the present invention allow for graft designs and material selection to take into consideration the detailed strain response for a particular anatomical structure, such as an ACL. Both graft design and material selection are very important for long term clinical success and involve a balance between scaffold
degradability, structural integrity needed for cell structure, overall scaffold mechanical properties, and the rate of cell matrix production.
[00050] Exemplary materials include those comprising poly-alpha-hydroxyesters such as polyglycoKc acid (PGA), poly-L-lactic acid (PLLA), and polylactic-co-glycolic copolymer (PLAGA), all of which have been approved by the FDA for these purposes. These types of degradable polymers do not elicit a permanent foreign body reaction and are gradually reabsorbed and replaced by natural tissue.
[00051] Protein modification of biomaterials may be used to improve cell adhesion and control the subsequent cellular response to material surfaces. Current strategies in improving cell attachment and augmenting subsequent cellular response include pre-coating these surfaces with molecules such as laminin, fϊbronectin (Fn) or grafting the Fn-related arginine- glycine-aspartic acid (RGD) tripeptide on biomaterials.
[00052] Tn addition to detailed characterization of ACL, the present invention also allows for high resolution ultrasound elastography of articular cartilage using the techniques described above and as discussed further in the Examples section below.
EXAMPLES
Example 1
Materials and Methods Tissue Isolation and Joint Preparation
[00053] Neonatal bovine calf (up to one week old) tibiofemoral joints obtained from an abattoir (Fresh Farm Beef, Vermont) were used for this Example. After removal of surrounding muscle and adipose tissue, the joint capsule was opened. Fascia lata and connective tissue were removed from the joint capsule with the ACL and posterior cruciate ligament (PCL) undisturbed. The PCL was maintained intact until immediately prior to testing in order to maintain joint stability and prevent premature damage to the ACL. During all joint preparation procedures, the ACL and surrounding tissues were kept hydrated with physiologic saline. The femur and tibia were cut to approximately 12 cm from the joint with a hacksaw, the periosteum removed, and bone marrow extracted from the intramedullary cavity to improve cement fixation of the joint. Subsequently, the tibia and femur were secured with custom anchors and cement to prevent slippage during testing. The joint was then mounted on a uniaxial material testing system (MTS 858 Bionix Testing System; MTS, Eden Prairie, MN) fitted with a custom cylindrical polycarbonate tank, the PCL was severed, and the medial femoral condyle was removed with a hacksaw to improve line-of-sight access to the ACL and insertions for the ultrasound transducer (Figures 3 A and 3B). These procedures may be applied in vivo for veterinary medicine purposes or in humans for assessing injury or age-related diseases such as osteoporosis.
Materials and Methods Tensile Testing
[00054J Tensile testing was performed with the femur-ACL-tibia complex (hereinafter
"FATC") in a tibial alignment following the methods of Woo et al. with modifications to accommodate ultrasound imaging. The femur and tibia were aligned along the tensile axis with 0° of knee flexion, and the sample was submerged in degassed physiologic saline. In addition to preventing tissue dehydration, the saline provided a medium for ultrasound propagation. A preload of 2 N was applied for one minute, and the joint was preconditioned by cyclic sawtooth loading from 0— 0.75 mm for 10 cycles at 20 rnm/min followed by a rest of 1 min. Three load regimens were applied to each sample (n = 3). First, the joint was cyclically loaded from 0 - 2 mm at 20 mm/min, with 0 mm being the displacement during the preload. Following a 30 minute rest, the joint was cyclically reloaded from 0 - 3 mm, with additional displacement applied during this testing regimen to ensure a detectable amount of deformation occurred across the insertions. Finally, after an additional 30 minute rest, the joint was loaded to failure at 10 mm/min.
Materials and Methods for Ultrasound Data Collection and Processing
[00055] While the joints were loaded in tension, an ultrasound scanner (Terason 2000;
Teratech, Inc., Rockville, MD) acquired Radio Frequency (RF) data at 5 MHz using a linear array. The ultrasound transducer was mounted inside the saline tank and positioned to image the ACL and insertions. Sequences of ultrasound RF data were acquired continuously during the applied loading repeatedly for periods of 3 seconds at 54 frames/s (128 RF lines, sampling frequency: 10 MHz). The axial displacement between a reference and successive frames was estimated offline and imaged using cross-correlation and recorrelation techniques with a window size of 3 mm and a window overlap of 80%.
[00056] Axial displacements, or displacements occurring in the direction orthogonal to the face of the transducer and parallel to the direction of ultrasound propagation, were
estimated for each RF frame with respect to a reference frame using a ID cross-correlation algorithm. In this algorithm, time-shifts between two backscattered signals are determined by the cross-correlation of small sliding windows over the entire 2D ultrasound image. At high decorrelation noise, recorrelation techniques were employed. Finally, the strain distribution was computed by differentiating the displacement map along the axial direction. For the numerical differentiation, a least-squares regression method was used. Displacement and strain were estimated relative to a reference frame, which was captured at the beginning of the application of tensile load, in order to obtain a temporal profile and map of the cumulative deformation at the ligament and the insertion.
Results Obtained
[00057] In all specimens tested, the ACL and the interface between the ACL and the femoral or tibial bone (Figure 4A) were readily identifiable on the ultrasound images (Figure 4B), as well as in all subsequent displacement maps and strain elastograms. Each tissue type exhibited a distinct ultrasound signature based on the structure and density of the tissue, which varies from soft to hard tissues and results in distinct speckle patterns distinguishing ligament, cartilage, and bone. As a highly dense tissue, bone is strongly echogenic resulting in bright ultrasound reflection and poor propagation of ultrasound waves through bone.
Ligament and bone are less dense and therefore less echogenic, enabling the structure of these tissues to be discernable on the B-mode images. Within the ACL, a narrow band of high strain in the middle and along the length of ACL was noted that also corresponded to a highly echogenic area on the B-scan images. This may reflect the parallel bundle- organization of the ACL. Distinctions between soft and hard tissue signatures were used to identify the ACL insertions into bone.
[00058] After establishing the ability to image the ACL and insertions with the experimental setup, mechanical testing was performed in order to obtain the mechanical
properties of the FATC using traditional mechanical testing means, as "well as to determine the localized mechanical behavior of the ACL and tibial insertion using ultrasound elastography. For neonatal bovine FATCs tested in a tibial orientation with 0° of knee flexion (Figure 3 A and B), the average stiffness was 59 ± 15 N/mm and the average modulus was 100 ± 30 MPa. Elastographic analysis yielded maps of cumulative deformation and strain in the ACL and tibial insertion during tensile loading. Figure 4C shows the distribution of deformation throughout the FATC, with magnitudes of deformation represented according to a colormap, with small deformations blue and large deformations red. The magnitude of displacement was found to be the highest (red in Figure 4C) within the ACL proper and decreased in value in a gradual transition (orange, yellow, and green) from ligament (red) to bone (blue). In addition, elastographic analysis revealed through strain maps that the strain profile at the tibial insertion was highly complex as the FATC was loaded in tension (Figure 4D). Both compressive and tensile strains were visualized at the tibial insertion site, indicated by the green-blue and yellow-red regions, respectively, on the elastogram in Figure 4D. In addition, compressive strains (blue) were found in the ACL itself, most likely because in this experiment, the ultrasound transducer was aligned with respect to the insertion (region of interest) during loading. A narrow region of blue appears along the entire upper tissue- saline interface in the elastogram in Figure 4D, which is an artifact of the high acoustical impedance difference between the FATC and saline and should not be confused with compressive strain.
[00059] Point-wise temporal displacement and strain analyses at the tibial ACL insertion demonstrate the deformation at the insertion over time with the development of both tensile and compressive strain. As shown in Figures 5A and 5B, which depicts the displacement and strain response, respectively, at a single point at the tibial insertion undergoing compressive strain (Figure 4D, arrow), a constant rate of deformation was
observed, with a corresponding non-constant evolution of compressive strain. The accumulation of deformation and strain throughout the entire FATC can be seen in Figure 6, which shows deformation (Figure 6A) and strain (Figure 6B) for a sequence of 16 frames with respect to a single reference frame. While the magnitude of deformation remains relatively constant in the bone (blue) over time from frames 1-16, the ACL proper undergoes high deformation, as shown be the change in color of the ACL from blue to green, yellow, orange, and finally red (Figure 6A). Additionally, Figure 6B reveals the evolution of strain over time, indicating the development of complex strain with both compressive (blue) and tensile (red) components at the interface between the ACL and bone.
[00060] In the elastographic analysis method used in this study, displacement and strains are measured along the axis parallel to the direction of ultrasound beam propagation. This may introduce artifacts corresponding to the orientation of the transducer with respect to the FATC. To ensure that results were not dependent on the orientation of the transducer, an additional trial was performed with the transducer rotated such that the face of the transducer was aligned along the principal axis of the ACL (Figure 7A). It is important to note that in this sample, a gradual decrease in the magnitude of deformation was again observed across the insertion, with high deformation (green) in the ACL, a low degree of deformation in bone (dark red), and a transition from green to yellow, orange and red at the interface between the ACL and bone (Figure 7B). Additionally, compressive strain (blue regions, Figure 7C) is again observed at the interface while the joint is loaded in tension. Tensile strains (red regions) are found at the same interface. These results are consistent with those presented in Figure 4 and combined collectively demonstrate the complexity of the strain profile at the tibial insertion.
[00061] These results demonstrate that displacement is non-uniformly distributed throughout the FATC and that strain in the ACL insertions is complex. Displacement images
reveal that deformation is higher in the ACL imidsubstance compared to the tibial insertion, and that a gradual transition exists in the degree of deformation from the ligament proper, through the tibial insertion, into bone. This distribution indicates a tissue type-dependent increase in stiffness progressing from ligament to interface and then to bone. Strain elastograms revealed that the strain distribution at the tibial insertion is highly complex, with both tensile and compressive strain components localized at the tibial insertion site. The complexity of strain distribution within the insertions may be due to the transfer of tensile strain from the ligament to bone through the intcrfacial fibrocartilagc tissue found at direct ligament and tendon insertions.
[00062] The results obtained by methods of the present invention and described herein constitute the first experimental determination of the complex strain distribution at ACL insertion sites and allow for the following observations: First, the presence of a
fibrocartilaginous transitional tissue between ligament and bone demonstrates that a compressive strain component exists in that region during physiological loading. Second, collagen fibers extending from ligament into bone at the insertions, when loaded in tension, transmit shear and compressive stresses through the fibrocartilage zones of the insertions.
[00063] These results allow for the determination of mechanical properties of the ACL and ACL-bone interface and represent a milestone in the understanding of the localized functionality of orthopaedic tissues, specifically at soft to hard tissue interfaces. The results obtained from the above-described methods also inform the engineering of interfacial scaffolds based on data derived from the analysis of the mechanical properties of healthy ACL insertions. In addition these methods allow for quantitative and noninvasive evaluation of the success of efforts to improve graft to bone healing.
[00064] These results demonstrate that these methods of ultrasound elastography provide valuable information on the mechanical behavior of the ACL insertions upon
physiological loading. This information allows for a better understanding of the structure- function relationship inherent at the ACL-bone interface, as well as tendon graft healing and re-establishment of a functional tendon-bone interface. In addition, mechanical parameters, such as Young's modulus, shear modulus and Poisson's ratio can be estimated and imaged based on the elastographic measurements.
Example 2
[00065] Techniques applied in Example 1 to anterior cruciate ligaments in vitro, arc applied in situ to characterize ACL ligaments. Highly detailed data characterizing strain responses of ACL ligaments arc obtained.
Example 3
[00066] Techniques applied in Example 1 to anterior cruciate ligaments are applied to other structures of the tibiofemoral joint, in vitro and in situ, including the posterior cruciate ligament, cartilage, and medial and lateral menisci. Highly detailed data characterizing strain responses of these structures are obtained.
Example 4
[00067] Techniques applied in Example 1 are applied to other joints of the upper and lower extremities including the foot, ankle, hip, temporomandibular joint (TMJ), shoulder, elbow, hand and wrist, in vitro and in situ. Highly detailed data characterizing strain responses of structures associated with these joints are obtained.
Example 5
[00068] Techniques applied in Example 1 are applied to tendons, such as the achilles tendon and flexor and extensor tendons of mammalian extremities, in vitro and in situ.
Highly detailed data characterizing strain responses of tendons are obtained.
Example 6
Materials and Methods
[00069J Full-thickness, 1-cm diameter cylindrical samples of articular cartilage (n=3) were obtained from bovine, femoral condyles with an average thickness of 4.71 ± 0.47 mm and femoral head with an average thickness of 3.39 ± 0.86 mm in immature, healthy bovine. The samples were immersed in PBS within a custom-made loading device, illustrated in Figure 8, mounted onto an Instron 5848 MicroTcstcr (Norwood, MA) equipped with a 100 N load cell (accuracy +/- 0.5N, with a 1 -micron displacement resolution of the actuator).
[00070] The specimens were oriented such that the deep portion of the cartilage contacted an aluminum loading plate and the articular surface rested upon another rigid, impermeable surface containing a 3 mm opening to serve as the acoustic window for the high-resolution ultrasound transducer (f/2, 8mm focus, 55 MHz, 46 Hz frame rate, Vevo 770, Visualsonics, Toronto, Canada). The ultrasound probe was separated by 3-mm from the surface of the articular cartilage.
[00071] A tare strain of 0.1% based on the measurement of the undeformed cartilage plugs was sustained for 30 seconds, followed by a ramp to strains ranging from 0.5 to 4.0% strain at 0.1 mm/sec for two femoral condyle samples and one femoral head sample, the results of which are graphically described in Figure 9. For these samples, B-mode ultrasound scans were acquired immediately after the tare strain and immediately after ramped compression of the cartilage. Subsequently, a pre-strain of 10% based on the measurement of the undeformed cartilage plugs was applied for 30 seconds, followed by a ramp to strains ranging from 2% at 0.1 mm/sec for two femoral head samples. Radio frequency (RF) ultrasound signals of these samples were acquired once equilibrium was attained.
Displacement images and elastograms were generated using ID crosscorrelation techniques and gradient operators on the RF signals (window size 0.3 mm, 85% overlap), respectively.
Median filtering of the displacement data was also implemented. This arrangement simulates a device in which an ultrasound transducer is incorporated into an arthroscopic indentation device.
Results Obtained
[00072] As shown in Figure 9, ramp compressions with high resolution displacement were applied and rapid stress relaxation of the articular cartilage samples recorded. As shown in Figures 1OA and 1OB, high quality RF signals and B-modc images were obtained of the cartilage plugs using the high resolution ultrasound with a 1/16" (Figure 10A) and 1/32" (Figure 10B) diameter opening in the loading plate. The full thickness of the cartilage sample can be observed and the articular layer can be discerned from zone 2. Abundant speckle is also visible within the cartilage, thus facilitating displacement and strain imaging.
[00073] As shown in Figure 1 IA, displacement images for all the samples imaged using the 1/8" opening in the loading plate revealed that while displacement of the region of interest did occur, a gradient of the displacements did not result. Rather, it appears that the entire region of the sample underwent rigid motion, displacing upwards through the opening as evidenced by circular indentations which were observed on the cartilage surfaces following the loading protocol. As shown in Figure 1 IB, the corresponding elastogram therefore displayed zero strain within most of the sample. Only the articular surface and the interface exhibited local strain. Similar maps were obtained regardless of the extent of post- compression strain.
[00074] Consistent with Figure 12A, the majority of displacement images for all the samples imaged using the 1/32" opening in the loading plate revealed that while displacement of the region of interest did occur, local displacements did not result. As shown in Figure 12B, the corresponding elastogram displayed about 0.5% strain within a portion of zone 2, while the articular surface and the interface exhibited strains on the order of 0.2%.
[00075] Whether the RF signals were captured immediately after the compressions or once equilibrium was achieved did not significantly impact the appearance of the cartilage on the elastograms. The elastograms were also similar between femoral condyle specimens and femoral head samples.
[00076] This Example confirms the usefulness of high resolution ultrasound elastographic imaging of articular cartilage, for example, for the early diagnosis and monitoring of treatment for articular cartilage pathologies, such as osteoarthritis.
Example 7
[00077] Techniques applied in Example 6 to articular cartilage in vitro, arc applied in situ. Highly detailed data images characterizing local strains to articular cartilage are obtained.
Example 8
[00078] Techniques applied in Example 6 are applied to articular cartilage of the joints of the upper and lower extremities including the ankle, hip, shoulder, elbow and wrist, in vitro and in situ. Highly detailed data images characterizing local strains of articular cartilage associated with these joints are obtained.
[00079] It will be understood that the foregoing description is illustrative of the principles of the invention, and that various modifications can be made by those skilled in the art without departing from the scope and spirit of the invention.
Claims
1. A method for obtaining information about the mechanical behaviour of structures associated with mammalian joints comprising:
(a) creating deformation in a joint structure of interest;
(b) using an ultrasound scanner in a linear array to acquire sequences of ultrasound data of the joint structure; and
(c) estimating a displacement component or displacement distribution between a reference frame of backscattcrcd signals and a successive frame of backscattcrcd signals, wherein said displacement component is estimated with a matching algorithm.
2. The method of claim 1, further comprising repeating step (c) until sufficient data is obtained to estimate strain distribution in the joint structure of interest.
3. The method of claim 2, wherein the joint structure of interest is selected from the group consisting of: ligaments and articular cartilage and any combination thereof.
4. The method of claim 3, wherein the matching algorithm includes determining time- shifts between two RF signals by cross-correlating sliding windows over a 2D or 3D ultrasound image to provide an estimation of axial, lateral or elevational displacement components.
5. The method of claim 4, further comprising computing a strain or strain rate distribution.
6. The method of claim 5, wherein computing a strain or strain rate distribution comprises differentiating a displacement map along one of the principle directions.
7. The method of claim 6 wherein the said differentiating comprises numerical differentiation, wherein the numerical differentiation includes least-square regression.
8. The method of claim 1 , wherein said estimation of the displacement component further includes a recorrelation algorithm.
9. The method of claim 1 or 8, further comprising the use of a window size between about 1 and about 5 mm and a window overlap between about 50 and about 99%.
10. The method of claim 1 , wherein the joint structure of interest is associated with the tibiofemoral joint.
11. The method of claim 10, wherein the ligament is an anterior cruciate ligament.
12. The method of claim 10, wherein the ligament is a posterior cruciate ligament.
13. The method of claim 1 , wherein the information is obtained in vitro.
14. The method of claim 1, wherein the information is obtained in vivo.
15. The method of claim 1 , wherein the information is obtained in situ noninvasively .
16. The method of claim 1, wherein the information is obtained in situ during a minimally invasive procedure such as an arthroscopy.
17. The method of claim 1, wherein the information is obtained in situ during an invasive procedure such as hip surgery.
18. The method of claim 1 , wherein the joint structure of interest is associated with a joint selected from the group consisting of: foot, ankle, knee, hip, hand, wrist, elbow, shoulder and temporomandibular joint (TMJ).
19. The method of claim 1, wherein the mammal is selected from the group consisting of: human, equine, canine and feline.
20. The method of claim 1 or 2, wherein said information is used to inform the design of tissue grafts.
21. A tissue graft produced using the information obtained by the method of claim 1 or 2.
22. A method for obtaining information about the mechanical behaviour of mammalian tendons comprising:
(a) creating deformation in a tendon of interest;
(b) using an ultrasound scanner with a piezoelectric element or an array of piezoelectric elements to acquire sequences of ultrasound data of the tendon; and
(c) estimating the displacement between a reference frame of the ultrasound data and successive frames of the ultrasound data, wherein said displacement or displacement distribution is estimated with a cross-correlation algorithm.
23. The method of claim 22, further comprising repeating step (c) until sufficient data is obtained to estimate strain distribution in the tendon of interest.
24. The method of claim 22, wherein the cross-correlation algorithm includes determining time-shifts between two backscattered RF signals by cross-correlating sliding windows over a 2D or 3D ultrasound image to provide an estimation of axial, lateral or elevational displacement components ultrasound image.
25. The method of claim 24, wherein the cross-correlation algorithm includes determining time-shifts between two backscattered RF signals by cross-correlating sliding windows over a 2D or 3D ultrasound image to provide an estimation of axial, lateral or elevational displacement components ultrasound image.
26. The method of claim 24, wherein the cross-correlation algorithm includes determining time-shifts between two backscattered RF signals by cross-correlating sliding windows over a 3D ultrasound image to provide an estimation of axial, lateral or elevational displacement components ultrasound image.
27. The method of claim 23, further comprising computing a strain or strain rate distribution.
28. The method of claim 27, wherein computing a strain or strain rate distribution comprises differentiating a displacement map along the axial direction.
29. The method of claim 28 wherein said differentiating comprises numerical
differentiation and wherein the numerical differentiation includes least-square regression.
30. The method of claim 22, wherein said estimation of axial displacement further includes a recorrelation algorithm.
31. The method of claim 23 or 28, further comprising the use of a window size between about 1 and 5mm and a window overlap between about 50 and 99%.
32. The method of claim 1 or 22, wherein said deformation is selected from the group consisting of tension, compression and relaxation, or any combination thereof.
33. The method of claim 1 or 22, wherein said deformation is generated by the ultrasound probe itself.
34. The method of claim 1 or 22, wherein said deformation is generated by the scanned subject itself.
35. The method of claim 22, wherein said piezoelectric element is embedded or is part of a surgical tool.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP06840170A EP1963805A4 (en) | 2005-12-09 | 2006-12-08 | Systems and methods for elastography imaging |
US12/096,254 US20090221916A1 (en) | 2005-12-09 | 2006-12-08 | Systems and Methods for Elastography Imaging |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US74924505P | 2005-12-09 | 2005-12-09 | |
US60/749,245 | 2005-12-09 | ||
US75112205P | 2005-12-16 | 2005-12-16 | |
US60/751,122 | 2005-12-16 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2007067987A2 true WO2007067987A2 (en) | 2007-06-14 |
WO2007067987A3 WO2007067987A3 (en) | 2007-12-06 |
Family
ID=38123656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2006/061809 WO2007067987A2 (en) | 2005-12-09 | 2006-12-08 | Systems and methods for elastography imaging |
Country Status (3)
Country | Link |
---|---|
US (1) | US20090221916A1 (en) |
EP (1) | EP1963805A4 (en) |
WO (1) | WO2007067987A2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010068450A1 (en) * | 2008-11-25 | 2010-06-17 | Mayo Foundation For Medical Education And Research | System and method for analyzing carpal tunnel using ultrasound imaging |
CN102361597A (en) * | 2010-01-20 | 2012-02-22 | 松下电器产业株式会社 | Displacement estimation method, displacement estimation device |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006044996A2 (en) * | 2004-10-15 | 2006-04-27 | The Trustees Of Columbia University In The City Of New York | System and method for automated boundary detection of body structures |
US10687785B2 (en) | 2005-05-12 | 2020-06-23 | The Trustees Of Columbia Univeristy In The City Of New York | System and method for electromechanical activation of arrhythmias |
WO2007035721A2 (en) * | 2005-09-19 | 2007-03-29 | The Trustees Of Columbia University In The City Of New York | Ultrasound method to open blood brain barrier |
GB0708358D0 (en) * | 2007-05-01 | 2007-06-06 | Cambridge Entpr Ltd | Strain image display systems |
US8502821B2 (en) * | 2008-02-04 | 2013-08-06 | C Speed, Llc | System for three-dimensional rendering of electrical test and measurement signals |
WO2011035312A1 (en) | 2009-09-21 | 2011-03-24 | The Trustees Of Culumbia University In The City Of New York | Systems and methods for opening of a tissue barrier |
WO2010014977A1 (en) | 2008-08-01 | 2010-02-04 | The Trustees Of Columbia University In The City Of New York | Systems and methods for matching and imaging tissue characteristics |
WO2010030819A1 (en) | 2008-09-10 | 2010-03-18 | The Trustees Of Columbia University In The City Of New York | Systems and methods for opening a tissue |
NZ616660A (en) * | 2011-03-25 | 2016-05-27 | Trackside Technologies Pty Ltd | Connective tissue monitoring, compositions for connective tissue treatment and methods for treating connective tissue |
US11000533B2 (en) | 2011-03-25 | 2021-05-11 | Trackside Technologies Pty Ltd. | Connective tissue monitoring, compositions for connective tissue treatment and methods for treating connective tissue |
US9320491B2 (en) | 2011-04-18 | 2016-04-26 | The Trustees Of Columbia University In The City Of New York | Ultrasound devices methods and systems |
WO2012162664A1 (en) | 2011-05-26 | 2012-11-29 | The Trustees Of Columbia University In The City Of New York | Systems and methods for opening of a tissue barrier in primates |
US9420997B2 (en) | 2012-06-14 | 2016-08-23 | Siemens Medical Solutions Usa, Inc. | Motion artifact suppression in ultrasound diagnostic imaging |
TWI584787B (en) * | 2012-09-10 | 2017-06-01 | 國立成功大學 | Living tendon detection apparatus and living tendon detection method |
WO2014059170A1 (en) | 2012-10-10 | 2014-04-17 | The Trustees Of Columbia University In The City Of New York | Systems and methods for mechanical mapping of cardiac rhythm |
US9247921B2 (en) | 2013-06-07 | 2016-02-02 | The Trustees Of Columbia University In The City Of New York | Systems and methods of high frame rate streaming for treatment monitoring |
US10322178B2 (en) | 2013-08-09 | 2019-06-18 | The Trustees Of Columbia University In The City Of New York | Systems and methods for targeted drug delivery |
US10028723B2 (en) | 2013-09-03 | 2018-07-24 | The Trustees Of Columbia University In The City Of New York | Systems and methods for real-time, transcranial monitoring of blood-brain barrier opening |
FR3013956B1 (en) * | 2013-12-04 | 2017-09-15 | Univ Rabelais Francois | IN VIVO EVALUATION METHOD OF THE PHYSIOPATHOLOGICAL STATE OF BIOLOGICAL TISSUE AND ASSOCIATED DEVICE |
US11202618B2 (en) * | 2015-09-21 | 2021-12-21 | Edan Instruments, Inc. | SNR improvement and operator-independence using time-varying frame-selection for strain estimation |
WO2018227088A1 (en) | 2017-06-08 | 2018-12-13 | Arizona Board Of Regents On Behalf Of Arizona State University | Systems and methods for generating an elastogram of brain tissue using mesoscopic wavelength ultrasound |
Family Cites Families (101)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3598111A (en) * | 1968-12-09 | 1971-08-10 | Health Technology Corp | Technique and apparatus for measuring and monitoring the mechanical impedance of body tissues and organ systems |
US4463608A (en) * | 1979-05-07 | 1984-08-07 | Yokogawa Hokushin Electric Corp. | Ultrasound imaging system |
US4777599A (en) * | 1985-02-26 | 1988-10-11 | Gillette Company | Viscoelastometry of skin using shear wave propagation |
US4822679A (en) * | 1985-08-26 | 1989-04-18 | Stemcor Corporation | Spray-applied ceramic fiber insulation |
US4858613A (en) * | 1988-03-02 | 1989-08-22 | Laboratory Equipment, Corp. | Localization and therapy system for treatment of spatially oriented focal disease |
US5107837A (en) * | 1989-11-17 | 1992-04-28 | Board Of Regents, University Of Texas | Method and apparatus for measurement and imaging of tissue compressibility or compliance |
US5457754A (en) * | 1990-08-02 | 1995-10-10 | University Of Cincinnati | Method for automatic contour extraction of a cardiac image |
US6405072B1 (en) * | 1991-01-28 | 2002-06-11 | Sherwood Services Ag | Apparatus and method for determining a location of an anatomical target with reference to a medical apparatus |
JP3109749B2 (en) * | 1991-04-17 | 2000-11-20 | 株式会社東芝 | Ultrasound imaging device |
WO1992020290A1 (en) * | 1991-05-17 | 1992-11-26 | Innerdyne Medical, Inc. | Method and device for thermal ablation |
JPH06511481A (en) * | 1991-07-05 | 1994-12-22 | ユニバーシティ オブ ロチェスター | Ultra-fine non-agglomerated porous particles that incorporate air bubbles |
EP0627206B1 (en) * | 1993-03-12 | 2002-11-20 | Kabushiki Kaisha Toshiba | Apparatus for ultrasound medical treatment |
US5435310A (en) * | 1993-06-23 | 1995-07-25 | University Of Washington | Determining cardiac wall thickness and motion by imaging and three-dimensional modeling |
US5601084A (en) * | 1993-06-23 | 1997-02-11 | University Of Washington | Determining cardiac wall thickness and motion by imaging and three-dimensional modeling |
US6259943B1 (en) * | 1995-02-16 | 2001-07-10 | Sherwood Services Ag | Frameless to frame-based registration system |
US5662113A (en) * | 1995-06-30 | 1997-09-02 | Siemens Medical Systems, Inc | Edge enhancement system for ultrasound images |
US6351659B1 (en) * | 1995-09-28 | 2002-02-26 | Brainlab Med. Computersysteme Gmbh | Neuro-navigation system |
US5606971A (en) * | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
US5810731A (en) * | 1995-11-13 | 1998-09-22 | Artann Laboratories | Method and apparatus for elasticity imaging using remotely induced shear wave |
WO1997032277A1 (en) * | 1996-02-29 | 1997-09-04 | Acuson Corporation | Multiple ultrasound image registration system, method and transducer |
JP3652791B2 (en) * | 1996-06-24 | 2005-05-25 | 独立行政法人科学技術振興機構 | Ultrasonic diagnostic equipment |
US6026173A (en) * | 1997-07-05 | 2000-02-15 | Svenson; Robert H. | Electromagnetic imaging and therapeutic (EMIT) systems |
US5752515A (en) * | 1996-08-21 | 1998-05-19 | Brigham & Women's Hospital | Methods and apparatus for image-guided ultrasound delivery of compounds through the blood-brain barrier |
US8003705B2 (en) * | 1996-09-23 | 2011-08-23 | Incept Llc | Biocompatible hydrogels made with small molecule precursors |
US6090800A (en) * | 1997-05-06 | 2000-07-18 | Imarx Pharmaceutical Corp. | Lipid soluble steroid prodrugs |
WO1998048783A1 (en) * | 1997-04-30 | 1998-11-05 | Point Biomedical Corporation | Microparticles useful as ultrasonic contrast agents and for delivery of drugs into the bloodstream |
US20020039594A1 (en) * | 1997-05-13 | 2002-04-04 | Evan C. Unger | Solid porous matrices and methods of making and using the same |
US6095976A (en) * | 1997-06-19 | 2000-08-01 | Medinol Ltd. | Method for enhancing an image derived from reflected ultrasound signals produced by an ultrasound transmitter and detector inserted in a bodily lumen |
US5928151A (en) * | 1997-08-22 | 1999-07-27 | Acuson Corporation | Ultrasonic system and method for harmonic imaging in three dimensions |
US6106465A (en) * | 1997-08-22 | 2000-08-22 | Acuson Corporation | Ultrasonic method and system for boundary detection of an object of interest in an ultrasound image |
JP4260920B2 (en) * | 1998-05-13 | 2009-04-30 | 株式会社東芝 | Ultrasonic diagnostic equipment |
US6287765B1 (en) * | 1998-05-20 | 2001-09-11 | Molecular Machines, Inc. | Methods for detecting and identifying single molecules |
US6511426B1 (en) * | 1998-06-02 | 2003-01-28 | Acuson Corporation | Medical diagnostic ultrasound system and method for versatile processing |
US6241675B1 (en) * | 1998-06-09 | 2001-06-05 | Volumetrics Medical Imaging | Methods and systems for determining velocity of tissue using three dimensional ultrasound data |
US6425867B1 (en) * | 1998-09-18 | 2002-07-30 | University Of Washington | Noise-free real time ultrasonic imaging of a treatment site undergoing high intensity focused ultrasound therapy |
US6246895B1 (en) * | 1998-12-18 | 2001-06-12 | Sunnybrook Health Science Centre | Imaging of ultrasonic fields with MRI |
US6309355B1 (en) * | 1998-12-22 | 2001-10-30 | The Regents Of The University Of Michigan | Method and assembly for performing ultrasound surgery using cavitation |
US6547730B1 (en) * | 1998-12-31 | 2003-04-15 | U-Systems, Inc. | Ultrasound information processing system |
FR2791136B1 (en) * | 1999-03-15 | 2001-06-08 | Mathias Fink | IMAGING METHOD AND DEVICE USING SHEAR WAVES |
US6186402B1 (en) * | 1999-03-18 | 2001-02-13 | Axiohm Transaction Solutions, Inc. | Credit and smart card readers |
US6352507B1 (en) * | 1999-08-23 | 2002-03-05 | G.E. Vingmed Ultrasound As | Method and apparatus for providing real-time calculation and display of tissue deformation in ultrasound imaging |
US7520856B2 (en) * | 1999-09-17 | 2009-04-21 | University Of Washington | Image guided high intensity focused ultrasound device for therapy in obstetrics and gynecology |
US6447450B1 (en) * | 1999-11-02 | 2002-09-10 | Ge Medical Systems Global Technology Company, Llc | ECG gated ultrasonic image compounding |
US6312382B1 (en) * | 1999-11-15 | 2001-11-06 | Ronald Mucci | Method and apparatus for extracting cardiac information from acoustic information acquired with an ultrasound device |
US6514221B2 (en) * | 2000-07-27 | 2003-02-04 | Brigham And Women's Hospital, Inc. | Blood-brain barrier opening |
US6508768B1 (en) * | 2000-11-22 | 2003-01-21 | University Of Kansas Medical Center | Ultrasonic elasticity imaging |
WO2002043564A2 (en) * | 2000-11-28 | 2002-06-06 | Allez Physionix Limited | Systems and methods for making non-invasive physiological assessments |
US6537221B2 (en) * | 2000-12-07 | 2003-03-25 | Koninklijke Philips Electronics, N.V. | Strain rate analysis in ultrasonic diagnostic images |
US6689060B2 (en) * | 2001-02-28 | 2004-02-10 | Siemens Medical Solutions Usa, Inc | System and method for re-orderable nonlinear echo processing |
US6537217B1 (en) * | 2001-08-24 | 2003-03-25 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for improved spatial and temporal resolution in ultrasound imaging |
FR2830936B1 (en) * | 2001-10-16 | 2004-08-27 | Agronomique Inst Nat Rech | METHOD FOR MEASURING THE STATE OF TENSION OF A MATERIAL AND APPLICATIONS OF THIS METHOD |
US6855114B2 (en) * | 2001-11-23 | 2005-02-15 | Karen Drukker | Automated method and system for the detection of abnormalities in sonographic images |
US7166075B2 (en) * | 2002-03-08 | 2007-01-23 | Wisconsin Alumni Research Foundation | Elastographic imaging of in vivo soft tissue |
CA2421468A1 (en) * | 2002-03-14 | 2003-09-14 | Matsushita Electric Industrial Co., Ltd. | Image processing device and ultrasonic diagnostic device |
US6683454B2 (en) * | 2002-03-28 | 2004-01-27 | Ge Medical Systems Global Technology Company, Llc | Shifting of artifacts by reordering of k-space |
US7549985B2 (en) * | 2002-06-26 | 2009-06-23 | The Regents Of The University Of Michigan | Method and system to create and acoustically manipulate a microbubble |
US6780152B2 (en) * | 2002-06-26 | 2004-08-24 | Acuson Corporation | Method and apparatus for ultrasound imaging of the heart |
US20040049134A1 (en) * | 2002-07-02 | 2004-03-11 | Tosaya Carol A. | System and methods for treatment of alzheimer's and other deposition-related disorders of the brain |
US7314446B2 (en) * | 2002-07-22 | 2008-01-01 | Ep Medsystems, Inc. | Method and apparatus for time gating of medical images |
US6749571B2 (en) * | 2002-09-19 | 2004-06-15 | Wisconsin Alumni Research Foundation | Method and apparatus for cardiac elastography |
US7103400B2 (en) * | 2002-11-08 | 2006-09-05 | Koninklijke Philips Electronics, N.V. | Artifact elimination in time-gated anatomical imaging |
US7697972B2 (en) * | 2002-11-19 | 2010-04-13 | Medtronic Navigation, Inc. | Navigation system for cardiac therapies |
AU2003303525A1 (en) * | 2002-12-31 | 2004-07-29 | Ultra-Sonic Technologies, L.L.C. | Transdermal delivery using encapsulated agent activated by ultrasound and/or heat |
US6994673B2 (en) * | 2003-01-16 | 2006-02-07 | Ge Ultrasound Israel, Ltd | Method and apparatus for quantitative myocardial assessment |
US7257244B2 (en) * | 2003-02-24 | 2007-08-14 | Vanderbilt University | Elastography imaging modalities for characterizing properties of tissue |
US20040172081A1 (en) * | 2003-02-28 | 2004-09-02 | Dai-Yuan Wang | Intracardiac pressure guided pacemaker |
US7344509B2 (en) * | 2003-04-17 | 2008-03-18 | Kullervo Hynynen | Shear mode therapeutic ultrasound |
US7175599B2 (en) * | 2003-04-17 | 2007-02-13 | Brigham And Women's Hospital, Inc. | Shear mode diagnostic ultrasound |
AU2004251360A1 (en) * | 2003-06-25 | 2005-01-06 | Siemens Medical Solutions Usa, Inc. | Automated regional myocardial assessment for cardiac imaging |
US6984209B2 (en) * | 2003-07-02 | 2006-01-10 | The Brigham And Women's Hospital, Inc. | Harmonic motion imaging |
US7055378B2 (en) * | 2003-08-11 | 2006-06-06 | Veeco Instruments, Inc. | System for wide frequency dynamic nanomechanical analysis |
US7358226B2 (en) * | 2003-08-27 | 2008-04-15 | The Regents Of The University Of California | Ultrasonic concentration of drug delivery capsules |
US20050054930A1 (en) * | 2003-09-09 | 2005-03-10 | The University Court Of The University Of Dundee | Sonoelastography using power Doppler |
EP1520517B1 (en) * | 2003-09-30 | 2005-12-14 | Esaote S.p.A. | Method for estimating tissue velocity vectors and tissue deformation from ultrasonic diagnostic imaging data |
US7421101B2 (en) * | 2003-10-02 | 2008-09-02 | Siemens Medical Solutions Usa, Inc. | System and method for local deformable motion analysis |
US7753847B2 (en) * | 2003-10-03 | 2010-07-13 | Mayo Foundation For Medical Education And Research | Ultrasound vibrometry |
US7896821B1 (en) * | 2003-11-14 | 2011-03-01 | Perfusion Technology, LLC | Low intensity directed ultrasound (LODUS) mediated blood brain barrier disruption |
EP1684724A4 (en) * | 2003-11-19 | 2008-04-02 | Barnes Jewish Hospital | Enhanced drug delivery |
EP1693005A4 (en) * | 2003-12-10 | 2010-09-01 | Panasonic Corp | Ultrasonograph and ultrasonography |
US7951083B2 (en) * | 2004-02-05 | 2011-05-31 | Siemens Medical Solutions Usa, Inc. | Motion analysis improvements for medical diagnostic ultrasound |
CN1997999B (en) * | 2004-03-29 | 2010-09-08 | 彼德·T·杰尔曼 | Systems and methods to determine elastic properties of materials |
US7809426B2 (en) * | 2004-04-29 | 2010-10-05 | The Cleveland Clinic Foundation | Acquiring contrast-enhanced, T1 weighted, cine magnetic resonance images |
EP1762180B1 (en) * | 2004-06-09 | 2015-04-15 | Hitachi Medical Corporation | Elastic image display method and ultrasonographic device |
CN101389273B (en) * | 2004-08-05 | 2012-09-05 | 贝勒研究院 | Gene or drug delivery system |
US20060058651A1 (en) * | 2004-08-13 | 2006-03-16 | Chiao Richard Y | Method and apparatus for extending an ultrasound image field of view |
US7678050B2 (en) * | 2004-08-24 | 2010-03-16 | General Electric Company | Method and apparatus for detecting cardiac events |
US20060074315A1 (en) * | 2004-10-04 | 2006-04-06 | Jianming Liang | Medical diagnostic ultrasound characterization of cardiac motion |
US7223241B2 (en) * | 2004-12-16 | 2007-05-29 | Aloka Co., Ltd. | Method and apparatus for elasticity imaging |
US8858441B2 (en) * | 2005-05-12 | 2014-10-14 | The Trustees Of Columbia University In The City Of New York | System and method for electromechanical wave imaging of body structures |
WO2007035721A2 (en) * | 2005-09-19 | 2007-03-29 | The Trustees Of Columbia University In The City Of New York | Ultrasound method to open blood brain barrier |
CA2628100C (en) * | 2005-11-02 | 2016-08-23 | Visualsonics Inc. | High frequency array ultrasound system |
EP1991129B1 (en) * | 2006-02-24 | 2016-06-29 | NanoVibronix Inc. | System for surface acoustic wave treatment of skin |
US8574157B2 (en) * | 2007-02-14 | 2013-11-05 | General Electric Company | Method and apparatus for generating an ultrasound image of moving objects using deformable models |
WO2008118992A1 (en) * | 2007-03-26 | 2008-10-02 | Boston Scientific Scimed, Inc. | High resolution electrophysiology catheter |
US20090191244A1 (en) * | 2007-09-27 | 2009-07-30 | Children's Medical Center Corporation | Microbubbles and methods for oxygen delivery |
US20090247911A1 (en) * | 2008-03-25 | 2009-10-01 | Petr Novak | Multiple-angle switched high intensity focused ultrasound |
EP2309931B1 (en) * | 2008-07-10 | 2013-12-11 | Koninklijke Philips N.V. | Ultrasonic assessment of cardiac synchronicity and viability |
WO2010014977A1 (en) * | 2008-08-01 | 2010-02-04 | The Trustees Of Columbia University In The City Of New York | Systems and methods for matching and imaging tissue characteristics |
US8380305B2 (en) * | 2010-07-01 | 2013-02-19 | DynaDx Corporation | System and method for predicting successful defibrillation for ventricular fibrillation cardiac arrest |
EP2600771A1 (en) * | 2010-08-06 | 2013-06-12 | The Trustees of Columbia University in the City of New York | Medical imaging contrast devices, methods, and systems |
WO2012162664A1 (en) * | 2011-05-26 | 2012-11-29 | The Trustees Of Columbia University In The City Of New York | Systems and methods for opening of a tissue barrier in primates |
-
2006
- 2006-12-08 EP EP06840170A patent/EP1963805A4/en not_active Withdrawn
- 2006-12-08 US US12/096,254 patent/US20090221916A1/en not_active Abandoned
- 2006-12-08 WO PCT/US2006/061809 patent/WO2007067987A2/en active Application Filing
Non-Patent Citations (1)
Title |
---|
REVELL ET AL.: "Computer Vision Elastography: Speckle adaptive Motion Estimation for Elastography Using Ultrasound Sequences", IEEE TRANSACTIONS ON MEDICAL IMAGING, vol. 24, no. 6, June 2005 (2005-06-01) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010068450A1 (en) * | 2008-11-25 | 2010-06-17 | Mayo Foundation For Medical Education And Research | System and method for analyzing carpal tunnel using ultrasound imaging |
US8795181B2 (en) | 2008-11-25 | 2014-08-05 | Mayo Foundation For Medical Education And Research | System and method for analyzing carpal tunnel using ultrasound imaging |
CN102361597A (en) * | 2010-01-20 | 2012-02-22 | 松下电器产业株式会社 | Displacement estimation method, displacement estimation device |
Also Published As
Publication number | Publication date |
---|---|
US20090221916A1 (en) | 2009-09-03 |
EP1963805A4 (en) | 2010-01-06 |
EP1963805A2 (en) | 2008-09-03 |
WO2007067987A3 (en) | 2007-12-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20090221916A1 (en) | Systems and Methods for Elastography Imaging | |
Spalazzi et al. | Elastographic imaging of strain distribution in the anterior cruciate ligament and at the ligament–bone insertions | |
Arndt et al. | Non-uniform displacement within the Achilles tendon during passive ankle joint motion | |
Olson et al. | Biomechanical consequences of fracture and repair of the posterior wall of the acetabulum. | |
Georgiou et al. | Accurate diagnosis of hip prosthesis loosening using a vibrational technique | |
Beaulé et al. | Surgical correction of cam deformity in association with femoroacetabular impingement and its impact on the degenerative process within the hip joint | |
Haen et al. | Shear waves elastography for assessment of human Achilles tendon's biomechanical properties: an experimental study | |
Kimura et al. | Reconstruction of a defect of the rotator cuff with polytetrafluoroethylene felt graft: recovery of tensile strength and histocompatibility in an animal model | |
Saarakkala et al. | Ultrasound indentation of normal and spontaneously degenerated bovine articular cartilage | |
JP2005504563A (en) | Methods and compositions for resurfacing joints | |
Huang et al. | A new method to evaluate the elastic modulus of cortical bone by using a combined computed tomography and finite element approach | |
Glaser et al. | A non-invasive acoustic and vibration analysis technique for evaluation of hip joint conditions | |
Laasanen et al. | Ultrasound indentation of bovine knee articular cartilage in situ | |
Gao et al. | Ultrasound elasticity imaging for determining the mechanical properties of human posterior tibial tendon: a cadaveric study | |
Nieminen et al. | Quantitative assessment of articular cartilage using high-frequency ultrasound: research findings and diagnostic prospects | |
Chan et al. | Functional MRI can detect changes in intratissue strains in a full thickness and critical sized ovine cartilage defect model | |
Khokhlova et al. | Assessment of hip and knee joints and implants using acoustic emission monitoring: A scoping review | |
Jeng et al. | Instrumented indentation investigation on the viscoelastic properties of porcine cartilage | |
Suydam et al. | Is echogenicity a viable metric for evaluating tendon properties in vivo? | |
Konofagou et al. | Elastographic imaging of the strain distribution at the anterior cruciate ligament and ACL-bone insertions | |
Khokhlova et al. | Non-invasive Assessment of Cartilage Damage of the Human Knee using Acoustic Emission Monitoring: a Pilot Cadaver Study | |
Thomas et al. | The biomechanical properties of fascia lata grafts: a preliminary study | |
Gielen et al. | Computed Tomography (CT) of the stifle | |
Grantham et al. | Ligaments | |
Spalazzi et al. | Elastographic imaging of strain distribution within the anterior cruciate ligament and at the acl-bone insertions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2006840170 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 12096254 Country of ref document: US |