Next Article in Journal
Initial Boundary Value Problem for a Fractional Viscoelastic Equation of the Kirchhoff Type
Next Article in Special Issue
Fault Diagnosis of Hydroelectric Units Based on a Novel Multiscale Fractional-Order Weighted Permutation Entropy
Previous Article in Journal
Mild Solution for the Time-Fractional Navier–Stokes Equation Incorporating MHD Effects
Previous Article in Special Issue
A Study of Adaptive Fractional-Order Total Variational Medical Image Denoising
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Adaptive Fractional Differential Active Contour Image Segmentation Method

1
School of Automation and Electrical Engineering, Shenyang Ligong University, Shenyang 110159, China
2
Science and Technology on Electro-Optical Information Security Control Laboratory, Tianjin 300308, China
*
Author to whom correspondence should be addressed.
Submission received: 23 August 2022 / Revised: 4 October 2022 / Accepted: 6 October 2022 / Published: 10 October 2022

Abstract

:
When the image is affected by strong noise and uneven intensity, the traditional active contour models often cannot obtain accurate results. In this paper, a novel adaptive fractional differential active contour image segmentation method is proposed to solve the above problem. At first, in order to extract more texture parts of the image, an adaptively fractional order matrix is constructed according to the gradient information of the image, varying the fractional order of each pixel. Then, the traditional edge-stopping function in the regularization term is susceptible to noise, and a new fractional-order edge-stopping function is designed to improve noise resistance. In this paper, a fitting term based on adaptive fractional differentiation is introduced to solve the problem of improper selection of the initial contour position leading to inaccurate segmentation results so that the initial contour position can be selected arbitrarily. Finally, the experimental results show that the proposed method can effectively improve the segmentation accuracy of noise images and weak-edge images and can arbitrarily select the position selection of the initial contour.

1. Introduction

The active contour model (ACM) [1] has the advantages of sub-pixel accuracy, topology adaptability, etc. Therefore, it has been widely used in the fields of image segmentation [2,3,4,5], video clipping [6,7], scene understanding [8], and object tracking [9,10]. Image segmentation based on ACM is a nonlinear segmentation method. Its initial contour is induced to move towards the image boundary, and the target area in the image domain is finally extracted [11]. However, in practical applications, images are often affected by noise and nonuniformity of intensity, resulting in over-segmentation or local minimum problems. Therefore, it is necessary to optimize the image segmentation method based on ACM to improve the image segmentation accuracy.
ACMs are roughly divided into two categories: edge-based models [12,13,14,15] and region-based models [16,17,18,19,20]. The edge-based model mainly defines the edge-stopping function according to the gradient information near the active contour curve. This method can handle the image with a clear object boundary well, but the segmentation results of the image with blurred boundary are generally unsatisfactory. In addition, it is easy to fall into a local minimum due to factors such as noise and texture. The corresponding region-based model pays more attention to global statistical information so that it can segment images with weak edges and is robust to the initial contour position and noise. However, it is difficult to obtain ideal segmentation results for some images with non-uniform gray levels in the target or background area.
After years of research, scholars at home and abroad have proposed numerous ACMs. The following were representative ACMs. In 1993, Caselles et al. [21] first introduced the level set method into ACM and proposed the geodesic active contour (GAC) model. The basic idea of this model is to transform the problem of image segmentation into the problem of finding the shortest geodesic, in which the weight is determined by the image gradient at the position of the curve. Li et al. [22] designed a distance regularization term to be added to the energy function of the GAC model to solve the problem of re-initializing the level set function in the iterative process, thus ensuring the stability of the level set function evolution. Chan and Vese [23] proposed the Chan–Vese (CV) model, which is insensitive to noise as well as the position of the initial contour. However, it is difficult for this model to obtain ideal segmentation results for some images with uneven gray levels in the target or background areas. Li et al. [24] proposed a variable region fitting (RSF) model by embedding local image information. The RSF model greatly improves the segmentation performance of images with uneven intensity.
To further improve the segmentation performance of noise and weak edge images, scholars at home and abroad have found that fractional differential plays a special role in signal singularity detection and extraction [25,26,27,28]. It can improve the high-frequency component of the signal while retaining the low-frequency component of the signal nonlinearly. Ref. [12] proposed a fractional edge-stopping function using a fractional gradient instead of a conventional gradient. Ref. [18] proposed a fractional active contour model by applying different types of noise to images. Ref. [29] introduced new edge energy, which is driven by the difference between the inward and outward fractional differential of the evolution profile. The above fractional order models have advantages over traditional ACM. However, these models perform fixed order differentiation on the image, which means that the texture part of the image is easy to be smoothed or destroyed due to improper order, and the experiment needs to be repeated to obtain the desired order. Refs. [30,31] proposed a fractional differential model with adaptive fractional order, and its segmentation performance is improved to a certain extent compared with the fixed order fractional ACM. In the research process of adaptive fractional order, Refs. [32,33,34] provides research strategies for nonlinear parameter optimization, to effectively and quickly determine the optimal value of parameters. Refs. [35,36] adaptively selects fractional order according to the noise standard deviation of the source image, which not only retains the gradient feature but also makes up for the lack of ignoring pixels. At present, ACMs based on adaptive fractional order show improvement, but they only rely on the edge or region information of the image, and still have the problem of instability. To solve these problems, this paper proposes a new active contour image segmentation method based on the edge and region information of the image to achieve accurate segmentation of noisy images and weak edge images.
Based on these previous works, we propose an adaptive fractional differential active contour image segmentation method. The main contributions of the paper are summarized as follows:
  • According to the gradient information of the image, an adaptively changing order matrix is constructed, and the fractional order of each pixel of the image is adaptively changed. Different differential orders are used in the edge texture and smooth area of the image, which preserves richer image information.
  • A fractional-order edge-stopping function is proposed, which can effectively improve the performance of the segmentation of noisy images. It solves the problem that the traditional edge-stopping function is trapped in a local minimum due to the influence of noise points.
  • The fitting term is designed by combining adaptive fractional differential and Gaussian kernel function, which makes the new method effectively overcome the problem of sensitivity to the initial contour position.
The rest of the structure of the paper is as follows. Section 2 introduces the related knowledge of the classical active contour model. Section 3 proposes an adaptive fractional differentiation and fractional edge-stopping function. Section 4 proposes an adaptive fractional derivative active contour image segmentation method. The experimental results are given in Section 5 to illustrate the effectiveness of the algorithm. Finally, the conclusions are given in Section 6.

2. Active Contour Classic Model

In this section, we briefly introduce two classical active contour models, the edge-based distance-regularized level set evolution (DRLSE) model [37] and the region-based region scalable fitting (RSF) model [24].

2.1. Distance Regularized Level Set Evolution (DRLSE) Model

Let Ω be the image area, and ϕ ( x , y ) be the level set function defined on Ω . The energy function of the DRLSE model in [37] is defined as follows,
E ( ϕ ) = μ Ω R ( | ϕ | ) d x + λ Ω g ( x , y ) δ ( ϕ ) | ϕ | d x d y + v Ω g ( x ) H ( ϕ ) d x d y ,
where μ > 0 , λ > 0 , and ν R are the weight coefficients of the respective energy terms. R ( | ϕ | ) is a double-well potential function. g ( x , y ) is the edge-stopping function. δ ( ) and H ( ) are the Dirac function and the Heaviside function respectively.
The first term in the Equation (1) is a penalty term, which essentially maintains the regularity of the level set function in the evolution process; the second term is the long term, and the line integral of the edge-stopping function on the contour of the zero level set is calculated; the third term is the area term, which calculates the weighted area of the area within the active contour. The second and third terms are used to drive the curve to evolve towards the target boundary.
The integrand used in this model is,
R ( | ϕ | ) = { 1 ( 2 π ) 2 [ 1 - cos ( 2 π | ϕ | ) ] , | ϕ | < 1 1 2 ( | ϕ | 1 ) 2 , | ϕ | 1 .
In addition, for most edge-based active contour models [12,13,14,15], the edge-stopping function g ( x , y ) is defined as follows,
g ( x , y ) = 1 1 + | G σ I ( x , y ) | 2 ,
where G σ is a Gaussian kernel function with standard deviation σ , which is used to smooth the image to reduce the effect of noise. I ( x , y ) is image information.
δ ( ) and H ( ) are defined in [37] as follows,
δ ε ( x ) = { 1 2 ε [ 1 + cos ( π x ε ) ] , | x | < ε 0 , | x | > ε ,
H ε ( x ) = { 1 2 ( 1 + x ε + 1 π sin ( π x ε ) ) , | x | ε 1 , x > ε 0 , x < ε ,
where ε is the smoothing parameter, and its value is usually set to 1.5.
In Equation (1), E ( ϕ ) is a function of the level set function ϕ , and the optimal ϕ is found by minimizing the energy function E ( ϕ ) . Using the variational method and according to the gradient descent flow theory [38], the gradient descent flow equation of the level set function ϕ is obtained as follows,
ϕ t = μ div ( d R | ϕ | ϕ ) + λ δ ε ( ϕ ) div ( g ϕ | ϕ | ) + v g δ ε ( ϕ ) ,
where div is the divergence operator, and the function d R is defined as,
d R ( | ϕ | ) = R ( | ϕ | ) | ϕ | = { sin ( 2 π | ϕ | ) 2 π | ϕ | , | ϕ | < 1 1 1 | ϕ | , | ϕ | 1 .
Compared with the traditional edge active contour model, the DRLSE model can accurately and stably segment objects with a clear boundary. However, when segmenting a noisy image, the evolution curve tends to fall into a local minimum due to noise interference; relying only on edge information, the result is sensitivity to the position of the initial contour. In addition, at the weak boundary of the image, the gradient of the image is small, which makes the value of the edge-stopping function larger, and the evolution curve easily crosses the target boundary.

2.2. Region Scalable Fitting (RSF) Model

The RSF model [24] uses the Gaussian kernel function to define the energy functional, and the expression is as follows,
E ( ϕ ) = λ 1 Ω K σ ( x y ) ( I ( y ) f 1 ( x ) ) 2 H ( ϕ ( y ) ) d x d y + λ 2 Ω K σ ( x y ) ( I ( y ) f 2 ( x ) ) 2 ( 1 H ( ϕ ( y ) ) ) d x d y + μ Ω 1 2 ( | ϕ | 1 ) 2 d x + ν Ω | H ( ϕ ( x ) ) | d x
where K σ ( x ) is a Gaussian kernel function with a standard deviation of σ , f 1 ( x ) , f 2 ( x ) are the fitting means of the given pixel x in the local area inside and outside the contour curve respectively, and the expression is,
f i ( x ) = K σ ( x ) [ H i ( ϕ ( x ) ) I ( x ) ] K σ ( x ) H i ( ϕ ( x ) )         i = 1 , 2 ,
where H 1 ( ϕ ( x ) ) = H ( ϕ ( x ) ) , H 2 ( ϕ ( x ) ) = 1 H ( ϕ ( x ) ) .
The first term in the Equation (8) is a local fitting term, which is used to drive the evolution curve to stop at the target contour; the second term is the penalty term to prevent the re-initialization of the evolution curve; the third term is the regular term, which controls the length of the evolution curve.
The gradient descent flow equation of the level set function ϕ is obtained by minimizing the energy functional (8) as follows,
ϕ t = δ ( x ) ( λ 1 K σ ( x y ) ( I ( y ) f 1 ( x ) ) 2 d y + λ 2 K σ ( x y ) ( I ( y ) f 2 ( x ) ) 2 d y + μ ( 2 ϕ div ( ϕ | ϕ | ) ) + ν δ ( x ) div ( ϕ | ϕ | )
The RSF model greatly improves the segmentation performance of the image with uneven intensity. Because only the local information of the image is used, the model is very sensitive to noise. In addition, the model is prone to fall into local optimization and is sensitive to the initial position of the evolution curve. To obtain the ideal segmentation result, it is necessary to select the appropriate initial position according to the segmentation result, which limits its application.
Through the above brief analysis of DRLSE and RSF models, it can be concluded that both models are sensitive to the position of the initial contour and noise. In addition, this is not applicable to the target edge blurred image. Therefore, it is necessary to optimize the active contour image segmentation method.

3. Adaptive Fractional Differential

The fractional differential has the characteristic of increasing the high-frequency component of the signal while retaining the low-frequency component of the signal nonlinearly. Therefore, fractional differentiation can increase the signal size of the image edge and texture part, while retaining the regional information where the gray value changes slowly, and has good robustness to noise. Therefore, fractional differentiation has received more and more attention and applications in the field of image processing.

3.1. Fractional Differential

Fractional order differential is a branch of integer order differential, which extends the traditional integer order differential. There are three classical fractional differential definitions, namely Riemann–Liouville (R-L), Grüwald–Letnikov (G-L), and Caputo [25,26,27,28]. R-L and Caputo are both defined by the Cauchy integral formula, which is complex and is not conducive to large-scale data calculation. The definition referenced by G-L can be converted into convolution form, so the former has a better effect on image processing than the other two definitions.
G-L of function f ( t ) is defined as follows,
D t α a G L f ( t ) = lim h 0 h α j = 0 [ ( t a ) / h ] ( 1 ) j ( α j ) f ( t j h ) h α j = 0 [ ( t a ) / h ] ω j ( α ) f ( t j h )
where ( α j ) represents the combination number parameter.
For n × m a two-dimensional image f ( x , y ) , the definition of fractional differential in the x - and the y - directions are as follows:
α f ( x , y ) x α f ( x , y ) + ( α ) f ( x 1 , y ) + ( α ) ( α + 1 ) 2 f ( x 2 , y ) + + Γ ( α + 1 ) n ! Γ ( α + n + 1 ) f ( x n , y ) +
α f ( x , y ) y α f ( x , y ) + ( α ) f ( x , y 1 ) + ( α ) ( α + 1 ) 2 f ( x , y 2 ) + + Γ ( α + 1 ) n ! Γ ( α + n + 1 ) f ( x , y n ) +

3.2. Adaptive Fractional Differential

The traditional image processing method, based on the fractional differential, is to use the same order differential operation for each pixel in the image. This is often inconsistent with the characteristics of different pixel values so that the texture details of the image are smoothed or destroyed, and eventually, the image information is lost. It will seriously affect the image quality and further analysis and processing.
In this paper, we propose a new fractional-order varying method for each pixel and construct an adaptive fractional order matrix, P , according to the gradient information of the image. So, different differential orders are adopted in the smooth part and the edge texture part of the image to obtain more abundant image information. The expression of the adaptive order matrix P is,
P = | I ( x , y ) | + α ones | I ( x , y ) | + β ones ,
where I ( x , y ) is a two-dimensional image of n × m , | I ( x , y ) | is the gradient information of the image I ( x , y ) , and the dimension is n × m , the ones represents a matrix whose elements are all 1 and its dimension is n × m , and α and β are the correction parameters. The order interval of the fractional differential operation is controlled by the value of the correction parameter, and the value depends on the specific situation in different images.
The experimental object of this paper is mainly infrared images. Because the infrared image is obtained by “measuring” the heat radiated from the object, it usually has the characteristics of poor resolution, low contrast, large background noise, and blurred visual effect. Through many experiments and relevant literature references, combined with the characteristics of infrared images, the interval of differential order and the constraint conditions of modified parameters will be analysed later.
The fractional order matrix P can be adaptively adjusted according to the local statistical information and structural features of the image. Make it have a larger differential order at the strong edges of the image and a smaller differential order at the weak edges and textures of the image. Taking the template of image grayscale information 5 × 5 as an example, the correction parameters are α = 1 , β = 0.8 , and the corresponding order changes are shown in Figure 1.
In this paper, the fractional-order variable differential of the image I ( x , y ) can be obtained by using the definition of G-L, and the results of the fractional-order differential operation with different orders for each pixel point are as follows,
D P G L I ( x , y ) = ( D P 11 I ( 1 , 1 ) D P 12 I ( 1 , 2 ) D P 1 m I ( 1 , m ) D P 21 I ( 2 , 1 ) D P 22 I ( 2 , 2 ) D P 2 m I ( 2 , m ) D P n 1 I ( n , 1 ) D P n 2 I ( n , 2 ) D P n m I ( n , m ) ) R n × m ,
where D P G L I ( x , y ) is fractional order changing image, P n m is the element in the fractional order matrix P .
Since the fractional derivative is a linear operator, in this paper, the absolute value operations in x , y directions are used to ensure that the linear change of the image grayscale information is preserved. The image I ( x , y ) processed by adaptive fractional differentiation can be updated as,
D P G L I ( x , y ) = ( D x P I ( x , y ) , D y P I ( x , y ) ) ,
| D P G L I ( x , y ) | = | D x P I ( x , y ) | + | D y P I ( x , y ) | .

3.3. Fractional Edge-Stopping Function

In the edge-based active contour model, the edge stop function plays an important role in image segmentation. The traditional edge stop function can be a non-negative function with monotone decreasing, and its main function is to control the evolution speed of the curve. However, because it only depends on the gradient information of the image, it is easy to be affected by noise and can not accurately segment the noisy image and the target edge blurred image.
To improve the robustness of the traditional edge-stopping function to noise, this paper introduces the adaptive fractional differential into the Equation (3), defines a new edge-stopping function, and proposes an edge-stopping function based on the adaptive fractional differential,
g ( x , y ) = 1 1 + G σ | D P G L I ( x , y ) | ,
where | D P G L I ( x , y ) | is the result of performing the fractional-order variable operation on each pixel point of the image I ( x , y ) . The calculation method is shown in Equation (17).

4. Image Segmentation Based on Adaptive Fractional Differentiation

In this section, combined with Section 3, an image segmentation method based on adaptive fractional differentiation is proposed. The energy functional of this method consists of three parts: penalty term E P , regularization term E R , and fitting term E F , which will be described in detail below,
E = E P + E R + E F .

4.1. Penalty Term E P

The first term in the Equation (19) is the penalty term, and its function is to avoid the re-initialization step of ACM. The energy function is as follows,
E P = μ Ω R ( | ϕ | ) d x ,
where μ > 0 is the weight coefficient, and ϕ is the level set function. As defined in Equation (2), this is a double-well potential function.
In Ref [22], the accumulated function is the symbol distance function (SDF), and its expression is p = ( | ϕ | 1 ) 2 / 2 . The diffusion rate of the integrand can be obtained from the gradient flow equation, and its expression is r 1 = 1 ( 1 / | ϕ | ) . When | ϕ | = 1 , the level set function is the target boundary, and the diffusion rate r 1 = 0 is stable; when | ϕ | = 0 , the level set function be far from the target boundary, and the diffusion rate is r 1 . This will make the level set function more and more irregularly in the iterative process, resulting in oscillation during the evolution process. As a result, the calculation accuracy will be seriously affected, and the stability of the level set evolution will be destroyed. Therefore, the integrand proposed in [37] is introduced in this paper, and the diffusion rate of the integrand (2) can be obtained from the gradient flow equation,
r 2 = { ( | ϕ | 1 ) ( 2 | ϕ | 1 ) , | ϕ | < 1 1 1 | ϕ | , | ϕ | 1 .
It can be seen from Equation (21) that when | ϕ | = 1 , r 1 = 0 is in a stable state; when | ϕ | = 0 , r 1 = 1 is also in a stable state. Therefore, under the action of the double-well potential function, the final level set function is relatively smoother, as shown in Figure 2.
The gradient descent flow equation of the level set function ϕ obtained by minimizing the energy functional (20) is as follows,
ϕ t = μ div ( d R | ϕ | ϕ ) ,
where the function d R is defined as Equation (7).

4.2. Regularization Term E R

The second term in the Equation (19) is a regular term, which is used to maintain the regularity of the evolution curve and calculate the line integral of the curve. The energy function is as follows,
E R = ν Ω g ( x , y ) ( | H ( ϕ ) | ) d x d y ,
where ν > 0 is the weight coefficient, and g ( x , y ) is the fractional edge-stopping function, as shown in Equation (18) and described in detail in Section 3.
The gradient downflow equation of the level set function ϕ is obtained by minimizing the energy functional (23) as follows,
ϕ t = ν g δ ( ϕ ) div ( ϕ | ϕ | ) .
Due to the introduction of a fractional-order edge-stopping function in this term, the model has a certain robustness to noise. The improved fractional-order edge-stopping function only relies on the gradient information of the image and cannot solve the problem that the model is sensitive to the initial contour position. Therefore, in Section 4.3, the adaptive fractional derivative is introduced into the region fitting term, and the use of the fractional derivative can enhance the advantages of the image texture details and optimize the segmentation method.

4.3. Fitting Term E F

The third term in the Equation (19) is a fitting term, which is used to drive the evolution curve to stop at the target contour and make the model robust to the initial contour position. The fitting term proposed in this paper includes two parts: one is the local region fitting term introduced from [24], and the other is the fitting term based on the adaptive fractional differential. The energy functional is as follows,
E F = λ [ Ω K σ ( x y ) ( I ( y ) f 1 ( x ) ) 2 H ( ϕ ( y ) ) d x d y + Ω K σ ( x y ) ( I ( y ) f 2 ( x ) ) 2 ( 1 H ( ϕ ( y ) ) ) d x d y ] + E F F
where f 1 ( x ) and f 2 ( x ) are respectively the fitting mean of the local area inside and outside the contour curve of a given pixel point x , and the expression is Equation (9). K σ ( x ) is the Gaussian kernel function with standard deviation σ , I ( y ) is the intensity of the image I ( x , y ) at point y , ϕ ( y ) is the level set function at a point y , and δ ( ) and H ( ) are Equations (4) and (5).
The energy function of the fitting term E F F , based on adaptive fractional differentiation is,
E F F = γ [ Ω K σ ( x y ) ( F I ( y ) b 1 ( x ) ) 2 H ( ϕ ( y ) ) d x d y + Ω K σ ( x y ) ( F I ( y ) b 2 ( x ) ) 2 ( 1 H ( ϕ ( y ) ) ) d x d y ]
where b 1 ( x ) and b 2 ( x ) are the local area fitting means of the given pixel x inside and outside the contour curve, respectively,
b i ( x ) = K σ ( x ) [ H i ( ϕ ( x ) ) F I ( x ) ] K σ ( x ) H i ( ϕ ( x ) )         i = 1 , 2 .
Among them, F I ( x , y ) is the image I ( x , y ) after adaptive fractional differential processing. The specific derivation can be obtained from Equation (17) and is described in detail in Section 3.2,
F I ( x , y ) = I ( x , y ) + | D P G L I ( x , y ) | .
The gradient descent flow equation of the level set function ϕ obtained by minimizing the energy functional (25) is as follows,
ϕ t = λ δ ( ϕ ) [ K σ ( x y ) ( I ( y ) f 1 ( x ) ) 2 d x K σ ( x y ) ( I ( y ) f 2 ( x ) ) 2 d x ] γ δ ( ϕ ) [ K σ ( x y ) ( F I ( y ) b 1 ( x ) ) 2 d x K σ ( x y ) ( F I ( y ) b 2 ( x ) ) 2 d x ]
Finally, the gradient descent flow equation by minimizing the level set function ϕ of the total energy functional (19) in this paper is as follows,
ϕ t = μ div ( d R | ϕ | ϕ ) + ν g δ ( ϕ ) div ( ϕ | ϕ | ) λ δ ( ϕ ) [ K σ ( x y ) ( I ( y ) f 1 ( x ) ) 2 d x K σ ( x y ) ( I ( y ) f 2 ( x ) ) 2 d x ] γ δ ( ϕ ) [ K σ ( x y ) ( F I ( y ) b 1 ( x ) ) 2 d x K σ ( x y ) ( F I ( y ) b 2 ( x ) ) 2 d x ]
Then the level set evolution equation in Equation (30) can be discretized as,
ϕ i , j n + 1 ϕ i , j n Δ t = L ( ϕ i , j n ) ,
where ( i , j ) is the spatial index, k is the time index, and L ( ϕ i , j n ) is the numerical approximation of the function on the right side of the evolution Equation (31).
The iterative equation of the level set function can be expressed as,
ϕ i , j n + 1 = ϕ i , j n + Δ t L ( ϕ i , j n ) .

4.4. Implementation and Algorithm

The implementation and algorithm steps of the active contour model of image segmentation, based on adaptive fractional differentiation proposed in this paper, are as follows:
Step1: Input the image I ( x , y ) , initialize the level set function, and set the initial parameter values: μ = 0.2 , ν = 1.5 , λ = γ = 1 , Δ t = 5 , σ = 3 , ε = 1.5 .
Step2: Calculate the result D P G L I ( x , y ) after the adaptive fractional differential operation from the Equation (15) in Section 3.2, and update the image I ( x , y ) according to the Equation (17).
Step3: Calculate the fractional-order edge-stopping function g ( x , y ) from the Equation (18) in Section 3.3.
Step4: Calculate the penalty term energy function E P from the Equation (20) in Section 4.1.
Step5: Calculate the energy function of the regular term E R by the Equation (23) in Section 4.2.
Step6: Calculate the energy function of the fitting term E F F based on the adaptive fractional differential by the Equation (26) in Section 4.3, and then calculate the energy function of the total fitting term E F by the Equation (25).
Step7: Calculate the gradient descent flow equation ϕ / t of the level set function ϕ of the total energy function E by the Equation (30) in Section 4.3.
Step8: Calculate from the Equation (32) in Section 4.3, and update the level set function ϕ n + 1 ( x ) and ϕ n ( x ) .
Step9: Check whether the set number of iterations is reached. If the number of iterations is reached, the curve stops evolving to output the final target contour; otherwise, return to Step 2.
Algorithm The algorithm for solving the proposed model
Input:   Infrared   image   I ( x , y ) n × m .
Output: Image segmentation results.
1: Set initialization parameters: μ , ν , λ , γ , Δ t , σ , ε , Iter max ;
2: Calculate   D P G L I ( x , y ) by Equation (16);
3: for n = 1 : I t e r max do
4: Calculate   the   penalty   term   energy   function   E P by Equation (20);
5: Calculate   the   regularization   term   energy   function   E R by Equation (23);
6: Calculate the fitting term energy function E F by Equation (25);
7:   Calculate   the   gradient   descent   flow   equation   ϕ / t of the total energy function E by Equation (30);
8: Update the level set function ϕ n + 1 ( x )   and   ϕ n ( x ) by Equation (32);
9: if n = = Iter max then
10: break;
11: end if
12: end for

5. Experimental Results

In this section, the experimental contents are the robustness verification to noise, the robustness verification of the initial contour position, the comparison between fixed-order and adaptive differentiation, and the comparison with other methods (CV, RSF, Ref. [39], Ref. [40]). The computer configuration used in the experiment is an Intel CORE i5 processor, 8 GB memory, Windows 10 operating system, and the experimental platform MATLAB (R2016a).
In this paper, dice similarity coefficient (DSC) and intersection over union (IOU) are selected as performance evaluation indicators,
DSC ( A , B ) = 2 | R A R B | | R A | + | R B | ,
IOU ( A , B ) = | R A R B | | R A R B | .
The segmentation results were quantitatively compared with the expert manual segmentation results, and the area-based DSC index and the IOU overlap performance index were used for comparison. Among them, | | is the number of pixels, R A and R B are the expert segmentation area and the experimental segmentation result area respectively. The closer the DSC and IOU values are to 1, the higher the segmentation accuracy.

5.1. Verify the Robustness of Noise

To verify that the active contour image segmentation method based on adaptive fractional differentiation proposed in this paper has good noise resistance, this section will conduct verification experiments on different types of noisy images.
In Figure 3, the infrared thermal image of candle burning is taken as the experimental object. From top to bottom are the original image, the noise image with Gaussian noise, speckle noise, and mixed noise. The mixed noise includes Gaussian noise and speckle noise, and the noise density is 0.03. From left to right are the noisy image, CV, RSF, Ref. [39], Ref. [40], and our method. The red rectangle in Figure 3a represents the initial contour, and the green curve in (b–f) represents the evolution result of the contour curve. The number of iterations of this experiment is I t e r max = 20 .
As shown in Figure 3b–e, the segmentation result is affected by noise to varying degrees, and some noise points are mistakenly considered as edges of the image. From the segmentation results in Figure 3f, it can be seen that the method proposed in this paper obtain better segmentation results under three types of noise conditions. These results show that the Gaussian kernel in the local fitting term is effective in dealing with the noise of the evolution process.
In addition to visual effects, IOU and DSC performance index values under different models in Table 1 were also compared. It can be seen that our model can obtain better performance indicators and segmentation results.

5.2. Initial Contour Position Robustness Verification

In this section, test experiments are used to verify the robustness of the proposed model to the initial curve, and the infrared image of the aircraft is taken as the experimental object. Figure 4 shows the segmentation results of different models and the proposed method for non-uniform images with different initial shapes. In Figure 4, from left to right are the input image, CV, RSF, Ref. [39], Ref. [40], and the proposed method where the input image has initial contours at different locations. In addition, the evolution process diagram of the level set function of different models is shown inFigure 5, which verifies the stability of the level set function of the model in this paper during the evolution process. The red rectangle in Figure 4a represents the initial contour, and the green curve in (b–f) represents the evolution result of the contour curve. The number of iterations of this experiment is I t e r max = 40 .
The results in Figure 4c–e show that the three models are sensitive to the initial contour position and often fail to achieve ideal segmentation results due to improper initial conditions. In Figure 4b–e, the CV model and the method in this paper can obtain similar results under different initial conditions. It is verified that our model is robust to the position of the initialization contour. In Figure 5a–d, the level set functions of other models will become more and more irregular in the iteration process, leading to oscillation in the evolution process, and then the accuracy of calculation will be seriously affected. Finally, the stability of level set evolution will be destroyed. In Figure 5e under the action of the double-well potential function, the final level set function is relatively smoother, which is conducive to the stability of the model.

5.3. Comparison between Fixed Order and Adaptive Differential

In this subsection, we will segment images with blurred object edges, and verify the performance of combining local region fitting and adaptive fractional differentiation in image segmentation. Furthermore, we validate the constraints on the correction parameters in adaptive fractional differentiation.
In the segmentation results in Figure 6, there are five different infrared images Pic 1–Pic 5 from top to bottom. The image contents are aircraft, human arms, and daily living environment. From left to right are the input image and initial contour, ordinary integer-order differential, 0.8-order fixed-order fractional differential, and our method. The red rectangle in Figure 6a represents the initial contour, and the green curve in (b–d) represents the evolution result of the contour curve. The number of iterations of this experiment is I t e r max = 40 .
The above is the segmentation result and performance index after processing the infrared grayscale image. Figure 6b shows that for infrared images with large background noise and blurred visual effects, ordinary integer order differential is easy to fall into a local minimum. In Figure 6c, a 0.8-order fixed-order fractional order G-L differential is adopted, and the same order differential processing is performed in the low-frequency part and the high-frequency part of the image, resulting in smooth texture details of the image. Therefore, compared with ordinary integer order differentials, the segmentation result and performance index are poor. In Figure 6d, after adding the adaptive fractional local region fitting term, we observe it has certain robustness to noise, can effectively separate the background and the target, and is more noise resistant. At the same time, due to the adaptive fractional differentiation for different parts of texture details, more target details of the image can be extracted, and a better image segmentation effect is achieved. Table 2 also shows the effectiveness of this method in terms of evaluation index values.
In addition, [32,33,34] solved the optimal control problem of nonlinear fractional order systems through gradient-based optimization methods. In this paper, the nonlinear parameter optimization strategy is used to process the modified parameters α and β , taking the performance indicators IOU and DSC as the objective function, and the modified parameter α , β as the input variables, taking random values in [ 0 , 1 ] , and finally obtaining the optimal values of IOU and DSC. In Figure 7, the corresponding IOU and DSC values of the correction parameters under random values are shown respectively.
Figure 7 shows that for infrared images, when α > β , the image segmentation performance IOU and DSC values are better than α < β . Therefore, in the order matrix P with adaptive order change, the constraint condition α > β of its modified parameters is established; in combination with the highest points of Figure 7a,b, the optimal values of α and β are 1 and 0.8 respectively. Next, this paper selects those α and β values of the best and non-best in Figure 7, processes different infrared images, and compares their processing effects in the edge part. In Figure 8, from left to right, there are seven different infrared images, including human body movements, flames when candles burn, aircraft, and daily living environment. From top to bottom are the input image, the nonoptimal segmentation result, and the optimal segmentation result. The red rectangle in the first row of images represents the initial contour, and the green curve in the second and third rows of images represents the evolution result of the contour curve.
The results in Figure 8 show that, for infrared images, the best correction parameters enable the method in this paper to obtain finer edge contours and show better results.

5.4. Compare with Other Methods

In this subsection, we compare the proposed model with CV, RSF, Ref. [39], and Ref. [40]. The effectiveness of the proposed model in processing noisy and uneven grayscale infrared images is verified. In Figure 9, there are five different infrared images Pic6–Pic10 from top to bottom. The image content is the infrared thermal image of burning candles, human palms, and animals. From left to right are the input image and the segmentation results of the initial contour, CV, RSF, Ref. [39], Ref. [40], and the proposed model in this paper, respectively. The red rectangle in Figure 9a represents the initial contour, and the green curve in (b–d) represents the evolution result of the contour curve. The number of iterations of this experiment is I t e r max = 40 .
The above is the segmentation result and performance index after processing the infrared thermal image. Figure 9b–e show that for infrared images with large background noise and blurred visual effect, the four models are easily affected by noise and initial contour position, and finally fall into the local minimum value. In Figure 9f, the method proposed in this paper can effectively separate the background from the target, has more noise resistance, extracts more target details of the image, and achieves a better image segmentation effect. It is shown in Table 3 that when fractional variable order differential processing is adopted, our model can obtain finer edge contours and better handle some details compared with other models. At the same time, it can also reflect the effectiveness of this model in evaluating the index value.

6. Conclusions

In this paper, an adaptive fractional differential active contour image segmentation method is proposed. Adaptive fractional differentiation is adopted to extract more textured parts from the image by adaptive order differentiation for each pixel. The fractional edge-stopping function improves noise resistance. Additionally, the fitting term based on adaptive fractional differentiation solves the problem that improper selection of the initial contour position will lead to inaccurate segmentation results. The experimental results show that for the infrared image with strong noise and blurred edges, this method can solve the problems of being sensitive to noise and initial contour position, and effectively segment the image of the weak edge target. In future research work, we will explore adaptive weighting techniques to achieve faster calculations and use fewer parameters and construct a better-performing active contour model for image segmentation.

Author Contributions

Conceptualization, Y.Z.; methodology, Y.Z., L.Y. and Y.L.; software, Y.Z. and L.Y.; validation, Y.Z. and L.Y.; formal analysis, Y.Z., L.Y. and Y.L.; investigation, Y.Z. and L.Y.; resources, Y.Z.; data curation, Y.Z., L.Y. and Y.L.; writing—original draft preparation, Y.Z., L.Y. and Y.L.; writing—review and editing, Y.Z., L.Y. and Y.L.; visualization, Y.Z., L.Y. and Y.L.; supervision, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Liaoning Provincial Education Department Scientific Research Project (LJKZ0245) and National Key Laboratory Project (2021JCJQLB055006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data were computed using our algorithm.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kass, M.; Witkin, A.; Terzopoulos, D. Snakes: Active contour models. Int. J. Comput. Vis. 1988, 1, 321–331. [Google Scholar] [CrossRef]
  2. Lv, H.; Zhang, Y.; Wang, R. Active contour model based on local absolute difference energy and fractional-order penalty term. J. Appl. Math. Model. 2022, 107, 207–232. [Google Scholar] [CrossRef]
  3. Jin, R.; Weng, G. Active contours driven by adaptive functions and fuzzy c-means energy for fast image segmentation. J. Signal Process. 2019, 163, 1–10. [Google Scholar] [CrossRef]
  4. Lv, H.; Zhang, F.; Wang, R. Robust active contour model using patch-based signed pressure force and optimized fractional-order edge. J. IEEE Access 2021, 9, 8771–8785. [Google Scholar] [CrossRef]
  5. Wei, X.; Zheng, W.; Xi, C.; Shang, S. Shoreline extraction in SAR image based on advanced geometric active contour model. J. Remote Sens. 2021, 13, 642. [Google Scholar] [CrossRef]
  6. Karami, E.; Shehata, M.S.; Smith, A. Adaptive polar active contour for segmentation and tracking in ultrasound videos. J. IEEE Trans. Circuits Syst. Video Technol. 2018, 29, 1209–1222. [Google Scholar] [CrossRef]
  7. Hu, P.; Wang, G.; Kong, X.; Kuen, J.; Tan, Y.-P. Motion-guided cascaded refinement network for video object segmentation. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Salt Lake City, UT, USA, 18–22 June 2018; pp. 1400–1409. [Google Scholar]
  8. Zhang, P.; Liu, W.; Lei, Y.; Wang, H.; Lu, H. Deep multiphase level set for scene parsing. J. IEEE Trans. Image Process. 2020, 29, 4556–4567. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, S.; Wu, J.; Feng, Z. Dual-Fusion active contour model with semantic information for saliency target extraction of underwater images. J. Appl. Sci. 2022, 12, 2515. [Google Scholar] [CrossRef]
  10. Pan, Z.; Liu, S.; Fu, W. A review of visual moving target tracking. J. Multimed. Tools Appl. 2017, 76, 16989–17018. [Google Scholar] [CrossRef]
  11. Lobregt, S.; Viergever, M.A. A discrete dynamic contour model. J. IEEE Trans. Med. Imaging 1995, 14, 12–24. [Google Scholar] [CrossRef]
  12. Bendaoud, A.; Hachouf, F. Fractional calculus for improving edge-based active contour models. In Proceedings of the 2021 International Conference on Recent Advances in Mathematics and Informatics (ICRAMI), Tebessa, Algeria, 21–22 September 2021; IEEE: Piscataway, NJ, USA, 2021; pp. 1–5. [Google Scholar]
  13. Liu, C.; Liu, W.; Xing, W. An improved edge-based level set method combining local regional fitting information for noisy image segmentation. J. Signal Process. 2017, 130, 12–21. [Google Scholar] [CrossRef]
  14. Yu, H.; He, F.; Pan, Y. A novel segmentation model for medical images with intensity inhomogeneity based on adaptive perturbation. J. Multimed. Tools Appl. 2019, 78, 11779–11798. [Google Scholar] [CrossRef]
  15. Liu, C.; Liu, W.; Xing, W. A weighted edge-based level set method based on multi-local statistical information for noisy image segmentation. J. Vis. Commun. Image Represent. 2019, 59, 89–107. [Google Scholar] [CrossRef]
  16. Yang, Y.; Hou, X.; Ren, H. Efficient active contour model for medical image segmentation and correction based on edge and region information. J. Expert Syst. Appl. 2022, 194, 116436. [Google Scholar] [CrossRef]
  17. Li, M.; Zhan, Y.; Ge, Y. Fractional distance regularized level set evolution with its application to image segmentation. J. IEEE Access 2020, 8, 84604–84617. [Google Scholar] [CrossRef]
  18. Chen, B.; Huang, S.; Liang, Z.; Chen, W.; Pan, B. A fractional order derivative based active contour model for inhomogeneous image segmentation. J. Appl. Math. Model. 2019, 65, 120–136. [Google Scholar] [CrossRef]
  19. Yan, X.; Weng, G. Hybrid active contour model driven by optimized local pre-fitting image energy for fast image segmentation. J. Appl. Math. Model. 2022, 101, 586–599. [Google Scholar] [CrossRef]
  20. Fang, J.; Liu, H.; Zhang, L.; Liu, J.; Liu, H. Region-edge-based active contours driven by hybrid and local fuzzy region-based energy for image segmentation. J. Inf. Sci. 2021, 546, 397–419. [Google Scholar] [CrossRef]
  21. Caselles, V.; Kimmel, R.; Sapiro, G. Geodesic active contours. Int. J. Comput. Vis. 1997, 22, 61–79. [Google Scholar] [CrossRef]
  22. Li, C.; Xu, C.; Gui, C.; Fox, M.D. Level set evolution without re-initialization: A new variational formulation. In Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR’05), San Diego, CA, USA, 20–25 June 2005; IEEE: Piscataway, NJ, USA, 2005; pp. 430–436. [Google Scholar]
  23. Chan, T.F.; Vese, L.A. Active contours without edges. J. IEEE Trans. Image Process. 2001, 10, 266–277. [Google Scholar] [CrossRef] [PubMed]
  24. Li, C.; Kao, C.-Y.; Gore, J.C.; Ding, Z. Minimization of region-scalable fitting energy for image segmentation. J. IEEE Trans. Image Process. 2008, 17, 1940–1949. [Google Scholar]
  25. Zhang, X.; Liu, R.; Ren, J.; Gui, Q. Adaptive fractional image enhancement algorithm based on rough set and particle swarm optimization. J. Fractal Fract. 2022, 6, 100. [Google Scholar] [CrossRef]
  26. Yan, H.; Zhang, X. Adaptive fractional multi-scale edge-preserving decomposition and saliency detection fusion algorithm. J. ISA Trans. 2020, 107, 160–172. [Google Scholar] [CrossRef]
  27. Zhang, X.; Dai, L. Image enhancement based on rough set and fractional order differentiator. J. Fractal Fract. 2022, 6, 214. [Google Scholar] [CrossRef]
  28. Yan, H.; Zhang, J.-X.; Zhang, X. Injected infrared and visible image fusion via L1 decomposition model and guided filtering. J. IEEE Trans. Comput. Imaging 2022, 8, 162–173. [Google Scholar] [CrossRef]
  29. Khanna, S.; Chandrasekaran, V. Fractional differentiation-based edge energy driven active contours for robust image segmentation. ICTACT J. Image Video Process. 2021, 12, 2493–2501. [Google Scholar]
  30. Zhang, Y.; Tian, Y. A new image segmentation method based on fractional-varying-order differential. J. Beijing Inst. Technol. 2021, 30, 254–264. [Google Scholar]
  31. Li, M.-M.; Li, B.-Z. A novel active contour model for noisy image segmentation based on adaptive fractional order differentiation. J. IEEE Trans. Image Process. 2020, 29, 9520–9531. [Google Scholar] [CrossRef]
  32. Liu, C.; Gong, Z.; Yu, C.; Wang, S.; Teo, K.L. Optimal control computation for nonlinear fractional time-delay systems with state inequality constraints. J. Optim. Theory Appl. 2021, 191, 83–117. [Google Scholar] [CrossRef]
  33. Gong, Z.; Liu, C.; Teo, K.L.; Yi, X. Optimal control of nonlinear fractional systems with multiple pantograph-delays. J. Appl. Math. Comput. 2022, 425, 127094. [Google Scholar] [CrossRef]
  34. Gong, Z.; Liu, C.; Teo, K.L.; Wang, S.; Wu, Y. Numerical solution of free final time fractional optimal control problems. J. Appl. Math. Comput. 2021, 405, 126270. [Google Scholar] [CrossRef]
  35. Zhang, X.; Yan, H. Medical image fusion and noise suppression with fractional-order total variation and multi-scale decomposition. J. IET Image Process. 2021, 15, 1688–1701. [Google Scholar] [CrossRef]
  36. Zhang, X.; He, H.; Zhang, J.-X. Multi-focus image fusion based on fractional order differentiation and closed image matting. J. ISA Trans. 2022, in press. [Google Scholar] [CrossRef] [PubMed]
  37. Li, C.; Xu, C.; Gui, C.; Fox, M.D. Distance regularized level set evolution and its application to image segmentation. J. IEEE Trans. Image Process. 2010, 19, 3243–3254. [Google Scholar] [CrossRef]
  38. Aubert, G.; Kornprobst, P.; Aubert, G. Mathematical Problems in Image Processing: Partial Differential Equations and the Calculus of Variations; Springer: Berlin/Heidelberg, Germany, 2006; p. 147. [Google Scholar]
  39. Ding, K.; Xiao, L.; Weng, G. Active contours driven by local pre-fitting energy for fast image segmentation. J. Pattern Recognit. Lett. 2018, 104, 29–36. [Google Scholar] [CrossRef]
  40. Liu, H.; Fang, J.; Zhang, Z.; Lin, Y. A novel active contour model guided by global and local signed energy-based pressure force. J. IEEE Access 2020, 8, 59412–59426. [Google Scholar] [CrossRef]
Figure 1. Example of adaptive order matrix P with dimension 5 × 5 .
Figure 1. Example of adaptive order matrix P with dimension 5 × 5 .
Fractalfract 06 00579 g001
Figure 2. Level set function. (a) Level set function in [22]; (b) level set function in [37].
Figure 2. Level set function. (a) Level set function in [22]; (b) level set function in [37].
Fractalfract 06 00579 g002
Figure 3. Segmentation results for different types of noise. (a) Input image and initial contour; (b) CV; (c) RSF; (d) Ref. [39]; (e) Ref. [40]; (f) Ours.
Figure 3. Segmentation results for different types of noise. (a) Input image and initial contour; (b) CV; (c) RSF; (d) Ref. [39]; (e) Ref. [40]; (f) Ours.
Fractalfract 06 00579 g003
Figure 4. Segmentation results in different initial contours. (a) Input image and initial contour; (b) CV; (c) RSF; (d) Ref. [39]; (e) Ref. [40]; (f) Ours.
Figure 4. Segmentation results in different initial contours. (a) Input image and initial contour; (b) CV; (c) RSF; (d) Ref. [39]; (e) Ref. [40]; (f) Ours.
Fractalfract 06 00579 g004
Figure 5. Level set functions for different models. (a) CV; (b) RSF; (c) Ref. [39]; (d) Ref. [40]; (e) Ours.
Figure 5. Level set functions for different models. (a) CV; (b) RSF; (c) Ref. [39]; (d) Ref. [40]; (e) Ours.
Fractalfract 06 00579 g005
Figure 6. Segmentation results in a fixed order and variable order. (a) Input image Pic 1–Pic 5 and initial contour; (b) General derivative; (c) Fixed order; (d) Adaptive order.
Figure 6. Segmentation results in a fixed order and variable order. (a) Input image Pic 1–Pic 5 and initial contour; (b) General derivative; (c) Fixed order; (d) Adaptive order.
Fractalfract 06 00579 g006
Figure 7. Performance index value. (a) IOU value; (b) DSC value.
Figure 7. Performance index value. (a) IOU value; (b) DSC value.
Fractalfract 06 00579 g007
Figure 8. Comparison of infrared image segmentation results between the best and non-best. (a) Infrared image of human motion in simple background; (b) Infrared image of human motion in complex background; (c) Infrared image of candle burning flame; (d) Infrared image of aircraft; (e) Infrared image of fighter; (f) Infrared image under complex background; (g) Infrared image under complex background.
Figure 8. Comparison of infrared image segmentation results between the best and non-best. (a) Infrared image of human motion in simple background; (b) Infrared image of human motion in complex background; (c) Infrared image of candle burning flame; (d) Infrared image of aircraft; (e) Infrared image of fighter; (f) Infrared image under complex background; (g) Infrared image under complex background.
Fractalfract 06 00579 g008
Figure 9. Segmentation results of infrared images. (a) Input image Pic6–Pic10 and initial contour; (b) CV; (c) RSF; (d) Ref. [39]; (e) Ref. [40]; (f) Ours.
Figure 9. Segmentation results of infrared images. (a) Input image Pic6–Pic10 and initial contour; (b) CV; (c) RSF; (d) Ref. [39]; (e) Ref. [40]; (f) Ours.
Fractalfract 06 00579 g009
Table 1. Corresponding evaluation index values.
Table 1. Corresponding evaluation index values.
Noise ImagePerformance IndexModels
CVRSFRef. [39]Ref. [40]Ours
ClearIOU0.91460.90330.81530.71040.9856
DSC0.93680.94920.89830.83070.9927
GaussianIOU0.91830.86780.75400.81090.9814
DSC0.93350.92920.85980.89550.9898
SpeckleIOU0.92760.76960.80580.71250.9714
DSC0.94310.86980.89250.83210.9898
Gaussian + SpeckleIOU0.92220.78130.74280.73100.9687
DSC0.94020.87720.85240.84460.9930
Table 2. Corresponding evaluation index values are in Figure 6.
Table 2. Corresponding evaluation index values are in Figure 6.
ImagePerformance IndexModels
General DerivativeFixed OrderAdaptive Order
Pic1IOU0.98810.97980.9896
DSC0.99400.98980.9948
Pic2IOU0.97920.97220.9846
DSC0.98950.98590.9923
Pic3IOU0.84650.82000.9760
DSC0.91690.90010.9878
Pic4IOU0.95320.96800.9748
DSC0.97600.98380.9900
Pic5IOU0.96880.78470.9760
DSC0.98420.87940.9878
Table 3. Performance index values corresponding to Figure 9.
Table 3. Performance index values corresponding to Figure 9.
ImagePerformance IndexModels
CVRSFRef. [39]Ref. [40]Ours
Pic6IOU0.87810.90590.72930.57590.9832
DSC0.93510.95060.84350.73090.9915
Pic7IOU0.62450.51390.44820.37170.7101
DSC0.76890.67890.61900.54200.8579
Pic8IOU0.93310.65890.64100.58090.9510
DSC0.96540.79440.78120.73490.9749
Pic9IOU0.81250.86470.83090.82560.9929
DSC0.89660.92750.90770.90450.9965
Pic10IOU0.75270.69880.71350.66660.9539
DSC0.85890.82270.83280.79990.9764
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Yang, L.; Li, Y. A Novel Adaptive Fractional Differential Active Contour Image Segmentation Method. Fractal Fract. 2022, 6, 579. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6100579

AMA Style

Zhang Y, Yang L, Li Y. A Novel Adaptive Fractional Differential Active Contour Image Segmentation Method. Fractal and Fractional. 2022; 6(10):579. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6100579

Chicago/Turabian Style

Zhang, Yanzhu, Lijun Yang, and Yan Li. 2022. "A Novel Adaptive Fractional Differential Active Contour Image Segmentation Method" Fractal and Fractional 6, no. 10: 579. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6100579

Article Metrics

Back to TopTop