An optimization method for correction of ultrasound

0 downloads 0 Views 2MB Size Report
3D EMA (NDI WAVE) sensors on a ruler (to simulate speaker's head) and probe. - Manipulate ruler/probe movements. Pick initial guess of origin. - Optimize for ...
Oct 4~6, 2017

An optimization method for correction of ultrasound probe-related contours to head-centric coordinates Wei-rong Chen1, Mark Tiede1, Shuwen Chen2 1 Haskins 2

Laboratories

The Chinese University of Hong Kong

[email protected], [email protected], [email protected] Ver: oct-06-2017

1

Challenge • One well-known challenge for using ultrasound in phonetic studies is the correction of probe/head movements. Tongue held still, probe moving

Relative movements between probe and head induce errors in ultrasound images.

2

Two solutions 1 Constraint method Head stabilizer headset (helmet) Example: Restrict probe/head/jaw movements. How: Challenges: 1. Difficult to guarantee purely no movements 2. Potentially affects naturalness of speech.

2 Correction method Example: How:

HOCUS system (Whalen et al, 2005) Track and transfer probe/head movements to ultrasound contours as correction

Challenges: 1. Difficult to define the origin of coordinates. 2. Accuracy is difficult to assess. Whalen, D. H., et al. (2005). "The Haskins Optically Corrected Ultrasound System (HOCUS).” JSLHR 48(3): 543-553.

3

Correction method • 3 tasks for correction method: - Mutual scaling (easy task) - Transformation of the orientations of the planes in coordinates (solved in HOCUS) Not settled yet! - Determining the mutual origin.

4

Objectives of the current work • This works demonstrates a simple way for correction methods to : 1. Calibrate the origin of coordinate system for ultrasound probe. 2. Assess the accuracy of the correction of probe/head movements. 5

“Floating ruler” experiment1 - 3D EMA (NDI WAVE) sensors on a ruler (to simulate speaker’s head) and probe - Manipulate ruler/probe movements. Pick initial guess of origin. - Optimize for the origin by aligning contours in ultrasound image (estimation) and in EMA (truth). - Use the optimal origin for the correction of ruler/probe movements.

1. Thanks

to Doug Whalen for this cool name!

6

Floating ruler experiment U1-U3: EMA sensors on probe RR1-RR3: EMA sensors on the upper part of ruler R1, R2: EMA sensors on the bottom part of ruler Ruler bottom edge: position of ruler bottom edge relative to R1R2, traced from photo image. Probe top surface: position of probe top surface relative to U1U3, traced from photo image.

7

Floating ruler experiment 6 settings: 1. Baseline: aligning ruler and probe in straight, centered position. 2. Rot_ccw: Rotate ruler by 20~25° counterclockwise. 3. Rot_cw: Rotate ruler by 20~25° clockwise. 4. Shift_L: Shift ruler to the left by approx. 2cm 5. Shift_R: Shift ruler to the right by approx. 2cm 6. Shift_D: Shift ruler down by approx. 1cm

8

Contour tracking in ultrasound images baseline

GOAL of head/probe movement correction: Accurately recover/transform the contours in 2-6 to fit the baseline, by removing probe/ruler movements 9

TASK: Co-register the origins in both coordinate systems Where is the origin [0,0] ?

10

STEP 1: Initial guess for the coordinate origin (based on our understanding of how ultrasound works)

11

Pick our best guess for the coordinate origins Ultrasound

EMA

12

STEP 2: Transfer probe movements to ruler (head) movements (Remove probe movements w.r.t. baseline)

13

After probe-movements are removed If the origin is correct: - ruler bottom contours in ultrasound and in EMA should perfectly align. Ultrasound:

EMA: 14

Calculate the error of alignments with initial guess

US = ultrasound

15

STEP 3: Find the optimal origin: 1. Define cost function J(θ) = mean distance between ruler bottom contours in ultrasound and in EMA, where θ = the position of origin. 2. Solve 𝜃෠ = 𝑎𝑟𝑔𝑚𝑎𝑥𝜃 𝐽(𝜃) using gradient descent.

16

Compare optimization results Initial guess:

Optimized: 17

STEP 4: Correct ruler/head movement with optimized origin 1. Calculate and store the ruler/head movements w.r.t baseline 2. Transfer ruler/head movements to ultrasound contours in 2-6 as correction.

18

Calculate and store the ruler movements w.r.t baseline For each setting in 2-6: Calculate the rigid body transformation (rotation R and translation T) from the baseline (setting1). B = baseline ruler position (setting 1), A(i) = ruler position of the setting i. B = R(i) x A(i) + T(i) R(i) = rotation from A(i) to B T(i) = translation from A(i) to B

19

FINAL RESULTS of correction: Update the origin in ultrasound by the optimal value Apply R and T to ultrasound contours to correct for probe/ruler movements

‘x’: initial guess of origin (0,0) ‘x’: optimal origin (-2.8, -2) Mean error after correction = 0.45 mm

20

RESULT: APPLY optimal origin to tongue contour correction:

21

Discussion • How accurate was our initial guess of origin? Actually not too bad. Initial guess [0,0] vs. optimal [-2.8, -2] in millimeter is not a huge difference. In most settings, the improvement of accuracy by the optimal origin is less than 1 mm.

22

Discussion • Why optimal origin and initial guess are different? • What is the source of this discrepancy? Some possibilities: Machine-internal process. Errors in hand-picking initial guess

23

Discussion • Since the initial guess is not bad, then why bother to find the optimal origin? Uncertain errors occur when hand-picking the initial guess. Uncertainties about the internal processes in ultrasound machine. Good initial guess for this machine doesn’t guarantee the same accuracy for another machine. The proposed “floating ruler” is a simple and cheap method of calibration. Just collect 6 frames only once for each probe in each machine. 24

Discussion • Alternatives for EMA? The use of EMA is for tracking probe/ruler(head) movements. Alternatives: - Optical tracking as in HOCUS system - Video-based landmark tracking (e.g., ‘blue-dot-tracking’) - 3D video camera (e.g., Microsoft Kinect) - 3D facial recognition technology (e.g., OpenFace)

25

Conclusion • The proposed “Floating ruler” method is a simple and necessary procedure for both calibration and accuracy assessment if head-correction for ultrasound contours is performed.

26

Future directions • Simplify and standardize the procedure. • As the transformational correction is only affected by ‘ rotation’, more settings with different probe angles need to be tested. • Apply the optimal origin to the head-correction for speech data.

27

Acknowledgements • We thank Doug Whalen, Christine Shadle, Emily Phillips, Kevin Roon for comments and helps. • This works is supported by NIH grant DC-002717 to Haskins Laboratories.

28