뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow...

56
超音波影像實驗室 Ultrasonic Imaging Laboratory Dr. Pai-Chi Li (李百祺博士) Room 303 Department of Electrical Engineering National Taiwan University No.1, Sec.4, Roosevelt Road Taipei, Taiwan 106, R.O.C. TEL: + 886-2-23635251-303 Http:// land.ee.ntu.edu.tw

Transcript of 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow...

Page 1: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

超音波影像實驗室Ultrasonic Imaging Laboratory

Dr. Pai-Chi Li (李百祺博士)

Room 303Department of Electrical EngineeringNational Taiwan UniversityNo.1, Sec.4, Roosevelt RoadTaipei, Taiwan 106, R.O.C.TEL: + 886-2-23635251-303Http:// land.ee.ntu.edu.tw

Page 2: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Ultrasonic Imaging Laboratory Members

Advisor: Dr. Pai-Chi LiPh. D. Student: 4 Master Degree Student: 5Administrators: 1Research Assistant: 3Graduated Student: 7 (M.S.)

Page 3: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

PERSONAL EDUCATION

National Taiwan University 1983-1987 B.S. Electrical EngineeringUniversity of Michigan 1989-1990 M.S. EE: SystemsUniversity of Michigan 1991-1994 Ph.D. EE: Systems

RESEARCH INTEREST

Signal and imaging processing, Ultrasonic medical imaging

Pai-Chi Li (Associate Professor, Senior Member of IEEE )Electrical EngineeringNational Taiwan University+ [email protected]

WORK EXPERIENCE

Adjunct Associate Investigator Health Research Institutes, Taiwan 2001-present.Member of Technical Staff Acuson Corporation, U.S.A. 1994-1997.

Page 4: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

PROGRAM OVERVIEW

1. Ultrasonic Synthetic Aperture Imagingand Adaptive Imaging

2. Ultrasonic Nonlinear Imaging 3. Ultrasonic Elastic Imaging4. 3-D Ultrasonic Imaging5. High Frequency Ultrasonic Imaging6. Blood Flow Estimation Using Ultrasonic Contrast Agent

Page 5: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

1. Ultrasonic Synthetic Aperture Imagingand Adaptive Imaging

Page 6: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

• Ultrasonic Synthetic Aperture Imaging– Filter based synthetic focusing technique

Ran

ge

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

Baseband Demodulation

sinθ

Beam

Buffer

Image

Buffer

Filter Bank

Ran

ge

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

Baseband Demodulation

sinθ

Beam

Buffer

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

�������������������������������������������������

Baseband Demodulation

sinθ

Beam

Buffer

Image

Buffer

Filter BankFilter Bank

T r a n s d u c e r A /DB a s e b a n d

D e m o d u la t io n

B e a m B u f f e r

R a n g e - D e p e n d e n t

F i l te r B a n kI m a g e B u f f e r

S ig n a l P r o c e s s in g

S c a n C o n v e r s io n

D is p la y

B e a m f o r m e rT r a n s d u c e r A /DB a s e b a n d

D e m o d u la t io n

B e a m B u f f e r

R a n g e - D e p e n d e n t

F i l te r B a n kI m a g e B u f f e r

S ig n a l P r o c e s s in g

S c a n C o n v e r s io n

D is p la y

B e a m f o r m e r

M.-L. Li and P.-C. Li, "Filter Based Synthetic Transmit and Receive Focusing", Ultrasonic Imaging, Vol. 23, pp. 73-89, April, 2001.

Page 7: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Adaptive Imaging

• A new adaptive imaging technique using generalized coherence factor (GCF) is proposed.

• GCF is derived based on the spectrum of the received array data along the array direction.

• GCF is an index of beamforming quality • GCF is used as a weighting factor to the

reconstructed image.

Page 8: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

IdeaSpectrum (fx)

No aberration

10 dB/div

Range (R

)

Channel (x)

Aberrated

10 dB/div

More high frequency components

Page 9: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Generalized Coherence Factor

energy spectral totalrangefrequency low specified-prea hinenergy wit spectralGCF =

10 dB/div• High GCF corresponds to good focusing

quality and the image intensity should be maintained

No aberration

10 dB/div Lower GCF should be used to reduce the image data because significant beamforming errors are present

