Research article Special Issues

Residual based attention-Unet combing DAC and RMP modules for automatic liver tumor segmentation in CT


  • Received: 12 January 2022 Revised: 20 February 2022 Accepted: 24 February 2022 Published: 09 March 2022
  • Purpose: Due to the complex distribution of liver tumors in the abdomen, the accuracy of liver tumor segmentation cannot meet the needs of clinical assistance yet. This paper aims to propose a new end-to-end network to improve the segmentation accuracy of liver tumors from CT. Method: We proposed a hybrid network, leveraging the residual block, the context encoder (CE), and the Attention-Unet, called ResCEAttUnet. The CE comprises a dense atrous convolution (DAC) module and a residual multi-kernel pooling (RMP) module. The DAC module ensures the network derives high-level semantic information and minimizes detailed information loss. The RMP module improves the ability of the network to extract multi-scale features. Moreover, a hybrid loss function based on cross-entropy and Tversky loss function is employed to distribute the weights of the two-loss parts through training iterations. Results: We evaluated the proposed method in LiTS17 and 3DIRCADb databases. It significantly improved the segmentation accuracy compared to state-of-the-art methods. Conclusions: Experimental results demonstrate the satisfying effects of the proposed method through both quantitative and qualitative analyses, thus proving a promising tool in liver tumor segmentation.

    Citation: Rongrong Bi, Chunlei Ji, Zhipeng Yang, Meixia Qiao, Peiqing Lv, Haiying Wang. Residual based attention-Unet combing DAC and RMP modules for automatic liver tumor segmentation in CT[J]. Mathematical Biosciences and Engineering, 2022, 19(5): 4703-4718. doi: 10.3934/mbe.2022219

    Related Papers:

    [1] Nalin Fonseka, Jerome Goddard Ⅱ, Alketa Henderson, Dustin Nichols, Ratnasingham Shivaji . Modeling effects of matrix heterogeneity on population persistence at the patch-level. Mathematical Biosciences and Engineering, 2022, 19(12): 13675-13709. doi: 10.3934/mbe.2022638
    [2] Robert Stephen Cantrell, Chris Cosner, William F. Fagan . The implications of model formulation when transitioning from spatial to landscape ecology. Mathematical Biosciences and Engineering, 2012, 9(1): 27-60. doi: 10.3934/mbe.2012.9.27
    [3] Nazanin Zaker, Christina A. Cobbold, Frithjof Lutscher . The effect of landscape fragmentation on Turing-pattern formation. Mathematical Biosciences and Engineering, 2022, 19(3): 2506-2537. doi: 10.3934/mbe.2022116
    [4] Zhilan Feng, Robert Swihart, Yingfei Yi, Huaiping Zhu . Coexistence in a metapopulation model with explicit local dynamics. Mathematical Biosciences and Engineering, 2004, 1(1): 131-145. doi: 10.3934/mbe.2004.1.131
    [5] James T. Cronin, Jerome Goddard II, Amila Muthunayake, Ratnasingham Shivaji . Modeling the effects of trait-mediated dispersal on coexistence of mutualists. Mathematical Biosciences and Engineering, 2020, 17(6): 7838-7861. doi: 10.3934/mbe.2020399
    [6] Nariyuki Nakagiri, Hiroki Yokoi, Yukio Sakisaka, Kei-ichi Tainaka . Population persistence under two conservation measures: Paradox of habitat protection in a patchy environment. Mathematical Biosciences and Engineering, 2022, 19(9): 9244-9257. doi: 10.3934/mbe.2022429
    [7] Kunquan Lan, Wei Lin . Population models with quasi-constant-yield harvest rates. Mathematical Biosciences and Engineering, 2017, 14(2): 467-490. doi: 10.3934/mbe.2017029
    [8] Minjuan Gao, Lijuan Chen, Fengde Chen . Dynamical analysis of a discrete two-patch model with the Allee effect and nonlinear dispersal. Mathematical Biosciences and Engineering, 2024, 21(4): 5499-5520. doi: 10.3934/mbe.2024242
    [9] Robert Stephen Cantrell, Chris Cosner, William F. Fagan . Edge-linked dynamics and the scale-dependence of competitive. Mathematical Biosciences and Engineering, 2005, 2(4): 833-868. doi: 10.3934/mbe.2005.2.833
    [10] Kuang-Hui Lin, Yuan Lou, Chih-Wen Shih, Tze-Hung Tsai . Global dynamics for two-species competition in patchy environment. Mathematical Biosciences and Engineering, 2014, 11(4): 947-970. doi: 10.3934/mbe.2014.11.947
  • Purpose: Due to the complex distribution of liver tumors in the abdomen, the accuracy of liver tumor segmentation cannot meet the needs of clinical assistance yet. This paper aims to propose a new end-to-end network to improve the segmentation accuracy of liver tumors from CT. Method: We proposed a hybrid network, leveraging the residual block, the context encoder (CE), and the Attention-Unet, called ResCEAttUnet. The CE comprises a dense atrous convolution (DAC) module and a residual multi-kernel pooling (RMP) module. The DAC module ensures the network derives high-level semantic information and minimizes detailed information loss. The RMP module improves the ability of the network to extract multi-scale features. Moreover, a hybrid loss function based on cross-entropy and Tversky loss function is employed to distribute the weights of the two-loss parts through training iterations. Results: We evaluated the proposed method in LiTS17 and 3DIRCADb databases. It significantly improved the segmentation accuracy compared to state-of-the-art methods. Conclusions: Experimental results demonstrate the satisfying effects of the proposed method through both quantitative and qualitative analyses, thus proving a promising tool in liver tumor segmentation.



    A molecular graph is a graph such that its vertices correspond to the atoms and the edges to the bonds. Chemical Graph Theory is a branch of Mathematical Chemistry, which has an important effect on the development of the Chemical Sciences. A topological index is a numerical parameter mathematically derived from the graph structure. Several such topological indices have been considered in Theoretical Chemistry and have found some applications, especially in QSPR/QSAR study see [1,2].

    Let G be a finite, simple, connected graph with vertex set V(G) and edge set E(G). The degree dG(u) of a vertex u is the number of vertices adjacent to u. Let Δ(G) and δ(G) denote the maximum and minimum degree, respectively, among the vertices of G. The Revan vertex degree of a vertex uV(G) is defined as rG(u)=Δ(G)+δ(G)dG(u). The edge connecting the vertices u,vV(G) will be denoted by uv. We refer to [3] for undefined notations and terminologies.

    The first and second Revan indices of a graph G were introduced by Kulli in [4], and they are defined as

    R1(G)=uvE(G)(rG(u)+rG(v)),R2(G)=uvE(G)rG(u)rG(v).

    The F-Revan index of a graph is defined in [5] as

    FR(G)=uvE(G)(rG(u)2+rG(v)2).

    The Sombor index was introduced by Gutman in [6], defined as

    SO(G)=uvE(G)dG(u)2+dG(v)2.

    Recently, a large number of studies on Sombor indices have been published, see e.g., [7,8,9,10,11,12,13,14,15].

    The Revan Sombor index was proposed by Kulli and Gutman in [16] as

    RSO(G)=uvE(G)rG(u)2+rG(v)2.

    The first and second Revan (a,b)-KA indices of a graph G are defined (generalizing the Revan Sombor index), respectively, as

    RKA1a,b(G)=uvE(G)(rG(u)a+rG(v)a)b,RKA21,c(G)=uvE(G)(rG(u)rG(v))c.

    Recently, some (a,b)-KA indices have been used to analyze the properties of chemical compounds and nanostructures. In particular, they have been applied to the study of benzenoid systems and phenylenes, see for example [17,18,19,20].

    Therefore, motivated both by their practical applications as well as by their theoretical properties, mainly because of the close relationship with the well-known Sombor index [15], in this paper we perform analytical and statistical studies on the recently introduced Revan Sombor indices.

    Theorem 1. If G is a graph with maximum degree Δ and minimum degree δ, and cR{0}, then

    2c/2ΔcRKA21,c(G)RKA12,c/2(G)2c/2δcRKA21,c(G) if c>0,2c/2δcRKA21,c(G)RKA12,c/2(G)2c/2ΔcRKA21,c(G) if c<0,

    and the equality in each bound is attained for each c if and only if G is a regular graph.

    Proof. If c>0, then

    RKA12,c/2(G)=uvE(G)(rG(u)2+rG(v)2)c/2=uvE(G)(rG(u)rG(v))c(1rG(u)2+1rG(v)2)c/2uvE(G)(rG(u)rG(v))c(1Δ2+1Δ2)c/2=2c/2ΔcRKA21,c(G),

    and we obtain the converse inequality if c<0. The equality in each bound is attained if and only if rG(u)=rG(v)=Δ, i.e., dG(u)=dG(v)=δ for every uvE(G), that is, G is a regular graph.

    Also, if c>0, then

    RKA12,c/2(G)=uvE(G)(rG(u)rG(v))c(1rG(u)2+1rG(v)2)c/2uvE(G)(rG(u)rG(v))c(1δ2+1δ2)c/2=2c/2δcRKA21,c(G),

    and we obtain the converse inequality if c<0. The equality in each bound is attained if and only if rG(u)=rG(v)=δ, i.e., dG(u)=dG(v)=Δ for every uvE(G), that is, G is a regular graph.

    Notice that the limit as c approaches zero, the bounds, both upper and lower, are all equal to m (number of edges of G).

    If we choose c=1 in Theorem 1, then we obtain the following result.

    Corollary 2. If G is a graph, then

    2ΔR2(G)RSO(G)2δR2(G)

    and the equality in each bound is attained if and only if G is a regular graph.

    Theorem 3. If G is a graph with m edges and maximum degree Δ, and c>0, then

    RKA12,c/2(G)mΔcRKA11,c(G),

    and the equality in the bound is attained for each c if and only if G is a regular graph.

    Proof. Cauchy-Schwarz inequality gives

    RKA12,c/2(G)=uvE(G)(rG(u)2+rG(v)2)c/2uvE(G)Δc/2(rG(u)+rG(v))c/2(uvE(G)Δc)(uvE(G)(rG(u)+rG(v))c)=mΔcRKA11,c(G).

    If the equality in the bound is attained for some c>0, then rG(u)=rG(v)=Δ for every uvE(G) and so, G is a regular graph.

    If G is regular, then

    mΔcRKA11,c(G)=mΔcm2cΔc=m2c/2Δc=RKA12,c/2(G).

    If we choose c=1 in Theorem 3, then we obtain the following result.

    Corollary 4. If G is a graph with m edges and maximum degree Δ, then

    RSO(G)mΔR1(G)

    and the equality in the bound is attained if and only if G is a regular graph.

    We need the following well known inequalities for x1,x2>0 [21]:

    xα1+xα2<(x1+x2)α2α1(xα1+xα2)if α>1,2α1(xα1+xα2)(x1+x2)α<xα1+xα2if 0<α<1,(x1+x2)α2α1(xα1+xα2)if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if x1=x2.

    These inequalities allow to obtain the following result.

    Theorem 5. If G is a graph and a,b,cR{0}, then

    RKA1ab/c,c(G)<RKA1a,b(G)2bcRKA1ab/c,c(G) if b>c,bc>0,2bcRKA1ab/c,c(G)RKA1a,b(G)<RKA1ab/c,c(G) if b<c,bc>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b<0,c>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b>0,c<0,

    and the equality in each non-strict inequality is attained for each a,b,c if and only if each connected component of G is a regular graph.

    Proof. If α=b/c, x1=rG(u)a and x2=rG(v)a, then the previous inequalities give

    rG(u)ab/c+rG(v)ab/c<(rG(u)a+rG(v)a)b/c2b/c1(rG(u)ab/c+rG(v)ab/c)if b/c>1,2b/c1(rG(u)ab/c+rG(v)ab/c)(rG(u)a+rG(v)a)b/c<rG(u)ab/c+rG(v)ab/cif 0<b/c<1,(rG(u)a+rG(v)a)b/c2b/c1(rG(u)ab/c+rG(v)ab/c)if b/c<0,

    and the equality in each non-strict inequality is attained if and only if rG(u)=rG(v), i.e., dG(u)=dG(v).

    Hence, we obtain

    (rG(u)ab/c+rG(v)ab/c)c<(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif b/c>1,c>0,2bc(rG(u)ab/c+rG(v)ab/c)c(rG(u)a+rG(v)a)b<(rG(u)ab/c+rG(v)ab/c)cif b/c>1,c<0,2bc(rG(u)ab/c+rG(v)ab/c)c(rG(u)a+rG(v)a)b<(rG(u)ab/c+rG(v)ab/c)cif 0<b/c<1,c>0,(rG(u)ab/c+rG(v)ab/c)c<(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif 0<b/c<1,c<0,(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif b<0,c>0,(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif b>0,c<0,

    and the equality in each non-strict inequality is attained if and only if dG(u)=dG(v).

    If we sum on uvE(G) these inequalities, then we obtain the desired inequalities. Furthermore, the equality in each non-strict inequality is attained for each a,b,c if and only if dG(u)=dG(v) for every uvE(G); and this happens if and only if each connected component of G is a regular graph.

    Remark 6. The excluded case b=c in Theorem 5 is not interesting, since RKA1ab/c,c(G)=RKA1a,b(G) if b=c.

    If we take a=2, b=1/2 and c=1/α in Theorem 5, then we obtain the following result for the Revan Sombor index.

    Corollary 7. If G is a graph and αR{0}, then

    RKA1α,1/α(G)<RSO(G)21/21/αRKA1α,1/α(G) if α>2,21/21/αRKA1α,1/α(G)RSO(G)<RKA1α,1/α(G) if 0<α<2,21/21/αRKA1α,1/α(G)RSO(G) if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if each connected component of G is a regular graph.

    The above bounds (in Theorem 5 and Corollary 7) are achieved for several families of graphs, e.g., the complete graph Kn, the complete bipartite graph Km,n, the cycle graph Cn, etc.

    Recall that one of the most studied topological indices is the first Zagreb index, defined by

    M1(G)=uvE(G)(dG(u)+dG(v))=uV(G)dG(u)2.

    Our next result relates M1 and RKA1α,1/α.

    Theorem 8. If G is a graph with m edges, maximum degree Δ and minimum degree δ, and αR{0}, then

    2(Δ+δ)m211/αRKA1α,1/α(G)M1(G)<2(Δ+δ)mRKA1α,1/α(G) if α>1,2(Δ+δ)mRKA1α,1/α(G)<M1(G)2(Δ+δ)m211/αRKA1α,1/α(G) if 0<α<1,M1(G)2(Δ+δ)m211/αRKA1α,1/α(G) if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if each connected component of G is a regular graph.

    Proof. If we take a=b=1 and c=1/α in Theorem 5, then we obtain the following result for RKA11,1(G):

    RKA1α,1/α(G)<RKA11,1(G)211/αRKA1α,1/α(G)if α>1,211/αRKA1α,1/α(G)RKA11,1(G)<RKA1α,1/α(G)if 0<α<1,211/αRKA1α,1/α(G)RKA11,1(G)if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if each connected component of G is a regular graph.

    These inequalities and the identity

    RKA11,1(G)=uvE(G)(rG(u)+rG(v))=2(Δ+δ)mM1(G)

    finish the proof.

    The first and second (a,b)-KA indices of a graph G are defined in [17], respectively, as

    KA1a,b(G)=uvE(G)(dG(u)a+dG(v)a)b,KA2a,b(G)=uvE(G)(dG(u)adG(v)a)b.

    The following result relates the KA1a,b and RKA1a,b indices.

    Theorem 9. If G is a graph with maximum degree Δ and minimum degree δ, and a,bR{0}, then

    δabΔabKA1a,b(G)RKA1a,b(G)ΔabδabKA1a,b(G) if ab>0,ΔabδabKA1a,b(G)RKA1a,b(G)δabΔabKA1a,b(G) if ab<0,

    and the equality in each inequality is attained for each a,b if and only if G is a regular graph.

    Proof. Let us define the function f:[δ,Δ]×[δ,Δ](0,) by

    f(x,y)=(Δ+δx)a+(Δ+δy)axa+ya.

    Note that if a>0, the numerator of f is strictly decreasing on [δ,Δ] on each variable and its denominator is strictly increasing on [δ,Δ] on each variable. Hence, f is a strictly decreasing function on [δ,Δ] on each variable, and so,

    δaΔa=f(Δ,Δ)f(x,y)f(δ,δ)=Δaδa.

    If x=dG(u) and y=dG(v), then

    δaΔa(dG(u)a+dG(v)a)rG(u)a+rG(v)aΔaδa(dG(u)a+dG(v)a).

    Notice that if b>0, then the above inequality is preserved; thus

    δabΔab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bΔabδab(dG(u)a+dG(v)a)b.

    In contrast, if b<0, then the inequality above is not preserved; thus

    Δabδab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bδabΔab(dG(u)a+dG(v)a)b.

    If a<0, then f is a strictly increasing function on [δ,Δ] on each variable, and so,

    Δaδa(dG(u)a+dG(v)a)rG(u)a+rG(v)aδaΔa(dG(u)a+dG(v)a).

    Thus,

    δabΔab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bΔabδab(dG(u)a+dG(v)a)b

    if b<0, and

    Δabδab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bδabΔab(dG(u)a+dG(v)a)b

    if b>0.

    If we sum on uvE(G) these inequalities, then we obtain the desired inequalities. Furthermore, since f is a strictly monotone function, the equality in each inequality is attained for each a,b if and only if dG(u)=dG(v)=δ for every uvE(G) or dG(u)=dG(v)=Δ for every uvE(G); and this happens if and only if G is a regular graph.

    The Revan Sombor and the Sombor indices have the same value for every regular graph. The following result, which is a consequence of Theorem 9 with a=2 and b=1/2, relates these two indices for every graph.

    Corollary 10. If G is a graph with maximum degree Δ and minimum degree δ, then

    δΔSO(G)RSO(G)ΔδSO(G)

    and the equality in each inequality is attained if and only if G is a regular graph.

    Another remarkable topological index is the harmonic index, defined in [22] as

    H(G)=uvE(G)2du+dv.

    This index has attracted a great interest in the lasts years (see, e.g., [23,24,25,26,27]).

    The next result relates the RKA1a,1/a and the harmonic indices.

    Theorem 11. If G is a graph with maximum degree Δ and minimum degree δ, and aR, then

    21/aΔδH(G)RKA1a,1/a(G)<12(Δ+δ)2H(G) if a>1,2ΔδH(G)<RKA1a,1/a(G)22+1/a(Δ+δ)2H(G) if 0<a<1,RKA1a,1/a(G)22+1/a(Δ+δ)2H(G) if a<0.

    The equality in the first inequality is attained for each a if and only if each connected component of G is Δ-regular or δ-regular. The equality in the two last inequalities is attained for each a if G is a regular graph.

    Proof. Recall that

    (xa1+xa2)1/a<x1+x2211/a(xa1+xa2)1/aif a>1,211/a(xa1+xa2)1/ax1+x2<(xa1+xa2)1/aif 0<a<1,211/a(xa1+xa2)1/ax1+x2if a<0,

    for x1,x2>0, and the equality in each non-strict inequality is attained for each a if and only if x1=x2.

    Let us define the function g:[δ,Δ]×[δ,Δ](0,) by

    g(x,y)=((Δ+δx)a+(Δ+δy)a)1/a(x+y).

    If a>1, then

    g(x,y)<((Δ+δx)+(Δ+δy))(x+y)=(2(Δ+δ)(x+y))(x+y)(2(Δ+δ)(Δ+δ))(Δ+δ)=(Δ+δ)2,

    since x+y[2δ,2Δ]. Thus,

    (rG(u)a+rG(v)a)1/a<(Δ+δ)21dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)<12(Δ+δ)2H(G).

    Also,

    g(x,y)21+1/a((Δ+δx)+(Δ+δy))(x+y)=21+1/a(2(Δ+δ)(x+y))(x+y)21+1/a4Δδ

    since x+y[2δ,2Δ]. The equality in the previous bound is attained if and only if we have either x=y=δ or x=y=Δ. Therefore,

    (rG(u)a+rG(v)a)1/a21/aΔδ2dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)21/aΔδH(G).

    The equality in this bound is attained for each a>1 if and only if for each edge uvE(G) we have either dG(u)=dG(v)=Δ or dG(u)=dG(v)=δ, that is, each connected component of G is Δ-regular or δ-regular.

    Assume now a<1, then

    g(x,y)21+1/a((Δ+δx)+(Δ+δy))(x+y)=21+1/a(2(Δ+δ)(x+y))(x+y)21+1/a(2(Δ+δ)(Δ+δ))(Δ+δ)=21+1/a(Δ+δ)2,

    since x+y[2δ,2Δ]. Thus,

    (rG(u)a+rG(v)a)1/a22+1/a(Δ+δ)22dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)22+1/a(Δ+δ)2H(G).

    If G is a regular graph with m edges, then

    RKA1a,1/a(G)=21/aΔm=22+1/a(2Δ)2mΔ=22+1/a(Δ+δ)2H(G).

    If 0<a<1, then we have

    g(x,y)>((Δ+δx)+(Δ+δy))(x+y)=(2(Δ+δ)(x+y))(x+y)4Δδ

    since x+y[2δ,2Δ]. Therefore,

    (rG(u)a+rG(v)a)1/a>2Δδ2dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)>2ΔδH(G).

    If we choose a=2 in Theorem 11, then we obtain the following result for the Revan Sombor index.

    Corollary 12. If G is a graph with maximum degree Δ and minimum degree δ, then

    2ΔδH(G)RSO(G)<12(Δ+δ)2H(G)

    and the equality in the first inequality is attained if and only if each connected component of G is Δ-regular or δ-regular.

    Recently, within a statistical approach to degree-based topological indices (TIs) on random graphs (see e.g., [28,29,30]), it has been shown that the average values of TIs like the average Randic index, R(G), the average harmonic index, H(G), and the average of some Sombor indices are highly correlated with the average Shannon entropy of the eigenvectors of the corresponding adjacency matrix [31,32]. Here, denotes the average over an ensemble of random graphs. This is a notable result because it puts forward the application of degree-based TIs beyond mathematical chemistry (it is relevant to add that random graphs have also been studied by means of eigenvalue-based TIs such as the Estrada index, the Laplacian Estrada index, and Rodriguez-Velazquez indices, see e.g., [33,34].) Specifically, on the one hand, certain average TIs can provide equivalent information than traditional, spectrum-based, random matrix theory measures. On the other hand, TIs may be used to predict spectral properties of random graphs.

    Therefore, motivated by potential applications, in what follows we apply for the first time (to our knowledge) Revan degree-based indices on random graphs. Indeed, since the inequalities obtained in Section 2 are not restricted to any particular type of graph, we first quantify them for random graphs and later we extend some of them to index average values, as needed in statistical studies of random graphs.

    Below we consider two prominent models of random graphs: Erdös-Rényi (ER) graphs and random geometric (RG) graphs. ER graphs [35,36] GER(n,p) are formed by n vertices connected independently with probability p[0,1]. While RG graphs [37,38] GRG(n,r) consist of n vertices uniformly and independently distributed on the unit square, where two vertices are connected by an edge if their Euclidean distance is less or equal than the connection radius r[0,2].

    In the Theorems presented in Section 2 we know the graph properties needed for the equalities involved have to be attained, however we do not really know how strong the inequalities could be for arbitrary graphs. Thus, in this subsection we quantify the inequalities of those Theorems on random graphs.

    To ease the quantification of the six Theorems presented in Section 2 we:

    (i) write the left inequality of Corollary 2 as

    0RSO(G)2ΔR2(G), (3.1)

    (ii) write Corollary 4 as

    0mΔR1(G)RSO(G), (3.2)

    (iii) write the right inequality in the second line of Corollary 7, for α=1/2, as

    0<RKA11/2,2(G)RSO(G), (3.3)

    (iv) write the left inequality in the second line of Theorem 8, for α=1/2, as

    0<M1(G)2(Δ+δ)m+RKA11/2,2(G), (3.4)

    (v) write the left inequality of Corollary 10 as

    0RSO(G)δΔSO(G), (3.5)

    and

    (vi) write the left inequality of Corollary 12 as

    0RSO(G)2δΔH(G). (3.6)

    Then, in Figure 1 [in Figure 2] we plot the right hand side of Eqs (3.1)–(3.6) as a function of the probability p [the connection radius r] of ER [RG] graphs of size n. For each value of p [r] we present results for 10 graph realizations.

    Figure 1.  Right hand side of (a) Eq (3.1), (b) Eq (3.2), (c) Eq (3.3), (d) Eq (3.4), (e) Eq (3.5) and (f) Eq (3.6) as a function of the probability p of Erdős-Rényi graphs GER(n,p) of size n. Each symbol corresponds to a single graph realization. Dashed lines in (c) [(d)] correspond to Eq (3.7) [(3.8)] with dG=(n1)p.
    Figure 2.  Right hand side of (a) Eq (3.1), (b) Eq (3.2), (c) Eq (3.3), (d) Eq (3.4), (e) Eq (3.5) and (f) Eq (3.6) as a function of the connection radius r of random geometric graphs GRG(n,r) of size n. Each symbol corresponds to a single graph realization. Dashed vertical lines in (a, b) [(e, f)] mark the maximum data values which occur at r0.8 [r0.64]. Dashed lines in (c) [(d)] correspond to Eq (3.7) [Eq (3.8)] with dG=(n1)r2(π8r/3+r2/2).

    For both random graph models, ER and RG graphs, we observe two different behaviors for the right hand side of Eqs (3.1)–(3.6): On the one hand, the right hand side of Eqs (3.1), (3.2), (3.5) and (3.6) are close to zero for n0 and p0 or r0, then grow with p or r, approach a maximum value at p=pmax or r=rmax, and finally decrease abruptly and become zero at p=1 or r=2. For ER graphs pmax is very close to 1, so the decrease of the right hand side of Eqs (3.1), (3.2), (3.5) and (3.6) is not displayed in Figure 1; while for RG graphs rmax0.8 for Eqs (3.1) and (3.2) and rmax0.64 for Eqs (3.5) and (3.6), as can be seen in Figure 2. On the other hand, the right hand side of Eqs (3.3) and (3.4) is close to zero for n0 and p0 or r0, then grow with p or r approaching their maximum values at p=1 or r=2. That is, the right hand side of Eqs (3.3) and (3.4) are increasing functions of n and p or r.

    Moreover, we can estimate the right hand side of Eqs (3.3) and (3.4) as follows. As it is shown in Ref. [39], for ER and RG graphs in the dilute limit (i.e., when dG1) we have that δΔdGrG; where dG and rG are the average degree and average Revan degree, respectively. Therefore, for ER and RG graphs in the dilute limit, we can write the right hand side of Eq (3.3) as

    RKA11/2,2(G)RSO(G)=uvE(G)(rG(u)1/2+rG(v)1/2)2uvE(G)rG(u)2+rG(v)2uvE(G)(2rG1/2)2uvE(G)2rG2=4rGm2rGm4dGm2dGm=422ndG2; (3.7)

    while the right hand side of Eq (3.4) can be written as

    M1(G)2(Δ+δ)m+RKA11/2,2(G)=uvE(G)(dG(u)+dG(v))2(Δ+δ)m+uvE(G)(rG(u)1/2+rG(v)1/2)2uvE(G)2dG4dGm+uvE(G)(2rG1/2)2=2dGm+4rGm2dGm+4dGm=2dGm=ndG2. (3.8)

    Indeed, in panels (c) and (d) of Figures 1 and 2 we plot (as dashed lines) Eqs (3.7) and (3.8), respectively, and observe a very good correspondence with the computational data, even for relatively small values of p or r. In Figure 1(c), (d) we used dG=(n1)p for ER graphs; while in Figure 2(c), (d) we used dG=(n1)f(r), with f(r)=r2(π8r/3+r2/2), for RG graphs.

    Since there is a growing interest in the average values of TIs, in particular when applied to random graphs (see e.g., [28,29,30,31]), it should be desirable to state inequalities for them. Fortunately, due to linearity, Theorem 5 can be straightforwardly extended to index average values; it reads as:

    Corollary 13. If {G} is an ensemble of graphs and a,b,cR{0}, then

    RKA1ab/c,c(G)<RKA1a,b(G)2bcRKA1ab/c,c(G) if b>c,bc>0,2bcRKA1ab/c,c(G)RKA1a,b(G)<RKA1ab/c,c(G) if b<c,bc>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b<0,c>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b>0,c<0.

    Moreover, from Theorem 8 we make the following conjecture:

    Conjecture 14. If {G} is an ensemble of graphs, with average maximum degree Δ, average minimum degree δ, average number of edges m, and αR{0}, then

    2(Δ+δ)m211/αRKA1α,1/α(G)M1(G) <2(Δ+δ)mRKA1α,1/α(G) if α>1,2(Δ+δ)mRKA1α,1/α(G)<M1(G) 2(Δ+δ)m211/αRKA1α,1/α(G) if 0<α<1,M1(G)2(Δ+δ)m211/αRKA1α,1/α(G) if α<0.

    Now, to validate Conjecture 14 we split and write the second inequality, with α=1/2, as

    0<M1(G)2(Δ+δ)m+RKA11/2,2(G) (3.9)

    and

    02(Δ+δ)m21RKA11/2,2(G)M1(G). (3.10)

    Then, in Figure 3 [in Figure 4] we plot the right hand side of Eqs (3.9) and (3.10) as a function of the probability p [the connection radius r] of ER [RG] graphs of size n. The averages are computed over ensembles of 106 random graphs. As can be clearly observed in these figures, the right hand side of Eqs (3.9) and (3.10) is larger than zero for all the combinations of n and p [r] used in this work; thus we validate Conjecture 14 on both ER and RG graphs. In Figure 3(c), (d) we compute m as m=n(n1)p/2, while in Figure 4(c), (d) we use m=n(n1)f(r)/2 with f(r)=r2(π8r/3+r2/2).

    Figure 3.  Right hand side of (a) Eq (3.9) and (b) Eq (3.10) as a function of the probability p of Erdős-Rényi graphs GER(n,p) of size n. Each data value was computed by averaging over 106 random graphs GER(n,p). Dashed lines in (a) are Eq (3.11) with dG=(n1)p. The dashed vertical line in (b) marks the maximum of the curves which occur at p0.5.
    Figure 4.  Right hand side of (a) Eq (3.9) and (b) Eq (3.10) as a function of the connection radius r of random geometric graphs GRG(n,r) of size n. Each data value was computed by averaging over 106 random graphs GRG(n,r). Dashed lines in (a) are Eq (3.11) with dG=(n1)r2(π8r/3+r2/2). The dashed vertical line in (b) marks the maximum of the curves which occur at r0.64.

    Also, from Figure 3(a) [from Figure 4(a)] we observe that the right hand side of Eq (3.9) is close to zero for n0 and p0 [r0], then grows with p [r] approaching their maximum value at p=1 [r=2]. While from Figure 3(b) [from Figure 4(b)] we note that the right hand side of Eq (3.10) is close to zero for n0 and p0 [r0], then grows with p [r], approaches a maximum value at p0.5 [r0.64], and finally decrease abruptly and become zero at p=1 [r=2].

    By the use of the same approximations that allowed us to get Eqs (3.7) and (3.8) in the previous subsection, we can also estimate the right hand side of Eq (3.9) for both ER and RG graphs in the dilute limit. This results in

    M1(G)2(Δ+δ)m+RKA11/2,2(G)ndG2. (3.11)

    Therefore, in Figure 3(a) [from Figure 4(a)] we plot Eq (3.11), as dashed lines, and observe good correspondence with the computational data.

    In this work, we study the recently introduced Revan Sombor index and the first and second Revan (a,b)-KA indices. First, we present new relations to provide bounds on these Revan indices of the Sombor family which also serve to relate them with other Revan indices (such as R1,2, which are Revan versions of the first and second Zagreb indices) and with standard degree-based indices (such as the Sombor index, the first and second (a,b)-KA indices, the first Zagreb index and the Harmonic index). We also quantify our relations on two models of random graphs: Erdös-Rényi graphs and random geometric graphs. Then, motivated by the growing interest in the application of topological indices in the analysis of random graphs, we extend some of our relations to index average values, so they can be effectively used for the statistical study of ensembles of random graphs. It is important to stress that our results on index average values are essentially obtained by averaging over the graph ensemble in the statistical sense rather than in the rigorous probabilistic sense, which is a topic we plan to explore in the near future.

    We hope our results may motivate further analytical as well as computational studies of Revan Sombor indices. Indeed, since the average Sombor indices have shown to work as complexity measures of random graphs and networks (equivalent to random matrix theory measures) [32], we plan to explore wether Revan Sombor indices could also be used in this direction.

    J. A. M. -B. thanks support from CONACyT (Grant No. 286633), CONACyT-Fronteras (Grant No. 425854), and VIEP-BUAP (Grant No. 100405811-VIEP2022), Mexico. The research of J. M. R. and J. M. S. was supported by a grant from Agencia Estatal de Investigación (PID2019-106433GBI00/AEI/10.13039/501100011033), Spain. J. M. R. was supported by the Madrid Government (Comunidad de Madrid-Spain) under the Multiannual Agreement with UC3M in the line of Excellence of University Professors (EPUC3M23), and in the context of the V PRICIT (Regional Programme of Research and Technological Innovation).

    The authors declare there is no conflict of interest.



    [1] J. Ferlay, H. R. Shin, F. Bray, D. Forman, C Mathers, D. M. Parkin, Estimates of worldwide burden of cancer in 2008: Globocan 2008, Int. J. Cancer, 27 (2010), 2893–2917. https://doi.org/10.1002/ijc.25516 doi: 10.1002/ijc.25516
    [2] K. M. Ratheesh, L. K. Seah, V. M. Murukeshan, Spectral phase-based automatic calibration scheme for swept source-based optical coherence tomography systems, Phy. Med. Biol., 21 (2016), 7652. https://doi.org/10.1088/0031-9155/61/21/7652 doi: 10.1088/0031-9155/61/21/7652
    [3] R. K. Meleppat, M. V. Matham, L. K. Seah, An efficient phase analysis-based wavenumber linearization scheme for swept source optical coherence tomography systems, Laser Phys. Lett., 5 (2015), 055601. https://doi.org/10.1088/1612-2011/12/5/055601 doi: 10.1088/1612-2011/12/5/055601
    [4] R. K. Meleppat, M. V. Matham, L. K. Seah, Optical frequency domain imaging with a rapidly swept laser in the 1300nm bio-imaging window, in International Conference on Optical and Photonic Engineering (ICOPEN 2015), International Society for Optics and Photonics, (2015), 9524: 95242R. https://doi.org/10.1117/12.2190530
    [5] N. Mu, H. Wang, Y. Zhang, J. Jiang, J. Tang, Progressive global perception and local polishing network for lung infection segmentation of COVID-19 CT images, Pattern Recognit., 120 (2021), 108168. https://doi.org/10.1016/j.patcog.2021.108168 doi: 10.1016/j.patcog.2021.108168
    [6] F. Zhu, Z. Gao, C. Zhao, Z. Zhu, J. Tang, Y. Liu, et al., Semantic segmentation using deep learning to extract total extraocular muscles and optic nerve from orbital computed tomography images, Optik, 244 (2021), 167551. https://doi.org/10.1016/j.ijleo.2021.167551 doi: 10.1016/j.ijleo.2021.167551
    [7] C. Zhao, Y. Xu, Z. He, J. Tang, Y. Zhang, J Han, et al., Lung segmentation and automatic detection of COVID-19 using radiomic features from chest CT images, Pattern Recognit., 119 (2021), 108071. https://doi.org/10.1016/j.patcog.2021.108071 doi: 10.1016/j.patcog.2021.108071
    [8] O. Ronneberger, P. Fischer T. Brox, U-net: convolutional networks for biomedical image segmentation, in International Conference on Medical image computing and computer assisted intervention, (2015), 234–241. https://doi.org/10.1007/978-3-319-24574-4_28
    [9] O. Oktay, J. Schlemper, L. L. Folgoc, M. Lee, M. Heinrich, K. Misawa, et al., Attention u-net: learning where to look for the pancreas, preprint, arXiv: 1804.03999.
    [10] L. C. Chen, G. Papandreou, I. Kokkinos, K. Murphy, A. L Yuille, Deep lab: semantic image segmentation with deep convolutional nets, atrous convolution, and fully connected CRFs, IEEE Trans. Pattern Anal. Mach. Intell., 40 (2017), 834–848. https://doi.org/10.1109/TPAMI.2017.2699184 doi: 10.1109/TPAMI.2017.2699184
    [11] E. Shelhamer, J. Long, T. Darrell, Fully convolutional networks for semantic segmentation, IEEE Trans. Pattern Anal. Mach. Intell, 39 (2016), 640–651. https://doi.org/10.1109/TPAMI.2016.2572683 doi: 10.1109/TPAMI.2016.2572683
    [12] L. C. Chen, G. Papandreou, I. Kokkinos, K. Murphy, A. L. Yuille, Semantic image segmentation with deep convolutional nets and fully connected crfs, preprint, arXiv: 1412.7062.
    [13] M. Z. Alom, M. Hasan, C. Yakopcic1, T. M. Taha, V. K. Asari1, Recurrent residual convolutional neural network based on U-Net (R2U-Net) for nedical image segmentation, preprint, arXiv: 1802.06955.
    [14] J. Wang, P. Lv, H. Wang, C. Shi, SAR-U-Net: squeeze-and-excitation block and atrous spatial pyramid pooling based residual U-Net for automatic liver segmentation in Computed Tomography, Comput. Methods Programs Biomed., 208 (2021), 106268. https://doi.org/10.1016/j.cmpb.2021.106268 doi: 10.1016/j.cmpb.2021.106268
    [15] Z. Zhou, M. M. Rahman Siddiquee, N. Tajbakhsh, J. Liang, Unet++: a nested u-net architecture for medical image segmentation, in Deep Learning in Medical Image Analysis and Multimodal Learning for Clinical Decision Support, Springer, Cham, (2018), 3–11.
    [16] F. Milletari, N. Navab, S. A. Ahmadi, V-net: fully convolutional neural networks for volumetric medical image segmentation, in 2016 IEEE Fourth International Conference on 3D Vision (3DV), (2016), 565–571. https://doi.org/10.1109/3DV.2016.79
    [17] X. Li, H. Chen, X. Qi, Q. Dou, C. W. Fu, P. A. Heng, H-DenseUNet: hybrid densely connected UNet for liver and tumor segmentation from CT volumes, IEEE Trans. Med. Imaging, 37 (2018), 2663–2674. https://doi.org/10.1109/TMI.2018.2845918 doi: 10.1109/TMI.2018.2845918
    [18] R. Mehta, J. Sivaswamy, M-net: A convolutional neural network for deep brain structure segmentation, in 2017 IEEE 14th International Symposium on Biomedical Imaging (ISBI 2017), (2017), 437–440. https://doi.org/10.1109/ISBI.2017.7950555
    [19] C. Yu, J. Wang, C. Peng, C. Gao, G. Yu, N. Sang, Bisenet: bilateral segmentation network for real-time semantic segmentation, in Proceedings of the European Conference on Computer Vision (ECCV), (2018), 325–341. https://doi.org/10.1007/978-3-030-01261-8_20
    [20] Z. Gu, J. Cheng, H. Fu, K. Zhou, H. Hao, Y. Zhao, et al., Ce-net: Context encoder network for 2d medical image segmentation, IEEE Trans. Med. Imaging, 38 (2019), 2281–2292. https://doi.org/10.1109/TMI.2019.2903562 doi: 10.1109/TMI.2019.2903562
    [21] S. Wiesler, H. Ney, A convergence analysis of log-linear training, Adv. Neural Inf. Process. Syst., 24 (2011), 657–665.
    [22] E. Vorontsov, A. Tang, C. Pal, S Kadoury, Liver lesion segmentation informed by joint liver segmentation, in 2018 IEEE 15th International Symposium on Biomedical Imaging (ISBI 2018). IEEE, (2018), 1332–1335. https://doi.org/10.1109/ISBI.2018.8363817
    [23] L. Zhou, C. Zhang, M. Wu, D-LinkNet: LinkNet with pretrained encoder and dilated convolution for high resolution satellite imagery road extraction, in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition Workshops, (2018), 182–186. https://doi.org/10.1109/CVPRW.2018.00034
  • 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(3865) PDF downloads(152) Cited by(13)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog