Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

акустика / gan_ws_acoustical_imaging_techniques_and_applications_for_en

.pdf
Скачиваний:
83
Добавлен:
04.05.2023
Размер:
5.27 Mб
Скачать

130

 

Acoustical Imaging: Techniques and Applications for Engineers

0 mm

 

 

 

 

 

 

 

1

 

 

df

1

 

 

 

00

 

H1(f )

00

 

T1(f )

80

160

 

80

160

1

 

 

 

1

 

 

T2(f )

00

 

H2(f )

 

00

 

80

160

 

80

160

dt

 

 

 

1

 

 

 

1

 

H3(f )

 

 

 

T3(f )

00

 

 

00

 

80

160

 

80

160

1

 

 

 

1

 

 

 

 

 

H4(f )

 

 

 

T4(f )

00

80

 

0

 

 

160

 

0

80

160

1.6 mm

 

 

 

 

 

 

 

 

 

 

 

Figure 7.5 Attenuation compensation: (left) transfer functions of skin in four different depths; (right) depth-adapted pseudo-inversely prefiltered transmitter signals (Passmann and Ermert [5] © IEEE)

amplitudes have to be raised), the noise is high. This effect can be reduced by applying prefiltered transmitter signals Ti ( f ) (Figure 7.5) instead of a postfiltering procedure, which can be accomplished with no additional effort using the B/D-scanning procedure, as independent scans for different depths/areas are carried out in any case. Even with the prefiltering procedure, the SNR is degraded, because the available transducer frequency band is not optimally utilized. Therefore, a technique that increases the energy of the signal must be employed. Hence we use the well-known pulse compression technique, in combination with nonlinear frequencymodulated (NLFM) chirp signals, for that purpose. The basic principle of pulse compression is shown in Figure 7.6. Here, instead of a short pulse (top) with lower and upper bandwidth limits f1 and f2, respectively, a frequency-modulated (FM) chirp signal x(t ) (bottom) is transmitted with the same bandwidth f = f2 f1 as the pulse:

 

 

=

 

 

 

T

 

 

 

 

 

 

 

 

x (t )

 

rect

 

2t T

cos[ϕ (t )]

 

(7.3)

 

 

 

 

 

 

 

=

 

1

 

+

2

 

 

=

f2

T

 

 

d (t )

 

2π

f

t

 

1

kt2 k

 

f1

(7.4)

 

 

 

 

 

 

 

Signal

Spectrum

System

 

 

Echo

 

 

Result

H(ω, df)

H(ω, df)

Filter

Figure 7.6 Pulse compression: basic principle (Passmann and Ermert [5] © IEEE)

High-Frequencies Acoustical Imaging

 

 

 

 

 

131

The instantaneous frequency f (t ) can be found by:

 

 

f (t ) =

1 dϕ (t )

= f1 + kt

(7.5)

 

 

 

 

2π dt

Due to the larger signal duration T of the chirp, the signal energy is increased. The received echo is compressed by a digital (or analogue) all-pass filter HF ( f ) which depends on both the spectrum X ( f ) of the transmitter signal and the system transfer function H( f ):

HF ( f ) = ej(arg(X ( f ))+arg(H( f )))

(7.6)

yielding the output signal y(t ) by an inverse Fourier transformation:

 

y (t ) = F−1 {X ( f ) H ( f ) HF ( f )}

 

= F−1{|X ( f )| |H ( f )|}

(7.7)

Prefiltering, using the pulse compression method, can be done in different ways, as shown in Figure 7.7. If a particular spectrum X ( f ) of the chirp signal is desired for a certain depth, it can be generated by a linear frequency-modulated (LFM) chirp, which is also amplitudemodulated, or by a nonlinear-modulated (NLFM) chirp with a constant amplitude [7]. The latter transmits the higher frequency spectral components for a longer period, thus utilizing the full amplifier power for the entire transmission time. Therefore, the signal energy of the prefiltered NLFM chirp is higher than that of the LFM chirp. A relation between the group delay time Td ( f ) and the desired spectrum X ( f ) can be found in reference [8]:

dTd ( f )

= C |X ( f )|2

(7.8)

d f

 

 

Solving (7.8) by integration:

 

 

 

 

 

 

f

 

 

 

 

Td ( f ) =

C1 |X (ν )|2 dν + C2

(7.9)

 

0

 

 

 

 

The constants C1 and C2 can be determined using the condition:

Td ( f1 ) = 0 Td ( f2 ) = T

(7.10)

where T is the duration of the chirp signal.

Signal

Spectrum

System

Filter

Result

H(ω, df) Filter

LFM

H(ω, df) Filter

NLFM

Figure 7.7 Pulse compression: prefiltering (Passmann and Ermert [5] © IEEE)

132

Acoustical Imaging: Techniques and Applications for Engineers

1

0.5

0

 

 

s/μm

0

200

400

Figure 7.8 Pulse compressed echo of a strong scatterer (Passmann and Ermert [5] © IEEE)

The inverse function of Td ( f ) yields the instantaneous frequency f (t ):

f (t ) = inv {Td ( f )} =

1 dϕ (t )

(7.11)

 

 

 

2π dt

and

 

 

 

 

 

ϕ (t ) = 2π

t

 

 

 

 

inv {Td ( f )} (τ ) dτ

(7.12)

0

Equations (7.3), (7.9) and (7.12) yield an expression for a nonlinear frequency-modulated chirp signal with the desired spectrum X ( f ):

x (t )

 

rect

2t T

cos

2π

 

 

t

inv

f

C1

 

X (ν )

2 dν

 

C2

 

(τ ) dτ

(7.13)

=

·

 

 

|

+

T

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

Using NLFM, signals with a duration of 0.5 μs and a bandwidth of 100 MHz yield a gain of 12 dB.

As the pulse compression method assumes a linear imaging system, a problem will arise in tissue regions that have large inhomogeneities or strong scattering as these cause echoes that interfere with the dynamic range of the imaging system. In this situation, the system behaves nonlinearly and, consequently, large side lobes are generated during the compression procedure as shown in Figure 7.8. For comparison, Figure 7.9 shows the compression result of the behaviour in a linear system.

To solve this problem, an adaptive imaging procedure is proposed:

1.Acquisition of a pulse and a chirp image simultaneously.

2.Detection of regions of strong backscattering in the raw data set.

1

0.5

0

200

400

s/μm

0

Figure 7.9 Pulse compressed echo of a weak scatterer (Passmann and Ermert [5] © IEEE)

High-Frequencies Acoustical Imaging

133

Figure 7.10 Chirp-mode skin image – 10 mm × 1.6 mm (Passmann and Ermert [5] © IEEE)

3.Compression of the chirp signal data.

4.Replacement of chirp image data by pulse image data at all detected locations of strong backscatterer and within a certain region around them, which is determined by a measurement considering the worst case nonlinear behaviour.

To further improve the detectability of small regions, an adaptive filter for speckle suppression, presented by Loupas et al. [9], can be used. This filter has to be modified for high-frequency ultrasound applications by an adaptive estimation of the resolution cell size.

A chirp and a pulse image, obtained by Passmann and Ermert [5], is shown in Figures 7.10 and 7.11. The chirp image shows homogeneous resolution over the full depth area but also large side lobes can be observed in the skin-entry-echo region. The pulse image is free of artefacts but shows degrading resolution in the lower image regions. Figures 7.12 and

Figure 7.11 Pulse-mode skin image – 10 mm × 1.6 mm (Passmann and Ermert [5] © IEEE)

134

Acoustical Imaging: Techniques and Applications for Engineers

Figure 7.12 Combined skin image – 10 mm × 1.6 mm (Passmann and Ermert [5] © IEEE)

7.13 are also due to Passmann and Ermert [5]. Figure 7.11 shows the combined image and Figure 7.13 shows a speckle reduced version of it. The combined image shows homogeneous resolution over the full-depth area and is free of artefacts, therefore the diagnostic value of the image is considerably improved. The speckle suppression allows easier detection of image details. Hence, high-frequency ultrasound shows the potential for high-resolution skin and eye imaging in vivo, especially for glaucoma research.

The strongly focused transducers employing the B/D-scan conception yield more isotropic images. Attenuation compensation increases the image homogeneity and pulse compression provides the signal energy required for attenuation compensation. An adaptive combination of pulse and chirp mode avoids nonlinearity artefacts, and adaptive speckle processing allows for easier detection of image details.

Figure 7.13 Speckle reduced combined image – 10 mm × 1.6 mm (Passmann and Ermert [5] © IEEE)

High-Frequencies Acoustical Imaging

135

7.8Electronic Circuits of Acoustical Microscope

7.8.1Gated Signal and Its Use in Acoustical Microscope

As the theory of the acoustical microscope [10, 11] is built on a monochromatic or single frequency concept, it is therefore necessary to understand the relationship between the frequency domain and the time domain. An acoustic signal may be described in terms of the components of each frequency present, known as the frequency domain, or in terms of the time domain, which shows the signal being oscillatory, as seen on the oscilloscope. These two domains can be related by a Fourier transform as

 

 

F ( f ) =

f (t ) ei2π f t dt

(7.14)

−∞

From the properties of the Fourier transform, it follows that if a pulse is of finite duration, then it cannot have a single frequency but rather consist of a finite spread or spectrum of frequency components. An ideally monochromatic signal would have only a single frequency present, and therefore would have to be of infinite duration. For instance, if the envelope of a pulse is of rectangular form, then in the frequency domain it will have the Fourier transform of a rectangular function (or a sinc function, defined as sinc(x) = sinx/x). For a rectangular pulse of length to, the width of the central maximum in the frequency spectrum is 2/t0. This gives a measure of the spread of frequencies present.

The convolution of two functions of the same variable, G( f ) and H( f ) can be written as

G f H f f d f

 

F ( f ) =

(7.15)

−∞

The convolution operation is a way of describing the product of two overlapping functions, integrated over the whole of their overlap for a given value of their relative displacement. The convolution theorem states that the Fourier transform of the product of two functions is equal to the convolution of their separate transforms

F {G ( f ) × H ( f )} = F {G ( f )} F{H ( f )}

(7.16)

where F denotes the Fourier transform and denotes convolution.

If we take a pulse as the output of a monochromatic oscillator of amplitude A0 and apply to it a rectangular gate of width t0, the resulting spectrum will be the convolution of the spectrum

of the oscillator (a delta function f0) and the spectrum of the gate, producing

 

F ( f )

=

A0t0 sin {π ( f f0 ) t0}

(7.17)

π ( f f0 ) t0

 

 

If the pulses are repetitive, with a repetition frequency f1, then the gated signal is convolved in the time domain with a comb function (a comb function is a series of delta functions at constant spacing). This means that the Fourier transform of a single-gated pulse is multiplied by the Fourier transform of the pulse repetition which is, in turn, a comb function starting at the origin and repeating at intervals of f1. Thus, the frequency spectrum of a series of gated

136

Acoustical Imaging: Techniques and Applications for Engineers

pulses is the Fourier transform of the envelope of an individual pulse, centred at the frequency of the oscillator, and composed not of a continuous distribution but of a series of discrete spikes (that is, delta functions) within that envelope.

In the acoustic microscope, the required signal can be selected not only in the time domain but also in the frequency domain. It is then possible to select the specimen echo and separate it in time from the unwanted lens echo over the basis of the design consideration of the focal length of the lens, and hence the resolution is available. Smaller lens reverberations are always present that cannot be separated in time from the specimen’s echo. Moreover, even if such a reverberation could be eliminated, there could still be a problem from the fact that different parts of the acoustic wave from the specimen may be reflected with different time shifts. This could occur, if the specimen has two closely spaced layers, each of which reflects some energy, or if some of the energy is reflected back after coupling into surface waves. If the signal is measured by a peak detector circuit, a reasonable value will be obtained when the interference is constructive, but when it is destructive the value will be quite wrong, corresponding simply to the largest component present. If very long duration pulses could be used, it might be possible to measure the middle where the overlap is adequate, but this is not practicable when there are constraints on the pulse length – as there are in all high-frequency microscopes.

Acoustic microscopy involves the interference of an incident wave, a reflected wave and a refracted wave and the interference theory is built on the assumption that the microscope is perfectly monochromatic. This is especially so for the theory of the V (z) technique, studying how the contrast varies with defocus. The monochromatic theory could be summed over a frequency spectrum that is actually used in a given system, but would be almost impossible to invert when interpreting measured results. For much of qualitative imaging the narrow band detected signal may be adequate.

For more accurate measurements, the heterodyne circuit is used. An example of a heterodyne circuit is shown in Figure 7.14 [15].

If the bandwidth of the detector system can be narrowed to a fraction of the spacing between adjacent spikes, it is possible to select one of them and to reject all other signals at spacings of multiples of the pulse repetition frequency. In this way, it is possible to isolate and measure a single-frequency component and then achieve a measurement that is monochromatic, subject only to the phase noise of the signal source. Here the local oscillator is phase-locked to the signal source and has a frequency that differs from the signal source by some small amount [12, 13].

In the simple heterodyne circuit shown in Figure 7.14, and is used in most imaging microscopes, the pulse length is defined by the switch S1. As the speed of this switch determines the minimum pulse length and, hence, the minimum lens focal length, it also determines the highest frequency of the microscope. Thus the limit to the resolution is ultimately determined by the highest useful speed of the switch. S2 is the single-pole double-throw SPDT) switch and A1 the low-noise preamplifier. The mixer is a device that takes two inputs: the radio frequency signal and the local oscillator signal and gives outputs as the sum and different frequencies of these two [14]. The local oscillator frequency is chosen so that the different frequency is in the middle of the pass band of the narrow band filter, which is placed in the circuit at the intermediate frequency output of the mixer. Gating of the preamplifier output is achieved by gating the local oscillator of the heterodyne circuit.

High-Frequencies Acoustical Imaging

 

 

 

 

 

137

Signal

Power

 

 

 

 

i.f.

Detector

Video

Peak

generator

amp

S1

S2

Mixer

Filter

amp.

diode

amp.

def.

 

 

 

A1

R L I

 

 

 

 

 

Acoustic

Gated local

oscillator

 

(a)lens

40 kHz

Frequency

53 kHz

Lock-in

 

 

 

 

shifter

Reference

amplifier

Signal

 

 

 

 

 

 

 

 

 

 

Output

 

 

225.000 MHz

 

93 kHz

 

 

 

 

 

L.O. isolation

L.O.

 

 

L.O.

 

 

 

gate

 

 

 

 

 

 

 

 

 

T

α

α

 

 

 

Clock SYNC

L I

 

 

 

 

 

 

R

 

 

 

 

 

40 kHz

 

 

SPDT

RL I

 

 

to gate

 

 

 

 

T

 

T

 

 

R.F. osc.

 

r.f.

 

 

100 kHz

53 kHz

 

gate

 

 

225.093 MHz

 

 

 

L.P.

B.P.

 

 

 

 

Signal Acoustic monitor

lens

(b)

Figure 7.14 Schematic radio frequency (r.f.) systems: (a) simple heterodyne circuit, S1 determines the pulse length, S2 switches the lens from transmit to receive, and A1 amplifies the reflected signal;

(b) quasi-monochromatic circuit; the two oscillators and the pulse repetition frequency are phase-locked, and the final signal is lock-in detected (Weaver [15])

7.8.2Quasi-Monochromatic Systems

For quantitative acoustical microscopy, that is, the V (z) curve measurement, the quasimonochromatic circuit of Figure 7.14(b) is more accurate [15].

The basic principles of this circuit are similar to those of Figure 7.14(a) but there are some important differences. The radio frequency (r.f.) oscillator and the local oscillator are two frequency synthesizers that are phase-locked to one another to give a difference frequency of precisely defined phase. This difference frequency is very much lower than in the simple heterodyne circuit in order to make it comparable with the pulse repetition frequency.

In order for the system to behave as truly monochromatic, it is important that the receiver gate should not introduce any distortion of the frequency spectrum of the reflected pulse. The

138

Acoustical Imaging: Techniques and Applications for Engineers

reflected pulse will be multiplied in the time domain by the receiver gate, and that corresponds in the frequency domain to a convolution with the (repetitive) pulse gate transform.

After the output of the mixer has been amplified and low-pass filtered, it may be fed to a lock-in amplifier. This is a phase-sensitive detector performs at near-audio frequencies, rather as a mixer does at radio frequencies but with more sophisticated controls.

With a dual-phase lock-in amplifier, the in-phase and quadrature components, or equivalently the module and phase, can be measured simultaneously. If only the module of the detected component is required, then a readily available alternative to the circuit illustrated in Figure 7.14(b) is to use a spectrum analyzer with a measurement output in place of the heterodyne detector system in Figure 7.14(a), with its own tracking generator as the r.f. oscillator and the swept frequency range set to zero [16].

The quasi-monochromatic system has great advantages in terms of SNR, dynamic range and linearity over the quasi-continuous wave (qcw) system. It also allows the effect of the troublesome lens reverberation that occurs within the time gate to be at last eliminated. Since this also appears now as a monochromatic contribution, it can be removed by subtracting a uniform signal of compensating amplitude and phase at the reference frequency from the input to the lock-in amplifier. But most important of all, it allows measurements to be made at, effectively, a single frequency, so that they can be analyzed accurately in terms of the monochromatic theory of V (z) (the video signal as an explicit function of defocus).