Aberrated

Page 10: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

ResultsSimulation Experiment

Anechoic cyst

Aberrated

GCF Corrected

Page 11: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

2. Ultrasonic Nonlinear Imaging

Page 12: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Ultrasonic Tissue Harmonic Imaging• Tissue Harmonic

– the harmonic component generated from finite amplitude distortion

t

PressureBefore distortion

After distortion

MHz

Received Signal

MHz

Transmit Signal

MHz

SUM

Positive pulse

Negative pulse

• Pulse Inversion – Better fundamental rejection, lower frame rate

Page 13: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Ultrasonic Tissue Harmonic Imaging

Page 14: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Ultrasonic Tissue Harmonic Imaging• Harmonic spatial covariance analysis

• Effects of SNR • Effects of sound velocity inhomogeneities

Single Crystal Transmitter

Array Receiver

Image Plane

Speckle Target

Page 15: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Ultrasonic Tissue Harmonic Imaging• Motion artifacts of Pulse Inversion Technique

• Effects of SNR • Effects of sound velocity inhomogeneities

0 0.1 0.2 0.3 0.4 0.5 0.600.20.40.60.8

11.2

Axial Displacement (λ)

Inte

nsity

FundamentalTissue Harmonic

-6 -4 -2 0 2 4 6-40

-30

-20

-10

0

Lateral Position(mm)dB

No Motion 0.05 mm0.1 mm

Page 16: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

3. Ultrasonic Elastic Imaging

Page 17: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Ultrasonic elastic imaging

• Ultrasonic strain compounding image based on a fast speckle tracking algorithm

• Strain compounding technique– Improve contrast resolution of the image– Steps:

Obtain an uncompressed image as a basisApplying an external force on the object yields deformationModify the deformation in the image planeAverage the modified and original images

Page 18: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

2D fast speckle tracking algorithm• Block Sum Pyramid

– Take threshold: SADmin– Reduce the computations of SAD

• Multilevel Block matching– Reduce numbers of points to be searched– We use 2 levels– 1st level: window size= w*w, 9 points– 2nd level: window size= ½*w* ½ *w, all points

)2,2()12,2(

)2,12()12,12(),(1

jiXjiX

jiXjiXjiX

mm

mmm

+−+

−+−−=−

∑ ∑= =

−=m m

i j

mmm jiYjiXYXSAD2

1

2

1

),(),(),(

block: X

2m

2m

2mblock: Y

2m

Compare SADm

Pyramid structure

w

w1st level

Page 19: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Results• Algorithm performance For 121 pixels

12 : 1.5 : 10.998 (s)1.53 (s)12.3 (s)C3.6 : 2 : 17.64 (s)14.17 (s)27.14 (s)Matlab

RatioBSPA &Multilevel

BSPAFSALanguage BSP & Multilevelalgorithm is indeed not only faster than traditional algorithm, but also as accurateas traditional one.

• Compounding imageHeight (mm)

Width (mm)5

5

10

15

Width (mm)

Height (mm)

5

5

10

15

L: original image R: compounding image with BSP & Multilevel

Liver

SNR=0.1164

SNR

Computing speed

Page 20: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Height (mm)

Width (mm)5 10 15 20

5

10

15

20

23.5

Height (mm)

Width (mm)5 10 15 20

5

10

15

20

23.5

Cont’dThyroid

SNR=0.9625

SNR

Computing speed

L: original image R: compounded image with BSP & Multilevel

Height (mm)

Width (mm)5 10 15 20 25

5

10

15

20

25

30

Height (mm)

Width (mm)5 10 15 20 25

5

10

15

20

25

30

Breast

SNR=0.1692

SNR

Computing speed

Page 21: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Young’s Modulus Measurements of Human Liver and Correlation with Pathological Findings

T i m e ( s )

0 2 4 6 8 1 0 1 2 1 4 1 6 1 8

Re

ad

ing

fro

m e

lec

tro

nic

ba

1

2

3

4

5

6

7

S t r a i n

0 . 0 9 0 . 1 0 0 . 1 1 0 . 1 2 0 . 1 3 0 . 1 4 0 . 1 5 0 . 1 6

Str

es

s (

Ne

wto

n/m

2 )

1 5 0

2 0 0

2 5 0

3 0 0

3 5 0

4 0 0

4 5 0

5 0 0

5 5 0

S l o p e ( Y o u n g 's m o d u l u s )= 1 7 0 1 . 3 P a s c a l s

The Experimental Set-up The Readings from Electrical Balance as a Function of Time

The Stress- Strain Curve

Page 22: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

The Young’s modulus of normal liver, cirrhotic liver and hepatic tumors

Page 23: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Tissue characterization of Ultrasonic B-image

• Compare normal with cirrhotic liver– Statistical method– Conventional and non-separable wavelet

decomposition method

Normal liver Cirrhotic liver

Page 24: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

4. 3-D Ultrasound Imaging

Page 25: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Gyro Angle Detector for 3D UltrasoundGyro Angle Detector for 3D Ultrasound

• Hardware:– PG-03 Piezo GYRO(GWS)– 8254 for counting pulse width– 8051 with RS232 interface– IBM compatible PC

• Software:– Visual BASIC

Page 26: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

A Free-Hand 3D Ultrasound Imaging System

Page 27: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Correlation-Based Analysis for Complex Motion

Acquire 2 images

Find x’、z’ and α’

Images are spatial matched

Calculate 2D C.C.

Find the best matched combination

x’、y’、z’、α’、β’、γ’

Speckle tracking

Correct the x’ and z’ motions

Go through all β,γ, Y combinations

Page 28: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

A 3D System Integration

3D Rendering

• Platform:- Win NT、OpenGL

Frame Grabber

Baby Phantom ATL UM-9

Page 29: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

5. High Frequency Ultrasonic Imaging

Page 30: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High frequency ultrasonic imaging system

AxialfilterDemod.

10 20 30 40 50 60 70 80 90 100-60

-50

-40

-30

-20

-10

0

10

frequency(MHz)

Inte

nsity

(dB

)

PR5900 0.4xGCPuls0.6xGCPuls•Fully digital system architecture

•50 MHz center frequency•60 % fractional bandwidth

ADC

Lateralfilter

Codedexcitation

DAC

Hardware DesignTissue attenuationDynamic focus

Page 31: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Wire phantomDiameter = 52µm

• 52um nylon wire phantom• Gaussian Pulse Gaussian Chirp Pulse

60

0.4 0.8 1.2 1.6 2

10.01

5010.78

4011.55

dept

h(m

m)

30

12.32

20

13.0910

13.860.4 0.8 1.2 1.6 2

Position (mm)

Page 32: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Resolution test

•Lateral projection: axial projection•Spatial resolution is about 60 µm

0

0 0.2 0.4 0.6 0.8 1-60

-50

-40

-30

-20

-10

0GPS GCPS AFT PC

Position (mm)

Inte

nsity

(dB

)

-10

-20

-30

-40

-50

-60

GPS GCPS AFT PC

11.4 11.7 12.3 12.612Depth (mm)

Page 33: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

in-vitro pig eye image

0 mm

sclera

irislens

cornea

anterior chamberDep

th

6 mm16 mm0 mm

Position

Page 34: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Tissue harmonic imaging

• Pulse inversion technique cancels fundamental signal500

520

540

560

580

600

620

640

660

680

700

50

45

40

35

30

25

20

15

10

5

20 40 60 80 10020 40 60 80 100

Page 35: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Ultrasound Doppler

• 50MHz High Frequency Ultrasound : wideband transmitted signal (short transmitted pulse) and narrow lateral beamwidth→ better spatial and velocity resolution (down to mm/s), capable of estimating low velocities blood flow in small vessels.

D A C

A D C

P C

5 0 0 M H z

2 0 0 M H z

fc≒ 5 0 M H z

flow

3 D P o sitio n ing

T-P o rt P o w er A m p .

L N A

Page 36: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Ultrasound Doppler

In-Vitro Flow Estimation : Autocorrelation Technique

R t St S d( ) ( ) ()≡ + ∗

−∞

∫ τ τ τ

TTfπ

θ2

)(= , T : PRIR t R tej t( ) ( )( )= θ

Page 37: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Ultrasound Doppler

In-Vitro Flow Estimation : WMLE Technique

rk’(t) Delay Line of Length

(P-k)*T*(1+2v/c)

r1’(t)

r2’(t)

Σ h(t) 2 l(v)

Bank of delay lines, and filter h(t) matched to the expected demodulated echo signals which correspond to various velocities. The maximum likelihood velocity is then given by the filter with the largest output.

Page 38: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Ultrasound Doppler

In-Vitro 2D Flow Data : 500µm diameter cyst, with maximum velocity≒20mm/s

Page 39: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Ultrasound Experiment25~50MHz

Basic System Diagram Flow System

200MHz DA

500MHz AD

200MHz DA

500MHz AD500MHz AD

Experiment Condition

200MHz DA excite 25~50 MHz coded ultrasound wave500MHz AD receivingWide band ultrasound transducerUp to 20KHz PRF

Page 40: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Flow EstimationRF Butterfly Search (Multi-line)

Experimental -Flow Resultλ

Traditional Butterfly Search Line

10 20 30 40 50 60 70-10

-5

0

5

10

15

20

mm/sec

Multiple Butterfly Search Lines

Traditional RF Butterfly

Multi-line Butterfly

Page 41: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Flow EstimationColor Flow Image

-60

-40

-20

0

20

40

60

10 20 30 40 50 60 70

5

10

15

20

25

30

35

40

45

50

mm/sec

Page 42: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

High Frequency Harmonic Image25MHz 300 m Cyst imageµ

0 20 40 60 80 100-40

-30

-20

-10

0

10

dB

300µ m

Fundamental Image

5

10

15

20

25

30

35

40

45

50450

Fundamental Image500

Harmonic Image550

600

650

70020 40 60 80 100

Harmonic Image

5

10

15

20

25

30

35

40

45

50450

500

550

600

650

Frequency (MHz)700

20 40 60 80 100

Page 43: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

6. Blood Flow Estimation Using Ultrasonic Contrast Agent

Page 44: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Doppler Blood Flow Estimation (I)

- Doppler Angle Estimation- Blood Flow Estimation

Pulse/receiver

A/DControl and

Data Processing

Damper

Blood Mimicking Fluid

Tubing Pump

Flow Meter

trigger

Data Acquisition System

VXI mainframe

5MHz

PRI: 200 secµ

water tank60 65 70 75

0.20.4

0.60.81

Depth(mm)

Pow

er

0 20 40 60

0.20.4

0.60.81

Slow Time Index

Pow

er

0 2-4

-2

0

2

4

Mag

nitu

de

0 20 40 60-4

-2

0

2

Mag

nitu

de

0 08

Condition 1 Condition 2

Original Doppler signals AR Extrapolation Signals4

0 40 60Slow Time Index Slow Time Index

1. P.-C. Li, C.-J. Cheng and C.-C. Shen, “Doppler Angle Estimation Using Correlation”, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 47, No. 1, pp. 188-196, 2000.

2. P.-C. Li, C.-J. Cheng and C.-K Yeh, "On the Velocity Estimation Using Speckle Decorrelation", IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 48, No. 4, pp. 1084-1091, July, 2001.

3. C.-K. Yeh and P.-C. Li, "Doppler Angle Estimation Using AR Modeling", IEEE Transactions onUltrasonics, Ferroelectrics and Frequency Control. June, 2002

Page 45: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Blood Flow Estimation Using Ultrasonic Contrast Agent

- Blood Flow Estimation (Indicator-Dilution Theory)- Time-Vary Method

Reservoir

ReservoirPump

Syringe

Flow Meter

Damper

I

Mixing Chamber

7-MHz Linear Array Transducer

50µm

Time

Inte

nsity

1 sec

Time-Intensity Curve (TIC)

300 sec

1. C.-K. Yeh S.-W. Wang and P.-C. Li, Feasibility Study on the Time-Intensity Based Blood Flow Measurements Using Deconvolution,” Ultrasonic Imaging, vol. 23, pp. 90-105, April, 2001,

2. P.-C. Li, C.-K. Yeh and S.-W. Wang, "Time-Intensity Based Volumetric Flow Measurements: An In Vitro Study", Ultrasound in Medicine and Biology. vol. 28, no. 3, pp. 349–358, 2002

Page 46: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Doppler Blood Flow Estimation in Pulsatile Flow (II)

- Doppler Angle Estimation- Blood Flow Estimation

Dep

th

Time

0 0.2 0.4 0.635

40

45

50

55 Ave Flow: 11 ml/sec

Time (sec)

Dop

pler

Ang

le

Data Acquisition System

CompuFlow 1000(Blood Mimicking

Fluid)trigger

DAC 200ECG output

Pulse/receiver

Gage A/D Gage D/A

PRI: 400

triggerRF signals

PC

secµ

water tank

3.5 MHz

trigger

1. Chih-Kuang Yeh and P. C. Li, “Doppler Angle Estimation of Pulsatile Flow Using AR Modeling,” Ultrasonic Imaging, 2002 (submitted )

Page 47: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Blood Flow Estimation Using Ultrasonic Contrast Agent

- Shadowing Effect- Input and Output Time-Intensity Curves (IOTIC)

Flow Meter

Reservoir

AA

BB

CC

Reservoir

Pump

Syringe

260 ml

Damper

I

Water Tank

Mixing Chamber

Transducer (curve array)Dialysis Cartridge

1. Chih-Kuang Yeh and P. C. Li, “Contrast specific ultrasonic flow measurements based on both input and output time intensities,” Ultrasound in Medical & Biology, 2002 (Submitted)

Page 48: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Assessment of Parameters in Pulsatile Flow using Ultrasound Contrast Agent

• Provide a model for the assessment of perfusion characteristics

• Dilution theory– MTT : mean transient time– theory : V/Q (ideal)

• LTI system TV (time-varying)

∫∫

∞×

=

0

0

)(

)(

dttn

dttntMTT

o

o

no(t)

baseline = 0t = 0

Page 49: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Experimental setup

UHDC

IISyringe

B

Transducer

A

Page 50: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Simulation methods

Superposition theory :

0

),( /)( <

= −

ξξ ξ

tth VtQ

)(),()( ξξξ dIthtI Io ∫∞

∞−=

≥ξte

),( ξth)(ξQ

Page 51: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

Results: the theoretical values & MTT

simulation result experimental result

Page 52: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

EQUIPMENT (1) Panametrics Model 5900 PR (pulser / receiver) 1 200 MHz digital

GW Model GFG-813 (function generator) 1

GW Model GFG-8016D (function generator) 1

HP Model 54603B (oscilloscope) 1 60 MHz 附 probe X 4 power line X 1

Tektronix Model TDS 380 (oscilloscope)1 Two channel / digital

real-time / 400MHz / 2GS/s

TAIK Model TK-12001D (DC power supply) 1

EPE Model EP-3000 (DC power supply) 1

Cimarec Model SP46925 (stirrer/heater) 1 S/N:1069980822742

OHAUS Model IP12KS 1

Microtime Model 51/52-E (WINICE) 1 S/N:A00I955172

Cole Parmer Model 07596-20 (damper) 1

Cole Parmer Model 77021-60 (pumper) 1

GaGe Model CompuGen 1100 (AFG)1 ISA interface 1M

RAM on board S/N: G00086 80 MHz

GaGe Model CompuScpoe 12100 (A/D)1 PCI interface 1M

RAM on board S/N: P10243 100 MHz

Amplifier Research LN1000A(LNA) 1 DC transformer

Amplifier Research P25A250A(PA) 1

Signatec DAC2001 User manual ,

CDX1,BNC to SMD cableX3

Signatec PDA500 1 User manual , CDX1

Signatec PMP8-A 1 User manual , CDX1

Panametrics HF cable 3 1 ft., 3 ft., 6ft.

Page 53: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

EQUIPMENT (2): Transducer

Page 54: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

EQUIPMENT (3): Commercial Ultrasound Machine

(SonoSite)(hand-carried ultrasound system) (ATL Ultramark 9 HDI)

Page 55: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

EQUIPMENT (4): Phantom

Breast (II) BabyBreast (I)

Page 56: 뙗궵ꩩ뱶릳맪엧ꯇ - 國立臺灣大學home.ee.ntu.edu.tw/classnotes/LAB.pdfDoppler Blood Flow Estimation (I) - Doppler Angle Estimation - Blood Flow Estimation Pulse/receiver

EQUIPMENT (5)

(Digital Sonifier, BRANSON)Making Microbubbles

(UHDC Flow System)Simulation Physical Pulsatile Flow