Abstract:
An ultrasound machine processes ultrasonic data according to acoustoelastic properties of the materials to obtain strain information without specific assumptions with respect to the material properties of the measured material or a variety of different material properties normally not obtained by ultrasound machines.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
       [0001]    — 
       STATEMENT REGARDING FEDERALLY-SPONSORED RESEARCH OR DEVELOPMENT 
       [0002]    — 
       BACKGROUND OF THE INVENTION 
       [0003]    The present invention relates to ultrasonic quantitative measurements and, in particular, to an improved apparatus and method for making ultrasonic measurements in axial material properties when only lateral access to the material is possible. 
         [0004]    Conventional ultrasonic imaging provides a mapping of ultrasonic echo signals onto an image plane where the intensity of the echo, caused principally by relatively small differences in material properties between adjacent material types, is mapped to brightness of pixels on the image plane. While such images serve to distinguish rough structure within the body, they provide limited insight into the physical properties of the imaged materials. 
         [0005]    Ultrasound “elastography” may be used to produce data and images revealing material properties such as strain under an externally applied stress, Poisson&#39;s ratio, Young&#39;s modulus, and other common strain and strain-related measurements. Quasi-static elastography provides these measurements by taking two images of a material in two different states of compression, for example no compression and a given positive compression. The material may be compressed by a probe (including the transducer itself) or, for biological materials, by muscular action or movement of adjacent organs. Strain may be deduced from these two images by computing gradients of the relative shift of the material in the two images along the compression axis. Quasi-static elastography is analogous to a physician&#39;s palpation of tissue in which the physician determines stiffness by pressing the material and detecting the amount of material yield (strain) under this pressure. 
         [0006]    Ultrasound “acoustoelasticity” provides an alternative method of computing material properties using a mathematical description of the relationships between sound propagation and material properties to deduce the properties more directly. For example, a measurement of reflected energy as ultrasound passes through a front and rear boundary of a material in tension can be used to deduce the strain of the material along the axis of tension and the stiffness of the material along the direction of ultrasound propagation. This technique, pioneered by the present inventors, is described in U.S. patent application 2007/0089530 filed Oct. 16, 2006 and hereby incorporated by reference. 
         [0007]    The measurement of loads (stress) and stiffness in axial direction can be important in medicine, for example, in the assessment of disease or repair of ligaments and tendons. While the above described acoustoelastic technique allows the determination of axial strain (e.g. stretching along the length of the tendon) it does not measure the axial stress (e.g. tension along the length of the tendon) or axial material properties such as stiffness which may, particularly in tendons and ligaments, vary substantially between the axial and transverse directions. 
         [0008]    While it is generally possible to measure axial properties simply by shifting the direction of the transducers to be aligned along the axis of measurement, this is not practical with in situ tendon and ligaments. Many important applications where measurement of axial properties of tendons and ligaments is desired, occur when the tendons and ligaments are under tension in the body and only lateral access to the tendon or ligament is available, for example, through thin overlying tissue or a surgical incision. 
       BRIEF SUMMARY OF THE INVENTION 
       [0009]    The present inventors have developed a new apparatus and technique that can provide axial properties of the material when only lateral access is available for the ultrasound probes. Generally, the device provides a transverse directed ultrasonic transducer that determines transverse properties using the acoustoelastic techniques described above. A pair of angled transducers then provides an angled measurement of both transverse and axial properties. Mathematical combination of these two measurements can extract isolated axial properties. 
         [0010]    Specifically then, the present invention provides an ultrasound system for measuring axial material properties of material using ultrasound. The system includes an ultrasound transducer assembly transmitting an ultrasound signal into the material along a first transverse direction and receiving the ultrasound signal after reflection and modification by the material with a first and second tension applied along the axial extent of the material. In addition, the ultrasound transducer may transmit an ultrasound signal into the material along a first oblique direction at an angle between the transverse direction and the axial extent of the material and receive the ultrasound signal after reflection and modification by the material at the first and second tension. A processor receives and processes these signals along the transverse and angle directions according to a stored program to determine at least one of axial stress and axial stiffness of the material. 
         [0011]    It is thus one feature of at least one embodiment of the present invention to provide a measurement of axial properties of the material when only lateral access to the material can be obtained. 
         [0012]    The processor may determine transverse material properties from the ultrasonic signal along the transverse direction and may determine axial properties from the transverse properties by analysis of the time of flight of the ultrasound along the transverse and angle directions and the known geometry of the path of the ultrasound along the transverse and oblique direction. 
         [0013]    It is thus one feature of at least one embodiment of the present invention to deduce axial properties from transverse properties and a combination of transverse and axial time of flight of the ultrasonic signal. 
         [0014]    The processor may determine the transverse material properties from a time of flight of ultrasound between the first and second material positions and a reflection coefficient indicating reflected ultrasonic energy at the first and second material positions. The reflection coefficient may be a measure of a relative power reflected at the material positions in a predetermined frequency range. 
         [0015]    It is thus one feature of at least one embodiment of the present invention to make use of acoustoelastic analysis to determine the transverse properties of the material. 
         [0016]    The ultrasonic transducer assembly may include three ultrasonic transducers fitting within a housing having surfaces to align the housing with respect to an axial direction of the material wherein the first ultrasonic transducer directs ultrasonic energy substantially perpendicular to the tissue when the housing is aligned with the tissue and the second and third ultrasonic transducers direct energy at complementary obtuse and acute angles with respect to a perpendicular to the tissue when the housing is aligned with the tissue. 
         [0017]    It is thus one feature of at least one embodiment of the present invention to provide a simple apparatus for rapid determination of axial material properties, for example, in a surgical setting. 
         [0018]    The oblique direction is at an angle that is smaller than “critical angle” that will result in surface wave. 
         [0019]    It is thus one feature of at least one embodiment of the present invention to minimize interference resulting from a surface wave in the transverse direction. 
         [0020]    The ultrasound system may be surrounded by a transmission material and the first and second material positions are interfaces between the material and the transmission material. The interface material may be tissue or maybe a coupling elastomer of known material properties. 
         [0021]    It is thus one feature of at least one embodiment of the present invention to provide a system that may be used with ligaments and tendons surrounded by other tissue or isolated for measurement. 
         [0022]    These particular objects and advantages may apply to only some embodiments falling within the claims, and thus do not define the scope of the invention. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0023]      FIG. 1  is a simplified block diagram of an ultrasound scanner suitable for use with the present invention having a transducer system providing simultaneous transverse and oblique ultrasonic measurements; 
           [0024]      FIG. 2  is a simplified block diagram of a combination of signals from the two transducer directions to provide for axial measurements; 
           [0025]      FIG. 3  is an elevational cross-section of a probe for providing the simultaneous ultrasonic measurements of  FIG. 1 ; 
           [0026]      FIG. 4  is simplified representation of the tendon showing directions and variables used in the calculation of the present invention; 
           [0027]      FIG. 5  is a simplified representation of the transverse measurements used by the present invention; 
           [0028]      FIG. 6  is a flow chart of the steps of the present invention as executed by a processor following a stored program; 
           [0029]      FIGS. 7   a  and  7   b  are diagrams showing alternative transducer designs; and 
           [0030]      FIG. 8  is a diagram similar to that of  FIG. 7  showing an alternative embodiment of the transducer array of  FIG. 1  for evaluating two-dimensional or planar properties of materials. 
       
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
       [0031]    Referring now to  FIG. 1 , an acoustoelastographic ultrasound system  10  suitable for use with the present invention may employ an ultrasonic imaging machine  12  providing the necessary hardware and/or software to collect and process ultrasonic echo signals using a processor  14  executing a stored program  16  held within the ultrasonic imaging machine  12 . 
         [0032]    An ultrasonic transducer assembly  18  associated with the ultrasonic imaging machine  12  may have a housing  21  positionable about a ligament or tendon  20  extending an axial direction  26 . The ultrasonic transducer assembly  18  may include, in one embodiment, three ultrasonic transducers  24  in the housing to direct one ultrasonic beam generally perpendicular to the axial direction  26  of the tendon  20  and to project and receive a second ultrasonic beam at any angle that is smaller than “critical angle” (e.g., an angle above which the ultrasound passes predominately in a surface wave perpendicular to the to the axial direction  26 . Generally angle of 45° to the axial direction can be used. 
         [0033]    The echo signals from these beams may be received at the ultrasonic transducer assembly  18  to be converted to an electrical echo signal. 
         [0034]    Multiple acquisitions of echo signals may be obtained with the ligament  20  in different states of tension, for example, by muscular activity or by mechanical gripping of the tendon or similar technique. 
         [0035]    The electrical echo signals may be communicated along lead  28  to interface circuitry  30  of the ultrasonic imaging machine  12 . The interface circuitry  30  provides amplification, digitization, and other signal processing of the electrical signal as is understood in the art of ultrasonic imaging. The digitized echo signals are then transmitted to a memory  31  for storage (along with program  16 ) and subsequent processing by a processor  14 , as will be described below. 
         [0036]    After processing, the echo signals may be used to construct an image displayed on graphical display  32  and axial material properties may be displayed quantitatively on the graphical display  32  as well as conventional B-mode images. Input commands affecting the display of the echo signals and their processing may be received via a keyboard or cursor control device  34 , such as a mouse, attached to the processor  33  via interface circuitry  30 , as is well understood in the art. 
         [0037]    Referring now to  FIG. 2 , as noted above, a center transducer  24   b  of the ultrasonic transducer assembly  18  may be positioned within the housing  21  to direct a first beam perpendicularly to the axial direction  26 . This first beam passes through an outer transmissive material  40  into the tendon  20  through an interface between the two that causes a first reflection of ultrasonic energy back to the transducer  24   b . The wave continues until it strikes an interface between the tendon  20  and outer transmissive material  40  on the opposite side of the tendon  20  which also produces an echo reflecting altered sonic energy back to the transducer  24   b . The outer transmissive material  40  may be other tissue surrounding the tendon  20  or maybe compliant elastomeric pads when the apparatus is used surgically. In general, the signals from these echoes are processed to determine transverse properties of the tendon along an axis perpendicular to the axial direction  26 . Specifically these transfers properties as will be described below include a transverse stiffness value C 33  (where the 33 subscript indicates the transverse orientation) and a transverse strain value e 33  occurring between the two tensions applied to the tendon, as well as a transverse velocity value V 33  and a transverse time of flight T 33  being the time that it takes for the ultrasonic energy to pass between the two interfaces between outer transmissive material  40  and ligament  20 . The determination of these values will be described below and is indicated generally by process block  42 . 
         [0038]    Housing  21  also holds transducers  24   a  and  24   c  that flank transducer  24   b  and are tipped at approximate 45° true with respect to a normal to axial direction  26  to project a second ultrasonic beam from transducer  24   a  to be received by transducer  24   c  along this oblique path that has both transverse and axial dimensions. In this case at least a three echoes can be obtained from which it may deduce the time it takes the ultrasonic energy to cross into a ligament  20  from transducer  24   a  and to exit the ligament  20  proximate to transducer  24   c  after reflecting off the interface between outer transmissive material  40  and the ligament  20  opposite transducer  24   b . This time information is extracted as indicated by process block  44  and includes both transverse and axial information. Of principal interest is the time of flight T 13  (where the subscript  13  indicates both axial and transverse components) within the tendon  20  along this geometric path. 
         [0039]    In the present invention, information from process block  42  providing transverse material property may then be combined with the information from process block  44  providing transverse and axial information to produce axial information per process block  46  using a set of transducers that have only lateral access to the tendon  20 , that is from the side and not from the ends of the tendon as would be required using normal ultrasonic analysis for axial properties. 
         [0040]    Referring now to  FIG. 3 , in one embodiment it is envisioned that the housing  21  holding the ultrasonic transducer assembly  18  may provide a clamping mechanism  48  allowing the tendon  20  to be engaged, for example, by surgical opening and providing coupling outer transmissive material  40 , for example, a compliant elastomer, to cushion and surround the tendon  20  and provide ultrasonic coupling from the transducers  24  to the tendon  20 . 
       Obtaining Transverse Properties 
     Using Acoustoelastic Strain Gauge (ASG) 
       [0041]    Referring to  FIG. 6 , the program  16  initially collects the echo data along the transverse and lateral axes as indicated by process block  51 . Referring momentarily to  FIG. 5 , in the transverse direction transducer  24   b  will provide ultrasonic signal  50  directed transversely into the tendon  20  having a wave shape  52  at a first time t 0 . This ultrasonic signal  50  will produce two echoes, echo  54  from a closest interface between outer transmissive material  40  and the tendon  20  and echo  56  from a furthest interface between the outer transmissive material  40  and the tendon  20 . These echoes  54  and  56  will produce waveforms  58  and  60  at time t 1  and t 2  respectively. 
         [0042]    Referring now also to  FIG. 4 , the difference in time between t 1  and t 2 , termed (T 0     —     33 ) will be used to determine a transverse velocity V 33  as will be described below as well as transverse stiffness C 33  and transverse strain e 33  where the subscript  33  indicates the transverse direction. 
         [0043]    Similarly, in the oblique direction transducer  24   a  will provide an angled ultrasonic signal  64  having a wave shape  66  at time t 0  that will be received by transducer  24   c  at time t 1  as wave shape  68 . The difference between these two times t 1  and t 2 , termed (T 0     —     TOTAL ) will be used to deduce the axial velocity V 11  as will be described below as well as axial stiffness C 11  and axial strain e 11  where the “11” subscript indicates the axial direction. The actual time for the ultrasonic signals to pass through the tendon  20 , excluding the outer transmissive material  40 , may be either assumed to be the total time when the regions of the outer transmissive material  40  are small or may be deduced by looking at multiple echoes that indicate the interface crossings. 
         [0044]    Referring now to again to  FIG. 6  at process block  70  the transverse information described with respect to  FIG. 5  is used to extract time of flight and reflection coefficients for two states of tension of the ligament  20  as has previously been described with respect to the above referenced patent application 2007/0089530 now reproduced for completeness below. 
       Frequency Domain Analysis 
       [0045]    Strain e presented in sections hereafter represents the axial strain e 11 . Based on the assumption of tissue incompressibility (volume consistency through deformation), the strain in transverse (thickness) direction e 33  is related to axial strain e=e 11  by following relation. 
         [0000]    
       
         
           
             
               
                 
                   
                     e 
                     33 
                   
                   = 
                   
                     
                       1 
                       
                         
                           1 
                           + 
                           
                             e 
                             11 
                           
                         
                       
                     
                     = 
                     
                       1 
                       
                         
                           1 
                           + 
                           e 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   0 
                   ) 
                 
               
             
           
         
       
     
         [0046]    In addition, all material properties in all directions (axial and transverse) are treated as function of axial strain e. 
         [0047]    First, the transverse stiffness (e.g., along the transverse axis) is assumed to be some function of the applied axial strain e, for example, a general second order function such as: 
         [0000]        {tilde over (C)}   33 ( e )≈ c   2   ·e   2   +c   1   ·e+C   33 .  (1)
 
         [0048]    As will be discussed later, the order or the form of this function is not critical and is chosen so that it can be fit to the data. Further it is assumed that the tissue is nearly incompressible, that is, e 11  is a function of and can be deduced from e 33  and vice versa and so stiffness can be formulated with either axial or transverse strain as the dependent variable. This is true for other strain dependent phenomena, i.e. wave propagation velocity, wave travel time, and reflection coefficient which are defined below. From this point forward, it is assumed that all strain dependent phenomena are defined as a function of axial strain e with the subscripts dropped for simplicity. 
         [0049]    Because of this assumption of near incompressibility (i.e. the material density does not depend on the amount of strain or deformation), it follows that: 
         [0000]      ρ S ( e )≈ρ S ( e =0)=ρ 0   S .  (2)
 
         [0050]    Under these assumptions, the wave velocity through thickness (in transverse direction) is given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     V 
                      
                     
                       ( 
                       e 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                           
                             
                               C 
                               ~ 
                             
                             33 
                           
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                         
                           
                             ρ 
                             S 
                           
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                       
                     
                     ≈ 
                     
                       
                         
                           
                             
                               
                                 C 
                                 ~ 
                               
                               33 
                             
                              
                             
                               ( 
                               e 
                               ) 
                             
                           
                           
                             ρ 
                             0 
                             S 
                           
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0051]    The reflection coefficient at the interface of surrounding medium (outer transmissive material  40 ) and target medium (e.g. ligament  20 ) is defined by: 
         [0000]    
       
         
           
             
               R 
                
               
                 ( 
                 e 
                 ) 
               
             
             = 
             
               
                 
                   
                     ρ 
                     S 
                   
                    
                   
                     V 
                      
                     
                       ( 
                       e 
                       ) 
                     
                   
                 
                 - 
                 
                   
                     ρ 
                     W 
                   
                    
                   
                     V 
                     
                       W 
                        
                       
                           
                       
                     
                   
                 
               
               
                 
                   
                     ρ 
                     S 
                   
                    
                   
                     V 
                      
                     
                       ( 
                       e 
                       ) 
                     
                   
                 
                 + 
                 
                   
                     ρ 
                     W 
                   
                    
                   
                     V 
                     W 
                   
                 
               
             
           
         
       
     
         [0052]    or re-organized form 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ρ 
                       S 
                     
                     · 
                     
                       V 
                        
                       
                         ( 
                         e 
                         ) 
                       
                     
                   
                   = 
                   
                     
                       
                         1 
                         + 
                         
                           R 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                       
                       
                         1 
                         - 
                         
                           R 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                       
                     
                      
                     
                       
                         ρ 
                         W 
                       
                       · 
                       
                         
                           V 
                           W 
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
           
         
       
     
         [0053]    In these relations, superscripts S and W respectively represent the target tissue  20  (solid) and the surrounding medium (e.g., water equivalent) or known elastomer of outer transmissive material  40 . 
         [0054]    Combining equations (1)-(4), a relation between material properties (c 1 , c 2 , C 33 , ρ 0 ) and impedance square (right hand side of equation) for each stretched state is derived as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ρ 
                       0 
                       S 
                     
                      
                     
                       ( 
                       
                         
                           
                             c 
                             2 
                           
                           · 
                           
                             e 
                             2 
                           
                         
                         + 
                         
                           
                             c 
                             1 
                           
                           · 
                           e 
                         
                         + 
                         
                           C 
                           33 
                         
                       
                       ) 
                     
                   
                   = 
                   
                     
                       
                         [ 
                         
                           
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   R 
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                               
                               ) 
                             
                             
                               ( 
                               
                                 1 
                                 - 
                                 
                                   R 
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                            
                           
                               
                           
                            
                           
                             ρ 
                             W 
                           
                            
                           
                             V 
                             W 
                           
                         
                         ] 
                       
                       2 
                     
                     = 
                     
                       
                         IPS 
                          
                         
                           ( 
                           e 
                           ) 
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0055]    When the material properties of the outer transmissive material  40  are known, the reflection coefficient can be retrieved from measured echo signals, and the impedance square parameter on the right hand side is known. As will be noted below, however, knowledge of the material properties of the surrounding material is not required. 
         [0056]    Similarly, the relation between the material property and impedance square at non-stretched state can be derived as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ρ 
                       0 
                       S 
                     
                      
                     
                       C 
                       33 
                     
                   
                   = 
                   
                     
                       
                         [ 
                         
                           
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   R 
                                   0 
                                 
                               
                               ) 
                             
                             
                               ( 
                               
                                 1 
                                 - 
                                 
                                   R 
                                   0 
                                 
                               
                               ) 
                             
                           
                            
                           
                             ρ 
                             W 
                           
                            
                           
                             V 
                             W 
                           
                         
                         ] 
                       
                       2 
                     
                     = 
                     
                       
                         IPS 
                         0 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0057]    By the taking the ratio of reflection information equations (5) and (6), the unknown target tissue density cancels and is eliminated as a parameter required to determine the other material properties. This is important because in-vivo tissue density is very difficult to measure. 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             ρ 
                             0 
                             S 
                           
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                          
                         
                           
                             
                               C 
                               ~ 
                             
                             33 
                           
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                       
                       
                         
                           ρ 
                           0 
                           s 
                         
                          
                         
                           C 
                           33 
                         
                       
                     
                     ≈ 
                     
                       
                         
                           
                             
                               c 
                               2 
                             
                             · 
                             
                               e 
                               2 
                             
                           
                           + 
                           
                             
                               c 
                               1 
                             
                             · 
                             e 
                           
                         
                         
                           C 
                           33 
                         
                       
                       + 
                       1 
                     
                   
                   = 
                   
                     
                       
                         
                           [ 
                           
                             
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     R 
                                      
                                     
                                       ( 
                                       e 
                                       ) 
                                     
                                   
                                 
                                 ) 
                               
                               
                                 ( 
                                 
                                   1 
                                   - 
                                   
                                     R 
                                      
                                     
                                       ( 
                                       e 
                                       ) 
                                     
                                   
                                 
                                 ) 
                               
                             
                              
                             
                               ρ 
                               W 
                             
                              
                             
                               V 
                               W 
                             
                           
                           ] 
                         
                         2 
                       
                       
                         
                           [ 
                           
                             
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     R 
                                     0 
                                   
                                 
                                 ) 
                               
                               
                                 ( 
                                 
                                   1 
                                   - 
                                   
                                     R 
                                     0 
                                   
                                 
                                 ) 
                               
                             
                              
                             
                               ρ 
                               W 
                             
                              
                             
                               V 
                               W 
                             
                           
                           ] 
                         
                         2 
                       
                     
                     = 
                     
                       
                         
                           
                             [ 
                             
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     R 
                                      
                                     
                                       ( 
                                       e 
                                       ) 
                                     
                                   
                                 
                                 ) 
                               
                               
                                 ( 
                                 
                                   1 
                                   - 
                                   
                                     R 
                                      
                                     
                                       ( 
                                       e 
                                       ) 
                                     
                                   
                                 
                                 ) 
                               
                             
                             ] 
                           
                           2 
                         
                         
                           
                             [ 
                             
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     R 
                                     0 
                                   
                                 
                                 ) 
                               
                               
                                 ( 
                                 
                                   1 
                                   - 
                                   
                                     R 
                                     0 
                                   
                                 
                                 ) 
                               
                             
                             ] 
                           
                           2 
                         
                       
                       = 
                       
                         
                           IPS 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                         
                           IPS 
                           0 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
         [0058]    Here, the left hand side term represents the unknown normalized material property. Again, the right hand side is measured or known information. Since the material properties of the surrounding medium (ρ W V W ) are cancelled out in the same process, the material properties of the surrounding medium that were originally assumed to be known are also not essential to the analysis. 
         [0059]    Time Domain Analysis: 
         [0060]    The travel time at e=0 (non-stretched state or lesser stretched state) can be related to the unloaded tissue thickness D and wave velocity V 0  by: 
         [0000]    
       
         
           
             
               
                 
                   
                     T 
                     0 
                   
                   = 
                   
                     
                       
                         2 
                          
                         D 
                       
                       
                         V 
                         0 
                       
                     
                     = 
                     
                       2 
                        
                       D 
                        
                       
                         
                           
                             
                               ρ 
                               0 
                               s 
                             
                             
                               C 
                               33 
                             
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
           
         
       
     
         [0061]    Similarly, the travel time at e≠0 (stretched state or greater stretched state) is given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     T 
                      
                     
                       ( 
                       e 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         2 
                          
                         
                           d 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                       
                       
                         V 
                          
                         
                           ( 
                           e 
                           ) 
                         
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
           
         
       
     
         [0062]    Here, the wave velocity V and tissue thickness d at a stretched state are given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       V 
                       ≈ 
                       
                         
                           
                             
                               
                                 C 
                                 ~ 
                               
                               33 
                             
                              
                             
                                 
                             
                              
                             
                               ( 
                               e 
                               ) 
                             
                           
                           
                             ρ 
                             0 
                             s 
                           
                         
                       
                     
                     = 
                     
                       
                         
                           
                             
                               c 
                               2 
                             
                              
                             
                               e 
                               2 
                             
                           
                           + 
                           
                             
                               c 
                               1 
                             
                              
                             e 
                           
                           + 
                           
                             C 
                             33 
                           
                         
                         
                           ρ 
                           0 
                           s 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   and 
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
             
               
                 
                   
                     d 
                      
                     
                       ( 
                       e 
                       ) 
                     
                   
                   = 
                   
                     
                       D 
                       
                         
                           1 
                           + 
                           e 
                         
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
           
         
       
     
         [0063]    By taking the ratio of the wave travel time equations (10) and (11): 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             ( 
                             
                               
                                 T 
                                 0 
                               
                               
                                 T 
                                  
                                 
                                   ( 
                                   e 
                                   ) 
                                 
                               
                             
                             ) 
                           
                           2 
                         
                         = 
                           
                          
                         
                           
                             ( 
                             
                               
                                 2 
                                  
                                 
                                   D 
                                   / 
                                   
                                     V 
                                     0 
                                   
                                 
                               
                               
                                 2 
                                  
                                 
                                   d 
                                   / 
                                   V 
                                 
                               
                             
                             ) 
                           
                           2 
                         
                       
                     
                   
                   
                     
                       
                         = 
                           
                          
                         
                           
                             ( 
                             
                               
                                 
                                   
                                     
                                       
                                         C 
                                         ~ 
                                       
                                       33 
                                     
                                      
                                     
                                       ( 
                                       e 
                                       ) 
                                     
                                   
                                   
                                     C 
                                     33 
                                   
                                 
                               
                                
                               
                                 
                                   1 
                                   + 
                                   e 
                                 
                               
                             
                             ) 
                           
                           2 
                         
                       
                     
                   
                   
                     
                       
                         
                           = 
                             
                            
                           
                             
                               
                                 IPS 
                                  
                                 
                                   ( 
                                   e 
                                   ) 
                                 
                               
                               
                                 IPS 
                                 0 
                               
                             
                              
                             
                               ( 
                               
                                 1 
                                 + 
                                 e 
                               
                               ) 
                             
                           
                         
                         , 
                       
                     
                   
                 
               
               
                 
                   ( 
                   12 
                   ) 
                 
               
             
           
         
       
     
         [0064]    the strain e=e 11  applied to the tissue can be evaluated as: 
         [0000]    
       
         
           
             
               
                 
                   e 
                   = 
                   
                     
                       
                         
                           IPS 
                           0 
                         
                         
                           IPS 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                       
                        
                       
                         
                           ( 
                           
                             
                               T 
                               0 
                             
                             
                               T 
                                
                               
                                 ( 
                                 e 
                                 ) 
                               
                             
                           
                           ) 
                         
                         2 
                       
                     
                     - 
                     1. 
                   
                 
               
               
                 
                   ( 
                   13 
                   ) 
                 
               
             
           
         
       
     
         [0065]    Since the right hand side of this relation contains parameters that are measured directly from wave signals, the applied strain can be evaluated directly without any supplemental information. Once the applied strains are evaluated, they are used in equation (7) to compute the normalized stiffness: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             c 
                             2 
                           
                           · 
                           
                             e 
                             2 
                           
                         
                         + 
                         
                           
                             c 
                             1 
                           
                           · 
                           e 
                         
                       
                       
                         C 
                         33 
                       
                     
                     + 
                     1 
                   
                   = 
                   
                     
                       
                         
                           [ 
                           
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   R 
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                               
                               ) 
                             
                             
                               ( 
                               
                                 1 
                                 - 
                                 
                                   R 
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                           ] 
                         
                         2 
                       
                       
                         
                           [ 
                           
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   R 
                                   0 
                                 
                               
                               ) 
                             
                             
                               ( 
                               
                                 1 
                                 - 
                                 
                                   R 
                                   0 
                                 
                               
                               ) 
                             
                           
                           ] 
                         
                         2 
                       
                     
                     = 
                     
                       
                         
                           IPS 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                         
                           IPS 
                           0 
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
         [0066]    If the density of the target tissue ρ 0   S  is known prior to the testing, the coefficients in the transverse stiffness function can be directly evaluated by substituting the evaluated strain into relation (5): 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         c 
                         2 
                       
                       · 
                       
                         e 
                         2 
                       
                     
                     + 
                     
                       
                         c 
                         1 
                       
                       · 
                       e 
                     
                     + 
                     
                       C 
                       33 
                     
                   
                   = 
                   
                     
                       
                         
                           1 
                           
                             ρ 
                             0 
                             S 
                           
                         
                          
                         
                           [ 
                           
                             
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     R 
                                      
                                     
                                       ( 
                                       e 
                                       ) 
                                     
                                   
                                 
                                 ) 
                               
                               
                                 ( 
                                 
                                   1 
                                   - 
                                   
                                     R 
                                      
                                     
                                       ( 
                                       e 
                                       ) 
                                     
                                   
                                 
                                 ) 
                               
                             
                              
                             
                               ρ 
                               W 
                             
                              
                             
                               V 
                               W 
                             
                           
                           ] 
                         
                       
                       2 
                     
                     = 
                     
                       
                         1 
                         
                           ρ 
                           0 
                           S 
                         
                       
                        
                       
                         
                           IPS 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     5 
                     ′ 
                   
                   ) 
                 
               
             
           
         
       
     
         [0067]    Once the transverse stiffness is given as a function of applied axial strain e, the tissue thickness at any strain level can be calculated from relation: (9) 
         [0000]    
       
         
           
             
               
                 
                   
                     d 
                      
                     
                       ( 
                       e 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                           T 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                         2 
                       
                        
                       
                         V 
                          
                         
                           ( 
                           e 
                           ) 
                         
                       
                     
                     = 
                     
                       
                         
                           T 
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                         2 
                       
                        
                       
                         
                           
                             
                               
                                 
                                   c 
                                   2 
                                 
                                 · 
                                 
                                   e 
                                   2 
                                 
                               
                               + 
                               
                                 
                                   c 
                                   1 
                                 
                                 · 
                                 e 
                               
                               + 
                               
                                 C 
                                 33 
                               
                             
                             
                               ρ 
                               0 
                               S 
                             
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   14 
                   ) 
                 
               
             
           
         
       
     
         [0068]    The result of these calculations is to determine the transverse stiffness C 33 , and axial strain e (termed e in these equations) as indicated by process block  72  of  FIG. 6 . 
       Obtaining Axial Properties 
     Using Acoustoelastic Strain Gauge (ASG) Derived Values 
       [0069]    Using these derived transfer properties, the axial properties can now be determined. In the present invention, as noted above, we have previously determined e 1   
         [0000]    
       
         
           
             
               
                 
                   e 
                   = 
                   
                     
                       
                         
                           
                             C 
                             33 
                           
                           
                             
                               
                                 C 
                                 ~ 
                               
                               33 
                             
                              
                             
                               ( 
                               e 
                               ) 
                             
                           
                         
                          
                         
                           
                             ( 
                             
                               
                                 T 
                                 
                                   0 
                                    
                                   _ 
                                    
                                   33 
                                 
                               
                               
                                 
                                   T 
                                   33 
                                 
                                  
                                 
                                   ( 
                                   e 
                                   ) 
                                 
                               
                             
                             ) 
                           
                           2 
                         
                       
                       - 
                       1 
                     
                     = 
                     
                       
                         
                           
                             
                               [ 
                               
                                 
                                   ( 
                                   
                                     1 
                                     + 
                                     
                                       
                                         R 
                                         33 
                                       
                                        
                                       
                                           
                                       
                                        
                                       
                                         ( 
                                         e 
                                         ) 
                                       
                                     
                                   
                                   ) 
                                 
                                 
                                   ( 
                                   
                                     1 
                                     - 
                                     
                                       
                                         R 
                                         33 
                                       
                                        
                                       
                                         ( 
                                         e 
                                         ) 
                                       
                                     
                                   
                                   ) 
                                 
                               
                               ] 
                             
                             
                               2 
                                
                               
                                   
                               
                             
                           
                           
                             
                               [ 
                               
                                 
                                   ( 
                                   
                                     1 
                                     + 
                                     
                                       R 
                                       
                                         0 
                                          
                                         _ 
                                          
                                         33 
                                       
                                     
                                   
                                   ) 
                                 
                                 
                                   ( 
                                   
                                     1 
                                     - 
                                     
                                       R 
                                       
                                         0 
                                          
                                         _ 
                                          
                                         33 
                                       
                                     
                                   
                                   ) 
                                 
                               
                               ] 
                             
                             2 
                           
                         
                          
                         
                           
                             ( 
                             
                               
                                 T 
                                 
                                   0 
                                    
                                   _ 
                                    
                                   33 
                                 
                               
                               
                                 
                                   T 
                                   33 
                                 
                                  
                                 
                                   ( 
                                   e 
                                   ) 
                                 
                               
                             
                             ) 
                           
                           2 
                         
                       
                       - 
                       1 
                     
                   
                 
               
               
                 
                   
                     ( 
                     7 
                     ) 
                   
                    
                   
                     ( 
                     13 
                     ) 
                   
                 
               
             
           
         
       
     
         [0070]    and we have previously determined normalized C 33   
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             C 
                             ~ 
                           
                           33 
                         
                          
                         
                           ( 
                           e 
                           ) 
                         
                       
                       
                         C 
                         33 
                       
                     
                     ≈ 
                     
                       
                         
                           
                             c 
                             2 
                           
                           
                             C 
                             33 
                           
                         
                         · 
                         
                           e 
                           2 
                         
                       
                       + 
                       
                         
                           
                             c 
                             1 
                           
                           
                             C 
                             33 
                           
                         
                         · 
                         e 
                       
                       + 
                       1 
                     
                   
                   = 
                   
                     
                       
                         
                           [ 
                           
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   
                                     R 
                                     33 
                                   
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                               
                               ) 
                             
                             
                               ( 
                               
                                 1 
                                 - 
                                 
                                   
                                     R 
                                     33 
                                   
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                           ] 
                         
                         2 
                       
                       
                         
                           [ 
                           
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   R 
                                   
                                     0 
                                      
                                     _ 
                                      
                                     33 
                                   
                                 
                               
                               ) 
                             
                             
                               ( 
                               
                                 1 
                                 - 
                                 
                                   R 
                                   
                                     0 
                                      
                                     _ 
                                      
                                     33 
                                   
                                 
                               
                               ) 
                             
                           
                           ] 
                         
                         2 
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   15 
                   ) 
                 
               
             
           
         
       
     
         [0071]    This equation is a fundamental relationship in acoustoelastics: 
         [0000]    
       
         
           
             
               
                 
                   
                     V 
                     → 
                   
                   = 
                   
                     
                       
                         
                           V 
                           → 
                         
                         11 
                       
                       + 
                       
                         
                           V 
                           → 
                         
                         33 
                       
                     
                     = 
                     
                       
                         
                           
                             
                               
                                 
                                   
                                     
                                       C 
                                       ~ 
                                     
                                     11 
                                   
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                                 + 
                                 
                                   
                                     t 
                                     11 
                                   
                                    
                                   
                                     ( 
                                     e 
                                     ) 
                                   
                                 
                               
                                
                               
                                   
                               
                             
                             ρ 
                           
                            
                           
                             i 
                             → 
                           
                         
                       
                       + 
                       
                         
                           
                             
                               
                                 
                                   
                                     C 
                                     ~ 
                                   
                                   33 
                                 
                                  
                                 
                                   ( 
                                   e 
                                   ) 
                                 
                               
                               ρ 
                             
                              
                             
                               k 
                               → 
                             
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   16 
                   ) 
                 
               
             
           
         
       
     
         [0072]    In equation 16, t 11 (e) is axial stress as a function of axial strain. 
         [0073]    We don&#39;t know V directly because the geometry of the path may be affected by propagation velocity changes in the tissue and changes in stiffness according to the loading (the acoustoelastic effect) but we know normalized V 33  from perpendicular transducer and can deduce normalized V 11  as follows: 
         [0000]        T   0     —     TOTAL   2   =T   0     —     33   2   +T   0     —     11   2   −&gt;T   0     —     11 =√{square root over ( T   0     —     TOTAL   2   −T   0     —     33   2 )},  (17)
 
         [0000]        T   TOTAL ( e ) 2   =T   33 ( e ) 2   +T   11 ( e ) 2   −&gt;T   11 ( e )=√{square root over ( T   TOTAL ( e ) 2   −T   33 ( e ) 2 )}{square root over ( T   TOTAL ( e ) 2   −T   33 ( e ) 2 )}  (18)
 
         [0000]      and 
         [0000]        d   11 ( e )= D   11 ·(1 +e )  (19)
 
         [0074]    where the value of T 0     —     TOTAL  and T 0     —     33  is determined as described above. D 11  and d 11 (e) represent the axial length of a element of tissue under consideration at pre- and post-deformed states respectively. 
         [0075]    This gives us normalized V 11 , per process block  74 , as follows: 
         [0000]    
       
         
           
             
               
                 
                   e 
                   = 
                   
                     
                       
                         
                           
                             
                               
                                 V 
                                 11 
                               
                                
                               
                                 ( 
                                 e 
                                 ) 
                               
                             
                             
                               V 
                               
                                 0 
                                  
                                 _ 
                                  
                                 11 
                               
                             
                           
                            
                           
                             
                               
                                 T 
                                 11 
                               
                                
                               
                                 ( 
                                 e 
                                 ) 
                               
                             
                             
                               T 
                               
                                 0 
                                  
                                 _ 
                                  
                                 11 
                               
                             
                           
                         
                         - 
                         1 
                       
                       -&gt; 
                       
                         
                           
                             V 
                             11 
                           
                            
                           
                             ( 
                             e 
                             ) 
                           
                         
                         
                           V 
                           
                             0 
                              
                             _ 
                              
                             11 
                           
                         
                       
                     
                     = 
                     
                       
                         ( 
                         
                           e 
                           + 
                           1 
                         
                         ) 
                       
                        
                       
                           
                       
                        
                       
                         
                           
                             T 
                             
                               0 
                                
                               _ 
                                
                               11 
                             
                           
                           
                             
                               T 
                               11 
                             
                              
                             
                               ( 
                               e 
                               ) 
                             
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   20 
                   ) 
                 
               
             
           
         
       
     
         [0076]    Knowing C 33 , e and V 11  and V 33  we can determine C 11 +t 11  using equation (16). Separating axial stiffness and stress, C 11  and t 11  can be done for a linear material simply by knowing that C 11  e=t 11 . For a nonlinear material such as tendons it is necessary to integrate C 11  with respect to e to get t 11  as follows: 
         [0000]    
       
         
           
             
               
                 
                   
                     t 
                     11 
                   
                   = 
                   
                     
                       ∫ 
                       e 
                     
                      
                     
                       
                         C 
                         11 
                       
                        
                       
                         
                            
                           e 
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   21 
                   ) 
                 
               
             
           
         
       
     
         [0077]    In practice this integral can be resolved by defining C 11  as a function of e in terms of a, b, and c used as parameters for the nonlinear function, for example, coefficients of a second order polynomial. Stress t 11  is also defined by the same parameters plus the constant of integration which will be zero. Accordingly, knowledge of a, b, and c. which can be determined by empirical analysis of tendons, will permit the calculation of stiffness C 11  as a function of strain and also calculate stress t 11  as a function of strain per process block  76 . Stresses can be integrated over an area to compute axial force. 
         [0078]    Referring now to  FIG. 7   a  in an alternative embodiment of the ultrasonic transducer assembly  18 , a single ultrasonic transducer  24   a  may serve the role of both of transducers  24   a  and  24   b  described in  FIG. 2 . In this embodiment, transducer  24   a  is originally positioned at an angle less than the critical angle with respect to axial direction  26  for obtaining the oblique data and then pivoted as indicated by arrow  80  to a perpendicular direction to serve in the role of transducer  24   b  as described above. 
         [0079]    These methods could be used for materials other than biological tissues. The ultrasonic signals processed can be raw ultrasound or processed B-mode data. 
         [0080]    Alternatively, as shown in  FIG. 7   b , a phased array detector  82 , providing for multiple transmitting and receiving elements  84 , may be operated by appropriately phasing the transmit and receive signals among each of the elements  84  to direct a first beam  86  at an oblique angle and the second beam  87  at a perpendicular angle providing comparable data to that obtained by the ultrasonic transducer assembly  18  of  FIG. 2 . 
         [0081]    Referring now to  FIG. 8 , the present invention is not limited to deducing axial properties along a single dimension, that is, along a single dimension generally perpendicular to the lateral access of the transducer array  18 , but can be used to evaluate planar properties along axes at multiple angles generally perpendicular to the access direction of the transducer array  18 . 
         [0082]    In one embodiment, this may be done by establishing an array  18  of multiple oblique transducers  24   a ,  24 ′ a ,  24   c  and  24 ′ c  where transducers  24   a  and  24   c  are arranged as described above with respect to  FIG. 1  to direct ultrasound within a vertical plane aligned with axis  26 , and transducers  24 ′ a  and  24 ′ c  are arrayed similarly in a vertical plane but one that is generally perpendicular to the plane of transducers  24   a  and  24   c  and aligned with an axis  26 ′ perpendicular to axis  26 . 
         [0083]    This approach will permit the collection of oblique information along two perpendicular axes  26  and  26 ′ (with transducer  24   b  operating for both measurements to collect the perpendicular information), and by simple extension of the mathematics described above, to calculate axial properties in two dimensions within a plane of the material. 
         [0084]    It will be understood that this approach can be extended by using additional oblique transducer pairs at more angles than the two perpendicular angles of axes  26  and  26 ′ as shown in  FIG. 8 . Further it will be appreciated that a single transducer pair  24   a  and  24   c  can be rotated as indicated by arrow  90  about a perpendicular to the measured plane, while retaining their angles with respect to a planar surface of the material  20 . Alternatively, or in addition, the multi-transducer array having transducers  24   a ,  24 ′ a ,  24   c  and  24 ′ c  may be rotated, for example by increments of less than 90 degrees, for even faster data acquisition. 
         [0085]    Such a device as described above, may calculate planar stiffness in membranes like the pericardium, lung, heart valve leaflets, vocal folds, and fascia for compartment syndrome. The tissue  20  would need to be subject to loadings in multiple directions which may be accomplished by mechanical intervention or through the use of natural muscular motions by timing the measurements to take advantage of these periodic applied tensions. 
         [0086]    It is specifically intended that the present invention not be limited to the embodiments and illustrations contained herein, but include modified forms of those embodiments including portions of the embodiments and combinations of elements of different embodiments as come within the scope of the following claims.