ANISOTROPY AND OTHER CHARACTERISTICS

Related Applications

[0001] This application claims the benefit of provisional application serial number 62/336,437, filed May 1 3, 201 6, the disclosure of which is hereby incorporated herein by reference in its entirety.

Background of the Invention

[0002] Acoustic radiation force (ARF)-based ultrasound imaging

techniques are used to characterize the viscoelastic properties of tissue.

Some tissues are anisotropic, meaning their viscoelastic properties exhibit directional dependence. The ability to measure and image the degree of anisotropy may be diagnostically relevant because some pathologies can alter tissue structure and thereby the degree of anisotropy. Thus, the ratio of maximum elasticity to minimum elasticity is a diagnostically relevant tool to detect tissue pathology.

[0003] Conventional methods measure relative elasticity by measuring shear wave velocity. This is difficult because it requires taking multiple frames of an image (e.g., a movie) and trying to measure the propagation speed of the shear wave as it radiates out from the point of force. This computationally expensive process is made more difficult by the fact that the propagating wave can be channeled by structures such as arteries and veins and can be reflected at boundaries of changing elasticity such as the surface of bones.

[0004] Thus, there is a need for improved methods of taking

measurements of tissues and other materials, including methods and systems for assessing material anisotropy.

Summary of the Invention

[0005] This technology utilizes acoustic radiation force based ultrasound or other force-producing techniques to excite tissue or materials and uses the resulting displacements, which can be measured quite accurately, to estimate the degree of anisotropy of the tissue or material. The technology utilizes a standard ultrasound scanner and a one-dimensional or two-dimensional ultrasound transducer to perform these measurements by taking the ratio of the displacements in two directions. The degree of anisotropy may be relevant as a diagnostic metric for applications in imaging skeletal muscle, kidney, and cardiac muscle, all of which are known to be anisotropic, as well as other organ systems. To date, this technology has been simulated using finite element analysis in anisotropic material models with varying degrees of anisotropy, and has been clinically demonstrated in human muscle and kidney tissues.

[0006] We have previously shown that ARF-induced displacements allow for point measurements of degree of anisotropy, using finite element simulation. However, to make these measurements with a linear array transducer requires a 90 degree rotation of the imaging transducer or the sample, potentially leading to misalignment. Alternatively, we have demonstrated the feasibility of creating images of degree of anisotropy using ARF-induced displacements produced by a simulated two-dimensional (2D) matrix array transducer. With a 2D transducer, rotation is not needed. Other transducer geometries, such as a row-column design, would also support such anisotropy measurements without transducer rotation.

[0007] According to one aspect of the subject matter described herein, a method for assessing material anisotropy includes: applying a first force to a material sample and taking a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and calculating a material anisotropy based on the first and second measurements, wherein the first force and the second force are applied in any order.

[0008] In one embodiment, the angle difference between the first direction and the second direction is 90 degrees. [0009] In one embodiment, the angle difference between the first direction and the second direction is an angle other than 90 degrees.

[0010] In one embodiment, the first and second measurements comprise measurements of displacement, measurements of relaxation time constant for constant stress (Tau), measurements of recovery time, measurements of strain, measurements of elasticity, or measurements of viscosity.

[0011] In one embodiment, the material sample comprises a tissue sample.

[0012] In one embodiment, calculating a material anisotropy based on the first and second measurements comprises calculating the material anisotropy based on a ratio of the larger of two measurements to the smaller of two measurements.

[0013] In one embodiment, calculating a material anisotropy comprises calculating an elasticity anisotropy.

[0014] In one embodiment, calculating a material anisotropy comprises calculating a viscosity anisotropy.

[0015] In one embodiment, the coronal plane is at a focal depth within the material sample.

[0016] In one embodiment, the direction of anisotropy of the material is known and wherein the first force has an orientation at a known angle relative to the direction of anisotropy of the material.

[0017] In one embodiment, the direction of anisotropy of the material is unknown and wherein the method further comprises sweeping an angle of orientation of a long axis of force within the coronal plane until two different measurements are taken at a known difference in angle of orientation, and using the two different measurements and known difference in angle of orientation to calculate the material anisotropy.

[0018] In one embodiment, the direction of anisotropy of the material is unknown and wherein the method further comprises sweeping an angle of orientation of a long axis of force within the coronal plane until maximum and minimum measurements are taken, and using the maximum and minimum measurements to calculate the material anisotropy.

[0019] In one embodiment, at least one of the forces is applied using a mechanical indenter. [0020] In one embodiment, at least one of the forces is applied using a magnetic field.

[0021] In one embodiment, at least one of the forces is applied using an acoustic radiation force.

[0022] In one embodiment, the at least one force is applied using an ultrasound transducer.

[0023] In one embodiment, the first force is applied at a first focal depth, wherein the second force is applied at a second focal depth, and wherein the first focal depth is the same as the second focal depth.

[0024] In one embodiment, the first force is applied at a first focal depth, wherein the second force is applied at a second focal depth, and wherein the first focal depth is different than the second focal depth.

[0025] In one embodiment, the transducer's position and orientation are the same during application of the first force and the second force. (I.e., the position and orientation of the transducer does not change between application of the first force and application of the second force.)

[0026] In one embodiment, the ultrasound transducer comprises a lens.

[0027] In one embodiment, the ultrasound transducer comprises a plurality of emitter elements.

[0028] In one embodiment, the elements are arranged in a one- dimensional configuration.

[0029] In one embodiment, the one-dimensional configuration comprises a linear array.

[0030] In one embodiment, the direction of the long axis of the first force is controlled by rotation of the transducer about an axis normal to the coronal plane.

[0031] In one embodiment, the elements are arranged in a two- dimensional configuration.

[0032] In one embodiment, the two-dimensional configuration comprises one of: a fully-populated grid; a sparsely-populated grid; a checkerboard pattern; a tic-tac-toe board pattern; a cross shape; a square or rectangle; a circle, oval, or ellipse; or combinations of the above. [0033] In one embodiment, the direction of the long axis of the first force is controlled without rotation of the transducer about an axis normal to the coronal plane.

[0034] In one embodiment, at least one of the measurements is taken using ultrasound.

[0035] In one embodiment, at least one of the measurements is taken using magnetic resonance imaging (MRI).

[0036] In one embodiment, at least one of the measurements is taken using an optical method. In one embodiment, the optical method comprises optical coherence tomography (OCT).

[0037] In one embodiment, at least one of the measurements is taken using a mechanical micrometer.

[0038] In one embodiment, at least one of the measurements is taken using an optical microscope or other optical imaging or video method.

[0039] In one embodiment, the step of applying the first force and taking the first measurement is performed a plurality of times to produce a plurality of first measurements, wherein the step of applying the second force and taking the second measurement is applied a plurality of times, and wherein a material characteristic is calculated based on the plurality of first and second measurements. Examples of characteristics that may be so calculated include, but are not limited to, viscosity anisotropy, elasticity anisotropy, Tau values, and other metrics.

[0040] According to another aspect of the subject matter described herein, a system for assessing material anisotropy is adapted to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and calculate a material anisotropy based on the first and second measurements, wherein the first force and the second force are applied in any order.

[0041] According to another aspect of the subject matter described herein, a system for assessing material anisotropy comprises: means for applying a first force to a material sample and taking a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; means for applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and means for calculating a material anisotropy based on the first and second measurements, wherein the first force and the second force are applied in any order.

[0042] According to another aspect of the subject matter described herein, a system for assessing a material anisotropy comprises: a force applying module operable to apply, to a material sample, a plurality of forces having a direction and a coronal plane normal to the direction of the force; a

measurement module operable to measure, for each of the plurality of forces, a displacement of the material caused by the force; and a calculation module operable to calculate a material anisotropy based on the plurality of measured material displacements. The plurality of forces may comprise one or more forces in each of two or more directions. In one embodiment, for example, a force having one profile or orientation may be applied one or more times, and a force having another profile or orientation may be applied one or more times. In another embodiment, a force having a particular profile or orientation may be applied multiple times in a manner designed to detect, accommodate, or avoid a resonant property of the tissue.

[0043] According to another aspect of the subject matter described herein, a non-transitory computer readable medium storing software instructions that when executed by one or more processors of a system for assessing material anisotropy cause the system to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and calculate a material anisotropy based on the first and second measurements, wherein the first force and the second force are applied in any order.

[0044] According to another aspect of the subject matter described herein, a computer program comprising instructions which, when executed by at least one processor, cause the at least one processor to carry out any of the methods described herein.

[0045] According to another aspect of the subject matter described herein, a method for taking a material measurement comprises: applying a first force to a material sample and taking a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes (axes having equivalent length) within the coronal plane; and normalizing the first measurement using the second measurement to produce a first normalized measurement, wherein the first force and the second force are applied in any order.

[0046] In one embodiment, the angle difference between the first direction and the second direction is 90 degrees.

[0047] In one embodiment, the angle difference between the first direction and the second direction is an angle other than 90 degrees.

[0048] In one embodiment, the first and second measurements comprise: measurements of displacement; measurements of relaxation time constant for constant stress (Tau); measurements of recovery time; measurements of strain; measurements of elasticity; or measurements of viscosity.

[0049] In one embodiment, the material sample comprises a tissue sample.

[0050] In one embodiment, the method further comprises: applying a third force to a material sample and taking a third measurement, the third force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and normalizing the third measurement using the second measurement to produce a second normalized measurement, wherein the first force, the second force, and the third force are applied in any order.

[0051] In one embodiment, the method further comprises: applying a third force to a material sample and taking a third measurement, the third force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and applying a fourth force to the material sample and taking a fourth measurement, the fourth force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and normalizing the third measurement using the fourth measurement to produce a second normalized measurement, wherein the first force, the second force, the third force, and the fourth force are applied in any order.

[0052] In one embodiment, the method further comprises calculating a material anisotropy based on the first normalized measurement and the second normalized measurement.

[0053] In one embodiment, calculating a material anisotropy based on the first and second measurements comprises calculating the material anisotropy based on a ratio of the larger of two measurements to the smaller of two measurements.

[0054] In one embodiment, calculating a material anisotropy comprises calculating an elasticity anisotropy. [0055] In one embodiment, calculating a material anisotropy comprises calculating a viscosity anisotropy.