7.8.3Very Short Pulse Technique

For some purposes, it may be necessary to have accurate frequency definition; for others, a good time discrimination is useful. Because of the Fourier relationship between frequency and time, the more precisely the time of a signal is known, the greater is it necessary to know the bandwidth of frequencies. Approximately, the time resolution τ is the reciprocal of the

bandwidth Bw that their product Bwτ = 1.

The design of a system for working with short pulses follows the same principles as other pulsed ultrasonic systems, such as ultrasound flaw detectors, but it is necessary to achieve very much great stability and high bandwidth.

In Sections 8.3, 9.2 and 9.3, the techniques and applications of time-resolved quantitative measurements will be given in which impulse excitation is essential. But many uses of acoustic microscope call primarily for an image, and for that purpose conventional gated continuous wave electronics are quite adequate and indeed are to be preferred because of the good signal strength that they give for obtaining pictures at a reasonable speed with good contrast.

References

[1]Korpel, A., Kessler, L.W. and Palermo, P.R. (1871) Acoustic microscope operating at 100 MHz. Nature, 232(5306), 110–111.

[2]Hoss,¯ A., Ermert, H., el Gammal, S. and Altmeyer, P. (1989) A 50 MHz ultrasonic imaging system for dermatologic application. IEEE Ultrasonics Symposium Proceedings, pp. 849–852.

[3]Crocker, R.L., Gray, N.J., Phillips, R.B. and Chivers, R.C. (1989) High frequency ultrasonic B-scanner.

Ultrasonics International 89 Conference Proceedings, pp. 578–583.

[4]Nakamura, M.N., Obara, K. and Segawa, M. (1976) The strong piezoelectricity in polyvinylidene fluoride (pvdf). Ultrasonics, 14, 15–23.

High-Frequencies Acoustical Imaging

139

[5]Passmann, C. and Ermert, H. (1995) Adaptive 150 MHz ultrasound imaging of the skin and the eye using an optimal combination of short pulse mode and pulse compression mode. Proceedings IEEE Ultrasonics Symposium, Seattle, USA, 1995.

[6]Passmann, C. and Ermert, H. (1995) In vivo imaging of the skin in the 100 MHz region using the synthetic aperture concept. Proceedings IEEE Ultrasonics Symposium, Seattle, 1995.

[7]Pollakowski, M. (1993) Ein Beitrag zur Anwendung der Pulskom-pressionstechnik in der Zerstorungsfreien¨ Werkstoffprufung¨ mit Ultraschall. Ph.D. thesis, Ruhr-Universitat Bochum, Shaker, Aachen.

[8]Fowle, E.N. (1961) A general method for controlling the time and frequency envelopes of FM signals. Technical Report 41G-0008, Massachusetts Institute of Technology (MIT), Lincoln Laboratory, Lexington, Mass., June 1961.

[9]Loupas, T., McDicken, W.N. and Allan, P.L. (1989) An adaptive weighted median filter for speckle suppression in medical ultrasonic images. IEEE Trans. Ultrason. Ferroel. Freq. Control, UFFC-36(1): 129–135.

[10]Briggs, A. (1992) Acoustic Microscopy, Clarendon Press, Oxford, pp. 64–78.

[11]Atalar, A. and Hoppe, M. (1986) High-performance acoustic microscope. Rev. Sci. Instrum., 57, 2568–2576.

[12]Liang, K.K., Bennett, S.D., Khuri-Yakub, B.T. and Kino, G.S. (1985) Precise phase measurements with the acoustic microscope. IEEE Trans., SU-32, 266–273.

[13]Liang, K.K., Bennett, S.D. and Kino, G.S. (1986) Precise phase measurements with short tone burst signals in acoustic microscopy. Rev. Sci. Instrum., 57, 446–452.

[14]Henderson, B.C. (1990) Mixers in microwave systems. Watkins-Johnson Company Tech-notes, 17(1), 1–15; 17(2), 1–13.

[15]Weaver, J.M.R. (1991) An optimal R.F. system for quantitative acoustic microscopy. IEEE Trans. Ultrason. Ferroel. Freq. Control. UFFC.

[16]Kushibiki, J. and Chubachi, N. (1985) Material characterization by line-focus-beam acoustic microscope. IEEE Trans., SU-32, 189–212.