[go: up one dir, main page]

CN119169034A - Global image segmentation method based on global Gaussian distribution and statistical norm - Google Patents

Global image segmentation method based on global Gaussian distribution and statistical norm Download PDF

Info

Publication number
CN119169034A
CN119169034A CN202411179674.2A CN202411179674A CN119169034A CN 119169034 A CN119169034 A CN 119169034A CN 202411179674 A CN202411179674 A CN 202411179674A CN 119169034 A CN119169034 A CN 119169034A
Authority
CN
China
Prior art keywords
energy functional
global
level set
set function
image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202411179674.2A
Other languages
Chinese (zh)
Other versions
CN119169034B (en
Inventor
刘国军
刘相国
马月梅
岳靖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ningxia University
Original Assignee
Ningxia University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ningxia University filed Critical Ningxia University
Priority to CN202411179674.2A priority Critical patent/CN119169034B/en
Publication of CN119169034A publication Critical patent/CN119169034A/en
Application granted granted Critical
Publication of CN119169034B publication Critical patent/CN119169034B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

The invention provides a global image segmentation method based on global Gaussian distribution and statistical norms, which relates to the technical field of image processing and application and comprises the steps of constructing an energy functional based on the global Gaussian distribution and the statistical norms; the method comprises the steps of obtaining pixel coordinates of an image to be segmented, inputting the pixel coordinates of the image to be segmented into the energy functional, carrying out iterative operation on the level set function in the energy functional until the difference between the energy functional value obtained by iterative calculation of the current round and the energy functional value obtained by iterative calculation of the previous round is smaller than a preset precision threshold value, ending iterative calculation, outputting a final level set function of the image to be segmented, updating the current level set function after each iterative calculation, and carrying out iterative update of the next round by utilizing the updated level set function if the condition of ending the iterative calculation is not reached. The scheme can improve the precision of image segmentation.

Description

Global image segmentation method based on global Gaussian distribution and statistical norms
Technical Field
The invention relates to the technical field of image processing and application, in particular to a global image segmentation method based on global Gaussian distribution and statistical norms.
Background
The image is one of carriers for transmitting information in daily life of people, the research of image technology is extended to various aspects, and the image is widely applied to the fields of target detection, medical image processing, industry, face recognition, intelligent transportation and the like. Image processing has become very important with the continuous progress of technology, and image segmentation is an important field in image processing. The image segmentation mainly divides an image into a plurality of non-overlapping areas, and segments a target area in the image, thereby providing basis for subsequent image analysis and image understanding. Therefore, how to obtain better image segmentation results and improve segmentation accuracy and efficiency have become one of the key problems of image segmentation.
The current common image segmentation mode is a segmentation mode based on a local level set, has a good effect on segmenting an image with uneven intensity, and can improve segmentation accuracy under the condition of uneven intensity by using local statistical characteristics. However, since this approach loses global features, it is highly sensitive to the initial position of the contour. Meanwhile, when the algorithm is processing high noise images, the existing scheme may also perform poorly because noise may cause local intensity variations, which may be erroneously captured by contours, resulting in erroneous segmentation. Therefore, the existing image segmentation scheme has poor accuracy in segmenting the image.
Disclosure of Invention
In view of the foregoing, it is necessary to propose a global image segmentation method based on global gaussian distribution and statistical norms to improve the accuracy of image segmentation.
In a first aspect, the present invention provides a global image segmentation method based on global gaussian distribution and statistical norms, including:
Constructing an energy functional based on global Gaussian distribution and a statistical norm, wherein the energy functional comprises a level set function for representing an image segmentation result;
acquiring pixel coordinates of an image to be segmented;
Inputting pixel coordinates of the image to be segmented into the energy functional, carrying out iterative operation on a level set function in the energy functional based on a gradient descent method until the difference between the energy functional value obtained by the iterative calculation of the current round and the energy functional value obtained by the iterative calculation of the previous round is smaller than a preset precision threshold, and ending the cyclic iterative calculation;
Outputting a final level set function of the image to be segmented;
And after each iterative calculation, if the condition of ending the cyclic iterative calculation is not met, updating the current level set function, and carrying out the iterative updating of the next round by utilizing the updated level set function.
Preferably, the construction of the energy functional based on the global Gaussian distribution and the statistical norm comprises the steps of starting from global information of an image, defining a global Gaussian distribution fitting energy for each pixel in the image to be segmented by taking a global mean value and a global variance as variables, and integrating the fitting energy of all pixel points to construct the energy functional.
Preferably, constructing the energy functional specifically includes:
the fitting energy is described statistically by global gaussian distributions of different mean and variance, respectively, and the following primary energy functional is constructed:
Wherein, Is a probability density function defined as the global region intensity of a gaussian distribution, i=1, 2, Ω 1 and Ω 2 are the inner and outer regions, respectively, of the image contour C to be segmented, the following energy functions can be obtained:
Wherein, The energy functional is characterized, phi is a level set function, x is the pixel coordinates in the image I (x) to be segmented, mu i (x) is the global mean,As a global variance of the values of the variables,For the degree of deviation of the gray value at x from the average intensity value μ i (x), λ (x) is a multiplier that adjusts the deviation from μ i (x), α 1 and α 2 are both fixed balance parameters, and H (Φ (x)) is a Heaviside function.
Preferably, the iterative operation is performed on the level set function in the energy functional based on the gradient descent method, including:
Step S31, configuring initial iteration parameters, wherein the iteration parameters comprise initial energy functional values;
Step S32, respectively calculating the current energy functional related mu i (x), Partial derivative of λ (x);
Step S33, making each calculated partial derivative be 0, and calculating mu i (x) by using pixel coordinates of the image to be segmented, A calculated value of λ (x);
step S34, calculating mu i (x), The calculated value of lambda (x) is brought into the current energy functional to obtain an energy functional value;
And step 35, judging whether the difference between the energy functional value obtained by current calculation and the energy functional value obtained by previous iterative calculation is smaller than a preset precision threshold, if so, ending iterative operation, otherwise, returning to step 32 for execution after updating the level set function in the current energy functional.
Preferably, in step S35, the updating the level set function in the current energy functional includes:
Calculating the partial derivative of the current energy functional with respect to a level set function;
And updating the level set function in the current energy functional by utilizing the partial derivative of the current energy functional with respect to the level set function, and taking a new energy functional formed after the level set function is updated as the energy functional of the next iterative operation.
Preferably, the updating the level set function in the current energy functional with respect to the partial derivative of the level set function by using the current energy functional includes:
updating the level set function using the following calculation formula:
Where Δt is the time step.
Preferably, after constructing the energy functional, further comprising:
correcting the constructed energy functional by using the following regularization term to ensure the smoothness of the level set function in the curve evolution process;
wherein p (phi (x)) is the regularization term, Gradient as a function of level set.
Preferably, the preset precision threshold is 10 -3.
In a second aspect, the present invention provides a computer readable storage medium having stored thereon a computer program which, when executed in a computer, causes the computer to perform the method according to any of the first aspects.
In a third aspect, the present invention provides a computing device comprising a memory having executable code stored therein and a processor, the processor implementing a method as in any of the first aspects when executing the executable code.
According to the technical scheme, when the global image is segmented based on the scheme, firstly, an energy functional is built based on global Gaussian distribution and statistical norms, then, pixel coordinates of the image to be segmented are acquired and then input into the energy functional, further, iterative operation is carried out on a level set function in the energy functional based on a gradient descent method until the difference value between the energy functional value obtained by current iterative calculation and the energy functional value obtained by the previous round is smaller than a preset precision threshold, and then, loop iterative calculation is ended, and a final level set function of the image to be segmented is output. Therefore, in order to solve the problem of high sensitivity to the initial position of the outline and noise caused by the lack of global features, the method constructs the energy functional based on the probability density function of the global area intensity of the Gaussian distribution with statistical rule characteristics from the global information of the image, not only can consider the features of each pixel point of the image to be segmented by taking the global mean value and the variance as variables, but also can measure the similarity among the pixels through the Gaussian distribution statistical characteristics, so that the boundary of the target area can be identified more accurately under the complex background. Therefore, the scheme can improve the precision of image segmentation.
Drawings
Fig. 1 is a flowchart of a global image segmentation method based on global gaussian distribution and statistical norms according to an embodiment of the present invention.
Fig. 2 shows image segmentation contrast results of different initial contours.
Fig. 3 shows image segmentation comparison results of different noise and different levels.
Fig. 4 shows the segmentation comparison result of the natural image.
Detailed Description
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
As shown in fig. 1, the present invention provides a global image segmentation method based on global gaussian distribution and statistical norms, which may include the following steps:
Step 101, constructing an energy functional based on global Gaussian distribution and statistical norms, wherein the energy functional comprises a level set function for representing an image segmentation result;
102, acquiring pixel coordinates of an image to be segmented;
Step 103, inputting pixel coordinates of an image to be segmented into an energy functional, carrying out iterative operation on a level set function in the energy functional based on a gradient descent method until the difference between the energy functional value obtained by iterative computation of the current round and the energy functional value obtained by iterative computation of the previous round is smaller than a preset precision threshold value, and ending the cyclic iterative computation, wherein after each iterative computation, if the condition of ending the cyclic iterative computation is not reached, updating the current level set function, and carrying out iterative update of the next round by utilizing the updated level set function;
and 104, outputting a final level set function of the image to be segmented.
In order to solve the problem that the initial position of the contour and noise are highly sensitive due to the lack of global features, the embodiment of the scheme starts from the global information of the image, constructs the energy functional based on the probability density function of the global area intensity of the Gaussian distribution with statistical rule characteristics, not only considers the feature of each pixel point of the image to be segmented by taking the global mean value and the variance as variables, but also can measure the similarity among pixels through the Gaussian distribution statistical characteristics, so that the boundary of the target area can be identified more accurately under the complex background. Therefore, the scheme can improve the precision of image segmentation.
For step 101, constructing an energy functional based on the global Gaussian distribution and the statistical norm, wherein the energy functional comprises a level set function for representing an image segmentation result;
In this embodiment, in order to solve the problem of high sensitivity to the initial position and noise of the contour caused by the lack of global features, consider that, starting from the global information of the image, for each pixel in the image to be segmented, a new global gaussian distribution fitting energy is defined by taking the global mean and variance as variables, and the fitting energy of all pixel points is integrated to form a data item of an energy function, namely an energy functional. Specifically, when constructing the energy functional, the fitting energy can be described by using global gaussian distribution of different mean and variance from the statistical perspective, and the following primary energy functional is constructed:
Wherein, p i,x(I(x),μi (x), Λ (x)) is a probability density function defined as the intensity of the global region of the gaussian distribution, where,Omega 1 and omega 2 are the inner and outer regions of the image contour C to be segmented, respectively, so that the following energy functions can be obtained:
Wherein, The energy functional is characterized, phi is a level set function, x is the pixel coordinates in the image I (x) to be segmented, mu i (x) is the global mean,As a global variance of the values of the variables,For the degree of deviation of the gray value at x from the average intensity value μ i (x), λ (x) is a multiplier that adjusts the deviation from μ i (x), α 1 and α 2 are both fixed balance parameters, and H (Φ (x)) is a Heaviside function.
For step 102, obtaining pixel coordinates of an image to be segmented;
In the step, firstly, an image to be segmented is obtained, then a pixel coordinate system of the image to be segmented is constructed, and further the pixel coordinate of the image to be segmented is obtained. For example, a two-dimensional coordinate system u-v in units of pixels is established with the upper left corner of the image to be segmented as the origin. The abscissa u and ordinate v of a pixel are the number of columns and rows, respectively, in the image array. The unit of the pixel coordinate system u-v is a pixel, which is a discrete image coordinate or pixel coordinate, with the origin in the upper left corner of the picture.
For step 103, inputting pixel coordinates of the image to be segmented into an energy functional, and carrying out iterative operation on a level set function in the energy functional based on a gradient descent method until the difference between the energy functional value obtained by the iterative computation of the current round and the energy functional value obtained by the iterative computation of the previous round is smaller than a preset precision threshold value, ending the cyclic iterative computation;
In the step, the energy functional value obtained after each iterative operation is compared with the energy functional value obtained by the previous iterative operation, and whether the difference value between the energy functional value and the energy functional value is smaller than a preset precision threshold value is judged. If so, the energy functional is not reduced continuously, namely the iterative operation reaches the precision requirement, the output level set function can meet the segmentation precision requirement, and if not, the iterative operation does not reach the precision requirement, the iterative operation is further continued. Specifically, step 103 may be implemented by performing an iterative operation on the level set function in the energy functional based on the gradient descent method as follows:
Step S31, configuring initial iteration parameters, wherein the iteration parameters comprise initial energy functional values;
Step S32, respectively calculating the current energy functional related mu i (x), Partial derivative of λ (x);
Step S33, making each calculated partial derivative be 0, and calculating mu i (x) by using pixel coordinates of the image to be segmented, A calculated value of λ (x);
step S34, calculating mu i (x), The calculated value of lambda (x) is brought into the current energy functional to obtain an energy functional value;
and step S35, judging whether the difference between the energy functional value obtained by current calculation and the energy functional value obtained by previous iterative calculation is smaller than a preset precision value, if so, ending iterative operation, otherwise, returning to step S32 for execution after updating the level set function in the current energy functional.
In this embodiment, during the iterative operation, each iteration parameter is initialized, for example, an initial energy functional value is given, and then the current energy functional is calculated with respect to μ i (x),Lambda (x) partial derivative. Specifically, the gradient descent flow is: where δ (φ (x)) is the derivative of H (φ (x)), i.e In the formula mu i (x),Λ (x) is unknown, so considering the partial derivative of the energy function and let the derivative equal to zero, it is possible to obtain:
Wherein:
Thus, the calculated mu i (x), And further judging the magnitude relation between the difference value between the current energy functional value and the energy functional value obtained in the previous round and the preset precision value. Of course, if the first iterative operation is performed, the energy functional value of the previous round is the set initial energy functional value. If the difference value of the energy functional values is smaller than a preset precision threshold value, for example, the difference value is smaller than 10 -3, the preset precision requirement is met, along with the progress of iterative operation, the energy functional does not evolve to the minimum, and at the moment, the level set function in the output energy functional is the final level set function of the image to be segmented, and the optimal segmentation curve of the image to be segmented is obtained. If the difference value of the energy functional values is larger than a preset precision threshold, namely the difference value is larger than 10 -3, the precision requirement is not met, and iterative operation is continued after the level set function is updated.
Further, in step S35, when updating the level set function in the current energy functional, the partial derivative of the current energy functional with respect to the level set is calculated, then the partial derivative of the current energy functional with respect to the level set function is utilized to update the level set function in the current energy functional, and the new energy functional formed after the level set function is updated is used as the energy functional of the next iteration operation. Specifically, when the level set function in the current energy functional is updated with respect to the partial derivative of the level set function, the level set function may be updated using the following calculation formula:
Where Δt is the time step.
In addition, in order to avoid the problem of reinitialization and ensure that the curve keeps the level set function smooth in the evolution process, the constructed energy functional can be modified by introducing regularization terms. The specific regularization term may be obtained as follows:
Where p (phi (x)) is a regularization term, Gradient as a function of level set.
And when the regularization term is used for correcting the energy functional, the following new energy functional can be obtained:
in practice, iterative operation can be performed by using the regularized and corrected energy functional, so that a more accurate image segmentation result is obtained.
The effects of the present solution are further described below with reference to specific application examples.
In order to verify the effectiveness of the scheme on natural images, 6 natural images are selected in BSDS data sets to be tested, and the segmentation performance of the scheme is measured according to the following 4 objective evaluation indexes, namely, a Jaccard index (JSC) rate, a Dice Similarity Coefficient (DSC) rate, a Recall rate (Recall) and a Precision rate (Precision). Specifically, each evaluation index is:
in the above four objective evaluation indexes, G s is a divided image, G GT is a GT image, |·|indicates the number of pixels in the region, TP indicates a true example, FP indicates a false positive example, and FN indicates a false negative example.
(1) Verifying the robustness of an initial contour
Experiments were performed by setting initial contours of different sizes and shapes, and the segmentation results are shown in fig. 2. As can be seen from fig. 2, in (a 0)-(a3), a green rectangle represents an initial contour and red represents a division result, in (b 0)-(b3), a green triangle represents an initial contour and red represents a division result, in (c 0)-(c3), a green circle represents an initial contour and red represents a division result, and in (d 0)-(d3), a green trapezoid represents an initial contour and red represents a division result. The whole segmentation result shows that whether the shape of the initial contour is rectangular, round, triangular or trapezoidal, or the initial contour comprises one target area, two target areas, three target areas and even no target area, the scheme can accurately segment the image.
(2) Robustness to noise
In order to verify the robustness of the scheme to noise, the method is verified from two aspects of subjective evaluation and objective evaluation.
1) Subjective evaluation
Setting the initial contour as a circle, selecting a real image, sequentially adding a mean value of 0, gaussian noise with variances of 0.015,0.02,0.025 and 0.03 and spiced salt noise with variances of 0.01,0.012,0.014 and 0.016 into the image, and dividing the image into a segmentation result shown in figure 3. Wherein, (a 0)-(a3) adds the segmentation result of different levels of gaussian noise to the image respectively, and (b 0)-(b3) adds the segmentation result of different levels of salt-and-pepper noise to the image respectively. The segmentation result of fig. 3 shows that the scheme can realize image segmentation of Gaussian noise with different levels and salt and pepper noise with different levels, and has good robustness to Gaussian noise and salt and pepper noise.
The segmentation results of the natural image as shown in fig. 4, (a 1)-(a6) are six different scene segmentation results without noise, and (b 1)-(b6) are six different scene segmentation results with different noise. As can be seen from fig. 4, the present solution can realize the segmentation of natural images under different scenes, and also can realize the segmentation of noise natural images under different scenes. Thereby verifying the feasibility and effectiveness of the scheme on natural image segmentation.
2) Objective evaluation
For better analysis of the experimental results, the above-described fig. 4 was objectively evaluated by four objective evaluation indexes, and the evaluation results are shown in table 1 below:
TABLE 1
Table 1 shows objective evaluation results of the present embodiment on natural image segmentation. The four objective indexes have the value ranges of [0,1], and the closer the objective evaluation value is to 1, the better the segmentation effect is, the closer to 0 is, and the worse the segmentation effect is. From table 1, it can be seen that, for the DSC evaluation index, the segmentation accuracy is in the range of 97.79% -99.48%, for the Precision index, the segmentation accuracy is in the range of 98.67% -99.99%, for the Recall evaluation index, the segmentation accuracy is in the range of 97.16% -99.36%, and for the JSC index, the segmentation accuracy is in the range of 95.68% -99.14%, so that the effectiveness of the scheme in segmenting images is further verified.
The present specification also provides a computer readable storage medium having stored thereon a computer program which, when executed in a computer, causes the computer to perform the method of any of the embodiments of the specification.
The present specification also provides a computing device comprising a memory having executable code stored therein and a processor which when executing the executable code implements the method of any of the embodiments of the specification.
The modules or units in the device of the embodiment of the invention can be combined, divided and deleted according to actual needs. The foregoing disclosure is illustrative of the preferred embodiments of the present invention, and is not to be construed as limiting the scope of the invention, as it is understood by those skilled in the art that all or part of the above-described embodiments may be practiced with equivalents thereof, which fall within the scope of the invention as defined by the appended claims.

