A Novel Gradient Vector Flow Snake Model Based on Convex Function for Infrared Image Segmentation

Infrared image segmentation is a challenging topic because infrared images are characterized by high noise, low contrast, and weak edges. Active contour models, especially gradient vector flow, have several advantages in terms of infrared image segmentation. However, the GVF (Gradient Vector Flow) model also has some drawbacks including a dilemma between noise smoothing and weak edge protection, which decrease the effect of infrared image segmentation significantly. In order to solve this problem, we propose a novel generalized gradient vector flow snakes model combining GGVF (Generic Gradient Vector Flow) and NBGVF (Normally Biased Gradient Vector Flow) models. We also adopt a new type of coefficients setting in the form of convex function to improve the ability of protecting weak edges while smoothing noises. Experimental results and comparisons against other methods indicate that our proposed snakes model owns better ability in terms of infrared image segmentation than other snakes models.


Introduction
Infrared radiation is an invisible type of electromagnetic wave, whose wavelength ranges between the radio wave and the visible light. Any object in nature whose temperature is over absolute zero (−273 • C) is able to radiate the infrared ray. Compared with the visible light, the infrared wave has its unique characteristics. For example, compared with the visible light, its light quantum energy is much lower, the heat effect is stronger, it is more likely to be absorbed by a substance, and is less sensitive to the human eye.
The visible light between 0.4~0.75 µm can be sensed by human eyes. The light outside this range cannot be sensed without the aid of detectors. After in-depth studies, the interaction between the medium and the radiation source has been discovered and the infrared radiation laws have been summarized, which greatly stimulate the development of the infrared technology. The advent and development of the infrared thermal imaging system indirectly broadens the visual sensing scope of the eyes. The most commonly used detector is the infrared thermal detector. It can measure the infrared thermal radiation quantity in a non-contact manner, convert it into clearly visible images and display them on a screen.
The infrared technology was first applied for military purposes. In recent years, the technology has been widely used in transportation, medicine and other scientific areas. The on-board infrared scanning imaging technology can be used to monitor the location and affected regions of a forest fire, control a fire disaster and minimize losses. In medicine, the infrared technology can be used to detect inflamed organs and diagnose early symptoms of cancers. In the electronic equipment manufacturing

Traditional Snakes Model
In 1987, Kass and co-workers proposed an active contour model, which is also known as the Snakes model [17]. The traditional version of the active contour model is a continuous closed curve and represented with parameter curve c (q) = [x (q) , y (q)] , q ∈ [0, 1]. The energy functional is minimized by moving the curve in the image, as shown in Equation (1).
where α and β are the weighting coefficients that adjust the flexibility and rigidity of the curve in the active contour. The first integral term in the equation is the internal energy that ensures smoothness and continuity of the curve. The second integral term is the exterior energy, which contains the information of the image where the contour curve is located. It is a man-made constraint specifically introduced to guarantee the curve evolves towards the object contour more accurately and quickly. In the traditional Snakes model, the external energy is usually defined as the local features of the image where the controlling point or the connecting line is located. The gradient is usually used as the feature, as shown in Equation (2).
The limitations of the traditional Snakes model are as follows: 1 It is very sensitive to the location of the initial contour. During the practical segmentation process, the initial location of the contour must be manually put near the image edge of interest, resulting in poor interactivity. 2 It is prone to converge towards the false edge near the object, and is not robust to the noise. 3 Its convergence performance is poor for the object contour with sunken regions.

The GVF Snakes Model
The catching range of the traditional Snakes model is limited, and the external energy only exists in the regions near the object contour. Xu and Prince proposed a new external force model for the active contour, i.e., the gradient vector flow (GVF) Snakes model [18]. GVF refers to the vector field obtained by propagating the gradient vector of the edge graph for the given image. It is represented with a function and can be determined using the following dynamic evolution equation.
where µ is the parameter that can control the degree of smoothness of the external force field in GVF, and the value set to it increases with the noise intensity in the image. f is the edge graph of the input image. ∇ 2 is the Laplace operator. The largest advantage of GVF Snakes over the traditional Snakes model is its ability to broaden the catching range of the initial contour and to catch the high-curvature region in the object contour.

An Improved GVF Snakes Model
Although the GVF Snakes model has many advantages, it also has many limitations. Many variants have been presented to address the limitations.

GGVF Snakes Model
In 1998, Xu and co-workers introduced two weighting coefficients that can change in the image domain to the iteration equation of the GVF external force field. In this way, they obtained a new external force called the generic gradient vector flow (GGVF) external force field [19]. The evolution equation of this external force field is: where the parameter K determines the weight of the smooth term and the data term. As in the GVF Snakes model, the choice of K relates to the image noise. The larger the noise, the larger the value of K should be. This model provides an approach to the problem that GVF Snakes can hardly converge towards the long narrow sunken regions and is not very robust to the noise.
Afterwards, Qin proposed an improved model of the GGVF external force field in 2013, i.e., the component-based normalized GGVF model (CNGGVF) external force field [6]. CNGGVF addresses the problem of GGVF Snakes that it can hardly converge towards LTI (Long and Thin Indentation), which is an even number of pixels in width.

NGVF Snakes and NBGVF Snakes
The Laplace operator can be decomposed along the tangent and normal directions. Hence, the evolution equation of the GVF external force field can be rewritten as: where V TT and V NN denote the second-order derivative along the tangent and normal directions of V t . The parameters α and β determine the degree of image diffusion along the tangent and normal directions. Normally, the interpolation method yields the best results. Based on this, Ning et al. proposed a normal gradient vector flow (NGVF) external force field [20]. From Equation (7), we can know that: After their investigations, You et al. discovered that diffusion along the tangent direction of the image edge can protect the image edge, and the diffusion along the normal direction can smooth the noise. The NGVF external force field abandons the tangent diffusion, making it difficult for the NGVF Snakes model to protect the weak edge of images. In this context, Wang et al. proposed a normally biased gradient vector flow (NBGVF) external force field [21]. NBGVF completely retains the tangent diffusion and is capable of adapting normal diffusion to image structure.
To sum up, in the NBGVF Snakes model, the parameters and definitions in Equation (7) can be defined as: This improved version of the model has higher diffusion efficiency and is able to effectively protect the weak edges. The weaker the edge to be protected, the smaller the value of K.

Improved Version of the GVF Model
As discussed in the section above, the GGVF Snakes model enlarges the convergence range of the active contour, improves the LTI convergence performance and is more robust to the noise. Based on NGVF, which has higher diffusion efficiency, NBGVF provides a solution to the weak edge protection problem. Hence, this paper relies on the GVF external force model, and combines GGVF and NBGVF to propose a novel external force model.
The improved version of the external force is defined as a vector field, and it can be obtained by using the following energy functional: where V NN and V TT denote the second-order derivative along the normal and tangent directions. g (|∇ f |) and h (|∇ f |) denote the coefficients of the smooth and data terms in Equation (12). As defined in GGVF, the value of K increases with the noise intensity in the image, but this may lead to the weak edge being over-smoothed. Unlike the coefficients of the normal and tangent diffusion operators in NBGVF, both of the coefficients directly depend on the intensity rather than the gradient of the edge graph, thereby greatly reducing the computational complexity. Moreover, as shown in Figure 1, the variation of the parameters in Equations (15) and (16) with the intensity of the edge graph takes the form of convex function. Compared with the parameters α and β in NBGVF (Equations (10) and (11)), the coefficients of the proposed model change gradually when the value of f is high, and thus offer more protection to the weak edge in the infrared images. Hence, the proposed model is capable of segmenting the infrared images more accurately. Meanwhile, the coefficients of the proposed model fluctuate violently when the value of f is low. As a result, contour divergence is more efficient at a long distance from the edge. where and denote the second-order derivative along the normal and tangent directions. (|∇ |) and ℎ(|∇ |) denote the coefficients of the smooth and data terms in Equation (12). As defined in GGVF, the value of K increases with the noise intensity in the image, but this may lead to the weak edge being over-smoothed. Unlike the coefficients of the normal and tangent diffusion operators in NBGVF, both of the coefficients directly depend on the intensity rather than the gradient of the edge graph, thereby greatly reducing the computational complexity.
Moreover, as shown in Figure 1, the variation of the parameters in Equations (15) and (16) with the intensity of the edge graph takes the form of convex function. Compared with the parameters and in NBGVF (Equations (10) and (11)), the coefficients of the proposed model change gradually when the value of f is high, and thus offer more protection to the weak edge in the infrared images. Hence, the proposed model is capable of segmenting the infrared images more accurately. Meanwhile, the coefficients of the proposed model fluctuate violently when the value of f is low. As a result, contour divergence is more efficient at a long distance from the edge. After the parameter in the equation is set to 0.1, the variation of the two coefficients is shown in Figure 2. This figure shows that the two coefficients fluctuate violently, dwindle to zero and then jump to 1 when the intensity increases. This means that the diffusion of the image near the edge graph along the normal direction is inhibited quickly, and that only the tangent diffusion component is left finally. This type of variation is conducive to the protection of weak edges. As discussed above, the value of K should increases with the image noise to smooth the noise, but the weak edge is likely to be lost due to noise smoothing. Hence, the value of the parameter should be optimized to achieve a trade-off between noise smoothing and weak edge protection. After the parameter in the equation is set to 0.1, the variation of the two coefficients is shown in Figure 2. This figure shows that the two coefficients fluctuate violently, dwindle to zero and then jump to 1 when the intensity increases. This means that the diffusion of the image near the edge graph along the normal direction is inhibited quickly, and that only the tangent diffusion component is left finally. This type of variation is conducive to the protection of weak edges. As discussed above, the value of K should increases with the image noise to smooth the noise, but the weak edge is likely to be lost due to noise smoothing. Hence, the value of the parameter should be optimized to achieve a trade-off between noise smoothing and weak edge protection.
graph along the normal direction is inhibited quickly, and that only the tangent diffusion component is left finally. This type of variation is conducive to the protection of weak edges. As discussed above, the value of K should increases with the image noise to smooth the noise, but the weak edge is likely to be lost due to noise smoothing. Hence, the value of the parameter should be optimized to achieve a trade-off between noise smoothing and weak edge protection.

Numerical Implementation
Now, the external force field can be obtained by minimizing Equation (12). The Euler-Lagrange equation of the energy functional can be written as: In order to obtain the vector field in Equation (17), we introduce the parameter t and construct the following partial differential equation.
where V x V y is the first-order partial derivative with respect to x or y, V xx V yy is the second-order partial derivative with respect to x or y, and V xy is the result achieved by computing the partial derivative with respect to x and then to y. The equations above can be solved by finding the equilibrium solution to the following set of partial differential equations: Iterating Equation (20) yields the desired external force field. Hence, the evolution equation of this external force field can be written as: The algorithm steps are given in Figure 3.
where = , = , = , and = . Iterating Equation (20) yields the desired external force field. Hence, the evolution equation of this external force field can be written as: The algorithm steps are given in Figure 3.

Experimental Results and Analysis
In this section, the proposed GVF model will be compared with GVF [18], GGVF [19], NGVF [20], NBGVF [21], CN-GGVF [6], LIF [28] and SOAC [29] across different images. First, we apply these methods to standard images, including the U-shaped image and the LTI image. All of these images are the traditional images used to evaluate the basic performance of various Snakes models. Afterwards, we will evaluate the performance of the proposed model and other algorithms in terms of segmenting infrared images, such as the original infrared image and the infrared images corrupted with various types of noises. These segmentation results form the basis for detailed analysis and comparison. MATLAB R2014B is used as the development environment of the experiment programs in this paper. The computer configuration is Inter Core i5-4210M 2.6 GHz CPU and 8 GB RAM.
Subjective assessment has its limitations for the evaluation of segmentation performance. In our experiments, the segmentation results are evaluated using the following metrics: Precision, Recall and F1 measure [1]. Let M seg denote the actual segmentation results and G seg denote the segmentation baseline.
The metric Prevision can be expressed as: Similarly, Recall can be defined as: F1 measure provides an evaluation metric that combines Precision with Recall. It is defined as: A high value for any of these three metrics means that the segmentation is accurate and the result approximates to the ground truth.

Catching Range, Convergence for Convex and Concave Planes and Insensitivity to Initial Contours
In this set of experiments, we use the U-shaped and square images to test the performance of the proposed method. The contour of the proposed method evolves from a long distance away towards the target edge of the image. The parameter setting of the proposed method is {K, τ} = {0.1, 1} and the evolution is shown in Figures 4 and 5. It can be seen that the final contour is well matched with the target edge. The results in Figure 4a show the large catching range of the proposed model. The results in Figure 4a,b demonstrate the ability of the proposed method to obtain accurate segmentation results regardless of where the initial contour is placed and whether the contour is distant from the object or passes through the target edge.
F1 measure provides an evaluation metric that combines Precision with Recall. It is defined as: A high value for any of these three metrics means that the segmentation is accurate and the result approximates to the ground truth.

Catching Range, Convergence for Convex and Concave Planes and Insensitivity to Initial Contours
In this set of experiments, we use the U-shaped and square images to test the performance of the proposed method. The contour of the proposed method evolves from a long distance away towards the target edge of the image. The parameter setting of the proposed method is , = 0.1, 1 and the evolution is shown in Figures 4 and 5. It can be seen that the final contour is well matched with the target edge. The results in Figure 4a show the large catching range of the proposed model. The results in Figure 4a,b demonstrate the ability of the proposed method to obtain accurate segmentation results regardless of where the initial contour is placed and whether the contour is distant from the object or passes through the target edge.

Convergence for Long Narrow Edges
In this subsection, the proposed method will be evaluated and compared with other traditional active contour models in terms of LTI image segmentation. The parameter setting of the proposed method is , = 1, 0.5 . The classic LTI images can test the convergence performance of the active contour model in the case of long narrow edges. The experimental results are given in Figure 6. These results clearly show that only the proposed method has the ability to converge towards the bottom of the LTI images, while the other algorithms stop converging at the entrance to the LTI images. Hence, the proposed algorithm has remarkable superiority over the traditional models.   Figure 5 demonstrates the ability of the proposed method to converge for convex and concave planes and obtain the U-shaped edge accurately through segmentation.

Convergence for Long Narrow Edges
In this subsection, the proposed method will be evaluated and compared with other traditional active contour models in terms of LTI image segmentation. The parameter setting of the proposed method is {K, τ} = {1, 0.5}. The classic LTI images can test the convergence performance of the active contour model in the case of long narrow edges. The experimental results are given in Figure 6. These results clearly show that only the proposed method has the ability to converge towards the bottom of the LTI images, while the other algorithms stop converging at the entrance to the LTI images. Hence, the proposed algorithm has remarkable superiority over the traditional models.
active contour models in terms of LTI image segmentation. The parameter setting of the proposed method is , = 1, 0.5 . The classic LTI images can test the convergence performance of the active contour model in the case of long narrow edges. The experimental results are given in Figure 6. These results clearly show that only the proposed method has the ability to converge towards the bottom of the LTI images, while the other algorithms stop converging at the entrance to the LTI images. Hence, the proposed algorithm has remarkable superiority over the traditional models.

Parameter Settings Sensitivity
To demonstrate the parameter sensitivity of our proposed model, we change the parameter "τ" from 0.01 to 1 and obtain the experiment results. Next, the results are quantitative calculated by the F1 measure criterion. Then we obtain the average F1 measure value of the experiment. Finally, the following curve in Figure 7 is acquired. According to Figure 7, in the range of 0.1 to 0.5 and 0.8 to 1, the F1 measure values change is not so obvious. Thus, our proposed model can be insensitive to parameter settings in the certain range.

Parameter Settings Sensitivity
To demonstrate the parameter sensitivity of our proposed model, we change the parameter "τ" from 0.01 to 1 and obtain the experiment results. Next, the results are quantitative calculated by the F1 measure criterion. Then we obtain the average F1 measure value of the experiment. Finally, the following curve in Figure 7 is acquired. According to Figure 7, in the range of 0.1 to 0.5 and 0.8 to 1, the F1 measure values change is not so obvious. Thus, our proposed model can be insensitive to parameter settings in the certain range.

Segmentation Results for Common Real-World Infrared Images
In this subsection, we will use the infrared images to evaluate the comprehensive performance of the proposed model. We captured the infrared images of the airplane, ship and tank using the infrared camera at a resolution of 640 × 480. After being pre-processed, including grayscale conversion and edge map calculation, the images are segmented using the proposed model and other traditional algorithms. The segmentation results are then compared.

Segmentation Results for Common Real-World Infrared Images
In this subsection, we will use the infrared images to evaluate the comprehensive performance of the proposed model. We captured the infrared images of the airplane, ship and tank using the infrared camera at a resolution of 640 × 480. After being pre-processed, including grayscale conversion and edge map calculation, the images are segmented using the proposed model and other traditional algorithms. The segmentation results are then compared.
In the experiment, the parameter setting of the proposed model is {K, τ} = {0.2, 1}. The major influences that affect segmentation accuracy are the weak target edges and the interference from the edges of other objects near the target. Figure 8 shows the original infrared images used in the experiment and Figure 9 shows the segmentation results of various active contour models. The last column is the ground truth. From Figure 9, it can be seen intuitively that the propose model can segment the infrared images very accurately and is superior to other traditional models in terms of accuracy. As discussed at the beginning of this section, subjective evaluation has some limitations. Hence, we perform quantitative analysis of these results based on Equations (21)-(23).

Segmentation Results for Common Real-World Infrared Images
In this subsection, we will use the infrared images to evaluate the comprehensive performance of the proposed model. We captured the infrared images of the airplane, ship and tank using the infrared camera at a resolution of 640 × 480. After being pre-processed, including grayscale conversion and edge map calculation, the images are segmented using the proposed model and other traditional algorithms. The segmentation results are then compared. In the experiment, the parameter setting of the proposed model is , = 0.2, 1 . The major influences that affect segmentation accuracy are the weak target edges and the interference from the edges of other objects near the target. Figure 8 shows the original infrared images used in the experiment and Figure 9 shows the segmentation results of various active contour models. The last column is the ground truth. From Figure 9, it can be seen intuitively that the propose model can segment the infrared images very accurately and is superior to other traditional models in terms of  The data in Table 1 intuitively reveal the advantages and disadvantages of the proposed method over other algorithms in terms of infrared image segmentation. Consider the metric of F1 measure, which can reflect the segmentation performance overall. The value of this metric of the proposed method is higher than the other algorithms across the three images. This demonstrates the undisputed superiority of the proposed method.

Segmentation Results of Noise-Corrupted Infrared Images
The images used in the experiment of the previous subsection were captured in the experimental environment and were processed specifically. However, in the real-world applications, the images we obtain are mostly corrupted with the noise, thereby expecting the proposed algorithm to smooth these noises. In order to verify the proposed algorithm's insensitivity to the noise, we corrupt the original infrared images with the noise using the "imnoise" function in MATLAB. In this subsection, we will use these noise-corrupted images for experiments. We add the salt-pepper and The data in Table 1 intuitively reveal the advantages and disadvantages of the proposed method over other algorithms in terms of infrared image segmentation. Consider the metric of F1 measure, which can reflect the segmentation performance overall. The value of this metric of the proposed method is higher than the other algorithms across the three images. This demonstrates the undisputed superiority of the proposed method.

Segmentation Results of Noise-Corrupted Infrared Images
The images used in the experiment of the previous subsection were captured in the experimental environment and were processed specifically. However, in the real-world applications, the images we obtain are mostly corrupted with the noise, thereby expecting the proposed algorithm to smooth these noises. In order to verify the proposed algorithm's insensitivity to the noise, we corrupt the original infrared images with the noise using the "imnoise" function in MATLAB. In this subsection, we will use these noise-corrupted images for experiments. We add the salt-pepper and multiplicative noises to the original images. The images corrupted with the salt-pepper noise are named planeN, shipN and tankN (Figure 10), and the parameter setting of "D" (noise density) is 0.001. The images corrupted with the multiplicative noise are named planeN2, shipN2 and tankN2 (Figure 11), and the parameter setting of "V" (variance) is 0.01. The parameter setting of the proposed model is , = 0.1, 0.2 . In addition to the two major influences discussed in previous subsection, the noise is another factor that affects segmentation accuracy. For some models such as GVF and NBGVF, there is a trade-off between noise smoothing and weak edge protection. This is particularly true for images ship and tank. Some models tradeoff weak edge protection for noise smoothing for the purpose of obtaining better results. However, the proposed model is capable of protecting weak edges without sacrifice of noise smoothing, thereby The parameter setting of the proposed model is , = 0.1, 0.2 . In addition to the two major influences discussed in previous subsection, the noise is another factor that affects segmentation accuracy. For some models such as GVF and NBGVF, there is a trade-off between noise smoothing and weak edge protection. This is particularly true for images ship and tank. Some models tradeoff weak edge protection for noise smoothing for the purpose of obtaining better results. However, the proposed model is capable of protecting weak edges without sacrifice of noise smoothing, thereby The parameter setting of the proposed model is {K, τ} = {0.1, 0.2}. In addition to the two major influences discussed in previous subsection, the noise is another factor that affects segmentation accuracy. For some models such as GVF and NBGVF, there is a trade-off between noise smoothing and weak edge protection. This is particularly true for images ship and tank. Some models tradeoff weak edge protection for noise smoothing for the purpose of obtaining better results. However, the proposed model is capable of protecting weak edges without sacrifice of noise smoothing, thereby resulting in greater accuracy. Figures 12 and 13 shows the segmentation results of various active contour models. As in the previous subsection, we perform quantitative analysis of experimental results using the same three metrics (Precision, Recall, and F1 measure). The analysis results are given in the Tables 2 and 3.         The data show that after the noise is added to the image, the proposed algorithm can still segment images very accurately, as the value of each metric is above 0.9. The accuracy of the proposed method is very close or even higher than the best algorithm, and the proposed method leads the worse algorithms by a larger margin. These results demonstrate the ability of the proposed method to smooth noises in the infrared images more effectively.
To prove our proposed model can be applied to infrared images with different noise intensity, we add different multiplicative noise to the infrared image "plane" by changing the parameter "V" in the "imnoise" function from 0.01 to 0.03 (with the gap of 0.005), and then apply our proposed model to these noise-polluted images and obtain the test results.
According to Figure 14, we can clearly see that the segmentation results are almost the same in the different images. Thus, the results prove that our proposed model can process different noise intensity adaptively. Our proposed model is very robust to the noise. The data show that after the noise is added to the image, the proposed algorithm can still segment images very accurately, as the value of each metric is above 0.9. The accuracy of the proposed method is very close or even higher than the best algorithm, and the proposed method leads the worse algorithms by a larger margin. These results demonstrate the ability of the proposed method to smooth noises in the infrared images more effectively.
To prove our proposed model can be applied to infrared images with different noise intensity, we add different multiplicative noise to the infrared image "plane" by changing the parameter "V" in the "imnoise" function from 0.01 to 0.03 (with the gap of 0.005), and then apply our proposed model to these noise-polluted images and obtain the test results.
According to Figure 14, we can clearly see that the segmentation results are almost the same in the different images. Thus, the results prove that our proposed model can process different noise intensity adaptively. Our proposed model is very robust to the noise.

Discussion
The block diagrams in Figure 15 clearly show the segmentation accuracy of the proposed method and other algorithms. It can be seen that the proposed method is superior to other algorithms for most of the infrared images and is almost insensitive to the noise. The average CPU time and number of iterations in these experiments are shown in Table 4.

Discussion
The block diagrams in Figure 15 clearly show the segmentation accuracy of the proposed method and other algorithms. It can be seen that the proposed method is superior to other algorithms for most of the infrared images and is almost insensitive to the noise. The average CPU time and number of iterations in these experiments are shown in Table 4.

Discussion
The block diagrams in Figure 15 clearly show the segmentation accuracy of the proposed method and other algorithms. It can be seen that the proposed method is superior to other algorithms for most of the infrared images and is almost insensitive to the noise. The average CPU time and number of iterations in these experiments are shown in Table 4.  In order to prove that the proposed method is suitable for more applications, we apply the proposed method to natural images in several sets of experiments. The results are shown in Figure  16.  In order to prove that the proposed method is suitable for more applications, we apply the proposed method to natural images in several sets of experiments. The results are shown in Figure 16. In order to prove that the proposed method is suitable for more applications, we apply the proposed method to natural images in several sets of experiments. The results are shown in Figure  16.  We can see that the proposed method is still capable of segmenting the natural images satisfactorily. Hence, the high segmentation accuracy qualifies the proposed method for more applications.
To sum up, we first perform experiments to verify some basic properties of the proposed method and prove that the proposed method is vastly superior to other classic algorithms in terms of LTI convergence. Afterwards, we apply the proposed method to the infrared images and compare with classic algorithms. Results show that the proposed method has great advantages over other algorithms. Finally, we evaluate the segmentation performance of the proposed method for natural images. Results imply that the proposed method is suitable for more applications.

Conclusions
The infrared image segmentation technology is of great significance to real-world life and manufacturing. However, many issues have yet to be addressed. The research on the use of active contour model for infrared image segmentation is in the infancy, but it has attracted a lot of attention. In this paper, we adapt the active contour model to the infrared images by improving NBGVF. A series of experiments have been performed to prove segmentation accuracy superiority of the proposed method over other algorithms (GVF, GGVF, NGVF, NBGVF, CN-GGVF, LIF, and SOAC). Meanwhile, it is proven that the proposed method can smooth noises while protecting weak edges in the infrared images. Hence, the proposed method is vastly superior to other algorithms.