Calibration of Multiple Aperture Ultrasound Probes

• US 20140043933A1
• Filed: 08/12/2013
• Published: 02/13/2014
• Est. Priority Date: 08/10/2012
• Status: Active Grant

• Current Assignee: Maui Imaging, Inc.
• Original Assignee: Bruce R. Ritzi
• Inventors: Belevich, Artem, Call, Josef R., Ritzi, Bruce R., Osborn, Nathan W.
• Granted Patent: US 9,572,549 B2
• Time in Patent Office: Days
• US Class Current: 367/11
• CPC Class Codes: A61B 6/584 determining position of components of the apparatus or device using images of the phantom
• CPC Class Codes: A61B 8/12 in body cavities or body tracts, e.g. by using catheters catheters per se A61M25/00
• CPC Class Codes: A61B 8/4477 using several separate ultrasound transducers or probes
• CPC Class Codes: A61B 8/4488 the transducer being a phased array
• CPC Class Codes: A61B 8/4494 characterised by the arrangement of the transducer elements devices for short-range imaging using particular transducer elements arrangements under G01S15/8909; arrangements of transducers in generation of mechanical vibrations of infrasonic, sonic or ultrasonic frequency B06B1/0607
• CPC Class Codes: A61B 8/483 involving the acquisition of a 3D volume of data
• CPC Class Codes: A61B 8/587 Calibration phantoms models for medical purposes G09B23/28
• CPC Class Codes: G01S 7/52052 with simulation of echoes

Abstract

The quality of ping-based ultrasound imaging is dependent on the accuracy of information describing the precise acoustic position of transmitting and receiving transducer elements. Improving the quality of transducer element position data can substantially improve the quality of ping-based ultrasound images, particularly those obtained using a multiple aperture ultrasound imaging probe, i.e., a probe with a total aperture greater than any anticipated maximum coherent aperture width. Various systems and methods for calibrating element position data for a probe are described.

Claims

1. A method of calibrating an ultrasound probe, comprising the steps of:

placing a first array and a second array of the ultrasound probe in position to image a phantom, each of the first and second arrays having a plurality of transducer elements;

imaging the phantom with the first array to obtain a reference image, wherein imaging is dependent on data describing a position of each transducer element of the first array;

imaging the phantom with the second array to obtain a test image, wherein imaging is dependent on data describing a position of each transducer element of the second array;

quantifying a first error between the reference image and the test image;

iteratively optimizing the data describing the position of each transducer element of the second array until the first error is at a minimum.

2. The method of claim 1, further comprising imaging the phantom with a third array of the ultrasound probe to obtain a second test image, the third array having a plurality of transducer elements, quantifying a second error between the reference image and the second test image and iteratively optimizing data describing a position of each element of the third array until the second error is minimized.

3. The method of claim 1, further comprising storing raw echo data received while imaging the phantom with the second array.

4. The method of claim 3, wherein the iteratively optimizing step comprises:

adjusting the data describing the position of the transducer elements of the second array to create first adjusted position data;

re-beamforming the stored echo data using the first adjusted position data to form a second test image of the reflectors;

quantifying a second error between the second test image and the reference image; and

determining whether the second error is less than the first error.

5. The method of claim 4, wherein adjusting the data describing the position of the transducer elements of the second array includes adjusting a position of a reference point of the array and an angle of a surface of the array, but does not include adjusting a spacing between the elements of the second array.

6. The method of claim 5, further comprising, after a first iteratively optimizing step, performing a second iteratively optimizing step comprising:

adjusting the first adjusted position data, including adjusting a spacing between at least two transducer elements of the second array to create second adjusted position data;

re-beamforming the stored echo data using the second adjusted position data to form a third test image of the reflectors;

quantifying a third error between the third test image and the reference image; and

determining whether the third error is less than the second error.

7. The method of claim 1, wherein iteratively optimizing the transducer element position data comprises optimizing using a least squares optimization process.

8. The method of claim 1, wherein quantifying the first error comprises quantifying a distance between positions of reflectors in the reference image relative to positions of the same reflectors in the test image.

9. The method of claim 1, wherein quantifying the first error comprises quantifying a difference in brightness between reflectors in the reference image and reflectors in the test image.

10. The method of claim 1, wherein quantifying the first error comprises quantifying a difference between a pattern of reflectors and holes in the reference image compared with a pattern of holes and reflectors in the test image.

11. The method of claim 1, wherein the reference image and the test image are three-dimensional volumetric images of a three-dimensional pattern of reflectors, holes, or both reflectors and holes.

12. The method of claim 1, wherein the phantom comprises living tissue.

13. The method of claim 1, further comprising identifying positions of reflectors in the phantom and fitting a mathematically defined curve to a detected pattern of reflectors.

14. The method of claim 13, wherein the curve is a straight line.

15. The method of claim 12, wherein the step of quantifying a first error comprises calculating a coefficient of determination that quantifies a degree of fit of the curve to the pattern of reflectors.

16. A method of calibrating an ultrasound probe, comprising the steps of:

insonifying a plurality of reflectors of a phantom with the ultrasound probe;

receiving echo data with the ultrasound probe;

storing the echo data;

beamforming the stored echo data using first transducer element position data to form an image of the reflectors;

obtaining reference data describing the reflectors;

quantifying an error between the image and the reference data; and

iteratively optimizing the transducer element position data based on the quantified error.

17. The method of claim 16, wherein the iteratively optimizing step comprises iteratively optimizing the transducer element position data with a least squares optimization process.

18. The method of claim 16, wherein the iteratively optimizing step comprises:

adjusting the transducer element position data;

re-beamforming the stored echo data using the adjusted transducer element position data to form a second image of the reflectors;

quantifying a second error based on the second image; and

evaluating the second error to determine whether the adjusted transducer element position data improves the image.

19. The method of claim 18, wherein adjusting the transducer element position data comprises adjusting an array horizontal position variable, an array vertical position variable and an array angle variable.

20. The method of claim 19, wherein adjusting the transducer element position data does not comprise adjusting a spacing between adjacent transducer elements on a common array.

21. The method of claim 16, wherein the reference data is based on physical measurements of the phantom.

22. The method of claim 16, further comprising deriving the reference data from a reference image of the phantom.

23. The method of claim 22, wherein the reference image is obtained using a different group of transducer elements of the probe than a group of transducer elements used for the insonifying and receiving steps.

24. The method of claim 16, wherein the step of iteratively optimizing the transducer element position data comprises using a least squares optimization process.

25. The method of claim 16, further comprising identifying positions of reflectors in the phantom and fitting a mathematically defined curve to a detected pattern of reflectors.

26. The method of claim 25, wherein the curve is a straight line.

27. The method of claim 25, wherein the step of quantifying a first error comprises calculating a coefficient of determination that quantifies a degree of fit of the curve to the pattern of reflectors.

28-44. -44. (canceled)