Claims (10)

1.一种基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,包括:1. A global image segmentation method based on global Gaussian distribution and statistical norm, characterized by comprising: 基于全局高斯分布和统计范数构建能量泛函;其中,所述能量泛函中包含用于表征图像分割结果的水平集函数;An energy functional is constructed based on a global Gaussian distribution and a statistical norm; wherein the energy functional includes a level set function for characterizing image segmentation results; 获取待分割图像的像素坐标;Get the pixel coordinates of the image to be segmented; 将所述待分割图像的像素坐标输入所述能量泛函,并基于梯度下降法对所述能量泛函中的水平集函数进行迭代运算,直至当前轮次迭代计算得到的能量泛函值与上一轮次迭代计算得到的能量泛函值之间的差值小于预设精度阈值,结束循环迭代计算;Inputting the pixel coordinates of the image to be segmented into the energy functional, and iteratively calculating the level set function in the energy functional based on the gradient descent method, until the difference between the energy functional value calculated in the current round of iteration and the energy functional value calculated in the previous round of iteration is less than a preset accuracy threshold, and then terminating the loop iterative calculation; 输出所述待分割图像的最终水平集函数;Outputting the final level set function of the image to be segmented; 其中,在每次迭代计算后,若未达到结束循环迭代计算的条件,则对当前的水平集函数进行更新,并利用更新后的水平集函数进行下一轮次的迭代更新。Among them, after each iterative calculation, if the condition for ending the loop iterative calculation is not met, the current level set function is updated, and the updated level set function is used for the next round of iterative update. 2.根据权利要求1所述的基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,所述基于全局高斯分布和统计范数构建能量泛函,包括:从图像的全局信息出发,对于待分割图像中的每个像素,将全局均值和方差作为变量定义一个全局高斯分布拟合能量,并对所有像素点的拟合能量进行积分,构建得到能量泛函。2. According to claim 1, the global image segmentation method based on global Gaussian distribution and statistical norm is characterized in that the energy functional is constructed based on the global Gaussian distribution and statistical norm, including: starting from the global information of the image, for each pixel in the image to be segmented, a global Gaussian distribution fitting energy is defined using the global mean and variance as variables, and the fitting energy of all pixels is integrated to construct an energy functional. 3.根据权利要求2所述的基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,构建所述能量泛函具体包括:3. The global image segmentation method based on global Gaussian distribution and statistical norm according to claim 2, characterized in that constructing the energy functional specifically comprises: 从统计学角度分别利用不同均值和方差的全局高斯分布描述所述拟合能量,并构建如下初级能量泛函:From a statistical point of view, the fitting energy is described using global Gaussian distributions with different means and variances, and the following primary energy functional is constructed: 其中,是定义为高斯分布的全局区域强度的概率密度函数,i=1,2,Ω1和Ω2分别为待分割图像轮廓C的内部区域和外部区域,则可以得到如下能量泛函:in, is the probability density function of the global regional intensity defined as a Gaussian distribution, i = 1, 2, Ω 1 and Ω 2 are the inner and outer regions of the image contour C to be segmented, respectively, then the following energy functional can be obtained: 其中,表征能量泛函,φ为水平集函数,x为待分割图像I(x)中的像素坐标,μi(x)为全局均值,为全局方差,为x处的灰度值与平均强度值μi(x)的偏差程度,λ(x)为调整与μi(x)的偏差的乘数,α1和α2均为固定的平衡参数,H(φ(x))是Heaviside函数。in, represents the energy functional, φ is the level set function, x is the pixel coordinate in the image to be segmented I(x), μ i (x) is the global mean, is the global variance, is the degree of deviation between the grayscale value at x and the average intensity value μ i (x), λ(x) is the multiplier for adjusting the deviation from μ i (x), α 1 and α 2 are fixed balance parameters, and H(φ(x)) is the Heaviside function. 4.根据权利要求3所述的基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,所述基于梯度下降法对所述能量泛函中的水平集函数进行迭代运算,包括:4. The global image segmentation method based on global Gaussian distribution and statistical norm according to claim 3, characterized in that the iterative operation of the level set function in the energy functional based on the gradient descent method comprises: 步骤S31:配置初始迭代参数;其中,所述迭代参数包括初始能量泛函值;Step S31: configuring initial iteration parameters; wherein the iteration parameters include initial energy functional values; 步骤S32:分别计算当前能量泛函关于μi(x)、λ(x)的偏导数;Step S32: Calculate the current energy functional with respect to μ i (x), Partial derivatives of λ(x); 步骤S33:令计算得到的各偏导数为0,并利用待分割图像的像素坐标分别计算得到μi(x)、λ(x)的计算值;Step S33: Set the calculated partial derivatives to 0, and use the pixel coordinates of the image to be segmented to calculate μ i (x), The calculated value of λ(x); 步骤S34:将计算得到的μi(x)、λ(x)的计算值带入当前能量泛函,得到能量泛函值;Step S34: The calculated μ i (x), Substitute the calculated value of λ(x) into the current energy functional to obtain the energy functional value; 步骤S35:判断当前计算得到的能量泛函值与上一次迭代计算得到的能量泛函值之间的差值是否小于预设精度阈值;若是,则结束迭代运算;否则,对当前能量泛函中的水平集函数进行更新后,返回步骤S32执行。Step S35: determine whether the difference between the energy functional value currently calculated and the energy functional value calculated in the previous iteration is less than a preset accuracy threshold; if so, end the iterative operation; otherwise, after updating the level set function in the current energy functional, return to step S32 for execution. 5.根据权利要求4所述的基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,步骤S35中,所述对当前能量泛函中的水平集函数进行更新,包括:5. The global image segmentation method based on global Gaussian distribution and statistical norm according to claim 4, characterized in that in step S35, the updating of the level set function in the current energy functional comprises: 计算所述当前能量泛函关于水平集函数的偏导数;Calculating the partial derivative of the current energy functional with respect to the level set function; 利用当前能量泛函关于水平集函数的偏导数对当前能量泛函中的水平集函数进行更新,并水平集函数更新后所形成的新的能量泛函作为下一次迭代运算的能量泛函。The level set function in the current energy functional is updated using the partial derivative of the current energy functional with respect to the level set function, and the new energy functional formed after the level set function is updated is used as the energy functional for the next iterative operation. 6.根据权利要求5所述的基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,所述利用当前能量泛函关于水平集函数的偏导数对当前能量泛函中的水平集函数进行更新,包括:6. The global image segmentation method based on global Gaussian distribution and statistical norm according to claim 5, characterized in that the updating of the level set function in the current energy functional by using the partial derivative of the current energy functional with respect to the level set function comprises: 利用如下计算式对所述水平集函数进行更新:The level set function is updated using the following calculation formula: 其中,Δt为时间步长。Among them, Δt is the time step. 7.根据权利要求3所述的基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,在构建所述能量泛函之后,进一步包括:7. The global image segmentation method based on global Gaussian distribution and statistical norm according to claim 3, characterized in that after constructing the energy functional, it further comprises: 利用如下正则化项对所构建的能量泛函进行修正,以保证曲线演化过程中水平集函数的光滑;The constructed energy functional is modified using the following regularization term to ensure the smoothness of the level set function during the curve evolution process; 其中,p(φ(x))为所述正则化项,为水平集函数的梯度。Where p(φ(x)) is the regularization term, is the gradient of the level set function. 8.根据权利要求1所述的基于全局高斯分布和统计范数的全局图像分割方法,其特征在于,所述预设精度阈值为10-38 . The global image segmentation method based on global Gaussian distribution and statistical norm according to claim 1 , wherein the preset accuracy threshold is 10 −3 . 9.一种计算机可读存储介质,其上存储有计算机程序,当计算机程序在计算机中执行时,令计算机执行如权利要求1-8中任一所述的方法。9. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to execute the method according to any one of claims 1 to 8. 10.一种计算设备,包括存储器和处理器,存储器中存储有可执行代码,处理器执行可执行代码时,实现权利要求1-8中任一所述的方法。10. A computing device comprising a memory and a processor, wherein the memory stores executable codes, and when the processor executes the executable codes, the method according to any one of claims 1 to 8 is implemented.
CN202411179674.2A 2024-08-27 2024-08-27 Global image segmentation method based on global Gaussian distribution and statistical norm Active CN119169034B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411179674.2A CN119169034B (en) 2024-08-27 2024-08-27 Global image segmentation method based on global Gaussian distribution and statistical norm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411179674.2A CN119169034B (en) 2024-08-27 2024-08-27 Global image segmentation method based on global Gaussian distribution and statistical norm