[0056] In one embodiment, the coronal plane is at a focal depth within the material sample.

[0057] In one embodiment, the direction of anisotropy of the material is known and wherein the first force has an orientation at a known angle relative to the direction of anisotropy of the material.

[0058] In one embodiment, the direction of anisotropy of the material is unknown and wherein the method further comprises sweeping an angle of orientation of a long axis of force within the coronal plane until two different measurements are taken at a known difference in angle of orientation, and using the two different measurements and known difference in angle of orientation to calculate the material anisotropy.

[0059] In one embodiment, the direction of anisotropy of the material is unknown and wherein the method further comprises sweeping an angle of orientation of a long axis of force within the coronal plane until maximum and minimum measurements are taken, and using the maximum and minimum measurements to calculate the material anisotropy.

[0060] In one embodiment, at least one of the forces is applied using a mechanical indenter.

[0061] In one embodiment, at least one of the forces is applied using a magnetic field.

[0062] In one embodiment, at least one of the forces is applied using an acoustic radiation force.

[0063] In one embodiment, the at least one force is applied using an ultrasound transducer.

[0064] In one embodiment, the first force is applied at a first focal depth, wherein the second force is applied at a second focal depth, and wherein the first focal depth is the same as the second focal depth.

[0065] In one embodiment, the first force is applied at a first focal depth, wherein the second force is applied at a second focal depth, and wherein the first focal depth is different than the second focal depth.

[0066] In one embodiment, the transducer's position and orientation are the same during application of the first force and the second force. (I.e., the position and orientation of the transducer does not change between application of the first force and application of the second force.)

[0067] In one embodiment, the ultrasound transducer comprises a lens.

[0068] In one embodiment, the ultrasound transducer comprises a plurality of emitter elements.

[0069] In one embodiment, the elements are arranged in a one- dimensional configuration.

[0070] In one embodiment, the one-dimensional configuration comprises a linear array.

[0071] In one embodiment, the direction of the long axis of the first force is controlled by rotation of the transducer about an axis normal to the coronal plane.

[0072] In one embodiment, the elements are arranged in a two- dimensional configuration.

[0073] In one embodiment, the two-dimensional configuration comprises one of: a fully-populated grid; a sparsely-populated grid; a checkerboard pattern; a tic-tac-toe board pattern; a cross shape; a square or rectangle; a circle, oval, or ellipse; or combinations of the above.

[0074] In one embodiment, the direction of the long axis of the first force is controlled without rotation of the transducer about an axis normal to the coronal plane.

[0075] In one embodiment, at least one of the measurements is taken using ultrasound.

[0076] In one embodiment, at least one of the measurements is taken using MRI.

[0077] In one embodiment, at least one of the measurements is taken using an optical method. In one embodiment, the optical method comprises OCT.

[0078] In one embodiment, at least one of the measurements is taken using a mechanical micrometer.

[0079] In one embodiment, at least one of the measurements is taken using an optical microscope.

[0080] In one embodiment, the step of applying the first force and taking the first measurement is performed a plurality of times to produce a plurality of first measurements, wherein the step of applying the second force and taking the second measurement is applied a plurality of times, and wherein a material viscosity anisotropy is calculated based on the plurality of first and second measurements.

[0081] According to another aspect of the subject matter described herein, a system for taking a material measurement is adapted to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and normalize the first measurement using the second measurement to produce a first normalized measurement, wherein the first force and the second force are applied in any order.

[0082] According to another aspect of the subject matter described herein, a system for taking a material measurement comprises: means for applying a first force to a material sample and taking a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; means for applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and means for normalizing the first measurement using the second measurement to produce a first normalized measurement, wherein the first force and the second force are applied in any order.

[0083] According to another aspect of the subject matter described herein, a system for taking a material measurement comprises: a force applying module operable to apply, to a material sample, a plurality of forces having a direction and a coronal plane normal to the direction of the force; a

measurement module operable to take, for each of the plurality of forces, a material measurement; and a calculation module operable to normalize one of the plurality of measurements using another of the plurality of measurements to produce a normalized measurement.

[0084] According to another aspect of the subject matter described herein, a non-transitory computer readable medium storing software instructions that when executed by one or more processors of a system for taking a material measurement cause the system to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and normalize the first measurement using the second measurement to produce a first normalized measurement, wherein the first force and the second force are applied in any order.

[0085] According to another aspect of the subject matter described herein, a computer program comprising instructions which, when executed by at least one processor, cause the at least one processor to carry out any of the methods described herein.

[0086] According to another aspect of the subject matter described herein, a method for taking a material measurement comprises: applying a first force to a material sample and taking a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and assessing a material property based on a difference between the first measurement and the second measurement, wherein the first force and the second force are applied in any order. [0087] In one embodiment, the angle difference between the first direction and the second direction is 90 degrees.

[0088] In one embodiment, the angle difference between the first direction and the second direction is an angle other than 90 degrees.

[0089] In one embodiment, the first and second measurements comprise: measurements of displacement; measurements of relaxation time constant for constant stress (Tau); measurements of recovery time; measurements of strain; measurements of elasticity; or measurements of viscosity.

[0090] In one embodiment, the material sample comprises a tissue sample.

[0091] In one embodiment, assessing a material property based on a difference between the first measurement and the second measurement comprises determining whether a difference between the first and second measurements exceeds a threshold.

[0092] In one embodiment, the method comprises, in response to determining that the difference between the first and second measurements exceeds the threshold, determining a change in tissue pathology.

[0093] According to another aspect of the subject matter described herein, a system for taking a material measurement is adapted to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and determine a material property based on a difference between the first measurement and the second measurement, wherein the first force and the second force are applied in any order.

[0094] According to another aspect of the subject matter described herein, a system for taking a material measurement comprises: means for applying a first force to a material sample and taking a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; means for applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and means for assessing a material property based on a difference between the first measurement and the second measurement, wherein the first force and the second force are applied in any order.

[0095] According to another aspect of the subject matter described herein, a system for taking a material measurement comprises: a force applying module operable to apply, to a material sample, a plurality of forces having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; a measurement module operable to take, for each of the plurality of forces, a material measurement; and a calculation module operable to determine a material property based on the plurality of material

measurements.

[0096] According to another aspect of the subject matter described herein, a non-transitory computer readable medium storing software instructions that when executed by one or more processors of a system for taking a material measurement cause the system to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and determine a material property based on a difference between the first measurement and the second measurement, wherein the first force and the second force are applied in any order.

[0097] According to another aspect of the subject matter described herein, a computer program comprising instructions which, when executed by at least one processor, cause the at least one processor to carry out any of the methods described herein. [0098] According to another aspect of the subject matter described herein, a system for taking a material measurement comprises: a force producing unit for applying a force to a material sample; a measuring device for taking a material measurement corresponding to the forces; one or more processors; and memory storing instructions executable by the one or more processors, whereby the system is operable to determine a material property based on a plurality of applied forces and their corresponding plurality of material measurements.

[0099] In one embodiment, at least one of the material measurements comprise: a measurement of displacement; a measurement of relaxation time constant for constant stress (Tau); a measurement of recovery time; a measurement of strain; a measurement of elasticity; or a measurement of viscosity.

[0100] In one embodiment, the material sample comprises a tissue sample.

[0101] In one embodiment, via execution of instructions by the one or more processors, the system is operable to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and calculate a material anisotropy based on the first and second measurements, wherein the first force and the second force are applied in any order.

[0102] In one embodiment, via execution of instructions by the one or more processors, the system is operable to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and normalize the first measurement using the second measurement to produce a first normalized measurement, wherein the first force and the second force are applied in any order.

[0103] In one embodiment, via execution of instructions by the one or more processors, the system is further operable to: apply a third force to a material sample and take a third measurement, the third force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and normalize the third measurement using the second measurement to produce a second normalized measurement, wherein the first force, the second force, and the third force are applied in any order.

[0104] In one embodiment, via execution of instructions by the one or more processors, the system is further operable to: apply a third force to a material sample and take a third measurement, the third force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force; and apply a fourth force to the material sample and take a fourth measurement, the fourth force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and normalize the third measurement using the fourth measurement to produce a second normalized measurement, wherein the first force, the second force, the third force, and the fourth force are applied in any order.

[0105] In one embodiment, via execution of instructions by the one or more processors, the system is operable to: apply a first force to a material sample and take a first measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; apply a second force to the material sample and take a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane; and determine a material property based on a difference between the first measurement and the second measurement, wherein the first force and the second force are applied in any order.

[0106] In one embodiment, determining a material property based on a difference between the first measurement and the second measurement comprises determining whether a difference between the first and second measurements exceeds a threshold.

[0107] In one embodiment, the material sample comprises a tissue sample and wherein, in response to determining that the difference between the first and second measurements exceeds a threshold, determining a change in tissue pathology.

[0108] In one embodiment, via execution of instructions by the one or more processors, the system is operable to calculate a material anisotropy based on the first normalized measurement and the second normalized

measurement.

[0109] In one embodiment, calculating a material anisotropy comprises calculating an elasticity anisotropy.

[0110] In one embodiment, calculating a material anisotropy comprises calculating a viscosity anisotropy.

[0111] In one embodiment, the coronal plane is at a focal depth within the material sample.

[0112] In one embodiment, the direction of anisotropy of the material is known and wherein the first force has an orientation at a known angle relative to the direction of anisotropy of the material.

[0113] In one embodiment, the direction of anisotropy of the material is unknown and wherein, via execution of instructions by the one or more processors, the system is operable to sweep an angle of orientation of a long axis of force within the coronal plane until two different measurements are taken at a known difference in angle of orientation, and use the two different measurements and known difference in angle of orientation to calculate the material anisotropy. [0114] In one embodiment, the direction of anisotropy of the material is unknown and wherein, via execution of instructions by the one or more processors, the system is operable to sweep an angle of orientation of a long axis of force within the coronal plane until maximum and minimum

measurements are taken, and use the maximum and minimum

measurements to calculate the material anisotropy.

[0115] In one embodiment, the force producing unit comprises a mechanical indenter.

[0116] In one embodiment, the force producing unit produces a magnetic field.

[0117] In one embodiment, the force producing unit comprises an ultrasound transducer.

