CN107895373B - Image segmentation method based on local region consistency manifold constraint MRF model - Google Patents
Image segmentation method based on local region consistency manifold constraint MRF model Download PDFInfo
- Publication number
- CN107895373B CN107895373B CN201711123155.4A CN201711123155A CN107895373B CN 107895373 B CN107895373 B CN 107895373B CN 201711123155 A CN201711123155 A CN 201711123155A CN 107895373 B CN107895373 B CN 107895373B
- Authority
- CN
- China
- Prior art keywords
- mrf
- image
- manifold
- model
- local
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-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 discloses an image segmentation method based on a local region consistency manifold constraint MRF model, which is based on a PairwiseMRF model, an image segmentation model based on the local region MRF model is constructed on an expansion neighborhood of an MRF node, and the prior distribution of a local region effectively avoids the interference of noise or texture mutation; meanwhile, based on manifold learning theory, a manifold regular term under a probability framework is established, the prior of the complex local space geometric structure of the natural image is effectively described by utilizing the probability distribution of the local area, and the local space geometric structure described by the manifold learning is introduced into the local area MRF segmentation model. Experiments prove that compared with the prior method, the method disclosed by the invention not only avoids the over-smoothing punishment of local area prior, but also effectively keeps the local geometric structure information of image segmentation.
Description
Technical Field
The invention belongs to the technical field of image processing, and particularly relates to an image segmentation method.
Background
An image segmentation method based on a Markov Random Field (MRF) model is widely applied, and the segmentation method is based on local spatial correlation of image information, uses a two-dimensional Random Field model to describe feature information of an image, and enhances the correlation of structure information by the MRF-based image segmentation model, so that the speed and the precision of natural image segmentation are greatly improved, but for natural images with abundant statistical characteristics, the currently common point-to-point interaction structure-based MRF (pairwise MRF) model cannot sufficiently describe complex statistical information and prior knowledge of the natural images, and therefore, in the image segmentation problem, the point-to-point MRF-based image processing algorithm often has wrong segmentation in the regions or edges of the images.
The complex natural image usually has a high-dimensional non-gaussian statistical characteristic, and the prior knowledge of the image needs to be accurately constructed to model the dependency relationship of the image structure in an extended neighborhood. A brain MR image segmentation algorithm [ J ] automatic science report based on an image slice Markov random field, 2014(08): 1754-. Zhufeng et al [ Zhufeng, Luo Li Min, Song Yun Qing et al. ] image segmentation [ J ] based on adaptive spatial neighborhood information Gaussian mixture model computer research and development, 2011,48(11): 2000) plus 2007 ] establishes weighted distribution of neighborhood information through the median probability of pixels in a local region, can effectively suppress noise and maintain edges, but the median probability still cannot completely reflect the characteristic spatial relationship of the local region.
Although the conventional local area MRF-based method adopts a distance measurement mode, such as euclidean distance, K-L divergence, etc., to introduce spatial information of local areas of an image and to serve as a constraint on consistency of the local areas, in the euclidean space, a sample point which represents a short distance in the space is not necessarily the closest distance in the space where an object itself is located. Therefore, the local region MRF method using the distance measurement cannot effectively describe the global consistency characteristics of the complex high-dimensional data, resulting in the erroneous segmentation of the complex natural image.
Disclosure of Invention
The invention aims to provide an image segmentation method based on a local region consistency manifold constraint MRF model, which aims to solve the problem that the conventional region MRF image segmentation model based on the Euclidean space distance measurement technology cannot describe the global consistency characteristic of complex high-dimensional data; the invention can effectively improve the image segmentation effect.
In order to achieve the purpose, the invention adopts the following technical scheme:
the image segmentation method based on the local region consistency manifold constraint MRF model comprises the following steps:
step 1: inputting a natural image X to be segmented as X ═ X1,x2,…xN};
Step 2: initializing parameters: segmentation class number K, local region prior Potts model parameter beta, Lagrange multiplier lambda and Gibbs sampling algorithm initial temperature T0;
And step 3: constructing a nearest neighbor graph G which represents the local geometric structure of the image data manifold, wherein V is a vertex set, E is an edge set, and W is a weight matrix of the graph G;
and 4, step 4: searching low-dimensional representation of manifold through graph embedding, projecting sample data in high-dimensional space complex distribution to low-dimensional space, establishing a dimension reduction mapping relation in a local area by utilizing the homomorphism of the manifold in the local area of an image and a Euclidean space, establishing a geometric structure relation representing the manifold in the low-dimensional space, and establishing a manifold regular term under a probability frame;
and 5: based on a pairwise MRF model, modeling the dependency relationship of an image structure on an expansion neighborhood, expressing more complex image region characteristics, and establishing a region-based MRF energy segmentation model; introducing a consistent manifold constraint term of a local area of the image by utilizing a Lagrange multiplier method to obtain a consistent manifold constraint MRF image segmentation model of the local area;
step 6: method for estimating GMM parameter of local region consistency manifold constraint MRF image segmentation model provided by the invention by adopting expectation maximization algorithm
And 7, optimizing the consistent manifold constraint MRF segmentation model of the local area by adopting a Gibbs sampling algorithm, and outputting an image segmentation result.
Further, step 2 specifically includes:
2a) let Ω ═ {1,2, …, K } denote the pixel node label space, and the segmentation class number K is determined manually;
2b) the prior Potts model parameter beta of the local region belongs to [0.1,5], and the Lagrangian multiplier lambda of the manifold consistency region constraint term of the local region belongs to [10,100 ];
2c) initial temperature T of Gibbs sampling algorithm0=4.0。
Further, step 3 specifically includes:
given image observation sample set X ═ X1,x2,…xNConstructing a nearest neighbor graph G (V, E, W) representing a local geometric structure of an image data manifold in a sample space set X based on manifold learning theory, wherein V is a vertex set, E is an edge set, and W is a weight matrix of the graph G;
in the formula, Nk(xs) Is a sample point xsK neighbor set of (1), if xs∈Nk(xr) Or xr∈Nk(xs) Then sample point xsAnd xrIs two adjacent sample points with weight wsrAccording to the thermal kernel function representation method weight wsrThe definition is as follows:
wherein t is 5, which is a thermonuclear parameter.
Further, step 4 specifically includes:
4a) under the probabilistic framework, let the observation sample set X ═ { X1,x2,…xNThe corresponding implicit variable set is Y ═ Y1,y2,…yNIn which y iss∈{1,2,…K, wherein K is the total number of the segmentation categories;
4b) let the smooth insert fk(x) Is a conditional probability distribution p (y ═ k | x), i.e.:
4c) based on a spectrogram theory, utilizing Laplacian feature mapping, establishing a nearest neighbor graph describing the local geometric structure of the image based on the step 3, and then searching for low-dimensional representation of the manifold through graph embedding; local information measure S of the established conditional probability distribution p (y ═ k | x)kThe definition is as follows:
D(p(k|xs),p(k|xr))2representing two adjacent nodes xs,xrConditional probability distribution p (y) ofs|xs) And p (y)r|xr) The distance of (c).
Further, step 5 specifically includes the following steps:
5a) based on a pairwise MRF model, more complex image region characteristics are expressed by utilizing the dependency relationship of an image structure on an expansion neighborhood, and the established region-based MRF energy segmentation model is shown as the following formula:
in the formula: eDataRepresenting the likelihood energy term of the MRF model in the local region, p (x)s|ys,Θ)Indicates a given condition (y)sΘ) time node xsTheta is a parameter of the conditional likelihood distribution, deltasA set of neighborhood nodes representing a node s; eSmoothThe prior smoothing energy term of the MRF model is obtained, and beta is a prior parameter of the MRF; w is ar(ys,yr) As a weight function of the center node s and the neighborhood nodes r, a standard normal distribution function is described as follows:
in the formula:representing a local area deltasA set of labels of neighborhood nodes with equal label values of the central node s is neutralized; i S (yr) I represents the set S (y)r) A potential of (d);
5b) based on the MRF global energy segmentation model of the established local region model, introducing a consistent manifold constraint term of the local region of the image by utilizing a Lagrange multiplier method to obtain the MRF image segmentation model of the consistent manifold constraint of the local region as shown in the following formula:
6. the image segmentation method based on the local region consistency manifold constraint MRF model according to claim 1, wherein the step 6 specifically comprises:
wherein, | δsAnd | is the number of neighborhood nodes of the central node s.
And 7, optimizing the consistent manifold constraint MRF segmentation model of the local region by adopting a Gibbs sampling algorithm, and estimating an image segmentation result based on a Maximum A Posteriori (MAP) criterion. Further, step 7 specifically includes:
7a) setting the maximum iteration number t of the Gibbs sampling algorithmmaxInitial temperature T04.0, annealing speed 0.95;
7b) sampling a consistent manifold constraint MRF segmentation model of a local area by adopting a Gibbs sampling algorithm, wherein a node y corresponding to a central pixelsBased on maximum a posteriori criteria (MAP):accept a new tag k according to the following probability:
in the formula (I), the compound is shown in the specification,for estimated optimal label, ELocal(k) Local region energy given a central pixel label of k;
7c) calculating pointsGlobal energy of cut modelIf it is notOr a maximum number of iterations t is reachedmaxThen, thenOutputting a segmentation result; otherwise, the temperature T is reduced to 0.95T(t)Returning to step 7b) and continuing the iteration.
Compared with the prior art, the invention has the advantages that:
based on a Pairwise MRF model, an image segmentation model based on a local region MRF model is constructed in an extended neighborhood of an MRF node, and the prior distribution of a local region effectively avoids the interference of noise or texture mutation; meanwhile, based on manifold learning theory, manifold regular terms under a probability framework are established, the local area probability distribution is utilized to effectively describe the prior of the complex local space geometric structure of the natural image, and the local space geometric structure described by the manifold learning is introduced into the local area MRF segmentation model.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 shows experimental results of an embodiment of the present invention; wherein FIG. 2(a) is a beach, mountain, and windsurf artwork; FIG. 2(b) shows the results of the beach, mountain, and windsurfing based on Pairwise MRF segmentation model; FIG. 2(c) shows the results of the segmentation of beach, mountain, and windsurfing based on the Local Space Adaptive MRF (LSAMLF) model; FIG. 2(d) shows the segmentation results of beacon, mountain, and windsurfing models proposed by the present invention.
Detailed Description
The present invention is further described in detail below with reference to the attached drawings so that those skilled in the art can implement the invention by referring to the description text.
Referring to fig. 1, an image segmentation method based on a local region consistency manifold constraint MRF model according to the present invention includes the following steps:
step 1: inputting a natural image X to be segmented as X ═ X1,x2,…xs…xN},xsRepresenting a pixel;
step 2: initializing parameters: segmentation class number K, local region prior Potts model parameter beta, Lagrange multiplier lambda and Gibbs sampling algorithm initial temperature T0;
2a) Let Ω ═ {1,2, …, K } denote the pixel node label space, and the segmentation class number K is determined manually.
2b) In the embodiment of the invention, the prior Potts model parameter beta of the local region belongs to [0.1,5], and the Lagrangian multiplier lambda belongs to [10,100] of the constraint term of the manifold consistency region of the local region.
2c) Initial temperature T of Gibbs sampling algorithm in the embodiment of the invention0=4.0。
And step 3: creating nearest neighbor graph describing local geometry of image
Given image observation sample set X ═ X1,x2,…xNBased on manifold learning theory, constructing nearest neighbor graph G (V, E, W) representing local geometric structure of image data manifold in sample space set X, wherein V is vertex set, E is edge set, W is weight matrix of graph G,
in the formula, Nk(xs) Is a sample point xsK neighbor set of (1), if xs∈Nk(xr) Or xr∈Nk(xs) Then sample point xsAnd xrIs two adjacent sample points with weight wsrAccording to the thermal kernel function representation method weight wsrThe definition is as follows:
wherein t is 5-20, which is a thermonuclear parameter.
And 4, step 4: establishing a manifold regular term under a probability framework:
4a) under the probabilistic framework, let the observation sample set X ═ { X1,x2,…xNThe corresponding implicit variable set is Y ═ Y1,y2,…yNIn which y issE {1,2, … K }, where K is the total number of segmentation classes.
4b) Let the smooth insert fk(x) Is a conditional probability distribution p (y ═ k | x), i.e.:
according to the manifold assumption that points close to each other in the manifold are most likely to belong to the same class, the conditional probability distributions p (y) of the two sample pointss|xs) And p (y)r|xr) It is also approximate, that is, in the implicit geometry of the marginal probability distribution p (X) of the observation sample set X, the conditional probability distribution p (y | ·) is smooth along the geodesic line.
4c) And (3) establishing a nearest neighbor graph describing the local geometric structure of the image based on the spectrogram theory by utilizing Laplacian feature mapping and then searching for a low-dimensional representation of the manifold through graph embedding. Thus, the sample data is projected to a low-dimensional space in a complex distribution of a high-dimensional space, a dimensionality reduction mapping relation is established locally by utilizing the homomorphism of the manifold in a local and Euclidean space, and a geometric structure relation expressing the manifold is established in the low-dimensional space. Local information measure S of the established conditional probability distribution p (y ═ k | x)kThe definition is as follows:
D(p(k|xs),p(k|xr))2representing two adjacent nodes xs,xrConditional probability distribution p (y) ofs|xs) And p (y)r|xr) The distance of (c).
Due to local information measure SkIs calculated along the shortest distance on the manifold and is therefore based on the local information measure SkThe measure of (a) more accurately reflects the prior information of the spatial distribution of the data.
And 5: establishing a local region consistency manifold constraint MRF segmentation model
5a) Based on a pairwise MRF model, more complex image region characteristics are expressed by utilizing the dependency relationship of an image structure on an expansion neighborhood, and the established region-based MRF energy segmentation model is shown as the following formula:
in the formula: eDataRepresenting the likelihood energy term of the MRF model in the local region, p (x)s|ysΘ) represents a given condition (y)sΘ) time node xsTheta is a parameter of the conditional likelihood distribution, deltasA set of neighborhood nodes representing a node s; eSmoothIs the prior smoothing energy term of the MRF model, and beta is the prior parameter of the MRF. w is ar(ys,yr) As a weight function of the center node s and the neighborhood nodes r, a standard normal distribution function is described as follows:
in the formula:representing a local area deltasA set of labels of neighborhood nodes equal to the label value of the central node s; i S (yr) I represents the set S (y)r) The potential of (c).
5b) Based on the MRF global energy segmentation model of the established local region model, introducing a consistent manifold constraint term of the local region of the image by utilizing a Lagrange multiplier method to obtain the MRF image segmentation model of the consistent manifold constraint of the local region as shown in the following formula:
as can be seen from the above equation, the local region manifold constraint MRF image segmentation model includes three terms: (1) in the Gaussian distribution likelihood model with local area weight, the weight wr(ys,yr) Takes into account spatial information between similar pixels; (2) the prior distribution of the local area effectively avoids the interference of noise or texture mutation, and keeps some obvious local structures of the natural image; (3) regularization termMore complex high-dimensional manifold structure information is introduced for spatial geometry structure prior of local region probability distribution, and over-smooth punishment of local region prior is avoided. λ is a regularization term parameter that controls the complexity of the implicit geometric function in the probability distribution. Therefore, the model fully utilizes the gray information and the complex spatial structure information of the image, has good robustness on the noise and the texture mutation of the image, and keeps the detail characteristics of the image.
Step 6: parameter estimation for proposed models
wherein, | δsAnd | is the number of neighborhood nodes of the central node s.
Step 7, optimizing the proposed model by adopting a Gibbs sampling algorithm
7a) Setting the maximum iteration number t of the Gibbs sampling algorithmmaxInitial temperature T04.0, annealing speed 0.95.
7b) Sampling a local region consistency manifold constraint MRF segmentation model by adopting a Gibbs sampling algorithm, and performing MAP criterionFor the estimated optimal label, the node y corresponding to the central pixelsAccept a new tag k according to the following probability:
in the formula, ELocal(k) Local region energy given a central pixel label of k;
7c) computing global energy of a segmentation modelIf it is not(ε is a small threshold, typically 10-6) Or a maximum number of iterations t is reachedmaxThen, thenAnd outputting a segmentation result. Otherwise, the temperature T is reduced to 0.95T(t)Returning to step 7b) and continuing the iteration.
The effect of the present invention will be further described with reference to fig. 2.
FIG. 2(a) shows a natural image to be segmented, wherein three images from left to right are respectively segmented into 4 classes, 4 classes and 6 classes; FIG. 2(b) is the result of the Pairwise MRF based segmentation model; FIG. 2(c) is a segmentation result based on a Local Spatial Adaptive MRF (LSAMLF) model; FIG. 2(d) shows the segmentation result of the present invention.
From the comparison result of the segmentation experiment, the segmentation result based on the Pairwise MRF segmentation model is the worst, because the point pair interaction structure of the Pairwise MRF model is too simple, the region and the global prior of a complex image are difficult to effectively express, and the image segmentation result is greatly interfered by factors such as noise or texture mutation of the image. As shown in the segmentation result of the "beach" image in fig. 2(b), the segmentation regions such as "coconut leaf" and "sea water" have strong abrupt texture changes, so that the segmentation effect of the "coconut leaf" and "sea water" regions is poor, and partial erroneous segmentation occurs. In the segmentation result of the "mountain" image in fig. 2(b), the "grassland" and "rock" regions have significant segmentation noise. In the segmentation result of the "windsurfing" image in fig. 2(b), the segmentation of the "sea water" region has more segmentation speckle noise, and particularly, the upper half part of the "sail" has obvious segmentation errors. The LSAMRF segmentation model has the advantages that the region characteristics of the image are introduced and more image region information is contained, so that the model has better robustness to sudden changes of noise or texture signals, and the better image segmentation result reflects the superiority of the model. The segmentation results are respectively shown in fig. 2(c), and it can be seen from the segmentation results that the "coconut leaf" and "sea water" of the "beach" image, "grassland" and "rock" region of the "mountain" image, "sea water" region of the "windsurfing" image, etc. have a smoother segmentation effect, and the upper half of the "sail" in the "windsurfing" image can be segmented more correctly, but because of the solution process of the region MRF model, the energy minimization of the local region easily causes the "sitting competition" phenomenon, which causes the phenomenon of "edge band" at the segmentation edge of some images, as shown in the "mountain" image segmentation result of fig. 2(c), an edge band appears between the "sky" and "mountain". The model provided by the invention firstly considers the spatial information among similar pixels by using the local spatial adaptive GMM, effectively avoids the interference of noise or texture mutation, and avoids the prior over-smooth punishment of a conventional local area; then, a local manifold constraint term of the image is established, complex space geometric structure prior of a local area is effectively described, and some obvious local structures of a natural image are kept. From the image segmentation result, the model provided by the invention effectively solves the edge banding phenomenon of the image segmentation result of "mountain" in fig. 2(c), and in the image segmentation result of "windsurfing" in fig. 2(a), the model not only has better robustness in the segmentation of the "ocean" region, but also keeps some detail features of the image, such as the "horizontal handle" in the "windsurfing" image, and is effectively extracted. Therefore, the model provided by the invention fully utilizes the gray information and the complex spatial structure information of the image, has good robustness on the noise and the texture mutation of the image, and keeps the detail characteristics of the image.
Claims (7)
1. The image segmentation method based on the local region consistency manifold constraint MRF model is characterized by comprising the following steps of:
step 1: inputting a natural image X to be segmented as X ═ X1,x2,…xN};
Step 2: initializing parameters: number of segmentation class K, local area priorInitial temperature T of Potts model parameter beta, Lagrange multiplier lambda and Gibbs sampling algorithm0;
And step 3: constructing a nearest neighbor graph G which represents the local geometric structure of the image data manifold, wherein V is a vertex set, E is an edge set, and W is a weight matrix of the graph G;
and 4, step 4: searching low-dimensional representation of manifold through graph embedding, projecting sample data in high-dimensional space complex distribution to low-dimensional space, establishing a dimension reduction mapping relation in a local area by utilizing the homomorphism of the manifold in the local area of an image and a Euclidean space, establishing a geometric structure relation representing the manifold in the low-dimensional space, and establishing a manifold regular term under a probability frame;
and 5: based on a pairwise MRF model, modeling the dependency relationship of an image structure on an expansion neighborhood, expressing more complex image region characteristics, and establishing a region-based MRF energy segmentation model; introducing a consistent manifold constraint term of a local area of the image by utilizing a Lagrange multiplier method to obtain a consistent manifold constraint MRF image segmentation model of the local area;
step 6: method for estimating GMM parameter of local region consistency manifold constraint MRF image segmentation model provided by the invention by adopting expectation maximization algorithm
And 7, optimizing the consistent manifold constraint MRF segmentation model of the local area by adopting a Gibbs sampling algorithm, and outputting an image segmentation result.
2. The image segmentation method based on the local region consistency manifold constraint MRF model according to claim 1, wherein the step 2 specifically comprises:
2a) let Ω ═ {1,2, …, K } denote the pixel node label space, and determine the segmentation class number K;
2b) the prior Potts model parameter beta of the local region belongs to [0.1,5], and the Lagrangian multiplier lambda of the manifold consistency region constraint term of the local region belongs to [10,100 ];
2c) initial temperature of Gibbs sampling algorithmT0=4.0。
3. The image segmentation method based on the local region consistency manifold constraint MRF model according to claim 2, wherein the step 3 specifically comprises:
giving a natural image X to be segmented as X1,x2,…xNConstructing a nearest neighbor graph G (V, E, W) representing a local geometric structure of an image data manifold in a natural image X to be segmented based on manifold learning theory, wherein V is a vertex set, E is a side set, and W is a weight matrix of the graph G;
in the formula, Nk(xs) Is a sample point xsK neighbor set of (1), if xs∈Nk(xr) Or xr∈Nk(xs) Then sample point xsAnd xrIs two adjacent sample points with weight wsrAccording to the thermal kernel function representation method weight wsrThe definition is as follows:
wherein t is 5, which is a thermonuclear parameter.
4. The image segmentation method based on the local region consistency manifold constraint MRF model according to claim 3, wherein the step 4 specifically comprises:
4a) under a probability framework, a natural image X to be segmented is made to be { X ═ X1,x2,…xNThe corresponding implicit variable set is Y ═ Y1,y2,…yNIn which y iss∈{1,2,…K};
4b) Let the smooth insert fk(x) Is a conditional probability distribution p (y ═ k | x), i.e.:
4c) Based on a spectrogram theory, utilizing Laplacian feature mapping, establishing a nearest neighbor graph describing the local geometric structure of the image based on the step 3, and then searching for low-dimensional representation of the manifold through graph embedding; local information measure S of the established conditional probability distribution p (y ═ k | x)kThe definition is as follows:
D(p(k|xs),p(k|xr))2representing two adjacent nodes xs,xrConditional probability distribution p (y) ofs|xs) And p (y)r|xr) The distance of (c).
5. The image segmentation method based on the local region consistency manifold constraint MRF model according to claim 4, wherein the step 5 specifically comprises the following steps:
5a) based on a pairwise MRF model, more complex image region characteristics are expressed by utilizing the dependency relationship of an image structure on an expansion neighborhood, and the established region-based MRF energy segmentation model is shown as the following formula:
in the formula: eDataRepresenting the likelihood energy term of the MRF model in the local region, theta being a parameter of the conditional likelihood distribution, deltasA set of neighborhood nodes representing a node s; eSmoothIs the a priori smoothed energy term of the MRF model,beta is a prior parameter of MRF; w is ar(ys,yr) As a weight function of the center node s and the neighborhood nodes r, a standard normal distribution function is described as follows:
in the formula:representing a local area deltasA set of labels of neighborhood nodes with equal label values of the central node s is neutralized; i S (y)r) I denotes the set S (y)r) A potential of (d);
5b) based on the MRF global energy segmentation model of the established local region model, introducing a consistent manifold constraint term of the local region of the image by utilizing a Lagrange multiplier method to obtain the MRF image segmentation model of the consistent manifold constraint of the local region as shown in the following formula:
6. the image segmentation method based on the local region consistency manifold constraint MRF model according to claim 5, wherein the step 6 specifically comprises:
wherein, | δsAnd | is the number of neighborhood nodes of the central node s.
7. The image segmentation method based on the local region consistency manifold constraint MRF model according to claim 6, wherein the step 7 specifically comprises:
7a) setting the maximum iteration number t of the Gibbs sampling algorithmmaxInitial temperature T04.0, annealing speed 0.95;
7b) sampling a consistent manifold constraint MRF segmentation model of a local area by adopting a Gibbs sampling algorithm, wherein a node y corresponding to a central pixelsAccept a new tag k according to the following probability:
in the formula, ELocal(k) Local region energy given a central pixel label of k;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711123155.4A CN107895373B (en) | 2017-11-14 | 2017-11-14 | Image segmentation method based on local region consistency manifold constraint MRF model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711123155.4A CN107895373B (en) | 2017-11-14 | 2017-11-14 | Image segmentation method based on local region consistency manifold constraint MRF model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107895373A CN107895373A (en) | 2018-04-10 |
CN107895373B true CN107895373B (en) | 2021-04-27 |
Family
ID=61805199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711123155.4A Active CN107895373B (en) | 2017-11-14 | 2017-11-14 | Image segmentation method based on local region consistency manifold constraint MRF model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107895373B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108629134B (en) * | 2018-05-11 | 2022-05-03 | 浙江工业大学 | Similarity strengthening method for small fields in manifold |
CN109785331B (en) * | 2018-12-04 | 2023-01-06 | 山东大学 | Sonar image segmentation method based on self-adaptive pixel value constraint and MRF |
CN111008977B (en) * | 2019-12-02 | 2023-04-25 | 西安建筑科技大学 | Image segmentation method of high-order MRF model based on multi-node topological overlap measure |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657995A (en) * | 2015-02-12 | 2015-05-27 | 合肥工业大学 | Remote sensing image segmentation method adopting region splitting technology |
CN107292346A (en) * | 2017-07-05 | 2017-10-24 | 四川大学 | A kind of MR image hippocampus partitioning algorithms learnt based on Local Subspace |
CN107316306A (en) * | 2017-06-23 | 2017-11-03 | 电子科技大学 | A kind of diameter radar image fast partition method based on Markov model |
CN107330907A (en) * | 2017-06-28 | 2017-11-07 | 陕西师范大学 | A kind of MRF image partition methods of combination deep learning shape prior |
-
2017
- 2017-11-14 CN CN201711123155.4A patent/CN107895373B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657995A (en) * | 2015-02-12 | 2015-05-27 | 合肥工业大学 | Remote sensing image segmentation method adopting region splitting technology |
CN107316306A (en) * | 2017-06-23 | 2017-11-03 | 电子科技大学 | A kind of diameter radar image fast partition method based on Markov model |
CN107330907A (en) * | 2017-06-28 | 2017-11-07 | 陕西师范大学 | A kind of MRF image partition methods of combination deep learning shape prior |
CN107292346A (en) * | 2017-07-05 | 2017-10-24 | 四川大学 | A kind of MR image hippocampus partitioning algorithms learnt based on Local Subspace |
Non-Patent Citations (2)
Title |
---|
A MRF model-based segmentation approach to classification for multispectral imagery;Anjan Sarkar 等;《IEEE》;20020531;第40卷(第5期);全文 * |
视觉感受与Markov随机场相结合的高分辨率遥感影像分割法;许妙忠 等;《测绘学报》;20150228;第44卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107895373A (en) | 2018-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107895373B (en) | Image segmentation method based on local region consistency manifold constraint MRF model | |
Kohli et al. | Robust higher order potentials for enforcing label consistency | |
Lucchi et al. | Learning for structured prediction using approximate subgradient descent with working sets | |
US20100134487A1 (en) | 3d face model construction method | |
Ayed et al. | Graph cut segmentation with a global constraint: Recovering region distribution via a bound of the bhattacharyya measure | |
Li et al. | SOLD: Sub-optimal low-rank decomposition for efficient video segmentation | |
CN109934843B (en) | Real-time contour refinement matting method and storage medium | |
US9569699B2 (en) | System and method for synthesizing portrait sketch from a photo | |
CN109242019B (en) | Rapid detection and tracking method for optical small target on water surface | |
Li et al. | A fast level set algorithm for building roof recognition from high spatial resolution panchromatic images | |
CN104732551A (en) | Level set image segmentation method based on superpixel and graph-cup optimizing | |
CN112132862B (en) | Adaptive scale estimation target tracking algorithm based on unmanned aerial vehicle | |
Glocker et al. | Triangleflow: Optical flow with triangulation-based higher-order likelihoods | |
CN107909597B (en) | Multi-scale MRF model image segmentation method with edge preservation function | |
CN105590327A (en) | Motion estimation method and apparatus | |
Salzmann et al. | Beyond feature points: Structured prediction for monocular non-rigid 3d reconstruction | |
Zhong et al. | Robust image segmentation against complex color distribution | |
CN115619795A (en) | Brain MR and natural image segmentation method based on skewed distribution | |
CN111008977B (en) | Image segmentation method of high-order MRF model based on multi-node topological overlap measure | |
Tsechpenakis et al. | Deformable probability maps: Probabilistic shape and appearance-based object segmentation | |
CN110363777A (en) | A kind of sea image semantic segmentation method based on reducible space constraint mixed model | |
Cheng et al. | A novel cell structure‐based disparity estimation for unsupervised stereo matching | |
Han et al. | Quasi-regular facade structure extraction | |
Ding et al. | Adaptive Range guided Multi-view Depth Estimation with Normal Ranking Loss | |
Lin et al. | Fast Regularity-Constrained Plane Reconstruction |
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 |