Calibration of multiple aperture ultrasound probes

• US 9,572,549 B2
• Filed: 08/12/2013
• Issued: 02/21/2017
• Est. Priority Date: 08/10/2012
• Status: Active Grant

• Current Assignee: Maui Imaging, Inc.
• Original Assignee: Maui Imaging, Inc.
• Inventors: Belevich, Artem, Call, Josef R., Ritzi, Bruce R., Osborn, Nathan W.
• Primary Examiner(s): Hulka, James
• Assistant Examiner(s): ARMSTRONG, JONATHAN D
• Application Number:
  US13/964,701
• Publication Number:
  US 20140043933A1
• Time in Patent Office: 1,289 Days
• Field of Search: 367/13
• US Class Current: 1/1
• 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:

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

retrieving first transducer element position data explicitly defining a position of each transducer element of the first array relative to a common coordinate system;

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

retrieving second transducer element position data explicitly defining a position of each transducer element of the second array relative to the common coordinate system;

imaging the phantom with the second array to obtain a test image, wherein imaging is dependent on the second position data;

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

iteratively optimizing the second transducer element position data 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 third position 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 second transducer element position 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 second transducer element position 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 second 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 13, 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:

securing the ultrasound probe having an array of transducer elements in a consistent position relative to a phantom;

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

receiving echo data with the ultrasound probe array transducer elements;

storing the echo data received by each of the transducer elements;

retrieving first transducer element position data explicitly defining a position of each transducer element of the ultrasound probe array relative to a common coordinate system;

beamforming the stored echo data using the 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 second image relative to the first 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, further comprising identifying positions of reflectors in the phantom and fitting a mathematically defined curve to a detected pattern of reflectors.

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

26. The method of claim 24, 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.