Research article Special Issues

Sequential shape similarity for active contour based left ventricle segmentation in cardiac cine MR image

  • Delineation of the boundaries of the Left Ventricle (LV) in cardiac Magnetic Resonance Images (MRI) is a hot topic due to its important diagnostic power. In this paper, an approach is proposed to extract the LV in a sequence of MR images. In the proposed paper, all images in the sequence are segmented simultaneously and the shape of the LV in each image is supposed to be similar to that of the LV in nearby images in the sequence. We coined the novel shape similarity constraint, and it is called sequential shape similarity (SSS in short). The proposed segmentation method takes the Active Contour Model as the base model and our previously proposed Gradient Vector Convolution (GVC) external force is also adopted. With the SSS constraint, the snake contour can accurately delineate the LV boundaries. We evaluate our method on two cardiac MRI datasets and the Mean Absolute Distance (MAD) metric and the Hausdorff Distance (HD) metric demonstrate that the proposed approach has good performance on segmenting the boundaries of the LV.

    Citation: Ke Bi, Yue Tan, Ke Cheng, Qingfang Chen, Yuanquan Wang. Sequential shape similarity for active contour based left ventricle segmentation in cardiac cine MR image[J]. Mathematical Biosciences and Engineering, 2022, 19(2): 1591-1608. doi: 10.3934/mbe.2022074

    Related Papers:

    [1] Yuan-Nan Young, Doron Levy . Registration-Based Morphing of Active Contours for Segmentation of CT Scans. Mathematical Biosciences and Engineering, 2005, 2(1): 79-96. doi: 10.3934/mbe.2005.2.79
    [2] Xiangfen Song, Yinong Wang, Qianjin Feng, Qing Wang . Improved graph cut model with features of superpixels and neighborhood patches for myocardium segmentation from ultrasound image. Mathematical Biosciences and Engineering, 2019, 16(3): 1115-1137. doi: 10.3934/mbe.2019053
    [3] Chenkai Chang, Fei Qi, Chang Xu, Yiwei Shen, Qingwu Li . A dual-modal dynamic contour-based method for cervical vascular ultrasound image instance segmentation. Mathematical Biosciences and Engineering, 2024, 21(1): 1038-1057. doi: 10.3934/mbe.2024043
    [4] M Nisha, T Kannan, K Sivasankari . A semi-supervised deep neuro-fuzzy iterative learning system for automatic segmentation of hippocampus brain MRI. Mathematical Biosciences and Engineering, 2024, 21(12): 7830-7853. doi: 10.3934/mbe.2024344
    [5] Biao Cai, Qing Xu, Cheng Yang, Yi Lu, Cheng Ge, Zhichao Wang, Kai Liu, Xubin Qiu, Shan Chang . Spine MRI image segmentation method based on ASPP and U-Net network. Mathematical Biosciences and Engineering, 2023, 20(9): 15999-16014. doi: 10.3934/mbe.2023713
    [6] Anna Andrews, Pezad Doctor, Lasya Gaur, F. Gerald Greil, Tarique Hussain, Qing Zou . Manifold-based denoising for Ferumoxytol-enhanced 3D cardiac cine MRI. Mathematical Biosciences and Engineering, 2024, 21(3): 3695-3712. doi: 10.3934/mbe.2024163
    [7] Xinlin Liu, Viktor Krylov, Su Jun, Natalya Volkova, Anatoliy Sachenko, Galina Shcherbakova, Jacek Woloszyn . Segmentation and identification of spectral and statistical textures for computer medical diagnostics in dermatology. Mathematical Biosciences and Engineering, 2022, 19(7): 6923-6939. doi: 10.3934/mbe.2022326
    [8] Tongping Shen, Fangliang Huang, Xusong Zhang . CT medical image segmentation algorithm based on deep learning technology. Mathematical Biosciences and Engineering, 2023, 20(6): 10954-10976. doi: 10.3934/mbe.2023485
    [9] Yingying Xu, Songsong Dai, Haifeng Song, Lei Du, Ying Chen . Multi-modal brain MRI images enhancement based on framelet and local weights super-resolution. Mathematical Biosciences and Engineering, 2023, 20(2): 4258-4273. doi: 10.3934/mbe.2023199
    [10] Jianhua Song, Lei Yuan . Brain tissue segmentation via non-local fuzzy c-means clustering combined with Markov random field. Mathematical Biosciences and Engineering, 2022, 19(2): 1891-1908. doi: 10.3934/mbe.2022089
  • Delineation of the boundaries of the Left Ventricle (LV) in cardiac Magnetic Resonance Images (MRI) is a hot topic due to its important diagnostic power. In this paper, an approach is proposed to extract the LV in a sequence of MR images. In the proposed paper, all images in the sequence are segmented simultaneously and the shape of the LV in each image is supposed to be similar to that of the LV in nearby images in the sequence. We coined the novel shape similarity constraint, and it is called sequential shape similarity (SSS in short). The proposed segmentation method takes the Active Contour Model as the base model and our previously proposed Gradient Vector Convolution (GVC) external force is also adopted. With the SSS constraint, the snake contour can accurately delineate the LV boundaries. We evaluate our method on two cardiac MRI datasets and the Mean Absolute Distance (MAD) metric and the Hausdorff Distance (HD) metric demonstrate that the proposed approach has good performance on segmenting the boundaries of the LV.



    Cardiac Magnetic Resonance (MR) Imaging provides high spatial and temporal resolution images and has been proven to be a powerful imaging modality [1,2]. It can capture the anatomical and functional information of a heart within a cycle [3], and thus, be well suited for the clinical diagnosis and research of cardiovascular diseases (CVDs) for radiologists and physicians. Quantification of the Left Ventricle (LV) of the heart usually acts as index of cardiac health [4]. Therefore, segmentation of the LV from cardiac MR images is an essential prerequisite for the diagnosis. From accurate segmentation results, many clinical diagnosis indices such as wall thickness, myocardial strain, ejection fraction, and circumferential shortening of myocardial fibers can be evaluated [5].

    Clinically, delineating the boundaries of the LV is often performed manually. However, manual segmentation often takes a long time, and is tedious, and subject to high inter-observer variability. Thus, there is a flurry of works devoted to segmenting the LV from MRI in the past decades, such as classification [6,7,8], watershed [9], graph cuts [10], active contour models [11,12,13,14,15,16], level set [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31], active shape models [32,33,34,35,36], and other methods [37,38]. Petitjean [39] and Peng [40] presented a comprehensive review respectively on this topic, interested readers can refer to [39,40] for details and references, and we merely present a brief survey here.

    Active contour models, or the snake models, were one of the most popular tools for segmentation of the region of interesting (ROI), and found successful applications in extracting the LV from MRI. The basic principle remains the same as those presented in the seminal work on the snake model proposed by Kass et al. [41]. Wu et al. coined the circle shape energy to aid the snake model, which considered the LV as a circle, and extracted the endo- and epicardium of the LV by coining the shape similarity constraint [11]. However, the delineation merely relying on the shape constraint still results in mistakes when the discontinuities are found in the object boundaries. The discontinuities may be caused by the papillary muscle or the organs nearby such as the liver. Santarelli also used the snake model for this task [13], in which the Gradient Vector Flow [42] acts as an external force. However, there is no solution for artifacts and weak edges. Ranganath [14] modified the internal energy of the snake by taking into account the curvature of the snake contour to smooth the boundary of the LV. McInerney et al. coined an active surface model to extract the surface of the LV by representing the surface using finite element [15].

    Geometric active contour, or level set, has been proven to be very powerful in a variety of segmentation tasks including the LV extraction. Li et al. coined a bi-layer level set [17] to extract the inner and outer boundaries of the LV, then the bi-ventricle [18]. Ismail et al. manually extract the LV in the first frame in a sequence of images, and supposed the intensity distribution of the myocardium in the first image is similar to the rest of the ones in the sequence [25]. Pluempitiwiriyawej introduced the ellipse-fitting scheme into the level set model and proposed an integrated model for cardiac MR image segmentation [28].

    With the development of deep learning, deep learning-based methods have shown impressive performance in medical image processing tasks. For example, Zhang et al. proposed the BLU-Net [43] to segment LV, Wang et al. combined the DenseNet segmentation and convolution neural network(CNN) regression to quantify full LV metrics [44], Shen et al. proposed a pure dilated residual U-Net [45] to segment the femur and tibia bones. Besides, deep learning combined with traditional methods has also been proved to be useful. Ngo combined deep learning and level set [46] to segment the left ventricle of the heart from cardiac cine MR data automatically, in which the deep learning method is responsible for locating the LV, and the level set is responsible for refining the boundary of the LV. Avendi et al. combined deep learning algorithms with deformable models [47] to evaluate automatically the LV function from short-axis cardiac MRI datasets.

    In this study, we propose an extension of our previous work in [11] based on the snake model. For a MR image sequence, it is often sampled from a cardiac cycle, and it is noticed that the shape of the LV in one image is similar to that of the LV in nearby ones in the sequence. In the aforementioned works, segmentation is carried out in an image-by-image manner, while in this work, we propose to segment all images in the sequence simultaneously and the shape of the LV in each image is supposed to be similar to that of the LV in nearby images in the sequence. We coined the novel shape similarity constraint, and it is called sequential shape similarity (SSS in short). Our approach integrates the prior knowledge from the image sequence into the snake model.

    The organization of this paper is as follows. Section II presents a brief introduction of our previous work in [11]. In Section III, the proposed sequential shape similarity (SSS) is detailed. In Section IV, the experimental results are reported, and the conclusion marks are presented in the final Section.

    We first made a brief review of our previous work [11], which include the snake model, GVC external force, Circle Shape Constraint (CSC) for endocardium segmentation, shape similarity constraint (SSC), and novel edge map for epicardium.

    Active contour models or snakes, coined by Kass et al. in 1988 [41], are dynamic curves c(s)=(x(s),y(s)),s[0,1] that evolve by minimizing the following energy functional,

    Esnake=10[12(α|c'(s)|2+β|c''(s)|2)+Eext]ds (1)

    where c'(s) and c''(s) are the first and second derivatives of c(s) with respect to arc length s. α and β are the weights controlling the smoothness and rigidity of the curve, respectively. Eext is derived from the image data and takes the smallest values at boundaries. The snake models gained great success in image processing.

    However, the traditional snakes are limited by the external force based on image gradient. Therefore, the gradient vector flow (GVF) for snakes is proposed by Xu and Prince [42]. The GVF is effective, yet, with high computation costs. To overcome this shortcoming, we previously proposed a new external force called Gradient Vector Convolution(GVC) [48] based on convolution and fast Fourier transform (FFT). The GVC can be calculated as follows:

    {u(x,y)=Kconfxv(x,y)=Kconfy (2)

    where Kcon is the convolution kernel, is the convolution operator. The GVC model behaves similarly to the well-known GVF model on enlarging capture range and robustness to initialization. Moreover, using the FFT, the GVC can be implemented in real-time since it is expressed in convolution, and thus, the snake model can also benefit from this fact. In practice, we take Kcon=1/(rh)n, where rh=x2+y2+h hR+,nR+. Kcon always shows high performance in terms of extending and smoothing gradient vectors. Thereby, the GVC snake can preserve edges and dive into deep concavities.

    For the segmentation of the endo- and epicardium of the LV, the difficulties originate primarily from the gaps or the weak boundaries caused by the artifacts and papillary muscles. In the seminal work in Eq (1), the internal energy, which consists of first and second derivatives, controls the continuity and smoothness of the contour, respectively. However, the first and second derivatives are merely local geometrical properties of the contour and if there were weak edges or gaps, the snake contour would fail to bridge the gaps and leak out. Moreover, if there is a local minimum caused by the external force, the snake contour would be trapped. To address these issues, the information about the overall shape of the endo- and epicardium should be incorporated into the snake energies. The shape information will guide the snake contour to delineate the target object correctly. The circle shape constraint (CSC), which roughly regards the shape of the endocardium as a circle, is adopted for the segmentation. This shape constraint tries to keep the contour being a circle during the movement. The circle shape constraint can be represented as follows,

    ECSC=λ210(R(s)R)2ds (3)

    where R(s)=(x(s)xc)2+(y(s)yc)2, R=10R(s)ds, xc=10x(s)ds, yc=10y(s)ds. (xc,yc) is the center of the snake contour. The energy function in Eq (3) characterizes the deviation of the snake contour from a circle of radius R located at (xc,yc). Both radius R and (xc,yc) change along with the motion of the snake contour. If the snake contour is attracted by some nuisances such as the or papillary muscle, tit is still approximately a circle under this constraint, which is the global shape of the LV.

    For segmentation of the epicardium, we proposed the Shape Similarity Constraint (SSC) that is based on the observation of shape similarity between the endocardium and the epicardium. The segmented endocardium is referred to as a prior shape constraint to extract the epicardium. The SSC constraint is formulated as,

    ESSC=η210((R(s)R)(r(s)r))2ds (4)

    The parameters in Eq (4) are meaningly similar to those in Eq (3), but R and R are for epicardium and r and r are for segmented endocardium.

    After the endocardium is segmented, the result is employed as the initial contour to segment the epicardium contour. However, the endocardium is merely a local minimum of the external energy derived from GVC (the red circle approximates the endocardium in the upper-right of Figure 1), and it would prevent the contour from converging towards the epicardium (the yellow circle approximates the epicardium in the upper-right of Figure 1. To address the issue, the original edge map is corrected by setting to zero the original edge map around and within the endocardium. This new edge map is used to calculate a new external force by GVC, and the endocardium is no longer the local minimum, and the new external force can directly put the snake contour forward to the epicardium.

    Figure 1.  Comparison of external forces for epicardium segmentation. The upper row: the original edge map and force field. The lower row: the modified edge map and force field.

    Let I be an MR image sequence containing K images. Suppose that there are 20 samples in a cardiac cycle, and the number K is equal to 20. Generally speaking, the purpose of the LV segmentation is to track the endocardium Rendok and epicardium Repik in each image Ik, k=1,,K. The aforementioned SSC can be described as follows, Rendok and Repik in image Ik are generally similar in shape and when segmenting Repik, the segmented Rendok is referred to as the prior shape constraint. For SSC, only the similarity in a single image is integrated into the snake model and some useful information among the images in a sequence is discarded. For example in Figure 2, there are gaps in the myocardium in Figure 2(b), see the location in the red circle, and it is hard to extract the myocardium in an image-by-image manner. However, in Figure 2(a), (c), the myocardium is intact, therefore, the shape of the myocardium in Figure 2(a), (c) can be employed to construct the constraint to segment the endo- and epicardium in Figure 2(b).

    Figure 2.  Three samples in an image sequence.

    Based on the observation in Figure 2, we propose the sequential shape similarity (SSS in short). For an image Ik (say the one in Figure 2(b)) and its neighbors Ik1 (the one in Figure 2) and Ik+1 (the one in Figure 2(c)) in sequence I, the shape of the endocardium Rendok is similar to Rendok1 and Rendok+1. Thus when segmenting Rendok, besides the CSC in Eq (3), Rendok1 and Rendok+1 can be also referred to as the prior shape constraint. Furthermore, the similarity range can be extended from Rendok±1 to Rendok±m, m[1,K/2]. For Rendok, there are 2 m endocardium in the sequence to construct a more effective similarity constraint. For Repik, besides the SSC in Eq (4), Repik1 and Repik+1 can be also referred to as the prior shape constraint. The proposed constraint is called sequential shape similarity (SSS). For clarification, we neglect the upper indices (epi, endo) in the derivation of the SSS, and the SSS can be represented as follow

    EkSSS=mj=mρj2((Rk(s)Rk)(Rk+j(s)Rk+j))2ds (5)

    where m[1,K/2]. The parameters in Eq (5) are meaningly similar to those in Eq (4), but Rk+j(s) and Rk+j are for image Ik+j and Rk(s) and Rk are for image Ik. In Eq (5), it is clear that Rk(s)Rk characterize the deviation of the snake contour from a circle of radius Rk in Ik and Rk+j(s)Rk+j measures the deviation of the snake contour from a circle with a radius Rk+j in Ik+j. Minimizing EkSSS makes the snake contour resembling the shape of LV in all the images in sequence I.

    Suppose there are n discrete points on the snake contour, Rk and Rk in Eq (5) can be estimated by Rk=1nni=1Rki, Rki=(xkixkc)2+(ykiykc)2, and xkc=1nni=1xki, ykc=1nni=1yki, i=1,2,,n, (xkc,ykc) is the centroid of the contour in Ik. By calculus of variation, the discretized equation of Eq (5) reads,

    {mj=mρj(xkixk+ji(xkcxk+jc)(RkRk+j)cos(2πi/n))=0mj=mρj(ykiyk+ji(ykcyk+jc)(RkRk+j)sin(2πi/n))=0 (6)

    here, the following function is employed to determine the weight ρj,

    ρj=aej2/b2,j[m,m] (7)

    ρj is to weigh the degree of similarity of Rendok and Rendok+j. For the image Ik+j, it is obvious that when j=±1, Ik+j is the neighbor of Ik, and there should be a larger ρj to weigh the similarity; when j increases or decreases, the Ik+j is far from Ik and the weight ρj should be decreased accordingly. Figure 3 illustrates the weight for the image Ik+j with varyingj, where a=1,b=2.

    Figure 3.  Weights corresponding to Eq (7), a = 1, b = 2.
    Figure 4.  Demonstration of segmentation of synthetic images, (a) an image sequence, (b) segmentation results of GVC snake without SSS, (c) segmentation results of GVC snake with SSS.

    The solution of Eq (5) is based on the shape similarity of the LV in the sequence, and the contours are changing during the snake deforming process. Thus, the deforming process of all of the contours should be carried on simultaneously, and all the images in the sequence are supposed to be segmented simultaneously. With the SSS and the strategies in Section II, we summary the segmentation algorithm in Algorithm 1.

    Algorithm 1. Segmentation of the cardiac MR image sequence
    Input Data: a sequence of the cardiac MR images, suppose K images
    Parameters: α, β, λ, η, n, h, a, b, m
    // First, segment the endocardium of the LV in each image
    1). Calculate the edge map of each image;
    2). Calculate the GVC external force for each image as in Eq (2);
    3). Initialize the snake contour in the kth image ck,endo(S),k=1,,K;
    4). The Number of converged snakes, iter=0;
    While iter<K do
        for k=1,,K do
            if ck,endo(S) is not converged, //For the kth image,
                Calculate the CSC in Eq (3);
                Calculate the SSS in Eq (5);
                    Deform the snake contour ck,endo(S);
                    if ck,endo(S) is converged, iter=iter+1; end
                end
        end
    end
    Return: the converged contour ck,endo(S),k=1,,K.
    // Then, segment the epicardium of the LV in each image
    1). Modify the edge map of each image as described in Section II.C;
    2). Calculate the GVC external force for each image as in Eq (2);
    3). Initialize the snake in the kth image ck,epi(S) using ck,endo(S), k=1,,K;
    4). The Number of converged snakes, iter=0;
    While iter<K do
            for k=1,,K do
                    if ck,epi(S) is not converged, //For the kth image,
                            Calculate the internal force in Eq (1);
                            Calculate the SSC in Eq (4);
                            Calculate the SSS in Eq (5);
                            Deform the snake contour ck,epi(S);
                            if ck,epi(S) is converged, iter=iter+1; end
                    end
            end
    end
    Return: the converged contour ck,epi(S), k=1,,K.

     | Show Table
    DownLoad: CSV

    In order to verify the performance of the proposed strategy, the parameters of our model are α=0.5,β=0.5,λ=0.4,η=0.6,n=2.6,h=1,a=1,b=2,m=1. In our experiments, two datasets are employed for evaluation: MICCAI 2009 [43] and the dataset published by Shuo Li (Shuo Dataset in short) [43,44,49]. Cardiac cine MRI data of MICCAI 2009 is provided by the MICCAI 2009 cardiac MR LV segmentation challenge organizers [43], and there are 2900 2D short-axis cine MR images of 145 subjects in the Shuo Dataset, and interested readers can find a detailed description in [49].

    Two metrics are employed to evaluate the segmentation performance, i.e., the Mean Absolute Distance (MAD) and the Hausdorff Distance (HD). Suppose the snake contour S={s1,s2,,sp} has p points, and the ground truth M={m1,m,,mq} has q points, the MAD and HD are respectively defined as

    MAD(S,M)=12(1ppi=1d(si,M)+1qqj=1d(mj,S)) (8)
    HD=max(maxsiSd(si,M),maxmjMd(mj,S)) (9)

    where d(si,M)=minsiS,mjMsimj is the distance from point si to the closet point on contour M and d(mj,S)=minsiS,mjMmjsi is that from point mj to the closet point on contour S.

    The MICCAI 2009 dataset is first employed for comparison, and the radial GVF method (RGVF) [16], level set method [23], the max-flow method (MFM) [25], and our previous work [11] are applied for comparison on this dataset, the detailed description of this dataset has been presented in [43]. The endocardia in Figure 5, especially in images #6–#10, are challenging for the snake model. However, the work in [11] regards the shape of the endocardium as a circle and employs the CSC as constraint, while the proposed strategy supposes that the shape of the endocardia in the consecutive image is similar and employs the SSS and CSC as constraints. The two methods both get pleasant results. But for epicardium segmentation, the artifacts and noise make some of the boundaries 'missing', and the snake moves to a false contour shown just with SSC in Figure 5(a), the results are a little far from pleasant. Figure 5(b) shows that the extra SSS effectively pushes the snake contour to the expected solution. Figure 6 shows the MAD errors corresponding to Figure 5. Overall, Integrating SSS to the GVC snake model can effectively reduce the effect of artifacts, and prevent the snake from weak edge leakage. One can learn from [11] the average MAD of RGVF, LSM, and MFM and our previous method for the endocardium segmentation are 9.06 pixels, 7.21 pixels, 4.79 pixels and 5.06 pixels, respectively, and those of the epicardium are 7.85 pixels, 5.70 pixels, 5.23 pixels and 5.18 pixels, respectively. Comparing them with the results shown in Figure 6, we can see the proposed strategy is able to achieve much better results. More representative segmentation results on MICCAI 2009 dataset are presented in Figure 7.

    Figure 5.  Experiments one subject in MICCAI 2009 dataset, (a) segmentation results using the strategies in [11], (b) segmentation results using the proposed strategies.
    Figure 6.  The MAD errors corresponding to Figure 5 using CSC, SSC and SSS.
    Figure 7.  Segmentation results on MICCAI 2009 dataset.

    To further demonstrate the effectiveness of SSS, the Shuo dataset is also employed for test [49]. Besides, in this experiment, we adopt the method of component edge map [48] after extracting the endocardium. First, with the endocardial center as a reference, we coin two edge maps by selectively opting the gradient in direction of x and y according to the grey level characteristics of images. In the left half part, the septum links the blood pools of RV and LV while the right ventricle is brighter. Thus, fxl=|min(Ix,0)| is utilized as an x-direction gradient to neglect the edge between the septum and the LV blood. For the right half, the lateral is surrounded by fat so we adopt fxr=max(Ix,0) to erase the edge between the fat and the lung as well as that between lateral and the LV pool. Similarly, we construct the y-direction gradient as follows: fyu=|min(Iy,0)|, fyd=max(Iy,0) where u,d represent the upper and lower half of the center respectively. The new edge map can be obtained by setting the fx,fy to zero around and within the endocardium. Then, this new edge map is employed to calculate a new external force by GVC. Figure 8 shows our results on a representative sequence of Shuo dataset, it can be seen that SSS combined with this new external force can get satisfactory results. The corresponding MAD indexes of endocardium and epicardium in this sequence were 1.15 pixels and 1.35 pixels respectively. We also calculated that the HD metrics of endocardium and epicardium in this sequence were 2.8210 mm and 3.5118 mm separately, while according to [43], the same indexes obtained from Graph cuts and Level set on some images in this dataset were 9.290 mm, 5.739 mm on endocardium and 4.200 mm, 5.484 mm on the epicardium. We can see that our method performs more excellently. Another set of experimental results of this dataset are shown in Figure 9.

    Figure 8.  Segmentation results of on a representative sequence of Shuo dataset.
    Figure 9.  Segmentation results on another sequence in Shuo dataset.

    In this work, we proposed a new strategy for left ventricle segmentation in cardiac cine MRI based on GVC snake. In the proposed strategy, all the images in a sequence are segmented simultaneously and the shape of the LV in each image is supposed to be similar to that of the LV in nearby images in the sequence. We coined the novel shape similarity constraint, and it is called sequential shape similarity (SSS). For the endocardium, the SSS is combined with the GVC snake with a circle constraint. With the SSC and SSS, the epicardia are segmented successively. Experimental results on the MICCAI 2009 dataset and Shuo dataset shown that the proposed segmentation approach performs well on conquering papillary muscles, weak edges, artifacts, and noise in the LV segmentation.

    This work was supported in part by the National Science Foundation Program of China (NSFC) (grant number: 61976241), and the International Science and technology cooperation plan project of Zhenjiang (grant number: GJ2021008).

    The authors declare there is no conflict of interest.



    [1] P. Croisille, D. Revel, MR imaging of the heart: functional imaging, Eur. Radiol., 10 (2000), 7–11. doi: 10.1007/s003300050003. doi: 10.1007/s003300050003
    [2] J. P. Earls, V. B. Ho, T. K. Foo, E. Castillo, S. D. Flamm, Cardiac MRI: recent progress and continued challenges, J. Magn. Reson. Imaging: JMRI, 16 (2002), 111–127. doi: 10.1002/jmri.10154. doi: 10.1002/jmri.10154
    [3] A. F. Frangi, W. J. Niessen, M. A. Viergever, Three-dimensional modeling for functional analysis of cardiac images, a review, IEEE Trans. Med. Imaging, 20 (2001), 2–5. doi: 10.1109/42.906421. doi: 10.1109/42.906421
    [4] D. Nguyen, K. Masterson, J. P. Vallée, Comparative evaluation of active contour model extensions for automated cardiac MR image segmentation by regional error assessment, Magn. Reson. Mater. Phys., Biol. Med., 20 (2007), 69–82. doi: 10.1007/s10334-007-0069-z. doi: 10.1007/s10334-007-0069-z
    [5] A. Fernández-Caballero, J. M. Vega-Riesco, Determining heart parameters through left ventricular automatic segmentation for heart disease diagnosis, Expert Syst. Appl., 36 (2009), 2234–2249. doi: 10.1016/j.eswa.2007.12.045. doi: 10.1016/j.eswa.2007.12.045
    [6] H. Hu, H. Liu, Z. Gao, L. Huang, Hybrid segmentation of left ventricle in cardiac MRI using gaussian-mixture model and region restricted dynamic programming, Magn. Reson. Imaging, 31 (2013), 575–584. doi: 10.1016/j.mri.2012.10.004. doi: 10.1016/j.mri.2012.10.004
    [7] M. Lorenzo-Valdés, G. I. Sanchez-Ortiz, A. G. Elkington, R. H. Mohiaddin, D. Rueckert, Segmentation of 4D cardiac MR images using a probabilistic atlas and the EM algorithm, Med. Image Anal., 8 (2004), 255–265. doi: 10.1016/j.media.2004.09.005. doi: 10.1016/j.media.2004.09.005
    [8] W. Bai, W. Shi, C. Ledig, D. Rueckert, Multi-atlas segmentation with augmented features for cardiac MR images, Med. Image Anal., 19 (2015), 98–109. doi: 10.1016/j.media.2014.09.005. doi: 10.1016/j.media.2014.09.005
    [9] J. Cousty, L. Najman, M. Couprie, S. Clément-Guinaudeau, T. Goissen, J. Garot, Segmentation of 4D cardiac MRI: Automated method based on spatio-temporal watershed cuts, Image Vision Comput., 28 (2010), 1229–1243. doi: 10.1016/j.imavis.2010.01.001. doi: 10.1016/j.imavis.2010.01.001
    [10] M. G. Uzunbaş, S. Zhang, K. M. Pohl, D. Metaxas, L. Axel, Segmentation of myocardium using deformable regions and graph cuts, in 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI), (2012), 254–257. doi: 10.1109/ISBI.2012.6235532.
    [11] Y. Wu, Y. Wang, Y. Jia, Segmentation of the left ventricle in cardiac cine MRI using a shape-constrained snake model, Comput. Vision Image Understanding, 117 (2013), 990–1003. doi: 10.1016/j.cviu.2012.12.008. doi: 10.1016/j.cviu.2012.12.008
    [12] Z. Zhang, C. Duan, T. Lin, S. Zhou, Y. Wang, X. Gao, GVFOM: a novel external force for active contour based image segmentation, Inf. Sci., 506 (2020), 1–18. doi: 10.1016/j.ins.2019.08.003. doi: 10.1016/j.ins.2019.08.003
    [13] M. F. Santarelli, V. Positano, C. Michelassi, M. Lombardi, L. Landini, Automated cardiac MR image segmentation: theory and measurement evaluation, Med. Eng. Phys., 25 (2003), 149–159. doi: 10.1016/S1350-4533(02)00144-3. doi: 10.1016/S1350-4533(02)00144-3
    [14] S. Ranganath, Contour extraction from cardiac MRI studies using snakes, IEEE Trans. Med. Imaging, 14 (1995), 328–338. doi: 10.1109/42.387714. doi: 10.1109/42.387714
    [15] T. McInerney, D. Terzopoulos, A dynamic finite element surface model for segmentation and tracking in multidimensional medical images with application to cardiac 4D image analysis, Comput. Med. Imaging Graphics, 19 (1995), 69–83, 1995. doi: 10.1016/0895-6111(94)00040-9. doi: 10.1016/0895-6111(94)00040-9
    [16] J. Liang, G. Ding, Y. Wu, Segmentation of the left ventricle from cardiac MR images based on radial GVF snake, in 2008 International Conference on BioMedical Engineering and Informatics, IEEE, 2 (2008), 238–242. doi: 10.1109/BMEI.2008.188.
    [17] C. Feng, S. Zhang, D. Zhao, C. Li, Simultaneous extraction of endocardial and epicardial contours of the left ventricle by distance regularized level sets, Med. Phys., 43 (2016), 2741–2755. doi: 10.1118/1.4947126. doi: 10.1118/1.4947126
    [18] Y. Liu, G. Captur, J. C. Moon, S. Guo, X. Yang, S. Zhang, et al., Distance regularized two level sets for segmentation of left and right ventricles from cine-MRI, Magn. Reson. Imaging, 34 (2016), 699–706. doi: 10.1016/j.mri.2015.12.027. doi: 10.1016/j.mri.2015.12.027
    [19] J. Lu, C. Feng, J. Yang, W. Li, D. Zhao, C. Wan, Segmentation of the cardiac ventricle using two layer level sets with prior shape constraint. Biomedical Signal Processing and Control, 68 (2021), 102671. doi: 10.1016/j.bspc.2021.102671. doi: 10.1016/j.bspc.2021.102671
    [20] C. Feng, C. Li, D. Zhao, C. Davatzikos, H. Litt, Segmentation of the left ventricle using distance regularized two-layer level set approach, in International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, Berlin, Heidelberg, (2013), 477–484. doi: 10.1007/978-3-642-40811-3_60.
    [21] J. Lu, C. Feng, W. Li, D. Zhao, ROI localization and initialization method for left ventricle segmentation, in Proceedings of the Third International Symposium on Image Computing and Digital Medicine, (2019), 12–16.
    [22] J. Lu, C. Feng, D. Zhao, Segmentation of the cardiac left ventricle from cine magnetic resonance images using local inhomogeneous intensity clustering with prior shape constraint, J. Med. Imaging Health Inf., 9 (2019), 70–77. doi: 10.1166/jmihi.2019.2542. doi: 10.1166/jmihi.2019.2542
    [23] I. B. Ayed, S. Li, I. Ross, Embedding overlap priors in variational left ventricle tracking, IEEE Trans. Med. Imaging, 28 (2009), 1902–1913. doi: 10.1109/TMI.2009.2022087. doi: 10.1109/TMI.2009.2022087
    [24] I. B. Ayed, H. Chen, K. Punithakumar, I. Ross, S. Li, Max-flow segmentation of the left ventricle by recovering subject-specific distributions via a bound of the Bhattacharyya measure, Med. Image Anal., 16 (2012), 87–100. doi: 10.1016/j.media.2011.05.009. doi: 10.1016/j.media.2011.05.009
    [25] N. Paragios, A level set approach for shape-driven segmentation and tracking of the left ventricle, IEEE Trans. Med. Imaging, 22 (2003), 773–776. doi: 10.1109/TMI.2003.814785. doi: 10.1109/TMI.2003.814785
    [26] M. Lynch, O. Ghita, P. F. Whelan, Segmentation of the left ventricle of the heart in 3-D+t MRI data using an optimized nonrigid temporal model, IEEE Trans. Med. Imaging, 27 (2008), 195–203. doi: 10.1109/TMI.2007.904681. doi: 10.1109/TMI.2007.904681
    [27] M. Lynch, O. Ghita, P. F. Whelan, Left-ventricle myocardium segmentation using a coupled level-set with a priori knowledge, Comput. Med. Imaging Graphics, 30 (2006), 255–262. doi: 10.1016/j.compmedimag.2006.03.009. doi: 10.1016/j.compmedimag.2006.03.009
    [28] C. Pluempitiwiriyawej, J. M. F. Moura, Y. L. Wu, C. Ho, STACS: new active contour scheme for cardiac MR image segmentation, IEEE Trans. Med. Imaging, 24 (2005), 593–603. doi: 10.1109/TMI.2005.843740. doi: 10.1109/TMI.2005.843740
    [29] T. Chen, J. Babb, P. Kellman, L. Axel, D. Kim, Semiautomated segmentation of myocardial contours for fast strain analysis in cine displacement-encoded MRI, IEEE Trans. Med. Imaging, 27 (2008), 1084–1094. doi: 10.1109/TMI.2008.918327. doi: 10.1109/TMI.2008.918327
    [30] J. Woo, P. J. Slomka, C. C. J. Kuo, B. W. Hong, Multiphase segmentation using an implicit dual shape prior: Application to detection of left ventricle in cardiac MRI, Comput. Vision Image Understanding, 117 (2013), 1084–1094. doi: 10.1016/j.cviu.2012.11.012. doi: 10.1016/j.cviu.2012.11.012
    [31] H. Lee, N. C. F. Codella, M. D. Cham, J. W. Weinsaft, Y. Wang, Automatic left ventricle segmentation using iterative thresholding and an active contour model with adaptation on short-axis cardiac MRI, IEEE Trans. Biomed. Eng., 57 (2010), 905–913. doi: 10.1109/TBME.2009.2014545. doi: 10.1109/TBME.2009.2014545
    [32] R. Beichel, H. Bischof, F. Leberl, M. Sonka, Robust active appearance models and their application to medical image analysis, IEEE Trans. Med. Imaging, 24 (2005), 1151–1169. doi: 10.1109/TMI.2005.853237. doi: 10.1109/TMI.2005.853237
    [33] J. Montagnat, H. Delingette, 4D deformable models with temporal constraints: application to 4D cardiac image segmentation, Med. Image Anal., 9 (2005), 87–100. doi: 10.1016/j.media.2004.06.025. doi: 10.1016/j.media.2004.06.025
    [34] H. C. van Assen, M. G. Danilouchkine, A. F. Frangi, S. Ordás, J. J. M. Westenberg, J. H. C. Reiber, et al., SPASM: A 3D-ASM for segmentation of sparse and arbitrarily oriented cardiac MRI data, Med. Image Anal., 10 (2006), 286–303. doi: 10.1016/j.media.2005.12.001. doi: 10.1016/j.media.2005.12.001
    [35] A. F. Frangi, D. Rueckert, J. A. Schnabel, W. J. Niessen, Automatic construction of multiple-object three-dimensional statistical shape models: application to cardiac modeling, IEEE Trans. Med. Imaging, 21 (2002), 1151–1166. doi: 10.1109/TMI.2002.804426. doi: 10.1109/TMI.2002.804426
    [36] J. Lötjönen, S. Kivistö, J. Koikkalainen, D. Smutek, K. Lauerma, Statistical shape model of atria, ventricles and epicardium from short- and long-axis MR images, Med. Image Anal., 8 (2004), 371–386. doi: 10.1016/j.media.2004.06.013. doi: 10.1016/j.media.2004.06.013
    [37] M. P. Jolly, Automatic segmentation of the left ventricle in cardiac MR and CT images, Int. J. Comput. Vision, 70 (2006), 151–163. doi: 10.1007/s11263-006-7936-3. doi: 10.1007/s11263-006-7936-3
    [38] X. Artaechevarria, A. Munoz-Barrutia, C. Ortiz-de-Solorzano, Combination strategies in multi-atlas image segmentation: application to brain MR data, IEEE Trans. Med. Imaging, 28 (2009), 1266–1277. doi: 10.1109/TMI.2009.2014372. doi: 10.1109/TMI.2009.2014372
    [39] C. Petitjean, J. N. Dacher, A review of segmentation methods in short axis cardiac MR images, Med. Image Anal., 15 (2011), 169–184. doi: 10.1016/j.media.2010.12.004. doi: 10.1016/j.media.2010.12.004
    [40] P. Peng, K. Lekadir, A. Gooya, L. Shao, S. E. Petersen, A. F. Frangi, A review of heart chamber segmentation for structural and functional analysis using cardiac magnetic resonance imaging, Magn. Reson. Mater. Phys., Biol. Med., 29 (2016), 155–195. doi: 10.1007/s10334-015-0521-4. doi: 10.1007/s10334-015-0521-4
    [41] M. Kass, A. Witkin, D. Terzopoulos, Snakes: Active contour models, Int. J. Comput. Vision, 1 (1988), 321–331. doi: 10.1007/BF00133570. doi: 10.1007/BF00133570
    [42] C. Xu, J. L. Prince, Snakes, shapes, and gradient vector flow, IEEE Trans. Image Process., 7 (1998), 359–369. doi: 10.1109/83.661186. doi: 10.1109/83.661186
    [43] H. Zhang, W. Zhang, W. Shen, N. Li, Y. Chen, S. Li, et al., Automatic segmentation of the cardiac MR images based on nested fully convolutional dense network with dilated convolution, Biomed. Signal Process. Control, 68 (2021), 102684. doi: 10.1016/j.bspc.2021.102684. doi: 10.1016/j.bspc.2021.102684
    [44] W. Wang, Y. Wang, Y. Wu, T. Lin, S. Li, B. Chen, Quantification of full left ventricular metrics via deep regression learning with contour-guidance, IEEE Access, 7 (2019), 47918–47928. doi: 10.1109/ACCESS.2019.2907564. doi: 10.1109/ACCESS.2019.2907564
    [45] W. Shen, W. Wu, H. Zhang, Z. Sun, J. Ma, S. Guo, et al., Automatic segmentation of the femur and tibia bones from X-ray images based on pure dilated residual U-Net, Inverse Probl. Imaging, 15 (2021). doi: 10.3934/ipi.2020057. doi: 10.3934/ipi.2020057
    [46] T. A. Ngo, Z. Lu, G. Carneiro, Combining deep learning and level set for the automated segmentation of the left ventricle of the heart from cardiac cine magnetic resonance, Med. Image Anal., 35 (2017), 159–171. doi: 10.1016/j.media.2016.05.009. doi: 10.1016/j.media.2016.05.009
    [47] M. R. Avendi, A. Kheradvar, H. Jafarkhani, A combined deep-learning and deformable-model approach to fully automatic segmentation of the left ventricle in cardiac MRI, Med. Image Anal., 30 (2016), 108–119. doi: 10.1016/j.media.2016.01.005. doi: 10.1016/j.media.2016.01.005
    [48] Y. Wang, Y. Jia, External force for active contours: gradient vector convolution, in PRICAI 2008: Trends in Artificial Intelligence, Berlin, Heidelberg, (2008), 466–472. doi: 10.1007/978-3-540-89197-0_43.
    [49] W. Xue, G. Brahm, S. Pandey, S. Leung, S. Li, Full left ventricle quantification via deep multitask relationships learning, Med. Image Anal., 43 (2018), 54–65. doi: 10.1016/j.media.2017.09.005. doi: 10.1016/j.media.2017.09.005
  • This article has been cited by:

    1. Bingjie Li, Tiejun Yang, Xiang Zhao, NVTrans‐UNet: Neighborhood vision transformer based U‐Net for multi‐modal cardiac MR image segmentation, 2023, 1526-9914, 10.1002/acm2.13908
    2. Yudong Zhang, Juan Manuel Gorriz, Deepak Ranjan Nayak, Optimization Algorithms and Machine Learning Techniques in Medical Image Analysis, 2023, 20, 1551-0018, 5917, 10.3934/mbe.2023255
    3. G. Sucharitha, Veeramalai sankardass, R. Rani, Nagaraj Bhat, A. Rajaram, RETRACTED: Deep learning aided prostate cancer detection for early diagnosis & treatment using MR with TRUS images, 2024, 46, 10641246, 3395, 10.3233/JIFS-235744
    4. Mehreen Irshad, Mussarat Yasmin, Muhammad Imran Sharif, Muhammad Rashid, Muhammad Irfan Sharif, Seifedine Kadry, A Novel Light U-Net Model for Left Ventricle Segmentation Using MRI, 2023, 11, 2227-7390, 3245, 10.3390/math11143245
    5. Alan Cervantes-Guzmán, Kyle McPherson, Jimena Olveres, Carlos Francisco Moreno-García, Fabián Torres Robles, Eyad Elyan, Boris Escalante-Ramírez, Anas Bilal, Robust cardiac segmentation corrected with heuristics, 2023, 18, 1932-6203, e0293560, 10.1371/journal.pone.0293560
    6. Sun Zheng, Zhu Aoying, Hou Yingsa, Sun Meichen, A physics-based iterative learning framework for quantitative parametric imaging with application to photoacoustic imaging, 2025, 142, 09521976, 109920, 10.1016/j.engappai.2024.109920
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2781) PDF downloads(89) Cited by(6)

Figures and Tables

Figures(9)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog