TOE Principles & Doppler  

How echocardiography works

Echocardiography provides a picture assimilated from the faint echoes that arise when short pulses of high frequency sound (ultrasound) are passed through tissues. The computation needed to produce quality pictures is non-trivial, but the fundamental principles can be understood by anyone with a high-school education.

Fundamentals

Consider an ultrasound transducer. This is usually a piezoelectric crystal - a crystal with the peculiar property that, when voltage is put across it, it deforms ever so slightly, and that when it is deformed (ever so slightly!) the voltage across it changes. A little thought will lead one to the deduction that if one varies the voltage across the crystal, a sound wave will be generated, as all sound really is is a pressure wave that passes through a medium. A crystal that is repetitively thickening and thinning out again must cause local changes in pressure - a sound! What is really sneaky is that reflected sound will cause a minuscule distortion when it strikes the crystal, and the consequent change in voltage can be picked up and interpreted. The timing of this echo will give us an indication of the depth at which it occurred, in the same way that German submarines were detected and their distance estimated by wartime sonar (sound navigation and ranging).

In TOE, we commonly use sound waves with a frequency of 2.5 to 5 million cycles per second (megaHertz), or more. This is well above the upper range of normal human hearing, which in children is perhaps 20 kiloHz, and in senescent adults drops considerably below this.

What is the relationship between frequency and wavelength?

Why do we use such high-frequency sound in echocardiography? If you want to "see" a tiny object with sound, you will need sound that has short wavelengths. This is why bats use such high frequency sound. Now, within the heart, we are interested in resolutions of perhaps a millimetre or less, so it makes sense that we need short wavelengths. Consider a sound with a frequency of 5 MHz passing through water (or human tissue). Sound passes through soft tissues at almost exactly 1540 metres per second, so 5 Million cycles will fit into 1540 metres. From this we can work out the wavelength, which is simply the distance between one peak of sound intensity and the next - here 1540 metres / 5 000 000, or about 0.3 millimetres. You can also work out the general formula that relates velocity, frequency and wavelength:

v = f * lambda         .. Equation 1.1

Also note the simple relationship between the time interval from one crest to another (lets call this delta t) and frequency:

delta t = 1 / f         .. Equation 1.2

.. and the relationship between velocity and distance travelled (deltas):

delta s = v * delta t         .. Equation 1.3

Imagine a beam of ultrasound passing through a tissue. Part of the beam will be attenuated (absorbed by the tissue), part will be reflected, and part will be scattered. We call the reflected part specular echoes - these are the components that we are interested in, while the scattered beam merely serves to confuse us! Scatter is worse in tissues that are inhomogeneous at a microscopic level, while specular reflection occurs mainly at interfaces between tissues of different density. For example, a beam passing through the point where blood meets, say a papillary muscle, will be strongly reflected at this point, provided the size of the reflecting object is at least a quarter of the wavelength of the ultrasound.

There is a trade-off - we pick up smaller objects more easily with shorter wavelengths, but as we shorten the wavelength, scatter worsens, increasing our confusion! In addition, as we shorten the wavelength (raise the frequency) of an ultrasound beam, we impair its tissue penetration.

How is sound attenuated in various tissues?

Different tissues differ dramatically in their ability to transmit sound. We usually talk about the "half-power distance", which is the depth at which the sound is halved. Here are some examples:

Half-power distances
(2 MHz beam)
Tissue Distance (mm)
Water 3800
Blood 150
Most soft tissue 10 to 50
Muscle 6 to 10
Bone 2 to 7
Air, Lung under 1
Data from deBruijn & Clements, Intraoperative use of echocardiography, Ch 2 p 16.
JB Lippincott Co, 1991, ISBN 0 397 51128 0.

You can see that there will be minimal transmission of ultrasound through lung and bone, and that even muscle attenuates the beam dramatically. This is the main reason for the superiority of TOE over transthoracic echo - there is just too much attenation with TTE due to the presence of chest wall and especially, lung tissue.

Timing

A beam of ultrasound moves fairly rapidly through tissue - 1540 m/s is far from sluggish! If you wish to interrogate tissues to a depth of say 15cm, the beam can get there and back in 0.150/1540 seconds, that is, 195 microseconds. This gives us the opportunity to send out multiple pulses of sound every second, and obtain a dynamic picture of what is going on in the tissue. For example, we might send out just three pulses at a frequency of 5 MHz, and then listen for about 250 microseconds, mark time for 750 odd microseconds, and then send out another set of pulses, and so on.

If we display on a screen the echoes coming back from one set of impulses, the picture isn't very informative. The Y axis will tell us how deep a particular echo is, and more intense echoes will be brighter, but the single line will be difficult to interpret. But if we draw successive lines on the screen, as each packet of information comes in, we get a motion display, usually abbreviated to M-mode. From the M-mode image we can precisely determine depths and therefore dimensions (and changes in dimensions) of structures. It's often quite easy to interpret which structure we're looking at on M-mode, but of course our task is made easier by modern technology, where we usually look at a two-dimensional image, and then simply indicate where we want our M-mode 'slice'. How do we obtain two-dimensional images?

Two-dimensional Images

Early two-dimensional echo probes had a rotating head. This allowed beams of ultrasound to be fired off at various (continually changing) directions, and the returning echoes provided information that could be reconstructed into a two-dimensional picture of the tissue being interrogated.

Current imaging is more subtle. Probes consist of an array of parallel piezoelectric elements - often sixty four or more of these. We use trickery to produce a wave-front that moves in a certain direction - if all of the piezoelectric elements fire at the same time, then a wave front will be formed from all of these little elements, and it will move off parallel to the plane of these elements. But what happens if we fire first one element, then the next in line, and so on? If you sit down and draw this, then you'll soon see that the wave front will move at an angle to the line of the piezoelectric array. This trick can be used to move the angle of the wavefront through ninety degrees (practically, although the theoretical maximum is a full 180o), so multiple beams can be used to scan a ninety degree sector of tissue. Returning echoes are then re-constituted into a 2-D picture.

Because of the rapid return of the echoes, a full image can be captured about every sixtieth of a second, and images can then be displayed to give an impression of continuous motion (two images are usually interlaced, so the 'refresh rate' of the image you see on the screen is usually about thirty times per second).

What determines the resolution of a picture?

In a two-dimensional image, there must be two different resolutions, the "side-to-side" (lateral) resolution being determined mainly by how narrow the beam is (assuming the wavelength is appropriate), and the resolution along the length of the beam (axial resolution) by frequency and number of pulses (The higher the frequency, the better the resolution, the more pulses you send, the more 'things get confused' and the lower the resolution). The axial resolution is approximated by:

Raxial = 0.77 * Npulses / f

where R is resolution in millimetres, Npulses is the number of pulses, and f is the transducer frequency in MHz. For a 5 MHz transducer and 3 pulses, the resolution works out at about half a millimeter. Lateral resolution is far worse - generally about 2 mm.

How deeply can we look?

The major factor that limits the depth to which an ultrasound beam can interrogate a tissue (apart from attenuation) is related to beam divergence. Even the "flattest" wavefront that originates from a transducer tends to diverge with distance (To use an optical analogy, even the most tightly collimated laser beam is half a kilometer wide when it reaches the moon)! The "depth of field" (the distance at which the beam is no longer parallel) is estimated by:

L = (Ltransducer) / lambda

Where L is the depth of field, Ltransducer is the size of the transducer, and lambda is the wavelength. The larger the transducer, or the smaller the wavelength, the greater the depth of field.

Signal Manipulation and Image Processing

The raw data coming from an array of piezoelectric crystals at the end of a TOE probe needs extensive massaging before it can be recognised by even the most skilled eye. There are several points at which one can intervene:
  1. The frequency, intensity and number of impulses sent into the tissue can be varied. For example, the outgoing signal can be damped, with a shorter, weaker signal resulting in less noise despite the weaker returning signal.
  2. The returning signal can be modified before it is digitised (turned into a number that can be manipulated by digital computers). Such modification is called pre-processing.
  3. The digitised information can be extensively manipulated before it is displayed. This is post-processing.

Pre-processing

A variety of filters and compensatory adjustments can be made to the incoming signal. These include:
  • Dynamic range manipulation - low intensity noise is removed, and weaker signals are often boosted relative to more intense signals (for example, by using the logarithm of intensities).
  • Gain / attenuation - Altering gain (or attenuation) changes all signals proportionally. Setting the gain too high will enhance all signals to the point of saturation; too much attenuation will result in loss of detail.
  • Depth compensation ("time-gain compensation") boosts signals from deeper tissues (as they will have been more attenuated).
  • A variety of other electronic manipulations of the signal is possible, for example M-mode signals are often made more crisp using leading edge enhancement, where the signal is electronically differentiated (the first derivative is taken) to localise the timing (and thus depth) of an echo more precisely.