[0118] In one embodiment, the first force is applied at a first focal depth, wherein the second force is applied at a second focal depth, and wherein the first focal depth is the same as the second focal depth.

[0119] In one embodiment, the first force is applied at a first focal depth, wherein the second force is applied at a second focal depth, and wherein the first focal depth is different than the second focal depth.

[0120] In one embodiment, the transducer's position and orientation are the same during application of the first force and the second force. (I.e., the position and orientation of the transducer does not change between application of the first force and application of the second force.)

[0121] In one embodiment, the ultrasound transducer comprises a lens.

[0122] In one embodiment, the ultrasound transducer comprises a plurality of emitter elements.

[0123] In one embodiment, the elements are arranged in a one- dimensional configuration.

[0124] In one embodiment, the one-dimensional configuration comprises a linear array.

[0125] In one embodiment, the ultrasound transducer produces a force having an oval or other profile with long and short axes within the coronal plane and wherein the direction of the long axis of the first force is controlled by rotation of the transducer about an axis normal to the coronal plane. [0126] In one embodiment, the elements are arranged in a two- dimensional configuration.

[0127] In one embodiment, the two-dimensional configuration comprises one of: a fully-populated grid; a sparsely-populated grid; a checkerboard pattern; a tic-tac-toe board pattern; a cross shape; a square or rectangle; a circle, oval, or ellipse; or combinations of the above.

[0128] In one embodiment, the ultrasound transducer produces a force having an oval or other profile with long and short axes within the coronal plane and wherein the direction of the long axis of the first force is controlled without rotation of the transducer about an axis normal to the coronal plane.

[0129] In one embodiment, at least one of the measurements is taken using ultrasound.

[0130] In one embodiment, at least one of the measurements is taken using MRI.

[0131] In one embodiment, at least one of the measurements is taken using an optical method. In one embodiment, the optical method comprises OCT.

[0132] In one embodiment, at least one of the measurements is taken using a mechanical micrometer.

[0133] In one embodiment, at least one of the measurements is taken using an optical microscope.

[0134] In one embodiment, via execution of instructions by the one or more processors, the system is operable to: apply a first force and take a first measurement a plurality of times to produce a plurality of first measurements; apply a second force and take a second measurement a plurality of times to produce a plurality of second measurements; and calculate a material viscosity based on the plurality of first and second measurements.

[0135] Those skilled in the art will appreciate the scope of the present invention and realize additional aspects thereof after reading the following detailed description of the preferred embodiments in association with the accompanying drawing figures. Brief Description of the Drawing Figures

[0136] The accompanying drawing figures incorporated in and forming a part of this specification illustrate several aspects of the invention, and together with the description serve to explain the principles of the invention.

[0137] Figure 1 A illustrates an ultrasound transducer according to an embodiment of the subject matter described herein.

[0138] Figure 1 B illustrates the ultrasound transducer in use over a material sample, showing an oval force profile produced within the material sample.

[0139] Figure 2 illustrates a flow chart for an exemplary process for assessing material anisotropy according to an embodiment of the subject matter described herein.

[0140] Figures 3, 4, and 5 illustrate flow charts showing portions of an exemplary process for taking a material measurement according to an embodiment of the subject matter described herein.

[0141 ] Figure 6 illustrates a flow chart for an exemplary process for taking a material measurement according to an embodiment of the subject matter described herein.

[0142] Figure 7 is block diagram illustrating a system for taking a material measurement, which may include assessing material anisotropy, according to an embodiment of the subject matter described herein.

