accurate detection of cardiac quiescence is crucial for many cardiac imaging modalities including cardiac computed tomography ( ct ) and magnetic resonance imaging ( mri ) . to minimize motion artifacts and blurring in cardiac images , it is necessary to acquire image data during quiescent phases of the cardiac cycle , i.e. , when the heart is quasi - stationary . current methods for predicting quiescent phases , and in turn triggering image data acquisition , rely almost exclusively on the electrocardiogram ( ecg ) . as an indirect representation of the mechanical state of the heart , ecg has been shown to be imprecise in predicting quiescent phases due to subject and heart - rate variability , . it is therefore desirable to detect quiescence from a mechanically derived signal of cardiac activity . for this work , cardiac quiescence will be detected from b - mode echocardiographic cine data using two frame - to - frame deviation - based methods . accurate gating ensures that data acquisition occurs during periods of cardiac quiescence , consequently minimizing blurring in the resulting ct images . in addition , accurate gating is necessary to ensure that the heart is in the same position in the image plane for techniques requiring data acquisition over multiple cardiac cycles . ctca is a promising alternative to the much more invasive catheter coronary angiography ( cca ) . more than 1.1 million ccas are performed annually in the united states of america costing approximately$40 billion and resulting in more than 14,000 major complications . furthermore , nearly 40% of these tests reveal no coronary artery disease ( as defined by less than 20% vessel stenosis ) . conversely , ctca does not require catheterization and is thus a less invasive and cost - effective alternative . echocardiography , i.e. , ultrasound of the heart , will be used as the modality to investigate quiescence due to its real - time nature , high temporal resolution , accessibility , and lack of ionizing radiation . in particular , because echocardiography directly provides cardiac motion information , it allows for the accurate assessment of cardiac quiescence . b - mode echocardiographic data consists of a sequence of two - dimensional images , allowing for the application of image processing methods , and provides insight into the in - plane two - dimensional motion of the heart as a function of time . echocardiographic tissue doppler imaging ( tdi ) is often used for measuring the velocity of the myocardium in the direction of the ultrasound interrogation . tdi relies on the doppler shift of the received ultrasonic pulses to calculate tissue velocity . while tdi is an efficient , straightforward measurement of cardiac motion , it is only sensitive in the direction of the pulse and therefore can not capture the two - dimensional motion that can be derived from b - mode data . in order to better understand and more accurately predict quiescent phases , the interventricular septum ( ivs ) will be observed using active contour based tracking and frame - to - frame deviation measurement techniques . the ivs was chosen because it has been shown to be an accurate indicator of cardiac and , importantly for ctca , coronary vessel quiescence . the remainder of this paper is organized as follows . section  ii starts with a description of the frame - to - frame deviation measures used to identify quiescent phases . then , the active contour techniques used to perform feature segmentation and tracking are provided . in section iii , the quiescence detection methods are verified and compared using bland - altman analysis . lastly , a discussion on the application of the methods and the implications of this work on cardiac imaging is given in section  iv . for this work , ecg and echocardiographic data were synchronously acquired from 10 human subjects with no known cardiac conditions ( 4 male , 6 female , 2345 years ) using an ultrasonix sonixtouch ultrasound machine ( analogic , peabody , ma , usa ) . emory university 's institutional review board approved the subject evaluations , and full , written , informed consent was obtained from each subject . the ivs was observed from the apical four - chamber view . for all but the first three subjects , the first two subject datasets consist of two 30-second intervals recorded pre- and post - exercise without breath hold . this choice was originally made to observe the effect of heart rate on cardiac quiescence . ecg data were acquired at a rate of 200 hz while b - mode data were acquired at frame rates of either 30 or 50 frames per second ( fps ) . if the heart is stationary for some period of a b - mode sequence , the image content of each frame during that period will be very similar . in this sense , frame - to - frame dissimilarity , deviation , can be associated with movement and hence consecutive frames with low deviation correspond to cardiac quiescence . for this work , deviation between frames will be expressed as a negative function of the pearson correlation coefficient taken over the ivs as it has been shown to be a strong indicator of image similarity . here , the pearson correlation can be thought of as a measure of the similarity between two frames and is defined as \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ $ \eqalignno{&\rho\left(i_{i},i_{j}\right)\cr & = { { \sum\nolimits_{\left(x , y\right)\in i}{\left(i_{i}\left(x , y\right)\!-\!\bar{i_{i}}\right)\left(i_{j}\left(x , y\right)-\bar{i_{j}}\right)}}\over{\sqrt{\sum\nolimits_{\left(x , y\right)\in i}\left(i_{i}\left(x , y\right)\!-\!\bar{i_{i}}\right)^{2}\sum\nolimits_{\left(x , y\right)\in i}\left(i_{j}\left(x , y\right)\!-\!\bar{i_{j}}\right)^{2}}}},&{\hbox{(1)}}}$$\end{document } where \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i_{i}$\end{document } is frame \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i$\end{document } , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ \bar{i_{\imath}}$\end{document } is the mean of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i_{i}$\end{document } , and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ ( x , y)$\end{document } indexes each frame . the correlation - based frame - to - frame deviation measure was calculated from a user - selected static rectangular region containing the ivs . for this work the region the deviation calculation of the region of interest ( roi ) is summarized as \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ $ \eqalignno{&{{\bf d}}\left(i , j\right)=1-\rho ( i_{i},i_{j})\cr & = 1\!-\!{{\sum\limits_{\left(x , y\right)\in s}{\left(i_{i}\left(x , y\right)\!-\!\bar{i_{\imath}^{s}}\right)\left(i_{j}\left(x , y\right)\!-\!\bar{i_{j}^{s}}\right)}}\over{\sqrt{\sum\limits_{\left(x , y\right)\in s}\left(i_{i}\left(x , y\right)\!-\!\bar{i_{i}^{s}}\right)^{2}\sum\limits_{\left(x , y\right)\in s}\left(i_{j}\left(x , y\right)\!-\!\bar{i_{j}^{s}}\right)^{2}}}},&{\hbox{(2)}}}$$\end{document } where \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i_{i}$\end{document } is frame \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i$\end{document } , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ s$\end{document } is the constant rectangular image region over which correlation is calculated , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ \bar{i_{i}^{s}}$\end{document } is the mean of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i_{i}$\end{document } in \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ s$\end{document } , and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ x$\end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ y$\end{document } index image \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i_{i}$\end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i_{j}$\end{document}. for this work , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ s$\end{document } is chosen so that the cardiac feature of interest is contained in \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ s$\end{document } for all frames . this is accomplished by observing the beginning of the b - mode sequence for at least one cardiac cycle . frame - to - frame deviation measurements for all possible comparisons of frame pairs for a given b - mode sequence are calculated according to ( 2 ) and recorded as a deviation matrix , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}\left(i , j\right)$\end{document } , where \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i$\end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ j$\end{document } are the time indices of the frames . when viewed as an image , this allows for the frame - to - frame deviation of the entire sequence to be efficiently observed , as shown in fig . 1.graphical representation of quiescence detection from a deviation matrix , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { \bf d}$\end{document } ( a ) , and the corresponding quiescence plot ( b ) . one cardiac cycle is shown with a period of 1.2 seconds . quiescent phases are indicated by dotted white squares where the mean value of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { \bf d}$\end{document } was chosen as the threshold . graphical representation of quiescence detection from a deviation matrix , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { \bf d}$\end{document } ( a ) , and the corresponding quiescence plot ( b ) . one cardiac cycle is shown with a period of 1.2 seconds . quiescent phases are indicated by dotted white squares where the mean value of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { \bf d}$\end{document } was chosen as the threshold . , this corresponds to low deviation between neighboring frames in the cardiac sequence . as such , square regions of low deviation along the diagonal of the deviation matrix , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , j)$\end{document } , indicate cardiac quiescent phases . for each point along the diagonal of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}\left(i , j\right)$\end{document } , the mean deviation measure of a square neighborhood centered at that point is calculated . the size of this neighborhood is increased until the mean deviation is no longer below a specified threshold , the mean of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}$\end{document } in this case . here the quiescence , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ q(i)$\end{document } , for each time instance , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i$\end{document } , is defined as the duration of the gating window centered at \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i$\end{document } indicated by the width of the largest square region centered at \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , i)$\end{document } with mean deviation below a specified threshold ( fig .  1(b ) ) . this method can be seen as an extension of the one - dimensional , m - mode method presented in   to two - dimensional b - mode data . the deviation measure , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , j)$\end{document } , can be shown to be nearly linear with feature displacement assuming the feature is rigid and the motion linear between frames \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ i$\end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ j$\end{document } , as described in the appendix . given a b - mode sampling period of 20 to 33 ms , these assumptions are reasonable for neighboring frames and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , i+1)$\end{document } is a suitable approximation of the velocity magnitude of the ivs within a scaling factor and constant offset . if necessary , the approximation of the velocity from \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}$\end{document } , defined as \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ $ v_{{\bf d}}(i)={{\bf d}}\eqno{\hbox{(3)}}$$\end{document } can be fit to a reference velocity provided by the active contour tracking methods described by ( 6 ) . alternatively , the approximated signals can be normalized and the velocity magnitude expressed in arbitrary units . interpolated velocity maps expressed in arbitrary units were calculated ( see fig .  2 ) to assist in the visual identification of cardiac quiescent phases . the velocity of each cardiac cycle , segmented by the r - peaks of the synchronously acquired ecg signal , is normalized to have a minimum of zero and a mean of one . this choice was made to eliminate the dependence on the maximum velocity , which was observed to vary more than the minimum and mean between cycles and acquisitions . cubic interpolation was used to fit the data to a uniform grid , allowing velocity to be displayed as a function of heart rate and cardiac cycle phase . an example of the interpolated velocity map for subject 4 is shown in fig .  2 , with the blue regions corresponding to minimal velocity and cardiac quiescence . 2 , the systolic and diastolic quiescent phases can be readily identified and the relative durations of these phases with respect to heart rate can be observed . 2.surface plot of the interpolated velocity from \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document } for all cardiac cycles observed for subject 4 . cardiac cycle percentage is shown on the x - axis , instantaneous heart rate on the y - axis , and velocity magnitude on the z - axis . note that , as expected , the diastolic quiescent phase centered at approximately 75% of the cardiac cycle decreases in duration as heart rate increases . surface plot of the interpolated velocity from \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document } for all cardiac cycles observed for subject 4 . cardiac cycle percentage is shown on the x - axis , instantaneous heart rate on the y - axis , and velocity magnitude on the z - axis . note that , as expected , the diastolic quiescent phase centered at approximately 75% of the cardiac cycle decreases in duration as heart rate increases . for this work , active contour tracking is used to validate the frame - to - frame deviation measures in regards to approximating the velocity magnitude of the ivs . active contour methods provide a flexible framework for feature - based segmentation from image data . active contour evolution relies both on constraints placed on the shape of the contour by the user and on the image characteristics along the contour . the balance of these two effects is tuned based on application and content of the image or sequence . classical active contour methods rely on discretizing the contour as a series of vertices , often referred to as snaxels , resulting in a contour \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ c(p)=\langle x(p ) , y(p)\rangle$\end{document } where \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ p$\end{document } is the contour length parameterized and normalized to one . the solution contour is found by minimizing the energy of the contour expressed as \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ $ e_{contour}=\int_{0}^{1}\left(e_{int}\left(c\left(p\right)\right)+e_{img}\left(c\left(p\right)\right)\right)dp,\eqno{\hbox{(4)}}$$\end{document } where \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ e_{int}$\end{document } is the energy associated with the shape of the contour and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ e_{img}$\end{document } is the energy associated with the image content . \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ e_{img}$\end{document } is most often chosen to be a decreasing function of the gradient of the image in order to attract the contour to edges in the image . extending active contour models to a sequence of frames is generally accomplished by finding the active contour on a frame - by - frame basis with each frame being solved with knowledge of only the current and previous frame . the solution of the previous frame is used as an estimate for the contour position in the following frame . different schemes , such as optical flow   and kalman filtering , are often employed to improve the estimation of the contour position in the next frame . active contours for this work are solved using a method based on , where optical flow is calculated according to   for the b - mode sequence and used in conjunction with the previous active contour solution to estimate the position of the contour in the subsequent frame . optical flow attempts to estimate frame - to - frame movement by defining a small neighborhood around each pixel in the current frame and then locating this neighborhood in the following frame . the cardiac - feature tracking algorithm identifies the feature of interest in each frame of b - mode data . the algorithm is initialized with the user providing an approximate outline consisting of 10 points along the boundary of the cardiac feature in the first frame of the sequence . the rough contour defined by these points is upsampled by five resulting in the evolution of a contour defined by 50 discrete points . the numerical active contour technique presented in   is used to find the outline of the feature in each frame by minimizing ( 4 ) . for all frames after the first , the position of the previous contour , along with optical flow data of the sequence , is used to initialize the active contour algorithm . knowledge of the velocity magnitude of the cardiac anatomic feature of interest can be used to validate the approximation provided by the frame - to - frame deviation methods . the velocity magnitude can be calculated from the two - dimensional velocity of the feature obtained by differencing the position of the contour in each frame . the position is calculated for each frame by finding the centroid of the cardiac feature , defined as the average location of the pixels inside the contour , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ $ \bar{x}_{i}={{1}\over{n}}\sum_{x\in s_{i}}x,\quad\bar{y}_{i}={{1}\over{n}}\sum_{y\in s_{i}}y,\eqno{\hbox{(5)}}$$\end{document } where \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ s_{i}$\end{document } is the set of all pixel locations inside the contour \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf x}}_{i}$\end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ n$\end{document } is the number of elements in \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ s_{i}$\end{document}. the velocity magnitude of the feature can then be found by approximating the two - dimensional velocity obtained by taking the first difference of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ \bar{x}_{i}$\end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ \bar{y}_{i}$\end{document}. this process is summarized by \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ $ v_{ac}(i)=\sqrt{\left(\bar{x}(i)-\bar{x}\left(i-1\right)\right)^{2}+\left(\bar{y}(i)-\bar{y}\left(i-1\right)\right)^{2}},\eqno{\hbox{(6)}}$$\end{document } with periods of low velocity magnitude corresponding to phases of cardiac quiescence . to assess the efficacy of b - mode detection of cardiac quiescence as well as the characteristics of cardiac quiescent phases for a small healthy population ( 10 subjects , 2345 years old ) the following results are provided . first , active contour tracking was validated in the direction of ultrasonic interrogation using synchronously acquired tdi data . next , the velocity magnitude of the ivs from the deviation measure , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document } , was confirmed using active contour tracking . then , bland - altman analysis was used to compare the quiescent phases indicated by the square regions of the deviation matrix , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , j)$\end{document } , and the velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}\left(i\right)$\end{document}. lastly , the nature of the quiescent phases for the ten subjects was observed from \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document}. the ivs was successfully segmented and tracked for seven of the 10 subjects using the active contour algorithm presented in section  ii . for those not tracked , poor acoustic windows , acoustic shadowing , and out of plane motion an example of a solved active contour for a given frame is provided in fig . lv and rv are the left and right ventricle , and la and ra are the left and right atrium . the position of the ivs was found for each frame using the centroid calculation given in ( 5 ) . from the \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ x$\end{document}- and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ y$\end{document}-position vectors , the velocity magnitude was calculated using first - difference methods as described in ( 6 ) . 3.b - mode frame obtained from an apical four - chamber view of the heart , with solved contour shown around the ivs . the four chambers : right ventricle ( rv ) , right atrium ( ra ) , left ventricle ( lv ) , and left atrium ( la ) are shown . b - mode frame obtained from an apical four - chamber view of the heart , with solved contour shown around the ivs . the four chambers : right ventricle ( rv ) , right atrium ( ra ) , left ventricle ( lv ) , and left atrium ( la ) are shown . to verify the accuracy of the two - dimensional active contour tracking , the velocity calculated from the active contour centroids was compared to the velocity from synchronously acquired tdi data . as a direct measurement of tissue velocity based on the doppler shift of the received ultrasonic pulses , tdi is often used to measure the one - dimensional velocity of cardiac features in the direction of ultrasonic interrogation . active contour based velocity in the direction of ultrasonic interrogation , axial direction , was found to agree with that of tdi , indicating that active contours provide accurate velocity information . in addition , the active contour based velocity in the lateral direction , orthogonal to the direction of interrogation and not detected by tdi , was found to be significant . this suggests that b - mode derived two - dimensional velocity provides more complete ivs motion information than tdi . clearly , the magnitude of the tdi velocity will be substantially different than the actual magnitude of the velocity , which is approximated by the velocity derived by active contours . 4.velocity in the axial ( a ) and lateral ( b ) directions are shown with the active contour velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{ac}$\end{document } , and tdi velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{tdi}$\end{document } , in blue and green , respectively . velocity in the axial ( a ) and lateral ( b ) directions are shown with the active contour velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{ac}$\end{document } , and tdi velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{tdi}$\end{document } , in blue and green , respectively . the deviation matrices , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , j)$\end{document } , for each subject were calculated and the velocity of the ivs was approximated by the one - off diagonal of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , j)$\end{document } , as described in section  ii . note that frame - to - frame deviation methods could be used for all 10 subjects as opposed to active contour methods that could only be used for seven subjects . to verify the accuracy of the velocity approximation , a number of velocity signals , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document } from ( 3 ) , were compared to the velocity signals derived from active contour tracking , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{ac}(i)$\end{document } from ( 6 ) . 5 shows an example of this fitting and demonstrates that the velocity derived from the frame - to - frame deviation measure matches that derived from contour tracking . in addition , the velocity approximation based on correlation was observed to have less noise in the quiescent phases with low velocity magnitude . this is most likely due to noise that is accentuated by the differentiation of digital signals . 5.comparison of active contour velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{ac}$\end{document } , and correlation based velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document } , in blue and green , respectively . note that the correlation based velocity approximation appears to have less noise in the quiescent phases with low velocity magnitude . comparison of active contour velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{ac}$\end{document } , and correlation based velocity , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document } , in blue and green , respectively . note that the correlation based velocity approximation appears to have less noise in the quiescent phases with low velocity magnitude . quiescent phases were determined as regions of the velocity signal , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document } , less than the mean of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document}. quiescence plots , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ q(i)$\end{document } , were then calculated by finding square regions along the diagonal of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , j)$\end{document } with average deviation under a specified threshold . the mean of \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { \bf d}$\end{document } was used as the threshold for quiescence . the center and duration of the diastolic quiescent phase for each cardiac cycle were then determined for the velocity and corresponding quiescence signals , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}(i)$\end{document } and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ q(i)$\end{document}. for the quiescence signals the peak corresponding to the diastolic quiescent phase was found with the position and height indicating the center and duration of the quiescent phase , respectively . both methods for identifying quiescent phases quiescent phases can be identified from both ( a ) the velocity signal , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document } , and ( b ) the quiescence , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ q(i)$\end{document } , from square regions of the deviation matrix . quiescent phases can be identified from both ( a ) the velocity signal , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document } , and ( b ) the quiescence , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ q(i)$\end{document } , from square regions of the deviation matrix . the center and duration of each diastolic quiescent phase detected from the velocity and quiescence signals were then compared for all subjects using bland - altman analysis to verify agreement between the measurement methods . as shown in fig . 7.comparison of diastolic quiescent phase identification methods for all cardiac cycles identified from both the velocity signal , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document } , and the square regions of the deviation matrix . the difference between the two measures is plotted on the vertical axis and the mean of the two measures is plotted on the horizontal axis . the cardiac cycle percentage of the phase centers ( a ) and durations ( b ) are plotted . comparison of diastolic quiescent phase identification methods for all cardiac cycles identified from both the velocity signal , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document } , and the square regions of the deviation matrix . the difference between the two measures is plotted on the vertical axis and the mean of the two measures is plotted on the horizontal axis . the cardiac cycle percentage of the phase centers ( a ) and durations ( b ) are plotted . due to the similarity of the velocity and quiescence signals in determining quiescent phases , the velocity signal was chosen to display the analysis of the data due to computational efficiency , requiring only one correlation measure per time index . 8 where the left column of plots indicates all quiescent phases observed longer than 83 ms for four subjects and the right column of plots depicts the interpolated velocity maps of those subjects . the cutoff of 83 ms was chosen as it corresponds to the minimum imaging time for a dual - source ct machine with a gantry rotation time of 333 ms . as observed in practice , at high heart rates the systolic quiescent phase becomes proportionally longer than the diastolic phase , . 8(b ) where the duration of the blue region in systole corresponding to minimal velocity becomes longer than that in diastole as heart rate increases . 8.scatter plots ( a ) and velocity color maps ( b ) shown for four subjects . for the scatter plots the dots correspond to the quiescent phase centers and the bars correspond to the durations as determined using the velocity signal , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document}. for the color maps blue corresponds to low velocity magnitude and red corresponds to high velocity magnitude . scatter plots ( a ) and velocity color maps ( b ) shown for four subjects . for the scatter plots the dots correspond to the quiescent phase centers and the bars correspond to the durations as determined using the velocity signal , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ v_{{\bf d}}$\end{document}. for the color maps blue corresponds to low velocity magnitude and red corresponds to high velocity magnitude . the diastolic quiescent phase statistics for each of the 10 subjects are provided in table  i. because the center and duration of the diastolic quiescent phase varies with heart rate , care should be taken in interpreting results for subjects who had a wide range in observed heart rates , indicated by a high heart rate standard deviation . from the table , a significant amount of inter- and intra - subject variability can be observed , suggesting that gating based solely on the ecg is suboptimal for predicting the timing of cardiac quiescent phases and that personalized gating protocols could improve this prediction . table isummary of diastolic quiescent phase statisticsquiescent phase centerquiescent phase durationheart ratecyclesmeanstd devmeanstd devmeanstd devsubject 174 % 6.4 % 451 ms203 ms60 bpm10.7 bpm61subject 270%4.4%352 ms113 ms58 bpm5.7 bpm39subject 375 % 4.3 % 203 ms40 ms99 bpm5.9 bpm10subject 474%1.4%324 ms33 ms58 bpm2.3 bpm37subject 579%6.6%246 ms107 ms63 bpm1.5 bpm31subject 686 % 4.8 % 188 ms71 ms83 bpm2.1 bpm43subject 782%2.2%182 ms32 ms77 bpm2.8 bpm42subject 884 % 3.4 % 301 ms80 ms68 bpm2.0 bpm39subject 975%2.4%446 ms73 ms59 bpm2.3 bpm35subject 1072 % 8.1 % 150 ms87 ms66 bpm2.9 bpm23 these methods are based on frame - to - frame deviation measures of a user - specified image region and were verified using active contour tracking . the two - dimensional velocity derived from active contour tracking was validated in the axial direction using tdi . in addition , the lateral velocity component , not detected by tdi , was shown to be significant . this suggests that two - dimensional , image - based methods may prove more accurate in determining cardiac quiescence than the more straightforward tdi method of determining myocardial velocity . these two methods for determining cardiac quiescence were used to analyze data for 10 subjects with no known cardiac conditions . the first method exploits the fact that the correlation between two frames is linearly related to feature displacement , given an assumption of linear motion and feature rigidity between the frames in question as described in the appendix . the velocity calculated using the deviation of neighboring frames and verified using active contour tracking suggests that , at least for the ivs , these assumptions are reasonable . the second method for determining quiescence from the deviation attempts to identify quiescent phases as square regions of low deviation from the deviation matrix , \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } $ { { \bf d}}(i , j)$\end{document}. despite using more deviation measures to determine quiescence for each time index , it was observed that this method was actually more susceptible to noise . one possible explanation for this is that the deviation measures of a long sequence of frames allow more time for inadvertent transducer , subject , or breathing motion whereas with the velocity approximation technique , only neighboring frames are compared . despite this , the two methods were found to agree with one another . the proposed deviation - based methods were able to identify quiescent phases for the three subjects where active contour tracking was not possible due to loss of tracking . it is important to note that although quiescent phases can be identified for b - mode sequences with poor imaging , accuracy may be decreased and can not be verified by active contour tracking . in addition to robustness , these methods provide a less noisy velocity signal than that provided by active contour methods , as active contour tracking requires differentiation to obtain the velocity . lastly , these methods are also much more computationally efficient than the active contour methods , requiring only correlation measures instead of solving complicated minimization problems and performing optical flow calculations for each b - mode frame . although not required for the applications in this work , deviation - based velocity calculation could be done in real - time . the choice of a region - based correlation method was made to provide a robust and efficient technique to determine the magnitude of the velocity of a specified roi . speckle tracking is another commonly used method to ascertain motion using optical flow to calculate the frame - to - frame displacement of the pixels inside an roi . for the application presented in this work , the pixel displacements could be averaged over the roi to obtain the velocity magnitude of the cardiac feature in the roi . this could then be used to find quiescent phases . because pixel - level 2-d motion between frames is estimated , more descriptive motion quantification could be performed . however , speckle tracking is not without drawbacks . it requires high image quality and high b - mode frame rate ( 60110 fps ) , placing substantial constraints on the imaging hardware . to minimize the effect of random noise , speckle tracking is usually averaged over at least three cardiac cycles , making speckle tracking for subjects with arrhythmias problematic . lastly , speckle tracking is computationally expensive and as a result is generally performed offline after an examination . another recent technique for robustly quantifying motion leverages the fact that a translation of an object in the image plane corresponds to a phase shift in the fourier domain . by taking the 2-d fourier transform of a sequence of images the 2-d velocity of the feature can be obtained from the phase shifts . because this method relies on the fourier transform and peak detection , it is more computationally expensive than the proposed region - based correlation method one potential concern with determining quiescence from the frame - to - frame deviation is that only the relative velocity magnitude is provided . normalization of the velocity signals partially alleviates this concern because a constant threshold for quiescence can be employed . if needed , absolute velocity magnitude can be obtained either by fitting the deviation - based velocity to a reference velocity , e.g. , the velocity from the active contour tracking , as described by ( 14)(16 ) , or by leveraging the relationship in ( 10 ) . the latter method does not require a reference velocity and would rely on how well the relationship between the area of feature overlap between frames and the feature displacement is known . another potential limitation is that the subject population for this work is not necessarily representative of that expected to receive ctca studies . however , the performance of echocardiography should be independent of subject population because echocardiography is not known to be limited by the age or cardiac health of the patient . therefore , the proposed detection methods should be applicable for patients receiving ctca studies . though liu et al . showed that the ivs is a strong indicator of quiescence of the left anterior descending coronary artery , other arteries may have different quiescent phases . for low heart rates the optimal ctca imaging window ( and consequently the longest quiescent phase ) occurs during the same cardiac phase during diastole for each coronary vessel . for higher heart rates , above 80 bpm , the longest quiescent phase occurs during systole and was shown to vary by 10% depending on the coronary artery in question , while the diastolic quiescent phases were the same . this suggests that for high heart rates the longest quiescent phases determined by observing the ivs may not correspond exactly to those of some coronary vessels . although this may be the case , the relationship between the quiescent phases of the ivs and those of the vessel in question may be deterministic , allowing for the actual quiescent phase to be obtained from that of the ivs . , it is reasonable to suggest that the correlation of neighboring frames can be used to efficiently approximate the magnitude of the velocity of cardiac features , and hence cardiac quiescence . the applicability of this technique to other cardiac features is a topic for further research . in addition , the frame - to - frame deviation of multiple static regions corresponding to different cardiac features in the imaging plane could be used to analyze the relationship between the timing of the quiescent phases of the cardiac features in question . from table i , a significant amount of both inter- and intra - subject variability was observed as evidenced by the different means of the quiescent phase centers and significant standard deviations of those centers for each subject . this suggests that gating based on the mechanics of the heart rather than the ecg may provide better prediction of cardiac quiescent phases for the gating of ctca and other cardiac imaging modalities . the quiescence detection methods presented here will be used for further study into such gating methods . this work shows that the tdi misses a significant amount of velocity in the lateral direction . therefore , it is likely that 2-d b - mode echocardiography motion estimation misses some amount of out - of - plane motion . given 3-d echocardiographic data , the extension of the proposed deviation measures to 3-d is straightforward and a topic for future study . this work addresses the first half of the gating problem , i.e. , quiescence detection . however , the second half of the problem , quiescence prediction , while relying on accurate detection , is more complicated . by knowing the precise timing of cardiac quiescent phases , effective methods for quiescence prediction can be developed and analyzed . the phases detected using the methods proposed in this work can serve as a baseline to compare phases predicted using ct compatible techniques . two potential techniques for cardiac - motion - informed gating are offline and real - time methods . offline methods would focus on obtaining patient specific gating parameters prior to ctca for the standard ecg gating that is currently in use . in this case the amount of potential improvement in images using this approach is a topic for further study . as such , it is necessary for the sensing device used for gating to be ct compatible . unfortunately , standard echocardiographic transducers are not ct compatible because of the extensive streak artifacts that they would cause in the ct images . fortunately , two potential modalities have been shown to be ct compatible : tdi and seismocardiography ( scg ) , . there exist tdi specific transducers that are radiolucent and would not interfere with ct data acquisition . though tdi was shown in this work to be insensitive to a significant amount of two - dimensional velocity , the performance of tdi in predicting quiescent phases was not investigated . scg uses an accelerometer to measure the acceleration of the chest wall due to cardiac motion . modern accelerometers are small enough to not interfere with ct acquisition and as such scg could provide a cardiac - motion - based signal that could be used for predicting cardiac quiescent phases and gating ctca . for both tdi and scg methods of real - time prediction , the methods proposed in this work could be used to assess the predictive accuracy of these methods . in addition to being a direct indication of cardiac motion , signals derived from the mechanics of the heart would improve the quiescence predictive performance for patients with cardiac arrhythmias . this is true because the signal features that can be used for triggering are much nearer to the actual quiescent phase . this minimizes the effect of inaccurate cardiac cycle length prediction caused by arrhythmias . from initial investigations there exist detectable features of the scg that are within 10% of the start of quiescence in terms of cardiac cycle phase . contrasting this against the 70% from the r - peak of the ecg to the diastolic quiescent phase or 40% for the systolic quiescent phase , it is apparent that the predictive performance of mechanical based gating signals would be better than that of the ecg in the presence of cardiac arrhythmias .
two novel methods for detecting cardiac quiescent phases from b - mode echocardiography using a correlation - based frame - to - frame deviation measure were developed . 
 accurate knowledge of cardiac quiescence is crucial to the performance of many imaging modalities , including computed tomography coronary angiography ( ctca ) . 
 synchronous electrocardiography ( ecg ) and echocardiography data were obtained from 10 healthy human subjects ( four male , six female , 2345 years ) and the interventricular septum ( ivs ) was observed using the apical four - chamber echocardiographic view . 
 the velocity of the ivs was derived from active contour tracking and verified using tissue doppler imaging echocardiography methods . in turn , the frame - to - frame deviation methods for identifying quiescence of the ivs were verified using active contour tracking . 
 the timing of the diastolic quiescent phase was found to exhibit both inter- and intra - subject variability , suggesting that the current method of ctca gating based on the ecg is suboptimal and that gating based on signals derived from cardiac motion are likely more accurate in predicting quiescence for cardiac imaging . 
 two robust and efficient methods for identifying cardiac quiescent phases from b - mode echocardiographic data were developed and verified . 
 the methods presented in this paper will be used to develop new ctca gating techniques and quantify the resulting potential improvement in ctca image quality .