Post-processing

The digitised image obtained from even the most sophisticated transducer is actually rather unattractive, and full of holes. Extensive post-processing is used to render the picture more attractive and intelligible. The radial information is transformed into a raster image (rectangular image made up of dots), commonly made up of about 65 000 or more picture elements (pixels). The image is grey-scale (varying shades of grey), with up to 1024 shades of grey stored in the memory of the echo machine, although fewer shades are usually shown on the screen! Other tricks include:
  • Temporal averaging can be used to smooth the current data with previous data (variable persistence) or decrease noise (by taking the average of several frames).
  • "Histogram equalisation" - each picture element is adjusted up or down according to the average intensity of the whole image.
  • Many other filters can be applied to 'enhance' the image - smoothing images to eliminate noise, enhancing edges, tracing boundaries, and so on.

The principles of Doppler echocardiography

The basic idea behind the Doppler effect is simple - as usual, the devil is in the details. We all know how an ambulance siren appears more high pitched as the vehicle moves towards us, and then wheee..ooooh becomes more bass as the flashing red light vanishes into the distance. Clearly if an object that is stationary relative to you is emitting sound at a frequency f, the interval between pulses will be 1/f seconds. But if the object is moving directly towards you and emits a pulse, it will have caught up slightly with this pulse at the time it emits the next pulse, so the frequency of the sound will be boosted. Similarly, if an object moving towards you reflects a sound you have emitted, it will raise the frequency (pitch) of the sound it reflects. Let's look at this in a little more detail.. (or you may wish to skip straight to the Doppler equation)!

Consider a doppler probe emitting a continuous stream of ultrasound towards a column of blood heading directly towards the probe. The oncoming red blood cells will reflect some of the sound back towards the probe. Let's find out how we can calculate the velocity of the blood flow, knowing the change in frequency of the reflected beam.

Some basics