Publications (2)

Publication Number Publication Date
CN119169034A true CN119169034A (en) 2024-12-20
CN119169034B CN119169034B (en) 2025-06-17

Family

ID=93886314

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411179674.2A Active CN119169034B (en) 2024-08-27 2024-08-27 Global image segmentation method based on global Gaussian distribution and statistical norm

Country Status (1)

Country Link
CN (1) CN119169034B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106570882A (en) * 2016-10-28 2017-04-19 中国科学院苏州生物医学工程技术研究所 Active contour image segmentation method for Gaussian mixture model
EP3188127A1 (en) * 2015-12-29 2017-07-05 Laboratoires Bodycad Inc. Method and system for performing bone multi-segmentation in imaging data
CN107798684A (en) * 2017-11-07 2018-03-13 河南师范大学 A kind of active contour image partition method and device based on symbol pressure function
CN111145179A (en) * 2019-11-20 2020-05-12 昆明理工大学 Gray uneven image segmentation method based on level set
CN116843706A (en) * 2023-05-12 2023-10-03 广西电网有限责任公司贵港供电局 An energy-adaptive fast level set image segmentation method
WO2024055446A1 (en) * 2022-09-13 2024-03-21 深圳先进技术研究院 Image segmentation method and apparatus, device, and readable storage medium

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3188127A1 (en) * 2015-12-29 2017-07-05 Laboratoires Bodycad Inc. Method and system for performing bone multi-segmentation in imaging data
CN106570882A (en) * 2016-10-28 2017-04-19 中国科学院苏州生物医学工程技术研究所 Active contour image segmentation method for Gaussian mixture model
CN107798684A (en) * 2017-11-07 2018-03-13 河南师范大学 A kind of active contour image partition method and device based on symbol pressure function
CN111145179A (en) * 2019-11-20 2020-05-12 昆明理工大学 Gray uneven image segmentation method based on level set
WO2024055446A1 (en) * 2022-09-13 2024-03-21 深圳先进技术研究院 Image segmentation method and apparatus, device, and readable storage medium
CN116843706A (en) * 2023-05-12 2023-10-03 广西电网有限责任公司贵港供电局 An energy-adaptive fast level set image segmentation method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
SONG YANGYANG ET AL: "Active contours driven by Gaussian function and adaptive-scale local correntropy-based K-means clustering for fast image segmentation", SIGNAL PROCESSING, 30 September 2020 (2020-09-30) *
ZHANG DAN ET AL: "A Xanthoceras sorbifolium crack segmentation method based on an improved level set", IET IMAGE PROCESSING, vol. 17, no. 5, 30 April 2023 (2023-04-30) *
孙林;柯正友;冯小博;徐久成;薛占熬;: "基于双符号压力函数的活动轮廓图像分割方法", 计算机工程与应用, no. 20, 30 April 2018 (2018-04-30) *
李更生 等: "基于区域信息耦合的自适应图像分割", 激光与光电子学进展, vol. 59, no. 2, 31 January 2022 (2022-01-31), pages 1 - 3 *
李更生,刘国军: "融合统计范数度量的局部区域图像分割", 光学精密工程, vol. 32, no. 16, 25 August 2024 (2024-08-25), pages 1 - 3 *
王辉;吴永武;杜应琼;: "一种图像分割模型的自适应参数选择方法", 计算机应用与软件, no. 05, 31 May 2019 (2019-05-31) *

