Ultrasound Imaging
Yao Wang Polytechnic University, Brooklyn, NY 11201 Based on J. L. Prince and J. M. Links, Medical Imaging Signals and Systems, and lecture notes by Prince. Figures are from the textbook except otherwise noted.
Ultrasound Imaging Yao Wang Polytechnic University, Brooklyn, NY - - PowerPoint PPT Presentation
Ultrasound Imaging Yao Wang Polytechnic University, Brooklyn, NY 11201 Based on J. L. Prince and J. M. Links, Medical Imaging Signals and Systems, and lecture notes by Prince. Figures are from the textbook except otherwise noted. Lecture
Yao Wang Polytechnic University, Brooklyn, NY 11201 Based on J. L. Prince and J. M. Links, Medical Imaging Signals and Systems, and lecture notes by Prince. Figures are from the textbook except otherwise noted.
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 2
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 3
imaging)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 4
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 5
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 6
(damping)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 7
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 8
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 9
Product of a decaying envelop and a sinusoidal function
( )
t f t ne 2 cos ) ( π
( )
[ ] [ ]
) ( ) ( 2 1 ) ( ) ( * ) ( 2 1 ) ( 2 1 2 cos ) (
2 2
f f N f f N f f f f f N e e t n t f t n
e e e t f j t f j e e
+ + − = + + − ⇔ + =
−
δ δ π
π π
Narrow band (around f0) pulse
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 10
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 11
Vibrating in z direction
Each point (x0,y0) produces a pressure wave p(x,y,z,t;x0,y0)
A scatter at (x,y,z) reflects p(x,y,z;t); generating p_s(x0’,y0’;t;x,y,z)
Wave
p(x,y,z;t) is superposition of above waves
transducer face The scattered signals
to a voltage signal r(x,y,z;t), all due to a single scatter at (x,y,z)
) , , (
' ' 0 y
x
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 12
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 13
(each point acts as a dipole rather than a monopole, Signal is strongest in the direction orthogonal to the dipole)
face
wave, with signal at transducer position x0’,y0’
the transducer
( ) ( )
2 2 2 2 1 2
) ( ) , ; , , , ( z y y x x r r c t n r z y x t z y x p + − + − = − =
−
0, face; in ) , ( if , 1 ) , ( ) ( ) , ( ) , , , (
1 2
= = − = ∫∫
−
y x y x s dy dx r c t n r z y x s t z y x p
) ; , , ( 1 ) , , ( ) , , (
' 1 ' ' '
r c t z y x p r z y x R t y x ps
−
− =
' ' ' 1 1 2 2 ' ' ' ' '
) ( ) , ( ' ) ' , ' ( ) , , ( ) ; , ( ) ' , ' ( ) , , , ( dy dx dy dx r c r c t n r z y x s r z y x s z y x KR dy dx t y x p r z y x s K t z y x r
s
− − = =
− −
= dxdydz y z y x r t r ) , , , ( ) (
z/r0’ due to dipole pattern
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 14
This approximation enables us to separate integration over x0,y0 and that over x0’,y0’
) ' ( ) ( ) ' 2 ( 2 ) ' ( 2 1 1 1 ) ' ( 2 1 1 1 1 2
1 1 1 1 1
exponent in the 2 using Also ) 2 ( ) ' ( ' : ion Approximat ) ' ( ) ' ( ) ( ) (
r z jk r z jk r r z c f j r c r c t f j e e r c r c t f j j e t f j j e
e e e e z/c t z c t n r c r c t n z r r e e r c r c t n r c r c t n e e t n t n
− − − − − − − − − − − − − − − − − − − − −
= = = − = − − ≈ ≈ − − = − − =
− − − − −
π π π φ π φ
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 15
Basic pulse-echo signal equation (depends only on transducer face, not pulse)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 16
This is the same result that we would have got had we assumed that all points on the transducer act as spherical wave generators and receivers
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 17
Fresnel region Fraunhofer region (far field, beam spreading)
λ 4
2
D
Geometric approximation: Wave is confined in a cylinder
Vibrating flat plate Near field boundary (NFB)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 18
S(u,v): FT of s(x,y) (valid in Fresnel region) (further approximation, in far field)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 19
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 20
ct
a
µ −
2
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 21
Since n(t) and q(x,y,z) depend only on the source, e_c(t) is affected by the reflectivity R(x,y,z) in the imaged body, or e_c(t) ~ R(0,0,z=ct/2)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 22
Previous result assumes the transducer is located at (0,0,0). When the transducer is at arbitrary (x0,y0), q(x,y,z) is changed to q(x-x0,y-y0,z) at (x0,y0,ct/2)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 23
R(x,y,z) (z=tc/2) convolved with h
defines the resolution cell of the imaging system
dimension as the transducer face in (x,y) plane, and the extend in z- direction = cT/2, if T is the length of the transmit pulse.
) , ( ) , ( ~ 2 / ) , ( ~ ) , , ( ) , , ( * * * ) , , ( ) , , ( ˆ Then ) , ( ) , , ( ) , , ( ~ assume we ion, approximat geometric With
2 2
y x s y x s c z n y x s z y x h z y x h e z y x R K z y x R y x s z y x zq z y x q
e kz j
− − = = = ≈ =
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 24
plate square for function sinc ) , (
ansform Fourier tr : ) , ( 2 / , ) , , ( ) , , ( * * * ) , , ( ) , , ( ˆ
2 2
= = = y x s v u S c z n z y z x S z y x h z y x h e z y x R K z y x R
e kz j
λ λ
Note that the spatial extend of the resolution cell now increases with Z. That means the lateral resolution decreases with Z. The lateral extension is also not finite (sinc function)
λ λ λ λ λ / for ) , ( ) , ( ) , , ( ) , , ( ~ axis) near i.e., y, x, z hen (correct w 2 ) ( assuming and ion approximat Fraunhofer With
2 2 ) ( 2 2
2 2
D z z y z x S e z y z x S z y x zq z y x q z y x
z y x jk
≥ ≈ = = >> ≈ +
+
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 25
(\lamda) includes 2 cycles of the modulating wave. An ultra sound pulse typically includes several \lambda
2
kz j
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 26
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 27
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 28
– Measure the reflectivity at different depth below the transducer position The horizontal axis can be interpreted as z, with z=t c/2
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 29
– High freq is used to produce very high axial resolution – Attenuation due to high freq is not a problem as the desired imaging depth is small
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 30
separate column of a 2D image, for a fixed transducer position
– Motion of an object point along the transducer axis (z) is revealed by a bright trace moving up and down across the image – Often used to image motion of the heart valves, in conjunction with the ECG
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 31
new pulse after each movement
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 32
– No tomographic measurement and reconstruction is necessary!
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 33
direction perpendicular to the place of each B-mode scan
From [Webb2003]
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 34
– \alpha is linearly increasing with f: alpha= a f
(often 80dB)
– z is the round trip disiance, depth of penetration is half – Lower penetration at higher frequency
– L=80dB, a=1, d_p=40/f (Mhz) cm
z
a
µ −
10 10
z
a
a
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 35
2d_p/c
– N pulses are required in each frame – Total time =N T_R – Frame rate F =1/NT_R=acf/NL
– 10-100 frames/sec
– reducing field of view and consequently N – Increasing frequency (at the cost of reduced depth of penetration)
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 36
multiple transducers
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 37
– Much smaller transducer elements than in linear array – Use electronic steering/focusing to vary transmit and receive beam directions
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 38
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 39
travels than T1:
to arrive at a point at the same time as T0, T1 should be delayed by
fires at
θ sin d d = ∆
c d c d T / sin / θ = ∆ =
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 40
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 41
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 42
Ti should receive T sec later than T_i+1
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 43
T0 fires at direction \theta, and all Ti’s receive after a certain delay, so that they are all receiving signal from the same point at a particular time
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 44
\sin\theta, z(t)=ct \cos \theta (possible scatterer)
c z x id z x ti / ) (
2 2 2 2
+ − + + =
c z x id z x t t t
i i
/ ) (
2 2 2 2
+ − − + = − = ∆
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 45
– Velocity vs. frequency shift (Doppler Freq)
– Stenosis or narrowing of the arteries causes blood flow velocity change – Children at risk of stroke have cerebral blood velocities 3-4 times
– Continuous wave (CW) Doppler measurement
– Pulse mode Doppler measurement
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 46
motion of the source and receiver
– Source moving away, wavelength longer, lower freq – Source moving closer, wavelength shorter, higher freq.
receiver, towards moving source : 2 / receiver, from away moving source : 2 / cos cos cos : frequency Doppler cos tor) motion vec source and ector receiver v source- between (angle : direction a in moving source : case General is frequency temporal Equivalent h wavelengt equivalent the as
thought be can above The .
distance a moves crest motion, source With /
distance a moves source motion source without /
distance a moves in wave crest 1 period
: direction wave the
speed with moving ) (freq source > <= < > ≈ − = − = − = > + = = + = + = = = =
D D O O O T D O T
T T
f f c v f v c v f f f f v c c f f v c c c f vT cT vT cT f v vT f c cT /f T v f π θ π θ θ θ θ θ θ λ λ
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 47
– Transducer transmitting a wave at freq fs, wavelength =c/fs – Object is a moving receiver with speed v, with angle θ – Target moving away (θ>= π/2), sound moves slower – Target moving closer (θ<=π/2), sound moves faster
receiver, towards moving source : 2 / receiver, from away moving source : 2 / cos : frequency Doppler cos > <= < >= = − = + =
D D S S O D S O
f f f c v f f f f c v c f π θ π θ θ θ
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 48
– Transmit wave at freq fs to object – Object moves with velocity v at angle θ – Object receives a wave with freq fo =(c+vcosq)/c fs – The object (scatterer) reflects this wave (acting as a moving source) – Transducer receives this wave with freq – Doppler freq
speed from doppler freq.
–
S O T
S S S T D
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 49
– No depth information – Assumes a single moving object below the transducer
– Use AM modulation to deduce the frequency shift f_D – Use time-domain correlation
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 50
Frequency Counter Spectrum Analyzer
From Graber: lecture notes for BMI F05
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 51
From [Webb2003]
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 52
From: Graber Lecture notes for BMS F05
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 53
signals detected at two different times
(correspondingly distance traveled) that yields maximum correlation
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 54
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 55
– Transmits short pulses and receives backscattered signals a number of times
– See Figure 3.24 in [Webb]
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 56
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 57
– Interlaced B-scan and color encoded Doppler images ⇒ limits acquisition rate to 2 kHz (freezing of B-scan image possible) – Variation of depth window (delay) allows 2D mapping (4-18 pulses per volume)
From: Graber Lecture notes for BMS F05
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 58
From: Graber Lecture notes for BMS F05
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 59
– Obstetrics and gynecology
– Breast imaging – Musculoskeletal structure – Cardiac diseases
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 60
– Properties of piezoelectric crystal, needs for damping and matching layer
– Relation between excitation pulse n(t), transducer face s(x,y) and received signal r(t) – General equation and different approximation – Need for time gain compensation – Envelope of received signal after time gain compensation is an estimate
– Beam width increases with z
– A-mode (reflectivity in z for fixed (x,y) position), M-mode (motion trace in z for fixed (x,y)), B-mode (reflectivity in one cross section) – 3D imaging – Doppler imaging
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 61
– Time delay for transmit element for steering and focusing – Time delay for receive focusing and dynamic focusing
– CW vs. time gated
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 62
EL5823 Ultrasound Imaging Yao Wang, Polytechnic U., Brooklyn 63
– Prince and Links, Medical Imaging Signals and Systems, Sec.
– Handouts: Webb Sec 3.10-3.13
– P11.4 (note T should be T=2 \lambda/c) – P11.6 – P11.7 – P11.8 – P11.9 – P11.14 – P10.12 (note f_R should be (c+v) f_0/c ) – P10.13