Medical diagnostic systems (Orvosbiolσgiai kιpalkotσ rendszerek) Beamforming in ultrasound ( Nyalαbalkotαs az ultrahangban) Miklσs Gyφngy • Array-dependent scanning methods (linear, phased, 2-D,...) • Beamforming strategies (fixed & dynamic focus, advanced) • Speckle reduction techniques (compounding) • Sidelobes and sidelobe reduction techniques (apodization) • Typically higher frequencies • example: 5-10 MHz • high resolution (~0.2 mm) • short penetration depth (~10 cm) • Subset of elements (subaperture) used to form each A-line • Good for imaging organs with easy access (e.g. abdominal organs) • Elevational (in-plane) focusing achieved using acoustic lens • typical resolution 10 mm • Spacing between elements on the order of a wavelength (see “grating lobes” slides later on) entire aperture used to generate A-line “Nerve movement in forearm during wrist extension” http://images.wellcome.ac.uk/ B0004357 Copyright work under Creative Commons licence generating A-line • Typically lower frequency • example: 1-4 MHz • low resolution (~0.6 mm) A-line • large penetration depth (~30 cm) • Good for imaging deep, hard-to­access (limited acoustic window e.g. due to ribs) organs (e.g. heart) • Elevational (in-plane) focusing achieved using acoustic lens • typical resolution 14 mm ? • Spacing between elements less than half a wavelength (see “grating lobes” slides later on) “Ultrasound image of normal 24 week fetus” http://images.wellcome.ac.uk N0019385 Copyright work under Creative Commons licence Medical diagnostic systems – Beamforming in ultrasound Arrays for 3-D imaging Freehand 1D array (position feedback with e.g. optical markers) – difficult registration (may be aided by position sensing) + simple to use Mechanised 1D array (fixed, predictable motion e.g. inside casing) – inflexible + simple registration 2D array – electronic complexity – element spacing + real-time 3D x7-2 array from Philips with 50Χ50 elements Medical diagnostic systems – Beamforming in ultrasound 3-D visualisation: surface vs. volume rendering Courtesy of Zonare Medical Systems 3D US image of left atrium [Agricola et al. 2010] http://www.zonare.com/products/clinical-images/id_7/ Creative Commons Attribution 3.0 Licence http://www.pagepress.org/journals/index.php/hi/article/viewArticle/h i.2010.e6/2133 Other types of array • Annular (concentric rings) [Anderson 2006] • accurate focusing along axis to produce A-line • need to be moved to provide B-mode • Element(s) inside catheter [Cobbold pp. 580-593] • single element moving inside catheter OR ring of elements used as phase array • exciting applications for imaging inside vessels, e.g. intravascular ultrasound (IVUS) For now, concentrate discussion of beamforming on 1-D arrays (linear and phased) Medical diagnostic systems – Beamforming in ultrasound Beamforming strategies • Fixed delay (near-field, far-field) • Dynamic focusing: vary receive beam with focus • Focal zone splicing: several transmission depths • Parallel receive beamforming (access to pre-beamformed data) • Synthetic aperture imaging (access to pre-beamformed data) Medical diagnostic systems – Beamforming in ultrasound Fixed delay beamforming Example: two transceivers A, B; point scatterer S • Estimation of scatterer strength becomes minimum-variance (electrical noise, “spatial noise” from other scatterers) combination of two transmissions and two receptions (cf. beamforming as spatial filtering [Van Veen and Buckley 1988]) • Transmit: delay B by 1 ΅s yA(t)= w(t); yB(t)= w(t-1) 4 ΅s • Signals received in phase at 5 ΅s • Receive: delay A by 1 ΅s (relative to B); sum b(t)= yA(t-5) + yB(t-4) • Delay-and-sum beamforming on receive Medical diagnostic systems – Beamforming in ultrasound Fixed delay beamforming • Simplest beamforming method • Same focus for transmit and receive .transmit_delays = -.receive_delays • Sharp focusing at target depth, but blurring at other depths • Far-field limit • plane wave • plane wave • linear variation of delays elements Example: D=10 mm, c=1500 m/s, f = 1.5 MHz. Far-field from D2/4. = 25 mm Fixed delay-and-sum beamforming on Tx, Rx. Adapted from [Burns 2005] Medical diagnostic systems – Beamforming in ultrasound Phased array beamforming (far-field) • Reference time (t=0) when pulse is transmitted from central element d • Usually, number of elements is even... • Fix relative delays, then scale 0 time to distance accordingly • Delay given by .(i,r)=2r/c+ id sin. id Medical diagnostic systems – Beamforming in ultrasound Dynamic receive beamforming • Once electric pulse converted into acoustic pulse (transmission), no user control over pulse propagation (single beam is formed) • Beamforming on receive, however, occurs d electronically/digitally and any number of beams can be synthesised 0 • Observation: as echoes return, they come from deeper and deeper objects • Idea: dynamically vary receive focus with time id • Corresponds to “stretching” (frequency modulation) of signal for out-of-centre elements • Delay given by .(i,r)=2z/c + (id)2/8z Medical diagnostic systems – Beamforming in ultrasound Focal zone splicing image • Well-focussed transmit beam leads spliced from all to good local sharpness Tx Tx1 Tx2 Tx3 • Transmit several pulses with different foci • Splice resulting images together depth • Reduction in frame rate • See [Szabo2004, p.309] for images adapted from [B-K Medical 2003, p. 59] Medical diagnostic systems – Beamforming in ultrasound Parallel receive beamforming: zone focussing • Frame rate depends on number of transmissions • With parallel digitization of element signals, several receive beams can be synthesised for a single, broader, transmit beam • Zonare Medical Systems (www.zonare.com) uses this technology Medical diagnostic systems – Beamforming in ultrasound Parallel receive beamforming: multiple beam transmission • Several transmit beams synthesised in one transmission • Parallel acquisition allows receive focussing (and thus separation) of the beams • Again, increase in frame rate results simultaneous synthesis of two beams adapted from [Cobbold 2007, p.476] Medical diagnostic systems – Beamforming in ultrasound Synthetic aperture imaging [Jensen et al. 2007] • Each element transmits a pulse on its own, one after the other • Echoes recorded by all (or many elements) at once • Assuming linearity, principle of superposition applies • Both transmit and receive foci can be synthesised retrospectively! • Higher image resolution • See [Jensen et al. 2007] for images Medical diagnostic systems – Beamforming in ultrasound Speckle reduction techniques • Speckle arises from interference between sub-wavelength scatterers from within one resolution cell • Spatial averaging (blurring) reduces speckle, but also reduces spatial resolution • Compounding • generate several images with different parameters • speckle hopefully weakly correlated • sum images to reduce speckle “noise” • averaging frames (temporal) causes blurring • look at two popular methods: spatial; frequency Medical diagnostic systems – Beamforming in ultrasound Spatial/angular compounding [Cobbold 2007, p. 469] • Subject imaged from several orientations • Scatterers interfere with different phases when angle of insonation/reception is varied • speckle weakly correlated • Summation of registered images from several orientations reduces speckle • Orientations generated from consecutive transmissions: • array need not move • frames separated by less than 100 ms • simple registration • real-time compound images Medical diagnostic systems – Beamforming in ultrasound Frequency compounding • Scattering directivity changes with frequency, but not origins of scattering • Therefore, strong correlation between images • Modest speckle reduction 6 MHz 8.5 MHz compound (“C8”) Images of tissue-mimciking agar gel wih vessel-mimicking inclusion; acquired on a z.one ultrasound system (Zonare) Medical diagnostic systems – Beamforming in ultrasound Sidelobes and sidelobe reduction techniques Sidelobes arise due to • limited aperture • sampled aperture (grating lobes) • temporal quantization errors (quantization error grating lobes) Sidelobe reduction techniques Pressure field against angle created by transducer. (Note that purely angular dependence implies far-field) Medical diagnostic systems – Beamforming in ultrasound Grating lobes [Szabo 2004, pp. 182-185] • Next lecture: in continuous wave mode, pressure in focal plane is 2-D Fourier transform of pressure amplitude distribution over aperture (i.e. apodization) • More precisely, angular distribution of pressure is 2-D Fourier transform of apodozation in far-field, and focusing brings the far-field to the focal plane • Thus (in analogy with the Fourier sampling theorem), if the discrete element spacing “samples” at d.0.5., grating lobes (“aliasing”) appear • Grating lobes first appear at .=±. (d=0.5.) and move closer to region of interest as d is further increased (new grating lobes appear every time d increased by 0.5.) • Are grating lobes avoided in practice? Medical diagnostic systems – Beamforming in ultrasound Grating lobes [Szabo 2004, pp. 182-185] • Are grating lobes avoided in practice? Two examples: • L10-5 linear array: fc = 7.5 MHz, N=128, D=38 mm • d = 1.5 . • P4-1 phased array: fc = 2.5 MHz, N=128, D=28 mm • d = 0.37 . • In linear scanning, the angle of inspection is 0, and the grating lobes are at such high angles that the interference from them is minimal (cf. obliquity factor reducing pressure field to 0 at .=. for soft baffles such as ultrasonic transducers [Cobbold 2007, pp. 447-448]) • Therefore, the requirement for d<0.5. is relaxed for linear arrays • However, in phased arrays angle of inspection varies substantially (e.g. -./4 to ./4), so grating lobes can have an effect, and d<0.5. adhered to • Grating lobes arise naturally out of CW analysis (see next lecture). Simple way to reduce grating lobes: shorten the pulse! Medical diagnostic systems – Beamforming in ultrasound Apodization – beam shaping [Szabo 2004, p. 193] • Again consider statement that pressure in focal plane is 2-D Fourier transform of apodization function (weightings used for the elements) • Therefore, to provide a sharp beam (good imaging resolution) on transmit and receive, transducer should be as large as possible • In the limiting case of an infinite plane or enveloping hemisphere, the beam would be an impulse • However, in the case of a finite aperture, the beam smears • This is in analogy with estimating the power spectrum of a signal from a limited time window • As in power spectrum estimation, uniform apodization causes sinc beam • Borrowing from spectral estimation techniques, different apodization (windowing) functions [Harris 1978] can be used to reduce the amplitude of sidelobes, at the expense of increasing the main lobe width Medical diagnostic systems – Beamforming in ultrasound Dynamic receive apodization • Increase receiving subaperture with depth • constant f# on receive • blurring (PSF) more uniform with depth • easier to deconvolve blurring computationally • easier to “deconvolve” (interpret image) by eye active subaperture Medical diagnostic systems – Beamforming in ultrasound Adaptive beamforming [Holm et al. 2009] • Technique taken from radar and sonar (as so often with ultrasound!) • Original aim was to cancel out jamming signal from enemy • Vary element weights w (apodization) so as to reject signal from elsewhere • This allows placing of sidelobes at regions of low energy (or low scattering) while maintaining position of main lobe (i.e. focus) • Example: Capon beamformer – minimise beamformed signal energy while keeping |w|=1 • See [Holm et al. 2009] for illustrative images References [Agricola et al. 2010] Real-time three dimensional transesophageal echocardiography: technical aspects and clinical applications [Anderson 2006] A seminar on k-space applied to medical ultrasound. http://dukemil.bme.duke.edu/Ultrasound/k-space/index.htm [B-K Medical 2003] Users guide for casablanca engine interface [Burns 2005] Introduction to the physical principles of ultrasound imaging and Doppler. http://medbio.utoronto.ca/students/courses/mbp1007/Fall2009/MBP1007_Burns_Utras ound.pdf [Cobbold 2007] Foundations of biomedical ultrasound [Gyφngy 2010] Passive cavitation mapping for monitoring ultrasound therapy [Holm et al. 2009] Capon beamforming for active ultrasound systems. http://heim.ifi.uio.no/~sverre/papers/09_MinVar-DSP-workshop.pdf [Jensen et al. 2007] System architecture of an experimental synthetic aperture real-time ultrasound system. http://www.jp-embedded.com/download/press/saurus/preprint.pdf [Szabo 2004] Diagnostic ultrasound imaging: Inside out [Van Veen and Buckley 1988] Beamforming: A versatile approach to spatial filtering