Patent Grant
10653392
All publications and patent applications mentioned in this specification are herein incorporated by reference in their entirety to the same extent as if each individual publication or patent application was specifically and individually indicated to be incorporated by reference.
This application relates generally to the field of ultrasound imaging, and more particularly to ping-based ultrasound imaging using apparent point-source transmitters.
In conventional scanline-based ultrasonic imaging, a focused beam of ultrasound energy is transmitted into body tissues to be examined and echoes returned along the same line are detected and plotted to form a portion of an image along the scanline. A complete image may be formed by repeating the process and combining image portions along a series of scanlines within a scan plane. Any information in between successive scanlines must be estimated by interpolation.
The same process has been extended to obtaining ultrasonic images of three-dimensional volumes by combining images from multiple adjacent slices (where each slice is in a different scan plane). Again, any information from any space in between successive scan planes must be estimated by interpolation. Because time elapses between capturing complete 2D slices, obtaining 3D image data for a moving object may be significantly impaired. So-called “4D” imaging systems (in which the fourth dimension is time) strive to produce moving images (i.e., video) of 3D volumetric space. Scanline-based imaging systems also have an inherent frame-rate limitation which creates difficulties when attempting 4D imaging on a moving object.
As a result of these and other factors, some of the limitations of existing 2D and 3D ultrasonic imaging systems and methods include poor temporal and spatial resolution, imaging depth, speckle noise, poor lateral resolution, obscured tissues and other such problems.
Significant improvements have been made in the field of ultrasound imaging with the creation of multiple aperture imaging, examples of which are shown and described in Applicant's prior patents and applications referenced above. Multiple aperture imaging methods and systems allow for ultrasound signals to be both transmitted and received from physically and logically separate apertures.
The various embodiments of systems and methods herein provide the ability to perform high resolution three-dimensional ultrasound imaging at frame rates sufficient to capture details of moving objects. Traditional scanline-based ultrasound imaging methods are limited to relatively slow frame rates due to the need to transmit and receive many scanlines to obtain a single two-dimensional plane. Extending such techniques to obtain imaging data from a complete 3D volume results in even slower frame rates due to the need to image many 2D slices.
As an example, assume that one needs to collect data from a cube of tissue 10 cm on a side at a depth ranging from 5 cm to 15 cm. If scanlines are transmitted from a common center, the shape that would be explored would be a truncated pyramid instead of a shape with comparable thickness in the proximal and distal regions. The tissue may be sampled with beams that are 2 mm (or less) apart on the distal face of the cube. To cover the distal surface one would need at least 50×50 directed beams or 2500 directed pulses. With a maximum pulse rate of approximately 2500 pulses/sec (which may be constrained by the speed of sound in tissue, the expected signal attenuation, and the background noise level), all of the required data may be collected in about one second. This collection time may be adequate for non-moving tissue such as bone, liver, etc., but is not fast enough to capture motion in arteries, or organs such as kidneys and especially the heart, or in moving joints or muscles.
On the other hand, with ping-based imaging, a single ping, propagating substantially uniformly in three dimensions, can insonify the entire volume, and dynamic beamforming (focusing) can identify the sources of the echo returns. Using ping-based imaging techniques, a minimum of three pings may be needed to obtain data for a 3D volume, while a minimum of two pings may be needed to obtain data for a 2D slice. In practical terms, ten to fifty (or more) pings may be used to achieve a desired image quality. For example, the use of 25 pings at a rate of 2500 pings per second may require only 0.01 seconds to acquire all the data for the entire 10 cm cube of tissue. For this particular example, data collection may be 100 times faster than with the scanline-based method.
Using ping-based ultrasound imaging techniques, both 2D and 3D frame rates may be increased substantially so as to allow for imaging of 3D volumes in real-time. Furthermore, by applying multiple aperture imaging techniques (e.g., transmitting and receiving ultrasound signals through multiple, spatially or physically separated acoustic windows), the resolution of such real-time 3D images may be dramatically improved relative to single-aperture techniques.
The following disclosure provides various embodiments of apparent point-source transducers, as well as systems and methods for using such apparent point-source transducers to perform high-frame-rate high-resolution real-time 2D, 3D and so-called 4D ultrasound imaging.
In one embodiment, a method of imaging an object with ultrasound energy is provided, the method comprising the steps of transmitting an un-focused ultrasound signal into a target medium from apparent point-source transmit transducer comprising a shell of piezoelectric material shaped as a spherical section with a spherical center point, receiving echoes reflected by a reflector in the target medium with an omnidirectional receive element that is different than the apparent point-source transmit transducer, determining a position of the reflector within the target medium by obtaining element position data describing a position of the spherical center point of the apparent point-source transmit transducer and a position of the receive element, calculating a total path distance as a sum of a first distance between the spherical center point and the reflector and a second distance between the reflector and the receive element, and determining a locus of possible points at which the reflector may lie, and producing an image of the reflector.
In some embodiments, the receive element comprises a shell of piezoelectric material shaped as a spherical section with a second spherical center point and wherein the position of the receive element is a position of the second spherical center point.
In another embodiment, the position of the receive element lies on a surface of the receive element.
In one embodiment, the method further comprises repeating the receiving, determining, and producing steps with a plurality of receive elements in a common receive aperture.
In one embodiment, the method further comprises repeating the receiving, determining, and producing with elements of a plurality of receive apertures.
In one embodiment, the method further comprises repeating the transmitting step with a separate second, apparent point-source transmit transducer.
In some embodiments, a straight-line distance between the apparent point-source transmit transducer and the receive element is greater than a maximum coherent aperture length for an intended imaging application.
In other embodiments, calculating the total path distance comprises adding apparent path segment representing a distance from a convex transmit transducer surface of the apparent point-source transmit transducer to the spherical center point.
In some embodiments calculating the total path distance comprises subtracting apparent path segment representing a distance from a concave transmit transducer surface of the apparent point-source transmit transducer to the spherical center point.
In alternative embodiments, the receive element has a circular shape.
An ultrasound imaging system is also provided, comprising a first apparent point-source transmit transducer shaped as a spherical section having a spherical center point, the first apparent point-source transmit transducer configured to transmit a three-dimensional semi-spherical pulse into a target object to be imaged, a first plurality of receive transducer elements configured to receive echoes of the three-dimensional semi-spherical pulse, a second plurality of receive transducer elements configured to receive echoes of the three-dimensional semi-spherical pulse, a controller configured to control transmission of the three-dimensional semi-spherical pulse and to determine a position of reflectors within the object based on a known position of the spherical center point of the apparent point-source transmit transducer, known positions of the elements of the first and second pluralities of receive transducer elements, a time at which the three-dimensional semi-spherical pulse was transmitted, and times at which the echoes are received.
In some embodiments, the first apparent point-source transmit transducer is convex relative to the target object.
In one embodiment, the first apparent point-source transmit transducer is concave relative to the target object.
In alternative embodiments, the first apparent point-source transmit transducer is shaped as a spherical section that is greater than half a sphere.
In some embodiments, the first apparent point-source transmit transducer is shaped as a spherical section that is less than half a sphere.
In one embodiment, the first apparent point-source transmit transducer is shaped as a spherical section that is half a sphere.
In some embodiments, the first apparent point-source transmit transducer has a spherical radius of between 0.2 mm and 10 mm.
In one embodiment, the first apparent point-source transmit transducer is configured to transmit ultrasound signals at a first frequency range.
In other embodiments, the first apparent point-source transmit transducer comprises a shell of piezoelectric material with a constant thickness.
In one embodiment, the system further comprises a second apparent point-source transmit transducer with a spherical radius and configured to transmit ultrasound signals at a second frequency range, the second frequency range being different than the first frequency range.
In some embodiments, the apparent point-source transmit transducer comprises a shell having a constant-thickness made of a continuous piezoelectric material.
In another embodiment, the apparent point-source transmit transducer comprises a shell having a constant-thickness made of a segmented piezoelectric material.
In some embodiments, the apparent point-source transmit transducer comprises a plurality of segments arranged into the spherical shape, wherein all segments are configured to transmit ultrasound signals simultaneously.
In one embodiment, the system further comprises a computer readable memory containing data describing a position of the spherical center point of the apparent point-source transmit transducer relative to at least one element of the first plurality of receive transducer elements.
In one embodiment, the system further comprises a computer readable memory containing an adjustment factor representing apparent path segment equal to a distance from a surface of the first apparent point-source transmit transducer to the spherical center point.
In another embodiment, each of the first plurality of receive elements and the second plurality of receive elements has a circular shape.
An ultrasound probe comprising: an apparent point-source transmit transducer comprising a shell of piezoelectric material shaped as a spherical section with a constant wall thickness and a spherical center point; and a receive array comprising a plurality of omnidirectional receive transducer elements, the receive array having a total aperture greater than a coherence width for an intended imaging application.
In some embodiments, the plurality of receive transducer elements are grouped into separate arrays.
In another embodiment, the plurality of receive transducer elements are contained in a continuous array.
In some embodiments, the receive elements have a cylindrical shape.
In additional embodiments, the receive elements have a spherical section shape.
In alternative embodiments, the ultrasound probe is sized and configured for insertion into a body lumen or cavity.
In one embodiment, the ultrasound probe is sized to cover approximately half of a human patient's chest.
In some embodiments, the total aperture is at least twice a coherence width for the intended imaging application.
In other embodiments, the total aperture is at least three times a coherence width for the intended imaging application.
In one embodiment, the probe comprises an array of transducer elements with a width of about 8 cm to about 10 cm.
An apparent point-source ultrasound transducer element is provided, comprising a shell of piezoelectric material shaped as a spherical section with a constant wall thickness and a spherical center point, a convex surface, and a concave surface, an acoustic damping material surrounding and bonded to the convex surface of the shell, and an electrical lead extending through the acoustic damping material and connected to the convex surface of the shell.
In one embodiment, the element further comprises an acoustic matching material filling and bonded to the concave surface of the shell.
In some embodiments, the shell has a transmitting surface shaped as a spherical section that is greater than half a sphere.
In other embodiments, the shell has a transmitting surface shaped as a spherical section that is less than half a sphere.
In additional embodiments, the shell has a transmitting surface shaped as a spherical section that is half a sphere.
In some embodiments, the shell is made of a composite material comprising a piezoelectric ceramic and a polymer.
In other embodiments, the shell comprises lead zirconate titanate (PZT).
In one embodiment, the shell has a transmitting surface area of at least three square millimeters.
In another embodiment, the shell has a transmitting surface area of at least five square millimeters.
In one embodiment, the shell has a transmitting surface area of at least ten square millimeters.
An apparent point-source ultrasound transducer element is also provided, comprising a shell of piezoelectric material shaped as a spherical section with a constant wall thickness and a spherical center point, a convex surface and a concave surface, an acoustic damping material filling and bonded to the concave surface of the shell, and an electrical lead extending through the acoustic damping material and connected to the concave surface of the shell.
In one embodiment, the element further comprises an acoustic matching material bonded to the convex surface of the shell.
In some embodiments, the shell has a transmitting surface shaped as a spherical section that is greater than half a sphere.
In another embodiment, the shell has a transmitting surface shaped as a spherical section that is less than half a sphere.
In one embodiment, the shell has a transmitting surface shaped as a spherical section that is half a sphere.
In an alternative embodiment, the shell is made of a composite material comprising a piezoelectric ceramic and a polymer.
In some embodiments, the shell comprises lead zirconate titanate (PZT).
In other embodiments, the shell has a transmitting surface area of at least three square millimeters.
In one embodiment, the shell has a transmitting surface area of at least five square millimeters.
In another embodiment, the shell has a transmitting surface area of at least ten square millimeters.
The novel features of the invention are set forth with particularity in the claims that follow. A better understanding of the features and advantages of the present invention will be obtained by reference to the following detailed description that sets forth illustrative embodiments, in which the principles of the invention are utilized, and the accompanying drawings of which:
The various embodiments will be described in detail with reference to the accompanying drawings. References made to particular examples and implementations are for illustrative purposes, and are not intended to limit the scope of the invention or the claims.
The present disclosure provides systems and methods for improving the quality of 2D, 3D and 4D ultrasound images through the use of one or more apparent point-source ultrasound transmitters. In some embodiments, such apparent point-source transmitters may be used in combination or integrally with multiple aperture ultrasound imaging systems, multiple aperture ultrasound probes, and/or multiple aperture ultrasound beamforming techniques. Various embodiments of such systems, methods and combinations are provided herein.
Although the various embodiments are described herein with reference to ultrasound imaging of various anatomic structures, it will be understood that many of the methods and devices shown and described herein may also be used in other applications, such as imaging and evaluating non-anatomic structures and objects. For example, the various embodiments herein may be applied to non-destructive testing applications such as evaluating the quality, integrity, dimensions, or other characteristics of various structures such as welds, pressure vessels, pipes, structural members, beams, etc. The systems and methods may also be used for imaging and/or testing a range of materials including human or animal tissues, solid metals such as iron, steel, aluminum, or titanium, various alloys or composite materials, etc.
The following paragraphs provide useful definitions for some terms used frequently herein. Other terms may also be defined as they are used below.
As used herein the terms “ultrasound transducer” and “transducer” may carry their ordinary meanings as understood by those skilled in the art of ultrasound imaging technologies, and may refer without limitation to any single component capable of converting an electrical signal into an ultrasonic signal and/or vice versa. For example, in some embodiments, an ultrasound transducer may comprise a piezoelectric device. In other embodiments, ultrasound transducers may comprise capacitive micro-machined ultrasound transducers (CMUT) or any other transducing device capable of converting ultrasound waves to and from electrical signals.
Transducers are often configured in arrays of multiple individual transducer elements. As used herein, the terms “transducer array” or “array” generally refers to a collection of transducer elements mounted to a common backing plate. Such arrays may have one dimension (1D), two dimensions (2D), 1.X dimensions (1.XD) or three dimensions (3D) as those terms are used elsewhere herein and/or as they are commonly understood in the art. Other dimensioned arrays as understood by those skilled in the art may also be used. Annular arrays, such as concentric circular arrays and elliptical arrays may also be used. An element of a transducer array may be the smallest discretely functional component of an array. For example, in the case of an array of piezoelectric transducer elements, each element may be a single piezoelectric crystal or a single machined section of a piezoelectric crystal.
As used herein, the terms “transmit element” and “receive element” may carry their ordinary meanings as understood by those skilled in the art of ultrasound imaging technologies. The term “transmit element” may refer without limitation to an ultrasound transducer element which at least momentarily performs a transmit function in which an electrical signal is converted into an ultrasound signal. Transmitted ultrasound signals may be focused in a particular direction, or may be unfocused, transmitting in all directions or a wide range of directions. Similarly, the term “receive element” may refer without limitation to an ultrasound transducer element which at least momentarily performs a receive function in which an ultrasound signal impinging on the element is converted into an electrical signal. Transmission of ultrasound into a medium may also be referred to herein as “insonifying.” An object or structure which reflects ultrasound waves may be referred to as a “reflector” or a “scatterer.”
As used herein, the term “aperture” may refer to a conceptual “opening” through which ultrasound signals may be sent and/or received. In actual practice, an aperture is simply a single transducer element or a group of transducer elements that are collectively managed as a common group by imaging control electronics. For example, in some embodiments an aperture may be a grouping of elements which may be physically separate and distinct from elements of an adjacent aperture. However, adjacent apertures need not necessarily be physically separate or distinct. Conversely, a single aperture may include elements of two or more physically separate or distinct transducer arrays. For example, distinct groups of transducer elements (e.g., a “left aperture”) may be constructed from a left array, plus the left half of a physically distinct center array, while a “right aperture” may be constructed from a right array, plus the right half of a physically distinct center array).
It should be noted that the terms “receive aperture,” “insonifying aperture,” and/or “transmit aperture” are used herein to mean an individual element, a group of elements within an array, or even entire arrays, that perform the desired transmit or receive function from a desired physical viewpoint or aperture. In some embodiments, such transmit and receive apertures may be created as physically separate components with dedicated functionality. In other embodiments, any number of send and/or receive apertures may be dynamically defined electronically as needed. In other embodiments, a multiple aperture ultrasound imaging system may use a combination of dedicated-function and dynamic-function apertures.
As used herein, the term “total aperture” refers to the overall size of all imaging apertures in a probe. In other words, the term “total aperture” may refer to one or more dimensions defined by a maximum distance between the furthest-most transducer elements of any combination of send and/or receive elements used for a particular imaging cycle. Thus, the total aperture may be made up of any number of sub-apertures designated as send or receive apertures for a particular cycle. In the case of a single-aperture imaging arrangement, the total aperture, sub-aperture, transmit aperture, and receive aperture may all have the same dimensions. In the case of a multiple aperture imaging arrangement, the dimensions of the total aperture includes the sum of the dimensions of all send and receive apertures plus any space between apertures.
In some embodiments, two apertures may be located adjacent to one another on a continuous array. In still other embodiments, two apertures may overlap one another on a continuous array, such that at least one element functions as part of two separate apertures. The location, function, number of elements and physical size of an aperture may be defined dynamically in any manner needed for a particular application. Constraints on these parameters for a particular application will be discussed below and/or will be clear to the skilled artisan.
Elements and arrays described herein may also be multi-function. That is, the designation of transducer elements or arrays as transmitters in one instance does not preclude their immediate re-designation as receivers in the next instance. Moreover, embodiments of the control system herein include the capabilities for making such designations electronically based on user inputs, pre-set scan or resolution criteria, or other automatically determined criteria.
Introduction to Point-Source Transmission Ultrasound Imaging
In various embodiments, point-source transmission ultrasound imaging, otherwise referred to as ping-based ultrasound imaging, provides several advantages over traditional scanline-based imaging. Point-source transmission differs in its spatial characteristics from a “phased array transmission” which focuses energy in a particular direction from the transducer element array along a directed scanline. A point-source pulse (ping) may be transmitted so as to generate either a two-dimensional a circular wavefront or a three-dimensional spherical wavefront in the scanning plane, thereby insonifying as wide an area as possible. Echoes from scatterers in the region of interest may return to all of the elements of receive apertures. Those echo signals may be filtered, amplified, digitized and stored in short term or long term memory (depending on the needs or capabilities of a particular system).
Images may then be reconstructed from received echoes by assuming that the wavefronts emitted from the point-source are physically circular in the region of interest. In actuality, the wavefront may also have some penetration in the dimension orthogonal to the scanning plane (i.e., some energy may essentially “leak” into the dimension perpendicular to the desired two-dimensional scanning plane, reducing the effective imaging depth). Additionally, the “circular” wavefront may actually be limited to a semicircle or a fraction of a circle less than 180 degrees ahead of the front face of the transducer according to the unique off-axis properties of the transducing material used. Similarly, a “spherical” wavefront may have an actual shape of a hemisphere or less than a hemisphere within the medium to be imaged.
A software-based, firmware-based, or hardware-based dynamic beamforming technique, in which a beamformer's focus may be continuously changed to focus at a particular pixel position as that pixel is being imaged, may be used to plot the position of echoes received from a point-source pulse. In some embodiments, a dynamic beamformer may plot the locus of each echo signal based on a round-trip travel time of the signal from the transmitter to an individual receive transducer element.
The locus of a single reflector will lie along either a two-dimensional ellipse (in the case of two-dimensional imaging) or a three-dimensional ellipsoid (in the case of three-dimensional imaging). A first focus of the ellipse or ellipsoid will be at the position of the transmit transducer element and the second focus will be at the position of the receive transducer element. Although several other possible reflectors lie along the same ellipse or ellipsoid, echoes of the same reflector will also be received by each of the other receive transducer elements of a receive aperture. The slightly different positions of each receive transducer element means that each receive element will define a slightly different ellipse or ellipsoid for a given reflector. Accumulating the results by coherently summing the ellipses or ellipsoids for all elements of a common receive aperture will indicate an intersection of the ellipses or ellipsoids for a reflector, thereby converging towards a point at which to display or define a pixel or voxel representing the reflector. The echo amplitudes received by any number of receive elements may thereby be combined into each pixel or voxel value. In other embodiments the computation can be organized differently to arrive at substantially the same image.
Various algorithms may be used for combining echo signals received by separate receive elements. For example, some embodiments may process echo-signals individually, plotting each echo signal at all possible locations along its ellipse, then proceeding to the next echo signal. Alternatively, each pixel location may be processed individually, identifying and processing all echoes potentially contributing to that pixel location before proceeding to the next pixel location.
Image quality may be further improved by combining images formed by the beamformer from one or more subsequent transmitted pings, transmitted from the same or a different point-source (or multiple different point-sources). Still further improvements to image quality may be obtained by combining images formed by more than one receive aperture. An important consideration is whether the summation of images from different pings, different transmit point-sources or different receive apertures should be coherent summation (phase sensitive) or incoherent summation (summing magnitude of the signals without phase information).
The decision as to whether to use coherent or incoherent summation may be influenced by the lateral extent/size of the receive aperture(s) and/or the transmit aperture(s). In some embodiments, it may be convenient to confine the size of an aperture to conform to the assumption that the average speed of sound is substantially the same for every path from a scatterer to each element of the receive aperture. For narrow receive apertures this simplifying assumption is easily met. However, as the width of the receive aperture increases, an inflection point is reached (referred to herein as the “maximum coherent aperture width” or “maximum coherence width”), beyond which the paths traveled by returning echoes of a common reflector will necessarily pass though different types of tissue having intrinsically different speeds of sound. When this difference results in receive wavefront phase shifts approaching or exceeding 180 degrees, additional receive elements extended beyond the maximum coherence width will actually degrade the image rather than improve it. The same considerations may also apply to the size of transmit apertures, which may include a plurality of coherently combined transducer elements. In the case of two-dimensional transducer arrays used in three-dimensional imaging (or 3D data collection), it may be useful to define a maximum coherent aperture size in two dimensions. Thus, in various embodiments a maximum coherent aperture may be defined as a group of transducer elements in a square, circle, polygon or other two-dimensional shape with a maximum distance between any two elements such that phase cancellation will be avoided when echo data received at the elements of the aperture are coherently combined.
Therefore, in order to realize the inherent benefits (e.g., in terms of increased spatial resolution) of a wide probe with a total aperture width far greater than the maximum coherent aperture width, the full probe width may be physically or logically divided into multiple apertures, each of which may be limited to an effective width less than or equal to the maximum coherent aperture width, and thus small enough to avoid phase cancellation of received signals. The maximum coherence width can be different for different patients (or different test objects) and for different probe positions on the same patient. In some embodiments, a compromise width may be determined for a given probe system. In other embodiments, a multiple aperture ultrasound imaging control system may be configured with a dynamic algorithm to subdivide the available elements in multiple apertures into groups that are small enough to avoid significant image-degrading phase cancellation. In various embodiments, a particular coherent aperture size may be determined automatically by a control system, or manually through user input via a user control such as a dial or slider.
In some embodiments, coherent (phase sensitive) summation may be used to combine echo data received by transducer elements located on a common receive aperture resulting from one or more pings. In some embodiments, incoherent summation may be used to combine echo data or image data received by separate receive apertures if such receive apertures that could possibly contain phase-cancelling data. Such may be the case with receive apertures that have a combined total aperture that is greater than a maximum coherence width for a given imaging target.
Point-Source Transmission for 3D Ultrasound Imaging
When a three-dimensional pulse is initiated from a point-source transmit transducer, the resulting semi-spherical wavefront travels into the region of interest (ROI) where some of the ultrasound energy may be reflected by scatterers in the ROI. Some of the echoes from the scatterers may travel back towards receive transducer elements of the probe, where the echoes may be detected, amplified, digitized, and stored in a short-term or long-term memory device. Each digitized sample value may represent a scatterer from the ROI. As in the 2D case, the magnitude of each received sample, along with its time of arrival and the exact positions of the transmit and receive transducers used, may be analyzed to define a locus of points identifying potential positions of the scatterer. In the 3D case, such a locus is an ellipsoid having as its foci the positions of the transmit and receive transducers. Each unique combination of transmit and receive transducer elements may define a separate view of the same reflector. Thus, by combining information from multiple transmit-receive transducer combinations, the actual location of each reflector may be more accurately represented.
For example, in some embodiments an image in a 3D array of voxels may be assembled in computer memory by beginning with an evaluation of a selected digital sample. The selected digitized sample value may be written into every voxel indicated by the corresponding ellipsoid described above. Proceeding to do the same with every other collected sample value, and then combining all resulting ellipsoids may produce a more refined image. Real scatterers would be indicated by the intersection of many ellipsoids whereas parts of the ellipsoids not reinforced by other ellipsoids would have low levels of signal and may be treated as noise (i.e., eliminated or reduced by filters or other image processing steps).
In other embodiments, the order of computation may be changed by beginning with a selected voxel in a final 3D image volume to be produced. For example, for a selected voxel, the closest stored sample may be identified for each transmitter/receiver pair. All samples corresponding to the selected voxel may then be evaluated and summed (or averaged) to produce a final representation of the voxel. Closeness of a sample to a selected voxel may be determined by calculating the vector distance from the three-dimensional position of a transmitter (i.e., the transmitter used to produce the sample) to the selected voxel position plus the vector distance from the selected voxel position to the position of a receiver used to produce the sample. Such a linear distance may be related to the time-divided sample values by dividing the total path length by speed of sound through the imaged object. Using such a method, the samples corresponding to a calculated time may be associated with the selected voxel.
Techniques for determining the location for received echo samples are generally referred to herein as beamforming, while techniques for combining information obtained from multiple transmitter/receiver combinations or from multiple separate pings transmitted using the same transmitter/receiver combination may generally be referred to as image layer combining. In various embodiments, a frame may be made up of any number of combined image layers. Frames may be displayed sequentially at a desired frame-rate on a display to form a moving image or video. The above-described beamforming processes may beneficially also be used for evaluating pixel values in a 2D cross-sectional slice through a 3D volume using raw echo data. In various embodiments, such 2D slices may be taken at any arbitrary angle or along any curved path through the 3D volume. The same techniques may also be used to zoom in (i.e., increase the size of features) using raw echo data rather than enlarging processed pixels or voxels.
Apparent Point-Source Transmitters
As described above, a point-source transmitter may be approximated using a single small transducer element of a transducer array. When performing 2D ping imaging using a 1D array (an array of elements with parallel longitudinal axes, typically including a lens to focus the signal into a single imaging plane), a single element may be able to produce a ping with sufficient energy in the imaging plane to achieve imaging at a reasonable depth. However, when imaging is purposefully extended into the third dimension, a single small transmit element of a typical transducer array may be insufficient to produce a ping with enough energy to obtain a viable image at the desired depth due to insufficient signal power. This may be understood in view of the fact that the power of a transmitted ultrasound pulse is dispersed in three dimensions rather than two, so the log amplitude of the wavefront attenuates according to an inverse-square relation rather than linearly. Depending on the frequency of the transmitted pulse and the attenuation rate of the material under observation, a low energy ping may weaken beneath the background noise level before returning a usable signal at desired depth. One solution may be to transmit a “ping” from multiple adjacent elements, but the more elements used, the less the transmit aperture approximates a point-source, which may have the effect of distorting the semi-spherical shape of the transmitted waveform (or semi-circular shape in the 2D case), which may result in reduced image quality. Using multiple transmit elements also reduces precision in the determination of a point to use as the transmit-source ellipsoid focus during beamforming calculations, thereby further reducing image quality. Such reduced image quality may be acceptable in some applications, but in other applications a higher quality image may be desired.
In various embodiments, an “apparent point-source transmitter” transducer may be configured to produce a waveform that both approximates an actual point-source and has sufficient energy to produce high quality images at the desired depth. In some cases, such apparent point-source transmitters may be configured such that ultrasound power output may be limited only by safety considerations within the imaged medium.
As used herein, the phrase “point-source” refers to a point in 3D space that represents a center point of a transmitted 2D or 3D ultrasound waveform. In some embodiments, such a point is ideally an infinitely small point corresponding to a produced wavefront with a consistent semi-spherical shape. In embodiments in which such a waveform is produced by a single small element, such a point may lie on the surface of the transducer element. As used herein, the terms “semi-spherical pulse” and “semi-spherical wavefront” may refer to any ultrasound wavefront with a spherical-section shape, including wavefronts with approximately spherical-section shapes greater than or less-than an ideal semi-sphere. Similarly, the terms “semi-circular pulse” and “semi-circular wavefront” may refer to any ultrasound wavefront which appears in an imaging plane to have a circular-section shape, including wavefronts with approximately circular-section shapes greater than or less-than an ideal semi-circle.
In some cases, multiple (e.g., two, three, four or more) small transducer elements from a common transmit/receive array may be excited simultaneously to produce a ping with more energy than may be produced by a single element. As a practical matter, when using multiple small elements to approximate a point-source transmitter, the “point” may effectively be larger and more spread out, which may tend to cause a loss of precision in beamforming calculations using the center “point” as the location of the transmitted ping (and, by extension, as one of the foci—along with the location of the receive element—of the ellipsoid representing the locus of all points for any given time-of-travel). Such a spread out or “smeared” transmit point may also lead to potentially undesirable variation in the shape of the produced waveform from an ideal semi-sphere. Some variation may be inevitably present in any point-source transmitter, but better results may be achieved with a point-source transmitter that produces a waveform as close to the ideal as possible.
An alternate solution is to provide a large transducer shaped and configured to produce a relatively high-power waveform that “appears” to have originated from a point-source—in other words, apparent point-source. When performing beamforming calculations to determine the location of reflectors based on the timing of received echoes, the location of the apparent point-source may be used as the origin of the transmitted ping wavefront. In some embodiments, particularly suitable shapes for transmitting transducers may include concave and convex spherical caps. Convex spherical caps may generally be referred to herein as “dome-shaped,” while concave spherical caps may be referred to as “bowl-shaped.” Some examples of imaging probes incorporating examples of such transducer elements are provided below.
An apparent point-source may exist when the point defining the origin of the semi-spherical wavefront lies somewhere other than on the surface of the transducer producing the wavefront. If the medium is assumed to be below or in front of the transducer, apparent point-source located above or behind a transducer surface may be referred to herein as a “negative” apparent point-source. On the other hand, apparent point-source located below the transducer surface may be referred to herein as a “positive” apparent point-source. A transducer configured to produce a wavefront that appears to have originated from apparent point-source may be referred to herein as an “apparent point-source transducer.”
An apparent point-source transducer 10 may comprise a shell 14 of a material exhibiting piezoelectric properties. The shell 14 may have a substantially constant thickness throughout. The transducer 10 may further include one or more electrical conductors extending from an interior surface 18 of the transducer shell 14. In the case of a dome-shaped transducer, the concave volume within the shell 14 may be filled with an acoustic damping material. Examples of suitable acoustic damping materials include polyurethanes, acrylics, epoxies (e.g., doped epoxies, such as tungsten-doped epoxy) or any other suitable acoustic backing materials.
In theory, a transducer in the shape of a complete sphere may produce a perfectly spherical wavefront with an apparent origin at the center of the sphere. However, the need for mechanical and electrical control of a transducer in a practical imaging system necessitates truncating the sphere to some degree. Thus, in some embodiments, a convex dome-shaped apparent point-source transducer 10 such as that shown
The diagrams of
The intersection 28 of the truncation plane 24 and the sphere 23 will be referred to herein as a cut circle which has a cut radius that is mathematically related to the spherical radius according to the equation:
a=sqrt(R2−E2)
Where a is the cut radius, R is the spherical radius and E is the cut elevation. R-E is the height (h) of the spherical cap.
The surface area of the resulting spherical cap is also mathematically related to the spherical radius (R) and the cut elevation (E) according to the equation:
Acap=2*π*R*(R−E)
The spherical radius used in the above equations should be the radius to the intended active transducer surface. Thus, for a convex dome-shaped transducer made from a transducer shell with a thickness t, the transducer surface area may be calculated using the outer spherical radius, which is equal to the inner radius plus the thickness.
The perpendicular distance 26 between the truncation plane 24 and a parallel plane 27 through the spherical center 16 will be referred to herein as the cut elevation. Cut elevation may be expressed as an actual distance or as a percent of the spherical radius 29. A cut elevation of exactly zero corresponds to a perfectly semi-spherical cap, while a cut elevation of 99% would result in a very small cap section with a surface area of about 0.5% of a complete sphere. As used herein, a positive cut elevation refers to a spherical cap such as that shown in
In any event, when performing ping-based, non-focused ultrasound imaging using an ultrasound transducer having the shape of a spherical cap, the spherical center point 16 may be treated as the mathematical origin of a wavefront emitted by the transducer for purposes of triangulation. The same may also be applied to convex (bowl-shaped) transducers.
In cases where living human or animal tissue is to be imaged, it may be desirable to keep the apparent point-source of a bowl transducer (i.e., the spherical center point at which ultrasound waves converge) from occurring too near or inside the living tissue. In some embodiments, this may be achieved by selecting appropriate spherical cap dimensions and/or by assembling probes with bowl-shaped apparent point-source transducers with one or more matching layers or other materials with a thickness sufficient to include the spherical center point. Thus, in some embodiments, a concave region of a bowl-shaped apparent point-source transducer 30 may be filled with an acoustic-coupling material that may be selected to have an inherent speed-of-sound substantially matching the medium to be imaged. The coupling material, which may also be referred to as a matching layer, may comprise any suitable material, such as saline solutions, glycerine, propylene glycol, silicone (e.g., RTV silicone), ballistic gelatin or other matching layer or lensing layer materials known to be suitable for use in ultrasound imaging of humans or animals. Alternatively, materials such as acrylics, glass, metals, composites, etc. may be used in matching layers for NDT (Non-Destructive Test) imaging of mechanical, structural or other non-living objects. In some embodiments, such a matching material may extend beyond the ring edge of the transducer shell 14 sufficiently to include the spherical center point 16, thereby eliminating any potential risk that may be presented by ultrasound energy converging at that point within an imaged medium.
In some embodiments, a larger apparent point-source transducer may be capable of inducing a higher energy wavefront in an imaged medium. In general, the maximum energy or power that may be produced by apparent point-source transducer may be proportional to the surface area of the transducer. The actual power produced by a transducer in a particular application may be controlled by varying the magnitude, frequency, duration, duty cycle, or other characteristics of an applied voltage signal. As a result, larger apparent point-source transducers may be used to transmit 3D ultrasound pings with more energy than smaller apparent point-source transducers.
The exact size of apparent point-source transducer may be partially dependent on the application in which it is to be used. Ultrasound signals attenuate as they pass through materials being imaged. As a result, the transmitted signal must have enough power that it may travel into the medium, reflect off of structures to be imaged, and return to the receive transducers with sufficient power that the signal may be adequately distinguished from noise. Thus, on one hand, it is desirable to provide the capability of transmitting ultrasound signals with as much power as possible. On the other hand, practical factors may limit the power level that may be safely used before causing injury (e.g., to human or animal patients being imaged) or damage (e.g., to sensitive materials or equipment being imaged or tested).
Because a desired maximum transducer power may be proportional to the transducer's surface area, the spherical radius and/or the cut elevation of apparent point-source transmit transducer may be selected based on a desired transducer surface area. For example, a 1D transducer element with a length of 14 mm and a width of 0.25 mm has a surface area of 3.5 mm2. If it is desired to make an equivalent apparent point-source transmitter, the same surface area may be achieved with an embodiment of a spherical cap apparent point-source transducer having a cut elevation of zero and a spherical radius of about 0.75 mm. In another embodiment, the same 3.5 mm2 surface area may also be achieved with a spherical cap apparent point-source transducer having a spherical radius of about 0.8 mm and a cut elevation of about 10% of the spherical radius (i.e., a cap height of about 0.7 mm and a cut radius of about 0.78 mm).
In various embodiments, any of various attributes such as the transducer surface area, the cut radius, the spherical radius, cap height, cut elevation, etc. may be used as a design starting point. In some cases, a particular surface area may be desired so as to achieve a desired transmit power level. Various examples of apparent point-source geometries based on various surface areas are provided in Table 1 below.
Alternatively, factors such as probe geometry or intended imaging target may be more easily met by designing apparent point-source transmitters based on spherical radius, cap height, cut radius or other geometric factors.
In other embodiments, apparent point-source transducers of different sizes may be used for imaging at different depths. In some cases, larger transducers may also be more susceptible to manufacturing variation. Such variation may lead to transducers that create non-uniform wavefronts. In some cases, the degree to which transducer surface irregularities may negatively affect imaging performance may be a function of the ultrasound wavelength being used. For example, higher frequency ultrasound (often best suited for relatively shallow-depth imaging due to typically greater attenuation as a function of imaging depth) may require a more accurate spherical surface than lower frequencies which may be better suited to deeper imaging. The term “near field” may generally refer to a region of the image plane nearest to the transducer. Thus, in some embodiments, relatively larger transducers may be used for imaging mid-field and/or far-field regions of a patient or object, while relatively smaller transducers may be used for imaging near-field regions.
For example, a smaller spherical cap apparent point-source transducer with a spherical radius of up to about 0.75 mm may be well suited to imaging near-field regions, and may be configured to transmit at relatively high frequencies (e.g., between about 5 MHz and about 10 MHz or more) for imaging at relatively shallow depths (e.g., to about 5-10 cm in human tissue). In other embodiments, a relatively larger apparent point-source transducer (e.g., with a spherical radius between about 0.75 mm and about 6 mm) may be well suited to imaging somewhat deeper regions, and may be operated to transmit at relatively low frequencies (e.g., between about 1 MHz and about 5 MHz) for imaging relatively deeper regions (e.g., greater than 10 cm).
Thus, in various embodiments, apparent point-source probes for use with ping-based multiple aperture ultrasound imaging techniques may contain one or more spherical cap apparent point-source transducers with a spherical radius of between about 0.2 mm and about 10 mm or more.
In some embodiments, one or more apparent point-source transmit transducers within an ultrasound probe may be operated at different power levels and/or at different frequencies when operating in different imaging modes in order to optimally image at a wide range of depths. In some embodiments, such imaging modes may be manually selected by an operator of an imaging system, and in other embodiments, such modes may be automatically selected based on a pre-programmed imaging process for a chosen imaging scenario.
Examples of Piezoelectric Materials and Manufacturing
As described above, a dome-shaped or bowl-shaped transducer may be in the form of a thin shell of a piezoelectric material in the shape of a truncated spherical cap. Such a shell may be made of any material exhibiting piezoelectric properties. Many naturally occurring and synthetic materials are known to exhibit piezoelectric properties that may be of a character suitable for use in ultrasound imaging applications. In the case of ping-based multiple aperture ultrasound imaging, ultrasound ping signals may be transmitted at frequencies commonly used in diagnostic medical ultrasound, e.g., in the range of about 1 MHz to about 20 MHz or more. Thus, apparent point-source transducers with fundamental frequencies within this range may be suitable for use in ping-based multiple aperture imaging.
Naturally-occurring piezoelectric materials include quartz, topaz and tourmaline, while man-made piezoelectric ceramic materials include lead zirconate titanate (PZT), barium titanate, lead metaniobate, & polyvinylidene difluoride (PVF2—not naturally piezoelectric, but may be made so by heating in the presence of a strong electrical field). Some man-made piezoelectric ceramic materials may be combined with non-piezoelectric polymer materials to create piezo-composites.
The thickness of apparent point-source transducer shell (whether bowl-shaped or dome-shaped) may be directly related to the fundamental frequency of the transducer. In some cases (e.g., for some piezoelectric ceramic materials), the thickness of a transducer shell may be equal to about half a wavelength of its corresponding fundamental frequency, or an odd number of wavelength halves such as 3/2 wavelength or 5/2 wavelength. However, depending on the materials used, the shell thickness may be differently related to the transducer's fundamental frequency. Manufacturing processes may also vary depending on the piezoelectric material used and other factors.
In order to produce a spherical-section shell with a substantially constant thickness, if there is a requirement that the shell have a thickness of half a wavelength, then there may be a minimum size for apparent point-source transducer configured for a particular fundamental frequency. For example, apparent point-source transducer sized for a fundamental frequency of 3 MHz may have a shell thickness of approximately ¼ mm (assuming a speed-of-sound of about 1550 m/s) and may have a minimum external diameter of about 1 mm. In other cases, smaller external diameters may be possible by using a different material, designing for a different speed-of-sound application, etc.
In some cases, the speed-of-sound characteristics of the piezoelectric material of the transducer itself may have directional characteristics (e.g., sound may travel faster along one crystal axis than another). In such cases, the shape of apparent point-source transducer may be varied from an ideal physical sphere (and/or by varying the transducer material thickness in portions of the sphere) in order to create a transducer that produces a more uniform spherical-section wavefront in the material to be imaged.
For example, natural or man-made piezoelectric material may be machined using traditional techniques in order to form the desired shape directly from a block of material. Such machining may be performed using mechanical cutters, water jets or any other available machining technique. Alternatively, a block or sheet of piezoelectric material may be machined into a plurality of elements attached to a flexible substrate which may then be formed into the desired shape. For example, a plurality of concentric ring cuts 42 and radial cuts 44 may be made in a sheet of piezoelectric material (as shown for example in
In some embodiments, a desired shape may be molded (e.g., by injection molding, die casting, or other molding process) from a piezo-composite material. Examples of molding processes that may be adapted to forming spherical-cap elements are described in U.S. Pat. Nos. 5,340,510 and 5,625,149.
It is also possible to produce ultrasound transducers in a desired shape using additive manufacturing techniques (commonly known as 3D printing techniques). For example, US Patent Publication 2013/0076207 and US Patent Publication 2013/0088122 describe systems and methods for forming transducers in the shape of cylindrical posts. Similar techniques may also be adapted to form transducers with spherical-cap shapes. Additionally, other manufacturing techniques such as laser sintering, stereo lithography, chemical vapor deposition or any other suitable techniques may be used to produce transducers in the shapes and sizes described herein.
Capacitive Micromachined Ultrasound Transducer (CMUT) formation techniques may also be used to form transducers of desired shapes onto a pre-shaped substrate. WO 2012/112540 shows and describes some examples of structures and techniques that may be adapted to forming spherical-cap shaped transducers. Alternately, a dome-shaped transducer may be made by forming an array of CMUT transducers on a substrate pre-formed into a desired shape (e.g., a concave or convex spherical cap as described above). In such embodiments, the CMUT elements may be electrically connected so as to transmit simultaneously without phasing.
Multiple Aperture Probes with Apparent Point-Source Transmitters
Apparent point-source transmitter transducers may be integrated into multiple aperture ultrasound probes designed for 2D or 3D imaging. Ultrasound probes configured for 3D ping-based multiple aperture imaging may comprise transducers arranged in an array extending substantial lengths in at least two dimensions across the surface of the object to be imaged. In some embodiments, some 3D probes may be used for apparent 2D imaging by only displaying a 2D slice while capturing data for a complete 3D volume, as described in further detail below.
In some embodiments, the probe 70 of
As in the example of
In various embodiments, transducer elements configured for receiving echoes of three-dimensional pings transmitted from apparent point-source transducer (e.g., receive elements 84) may themselves have spherical-section shapes defining apparent point receive elements. Receive transducer elements with spherical-section shapes may be well-suited to receiving echoes from all directions (i.e., omnidirectionally, or from as many directions as possible based on other constraints) from within a three-dimensional volume of the imaged medium. In some embodiments, apparent point-source transducer probe may comprise an array of apparent point-source transducers, each of which may comprise an independent concave or convex spherical-section shell. Such an array of apparent point-source transducer elements may include some elements that may be substantially larger than other elements. Alternatively, an array of apparent point-source transducer elements may be substantially the same size as one another. When receive elements are small, other shapes may also provide sufficient omni-directionality for uniformly receiving echoes from a complete volume.
Other receive-element shapes may also be used, depending on the size of such receive elements, and trade-offs such as that between probe manufacturability and image quality. Thus, in some embodiments transducer elements configured for receiving echoes of three-dimensional pings transmitted from apparent point-source transducer may alternatively have substantially circular receiving surfaces that may be planar, convex or concave. In elevation, such elements may have cylindrical or otherwise shaped sides. In further embodiments, transducer elements configured for receiving echoes of three-dimensional pings transmitted from apparent point-source transducer may have planar, concave or convex receiving surfaces with other shapes such as polygonal (e.g., square, hexagonal, octagonal, etc.).
Therefore, some probes may include a combination of transducers with different shapes and sizes in a variety of configurations. For example, in some embodiments, apparent point receive transducer elements may be arranged in arrays made up of aligned rows and columns, offset rows and/or columns, radial lines, rings, random arrays, or any other configuration. Apparent point receive elements may generally be substantially smaller than apparent point-source transmitters in the same probe, while retaining substantially omnidirectional receiving ability.
In various embodiments, a probe configured for 2D or 3D ping-based multiple aperture imaging may include multiple sizes of apparent point-source transmit transducers. As described above, larger apparent point-source transmit transducers may be beneficial for deep-tissue imaging, while smaller apparent point-source transmit transducers may be beneficial for relatively shallower-tissue imaging. In some embodiments, a probe may include apparent point-source transmit transducers of various sizes, each size being optimized for imaging at a particular depth. In some embodiments, images obtained at various depths may be combined to form a continuous image extending from relatively shallow regions to relatively deep regions.
In other embodiments, apparent point-source transducers may be used in combination with ultrasound probes sized to be inserted into a blood vessel or other bodily lumen within a patient (e.g., trans-esophageal probes, trans-urethral, trans-rectal probes, trans-vaginal probes).
Alternatively, as shown in
In some embodiments, apparent point-source transmit transducer may be provided on a separate probe from a receive transducer. Separate probes may have separate cables and may be entirely independent of one another. In other cases, a mechanical linkage may join the transmit and receive probes. When there is no mechanical linkage between separate transmit and receive transducers, the location of the transmit transducer can be determined by triangulation using receive timing from multiple receivers in the receive probe. Once the location of the transmitting apparent point-source is determined, the location of each of the scatterers in the tissue can be determined as described above using the intersection of ellipsoids generated via multiple receive elements.
In some embodiments, a transmit transducer 110 may be mounted to a distal end of an elongate catheter body 111 sized and configured to be positioned within a body cavity or lumen as shown for example in
Various methods may be used for triangulating the position of a transmit transducer relative to one or more receive transducers. For example, U.S. Pat. No. 8,473,239, which is incorporated herein by reference, provides a method for triangulating the position of a transducer element transmitting an ultrasound pulse based on the time-of-flight of a signal received by two or more receive elements located a known distance from one another. Similar techniques may be used to determine a location of one or more apparent point-source transmitters relative to one or more mechanically independent receive arrays. For example, an ultrasound ping transmitted from apparent point-source transmit transducer located on an intra-cavity probe positioned within a patient may be received with three or more elements of a receive array (or elements of separate receive arrays) positioned on the surface of the patient's skin at precisely known distances to one another.
When triangulating in three-dimensions, at least one of three of the receive elements should lie on a different plane than the remaining two. In other words, all receive elements should not lie in a common plane when triangulating the origin point of a three-dimensional ping.
The position of the apparent point-source transmitter may be triangulated based on the known position of the receive elements relative to one another and the difference between the time for the ping to reach the first receive element, the time for the ping to reach the second receive element, and the time for the ping to reach the third (and/or subsequent) receive element. In various embodiments, the time-of-flight values used for position measurements may be adjusted as described above in order to determine the location of the apparent point-source (i.e., the spherical center point) rather than the convex or concave transducer surface.
In some embodiments, using the position calculation method described above (or any other method), an imaging system with mechanically independent transmit and receive transducers may be configured to perform a transmitter-locating step while imaging. In some embodiments, the system may transmit a separate locating ping prior to transmitting one or more imaging pings. The position of the transmitter relative to the one or more receiver arrays may be determined and used in beamforming calculations for the echoes resulting from the one or more imaging pings following the locating ping. In some embodiments, locating pings may be transmitted at a different frequency than imaging pings.
In other embodiments, signals received directly from the imaging pings may be used to locate the transmitter relative to the one or more receiver arrays. In such embodiments, an initial peak of each transmitted imaging ping may be received at two or more receive elements at known positions relative to one another, and the position of the transmitter relative to the receive elements may be calculated as described above.
In some cases, all elements or arrays of a probe may be formed on a common substrate (e.g., using machining, molding, additive manufacturing, CMUT or other methods), while in other embodiments, apparent point-source and other arrays and/or elements may be made individually and assembled into a common probe structure using various other assembly methods.
Any 2D or 3D probe to be used for ping-based multiple aperture ultrasound imaging should be constructed and calibrated so as to record position information sufficient to define the actual acoustic position of each transducer element of the probe relative to a common coordinate system. Examples of systems and methods for aligning transducers during construction and for detecting the position of transducer elements are shown and described in Applicants' prior applications referenced above.
In some embodiments, apparent point-source transducers may be electrically connected within a probe so as to perform only transmit functions. In other embodiments, apparent point-source transducers may also be used to receive echoes by including a TX/RX switch.
Multiple Aperture Beamforming in 3D
Pings produced with apparent point-source transmit transducers may work particularly well with multiple aperture imaging, since echoes of the wide spherical-section 3D wavefront pulse may be received by multiple receivers over a wide area, much wider than any maximum coherence width of an imaged object with a non-homogeneous material. Beamforming during 3D ping-based multiple aperture imaging may involve calculating the locus of points corresponding to a time-of-flight of each received echo sample. In the case of 3D imaging, such a set of points is an ellipsoid—a three-dimensional elliptical surface having two focal points. One of the focal points is the position of the element on which the echo sample is received, and the other focal point is the position of the source from which the ping pulse was transmitted. In the case of apparent point-source transmitter transducer, the apparent transmit point used in such beamforming calculations is the center point of the spherical-cap shaped transducer forming the apparent point-source transmitter.
In some embodiments, ping-based multiple aperture imaging may operate by transmitting a spherical-section ultrasound ping signal from a first apparent point-source transmit transducer and receiving echoes with elements of two or more receive apertures. A complete “sub-image” may be formed by a controller or control system from echoes of a single ping received by a single receive element by triangulating the position of scatterers based on delay times between ping transmission and echo reception and the known three-dimensional positions of the transmit point and the receive point.
When using apparent point-source transducers, a time-of-transmit used in beamforming calculations may be substantially different from the time an electrical pulse is sent to the transducer. Because the sound waves do not actually originate from the apparent point-source (even though calculations are performed as if they do), the ping-transmit time used in beamforming calculations should be adjusted from the actual known value to an “apparent” value, corresponding to the time at which an outgoing waveform would appear to have been transmitted had it actually originated at the apparent point-source.
In the case of dome-shaped apparent point-source transducers, the transmitted ping start time may be adjusted by effectively adding an adjustment factor to each received echo time. In some embodiments, the adjustment time to be added may be a time value equal to the time required for a sound wave to travel from the spherical center to the convex outer transducer surface at a chosen speed of sound. In some embodiments, the chosen speed of sound for such an adjustment may be the same speed of sound used by an imaging system as the speed-of-sound in the object being imaged. In other embodiments, the chosen speed of sound for a convex apparent point-source adjustment may be the speed-of-sound in a stand-off or acoustic gel material immediately adjacent the transducer surface.
In the case of concave bowl-shaped apparent point-source transducers, the transmitted ping start time may be adjusted by effectively subtracting an adjustment factor from each received echo time. In some embodiments, the adjustment time to be subtracted may be a time value equal to the time required for a sound wave to travel from the inner concave transducer surface to the spherical center at a chosen speed of sound. In some embodiments, the chosen speed of sound for such an adjustment may be the same speed of sound used by an imaging system as the speed-of-sound in the object being imaged. In other embodiments, the chosen speed of sound for a concave apparent point-source adjustment may be the speed-of-sound in a stand-off or acoustic gel material immediately adjacent the transducer surface. In other embodiments, instead of subtracting time from each calculation, the adjustment may be calculated in terms of a number of data samples (e.g., based on a known sample rate), and an appropriate number of data samples at the beginning of a group of received echoes may be omitted from the beamforming calculation for each ping transmitted from a concave apparent point-source transmitter.
Images obtained from different unique combinations of one ping and one receive aperture may be referred to herein as “sub-image layers.” Multiple sub-image layers may be combined coherently to improve overall image quality. Additional image layer combining may be performed to further improve the quality of a final image. In the context of image layer combining, the term “image” may refer to a single two-dimensional pixel, a single voxel of a three-dimensional volume or a collection of any number of pixels or voxels.
Image layer combining may be described in terms of three image layer levels. These three cases include first-level image layers, second-level image layers and third-level image layers. (1) A first-level image layer may be formed from echoes received at a single receive aperture resulting from a single ping from a single transmit aperture (where a “transmit aperture” can be a single apparent point-source transmit element, a single small-element transmitter, or a group of transmit elements). For a unique combination of a single ping and a single receive aperture, the delayed echoes received by all the receive elements in the receive aperture may be summed coherently to obtain a first-level image layer. (2) Multiple first-level image layers resulting from echoes of multiple transmit pings (from the same or different transmit apertures) received at a single receive aperture can be summed together to produce a second-level image layer. Second-level image layers may be further processed to improve alignment or other image characteristics. (3) Third-level images may be obtained by combining second-level image layers formed with data from multiple receive apertures. In some embodiments, third-level images may be displayed as sequential time-domain frames to form a moving image.
In some embodiments, pixels or voxels of a first-level image layer may also be formed by summing in-phase and quadrature echo data, that is by summing each echo with an echo ¼ wavelength delayed for each receive-aperture element. In some cases, echo data may be sampled and stored as an in-phase data set and as a separate quadrature data set. In other cases, if the digital sampling rate is divisible by four, then a quadrature sample corresponding to an in-phase sample may be identified by selecting a sample at an appropriate number of samples prior to the in-phase sample. If the desired quadrature sample does not correspond to an existing whole sample, a quadrature sample may be obtained by interpolation. Combining in-phase and quadrature data for a single image (pixel, voxel or collection of pixels or voxels) may provide the advantage of increasing the resolution of the echo data without introducing blurring effects. Similarly, samples at values other than ¼ wavelength may be combined with in-phase samples in order to improve various imaging characteristics.
Combination, summation or averaging of various image layers may be accomplished either by coherent addition, incoherent addition, or a combination of the two. Coherent addition (incorporating both phase and magnitude information during image layer summation) tends to maximize lateral resolution, whereas incoherent addition (summing magnitudes only and omitting phase information) tends to average out speckle noise and minimize the effects of image layer alignment errors that may be caused by minor variations in the speed of sound through the imaged medium. Speckle noise is reduced through incoherent summing because each image layer will tend to develop its own independent speckle pattern, and summing the patterns incoherently has the effect of averaging out these speckle patterns. Alternatively, if the patterns are added coherently, they reinforce each other and only one strong speckle pattern results. Incoherent addition may be thought of as akin to instantaneous compound imaging, which has long been known as a means to suppress speckle noise.
In most embodiments, echoes received by elements of a single receive aperture are typically combined coherently. In some embodiments, the number of receive apertures and/or the size of each receive aperture may be changed in order to maximize some desired combination of image quality metrics such as lateral resolution, speed-of-sound variation tolerance, speckle noise reduction, etc. In some embodiments, such alternative element-to-aperture grouping arrangements may be selectable by a user. In other embodiments, such arrangements may be automatically selected or developed by an imaging system.
Variations in the speed of sound may be tolerated by incoherent addition as follows: Summing two pixels coherently with a speed-of-sound variation resulting in only half a wavelength's delay (e.g., approximately 0.25 mm for a 3 MHz probe) results in destructive phase cancellation, which causes significant image data loss; if the pixels are added incoherently, the same or even greater delay causes only an insignificant spatial distortion in the image layer and no loss of image data. The incoherent addition of such image layers may result in some smoothing of the final image (in some embodiments, such smoothing may be added intentionally to make the image more readable).
At all three image layer levels, coherent addition can lead to maximum lateral resolution of a multiple aperture system if the geometry of the probe elements is known to a desired degree of precision and the assumption of a constant speed of sound across all paths is valid. Likewise, at all image layer levels, incoherent addition leads to the best averaging out of speckle noise and tolerance of minor variations in speed of sound through the imaged medium.
In some embodiments, coherent addition can be used to combine image layers resulting from apertures for which phase cancellation is not likely to be a problem, and incoherent addition can then be used where phase cancellation would be more likely to present a problem, such as when combining images formed from echoes received at different receive apertures separated by a distance exceeding some threshold.
In some embodiments, all first-level images may be formed by using coherent addition of all sub-image layers obtained from elements of a common receive aperture, assuming the receive aperture has a width less than the maximum coherent aperture width. For second and third level image layers, many combinations of coherent and incoherent summation are possible. For example, in some embodiments, second-level image layers may be formed by coherently summing contributing first-level image layers, while third-level image layers may be formed by incoherent summing of the contributing second-level image layers.
Time-domain frames may be formed by any level of image layer depending on the desired trade-off between processing time and image quality. Higher-level images will tend to be of higher quality, but may also require more processing time. Thus, if it is desired to provide real-time imaging, the level of image layer combination processing may be limited in order to display images without significant “lag” being visible to the operator. The details of such a trade-off will depend on the particular processing hardware in use as well as other factors.
An example may be understood with reference to
In use, a 3D ping in the form of a spherical-section wavefront (e.g., a perfectly hemispherical wavefront, or a wavefront with a shape greater or less than a hemisphere) may be transmitted from apparent point-source transmit transducer (e.g., T1). The wavefront may travel into the object 96 and may be reflected by a reflector such as that shown at point A. Reflected ultrasound signals may then be received by the elements of a receive aperture (e.g., elements R2.1 through R2.n of receive aperture R2). Three-dimensional echo data may be collected by receiving echoes with receive apertures that do not lie on a common plane with each other and/or with the transmit aperture. For example, each of receive apertures R1 and R2 comprise elements that are spaced from one another in at least two dimensions relative to the imaged object 96 and thus the receive apertures R1 and R2 do not lie on the same plane as one another, nor on a common plane with the transmit aperture T1. Thus, in some embodiments, echoes may be received by all receive elements, and those elements may be grouped into apertures based on factors such as position and maximum coherence width.
The echo sample corresponding to point A may then be beamformed to determine the possible location of point A within the object 96. Beamforming may proceed by calculating a locus of possible positions for point A based on the known transmit time, the known echo receive time, the known location of a receive element (e.g., element R2.1 of receive aperture R2) and the known location of the apparent point-source (i.e., the spherical center of transmit element T1 in this example). The same calculation may be performed for echoes received at the remaining receive elements of the same receive aperture, each receive element defining a slightly different ellipsoid. The information from all of the receive elements may then be combined (e.g., coherently as described above) in order to converge on a small three-dimensional estimate of the location of point A. The same process may then be repeated for echo samples of point A received by elements of a second (or more) receive aperture from the same transmitted ping. Similarly, the same process may be repeated for a second (or more) 3D pings transmitted from the same or different transmit elements.
Because a high-quality 3D image may be obtained from echoes of a single transmitted ping using the above procedure, a ping-based multiple aperture imaging system with apparent point-source transmitters may be used for performing 4D imaging to display a real-time 3D video of a moving object.
Similarly, 3D pings transmitted from apparent point-source transmitter may also be used for ping-based multiple aperture Doppler, multiple aperture elastography or any other imaging technique that may benefit from high-frame-rate 3D imaging.
2D Imaging while Collecting 3D Data
In some embodiments, a form of 2D imaging may be performed using a probe configured for 3D imaging by simply beamforming and displaying only a 2D slice of received echo data from the received three-dimensional echo data. For example, such techniques may be used in order to reduce a beamform calculation and simplify display for real-time imaging using an imaging device with limited processing capability, while still retaining the full 3D echo data. The full 3D echo data may be beamformed and reviewed at a later time using a device with greater processing power. In some embodiments, the 2D slice to be beamformed and displayed may be automatically selected by an imaging device. Alternatively, the 2D slice to be beamformed and displayed may be selected or adjusted by an operator of the device.
Apparent Source for 2D Planar Imaging
In some embodiments, apparent source transmit transducer may be configured specifically for 2D planar imaging. While apparent point-source transmitter as described above will transmit ultrasound in three dimensions, a 2D apparent source transducer may be configured to transmit ultrasound signals that are confined to a single plane (or at least with minimal “leakage” out of the image plane). In some embodiments, such a 2D apparent source transducer 150 may have a shape of a shell 152 with a cylindrical section shape such as that shown in
As with the spherical-section transducers described above, cylindrical section transducers 150 may be made in various shapes, sizes, and circular cross-sections as needed (the description of
A cylindrical-section transducer 150 may also be configured such that ultrasound signals are transmitted into an imaged medium from either the convex surface or the concave surface. In various embodiments, a probe including a cylindrical section transducer 150 may include a focusing lens configured to focus transmitted ultrasound energy in the imaging plane. In use, time adjustments may be made to treat the cylindrical centerline (i.e., the circular center of the sphere) as the mathematical origin of pings transmitted from such a transducer.
In various embodiments, receive elements may also be formed from cylindrical section shell structures. Receive elements may typically have substantially smaller transducer surface areas, since echoes received by many receive transducers may be combined into to increase the received echo signal.
Multiple Aperture Ultrasound Imaging System Components
The block diagram of
The transmission of ultrasound signals from elements of the probe 202 may be controlled by a transmit control subsystem 204. In some embodiments, the transmit control subsystem 204 may include any combination of analog and digital components for controlling transducer elements of the probe 202 to transmit un-focused ultrasound pings at desired frequencies and intervals from selected transmit apertures according to a desired imaging algorithm. In some embodiments a transmit control system 204 may be configured to transmit ultrasound pings at a range of ultrasound frequencies. In some (though not all) embodiments, the transmit control subsystem may also be configured to control the probe in a phased array mode, transmitting focused (i.e., transmit beamformed) ultrasound scanline beams.
In some embodiments, a transmit control sub-system 204 may include a transmit signal definition module 206 and a transmit element control module 208. The transmit signal definition module 206 may include suitable combinations of hardware, firmware and/or software configured to define desired characteristics of a signal to be transmitted by an ultrasound probe. For example, the transmit signal definition module 206 may establish (e.g., based on user inputs or on predetermined factors) characteristics of an ultrasound signal to be transmitted such as a pulse start time, pulse length (duration), ultrasound frequency, pulse power, pulse shape, pulse direction (if any), pulse amplitude, transmit aperture location, or any other characteristics.
The transmit element control module 208 may then take information about the desired transmit pulse and determine the corresponding electrical signals to be sent to the appropriate transducer elements in order to produce this signal. In various embodiments, the signal definition module 206 and the transmit element control module 208 may comprise separate electronic components, or may include portions of one or more common components.
Upon receiving echoes of transmitted signals from a region of interest, the probe elements may generate time-varying electrical signals corresponding to the received ultrasound vibrations. Signals representing the received echoes may be output from the probe 202 and sent to a receive subsystem 210. In some embodiments, the receive subsystem may include multiple channels, each of which may include an analog front-end device (“AFE”) 212 and an analog-to-digital conversion device (ADC) 214. In some embodiments, each channel of the receive subsystem 210 may also include digital filters and data conditioners (not shown) after the ADC 214. In some embodiments, analog filters prior to the ADC 214 may also be provided. The output of each ADC 214 may be directed into a raw data memory device 220. In some embodiments, an independent channel of the receive subsystem 210 may be provided for each receive transducer element of the probe 202. In other embodiments, two or more transducer elements may share a common receive channel.
In some embodiments, an analog front-end device 212 (AFE) may perform certain filtering processes before passing the signal to an analog-to-digital conversion device 214 (ADC). The ADC 214 may be configured to convert received analog signals into a series of digital data points at some pre-determined sampling rate. Unlike most ultrasound systems, some embodiments of the ultrasound imaging system of
In order to convert the captured digital samples into an image, the data may be retrieved from the raw data memory 220 by an image generation subsystem 230. As shown, the image generation subsystem 230 may include a beamforming block 232 and an image layer combining (“ILC”) block 234. In some embodiments, a beamformer 232 may be in communication with a calibration memory 238 that contains probe calibration data. Probe calibration data may include information about the precise position, operational quality, and/or other information about individual probe transducer elements. The calibration memory 238 may be physically located within the probe, within the imaging system, or in location external to both the probe and the imaging system.
In some embodiments, after passing through the image generation block 230, image data may then be stored in an image buffer memory 236 which may store beamformed and (in some embodiments) layer-combined image frames. A video processor 242 within a video subsystem 240 may then retrieve image frames from the image buffer, and may process the images into a video stream that may be displayed on a video display 244 and/or stored in a video memory 246 as a digital video clip, e.g., as referred to in the art as a “cine loop”.
In some embodiments, the AFE 212 may be configured to perform various amplification and filtering processes to a received analog signal before passing the analog signal to an analog-to-digital conversion device. For example, an AFE 212 may include amplifiers such as a low noise amplifier (LNA), a variable gain amplifier (VGA), a bandpass or lowpass/anti-aliasing filter, and/or other amplification or filtering devices. In some embodiments, an AFE device 212 may be configured to begin passing an analog signal to an ADC 214 upon receiving a trigger signal. In other embodiments, an AFE device can be “free running”, continuously passing an analog signal to an ADC.
In some embodiments, each analog-to-digital converter 214 may generally include any device configured to sample a received analog signal at some consistent, predetermined sampling rate. For example, in some embodiments, an analog-to-digital converter may be configured to record digital samples of a time-varying analog signal at 25 MHz, which is 25 million samples per second or one sample every 40 nanoseconds. Thus, data sampled by an ADC may simply include a list of data points, each of which may correspond to a signal value at a particular instant. In some embodiments, an ADC 214 may be configured to begin digitally sampling an analog signal upon receiving a trigger signal. In other embodiments, an ADC device can be “free running”, continuously sampling a received analog signal.
In some embodiments, the raw data memory device 220 may include any suitable volatile or non-volatile digital memory storage device. In some embodiments, the raw data memory 220 may also comprise communication electronics for transmitting raw digital ultrasound data to an external device over a wired or wireless network. In such cases, the transmitted raw echo data may be stored on the external device in any desired format. In other embodiments, the raw data memory 220 may include a combination of volatile memory, non-volatile memory and communication electronics.
In some embodiments, the raw data memory device 220 may comprise a temporary (volatile or non-volatile) memory section, and a long-term non-volatile memory section. In an example of such embodiments, the temporary memory may act as a buffer between the ADC 214 and the beamformer 232 in cases where the beamformer 232 may be unable to operate fast enough to accommodate data at the full rate from the ADC 214. In some embodiments, a long-term non-volatile memory device may be configured to receive data from a temporary memory device or directly from the ADC 214. Such a long-term memory device may be configured to store a quantity of raw echo data for subsequent processing, analysis or transmission to an external device.
In some embodiments, the beamforming block 232 and the image layer combining block 234 may each include any digital signal processing and/or computing components configured to perform the specified processes (e.g., as described below). For example, in various embodiments the beamforming 232 and image layer combining 234 may be performed by software running on a single GPU, on multiple GPUs, on one or more CPUs, on combinations of CPUs & GPUs, on single or multiple accelerator cards or modules, on a distributed processing system, or a clustered processing system. Alternatively, these or other processes may be performed by firmware running on an FPGA (Field Programmable Gate Array) architecture or one or more dedicated ASIC (Application-Specific Integrated Circuit) devices.
In some embodiments, the video processor 242 may include any video processing hardware, firmware and software components that may be configured to assemble image frames into a video stream for display and/or storage.
Data Capture & Offline Analysis
In various embodiments, raw un-beamformed echo data resulting from a ping transmitted from apparent point-source transmit transducer and received by one or more arrays of receive transducer elements may be captured and stored in a raw data memory device for subsequent retrieval and analysis. In addition to such echo data, additional data may also be stored and/or transmitted over a network and retrieved for subsequent image generation and analysis. Such additional data may include calibration data describing the positions of the transmitting and receiving transducer elements, and transmit data describing the identity (or position) of transmitting transducers associated with specific echo data.
After retrieving such data, a clinician may use the data to reconstruct imaging sessions in a variety of ways while making adjustments that may not have been made during a live imaging session. For example, images of a series of 2D slices through the 3D volume may be generated and shown in succession in order to simulate a 2D transducer passing across a surface of the region of interest. Examples of such methods are described in Applicant's co-pending U.S. patent application Ser. No. 13/971,689 filed Aug. 20, 2013, which is incorporated herein by reference.
Some embodiments of a probe configured for imaging an entire patient's body or a substantial portion of a patient's body may comprise an array of apparent point-source transmitters and receive elements sized to cover a substantial portion of the desired region of a patient's body. For example, a probe may be sized to cover substantially half of a patient's chest or more. Such a probe may have a maximum dimension of about 8 cm to about 10 cm.
Alternatively, a much smaller probe capable of insonifying a conically-shaped volume of, for example, + or −30 degrees, can be placed on a patient's body such that an organ of interest may be included in the cone. Such a probe may be placed in more than one place to cover a larger volume of interest.
Although this invention has been disclosed in the context of certain preferred embodiments and examples, it will be understood by those skilled in the art that the present invention extends beyond the specifically disclosed embodiments to other alternative embodiments and/or uses of the invention and obvious modifications and equivalents thereof. Various modifications to the above embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, it is intended that the scope of the present invention herein disclosed should not be limited by the particular disclosed embodiments described above, but should be determined only by a fair reading of the claims that follow.
In particular, materials and manufacturing techniques may be employed as within the level of those with skill in the relevant art. Furthermore, reference to a singular item, includes the possibility that there are plural of the same items present. More specifically, as used herein and in the appended claims, the singular forms “a,” “and,” “said,” and “the” include plural referents unless the context clearly dictates otherwise. Also as used herein, unless explicitly stated otherwise, the term “or” is inclusive of all presented alternatives, and means essentially the same as the commonly used phrase “and/or.” It is further noted that the claims may be drafted to exclude any optional element. As such, this statement is intended to serve as antecedent basis for use of such exclusive terminology as “solely,” “only” and the like in connection with the recitation of claim elements, or use of a “negative” limitation. Unless defined otherwise herein, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
This application is a continuation of U.S. patent application Ser. No. 15/888,738, filed Feb. 5, 2018, which is a divisional of U.S. patent application Ser. No. 14/279,052, filed May 15, 2014, now U.S. Pat. No. 9,883,848, which application claims the benefit of U.S. Provisional Patent Application No. 61/877,555, filed Sep. 13, 2013, titled “Ultrasound Imaging Using Virtual Point-Source Transmission”, herein incorporated by reference in their entirety.
| Number | Name | Date | Kind |
|---|---|---|---|
| 3174286 | Erickson | Mar 1965 | A |
| 3895381 | Kock | Jul 1975 | A |
| 3974692 | Hassler | Aug 1976 | A |
| 4055988 | Dutton | Nov 1977 | A |
| 4072922 | Taner et al. | Feb 1978 | A |
| 4097835 | Green | Jun 1978 | A |
| 4105018 | Greenleaf et al. | Aug 1978 | A |
| 4180792 | Lederman et al. | Dec 1979 | A |
| 4205394 | Pickens | May 1980 | A |
| 4229798 | Rosie | Oct 1980 | A |
| 4259733 | Taner et al. | Mar 1981 | A |
| 4265126 | Papadofrangakis et al. | May 1981 | A |
| 4271842 | Specht et al. | Jun 1981 | A |
| 4325257 | Kino et al. | Apr 1982 | A |
| 4327738 | Green et al. | May 1982 | A |
| 4333474 | Nigam | Jun 1982 | A |
| 4339952 | Foster | Jul 1982 | A |
| 4452084 | Taenzer | Jun 1984 | A |
| 4501279 | Seo | Feb 1985 | A |
| 4511998 | Kanda et al. | Apr 1985 | A |
| 4539847 | Paap | Sep 1985 | A |
| 4566459 | Umemura et al. | Jan 1986 | A |
| 4567768 | Satoh et al. | Feb 1986 | A |
| 4604697 | Luthra et al. | Aug 1986 | A |
| 4662222 | Johnson | May 1987 | A |
| 4669482 | Ophir | Jun 1987 | A |
| 4682497 | Sasaki | Jul 1987 | A |
| 4694434 | Vonn Ramm et al. | Sep 1987 | A |
| 4781199 | Hirama et al. | Nov 1988 | A |
| 4817434 | Anderson | Apr 1989 | A |
| 4831601 | Breimesser et al. | May 1989 | A |
| 4893284 | Magrane | Jan 1990 | A |
| 4893628 | Angelsen | Jan 1990 | A |
| 4990462 | Sliwa, Jr. | Feb 1991 | A |
| 5050588 | Grey et al. | Sep 1991 | A |
| 5062295 | Shakkottai et al. | Nov 1991 | A |
| 5141738 | Rasor et al. | Aug 1992 | A |
| 5161536 | Vilkomerson et al. | Nov 1992 | A |
| 5197475 | Antich et al. | Mar 1993 | A |
| 5226019 | Bahorich | Jul 1993 | A |
| 5230339 | Charlebois | Jul 1993 | A |
| 5269309 | Fort et al. | Dec 1993 | A |
| 5278757 | Hoctor et al. | Jan 1994 | A |
| 5293871 | Reinstein et al. | Mar 1994 | A |
| 5299576 | Shiba | Apr 1994 | A |
| 5301674 | Erikson et al. | Apr 1994 | A |
| 5305756 | Entrekin et al. | Apr 1994 | A |
| 5339282 | Kuhn et al. | Aug 1994 | A |
| 5340510 | Bowen | Aug 1994 | A |
| 5345426 | Lipschutz | Sep 1994 | A |
| 5349960 | Gondo | Sep 1994 | A |
| 5355888 | Kendall | Oct 1994 | A |
| 5381794 | Tei et al. | Jan 1995 | A |
| 5398216 | Hall et al. | Mar 1995 | A |
| 5409010 | Beach et al. | Apr 1995 | A |
| 5442462 | Guissin | Aug 1995 | A |
| 5454372 | Banjanin et al. | Oct 1995 | A |
| 5503152 | Oakley et al. | Apr 1996 | A |
| 5515853 | Smith et al. | May 1996 | A |
| 5515856 | Olstad et al. | May 1996 | A |
| 5522393 | Phillips et al. | Jun 1996 | A |
| 5526815 | Granz et al. | Jun 1996 | A |
| 5544659 | Banjanin | Aug 1996 | A |
| 5558092 | Unger et al. | Sep 1996 | A |
| 5564423 | Mele et al. | Oct 1996 | A |
| 5568812 | Murashita et al. | Oct 1996 | A |
| 5570691 | Wright et al. | Nov 1996 | A |
| 5581517 | Gee et al. | Dec 1996 | A |
| 5625149 | Gururaja et al. | Apr 1997 | A |
| 5628320 | Teo | May 1997 | A |
| 5673697 | Bryan et al. | Oct 1997 | A |
| 5675550 | Ekhaus | Oct 1997 | A |
| 5720291 | Schwartz | Feb 1998 | A |
| 5720708 | Lu et al. | Feb 1998 | A |
| 5744898 | Smith et al. | Apr 1998 | A |
| 5769079 | Hossack | Jun 1998 | A |
| 5784334 | Sena et al. | Jul 1998 | A |
| 5785654 | Iinuma et al. | Jul 1998 | A |
| 5795297 | Daigle | Aug 1998 | A |
| 5797845 | Barabash et al. | Aug 1998 | A |
| 5798459 | Ohba et al. | Aug 1998 | A |
| 5820561 | Olstad et al. | Oct 1998 | A |
| 5838564 | Bahorich et al. | Nov 1998 | A |
| 5850622 | Vassiliou et al. | Dec 1998 | A |
| 5862100 | VerWest | Jan 1999 | A |
| 5870691 | Partyka et al. | Feb 1999 | A |
| 5876342 | Chen et al. | Mar 1999 | A |
| 5891038 | Seyed-Bolorforosh et al. | Apr 1999 | A |
| 5892732 | Gersztenkorn | Apr 1999 | A |
| 5916169 | Hanafy et al. | Jun 1999 | A |
| 5919139 | Lin | Jul 1999 | A |
| 5920285 | Benjamin | Jul 1999 | A |
| 5930730 | Marfurt et al. | Jul 1999 | A |
| 5940778 | Marfurt et al. | Aug 1999 | A |
| 5951479 | Holm et al. | Sep 1999 | A |
| 5964707 | Fenster et al. | Oct 1999 | A |
| 5969661 | Benjamin | Oct 1999 | A |
| 5999836 | Nelson et al. | Dec 1999 | A |
| 6007499 | Martin et al. | Dec 1999 | A |
| 6013032 | Savord | Jan 2000 | A |
| 6014473 | Hossack et al. | Jan 2000 | A |
| 6048315 | Chiao et al. | Apr 2000 | A |
| 6049509 | Sonneland et al. | Apr 2000 | A |
| 6050943 | Slayton et al. | Apr 2000 | A |
| 6056693 | Haider | May 2000 | A |
| 6058074 | Swan et al. | May 2000 | A |
| 6077224 | Lang et al. | Jun 2000 | A |
| 6092026 | Bahorich et al. | Jul 2000 | A |
| 6122538 | Sliwa, Jr. et al. | Sep 2000 | A |
| 6123670 | Mo | Sep 2000 | A |
| 6129672 | Seward et al. | Oct 2000 | A |
| 6135960 | Holmberg | Oct 2000 | A |
| 6138075 | Yost | Oct 2000 | A |
| 6148095 | Prause et al. | Nov 2000 | A |
| 6162175 | Marian, Jr. et al. | Dec 2000 | A |
| 6166384 | Dentinger et al. | Dec 2000 | A |
| 6166853 | Sapia et al. | Dec 2000 | A |
| 6193665 | Hall et al. | Feb 2001 | B1 |
| 6196739 | Silverbrook | Mar 2001 | B1 |
| 6200266 | Shokrollahi et al. | Mar 2001 | B1 |
| 6210335 | Miller | Apr 2001 | B1 |
| 6213958 | Winder | Apr 2001 | B1 |
| 6221019 | Kantorovich | Apr 2001 | B1 |
| 6222304 | Bernstein | Apr 2001 | B1 |
| 6231511 | Bae | May 2001 | B1 |
| 6238342 | Feleppa et al. | May 2001 | B1 |
| 6246901 | Benaron | Jun 2001 | B1 |
| 6251073 | Imran et al. | Jun 2001 | B1 |
| 6264609 | Herrington et al. | Jul 2001 | B1 |
| 6266551 | Osadchy et al. | Jul 2001 | B1 |
| 6278949 | Alam | Aug 2001 | B1 |
| 6289230 | Chaiken et al. | Sep 2001 | B1 |
| 6299580 | Asafusa | Oct 2001 | B1 |
| 6304684 | Niczyporuk et al. | Oct 2001 | B1 |
| 6309356 | Ustuner et al. | Oct 2001 | B1 |
| 6324453 | Breed et al. | Nov 2001 | B1 |
| 6345539 | Rawes et al. | Feb 2002 | B1 |
| 6361500 | Masters | Mar 2002 | B1 |
| 6363033 | Cole et al. | Mar 2002 | B1 |
| 6370480 | Gupta et al. | Apr 2002 | B1 |
| 6374185 | Taner et al. | Apr 2002 | B1 |
| 6394955 | Perlitz | May 2002 | B1 |
| 6423002 | Hossack | Jul 2002 | B1 |
| 6436046 | Napolitano et al. | Aug 2002 | B1 |
| 6449821 | Sudol et al. | Sep 2002 | B1 |
| 6450965 | Williams et al. | Sep 2002 | B2 |
| 6468216 | Powers et al. | Oct 2002 | B1 |
| 6471650 | Powers et al. | Oct 2002 | B2 |
| 6475150 | Haddad | Nov 2002 | B2 |
| 6480790 | Calvert et al. | Nov 2002 | B1 |
| 6487502 | Taner | Nov 2002 | B1 |
| 6499536 | Ellingsen | Dec 2002 | B1 |
| 6508768 | Hall et al. | Jan 2003 | B1 |
| 6508770 | Cai | Jan 2003 | B1 |
| 6517484 | Wilk et al. | Feb 2003 | B1 |
| 6526163 | Halmann et al. | Feb 2003 | B1 |
| 6543272 | Vitek | Apr 2003 | B1 |
| 6547732 | Jago | Apr 2003 | B2 |
| 6551246 | Ustuner et al. | Apr 2003 | B1 |
| 6565510 | Haider | May 2003 | B1 |
| 6585647 | Winder | Jul 2003 | B1 |
| 6597171 | Hurlimann et al. | Jul 2003 | B2 |
| 6604421 | Li | Aug 2003 | B1 |
| 6614560 | Silverbrook | Sep 2003 | B1 |
| 6620101 | Azzam et al. | Sep 2003 | B2 |
| 6652461 | Levkovitz | Nov 2003 | B1 |
| 6668654 | Dubois et al. | Dec 2003 | B2 |
| 6672165 | Rather et al. | Jan 2004 | B2 |
| 6681185 | Young et al. | Jan 2004 | B1 |
| 6690816 | Aylward et al. | Feb 2004 | B2 |
| 6692450 | Coleman | Feb 2004 | B1 |
| 6695778 | Golland et al. | Feb 2004 | B2 |
| 6702745 | Smythe | Mar 2004 | B1 |
| 6704692 | Banerjee et al. | Mar 2004 | B1 |
| 6719693 | Richard | Apr 2004 | B2 |
| 6728567 | Rather et al. | Apr 2004 | B2 |
| 6752762 | DeJong et al. | Jun 2004 | B1 |
| 6755787 | Hossack et al. | Jun 2004 | B2 |
| 6780152 | Ustuner et al. | Aug 2004 | B2 |
| 6790182 | Eck et al. | Sep 2004 | B2 |
| 6835178 | Wilson et al. | Dec 2004 | B1 |
| 6837853 | Marian | Jan 2005 | B2 |
| 6843770 | Sumanaweera | Jan 2005 | B2 |
| 6847737 | Kouri et al. | Jan 2005 | B1 |
| 6854332 | Alleyne | Feb 2005 | B2 |
| 6865140 | Thomenius et al. | Mar 2005 | B2 |
| 6932767 | Landry et al. | Aug 2005 | B2 |
| 7033320 | Von Behren et al. | Apr 2006 | B2 |
| 7087023 | Daft et al. | Aug 2006 | B2 |
| 7104956 | Christopher | Sep 2006 | B1 |
| 7217243 | Takeuchi | May 2007 | B2 |
| 7221867 | Silverbrook | May 2007 | B2 |
| 7231072 | Yamano et al. | Jun 2007 | B2 |
| 7269299 | Schroeder | Sep 2007 | B2 |
| 7283652 | Mendonca et al. | Oct 2007 | B2 |
| 7285094 | Nohara et al. | Oct 2007 | B2 |
| 7293462 | Lee et al. | Nov 2007 | B2 |
| 7313053 | Wodnicki | Dec 2007 | B2 |
| 7366704 | Reading et al. | Apr 2008 | B2 |
| 7402136 | Hossack et al. | Jul 2008 | B2 |
| 7410469 | Talish et al. | Aug 2008 | B1 |
| 7415880 | Renzel | Aug 2008 | B2 |
| 7443765 | Thomenius et al. | Oct 2008 | B2 |
| 7444875 | Wu et al. | Nov 2008 | B1 |
| 7447535 | Lavi | Nov 2008 | B2 |
| 7448998 | Robinson | Nov 2008 | B2 |
| 7466848 | Metaxas et al. | Dec 2008 | B2 |
| 7469096 | Silverbrook | Dec 2008 | B2 |
| 7474778 | Shinomura et al. | Jan 2009 | B2 |
| 7481577 | Ramamurthy et al. | Jan 2009 | B2 |
| 7491171 | Barthe et al. | Feb 2009 | B2 |
| 7497828 | Wilk et al. | Mar 2009 | B1 |
| 7497830 | Li | Mar 2009 | B2 |
| 7510529 | Chou et al. | Mar 2009 | B2 |
| 7514851 | Wilser et al. | Apr 2009 | B2 |
| 7549962 | Dreschel et al. | Jun 2009 | B2 |
| 7574026 | Rasche et al. | Aug 2009 | B2 |
| 7625343 | Cao et al. | Dec 2009 | B2 |
| 7637869 | Sudol | Dec 2009 | B2 |
| 7668583 | Fegert et al. | Feb 2010 | B2 |
| 7674228 | Williams et al. | Mar 2010 | B2 |
| 7682311 | Simopoulos et al. | Mar 2010 | B2 |
| 7699776 | Walker et al. | Apr 2010 | B2 |
| 7722541 | Cai | May 2010 | B2 |
| 7744532 | Ustuner et al. | Jun 2010 | B2 |
| 7750311 | Daghighian | Jul 2010 | B2 |
| 7764984 | Desmedt et al. | Jul 2010 | B2 |
| 7785260 | Umemura et al. | Aug 2010 | B2 |
| 7787680 | Ahn et al. | Aug 2010 | B2 |
| 7806828 | Stringer | Oct 2010 | B2 |
| 7819810 | Stringer et al. | Oct 2010 | B2 |
| 7822250 | Yao et al. | Oct 2010 | B2 |
| 7824337 | Abe et al. | Nov 2010 | B2 |
| 7833163 | Cai | Nov 2010 | B2 |
| 7837624 | Hossack et al. | Nov 2010 | B1 |
| 7846097 | Jones et al. | Dec 2010 | B2 |
| 7850613 | Stribling | Dec 2010 | B2 |
| 7862508 | Davies et al. | Jan 2011 | B2 |
| 7876945 | Lötjönen | Jan 2011 | B2 |
| 7880154 | Otto | Feb 2011 | B2 |
| 7887486 | Ustuner et al. | Feb 2011 | B2 |
| 7901358 | Mehi et al. | Mar 2011 | B2 |
| 7914451 | Davies | Mar 2011 | B2 |
| 7919906 | Cerofolini | Apr 2011 | B2 |
| 7926350 | Kröning et al. | Apr 2011 | B2 |
| 7927280 | Davidsen | Apr 2011 | B2 |
| 7972271 | Johnson et al. | Jul 2011 | B2 |
| 7984637 | Ao et al. | Jul 2011 | B2 |
| 7984651 | Randall et al. | Jul 2011 | B2 |
| 8002705 | Napolitano et al. | Aug 2011 | B1 |
| 8007439 | Specht | Aug 2011 | B2 |
| 8057392 | Hossack et al. | Nov 2011 | B2 |
| 8057393 | Yao et al. | Nov 2011 | B2 |
| 8079263 | Randall et al. | Dec 2011 | B2 |
| 8079956 | Azuma et al. | Dec 2011 | B2 |
| 8088067 | Vortman et al. | Jan 2012 | B2 |
| 8088068 | Yao et al. | Jan 2012 | B2 |
| 8088071 | Hwang et al. | Jan 2012 | B2 |
| 8105239 | Specht | Jan 2012 | B2 |
| 8135190 | Bae et al. | Mar 2012 | B2 |
| 8157737 | Zhang et al. | Apr 2012 | B2 |
| 8182427 | Wu et al. | May 2012 | B2 |
| 8202219 | Luo et al. | Jun 2012 | B2 |
| 8265175 | Barsoum et al. | Sep 2012 | B2 |
| 8277383 | Specht | Oct 2012 | B2 |
| 8279705 | Choi et al. | Oct 2012 | B2 |
| 8412307 | Willis et al. | Apr 2013 | B2 |
| 8414564 | Goldshleger et al. | Apr 2013 | B2 |
| 8419642 | Sandrin et al. | Apr 2013 | B2 |
| 8473239 | Specht et al. | Jun 2013 | B2 |
| 8478382 | Burnside et al. | Jul 2013 | B2 |
| 8483804 | Hsieh et al. | Jul 2013 | B2 |
| 8532951 | Roy et al. | Sep 2013 | B2 |
| 8582848 | Funka-Lea et al. | Nov 2013 | B2 |
| 8602993 | Specht et al. | Dec 2013 | B2 |
| 8627724 | Papadopoulos et al. | Jan 2014 | B2 |
| 8634615 | Brabec | Jan 2014 | B2 |
| 8672846 | Napolitano et al. | Mar 2014 | B2 |
| 8684936 | Specht | Apr 2014 | B2 |
| 9036887 | Fouras et al. | May 2015 | B2 |
| 9072495 | Specht | Jul 2015 | B2 |
| 9146313 | Specht et al. | Sep 2015 | B2 |
| 9152761 | Bhatia et al. | Oct 2015 | B2 |
| 9176078 | Flohr et al. | Nov 2015 | B2 |
| 9192355 | Specht et al. | Nov 2015 | B2 |
| 9217660 | Zlotnick et al. | Dec 2015 | B2 |
| 9220478 | Smith et al. | Dec 2015 | B2 |
| 9247874 | Kumar et al. | Feb 2016 | B2 |
| 9247926 | Smith et al. | Feb 2016 | B2 |
| 9265484 | Brewer et al. | Feb 2016 | B2 |
| 9268777 | Lu et al. | Feb 2016 | B2 |
| 9271661 | Moghari et al. | Mar 2016 | B2 |
| 9277861 | Kowal et al. | Mar 2016 | B2 |
| 9282945 | Smith et al. | Mar 2016 | B2 |
| 9339239 | Wang et al. | May 2016 | B2 |
| 9339256 | Specht et al. | May 2016 | B2 |
| 9392986 | Ning et al. | Jul 2016 | B2 |
| 9420994 | Specht | Aug 2016 | B2 |
| 9510806 | Smith et al. | Dec 2016 | B2 |
| 9526475 | Specht et al. | Dec 2016 | B2 |
| 9526485 | Yang | Dec 2016 | B2 |
| 9572549 | Belevich et al. | Feb 2017 | B2 |
| 9576354 | Fouras et al. | Feb 2017 | B2 |
| 9582876 | Specht | Feb 2017 | B2 |
| 9606206 | Boernert et al. | Mar 2017 | B2 |
| 9659152 | Mueller | May 2017 | B2 |
| 9668714 | Call et al. | Jun 2017 | B2 |
| 9775511 | Kumar et al. | Oct 2017 | B2 |
| 9788813 | Adam et al. | Oct 2017 | B2 |
| 9883848 | Specht et al. | Feb 2018 | B2 |
| 9901407 | Breisacher et al. | Feb 2018 | B2 |
| 9986975 | Specht et al. | Jun 2018 | B2 |
| 10064605 | Belevich et al. | Sep 2018 | B2 |
| 10130333 | Specht | Nov 2018 | B2 |
| 10206662 | Smith et al. | Feb 2019 | B2 |
| 10226234 | Specht et al. | Mar 2019 | B2 |
| 10267913 | Smith et al. | Apr 2019 | B2 |
| 10342518 | Specht et al. | Jul 2019 | B2 |
| 10380399 | Call et al. | Aug 2019 | B2 |
| 10401493 | Call et al. | Sep 2019 | B2 |
| 20020035864 | Paltieli et al. | Mar 2002 | A1 |
| 20020087071 | Schmitz et al. | Jul 2002 | A1 |
| 20020111568 | Bukshpan | Aug 2002 | A1 |
| 20020138003 | Bukshpan | Sep 2002 | A1 |
| 20020161299 | Prater et al. | Oct 2002 | A1 |
| 20030013962 | Bjaerum et al. | Jan 2003 | A1 |
| 20030028111 | Vaezy et al. | Feb 2003 | A1 |
| 20030040669 | Grass et al. | Feb 2003 | A1 |
| 20030228053 | Li et al. | Dec 2003 | A1 |
| 20040015079 | Berger et al. | Jan 2004 | A1 |
| 20040054283 | Corey et al. | Mar 2004 | A1 |
| 20040068184 | Trahey et al. | Apr 2004 | A1 |
| 20040100163 | Baumgartner et al. | May 2004 | A1 |
| 20040111028 | Abe et al. | Jun 2004 | A1 |
| 20040122313 | Moore et al. | Jun 2004 | A1 |
| 20040122322 | Moore et al. | Jun 2004 | A1 |
| 20040127793 | Mendlein et al. | Jul 2004 | A1 |
| 20040138565 | Trucco | Jul 2004 | A1 |
| 20040144176 | Yoden | Jul 2004 | A1 |
| 20040215075 | Zagzebski et al. | Oct 2004 | A1 |
| 20040236217 | Cerwin et al. | Nov 2004 | A1 |
| 20040236223 | Barnes et al. | Nov 2004 | A1 |
| 20040267132 | Podany | Dec 2004 | A1 |
| 20050004449 | Mitschke et al. | Jan 2005 | A1 |
| 20050053305 | Li et al. | Mar 2005 | A1 |
| 20050054910 | Tremblay et al. | Mar 2005 | A1 |
| 20050061536 | Proulx | Mar 2005 | A1 |
| 20050090743 | Kawashima et al. | Apr 2005 | A1 |
| 20050090745 | Steen | Apr 2005 | A1 |
| 20050111846 | Steinbacher et al. | May 2005 | A1 |
| 20050113689 | Gritzky | May 2005 | A1 |
| 20050113694 | Haugen et al. | May 2005 | A1 |
| 20050124883 | Hunt | Jun 2005 | A1 |
| 20050131300 | Bakircioglu et al. | Jun 2005 | A1 |
| 20050147297 | McLaughlin et al. | Jul 2005 | A1 |
| 20050165312 | Knowles et al. | Jul 2005 | A1 |
| 20050203404 | Freiburger | Sep 2005 | A1 |
| 20050215883 | Hundley et al. | Sep 2005 | A1 |
| 20050240125 | Makin et al. | Oct 2005 | A1 |
| 20050252295 | Fink et al. | Nov 2005 | A1 |
| 20050281447 | Moreau-Gobard et al. | Dec 2005 | A1 |
| 20050288588 | Weber et al. | Dec 2005 | A1 |
| 20060058664 | Barthe et al. | Mar 2006 | A1 |
| 20060062447 | Rinck et al. | Mar 2006 | A1 |
| 20060074313 | Slayton et al. | Apr 2006 | A1 |
| 20060074315 | Liang et al. | Apr 2006 | A1 |
| 20060074320 | Yoo et al. | Apr 2006 | A1 |
| 20060079759 | Vaillant et al. | Apr 2006 | A1 |
| 20060079778 | Mo et al. | Apr 2006 | A1 |
| 20060079782 | Beach et al. | Apr 2006 | A1 |
| 20060094962 | Clark | May 2006 | A1 |
| 20060111634 | Wu | May 2006 | A1 |
| 20060122506 | Davies et al. | Jun 2006 | A1 |
| 20060173327 | Kim | Aug 2006 | A1 |
| 20060262961 | Holsing et al. | Nov 2006 | A1 |
| 20060270934 | Savord et al. | Nov 2006 | A1 |
| 20070016022 | Blalock et al. | Jan 2007 | A1 |
| 20070016044 | Blalock et al. | Jan 2007 | A1 |
| 20070036414 | Georgescu et al. | Feb 2007 | A1 |
| 20070055155 | Owen et al. | Mar 2007 | A1 |
| 20070073781 | Adkins et al. | Mar 2007 | A1 |
| 20070078345 | Mo et al. | Apr 2007 | A1 |
| 20070088213 | Poland | Apr 2007 | A1 |
| 20070138157 | Dane et al. | Jun 2007 | A1 |
| 20070161898 | Hao et al. | Jul 2007 | A1 |
| 20070161904 | Urbano | Jul 2007 | A1 |
| 20070167752 | Proulx et al. | Jul 2007 | A1 |
| 20070167824 | Lee et al. | Jul 2007 | A1 |
| 20070232914 | Chen et al. | Oct 2007 | A1 |
| 20070238985 | Smith et al. | Oct 2007 | A1 |
| 20070242567 | Daft et al. | Oct 2007 | A1 |
| 20080110261 | Randall et al. | May 2008 | A1 |
| 20080110263 | Klessel et al. | May 2008 | A1 |
| 20080112265 | Urbano et al. | May 2008 | A1 |
| 20080114241 | Randall et al. | May 2008 | A1 |
| 20080114245 | Randall et al. | May 2008 | A1 |
| 20080114246 | Randall et al. | May 2008 | A1 |
| 20080114247 | Urbano et al. | May 2008 | A1 |
| 20080114248 | Urbano et al. | May 2008 | A1 |
| 20080114249 | Randall et al. | May 2008 | A1 |
| 20080114250 | Urbano et al. | May 2008 | A1 |
| 20080114251 | Weymer et al. | May 2008 | A1 |
| 20080114252 | Randall et al. | May 2008 | A1 |
| 20080114253 | Randall et al. | May 2008 | A1 |
| 20080114255 | Schwartz et al. | May 2008 | A1 |
| 20080125659 | Wilser et al. | May 2008 | A1 |
| 20080181479 | Yang et al. | Jul 2008 | A1 |
| 20080183075 | Govari et al. | Jul 2008 | A1 |
| 20080188747 | Randall et al. | Aug 2008 | A1 |
| 20080188750 | Randall et al. | Aug 2008 | A1 |
| 20080194957 | Hoctor et al. | Aug 2008 | A1 |
| 20080194958 | Lee et al. | Aug 2008 | A1 |
| 20080194959 | Wang et al. | Aug 2008 | A1 |
| 20080208061 | Halmann | Aug 2008 | A1 |
| 20080242996 | Hall et al. | Oct 2008 | A1 |
| 20080249408 | Palmeri et al. | Oct 2008 | A1 |
| 20080255452 | Entrekin | Oct 2008 | A1 |
| 20080269604 | Boctor et al. | Oct 2008 | A1 |
| 20080269613 | Summers et al. | Oct 2008 | A1 |
| 20080275344 | Glide-Hurst et al. | Nov 2008 | A1 |
| 20080285819 | Konofagou et al. | Nov 2008 | A1 |
| 20080287787 | Sauer et al. | Nov 2008 | A1 |
| 20080294045 | Ellington et al. | Nov 2008 | A1 |
| 20080294050 | Shinomura et al. | Nov 2008 | A1 |
| 20080294052 | Wilser et al. | Nov 2008 | A1 |
| 20080306382 | Guracar et al. | Dec 2008 | A1 |
| 20080306386 | Baba et al. | Dec 2008 | A1 |
| 20080319317 | Kamiyama et al. | Dec 2008 | A1 |
| 20090010459 | Garbini et al. | Jan 2009 | A1 |
| 20090012393 | Choi | Jan 2009 | A1 |
| 20090015665 | Willsie | Jan 2009 | A1 |
| 20090016163 | Freeman et al. | Jan 2009 | A1 |
| 20090018445 | Schers et al. | Jan 2009 | A1 |
| 20090024039 | Wang et al. | Jan 2009 | A1 |
| 20090036780 | Abraham | Feb 2009 | A1 |
| 20090043206 | Towfiq et al. | Feb 2009 | A1 |
| 20090048519 | Hossack et al. | Feb 2009 | A1 |
| 20090069681 | Lundberg et al. | Mar 2009 | A1 |
| 20090069686 | Daft et al. | Mar 2009 | A1 |
| 20090069692 | Cooley et al. | Mar 2009 | A1 |
| 20090079299 | Bradley et al. | Mar 2009 | A1 |
| 20090099483 | Rybyanets | Apr 2009 | A1 |
| 20090112095 | Daigle | Apr 2009 | A1 |
| 20090131797 | Jeong et al. | May 2009 | A1 |
| 20090143680 | Yao et al. | Jun 2009 | A1 |
| 20090148012 | Altmann et al. | Jun 2009 | A1 |
| 20090150094 | Van Velsor et al. | Jun 2009 | A1 |
| 20090182233 | Wodnicki | Jul 2009 | A1 |
| 20090182237 | Angelsen et al. | Jul 2009 | A1 |
| 20090198134 | Hashimoto et al. | Aug 2009 | A1 |
| 20090199392 | Singh et al. | Aug 2009 | A1 |
| 20090203997 | Ustuner | Aug 2009 | A1 |
| 20090208080 | Grau et al. | Aug 2009 | A1 |
| 20090259128 | Stribling | Oct 2009 | A1 |
| 20090264760 | Lazebnik et al. | Oct 2009 | A1 |
| 20090306510 | Hashiba et al. | Dec 2009 | A1 |
| 20090326379 | Daigle et al. | Dec 2009 | A1 |
| 20100010354 | Skerl et al. | Jan 2010 | A1 |
| 20100016725 | Thiele | Jan 2010 | A1 |
| 20100036258 | Dietz et al. | Feb 2010 | A1 |
| 20100063397 | Wagner | Mar 2010 | A1 |
| 20100063399 | Walker et al. | Mar 2010 | A1 |
| 20100069751 | Hazard et al. | Mar 2010 | A1 |
| 20100069756 | Ogasawara et al. | Mar 2010 | A1 |
| 20100085383 | Cohen et al. | Apr 2010 | A1 |
| 20100106431 | Baba et al. | Apr 2010 | A1 |
| 20100109481 | Buccafusca | May 2010 | A1 |
| 20100121193 | Fukukita et al. | May 2010 | A1 |
| 20100121196 | Hwang et al. | May 2010 | A1 |
| 20100130855 | Lundberg et al. | May 2010 | A1 |
| 20100145195 | Hyun | Jun 2010 | A1 |
| 20100168566 | Bercoff et al. | Jul 2010 | A1 |
| 20100168578 | Garson, Jr. et al. | Jul 2010 | A1 |
| 20100174194 | Chiang et al. | Jul 2010 | A1 |
| 20100174198 | Young et al. | Jul 2010 | A1 |
| 20100191110 | Insana et al. | Jul 2010 | A1 |
| 20100217124 | Cooley | Aug 2010 | A1 |
| 20100228126 | Emery et al. | Sep 2010 | A1 |
| 20100240994 | Zheng | Sep 2010 | A1 |
| 20100249570 | Carson et al. | Sep 2010 | A1 |
| 20100249596 | Magee | Sep 2010 | A1 |
| 20100256488 | Kim et al. | Oct 2010 | A1 |
| 20100262013 | Smith et al. | Oct 2010 | A1 |
| 20100266176 | Masumoto et al. | Oct 2010 | A1 |
| 20100286525 | Osumi | Nov 2010 | A1 |
| 20100286527 | Cannon et al. | Nov 2010 | A1 |
| 20100310143 | Rao et al. | Dec 2010 | A1 |
| 20100317971 | Fan et al. | Dec 2010 | A1 |
| 20100324418 | El-Aklouk et al. | Dec 2010 | A1 |
| 20100324423 | El-Aklouk et al. | Dec 2010 | A1 |
| 20100329521 | Beymer et al. | Dec 2010 | A1 |
| 20110005322 | Ustuner | Jan 2011 | A1 |
| 20110016977 | Guracar | Jan 2011 | A1 |
| 20110021920 | Shafir et al. | Jan 2011 | A1 |
| 20110021923 | Daft et al. | Jan 2011 | A1 |
| 20110033098 | Richter et al. | Feb 2011 | A1 |
| 20110044133 | Tokita | Feb 2011 | A1 |
| 20110066030 | Yao | Mar 2011 | A1 |
| 20110098565 | Masuzawa | Apr 2011 | A1 |
| 20110112400 | Emery et al. | May 2011 | A1 |
| 20110112404 | Gourevitch | May 2011 | A1 |
| 20110125017 | Ramamurthy et al. | May 2011 | A1 |
| 20110178441 | Tyler | Jul 2011 | A1 |
| 20110270088 | Shiina | Nov 2011 | A1 |
| 20110301470 | Sato et al. | Dec 2011 | A1 |
| 20110306886 | Daft et al. | Dec 2011 | A1 |
| 20110319764 | Okada et al. | Dec 2011 | A1 |
| 20120004545 | Ziv-Ari et al. | Jan 2012 | A1 |
| 20120035482 | Kim et al. | Feb 2012 | A1 |
| 20120036934 | Kröning et al. | Feb 2012 | A1 |
| 20120085173 | Papadopoulos et al. | Apr 2012 | A1 |
| 20120101378 | Lee | Apr 2012 | A1 |
| 20120114210 | Kim et al. | May 2012 | A1 |
| 20120121150 | Murashita | May 2012 | A1 |
| 20120137778 | Kitazawa et al. | Jun 2012 | A1 |
| 20120140595 | Amemiya | Jun 2012 | A1 |
| 20120141002 | Urbano et al. | Jun 2012 | A1 |
| 20120165670 | Shi et al. | Jun 2012 | A1 |
| 20120179044 | Chiang et al. | Jul 2012 | A1 |
| 20120226201 | Clark et al. | Sep 2012 | A1 |
| 20120235998 | Smith-Casem et al. | Sep 2012 | A1 |
| 20120243763 | Wen et al. | Sep 2012 | A1 |
| 20120253194 | Tamura | Oct 2012 | A1 |
| 20120265075 | Pedrizzetti et al. | Oct 2012 | A1 |
| 20120277585 | Koenig et al. | Nov 2012 | A1 |
| 20130070062 | Fouras et al. | Mar 2013 | A1 |
| 20130076207 | Krohn et al. | Mar 2013 | A1 |
| 20130079639 | Hoctor et al. | Mar 2013 | A1 |
| 20130083628 | Qiao et al. | Apr 2013 | A1 |
| 20130088122 | Krohn et al. | Apr 2013 | A1 |
| 20130116561 | Rothberg et al. | May 2013 | A1 |
| 20130131516 | Katsuyama | May 2013 | A1 |
| 20130144165 | Ebbini et al. | Jun 2013 | A1 |
| 20130204136 | Duric et al. | Aug 2013 | A1 |
| 20130204137 | Roy et al. | Aug 2013 | A1 |
| 20130258805 | Hansen et al. | Oct 2013 | A1 |
| 20130261463 | Chiang et al. | Oct 2013 | A1 |
| 20140058266 | Call | Feb 2014 | A1 |
| 20140073921 | Specht et al. | Mar 2014 | A1 |
| 20140086014 | Kobayashi | Mar 2014 | A1 |
| 20140147013 | Shandas et al. | May 2014 | A1 |
| 20140243673 | Anand et al. | Aug 2014 | A1 |
| 20150045668 | Smith et al. | Feb 2015 | A1 |
| 20150172878 | Luna et al. | Jun 2015 | A1 |
| 20160095579 | Smith et al. | Apr 2016 | A1 |
| 20160135783 | Brewer et al. | May 2016 | A1 |
| 20160256134 | Specht et al. | Sep 2016 | A1 |
| 20170209121 | Davies et al. | Jul 2017 | A1 |
| 20170224312 | Call et al. | Aug 2017 | A1 |
| 20180049717 | Adam et al. | Feb 2018 | A1 |
| 20180153511 | Specht et al. | Jun 2018 | A1 |
| 20180279991 | Call et al. | Oct 2018 | A1 |
| 20190008487 | Belevich et al. | Jan 2019 | A1 |
| 20190083058 | Specht | Mar 2019 | A1 |
| 20190175152 | Smith et al. | Jun 2019 | A1 |
| 20190200961 | Specht et al. | Jul 2019 | A1 |
| 20190328367 | Specht et al. | Oct 2019 | A1 |
| 20190370522 | Call et al. | Dec 2019 | A1 |
| Number | Date | Country |
|---|---|---|
| 1535243 | Oct 2004 | CN |
| 1781460 | Jun 2006 | CN |
| 101103927 | Jan 2008 | CN |
| 101116622 | Feb 2008 | CN |
| 101190134 | Jun 2008 | CN |
| 101453955 | Jun 2009 | CN |
| 100545650 | Sep 2009 | CN |
| 101609150 | Dec 2009 | CN |
| 101843501 | Sep 2010 | CN |
| 101912278 | Dec 2010 | CN |
| 102018533 | Apr 2011 | CN |
| 102112047 | Jun 2011 | CN |
| 102123668 | Jul 2011 | CN |
| 102599930 | Jul 2012 | CN |
| 102011114333 | Mar 2013 | DE |
| 1949856 | Jul 2008 | EP |
| 2058796 | May 2009 | EP |
| 2101191 | Sep 2009 | EP |
| 2182352 | May 2010 | EP |
| 2187813 | May 2010 | EP |
| 2198785 | Jun 2010 | EP |
| 1757955 | Nov 2010 | EP |
| 2325672 | May 2011 | EP |
| 1462819 | Jul 2011 | EP |
| 2356941 | Aug 2011 | EP |
| 1979739 | Oct 2011 | EP |
| 2385391 | Nov 2011 | EP |
| 2294400 | Feb 2012 | EP |
| 2453256 | May 2012 | EP |
| 1840594 | Jun 2012 | EP |
| 2514368 | Oct 2012 | EP |
| 1850743 | Dec 2012 | EP |
| 1594404 | Sep 2013 | EP |
| 2026280 | Oct 2013 | EP |
| 2851662 | Aug 2004 | FR |
| 49-11189 | Jan 1974 | JP |
| 54-44375 | Apr 1979 | JP |
| 55-103839 | Aug 1980 | JP |
| 57-31848 | Feb 1982 | JP |
| 58-223059 | Dec 1983 | JP |
| 59-101143 | Jun 1984 | JP |
| 59-174151 | Oct 1984 | JP |
| 60-13109 | Jan 1985 | JP |
| 60-68836 | Apr 1985 | JP |
| 01164354 | Jun 1989 | JP |
| 02501431 | May 1990 | JP |
| 03015455 | Jan 1991 | JP |
| 03126443 | May 1991 | JP |
| 04017842 | Jan 1992 | JP |
| 04067856 | Mar 1992 | JP |
| 05042138 | Feb 1993 | JP |
| 06125908 | May 1994 | JP |
| 06254092 | Sep 1994 | JP |
| 07051266 | Feb 1995 | JP |
| 07204201 | Aug 1995 | JP |
| 08154930 | Jun 1996 | JP |
| 08252253 | Oct 1996 | JP |
| 09103429 | Apr 1997 | JP |
| 09201361 | Aug 1997 | JP |
| 2777197 | May 1998 | JP |
| 10216128 | Aug 1998 | JP |
| 11089833 | Apr 1999 | JP |
| 11239578 | Sep 1999 | JP |
| 2001507794 | Jun 2001 | JP |
| 2001245884 | Sep 2001 | JP |
| 2002209894 | Jul 2002 | JP |
| 2002253548 | Sep 2002 | JP |
| 2002253549 | Sep 2002 | JP |
| 2003235839 | Aug 2003 | JP |
| 2004167092 | Jun 2004 | JP |
| 2004215987 | Aug 2004 | JP |
| 2004337457 | Dec 2004 | JP |
| 2004340809 | Dec 2004 | JP |
| 2004351214 | Dec 2004 | JP |
| 2005046192 | Feb 2005 | JP |
| 2005152187 | Jun 2005 | JP |
| 2005523792 | Aug 2005 | JP |
| 2005526539 | Sep 2005 | JP |
| 2006051356 | Feb 2006 | JP |
| 2006061203 | Mar 2006 | JP |
| 2006122657 | May 2006 | JP |
| 2006130313 | May 2006 | JP |
| 2006204923 | Aug 2006 | JP |
| 2007325937 | Dec 2007 | JP |
| 2008122209 | May 2008 | JP |
| 2008513763 | May 2008 | JP |
| 2008515557 | May 2008 | JP |
| 2008132342 | Jun 2008 | JP |
| 2008522642 | Jul 2008 | JP |
| 2008259541 | Oct 2008 | JP |
| 2008279274 | Nov 2008 | JP |
| 2008307087 | Dec 2008 | JP |
| 2009240667 | Oct 2009 | JP |
| 2010005375 | Jan 2010 | JP |
| 2010124842 | Jun 2010 | JP |
| 2010526626 | Aug 2010 | JP |
| 2011529362 | Dec 2011 | JP |
| 2013121493 | Jun 2013 | JP |
| 2014087448 | May 2014 | JP |
| 100715132 | Apr 2007 | KR |
| 1020080044737 | May 2008 | KR |
| 1020090103408 | Oct 2009 | KR |
| WO9218054 | Oct 1992 | WO |
| WO9800719 | Jan 1998 | WO |
| WO0164109 | Sep 2001 | WO |
| WO02084594 | Oct 2002 | WO |
| WO2005009245 | Feb 2005 | WO |
| WO2006114735 | Nov 2006 | WO |
| WO2007127147 | Nov 2007 | WO |
| WO2008097479 | Aug 2008 | WO |
| WO2009060182 | May 2009 | WO |
| WO2010095094 | Aug 2010 | WO |
| WO2010137453 | Dec 2010 | WO |
| WO2010139519 | Dec 2010 | WO |
| WO2011004661 | Jan 2011 | WO |
| WO2011057252 | May 2011 | WO |
| WO2011064688 | Jun 2011 | WO |
| WO2011100697 | Aug 2011 | WO |
| WO2011123529 | Oct 2011 | WO |
| WO2012028896 | Mar 2012 | WO |
| WO2012049124 | Apr 2012 | WO |
| WO2012049612 | Apr 2012 | WO |
| WO2012078639 | Jun 2012 | WO |
| WO2012091280 | Jul 2012 | WO |
| WO2012112540 | Aug 2012 | WO |
| WO2012131340 | Oct 2012 | WO |
| WO2012160541 | Nov 2012 | WO |
| WO2013059358 | Apr 2013 | WO |
| WO2013109965 | Jul 2013 | WO |
| WO2013116807 | Aug 2013 | WO |
| WO2013116809 | Aug 2013 | WO |
| WO2013116851 | Aug 2013 | WO |
| WO2013116854 | Aug 2013 | WO |
| WO2013116866 | Aug 2013 | WO |
| WO2013128301 | Sep 2013 | WO |
| Entry |
|---|
| Abeysekera et al.; Alignment and calibration of dual ultrasound transducers using a wedge phantom; Ultrasound in Medicine and Biology; 37(2); pp. 271-279; Feb. 2011. |
| Arigovindan et al.; Full motion and flow field recovery from echo doppler data; IEEE Transactions on Medical Imaging; 26(1); pp. 31-45; Jan. 2007. |
| Cai et al.; Off-axis directional acoustic wave beaming control by an asymmetric rubber heterostructures film deposited on steel plate in water; IEEE Intl.; 2009 Ultrasonics Symposium (IUS); pp. 1552-1554; Rome; Sep. 2009. |
| Capineri et al.; A doppler system for dynamic vector velocity maps; Ultrasound in Medicine & Biology; 28(2); pp. 237-248; Feb. 28, 2002. |
| Carson et al.; Measurement of photoacoustic transducer position by robotic source placement and nonlinear parameter estimation; Biomedical Optics (BiOS); International Society for Optics and Photonics (9th Conf. on Biomedical Thermoacoustics, Optoacoustics, and Acousto-optics; vol. 6856; 9 pages; Feb 28, 2008. |
| Chen et al.; Maximum-likelihood source localization and unknown sensor location estimation for wideband signals in the near-field; IEEE Transactions on Signal Processing; 50(8); pp. 1843-1854; Aug. 2002. |
| Chen et al.; Source localization and tracking of a wideband source using a randomly distributed beamforming sensor array; International Journal of High Performance Computing Applications; 16(3); pp. 259-272; Fall 2002. |
| Cristianini et al.; An Introduction to Support Vector Machines; Cambridge University Press; pp. 93-111; Mar. 2000. |
| Dunmire et al.; A brief history of vector doppler; Medical Imaging 2001; International Society for Optics and Photonics; pp. 200-214; May 30, 2001. |
| Du et al.; User parameter free approaches to multistatic adaptive ultrasound imaging; 5th IEEE International Symposium; pp. 1287-1290, May 2008. |
| Feigenbaum, Harvey, M.D.; Echocardiography; Lippincott Williams & Wilkins; Philadelphia; 5th Ed.; pp. 482, 484; Feb. 1994. |
| Fernandez et al.; High resolution ultrasound beamforming using synthetic and adaptive imaging techniques; Proceedings IEEE International Symposium on Biomedical Imaging; Washington, D.C.; pp. 433-436; Jul. 7-10, 2002. |
| Gazor et al.; Wideband multi-source beamforming with array location calibration and direction finding; Conference on Acoustics, Speech and Signal Processing ICASSP-95; Detroit, MI; vol. 3 IEEE; pp. 1904-1907; May 9-12, 1995. |
| Hasegawa et al.; High frame rate ultrasonic imaging of the heart by placing virtual point sources in front of array; IEEE Intl.; 2013 Ultrasonics Symposium (IUS); pp. 581-584; Prague; Jul. 2013. |
| Haun et al.; Efficient three-dimensional imaging from a small cylindrical aperture; IEEE Trans. on Ultrasonics, Ferroelectrics, and Frequency Control; 49(7); pp. 861-870; Jul. 2002. |
| Haykin, Simon; Neural Networks: A Comprehensive Foundation (2nd Ed.); Prentice Hall; pp. 156-187; Jul. 16, 1998. |
| Heikkila et al.; A four-step camera calibration procedure with implicit image correction; Proceedings IEEE Computer Scociety Conference on Computer Vision and Pattern Recognition; San Juan; pp. 1106-1112; Jun. 17-19,1997. |
| Hendee et al.; Medical Imaging Physics; Wiley-Liss, Inc. 4th Edition; Chap. 19-22; pp. 303-353; (year of pub. sufficiently earlier than effective US filing date and any foreign priority date) © 2002. |
| Hsu et al.; Real-time freehand 3D ultrasound calibration; CUED/F-INFENG/TR 565; Department of Engineering, University of Cambridge, United Kingdom; 14 pages; Sep. 2006. |
| Jeffs; Beamforming: a brief introduction; Brigham Young University; 14 pages; retrieved from the internet (http://ens.ewi.tudelft.nl/Education/courses/et4235/Beamforming.pdf); Oct. 2004. |
| Khamene et al.; A novel phantom-less spatial and temporal ultrasound calibration method; Medical Image Computing and Computer-Assisted Intervention—MICCAI (Proceedings 8th Int. Conf.); Springer Berlin Heidelberg; Palm Springs, CA; pp. 65-72; Oct. 26-29, 2005. |
| Kramb et al,.; Considerations for using phased array ultrasonics in a fully automated inspection system. Review of Quantitative Nondestructive Evaluation, 2004 Edition, ed. D. O. Thompson and D. E. Chimenti, American Inst. of Physics, pp. 817-825, Mar. 2004. |
| Ledesma-Carbayo et al.; Spatio-temporal nonrigid registration for ultrasound cardiac motion estimation; IEEE Trans. on Medical Imaging; vol. 24; No. 9; Sep. 2005. |
| Leotta et al.; Quantitative three-dimensional echocardiography by rapid imaging . . . ; J American Society of Echocardiography; vol. 10; No. 8; ppl 830-839; Oct. 1997. |
| Li et al.; An efficient speckle tracking algorithm for ultrasonic imaging; 24; pp. 215-228; Oct. 1, 2002. |
| Montaldo et al.; Building three-diminsional images using a time-reversal chaotic cavity; IEEE Trans. on Ultrasonics, Ferroelectrics, and Frequency Control; 52(9); pp. 1489-1497; Sep. 2005. |
| Morrison et al.; A probabilistic neural network based image segmentation network for magnetic resonance images; Proc. Conf. Neural Networks; Baltimore, MD; vol. 3; pp. 60-65; Jun. 1992. |
| Nadkarni et al.; Cardiac motion synchronization for 3D cardiac ultrasound imaging; Ph.D. Dissertation, University of Western Ontario; Jun. 2002. |
| Opretzka et al.; A high-frequency ultrasound imaging system combining limited-angle spatial compounding and model-based synthetic aperture focusing; IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, IEEE, US; 58(7); pp. 1355-1365; Jul. 2, 2011. |
| Press et al.; Cubic spline interpolation; §3.3 in “Numerical Recipes in FORTRAN: The Art of Scientific Computing”, 2nd Ed.; Cambridge, England; Cambridge University Press; pp. 107-110; Sep. 1992. |
| Saad et al.; Computer vision approach for ultrasound doppler angle estimation; Journal of Digital Imaging; 22(6); pp. 681-688; Dec. 1, 2009. |
| Sakas et al.; Preprocessing and volume rendering of 3D ultrasonic data; IEEE Computer Graphics and Applications; pp. 47-54, Jul. 1995. |
| Sapia et al.; Deconvolution of ultrasonic waveforms using an adaptive wiener filter; Review of Progress in Quantitative Nondestructive Evaluation; vol. 13A; Plenum Press; pp. 855-862; Jan. 1994. |
| Sapia et al.; Ultrasound image deconvolution using adaptive inverse filtering; 12 IEEE Symposium on Computer-Based Medical Systems, CBMS, pp. 248-253; Jun. 1999. |
| Sapia, Mark Angelo; Multi-dimensional deconvolution of optical microscope and ultrasound imaging using adaptive least-mean-square (LMS) inverse filtering; Ph.D. Dissertation; University of Connecticut; Jan. 2000. |
| Slavine et al.; Construction, calibration and evaluation of a tissue phantom with reproducible optical properties for investigations in light emission tomography; Engineering in Medicine and Biology Workshop; Dallas, TX; IEEE pp. 122-125; Nov. 11-12, 2007. |
| Smith et al.; High-speed ultrasound volumetric imaging system. 1. Transducer design and beam steering; IEEE Trans. Ultrason., Ferroelect., Freq. Contr.; vol. 38; pp. 100-108; Mar. 1991. |
| Specht et al.; Deconvolution techniques for digital longitudinal tomography; SPIE; vol. 454; presented at Application of Optical Instrumentation in Medicine XII; pp. 319-325; Jun. 1984. |
| Specht et al.; Experience with adaptive PNN and adaptive GRNN; Proc. IEEE International Joint Conf. on Neural Networks; vol. 2; pp. 1203-1208; Orlando, FL; Jun. 1994. |
| Specht, D.F.; A general regression neural network; IEEE Trans. on Neural Networks; vol. 2.; No. 6; Nov. 1991. |
| Specht, D.F.; Blind deconvolution of motion blur using LMS inverse filtering; Lockheed Independent Research (unpublished); Jun. 23, 1975. |
| Specht, D.F.; Enhancements to probabilistic neural networks; Proc. IEEE International Joint Conf. on Neural Networks; Baltimore, MD; Jun. 1992. |
| Specht, D.F.; GRNN with double clustering; Proc. IEEE International Joint Conf. Neural Networks; Vancouver, Canada; Jul. 16-21, 2006. |
| Specht, D.F.; Probabilistic neural networks; Pergamon Press; Neural Networks; vol. 3; pp. 109-118; Feb. 1990. |
| UCLA Academic Technology; SPSS learning module: How can I analyze a subset of my data; 6 pages; retrieved from the internet (http://www.ats.ucla.edu/stat/spss/modules/subset_analyze.htm) Nov. 26, 2001. |
| Urban et al; Implementation of vibro-acoustography on a clinical ultrasound system; IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control; 58(6); pp. 1169-1181 (Author Manuscript, 25 pgs.); Jun. 2011. |
| Urban et al; Implementation of vibro-acoustography on a clinical ultrasound system; IEEE Ultrasonics Symposium (IUS); pp. 326-329; Oct. 14, 2010. |
| Von Ramm et al.; High-speed ultrasound volumetric imaging—System. 2. Parallel processing and image display; IEEE Trans. Ultrason., Ferroelect., Freq. Contr.; vol. 38; pp. 109-115; Mar. 1991. |
| Wang et al.; Photoacoustic tomography of biological tissues with high cross-section resolution: reconstruction and experiment; Medical Physics; 29(12); pp. 2799-2805; Dec. 2002. |
| Wells, P.N.T.; Biomedical ultrasonics; Academic Press; London, New York, San Francisco; pp. 124-125; Mar. 1977. |
| Widrow et al.; Adaptive signal processing; Prentice-Hall; Englewood Cliffs, NJ; pp. 99-116; Mar. 1985. |
| Wikipedia; Point cloud; 2 pages; retrieved Nov. 24, 2014 from the internet (https://en.wikipedia.org/w/index.php?title=Point_cloud&oldid=472583138). |
| Wikipedia; Curve fitting; 5 pages; retrieved from the internet (http:en.wikipedia.org/wiki/Curve_fitting) Dec. 19, 2010. |
| Wikipedia; Speed of sound; 17 pages; retrieved from the internet (http:en.wikipedia.org/wiki/Speed_of_sound) Feb. 15, 2011. |
| Yang et al.; Time-of-arrival calibration for improving the microwave breast cancer imaging; 2011 IEEE Topical Conf. on Biomedical Wireless Technologies, Networks, and sensing Systems (BioWireleSS); Phoenix, AZ; pp. 67-70; Jan. 16-19, 2011. |
| Zhang et al.; A high-frequency high frame rate duplex ultrasound linear array imaging system for small animal imaging; IEEE transactions on ultrasound, ferroelectrics, and frequency control; 57(7); pp. 1548-1567; Jul. 2010. |
| Number | Date | Country | |
|---|---|---|---|
| 20190021697 A1 | Jan 2019 | US |
| Number | Date | Country | |
|---|---|---|---|
| 61877555 | Sep 2013 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 14279052 | May 2014 | US |
| Child | 15888738 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 15888738 | Feb 2018 | US |
| Child | 16137221 | US |
