UBC Faculty Research and Publications

Single wall closed-form differential ultrasound calibration. Najafi, Mohammad 2012

You don't seem to have a PDF reader installed, try download the pdf

Item Metadata


Rohling_SPIE_8316_83162A.pdf [ 328.71kB ]
JSON: 1.0107564.json
JSON-LD: 1.0107564+ld.json
RDF/XML (Pretty): 1.0107564.xml
RDF/JSON: 1.0107564+rdf.json
Turtle: 1.0107564+rdf-turtle.txt
N-Triples: 1.0107564+rdf-ntriples.txt

Full Text

 Single Wall Closed-form Differential Ultrasound Calibration  Mohammad Najafi, Narges Afsham, Purang Abolmaesumi, Robert Rohling Department of Electrical and Computer Engineering, University of British Columbia, Vancouver, BC, Canada ABSTRACT In freehand 3D ultrasound, images are acquired while the position of the transducer is recorded with a tracking device. Calibration is essential in this technique to find the transformation from the image coordinates to the reference coordinate system. The single wall technique is a common calibration method because a simple plane phantom is used. Despite its advantages, such as ease of phantom construction and image analysis, this method requires large number of images to converge to the solution. One reason is a lack of a closed-form solution. Also, the technique uses slightly ill- conditioned sets of equations with a high condition number due to limited range of scanning motions that produce clear images of the plane. Here, a novel closed-form formulation has been proposed for the single wall calibration technique. Also, differential measurements of the plane image are used instead of absolute plane detection to improve accuracy. The closed-form solution leads to more accurate and robust results while providing an insight into understanding error propagation and finding the optimal set of transducer poses. Results have been compared to the conventional single wall technique. A residual error of 0.14 mm is achieved for the proposed method compared to 0.91 mm in the conventional approach. Keywords: Free-hand ultrasound, calibration, single wall phantom, closed-form, differential.  1. INTRODUCTION In freehand 3D ultrasound, which is a technique of acquiring ultrasound images with a tracked transducer, one important step is transducer calibration. It provides the transformation that relates any pixel in the ultrasound image to the corresponding point in a 3D world coordinate system. A variety of approaches for calibration of ultrasound have been investigated and they have been compared in terms of reconstruction accuracy, reproducibility, and acquisition time1. One typical method of calibration is to image an artificial object, known as a phantom, with known geometrical parameters, and combine the prior knowledge of the phantom with its ultrasound images to solve for the calibration parameters. Image features of the phantom can be points, lines or more complex shapes. Point-based methods are performed by repetitive scanning of a point target such as a spherical fiducial, intersection of two wires, or the tip of a stylus2. Since aligning the ultrasound beam at a single point is not easy, scanning along strings or wires can be used instead of single location. The three-wire phantom1, Hopkins phantom3, and Z or N shaped phantoms4, are examples of these methods. The final accuracy of point-based or wire-based methods are highly dependent on the accuracy of manual or semi-automatic selection of the point or line centroid in ultrasound images1. Given the presence of noise and artifacts, it is argued that a line is segmented easier and more accurately than a point1. The single wall phantom and its enhanced version, the Cambridge phantom5, take advantage of this idea and facilitate the automatic segmentation of the line. The Cambridge phantom is very accurate but requires custom apparatus, accurate mounting of the probe, and is difficult to use for a novice user1. The single wall phantom is simply a plate, made of e.g. plexiglass, and does not require construction of a phantom but is less accurate especially for low numbers of acquired images acquired approximately perpendicular to the plate (when the plate appears as a clear sharp line). Here it is worth mentioning that wedge phantoms6,7 also use planar surfaces but they require accurate and specific phantom design in comparison with the simple design of the wall phantom. One idea to improve the accuracy of the calibration is to use differential measurements instead of absolute position of the features in the image. In fact, in the single wall technique, the slope of the line is a differential measurement whose accuracy is less dependent on the sharpness of the line. The solution of the conventional single wall technique5 is based Medical Imaging 2012: Image-Guided Procedures, Robotic Interventions, and Modeling, edited by David R. Holmes III, Kenneth H. Wong, Proc. of SPIE Vol. 8316, 83162A © 2012 SPIE · CCC code: 1605-7422/12/$18 · doi: 10.1117/12.911459 Proc. of SPIE Vol. 8316  83162A-1 Downloaded from SPIE Digital Library on 05 Mar 2012 to Terms of Use:  http://spiedl.org/terms  on measuring the absolute position of points along the line in the image and solves for the calibration parameters by minimizing the residual error. The method is iterative due to lack of a closed-form solution. It also has slightly ill- conditioned sets of equations due to the limited range of scanning motions that produce clear images5. Here, a closed-form method is proposed that first solves for the rotation components of the calibration matrix using the differential measurements (line slopes) and then, solves for the translation components. This improves the accuracy, and at the same time provides an insight into finding the optimal set of poses for the transducer. 2. METHODS 2.1 Mathematical Framework The calibration goal is to find the six degree-of-freedom transformation from the image to the transducer coordinate system ( I TT ). The transformation from the transducer to the reference coordinate system ( T RT ) is known by the readings from the tracker (Figure 1).  Figure 1. Showing coordinates systems of the Tracker (reference), transducer (T) and the ultrasound image (I) in the calibration setup. The method is based on scanning a surface with the tracked ultrasound transducer in several different poses. The rotation part of the I TT  transformation ( xR ) is first calculated using slopes of the lines in the images. The transformation from the image to the reference, I RT , can be defined from two free vectors iU r  and iV r  and a point iP0  as follows. iU r  is a unit vector in the direction passing through the center of array elements (lateral) and iV r  is a unit vector in the direction of ultrasound beam (axial).  iP0  is the origin of imaging plane in the reference coordinate system, which is the same as the translation vector in I RT . Pixel to millimeter ratios in axial and lateral direction are defined as xS  and yS , respectively. Here, it has been assumed that xS  is known from the transducer specifications provided by the manufacturer, but yS  is unknown and is solved in the algorithm. Using the above notations for each pose of the transducer we have:  ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ = 0 0 1 I T T R i RRU r  , ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ = 0 1 0 I T T R i RRV r . (1) Proc. of SPIE Vol. 8316  83162A-2 Downloaded from SPIE Digital Library on 05 Mar 2012 to Terms of Use:  http://spiedl.org/terms  Here I T R (= xR ) is the rotation part of the I TT  transform. And T R R (= iR ) is the rotation part of the transducer to reference transform, T RT , which is known and comes from the tracker for each pose (i). Considering the value of iU r  and iV r  at the initial pose of the transducer ( 0U r  and 0V r ), we can define them at pose i using Eq. 1.  00 1 0 )( URURRU d iii rrr == −  , 0010 )( VRVRRV diii rrr == − . (2) d iR is the relative rotation from the initial pose to the ith pose, which is known since both 0R  and iR  are determined from the tracker. Also the origin of the image ( iP0 ) can be defined as:  iiT R I T T I I T T R i ttRttRTTP rrr +=+= ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = 1 0 0 0 0 , (3) where t r  is the unknown translation and it r  is the transducer to reference translation. Therefore each pixel of the image (x,y) can be described in the reference coordinates as:  iyixi VySUxSPP rr ++= 0 . (4) 2.2 Closed-form Solution Considering the pixels on the line appeared in the ultrasound image from the intersection of the wall and the ultrasound image plane, they should satisfy the plane equation (Figure 2):  0].)[( 110 =−++ nQVySUxSP iyixi r rr . (5)  Figure 2. A line will appear in the ultrasound image from the intersection of the wall. Any pixel on the line P(x,y) should satisfy the plane equation. Assume another point of image on the line (x2,y2); by subtracting Eq. 5 for the two points and then dividing by xS , and then dividing by yΔ  we have:.  0.. =Δ+Δ nVyknUx ii r rrr    ⇒   0.. =+ nVkmnU ii r rrr , (6) Proc. of SPIE Vol. 8316  83162A-3 Downloaded from SPIE Digital Library on 05 Mar 2012 to Terms of Use:  http://spiedl.org/terms  where x y S S k = , 12 xxx −=Δ , x ym Δ Δ=  and 12 yyy −=Δ . Now by substituting Eq. 2 into Eq. 6:  0.. 00 =+ nVRkmnUR diidi r rrr . (7) In the above equation 130 ×U r , 130 ×V r  and 13×n r  are unknown and im  and 33× d iR  are known. To solve this equation now we use this property of the Kronecker product:  )()()( XvecABAXBvec t ⊗= . (8) The above equation means that taking vec  from the product of three matrixes can be written as the Kronecker product of the last and the first multiplied by the vec of the middle matrix. Here by taking vec  from two sides of Eq. 7 we have:  0)()()()( 00 =⊗+⊗ dittiditt RvecnVkmRvecnU r rrv . (9) Now by putting the above linear equation for all poses (i=1 to N) into matrix form, we have:  [ ] 1 1189 1 9 1 18 111111 0 33113311 33113311 × × × = ⎥⎥ ⎥⎥ ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎢⎢ ⎢⎢ ⎣ ⎡ ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ N N d NN d NN d N d N dddd ky ky x x RmRmRR RmRmRR M M LL MOMOMM LL , (10) where ttt nUx )( 019 rvr ⊗=×  and ttt nVy )( 019 r rr ⊗=× . The goal is to find [ ]9191 kykyxxX LL= , the null space of the left side matrix and the answer should be scaled to satisfy:  1)9:1( 19 == ×xX r . (11) which is obvious since 0U r , 0V r  and nr  are unit vectors. To find xR we can use Eq. 1 for any pose (i) by substituting calculated values of 0U r  and 0V r from the above solution.  We can easily find k using the below equation:  kykX == ×19)18:10( r . (12) So far we have solved the rotation part ( xR ) of the calibration matrix ( I TT ). To solve the translation part ( t r ) we substitute Eq. 3 into Eq. 5:  dnQnVySnUxSnttR tiiyiixii ==+++ rr rrrrrr ...).( , (13) where d is a parameter of the plane equation which is also unknown. Proc. of SPIE Vol. 8316  83162A-4 Downloaded from SPIE Digital Library on 05 Mar 2012 to Terms of Use:  http://spiedl.org/terms  Eq. 13 is a linear equation in which all the parameters except t r and d are known. Now by considering Eq. 14 as:  nVySnUxSc iiyiixi rrrr .. += ,    iti Rnn rr = ,   ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ − − = ]1,[ ]1,[ 1 t N t n n G M . (14) Then we can rewrite Eq. 13 as:  ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ −=⎥⎦ ⎤⎢⎣ ⎡ Nc c d t G M r 1 . (15) Eq. 15 can now be easily solved in a least-squares sense to find t r  and d :  ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ −=⎥⎦ ⎤⎢⎣ ⎡ − N tt c c GGG d t M r 1 1)( . (16)  3. EXPERIMENTS AND RESULTS The calibration experimental setup consists of SonixTouch ultrasound machine (Ultrasonix Medical Corporation, Richmond, BC, Canada), L14-5 10 MHz linear 2D ultrasound transducer, Optotrak Certus optical tracker (Northern Digital Inc, Waterloo, Ontario, Canada) and an aluminum plate immersed in a water tank.  The tracker identifies the pose of the ultrasound transducer by locating the infrared markers attached to the transducer. Ultrasound images of the aluminum plate (wall) for 40 different poses of the transducer along with the pose of the transducer have been recorded. A combination of rotational and translational movements has been applied to the transducer to generate the poses as suggested in5. The automatic line detection algorithm in5 has been implemented to find the lines in the images. Using the line slopes and the described method, the rotation part of the calibration transform is first calculated. To make the solution more robust, the unity constraints of 0U r , 0V r  and nr  have been considered as regularization terms while solving the set of linear equations in Eq. (10) with an iterative solver. The condition number is less than 1000; indicating good conditioning. To solve the translation part using Eq. (14), first the values of ic  are found from calculated iU r , iV r  and nr from the previous part and using the measured y-intercept of the lines. The solution of Eq. (14) gives the translation part ( t r ) and the wall to origin distance ( d ). The accuracy of the proposed method has been compared to the conventional single wall technique using the same data set from the line extraction software. As suggested in5, the equations are solved using the Levenberg-Marquardt method. The residual error of the calibration (in Eq. 5) for the proposed method is 0.14 mm while the residual error for the conventional wall calibration technique is 0.91 mm. Standard deviation of the calibration parameters for 20 trials for both methods have been summarized Table 1. Table 1.  Standard deviation of the calibration parameters for 20 trials, each with 40 ultrasound images of the phantom.  )(mmtx  )(mmty  )(mmtz  )(°xr  )(°yr  )(°zr Conventional single wall technique 1.3949    2.1482    1.0877 1.1510     2.2580 0.6764 Closed-form differential single wall technique 0.0714 0.0998 0.2393 0.0507 0.118 0.4412  Proc. of SPIE Vol. 8316  83162A-5 Downloaded from SPIE Digital Library on 05 Mar 2012 to Terms of Use:  http://spiedl.org/terms  4. CONCLUSION For the first time, a novel closed-form solution for the conventional single wall technique is proposed. Another novel idea is using differential measurements to perform calibration. By performing calculations on differential measurements (line slope) using the closed-form solution more accurate results are obtained compared to the conventional single wall technique5. The closed-form solution can give insight for understanding error propagation and for finding the optimal set of transducer poses. Also the closed-form solution provides fast convergence to the solution and is useful to avoid wrong sub-optimal solutions. Accuracy can be improved further by using RF data instead of B-mode images especially since it is easy to find line slopes by cross correlation of any two RF echo pulses. Accurate measurement is possible by using new methods of ultrasound motion tracking8. This is the future goal of this work. 5. ACKNOWLEGMENTS This work is supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) and Canadian Institutes of Health Research (CIHR). REFERENCES [1]  Hsu, P., Prager, R., Gee, A., and Treece, G., “Freehand 3D Ultrasound Calibration: A Review,” Advanced Imaging in Biology and Medicine, 47-84 (2009). [2]  Mercier, L., Langø, T., Lindseth, F., and Collins, D.L., “A review of calibration techniques for freehand 3-D ultrasound systems,” Ultrasound in Medicine & Biology 31(4), 449-471 (2005). [3]  Viswanathan, A., “Immediate Ultrasound Calibration from Three Poses and Minimal Image Processing,” Proceedings in Lecture Notes in Computer Science Medical Image Computing and Computer-Assisted Intervention (MICCAI) 3217, 446-454 (2004). [4]  Chen, T.K., Thurston, A.D., Ellis, R.E., and Abolmaesumi, P., “A Real-Time Freehand Ultrasound Calibration System with Automatic Accuracy Feedback and Control,” Ultrasound in Medicine & Biology 35(1), 79-93 (2009). [5]  Prager, R., “Rapid calibration for 3-D freehand ultrasound,” Ultrasound in Medicine & Biology 24(6), 855-869 (1998). [6]  Boctor, E.M., Iordachita, I., Choti, M.A., Hager, G., and Fichtinger, G., “Bootstrapped ultrasound calibration,” Studies in Health Technology and Informatics 119, 61-66 (2006). [7]  Afsham, N., Chan, K., Pan, L., Tang, S., and Rohling, R.N., “Alignment and calibration of high frequency ultrasound (HFUS) and optical coherence tomography (OCT) 1D transducers using a dual wedge-tri step phantom,” Proceedings of SPIE 7964(1), 796428-796428-8 (2011). [8]  Zahiri-Azar, R., and Salcudean, S.E., “Motion Estimation in Ultrasound Images Using Time Domain Cross Correlation With Prior Estimates,” IEEE Transactions on Biomedical Engineering 53(10), 1990-2000 (2006).   Proc. of SPIE Vol. 8316  83162A-6 Downloaded from SPIE Digital Library on 05 Mar 2012 to Terms of Use:  http://spiedl.org/terms


Citation Scheme:


Usage Statistics

Country Views Downloads
China 11 3
United States 2 0
City Views Downloads
Beijing 11 0
Mountain View 1 0
Sunnyvale 1 0

{[{ mDataHeader[type] }]} {[{ month[type] }]} {[{ tData[type] }]}


Share to:


Related Items