First of all, we need to set some ground rules. Don't be intimidated by the symbols - they're all pretty well self-explanatory. We will often refer to the following, so let's make them quite clear:

  1. Velocities:
    • vcell - the velocity of an oncoming red cell;
    • vsound - the velocity of sound in water (or blood);

  2. Distances and times:
    • delta t - a time interval
    • delta s - a distance travelled (in that time interval)

  3. Wavelengths and frequencies:
    • lambdaout - the wavelength of the beam emitted by the probe;
    • lambdar - the wavelength of the reflected beam
    • Likewise, fout and fr are simply the frequencies of the two beams.

  4. Simple equations (rules -we've already introduced these above):

    • v = f * lambda         .. Equation 1.1

    • delta t = 1 / f         .. Equation 1.2

    • delta s = v * delta t         .. Equation 1.3

Sorting out the Doppler Equation

Now let's attack the Doppler equation. Our goal is to find out the velocity of a column of blood moving towards a Doppler probe, if we know the change in frequency of the reflected ultrasound beam. Consider the following picture:

animation of doppler effect - see text

The stationary black block in the above diagram simply reflects the sound, and the frequency of the reflected sound is the same as that of the incoming beam, that is, fout. The moving magenta block (or red cell, if you wish) is quite different. At some point, the moving cell will encounter a wave crest, and reflect the wave. The cell will then carry on moving, and therefore will encounter the next wave crest sooner than if it were stationary. Can we calculate how long it will take for the moving cell to meet the next wave crest? Certainly..

In the time that it took the cell to meet the next wave crest, the cell has moved towards the wave crest with velocity vcell, and the wave crest has moved towards the cell with velocity vsound. The total distance the two have moved before they collide is clearly the distance between oncoming two wave crests, which is lambdaout. In other words:

  deltascell + deltaswavefront = lambdaout         .. Equation 2.1

Now, we know that delta s = v * delta t , so we can substitute in Equation 2.1 to get:

  vcell * delta t + vsound * delta t = lambdaout         .. Equation 2.2

Technically, the velocities vcell and vsound are actually vectors, and we are only interested in their magnitudes, so we should write the above equation as:
|vcell| * delta t + |vsound| * delta t = lambdaout
If we're using vectors, we should subtract the second from the first as they're moving in opposite directions .. but in all of the following, for clarity we will assume that we are talking about the magnitude |v| where we simply write "v". Also carefully note the assumption under which we are working - that the cell is moving directly towards the wavefront coming from the ultrasound probe.

Rearranging, we find that delta t, the time taken for the moving cell to meet the next wave crest, is simply:

  delta t = lambdaout / (vcell + vsound)         .. Equation 2.3

Think carefully now..

As we've mentioned, the red cell reflects first one and then a second crest. Consider the distance between the first reflected crest and the second one. In the time (delta t) it takes for the cell to encounter the second pulse, the first pulse has moved a certain distance back towards the ultrasound probe. But wait a bit! Because the cell is moving in the same direction as the reflected sound, the cell partially catches up with the beam it reflected! How do we take this into account? Easy. At the instant we are considering, the distance between the cell and the reflected wave (which is the same as the wavelength of the reflected beam) is simply the distance travelled by this wave front minus the distance travelled by the red cell in the same time. Knowing that distance equals velocity times delta t, we see that:

  lambdar = ( vsound * delta t ) - ( vcell * delta t )         .. Equation 3.1

In other words, we need to subtract the distance travelled by the red cell in the same direction as the reflected crest. We can now easily work out the frequency of the reflected sound. First, we rearrange Equation 3.1, thus:

  lambdar = delta t * (vsound - vcell)         .. Equation 3.2

Then we substitute the value we obtained for deltat from Equation 2.3 into Equation 3.2, so:

  lambdar = lambdaout * (vsound - vcell) / (vcell + vsound)         .. Equation 4.1

Using our knowledge of the relationship between frequency and wavelength, we can easily rearrange this to:

  vsound / fr = vsound / fout * (vsound - vcell) / (vcell + vsound)         .. Equation 4.2

Cancelling the vsound and rearranging, we get a formula for the frequency of the reflected beam..

  fr = fout * (vsound + vcell) / (vsound - vcell)         .. Equation 4.3

Look at the above equation

Does Equation 4.3 make sense? Yes. The reflected frequency will be boosted, as we expect. (The police siren is more shrill as it approaches us). Conversely, if the magnitude of vcell is negative, the frequency will be lowered. It's now easy to rearrange our equation to solve for vcell:

  fr * (vsound - vcell) = fout * (vsound + vcell)

  vcell = vsound (fr - fout) / (fout + fr)

If we abbreviate (fr - fout) to equal " delta f ", then we can state our equation even more succinctly..

  vcell = vsound * delta f / (2 * fout + delta f)         .. Equation 5.1

When we're dealing with physiological blood flow, the change in frequency is small compared to fout (This is because the velocity of blood flow is maximally a few metres per second, which is trivial compared with the 1540 m/s speed of sound in water). So we conveniently ignore the delta f in the denominator, to get our penultimate equation:

  vcell approx equal vsound * delta f / (2 * fout)         .. Equation 5.2

An important complication

Note that in the above we assumed that the column of blood was headed directly towards the Doppler probe. In real life, this is by no means always the case. We need to be very careful, because gross errors (underestimates) can occur if we don't realise that the above principles only apply where blood flow is in the line of the probe - either towards or away from the probe. The Doppler shift of blood flowing at right angles to a probe is always zero! Fortunately for us, it's easy to compensate for small errors in probe angulation. Consider blood flowing not directly towards a probe, but at an angle theta to the orientation of the probe:

blood flow at an angle to the probe..

Only the component of the blood flow directed towards the probe will be measured. If the true vector is away from the direction of the ultrasound beam, then:

cos theta = |observed vector| / |true vector|

We thus compensate by saying:

The above is the 'true' Doppler equation for sound reflected from a moving object, where vcell is the velocity of the object moving towards the probe, fout is the sound frequency emitted by the probe, vsound is the velocity of sound in water (1540ms-1), and theta is the angle between the ultrasound beam and the direction of motion of the moving object. delta f is the change in frequency.


Applications

Here we plan to have a section on actual doppler images, and Doppler assessment. It will be a long time before we can achieve this goal!

< Previous ^ Home Next >