[0143] Figure 8A shows the group velocity [ms ^'1 ] of a shear wave travelling in a transversely isotropic (Tl) material.

[0144] Figure 8B shows the XYZ coordinate system used to describe the material's axis of symmetry (AoS) orientation with respect to the acoustic radiation force impulse (ARFI) excitation 2D point spread function (PSF).

[0145] Figures 9A through 9D show the spatial distribution of the ARFI-2D

PSF in the lateral-elevational plane at the focal depth (36 mm) for focal configurations of F/1 .5, 3, 4, and 5, respectively with their respective degree of asymmetry values.

[0146] Figure 10A shows elliptical shear wave propagation in a

representative modeled Tl material. [0147] Figure 10B shows both expected (True) and measured (Measured) shear wave velocities (SWVs) along and across the AoS for the seven material-PSF orientations.

[0148] Figure 10C shows a plot of measured AoS orientations from the fitted ellipse versus true material AoS orientations.

[0149] Figures 1 1 A through 1 1 C show peak displacement (PD, left y-axis) and peak strains (PS, right y-axis) versus PSF-material orientation for ARFI focal configurations of F/1 .5, 3, and 5, respectively.

[0150] Figure 12A shows PD as a function of PSF-material orientation for the F/1 .5 2D ARFI PSF focal configuration in eight Tl materials.

[0151] Figure 12B shows PD as a function of material-PSF orientation for angles up to 360° for the F/1 .5 2D ARFI PSF focal configuration in three Tl materials with fitted sinusoidal function and root mean square error of the fit.

[0152] Figures 13A through 1 3D show ratios of PDs (PD ₉₀ /PD ₀ ) versus degree of material anisotropy (ratio of Young's moduli or ratio of shear moduli) for F/1 .5, 3, 4, and 5, respectively, 2D ARFI PSF focal configurations.

[0153] Figure 14A shows finite element method (FEM)-derived and ultrasonically-tracked PD in two representative Tl materials.

[0154] Figure 14B shows FEM-derived and ultrasonically-tracked ratio of PDs (PDgo/PDo) versus the ratio of shear moduli (μ_/μτ) for F/1 .5 2D ARFI

PSF focal configuration.

[0155] Figure 15A shows FEM-derived and ultrasonically-tracked PD in a representative Tl material (E _L =27.36 kPa, E _T = 27.09 kPa, μι = 36 kPa, and μ _τ = 9 kPa) as a function of PSF-material orientation for the F/1 .5 2D ARFI PSF focal configuration.

[0156] Figure 15B shows FEM-derived (dashed lines) and ultrasonically- tracked (solid lines) ratio of PDs (PD ₉₀ /PD ₀ ) in the representative material versus ARF duration for F/1 .5 2D ARFI PSF focal configuration.

[0157] Figure 16 is a graph of the degree of asymmetry of an ARF excitation as a function of depth.

[0158] Figure 17 is a graph of the degree of asymmetry of an ARF excitation as a function of depth for five different excitations, each having a different focal depth. Detailed Description of the Preferred Embodiments

[0159] The embodiments set forth below represent the necessary information to enable those skilled in the art to practice the invention and illustrate the best mode of practicing the invention. Upon reading the following description in light of the accompanying drawing figures, those skilled in the art will understand the concepts of the invention and will recognize applications of these concepts not particularly addressed herein. It should be understood that these concepts and applications fall within the scope of the disclosure and the accompanying claims.

[0160] Figure 1 A illustrates an ultrasound transducer according to an embodiment of the subject matter described herein. In the embodiment illustrated in Figure 1 A, a transducer 10 includes an ultrasonic emitter array 12, which produces an acoustic radiation force (ARF) having a direction of force 14. A coronal plane 1 6 is defined as the plane normal to the direction of force 14. In the embodiment illustrated in Figure 1 A, the ultrasonic emitter array 12 comprises a linear array of emitter elements, and thus may alternatively be referred to as "linear array 12".

[0161] Figure 1 B illustrates the ultrasound transducer in use over a material sample, showing an oval force profile produced within the material sample. In the embodiment illustrated in Figure 1 B, the transducer 1 0 produces a force having a direction 14. The force is directed into a material sample 18, which may be tissue. In one mode of operation, the transducer 10 produces a force having an oval-shaped profile. In another mode of operation, the transducer produces a force having a circular or other profile with equi-length axes.

[0162] In the embodiment illustrated in Figure 1 B, the transducer 10 is producing a force having an oval or other profile with a long axis 20 and a short axis 22. The shape delineated using dotted lines within the material samplel 8 is intended to illustrate the point that the acoustic force produced by the transducer 10 has a point of highest energy 24 - represented by a shaded oval - located at a focal depth 26 within the material 1 8. Above and below the focal depth 26, the size of the focal area expands and thus has a smaller energy per volume. [0163] In the example illustrated in Figure 1 B, the shape of the force profile at focal depth is an oval. The long axis of the oval force profile 24 has an angle theta (Θ) relative to an arbitrarily chosen reference vector 28 associated with the material sample 1 8. For a transducer with a linear array of elements, the long axis of force 20 is perpendicular to the long axis of the array. To change the angle theta (Θ), it is necessary to physically rotate the linear array 12 around an axis of rotation aligned with the direction of force 14. In another mode of operation (not shown in Figure 1 B) the shape of the force profile at focal depth is circular.

[0164] Figure 2 illustrates a flow chart for an exemplary process for assessing material anisotropy according to an embodiment of the subject matter described herein. In the embodiment illustrated in Figure 2, the method includes:

• applying a first force to a material sample and taking a first

measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane (step 100);

• applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force (step 1 02); and

· calculating a material anisotropy based on the first and second

measurements (step 104), wherein the first force and the second force are applied in any order.

[0165] It is noted that the first and second forces may, but need not always, occupy the same coronal plane. For example, a first oval force may have one orientation within a coronal plane and a second oval force may have a different orientation within the same coronal plane, e.g., at the same axial depth. In another example, a first oval force may occupy one coronal plane at a first axial depth and a second oval force may occupy another coronal plane at a second axial depth different from the first axial depth. Other

combinations of orientation, axial depth, focal points, or other characteristics, are within the subject matter of the present disclosure.

[0166] It is also noted that the forces may be applied at a wide range of depths along the axis of force. For example, one pair of forces having different profiles or orientations may be applied at one depth, another pair of forces may be applied at a different depth, and so on. Likewise a pair of forces may be applied at one set of respective depths, another pair of forces may be applied at another set of respective depths, and so on.

[0167] It is also noted that while most of the examples presented herein describe operations that involve a pair of forces, the concepts and techniques described herein may also be implemented using more than two forces and/or using more than two force profiles or orientations.

[0168] In one embodiment, multiple measurements may be made at a variety of different angles. For example, rather than looking at the ratio of values measured at 0 and 90 degrees to reflect anisotropy, the method may include looking at the values achieved at 0, 1 0, 20, 30, 40, ... 90 and larger degree angles and evaluate the measured value versus orientation angle. In alternative embodiments, more than two orientation angles could be measured and the method may include evaluation of more than just the ratio of values measured at these orientations.

[0169] Figures 3, 4, and 5 illustrate flow charts showing portions of an exemplary process for taking a material measurement according to an embodiment of the subject matter described herein. In the embodiment illustrated in Figure 3, the method includes:

• applying a first force to a material sample and taking a first

measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a first direction within the coronal plane (step 200);

• applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane (step 202); and

• normalizing the first measurement using the second measurement to produce a first normalized measurement, wherein the first force and the second force are applied in any order (step 204).

[0170] In the embodiment illustrated in Figure 4, the process described in Figure 3 continues with the following:

• applying a third force to a material sample and taking a third

measurement, the third force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force (step 300); and

• normalizing the third measurement using the second measurement to produce a second normalized measurement (step 302), wherein the first force, the second force, and the third force are applied in any order.

[0171] In the embodiment illustrated in Figure 5, the process described in Figure 3 continues with the following alternative steps to Figure 4:

· applying a third force to a material sample and taking a third

measurement, the third force having a direction, having a coronal plane normal to the direction of the force, and having an oval or other profile with long and short axes within the coronal plane, the long axis being oriented in a second direction within the coronal plane different from the first direction of the first force (step 400);

• applying a fourth force to the material sample and taking a fourth

measurement, the fourth force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane (step 402); and · normalizing the third measurement using the fourth measurement to produce a second normalized measurement (step 404), wherein the first force, the second force, the third force, and the fourth force are applied in any order. [0172] Figure 6 illustrates a flow chart for an exemplary process for taking a material measurement according to an embodiment of the subject matter described herein. In the embodiment illustrated in Figure 6, the method includes:

· applying a first force to a material sample and taking a first

measurement, the first force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane (step 500);

• applying a second force to the material sample and taking a second measurement, the second force having a direction, having a coronal plane normal to the direction of the force, and having a circular or other profile with equi-length axes within the coronal plane (step 502); and

• determining a material property based on a difference between the first measurement and the second measurement (step 504), wherein the first force and the second force are applied in any order.

[0173] Figure 7 is block diagram illustrating a system for taking a material measurement, which may include determining material anisotropy, according to an embodiment of the subject matter described herein. In the embodiment illustrated in Figure 7, the system includes a subsystem for force generation 30, a subsystem for motion dynamics observation 32, a subsystem for displacement calculation 34, a subsystem for parameterization 36, a subsystem for mathematical calculation 38, and a subsystem for producing system output 40.

[0174] The force generation subsystem 30 may produce ARF, a magnetic field, and/or a mechanical indentation. In the embodiment illustrated in Figure 2, the force generation subsystem 30 produces two forces directed towards a material sample, one force having a first orientation (orientation 1 ) in the coronal plane and the other force having a second orientation (orientation 2) in the coronal plane.

[0175] The motion dynamics observation subsystem 32 may make measurements based on ultrasound, magnetic resonance imagery, optical input (such as but not limited to pictures, videos, etc., including from highspeed cameras), optical coherence tomography, or using a mechanical means, such as but not limited to a micrometer. It should be noted that the angles of the forces and measurements with respect to the sample as shown in Figure 7 were chosen for illustration purposes only and are not intended to convey any particular orientation.

[0176] The displacement calculation subsystem 34 receives data produced by the motion dynamics observation subsystem 32 and calculates

displacement of the sample. In the embodiment illustrated in Figure 7, the displacement calculation subsystem 34 receives data for each measurement.

[0177] The parameterization subsystem 36 uses displacement information received from the displacement calculation subsystem 34 to determine certain parameters or material properties / material characteristics, such as PD, a value for TAU, time to peak (TTP), recovery time, Viscoelastic Response (VisR) relative elasticity, VisR relative velocity, etc.

[0178] The mathematical calculation subsystem 38 uses the parameters produced by the parameterization subsystem 36 to produce an output value, which is sent to system output subsystem 40. The calculations that may be performed by the mathematical calculation subsystem 38 include, but are not limited to, normalization of one data using another data; calculation of ratios of the two output values; and coefficient adjustment.

[0179] It will be understood that applying a force having a particular orientation or force profile may involve applying the force once or may involve applying the force a plurality of times. Applying a force a plurality of times may involve applying a force for a first duration of time, ceasing to apply the force for a second duration of time (e.g., pausing), and repeating the apply/pause cycle multiple times. The duration of time during which the force is applied, and/or the duration of the pause between one application of the force and a subsequent application of the force, does not need to be constant during the course of the repetitions. Likewise, the force that is applied need not be unchanged over the course of the repetitions.

[0180] In one embodiment, a plurality of applications of force and/or measurements is performed at different angles from each other. The angle difference between any two applications of force / measurements does not need to be 90 degrees. For example, a first application of force and measurement may be taken in one orientation and a second application of force and measurement may be taken at another orientation, where the difference in orientation may be 75 degrees, 50 degrees, or any other difference in orientation. Methods

[0181] Finite element simulations of a 1 .25 mm radius, Tl elastic, cylindrical inclusion embedded in an isotropic elastic material were performed using a commercially available Finite Element Method (FEM) solver, "LS- DYNA". A simulated 2D, 38 mm square matrix array transducer with 1922 elements, simulated with Field I I, was used to model the ARF excitation geometry. Each phantom was excited (for 70 us @ 4.21 MHz) by an asymmetrical ARF point spread function, with the long axis of the excitation oriented first along and then across the Tl material's axis of symmetry. The aperture was then swept laterally in 0.35 mm increments to acquire data along a 2D field of view. The degree of anisotropy was then computed as the ratio of the resultant peak displacements, and parametric images of degree of anisotropy and peak ARFI-induced displacement were rendered.

[0182] The imaging simulations were repeated for a total of 12 materials with differing material properties. The shear modulus of the isotropic material was held static at 1 .7 kPa, while the longitudinal and transverse shear moduli of the Tl material were varied from 5 to 1 0 and 5 to 30 kPa, in steps of 5 kPa, respectively. Contrast ratios were computed for each inclusion in the resulting parametric anisotropy images. Contrast-to-noise ratios (CNR) were also computed for the anisotropy and peak displacement parametric images.

[0183] The axial range over which anisotropy measurements could be performed was also investigated using the same 12 material combinations, in a homogeneous Tl material model. The axial range was determined by measuring the depth over which the ratio of peak displacements in the homogenous Tl material varied less than 5% for all 12 materials. The range was reported as mean ± standard deviation.

[0184] Using the same methods, a three-dimensional (3D) image of anisotropy was also simulated in a single material model with background shear modulus 1 .7 kPa and inclusion shear moduli of 5 and 30 kPa in the transverse and longitudinal directions, respectively. The 3D image was formed by sweeping the aperture in both the lateral and elevational directions and taking the ratio of the peak displacements at each position.

Results

[0185] For each of the six phantoms with ratios of longitudinal to transverse shear moduli ranging from 1 .0 to 6.0 in steps of 1 .0, the ratios of peak displacements increased proportionally with an increase in the ratio of shear moduli (1 .00 ± 0.08; 1 .06 ± 0.08; 1 .10 ± 0.07; 1 .1 3 ± 0.07; 1 .15 ± 0.07; 1 .17 ± 0.07). CNR measurements for peak displacement with the long axis of the force aligned with the axis of symmetry ranged from -25.7 dB to -16.9 dB from the least to most anisotropic lesion. With the long axis of the force normal to the axis of symmetry, the CNR ranged from -4.3 dB to +3.4 dB, from the least to most anisotropic lesion. Anisotropy images, created from the ratio of peak displacements yielded CNR values ranging from +18.7 dB to +20.9 dB in the same lesions. The axial extent of the anisotropy measurement field set was found to be 12.47 ± 2.35 mm for the materials used in this work. Using the 3D imaging technique, the degree of anisotropy measured throughout the 3D lesion was 1 .12 ± 0.04.

Conclusions

[0186] The results of the phantom simulations indicate that the measured degree of anisotropy is directly proportional to the modeled ratio of shear moduli. This work has also shown that the degree of anisotropy

measurement is valid over an axial range larger than the simulated inclusion, confirming the validity of phantom simulations. These data indicate that degree of anisotropy can be imaged using a ratio of peak displacements induced by a 2D matrix transducer.

[0187] To date, no other known ultrasound-based system is able to make images of anisotropy directly from the ratio of peak displacements. Other technologies that are able to make such images are dependent on the propagation of shear waves, which leads to decreased spatial and temporal resolution. Our method of creating anisotropy images depends only on the maximal displacement at each position through depth, which makes the resolution of our technique finer than that of shear-wave based methods. Our technique uses a two-dimensional matrix array transducer to create the appropriate forcing geometries in each orthogonal direction. Currently, two- dimensional matrix array transducers are expensive to build and have been demonstrated for producing radiation force in limited cases. It is noted that a 2D matrix is not necessarily needed to eliminate rotation. A row-column or other configuration of transducer could be used. In addition, as discussed below, it is possible to eliminate rotation by adjusting the focal depth.

In Silico Verification

[0188] Many soft tissues exhibit anisotropic mechanical properties. The simplest anisotropic material model is that of a Tl material, which is defined by a plane of isotropy perpendicular to an AoS. The Tl material model has been used to represent the strong anisotropy in mechanical properties observed in muscles, tendons, kidney, brain white matter, bones, and breast.

[0189] The mechanical properties of Tl materials may be described in terms of Young's and shear moduli; however, the relationship between these moduli is more complex in Tl than in isotropic materials. In linear, elastic, isotropic, incompressible (Poison's ratio v = 0.5) materials, the Young's modulus (£) is related to the shear modulus (μ) according to E = 3μ. The shear modulus may be estimated from the shear wave velocity { V) as , μ = pV ² , where p is the density of the material.

[0190] In contrast, for a Tl material, the Young's and the shear moduli differ along versus across the AoS. In the limit of small displacements, the constitutive equation for a linear elastic solid is defined by Hooke's law as aij ^{= c} ijkl ^£ kl (1)

[0191 ] where, σ _ί; is a stress tensor, s _kl is an infinitesimal strain tensor, and c _ijkl is a fourth-order stiffness tensor representing the properties of a material. Due to symmetry in Tl materials, c _ijkl contains five independent elastic constants (instead of 81 ). The reduced c _ijkl with five independent elastic constants will be denoted as c. The compliance matrix s = c ¹ can be expressed as follows

where

MT = (3)

2 (1 + v _TT )

[0192] ^ and E _T are the longitudinal and the transverse Young's moduli, jtt _L and jtt _T are the longitudinal and transverse shear moduli, and v _TT and v _LT are Poisson's ratio in the plane of isotropy and in the plane of symmetry, respectively. Longitudinal (L) and transverse (T) directions are defined relative to the AoS. The fractional volume change of an infinitesimal volume subjected to stress (i.e., dilatation) must be zero for an incompressible material. So, Poisson's ratio for an incompressible Tl material is given as.

E _T

[0193] Note that, incompressible Tl materials are defined by three material constants: μ _τ , μ _ι , and . Rouze et al. showed an elliptical propagation of shear waves induced by ARFI excitation in incompressible Tl materials. More specifically, the components of the shear wave group velocities form an ellipse, and they can be found from the following relationship in the Cartesian coordinate system

( _x cos Θ + V _y sin #) ( _x sin Θ— V _y cos #)

(5)

[0194] where, (V _x , V _y ) defines the path of a point on the ellipse, V _L and V _T are major and minor axis of the ellipse, and Θ is the angle between the major axis and the X-axis. Q, V _L = ^ , ard V _T = ^ also define the angle between the AoS and X-axis, SWV along the AoS, and SWV across the AoS, respectively. The orientation of the ellipse gives the orientation of the material's AoS.

[0195] Figure 8A shows the group velocity [ms ^'1 ] of a shear wave travelling in a Tl material with E _T = 34.56 kPa, E _L =216 kPa, μ _τ = 9 kPa, μι = 25 kPa and p = 1 gem ^"3 . The major axis length, minor axis length, and orientation of the ellipse give the SWV along the AoS (5 ms ^"1 ), across the AoS (3 ms ^"1 ), and orientation of the material's AoS (60°), respectively.

[0196] The anisotropy of Tl materials has been assessed as the ratio of SWV measured along versus across the AoS, or the "fractional anisotropy". However, SWV measurement is subject to error from tissue inhomogeneities within the measurement region, and the requirement for observations of displacement over time at multiple lateral locations reduces frame rate unless a plane wave imaging method is implemented. The drawback of plane wave imaging is a lack of transmit-focusing, which leads to loss of SNR and resolution and challenges applications in deep tissues and obese patients. To increase the SNR, coherent spatial compounding may be employed.

However, coherent compounding reduces the frame rate relative to conventional plane wave imaging. Assessing the degree of anisotropy in a Tl material in a single, focal spatial location could help to reduce error from inhomogeneities and increase frame rates without sacrificing lateral resolution or SNR.

[0197] The subject matter disclosed herein includes ARFI imaging, in which induced displacements are observed in only the ARFI region of excitation (ROE), as a novel approach to describing the degree of anisotropy in Tl materials. The rationale for this approach is that the mechanical response of Tl materials depends on the direction of loading with respect to the AoS. For example, it has been shown that displacements were higher when the long axis of an asymmetrical indenter was aligned parallel versus perpendicular to the direction of the AoS in Tl elastic materials. The inventors have determined that ARFI-induced PDs will be larger when the long axis of an asymmetrical ARFI2D PSF is aligned along versus across the AoS. The inventors have determined that the ratio of PDs achieved with the long axis of an asymmetrical ARFI- 2D PSF aligned along versus across the AoS reflects the degree of anisotropy in Tl materials.

[0198] The concept was first proven in silico, including an expanded description of the employed in silico model, validation of the Tl nature of the modeled materials, new results and analyses for 52 different simulated Tl materials, implementation of simulated ultrasonic displacement tracking, variation in the ARF duration, and a comprehensive discussion of the significance of the reported results. Table 1 , below, lists the used symbols.

Table 1

E _L Longitudinal Young's modulus

E _T Transverse Young' s modulus

L Longitudinal shear modulus

μ _τ Transverse shear modulus

V-L Ratio of the Longitudinal and Transverse shear moduli

μ _τ

EL Ratio of the Longitudinal and Transverse Young's moduli

PD _ang i _e Peak displacement when the AoS is oriented at angle (0°, 15°,

30°, 45°, 60°, 75°, and 90°)

PD ₉₀ Ratio of peak displacement (PD) when the AoS is oriented

PD ₀ along versus across the long axis of the ARFI-2D PSF

Methods

[0199] Figure 8B shows the XYZ coordinate system used to describe the material's AoS orientation with respect to the ARFI excitation 2D PSF, where lateral axis = Y axis, elevational axis = X axis, and axial axis = Z-axis. The Tl material's AoS is oriented along the X-axis, with the plane of isotropy in the YZ-plane and the plane of symmetry in the XZ-plane. Mechanical properties are similar in the plane of isotropy. An asymmetrical ARFI 2D PSF in the lateral-elevational plane, produced at the focal depth by a linear array transducer, is represented by a gray ellipse. For this material-PSF

arrangement, the long axis of the PSF is aligned with the AoS, for a "90°" orientation [0200] A FEM model was used to simulate the dynamic responses of an elastic Tl solid to an ARFI excitation. The 3D, dynamic response of the elastic Tl solid was solved through the balance of linear momentum using the software application "LS-DYNA3D" (Livermore Software Technology corp. Livermore, CA) with an explicit time-domain method. The axial, lateral, and elevational extents of the Finite Element (FE) mesh were each 25 mm for the purpose of measuring SWV to validate the material model; however, these values were adjusted to 25 mm, 12 mm, and 12 mm, respectively, for greater computational efficiency while estimating anisotropy based on the ARFI PD at the focal depth. The axial, lateral, and elevational dimensions spanned 25 mm to 50 mm, -6 mm to 6 mm, and -6 mm to 6 mm, respectively. The FE mesh element size was 0.2 χ 0.2 χ 0.2 mm ³ _.

[0201] The center of the transducer was aligned to the origin of the XYZ coordinate system. A built-in material model

(MAT_ORTHOTROPICTROPIC_ELASTIC in LS-DYNA3D) was used to simulate 52 different homogeneous, linear, elastic Tl materials with the following ranges of elastic constants: E _T = 6.53 - 89.69 kPa; E _L = 1 5.4 - 486.34 kPa; μ _τ = 1 .67 - 29.7 kPa; μι = ~\ 0 - 90 kPa. In the Tl elastic material, compressibility is controlled by the value of one Poisson's ratio, V _L T that must equal 0.5 to retain perfect incompressibility according to equation (4).

However, a perfectly incompressible material cannot be modeled in DYNA. To overcome this limitation, we set the simulated value of VL-r to 0.499, which approaches incompressibility conditions. Then, the second Poisson's ratio (VTT) for the simulated material was calculated using equation (4). Note that a Poisson's ratio of 0.499 has previously been implemented in DYNA

simulations for isotropic materials and for Tl materials. The density of all Tl materials was fixed to p = 1 gem ^"3 .

[0202] Table 2 contains the elastic constants with ratio of shear and Young's moduli of all simulated materials. Note that these elastic constants were generally selected to reflect known muscle properties (μι_ = 16-55 kPa; μτ = 1 -9 kPa>). For example, values of the shear moduli were selected such that the SWV was faster along versus across the muscle fibers (i.e., the AoS), and Young's moduli were selected such that the materials were stiffer along than across the AoS. Further, to explore the utility of our technique in more extreme scenarios, we also created one highly and one minimally anisotropic Tl material. A seven-element thick perfect matching layer (PML) was implemented around the mesh using the built-in material model

( M AT P M L_0 RTHOTROPIC_ELASTIC in LS-DYNA3D) with matching to the elastic constants used in the MAT_ORTHOTROPICTROPIC_ELASTIC material model. The purpose of the PML layer was to simulate an infinite medium and remove spurious wave reflections from the boundaries of the mesh.

[0203] Using a linear acoustics modeling package, "Field II", the ARFI field was simulated by first computing the acoustic intensity at each nodal location based on the method described in. The Siemens VF7-3 linear array transducer parameters were used to calculate the acoustic intensity in Field II with a transmitting frequency of 4.21 MHz, a fixed elevational focus at 37.5 mm, and lateral focal configurations of F/1 .5, 3, 4, and 5. The axial focal depth was fixed to 36 mm for all simulations. Field ll-derived intensity values were scaled to a peak intensity of 10,000 W/cm ² , and the volumetric acoustic radiation force magnitude was calculated using the expression:

where a is the absorption coefficient of the medium (assumed to be 0.5 dB/cm/MHz), / is the temporal-average beam intensity over a volume, and c is the speed of sound through the medium (assumed to be 1540 m/s).

[0204] The point load forces were calculated by spatially sampling the volumetric force. Point forces were then superimposed on the FE mesh, and the force magnitude was kept constant for 70 με for all simulated materials. Additionally, for one representative material, the force magnitude was kept constant for 25 με, 50 με, 70 με, and 100 με. The simulation was allowed to run for 3.5 ms total, and data were sampled every 0.1 ms for all materials, approximating a conventional ARFI ensemble with 10 kHz pulse repetition frequency (PRF). Additionally, for one representative material, data were sampled every 25 με, approximating a 40 kHz PRF.

[0205] The dimension of the generated displacements matrix with each simulation is four-dimensional (4D) (axial ^χ lateral ^χ elevational ^χ time). PD over time at the center of the push was calculated. A 4D matrix of strain was calculated by taking the spatial gradient of displacements at each time instant. Similar to PD, PS was calculated over time at the center of the push. To study the effect of PSF-AoS orientation on PD and PS, the AoS of the Tl materials was oriented at seven different angles: 0°, 15°, 30°, 45°, 60°, 75°, and 90°, where angle 0° and 90° correspond to the long axis of the ARFI-2D PSF being aligned across and along the material's AoS, respectively.

Additionally, three representative materials were orientated from 0° to 360° at an interval of 15°.

[0206] To evaluate the impact of ultrasonic tracking on PDs in Tl materials, simulation of ultrasonic tracking of ARF-induced displacement was adapted from the methods developed by Palmeri et al. First, a three-dimensional Field II scatterer phantom with fully-developed speckle (31 scatterers per resolution cell) was defined to span the volume of the FEM mesh. Next, the FEM- derived displacements were used to linearly-interpolate scatterer positions for every time step in the ARFI ensemble using MATLAB (Mathworks Inc., Natick, MA). After generating the scatterer position matrices for each time step in the ARFI ensemble, the corresponding radio frequency (RF) lines were simulated using Field I I. White Gaussian noise was added to each RF line using the awgn function in Matlab to simulate system signal to noise ratio (SNR) of 40dB. The Siemens VF7-3 linear array transducer parameters were used to simulate tracking pulses with a transmitting frequency of 6.15 MHz and lateral focal configurations of F/0.75. Tracking pulses were focused at 36 mm axially on transmit, and dynamic focusing was used on receive. Motion tracking was performed using one-dimensional normalized cross correlation with a kernel length of 2λ. After motion tracking, PD was calculated. For each simulated Tl material in a representative subset of six, ten unique phantoms with independent speckle realizations were generated.

[0207] To confirm that the modeled material was behaving as a Tl material, ARFI-induced SWV along and across the AoS was measured for seven different material-PSF orientations, i.e. the long axis of the PSF was aligned at 0°, 1 5°, 30°, 45°, 60 °, 75°, and 90° relative to the material's AoS. For each material-PSF orientation, the shear wave group velocities at angles from 0° - 360° in the lateral-elevational plane were measured at the focal depth using the method described by Wang et al. Then, an ellipse was fit through the measured SWVs to determine the SWV along the AoS, the SWV across the AoS, and the orientation of the material's AoS. The velocities were related to Tl material property according to equation (5).

D. Symmetry of ARFI Excitations

[0208] The spatial distribution of the ARFI 2D PSF was varied by changing the focal configuration, or F/#, defined as: /# = , where z is the axial focal depth and D is the width of the active aperture. For a standard 1 -D linear array, the lateral F/# is adjusted by electronically selecting the number of active elements, while the elevational F/# = is set by the height (h) of the transducer elements and the fixed-focused lens. The lens of the Siemens VF7-3 linear array transducer has a fixed elevational focus at 37.5 mm. We used "xdc_focused_multirow" function in Field II to model the fixed elevational focus at 37.5 mm. To manipulate the degree of asymmetry of ARFI excitations at the desired focal depth, lateral focal configuration was adjusted, and degree of asymmetry was measured as the ratio of the lateral versus elevational extent of the ARFI 2D PSF. For example, for an F/1 .5 lateral focal configuration at a focal depth of 36 mm, the lateral and elevational extents of the ARFI PSF were 1 mm and 3 mm, respectively, yielding a degree of asymmetry of 0.33. The degree of asymmetry of the ARFI-2D PSF at 36 mm for focal configurations of F/3, 4 and 5 was 0.69, 0.85, and 1 .15, respectively.

[0209] Figures 9A through 9D show the spatial distribution of the ARFI-2D PSF in the lateral-elevational plane at the focal depth (36 mm) for focal configurations of F/1 .5 (Figure 9A), F/3 (Figure 9B), F/4 (Figure 9C), and F/5 (Figure 9D) with their respective degree of asymmetry values, degree of asymmetry was measured as the ratio of the lateral versus elevational extent of the PSF. The intensity bar indicates normalized nodal force at the focal depth, with normalization relative to the maximum nodal force. The F/1 .5 lateral focal configuration yielded the most asymmetric ARFI PSF, while a nearly symmetrical ARFI PSF was achieved for the F/5 lateral focal configuration.

Results

[0210] Figure 10A shows elliptical shear wave propagation in a

representative modeled Tl material with E _T = 34.56 kPa, E _L = 216 kPa, μτ = 9 kPa, μι_ = 25 kPa, and p = 1 gem ^"3 . The SWVs along and across the AoS are expected (from equation (5)) to be 5 ms ^"1 and 3 ms ^'1 _, respectively. The black dots represent the measured SWV at angles from 0° - 360° and the fitted ellipse through these dots is shown as the black line.

[0211] Figure 10B shows both expected (True) and measured (Measured) SWVs along and across the AoS for the seven material-PSF orientations. The true SWV was calculated from the analytical equation (equation (5)), whereas the measured SWV was calculated from the fitted ellipse. The maximum percentage of error in the measured SWV along and across the AoS was < 2.5 %.

[0212] Figure 10C shows a plot of measured AoS orientations from the fitted ellipse versus true material AoS orientations, indicating the Pearson correlation coefficient was r = 1 .

[0213] Figures 1 1 A through 1 1 C show PD (left y-axis) and PS (right y-axis) versus PSF-material orientation for ARFI focal configurations of F/1 .5, 3, and 5. Figures 1 1 A and 1 1 B show results for two different simulated Tl materials, with the material in Figure 1 1 A having a higher degree of anisotropy. Metrics for material anisotropy were the ratio of the longitudinal and the transverse shear moduli (μι_/μτ) and the ratio of the longitudinal and transverse Young's moduli (E _L / E _T ). The transverse and longitudinal Young's moduli and the transverse shear modulus were kept fixed for both materials (E _T = 30.24 kPa, E _L = 47.25 kPa, and μτ = 9 kPa) but the longitudinal shear modulus was varied (μι_ = 36 kPa in Figure 1 1 A and μι_ = 16 kPa in Figure 1 1 B). Note that the PD and PS are higher in Figure 1 1 B than in Figure 1 1 A for the same focal configuration and orientation of the AoS because i _L is lower in Figure 1 1 B than in Figure 1 1 A. For both materials, PD and PS spanned minimum to maximum values when the AoS was oriented from 0° to 90° for F/1 .5 and F/3 ARFI focal configurations. However, PD and PS were constant across material orientations for the F/5 ARFI focal configuration.

[0214] Figure 1 1 C illustrates cross-sections of the 2D ARFI excitation PSFs in the lateral-elevational plane at the focal depth with respect to the alignment of the material's AoS. For asymmetrical ARFI excitations (F/1 .3 and F/3), the long axes of the ARFI PSFs were observed to rotate with respect to the AoS for different orientations. However, the interrogating ARFI region was the same in all orientations for the symmetrical focal configuration (F/5).

[0215] Our data presentation from here forward will be limited to PD because PD and PS were observed to exhibit similar trends. To quantify the variability in PD with material-PSF orientation, we measured the ratio of PDs when the material's AoS was at 90° versus 0° (i.e., PD ₉₀ / PD ₀ ). The ratio of PDs was 1 .56 and 1 .24 in the materials of panels (a) and (b), respectively when the ARFI focal configuration was F/1 .5. The ratios of shear moduli were 4.0 and 1 .78 in the materials of the panel (a) and (b), respectively. The ratio of Young's moduli was 1 .56 in both materials

[0216] Figure 12A shows PD as a function of PSF-material orientation for the F/1 .5 2D ARFI PSF focal configuration in eight Tl materials. For the material with μι_ / μτ and E _L / E _T closest to 1 , PD did not depend on the material orientation (dotted square line in Figure 12A). On the contrary, variability in PD with PSF-AoS orientation was greatest for the material with the highest degree of anisotropy (solid ' ^* ' line in Figure 12A).

[0217] Figure 12B shows PD as a function of material-PSF orientation for angles up to 360° for the F/1 .5 2D ARFI PSF focal configuration in three Tl materials with fitted sinusoidal function (Peak Displacement = B ₀ +

B ₁ sin B ₂ * Orientation + B ₃ )) and root mean square error (RMSE) of the fit. E _L, E _t , and μτwere fixed to 47.16 kPa, 30.231 kPa, and 9 kPa, respectively, whereas μι_ was varied to 16 (triangle), 36 (circle), and 90 (square) kPa in these three materials. The parameters of the fitted sinusoids were B ₀ = 14.34 (triangle), 1 1 .81 (circle), and 8.0 (square); B- = 1 .37 (triangle), 2.40(circle), and 2.58 (square); B ₂ = -2.0 (triangle), 2.02 (circle), and 2.04 (square); B ₃ = 4.74 (triangle), -1 .65(circle), and -1 .69 (square). Figure 12B similarly shows PD as a function of PSF-AoS orientation for angles up to 360° in three representative materials with increasing degrees of anisotropy. The parameters of sinusoidal functions fit to each material are listed in the caption. RMSE in the sinusoidal fit (see legend) shows that RMSE increases with increasing degree of shear modulus anisotropy. For the least anisotropic material (μι_ / μτ = 1 -78 and E _L / E _T = 1 .56), PD follows an approximately sinusoidal variation pattern as the PSF-AoS orientation angle increases. However, as the material increases in the degree of shear anisotropy, PD variation deviates from a sinusoidal profile when the long axis of the 2D ARFI PSF is aligned along the AoS (90° and 270° orientations). This deviation arises from large changes in PD associated with small changes in orientation angle when the asymmetric 2D ARFI PSF is primarily pushing along the AoS.

[0218] Rather than interrogating seven or more different material-PSF orientations, a more practical imaging scenario would entail measuring PD in only two orientations, i.e. with the long axis of the ARFI 2D PSF aligned across (0°) and along (90°) the material's AoS. Then, the ratio of these two PDs can be related to the degree of material anisotropy.

[0219] Figures 13A through 1 3D show ratios of PDs (PD ₉₀ /PD ₀ ) versus degree of material anisotropy (ratio of Young's moduli or ratio of shear moduli) for F/1 .5 and F/3 (Figures 1 3A and 12B) and for F/4 and F/5 (Figure 13C and 13D) 2D ARFI PSF focal configurations. Figures 13A and 13C show the impact of anisotropy in Young's modulus on the ratio of PDs when degree of anisotropy in shear modulus is fixed (μι_/μτ = 1 .01 , 1 .40, 1 .78, and 6.0).

Figures 13B and 13D show the impact of anisotropy in shear modulus on the ratio of PDs when degree of anisotropy in Young's modulus is fixed (E _L /E _T = 1 .01 , 1 .56, 1 .78, and 6.0). Note that Figures 13A and 13B use a different Y- axis scale from that used in Figures 13C and 13D.

[0220] First, note that for all four 2D ARFI PSF focal configurations, variation in the degree of Young's modulus anisotropy has little impact on PD ₉₀ / PD ₀ . However, PD ₉₀ / PD ₀ increases with increasing degree of anisotropy in shear modulus. Second, for all examined materials, PD ₉₀ / PD ₀ is larger when more asymmetric 2D ARFI PSF focal configurations are used.

[0221 ] Figure 14A shows FEM-derived (dashed lines) and ultrasonically- tracked (solid lines) PD in two representative Tl materials as a function of PSF-material orientation for the F/1 .5 2D ARFI PSF focal configuration. EL, E _T , and μτ were fixed to 27.36 kPa, 27.09 kPa, and 9 kPa, respectively, whereas μι_ was 49 kPa (triangle, gray) or 1 6 kPa (diamond, black).

Ultrasonically-tracked PDs are shown as mean ± standard deviation over 10 independent speckle realizations.

[0222] Figure 14B shows FEM-derived (dashed lines) and ultrasonically- tracked (solid lines) ratio of PDs (PD ₉₀ /PD ₀ ) versus the ratio of shear moduli (μι_/μτ) for F/1 .5 2D ARFI PSF focal configuration. Ultrasonically-tracked (gray, solid line, triangle) ratio of PDs is shown as mean ± standard deviation over 1 0 independent speckle realizations. E _L , E _T , and μτ were fixed to 27.36 kPa, 27.09 kPa, and 9 kPa, respectively, whereas μι_ was varied to values of 1 1 , 16, 25, 36, 49, and 90 kPa.

[0223] Two observations are notable. First, for both FEM-derived and ultrasonically-tracked cases, PDs follow a similar trend of spanning minimum to maximum values when the AoS is oriented from 0° to 90°. Second, in both materials, ultrasonically-tracked PDs are underestimated, consistent with expectation. However, this displacement underestimation does not substantially influence the ratio of PDs (PD ₉₀ /PD ₀ ), as shown in Figure 14B, for six simulated Tl materials with different ratios of shear moduli (μι_ / μτ )■ The FEM-derived and ultrasonically-tracked ratio of PDs increase similarly with increasing degree of anisotropy in shear modulus.

[0224] Figure 15A shows FEM-derived (dashed lines) and ultrasonically- tracked (solid lines) PD in a representative Tl material (E _L =27.36 kPa, E _T = 27.09 kPa, μι = 36 kPa, and μ _τ = 9 kPa) as a function of PSF-material orientation for the F/1 .5 2D ARFI PSF focal configuration with ARF durations of 25 με (diamond, black) or 1 00 με (triangle, gray). Ultrasonically-tracked PDs are shown as mean ± standard deviation over 10 different speckle realizations.

[0225] Figure 15B shows FEM-derived (dashed lines) and ultrasonically- tracked (solid lines) ratio of PDs (PD ₉₀ /PD ₀ ) in the representative material versus ARF duration for F/1 .5 2D ARFI PSF focal configuration. FEM-derived ratio of PDs is shown for two PRFs, 40 KHz (time sampling interval = 25 με, diamond) and 10 KHz (time sampling interval = 100 με, square).

Ultrasonically-tracked (gray, solid line, triangle) ratio of PDs is shown as mean ± standard deviation over 10 different speckle realizations. Tracked (gray, solid line, triangle) ratio of PDs is shown as mean ± standard deviation over 10 different speckle realizations.

[0226] As expected, the longer ARF duration induces larger PDs. But closer examination reveals that, for both FEM-derived and ultrasonically- tracked cases, the ratio of PDs (PD ₉₀ /PD ₀ ) is larger for the longer duration force. As shown in Figure 1 5B, PD ₉₀ /PD ₀ is lower for the 25 με and 50 με ARF durations than for the 75 με and 100 με ARF durations when the PRF is 10 KHz. However, when the PRF is 40 KHz, PD ₉₀ /PD ₀ does not vary appreciably with ARF duration.

Discussion

[0227] A novel method that exploits ARFI-induced PD to interrogate the degree of anisotropy in Tl materials has been developed and evaluated in silico in this manuscript. In this approach, the degree of anisotropy is assessed as the ratio of PD achieved when the long axis of an asymmetric ARFI 2D PSF is aligned along versus across the material's AoS. This new method differs from previous techniques for assessing material anisotropy in that the new approach does not observe shear wave propagation. Therefore, local inhomogeneity and shear wave reflections may not be as confounding to estimates of the degree of anisotropy.

[0228] A key finding of this work is that, for asymmetric ARFI 2D PSFs (F/1 .5 and F/3.0), the induced PDs and PSs are maximized when the long axis of the PSF is oriented along (90°) and minimized when the long axis is oriented across (0°) the material's AoS, as demonstrated in Figure 10C.

However, for the symmetrical ARFI-2D PSF (F/5.0), induced PD and PS are independent of PSF-AoS orientation because the interaction of the ARFI 2D PSF with the material is consistent (Figure 1 0C). This result is meaningful because it suggests that material anisotropy may be selectively interrogated or obviated by implementing asymmetric or symmetric ARFI 2D PSFs, respectively. Interrogating degree of anisotropy is clinically relevant. For example, muscle fiber fragmentation and disordered fatty/fibrous deposition with dystrophic degeneration degrade the SWV derived degree of anisotropy in affected muscles. However, if the material-PSF orientation is not precisely controlled, anisotropy may be confounding in the context of longitudinal and/or cross-sectional study designs. In such scenarios, differences in measurement orientation would confuse true difference in tissue structure or composition. Thus, a symmetric ARFI 2D PSF would be useful to ensure consistent measurement of material property.

[0229] Another critical outcome of this work is the demonstration that, for an asymmetrical ARFI 2D PSF, the degree of difference in PDs achieved at 90° and 0° PSF-AoS orientations increases as the degree of anisotropy of the material increases, with consistent PD at 1 80° increments, as shown in Figures 12A and 12B. This result is significant because it suggests that evaluating the difference in achieved PDs, as in the ratio of the PDs for example, reflects true degrees of difference in the directionally-dependent material properties.

[0230] It is important to consider that the ratio of PDs is more strongly influenced by anisotropy in shear moduli rather than Young's moduli, as shown in Figures 1 3A through 13D. For focal configurations of F/1 .5 and F/3, percent change in the ratio of PD with respect to the percent change in ratio of Young's moduli was 0.83% and 0.08%, respectively, when the ratio of shear moduli was fixed to 1 .01 . However, for the same focal configurations, percent change in the ratio of PDs with respect to the percent change in ratio of shear moduli (-¾j - ) was 10.5% and 3.72%, respectively, when the ratio of

Young's moduli was fixed to 1 .01 . For F/4 and F/5 focal configurations, the ratio of PDs did not depend on the ratio of Young's moduli (Figure 13C), but they did increase slightly with increasing ratio of shear moduli when the ratio of Young's moduli was fixed. Percentage change in the ratio of PDs with respect to the change in the ratio of shear moduli (-^- ) was 1 .35% and

0.31 % for F/4 and F/5 focal configurations, respectively when ratio of Young's moduli was fixed to 1 .01 . These results suggest that our approach is more sensitive to degree of anisotropy in shear modulus than Young's modulus.

[0231] Although the data presented herein show that PD ₉₀ / PD ₀ reflects the ratio of shear and, to a lesser extent, Young's moduli, a quantitative relationship between these terms has not been developed. Rather, the described approach is relevant to distinguishing materials in terms of their relative degrees of anisotropy.

[0232] The analyses presented up to Figures 13A and 13B were applied directly to the FEM simulated material deformations, without accounting for ultrasonic displacement tracking. It is well established that ultrasonically tracked tissue deformation is subject to jitter and displacement underestimation errors. Indeed, the ultrasonically-tracked displacements have variance (jitter) and are underestimated. However, despite these errors, the ultrasonically-tracked PDs increase with increasing material orientation angle in a manner similar to the FEM-derived PDs. Moreover, the ratios of PDs (PDgo/PDo) in six different Tl materials do not substantially differ between FEM-derived and ultrasonically-tracked measures. These data suggest that, even when ultrasonically-tracked, ARFI-induced PDs reflect the degree of anisotropy in Tl materials.

[0233] One area requiring further investigation is the potential impact of orientation-dependent differences in displacement magnitude that could disproportionately influence the magnitude of jitter error in PD ₉₀ versus PD ₀ measures. Jitter magnitude is inversely related to the correlation coefficient, which decreases with increasing displacement. Therefore, jitter could be larger when the long-axis of the asymmetrical ARFI 2D PSF is aligned along versus across the AoS, which could corrupt PD ratios. Jitter magnitude is also inversely related to center frequency, fractional bandwidth, SNR, and kernel length. Implementing methods to control for jitter error may be necessary when extending this work to experimental realization. For example, averaging displacement profiles across space and time, maximizing SNR, and minimizing displacement magnitude will serve to drive down jitter to support experimental assessment of the degree of anisotropy.

[0234] Another topic for deeper study is the possible impact of orientation- dependent differences in shear wave speed that could disproportionately influence the degree of displacement underestimation in PD ₉₀ versus PD ₀ measures. Displacement underestimation will be less severe when the long axis of the ARFI 2D PSF is oriented along versus across the material's AoS due to faster shear wave propagation and, therefore, resolution of shearing artifacts under the tracking PSF. Such displacement underestimation could be minimized by tracking with tighter focal configurations than the ARFI excitation, which would be achievable using matrix array transducers capable of elevational focusing.

[0235] In addition to ultrasonically tracking ARFI-induced displacements, another confounding factor is misalignment. During in vivo imaging, it may be difficult to align the long axis of the asymmetrical ARFI 2D PSF along and across the AoS, and misalignment will introduce error in the assessment of anisotropy degree. In Figure 1 2A, for the most anisotropic material (solid line Ό'), a 1 5° error in alignment (75° instead of 90 °) introduces roughly 1 0% error in the perceived peak displacement at 90° and an associated 9% error in PD ₉₀ /PD ₀ . A potential solution is to rotate the imaging transducer through multiple orientation angles to identify the orientations that achieve the largest (90°) and smallest (0°) displacements.

[0236] Figures 1 5A and 1 5B point out another relevant consideration when relating ARFI-induced PDs to the degree of anisotropy in Tl materials. While the duration of the ARF excitation does not directly impact the ratios of PDs (PD ₉₀ /PD ₀ ) achieved in linearly elastic materials, measurements of PD ₉₀ /PD ₀ are effected by the temporal relationship of the ARF duration to the PRF. To understand this, consider that for the shorter (25 με and 50 με) ARF durations with a 1 0 KHz PRF, more time elapses for elastic recovery between ARF cessation and temporal sampling of induced-displacement than for the longer ( 75 [is and 1 00 με) ARF durations. Thus, peak displacement is more severely underestimated when using the shorter versus the longer ARF durations due to greater temporal sampling error. When the simulated PRF is increased to 40 KHz, temporal sampling error is reduced, and PD ₉₀ /PD ₀ measures are more consistent over ARF durations. These results show the importance of the temporal sampling period, and its optimization is a topic of ongoing investigation. Finally, the simulated materials that were evaluated in this study were purely elastic. The impact of viscosity on the estimation of anisotropy as a ratio of ARFI-derived PDs is yet to be determined.

[0237] In this in silico study, a novel method for assessing the degree of anisotropy in elastic Tl materials using ARFI-induced PDs is presented. The results demonstrate that for asymmetrical ARFI 2D PSFs, that the ratio of PDs (and PSs) achieved with the long axis of the PSF aligned along versus across the material's AoS (PD ₉₀ / PD ₀ ) qualitatively reflects the degree of anisotropy in the shear and, to a lesser extent, Young's moduli of the material. This may be useful for diagnostic purposes, particularly in tissues and organs such as muscle in which degree of anisotropy degrades with disease progression. Further, this work shows that it is possible to selectively obviate material anisotropy by using symmetrical ARFI excitations. Thus, symmetrical ARFI focal configurations may be desirable when the orientation of the ARFI excitation to the material's AoS is not specifically known and standardization of measurement is important, such as for longitudinal or cross-sectional studies of anisotropic tissues and organs. The feasibility of extension to experimental realization is demonstrated by simulating ultrasonic tracking of ARFI-induced displacements and by varying the ARF duration. Overall, this study demonstrates that ARFI-derived PD is relevant to interrogating anisotropy in elastic Tl materials.

Table 2

Material E _L /E _T μ_/μτ E _L E _T ML MT _π

Number

1 1.01 1.22 27.36 27.09 11.0 9.0 0.505

2 1.01 1.78 27.36 27.09 16.0 9.0 0.505

3 1.01 2.78 27.36 27.09 25.0 9.0 0.505

4 1.01 4.0 27.36 27.09 36.0 9.0 0.505

5 1.01 5.44 27.36 27.09 49.0 9.0 0.505

6 1.01 10.0 27.36 27.09 90.0 9.0 0.505

7 1.56 1.22 47.16 30.23 11.0 9.0 0.679

8 1.56 1.78 47.16 30.23 16.0 9.0 0.679

9 1.56 2.78 47.16 30.23 25.0 9.0 0.679

10 1.56 4.0 47.16 30.23 36.0 9.0 0.679

11 1.56 5.44 47.16 30.23 49.0 9.0 0.679

12 1.56 10.0 47.16 30.23 90.0 9.0 0.679

13 1.78 1.22 55.08 30.94 11.0 9.0 0.719

14 1.78 1.78 55.08 30.94 16.0 9.0 0.719

15 1.78 2.78 55.08 30.94 25.0 9.0 0.719

16 1.78 4.0 55.08 30.94 36.0 9.0 0.719

17 1.78 5.44 55.08 30.94 49.0 9.0 0.719

18 1.78 10.0 55.08 30.94 90.0 9.0 0.719

19 6.0 1.22 207.0 34.50 11.0 9.0 0.917

20 6.0 1.78 207.0 34.50 16.0 9.0 0.917

21 6.0 2.78 207.0 34.50 25.0 9.0 0.917

22 6.0 4.0 207.0 34.50 36.0 9.0 0.917 6.0 5.44 207.0 34.50 49.0 9.0 0.917

6.0 10.0 207.0 34.50 90.0 9.0 0.917

12.53 1.78 275.95 22.02 10.0 5.62 0.960

6.25 1.78 188.76 30.20 14.0 7.87 0.920

4.17 1.78 158.56 38.02 18.0 10.11 0.880

2.78 1.78 125.08 44.99 22.0 12.36 0.820

1.56 1.78 76.54 49.06 26.0 14.61 0.679

1.02 1.78 51.91 50.89 30.0 16.85 0.510

12.53 1.40 350.86 28.01 10.0 7.14 0.960

6.25 1.40 240.00 38.40 14.0 10.0 0.920

4.17 1.40 201.60 48.35 18.0 12.86 0.880

2.78 1.40 159.03 57.21 22.0 15.71 0.820

1.56 1.40 97.31 62.38 26.0 18.57 0.679

1.02 1.40 66.00 64.71 30.0 21.43 0.510

12.53 1.01 486.34 38.81 10.0 9.90 0.960

6.25 1.01 332.67 53.23 14.0 13.86 0.920

4.17 1.01 279.45 67.01 18.0 17.82 0.880

2.78 1.01 220.44 79.29 22.0 21.78 0.820

1.56 1.01 134.89 86.47 26.0 25.74 0.679

1.02 1.01 91.49 89.69 30.0 29.70 0.510

12.53 6.0 81.87 6.53 10.0 1.67 0.960

6.25 6.0 56.00 8.96 14.0 2.33 0.920

4.17 6.0 47.04 11.28 18.0 3.0 0.880

2.78 6.0 37.11 13.35 22.0 3.67 0.820

1.56 6.0 22.71 14.56 26.0 4.33 0.679

1.02 6.0 15.40 15.10 30.0 5.0 0.510

8.33 5.70 91.0 10.92 16.0 2.81 0.94

1.67 1.01 91.0 54.0 16.0 15.85 0.703

5.88 3.99 91.0 15.35 16.0 4.01 0.92

2.63 2.13 72.29 27.18 16.0 7.50 0.813 Partial Angle Measurements

[0238] Although the examples described above generally involve comparing measurements taken at 0° and 90° relative to each other, the techniques described herein also work with measurements that are not a right angles to each other, and also work with measurements that are taken at angles other than 0° and 90° relative to the direction of anisotropy to the material being measured.

[0239] In one experiment, PD as a function of material-PSF orientation was measured in two representative Tl materials with ratio of shear moduli (SR = μΙ_ / μΤ) 1 .72 and 5.9. The degree of anisotropy was successfully derived using PD and SWV in three simulated material using rotation angle = 90°, 30°, 1 5°. The ratios of SWV at 0° to SWVs at 15°, 30°, and 90° were used to calculate SWV derived degree of anisotropy because PD and SWV are inversely related. Degree of anisotropy of material with SR=1 .22 was compared against degree of anisotropy of materials with SR=1 .78 and 2.78. PD was measured as a function of material orientation in the swine psoas major muscle, and the degree of anisotropy was derived using the measured PD and SWV in an isotropic elastic phantom and muscle.

Increasing Axial Measurement Range

[0240] Figure 16 is a graph of the degree of asymmetry of an ARF excitation as a function of depth, where the focal depth is 35mm. As Figure 16 illustrates, the degree of asymmetry of the ARF excitation is relatively constant through depth. For example, in Figure 16, the degree of asymmetry remains at approximately 0.25 from a depth of 5mm to a depth of 33mm. This allows for anisotropy imaging in a large region above the focal depth.

Rotation-less Measurements

[0241] In one embodiment, the variation in degree of asymmetry of an ARF excitation as a function of depth is exploited so as to perform measurements at different relative orientations without the need to physically rotate the ultrasonic transducer. Because the transducer's position is not changed between the two measurements, the two measurements are always properly registered to each other.

[0242] Figure 17 is a graph of the degree of asymmetry of an ARF excitation as a function of depth for five different excitations, each having a different focal depth. At an observation depth of 33mm, setting a focal depth of 33mm creates an ARF having a "tall oval" shape at the observation depth of 33mm; setting a focal depth to 38mm creates an ARF having a "wide oval" shape at the observation depth of 33mm. Thus, by changing the focal depth only, forces having long axes 90° apart can be created without having to rotate the transducer. This technique can provide very accurate results by eliminating errors due to unintended translation of the transducer during rotation, for example.

Experimental and Clinical Results

[0243] The techniques above have been shown to detect skeletal and muscle anisotropy in the presence of Duchenne Muscular Dystrophy (DMD). It has been demonstrated that the ratio of Relative Elasticity (RE) at the focal depth in both Rectus Femorus and Gastrocnemius Medial Head (GM) muscles. The DMD Rectus Femorus was significantly more anisotropic than the control muscles. Also, the transverse shear modulus was higher than the longitudinal shear modulus in the DMD Rectus Femorus, yielding a ratio greater than 1 .0, which is the opposite to that of healthy muscle. In the GM, the anisotropy ratio was not significantly different between DMD and control muscles.

[0244] In another study, the ratios of PDs (PDalong / PDacross) achieved using an asymmetrical F/1 .5 ARF focal configuration were linearly correlated (R2=1 ) to the ratios of SWVs (SWValong / SWVacross) in surgically exposed pig kidneys, in vivo. PD ratios measured in vivo in 3 pig kidneys and in 1 1 renal transplant patients using an asymmetrical F/1 .5 ARF focal configuration is greater than 1 , while the ratios achieved using a symmetrical F/5.0 ARF focal configuration are approximately 1 . In patients, the F/1 .5 and F/5.0 PD ratios were statistically different. PD ratios measured in vivo in 8 renal transplant patients using an asymmetrical F/1 .5 ARF focal configuration at 1 month and three months after transplantation did not statistically differ. The transplants had stable creatinine levels and renal function at one and three months post-transplant; the transplants were excluded if rejection was detected clinically. [0245] Those skilled in the art will recognize improvements and modifications to the preferred embodiments of the present invention. All such improvements and modifications are considered within the scope of the concepts disclosed herein and the claims that follow.