Also Published As

Publication number Publication date
CN119169034B (en) 2025-06-17

Similar Documents

Publication Publication Date Title
CN111795704B (en) Method and device for constructing visual point cloud map
JP5506785B2 (en) Fingerprint representation using gradient histogram
CN108388896B (en) License plate identification method based on dynamic time sequence convolution neural network
CN110334762B (en) Feature matching method based on quad tree combined with ORB and SIFT
US11017210B2 (en) Image processing apparatus and method
US12159405B2 (en) Method for detecting medical images, electronic device, and storage medium
CN110807473A (en) Target detection method, device and computer storage medium
CN109658406A (en) Recognition methods, device, computer equipment and the storage medium of blood-vessel image
CN114596320A (en) A kind of image segmentation method and device based on ALSLCV model
CN109101985A (en) It is a kind of based on adaptive neighborhood test image mismatch point to elimination method
CN106485651A (en) The image matching method of fast robust Scale invariant
CN110019895A (en) A kind of image search method, device and electronic equipment
CN114972821B (en) An image matching method, apparatus and electronic device
CN114723791B (en) Feature point matching method, device, electronic device and storage medium
CN117435934A (en) Matching method, device and storage medium for moving target trajectories based on bipartite graph
CN107844803A (en) The method and apparatus that a kind of picture compares
KR101725041B1 (en) Method and apparatus for detecting forged image
CN109002758B (en) Face feature point positioning method, device, equipment and storage medium
CN113379788A (en) Target tracking stability method based on three-element network
CN119169034B (en) Global image segmentation method based on global Gaussian distribution and statistical norm
CN114519729B (en) Image registration quality assessment model training method, device and computer equipment
US11074427B2 (en) Method for reconstructing an imprint image from image portions
CN116645525B (en) Game image recognition method and processing system
CN119579990B (en) A method, apparatus and readable storage medium for small target recognition in infrared images
CN113658235B (en) Accurate registration method of optical remote sensing image based on VGG network and Gaussian difference network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant