AU759136B2 - Improvements in ultrasound techniques - Google Patents

Description

WO 99/45839 -1- PCT/AU99/00157 Improvements in Ultrasound Techniques Field of the Invention The present invention relates to the field of ultrasound imaging and, in particular, to the utilisation of ultrasound data for producing volume flow information.

Background of the Invention In recent years, ultrasound imaging has become increasingly important, especially in the area of medical imaging of the internal portions of a person or animal's anatomy. The principles of ultrasound imaging are well known. However, recently the utilisation of phase data from an ultrasound signal has resulted in the production of flow information associated with an ultrasound scan. A number of manufacturers now offer "Colour Doppler Imaging" which provides information in the form of colour indicators showing the direction and speed of flow of blood within vessels.

It would be desirable to provide for non-invasive real time ultrasonic measurement of the blood volume flow rate in blood vessels with a high accuracy and high degree of robustness ideally utilizing standard hardware available with only software modification.

Summary of the Invention It is an object of the present invention to provide for the accurate measurement of flow volume within the vessels utilising ultrasound techniques.

In accordance with a first aspect of the present invention, there is provided a method of estimating the volume flow in a vessel utilising ultrasound techniques, the method comprising the steps of imaging the vessel utilising an ultrasonic transducer array so as to produce fluid velocity information for the vessel, the plane formed by the beams emitted by the transducers having a direction which is offset from the axis of the vessel such that the it forms an elliptical section within the vessel; and utilising the velocity information for points within the WO 99/45839 2 PCT/AU99/00157 vessel and the elliptical section to calculate a mean velocity for the fluid flow through the elliptical section.

Preferably, the mean velocity is calculated using a power weighted average of the velocity information, the power weighting being correlated to the intensity of the received signal at each point corresponding to the velocity information.

Preferably, the method further comprises the step of fitting an elliptical section to the fluid velocity information and determining the elliptical section parameters for the vessel from the fitted elliptical section.

The method is ideally suited to use on the internal portion of a human or animal's anatomy and can be implemented on a standard ultrasound machine by reprogramming of the software of the machine.

Brief Description of the Drawings Notwithstanding any other forms which may fall within the scope of the present invention, preferred forms of the invention will now be described, by way of example only, with reference to the accompanying drawings in which: Fig. 1 illustrates the standard longitudinal scan process; Fig. 2 illustrates a first rotation of the scan plane relative to a vessel; Fig. 3 illustrates a second rotation of the scan plane relative to a measured vessel; Fig. 4 6 illustrate the process of coordinate transforms utilised in the preferred embodiment; Fig. 7 illustrates the relationship between variables utilised in the preferred embodiment, and Fig. 8 illustrates the steps in the preferred embodiment.

Description of Preferred and Other Embodiments In the preferred embodiment, an oblique scan method is utilised to extract the details of an elliptical cross WO 99/45839 3 PCT/AU99/00157 sectional structure of a vessel and to utilise the elliptical structure, in conjunction with fluid movement information within the cross sectional structure, so as to derive a volume flow.

Turning initially to Fig. I, there is illustrated a schematic of the normal ultrasound processing technique.

The normal arrangement includes a transducer 1 which can include an array of transmitters and receivers arranged in a line. The transducer 1 emits a beam 2 which is of a substantially planar form having a predetermined thickness.

The beam 2 is projected through a vessel 3 with reflections being received back to the transducer 1 and subsequently processed so as to determine the structure of the vessel.

In Colour Doppler Imaging, a series of emissions are made in a periodic manner in rapid succession and the return reflections recorded. The phase difference between multiple returns is determined and utilised to calculate flow velocity information within the vessel 3 in a known manner.

In the preferred embodiment, an oblique scan is utilised where the orientation of the beam 2 does not line up with the vessel 3. In particular, as illustrated in Figs. 2 and 3, two forms of oblique scan beam rotation are provided. Firstly, as shown in Fig. 2 a rotation a of the transducer probe about an axis 7 that is perpendicular to the vessel longitudinal axis (ie. about the probe longitudinal axis) is utilised resulting in a first cross sectional volume 4. Secondly, as illustrated in Fig. 3, a rotation P around an axis 8 perpendicular to the scan plane is provided. This results in a final cross section 4 of the vessel 3 taking an elliptical sectional form being of a predetermined volume which is determined by the thickness characteristics of the transducer beam.

In the preferred embodiment a sequence of frames is imaged utilising the angled oblique scan of Fig. 2 and Fig.

3 and subsequently processed as discussed hereinafter to WO 99/45839 -4- PCT/AU99/00157 determine an actual volume flow rate through a vessel of interest.

Turning now to Fig. 8, there is illustrated the steps in the method of the preferred embodiment. The method proceeds 11 by initially utilizing standard modern ultrasound equipment to capture a series of frames using an oblique scan method.

From the captured sequence of frames, corresponding velocity information is obtained 12 utilising the known Colour Doppler techniques, the method of obtaining the colour flow data being entirely standard.

Next, the colour flow images are examined so as to identify an elliptical cross section and from this, the geometric parameters of the scan (a and P) are obtained 13.

The ratio of the ellipsoid minor to major axes is utilised to calculate the angle a whilst the orientation of the major axis relative to the horizontal in the scan and the colour Doppler steering angle can be used to determine f.

Once the overall elliptical structure has been determined, it is subsequently utilised to compute the cross sectional area of the vessel. Subsequently, with the knowledge of the geometry and the mean velocity, a volume flow rate estimate can be calculated 14.

A more detailed understanding of the process can be gained by examining the steps used to mathematically simulate the scanning and signal processing.

In order to determine and identify the ellipse parameters and to calculate the flow volume, the image plane received by the transducer element is scanned sequentially and for each scan coordinate, the corresponding vessel coordinate is obtained for processing so as to determine whether the sample coordinate (which represents a sample voxel) is within the vessel lumen or otherwise.

It is therefore necessary to transform between the WO 99/45839 5 PCT/AU99/00157 scan plane coordinate space and the vessel coordinate space. The transformation process will now be discussed with reference to Figs. 4 to 6. In Fig. 4, there is illustrated an x, y, z coordinate system of the vessel.

Fig. 4 also shows the transducer 1 aligned in the z-y plane. In the angle oblique scan process the vessel coordinates remain in the x, y, z, frame of reference whilst the scanning plane coordinates are transformed into frame of reference utilising two rotations.

As shown in Fig. 5, the first rotation is a rotation of the transducer 1 by a degrees around the z axis as illustrated in Fig. 5. Utilising a matrix notation, we can write this translation as: cos Z cos Z cos Z y' cos Z cos Z cos Z y (1) L cos Z cos Z (y,z')cos Z zJ Using the nomenclature from Fig. 5, the direction cosines are as follows: x' cos cos (7 2 a) cos 2) y' cos(/2-a) cos cos() os 2) y (2) z' cos (zr/2) cos cos zJ ie.

cos sin 0 y' sin cos 0 y (3) z' 0 0 1 z In the second transformation illustrated in Fig. 6, we go from the frame of reference to the frame of reference through the rotation of the transducer by P degrees. The transformation is expressed as follows: cos Z cos Z cosZ r' 1 y" cos Z cos Z cosZ y' (4) z"j cos z cos cosZ z' Using the nomenclature illustrated in Fig. 6 for the direction cosines, one obtains: WO 99/45839 6 PCT/AU99/00157 x1 cos cos cos( r/2) lrx' y" ]cos{ir 2) cos Cos (jz 2 y'( z cos cos (r 2-f8) cos(f) L' ie.

x 1 0 0 x' y" 0 cos -sin y' (6) Lz"J 0 +sin cos (ft) From equations and it is possible to combine the transformations to obtain a single transformation from z-y-z to X 1 0 0 ]cos -sin oTl cos(f) -sin (fi) sin cos 0Y (7) Lz" 0 sin cos 0 1I which gives: [x cos sin 0 y" cos sin cos (f)cos sin y (8) sin sin sin ()cos cos z or, through matrix inversion, [cos cos sin sin sin Fx.:1 -sin cos cos sin cos (9) LZ- 0 -sin cos 1zI Interestingly, the matrices in equations and (9) are the transposes of one another.

In the "idealised" simulation, the position of the sample volume is incremented through the extent of the scan-plane and the corresponding vessel coordinates calculated using equation If the sample volume lies within the vessel lumen (ie, if x 2 z 2 Radius )then a corresponding power value is assigned a value of 1.0 and the velocity assigned a value based on a radial distance from the vessel axis.

Unfortunately, the aforementioned analysis ignores the fact the transducer beam has a finite width and is likely to have an approximately Gaussian intensity cross section.

WO 99/45839 -7- PCT/AU99/00157 Hence, it is highly desirable to include the ultrasound parameters in the calculation of both power and velocity values for each sample volume or voxel location.

Geometry Determination (Ellipse Fitting) (Step 13 of Fig.

8) The first step in this process is determination of the geometry, which requires the fitting of an ellipse to the power or velocity data. The data utilised to fit the ellipse can include a combination of the power data, the velocity data and the grey scale information. In one example, the power data is utilized to fit the ellipse, suppressing power values below a threshold value and suppressing areas having significant grey scale values.

The algorithm used for ellipse fitting is based on a moment of inertia formula. For velocity data, only values above a specified velocity threshold are used and each value carries equal weight in the calculation of the moment of inertia. On the other hand, when fitting an ellipse to power data, values above the specified power threshold may be taken as weights in calculating the moment of inertia.

Thus, higher power regions may contribute more to the calculation than lower power areas.

Having fitting the ellipse, the angles a and P are calculated as follows: a sin- 1 (12) where a length of major axis of ellipse and b length of minor axis of ellipse Sdiameter of vessel S= angle between the major axis of the ellipse and the horizontal in the scan.

From a and P, the "Doppler angle" correction can be calculated as follows: cos (ODoppler)= COS cos((P-) where p is the Colour Doppler steering angle which can typically be say 70*, 90' or 1100 but is not limited WO 99/45839 8 PCT/AU99/00157 thereto. The relationship between a, b, p and p is as illustrated in Fig. 7.

The true area of the cross section (perpendicular to the axis of the vessel) can be calculated using angle a as follows: A.ucu, A ,,,ered. sin(a) (14) Mean Velocity Calculation (Step 14 of Fig. 8) As noted previously, once the geometric parameters have been measured from the elliptical cross section, the next step in the preferred embodiment is to compute the mean velocity over the volume bounded by the ellipse.

Preferably, a power weighted summation process as described by the following formula is used: v Pv" z" SP(y" z" wihin ellipse where P represents the power values received and v represents the received velocity values.

This scheme gives reduced weight to voxels with lower power, for example those which are only partially within the vessel. The power-weighted approach, therefore, addresses the "partial volume" effect related to the finite voxel size.

Naturally, angle correction is necessary to obtain the true mean-velocity from the mean measured velocity. Such a correction can be applied as follows: V measured V (16) cos(Dop,,cer) Volume Flow Rate Estimation (Step 14 of Fig. 8) Given the geometric parameters a, P, the measured area and the mean velocity, the volume flow rate Q can be calculated with the aid of equations (14) and (16) as follows: WO 99/45839 9 PCT/AU99/00157 V neualured Q A= A (17) cos (0 ,,ppler Vnea.,ure A ea.,,red sin tan (a) cos(a) cos((o cos(p -fl) It can be seen from the foregoing that with the utilisation of a power weighting system a more accurate velocity and flow volume result is provided.

Of course, many modifications are possible. One such modification would be to utilize a "predicted power weighting" of a modelled power weighting for a particular ellipsoid sectional structure. This alleviates the possible problems associated with noisy power data.

As a further modification, the positional data can be utilized to process the flow information to screen for tissue echoes etc. This modification can include the examination of the traditional grey scale data usually provided on ultrasound machine and using this information to suppress velocity data values in the area of significant tissue echoes due to their likely spurious data.

It would be appreciated by a person skilled in the art that numerous other variations and/or modifications may be made to the present invention as shown in the specific embodiments without departing from the spirit or scope of the invention as broadly described. The present embodiments are, therefore, to be considered in all respects to be illustrative and not restrictive.

Claims (8)

1. A method of estimating the volume flow in a vessel utilising ultrasound techniques, said method comprising the steps of: imaging said vessel utilising an ultrasonic transducer array emitting a beam so as to produce fluid velocity information for said vessel, the plane formed by the beams emitted by said transducers having a direction which is offset from the axis of said vessel such that said ultrasonic transducer array images an elliptical section within said vessel; and utilising said velocity information for points within said vessel and said elliptical section to calculate a mean velocity for said fluid flow through the elliptical section.

2. A method as claimed in claim 1 further comprising the step of fitting an elliptical section to said fluid velocity information and determining the elliptical section parameters for said vessel from said fitted elliptical section.

3. A method as claimed in claim 1 or 2 wherein said vessel comprises an internal portion of a human or animal's anatomy.

4. A method as claimed in any previous claim wherein said method is implemented on a standard ultrasound machine by reprogramming of the software of said machine.

A method as claimed in any previous claim wherein said utilising step comprises utilising a power weighted average of said velocity information, said power weighting being correlated to the intensity of said received signal at a point corresponding to said velocity information.

6. A method as claimed in any previous claim wherein said velocity information is preprocessed so as to suppress velocity data values in areas of significant grey-scale echoes from said ultrasonic transducers.

7. A method as claimed in claim 5 wherein power WO 99/45839 -11- PCT/AU99/00157 weightings less than a predetermined magnitude are discarded.

8. An apparatus when implementing the method of any of claims 1 to 8.