CN102163321A - Image segmentation method based on lattice Boltzman model - Google Patents
Image segmentation method based on lattice Boltzman model Download PDFInfo
- Publication number
- CN102163321A CN102163321A CN 201110060920 CN201110060920A CN102163321A CN 102163321 A CN102163321 A CN 102163321A CN 201110060920 CN201110060920 CN 201110060920 CN 201110060920 A CN201110060920 A CN 201110060920A CN 102163321 A CN102163321 A CN 102163321A
- Authority
- CN
- China
- Prior art keywords
- model
- boltzmann
- grid
- dimensional lattice
- 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
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
The invention discloses an image segmentation method based on a lattice Boltzman model, which comprises the following steps: performing the segmentation process to an image based on a microscopic lattice Boltzman model, and then realizing the solution to a macroscopic partial differential equation of the image segmentation so as to realize the automatic segmentation of the image. Compared with a typical image segmentation based on a partial differential equation method, the image segmentation method can realize an iterative computation in large step so as to effectively improve the efficiency of image segmentation with rapid segmentation speed and good segmentation effect.
Description
Technical field
The present invention relates to that a kind of (Lattice Boltzmann Model, image partition method LBM) belongs to image processing field based on grid Boltzmann model.
Technical background
Image segmentation refers to image is divided into the zone of each tool characteristic and proposes the technology and the process of interesting target, and it is a kind of important images analytical technology.Four ten years in the past, Study of Image Segmentation are subjected to the attention of people's height always.Image segmentation is the major issue in Flame Image Process and the analysis, also is a classic problem in the computer vision research.
The Partial Differential Equation method of Flame Image Process developed rapidly from the nineties in 20th century, had demonstrated the performance that is better than the traditional images disposal route
[1]In image processing field based on Partial Differential Equation method, because movable contour model has high efficiency, the simplicity of calculating and is specially adapted to modeling and the advantages such as deformable contour of extraction arbitrary shape, over nearly 20 years, movable contour model has had in rim detection, medical image segmentation and motion tracking to be used and very big development widely.Active contour in the image segmentation or " snake " (snake) model are at first proposed by M. Kass etc., its basic thought image segmentation problem is summed up as minimize a closed curve C(p) " energy " functional, realize that minimizing of " energy " functional promptly obtains short as far as possible and the closed curve of fairing as far as possible.But the major defect that the method exists is that it not only depends on the geometric configuration and the position of curve C, also depends on the parameter p of curve.In order to overcome this defective, V. Caselles, R. Sapiro proposed not contain free parameter geodesic line active contour (GAC) model in 1997, had eliminated the defective that classical snake model relies on free parameter.But being curve, the defective that GAC exists to approach to have in the image dark border of caving in.And for both there not being limbus also to lack the image of obvious textural characteristics, the GAC model is difficult to realize successful cutting apart.Based on this, T. Chan and L. Vese have proposed non-flanged active contour (C-V) model, have realized preferably segmentation effect automatically.
Digital picture is discontinuous, must be by numerical evaluation to obtain the approximate solution of partial differential equation.Many partial differential equation all are to disperse by explicit finite difference.Partial differential equation is found the solution complexity on the one hand, even does not exist; On the other hand, be subjected to the restriction of stability, the algorithm iteration step-length is very little, needs iteration many times just can produce a desired effect, and algorithm implementation procedure efficient is low.Along with developing rapidly of field robotization such as industry, agricultural, medical science, military affairs and intelligent demand, the requirement of image processing techniques is also improved day by day, therefore needing efficiently, implementation algorithm applies to the image real-time processing domain.
Grid Boltzmann model has physical thought clearly, the advantage of simple boundary treatment and fast parallel calculating, and it is the spatial model that disperses, and is particularly suitable for Digital Image Processing.From grid Boltzmann's microvisual model, design grid Boltzmann EVOLUTION EQUATION can be met macroscopical partial differential equation that Flame Image Process requires at last.Therefore grid Boltzmann method rapidly and efficiently provides the realization approach with accuracy for what realize Flame Image Process.
Summary of the invention
The objective of the invention is to the technical matters at the prior art existence, propose a kind of image partition method based on grid Boltzmann model, this method can not only improve the image segmentation quality, and can improve counting yield, is particularly useful for image and handles in real time.
In order to achieve the above object, design of the present invention is:
Image partition method based on grid Boltzmann model of the present invention comprises: set up two-dimensional lattice Boltzmann model, it is made of the computing grid of discretize, and the node of each grid is equivalent to a cellular, and the value of cellular is by the distribution function of particle
With the diffusion vector
Decision.
Finding the solution the EVOLUTION EQUATION that nonlinear diffusion equations realizes that gray level image is cut apart based on the method for grid Boltzmann model is:
Wherein
Be the vector of diffusion,
For iterations is
nIn time, be positioned at
xThe place has speed
The particle density distribution function, formula
Be relaxation factor,
Be the equilibrium state distribution function,
Be the Heaviside function derivative, its variable
uBe distance function,
c 1,
c 2Be respectively to cut apart the inside and outside average gray of evolution curve,
Be constant.
Based on the flow process of the image partition method of grid Boltzmann model as shown in Figure 4, the method that the present invention proposes adopts the variation level diversity method, introduce the Heaviside function, its derivative function is embedded in the relaxation factor of grid Boltzmann EVOLUTION EQUATION, find grid Boltzmann EVOLUTION EQUATION and macroscopical equation corresponding relation to find the solution nonlinear diffusion equations to realize image segmentation.
According to the foregoing invention design, the present invention adopts following technical proposals:
A kind of image partition method based on grid Boltzmann model is characterized in that operation steps is as follows:
(2) use two-dimensional lattice Boltzmann model, the initial balance state function of each action direction in the grid Boltzmann EVOLUTION EQUATION is set
(3) determine the iterations of grid Boltzmann EVOLUTION EQUATION
NAdopt the variation level diversity method, draw the Heaviside function, the initialization distance function
As the initial segmentation curve C;
(5) the inside and outside average gray value of computed segmentation curve
c 1With
c 2:
(6) calculating parameter
Derivative
, the relaxation factor in the traversal image calculation grid Boltzmann EVOLUTION EQUATION
(7) establish
nBe iterations, be according to two-dimensional lattice Boltzmann model modification balanced distribution function
(9) mechanism of calculating two-dimensional lattice Boltzmann model:
(10) upgrade particle distribution function
Obtain embedding the distanceization function
, get its zero level collection and be the new curve of cutting apart;
(11) judge whether to reach iterations
N, if reach
NWhen inferior, then output
uZero level assemble fruit and be segmentation contour; If do not reach
NInferior, then change step (4), repeating step (4)-(10) are up to reaching iteration
NZero level after output is handled behind the number of times is assembled fruit, is last segmentation result.
According to diffusion vector discrete in the two-dimensional lattice Boltzmann model
Direction number
qDifference, can be divided into D2Q5 and D2Q9 two classes to two-dimensional lattice Boltzmann model.
When the two-dimensional lattice Boltzmann's model in the above-mentioned steps (2) is the D2Q5 model, the initial balance state function
For:
When two-dimensional lattice Boltzmann model is the D2Q9 model in the above-mentioned steps (2), the initial balance state function
For:
,
n=1,2…,
N
The present invention compared with prior art, have following conspicuous outstanding substantive distinguishing features and remarkable advantage: above-mentioned image partition method based on grid Boltzmann model can not only be realized cutting apart automatically of image, obtain high-quality image segmentation effect, and guaranteeing to carry out the computing of big step-length under the stable situation of algorithm, thereby improve the efficient of calculating effectively.Being particularly useful for image handles in real time.
Description of drawings
Fig. 1 is two-dimensional lattice Boltzmann's model structure synoptic diagram that the computational grid by discretize constitutes.
Fig. 2 is based on the grid Boltzmann model structure synoptic diagram of D2Q5 among Fig. 1.
Fig. 3 is based on the grid Boltzmann model structure synoptic diagram of D2Q9 among Fig. 1.
Fig. 4 is the process flow diagram of the image partition method based on grid Boltzmann model of the present invention.
Fig. 5 is the former figure of gray level image to be split.
Fig. 6 is the initial curve figure of image local being cut apart setting.
Fig. 7 is based on the image segmentation design sketch of the grid Boltzmann model of D2Q5.
Fig. 8 is the initial curve figure of image overall being cut apart setting.
Fig. 9 is based on the image segmentation design sketch of the grid Boltzmann model of D2Q5.
Embodiment
Embodiment to the image partition method based on grid Boltzmann model of the present invention elaborates below: present embodiment is to implement under the prerequisite with technical scheme of the present invention; provided detailed embodiment, but protection scope of the present invention is not limited to following embodiment.
Embodiments of the invention are described with reference to the accompanying drawings as follows:
Shown in Fig. 1,2 and 3, set up two-dimensional lattice Boltzmann model, it is made of the computing grid of discretize, and the node of each grid is equivalent to a cellular, and its value is by the distribution function of particle
With the diffusion vector
Decision.Each renewal (iteration) of nodal value can be divided into two stages: migration phase and effect stage.Migration phase is transmitted particle by the field node to Centroid, and the effect stage then determines the quantity transmitted.Utilization grid Boltzmann models treated
M*
NDigital picture the time, each pixel and computing node can be mapped naturally, the gray-scale value of pixel is the population on the corresponding node then.Each pixel of node among Fig. 1 (round dot) representative image, arrow has shown the migration and the action direction of model.
According to discrete diffusion vector
Direction number
qDifference, can be two-dimensional lattice Boltzmann D
nQ
q(
nDimension
qDirection) be divided into two classes:
The D2Q5 model, as shown in Figure 2, it has 5 discrete speed:
The D2Q9 model, as shown in Figure 3, it has 9 discrete speed:
Based on the flow process of the image partition method of grid Boltzmann model as shown in Figure 4, the method that the present invention proposes embeds the Heaveside function derivative of image in the relaxation factor of grid Boltzmann EVOLUTION EQUATION, has found grid Boltzmann EVOLUTION EQUATION and macroscopical corresponding relation to find the solution nonlinear diffusion equations to realize image segmentation.
Embodiment 1: image local is cut apart
The segmentation effect of gray level image is shown in Fig. 5,6 and 7, and Fig. 5 is a head mri image to be split, and size is 388 * 388; Fig. 6 is the initial curve figure that image segmentation is provided with, and wherein white marking is initial circular curve; Fig. 7 is the image segmentation design sketch based on the grid Boltzmann model of D2Q5.
With the dividing method based on D2Q5 grid Boltzmann model is example, and its step is as follows:
(2) use two-dimensional lattice Boltzmann model, the initial balance state function of each action direction in the grid Boltzmann EVOLUTION EQUATION is set
:
(3) determine iterations
NBe 80, and initial curve C is set is:
(4) calculate the Heavside function
, parameter wherein
=1.0;
(5) the inside and outside average gray value c of computed segmentation curve
1And c
2:
(6) calculating parameter
Derivative
, the traversal image calculation is based on the relaxation factor in the grid Boltzmann EVOLUTION EQUATION of D2Q5
, wherein
=3;
(7) establish
nBe iterations, be according to two-dimensional lattice Boltzmann model modification balanced distribution function
(9) mechanism of calculating two-dimensional lattice Boltzmann model:
Wherein
=0.004;
(10) upgrade particle distribution function
Obtain embedding the distanceization function
, get its zero level collection and be the new curve of cutting apart;
(11) judge whether to reach iterations 80, if when reaching 80 times, then output
uZero level assemble fruit and be segmentation contour; If do not reach
NInferior, then change step (4), fruit is assembled up to the zero level that reaches after 80 number back output of iteration is handled in repeating step (4)-(10), is last segmentation result such as Fig. 7.At HP notebook (AMD Turion 64,2.01GHz CPU, 2G RAM), running environment is 5.98s for working time on the MATLAB 2008a.
Embodiment 2: image overall is cut apart
As Fig. 8 is the initial curve figure that head mri image segmentation to be split is provided with; Fig. 9 is the image segmentation design sketch based on the grid Boltzmann model of D2Q5.With the dividing method based on D2Q5 grid Boltzmann model is example, and its step is identical with embodiment one, and wherein the parameter setting changes to some extent: the initial curve of setting as shown in Figure 8; Iterations
NIt is 4 times.Last segmentation effect such as Fig. 9, be 1.28s working time.
Experiment shows, based on the image partition method of grid Boltzmann model, by different initial curves is set, can carries out target effectively and cut apart.For global segmentation, can cut apart fast automatically, segmentation result is good.
Claims (7)
1. image partition method based on grid Boltzmann model is characterized in that operation steps is as follows:
Use two-dimensional lattice Boltzmann model, the initial balance state function of each action direction in the grid Boltzmann EVOLUTION EQUATION is set
Determine the iterations of grid Boltzmann EVOLUTION EQUATION
NAdopt the variation level diversity method, introduce the Heaviside function, the initialization distance function
As initial curve C;
The average gray value that the computed segmentation curve is inside and outside
c 1With
c 2:
Calculating parameter
Derivative
, the relaxation factor in the traversal image calculation grid Boltzmann EVOLUTION EQUATION
If
nBe iterations, be according to two-dimensional lattice Boltzmann model modification balanced distribution function
Calculate the transition process of two-dimensional lattice Boltzmann model:
Calculate the mechanism of two-dimensional lattice Boltzmann model:
Upgrade particle distribution function
Obtain embedding the distanceization function
, get its zero level collection and be the new curve of cutting apart;
Judge whether to reach iterations
N, if reach
NWhen inferior, then output
uZero level is assembled fruit and is segmentation contour; If do not reach
NInferior, then change step (4), repeating step (4)-(10) are up to reaching iteration
NZero level after output is handled behind the number of times is assembled fruit, is last segmentation result.
4. the image partition method based on grid Boltzmann model according to claim 1 is characterized in that:
When the two-dimensional lattice Boltzmann's model in the above-mentioned steps (6) is the D2Q5 model, relevant relaxation factor
For:
5. the image partition method based on grid Boltzmann model according to claim 1 is characterized in that:
When two-dimensional lattice Boltzmann model is the D2Q9 model in the above-mentioned steps (6), relevant relaxation factor
For:
Wherein
Be parameter of curve, the smoothness of decision evolution curve.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110060920 CN102163321B (en) | 2011-06-15 | 2011-06-15 | Image segmentation method based on lattice Boltzman model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110060920 CN102163321B (en) | 2011-06-15 | 2011-06-15 | Image segmentation method based on lattice Boltzman model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102163321A true CN102163321A (en) | 2011-08-24 |
CN102163321B CN102163321B (en) | 2013-08-28 |
Family
ID=44464533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110060920 Expired - Fee Related CN102163321B (en) | 2011-06-15 | 2011-06-15 | Image segmentation method based on lattice Boltzman model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102163321B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102737378A (en) * | 2012-05-24 | 2012-10-17 | 中国科学院深圳先进技术研究院 | Method and system for partitioning image by using geodesic active contour |
CN103337097A (en) * | 2013-07-02 | 2013-10-02 | 上海大学 | Multiple Cartesian grid generation method applicable to LBM) |
CN104867110A (en) * | 2014-12-08 | 2015-08-26 | 上海大学 | Lattice Boltzmann model-based video image defect repairing method |
CN105241911A (en) * | 2015-09-23 | 2016-01-13 | 中国石油大学(北京) | LBM-based simulated low field nuclear magnetic resonance fluid analysis method and device |
CN106404730A (en) * | 2016-08-30 | 2017-02-15 | 上海大学 | Description method for propagating light in medium based on lattice Boltzmann model |
CN107240114A (en) * | 2017-05-27 | 2017-10-10 | 南京林业大学 | A kind of semi-automatic medical image cutting method based on a distance function shape constraining |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080107351A1 (en) * | 2006-11-07 | 2008-05-08 | Siemens Corporate Research, Inc. | Efficient segmentation of piecewise smooth images |
CN101673393A (en) * | 2009-09-25 | 2010-03-17 | 上海大学 | Image de-noising method based on lattice Boltzmann model |
-
2011
- 2011-06-15 CN CN 201110060920 patent/CN102163321B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080107351A1 (en) * | 2006-11-07 | 2008-05-08 | Siemens Corporate Research, Inc. | Efficient segmentation of piecewise smooth images |
CN101673393A (en) * | 2009-09-25 | 2010-03-17 | 上海大学 | Image de-noising method based on lattice Boltzmann model |
Non-Patent Citations (2)
Title |
---|
《Proceedings f the 29th Annual International Conference of the IEEE EMBS》 20070826 Yu Chen et al. Application of Lattice Boltzmann Method to Image Segmentation 全文 1-7 , * |
《电子测量技术》 20110331 张蕊等 图像修复的格子波尔兹曼方法 全文 1-7 第31卷, 第3期 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102737378A (en) * | 2012-05-24 | 2012-10-17 | 中国科学院深圳先进技术研究院 | Method and system for partitioning image by using geodesic active contour |
CN102737378B (en) * | 2012-05-24 | 2015-03-11 | 中国科学院深圳先进技术研究院 | Method and system for partitioning image by using geodesic active contour |
CN103337097A (en) * | 2013-07-02 | 2013-10-02 | 上海大学 | Multiple Cartesian grid generation method applicable to LBM) |
CN104867110A (en) * | 2014-12-08 | 2015-08-26 | 上海大学 | Lattice Boltzmann model-based video image defect repairing method |
CN105241911A (en) * | 2015-09-23 | 2016-01-13 | 中国石油大学(北京) | LBM-based simulated low field nuclear magnetic resonance fluid analysis method and device |
CN105241911B (en) * | 2015-09-23 | 2017-07-21 | 中国石油大学(北京) | The method and device that low-field nuclear magnetic resonance analyzes fluid is simulated based on LBM |
CN106404730A (en) * | 2016-08-30 | 2017-02-15 | 上海大学 | Description method for propagating light in medium based on lattice Boltzmann model |
CN106404730B (en) * | 2016-08-30 | 2019-10-11 | 上海大学 | The description method that light based on Lattice Boltzmann method model is propagated in the medium |
CN107240114A (en) * | 2017-05-27 | 2017-10-10 | 南京林业大学 | A kind of semi-automatic medical image cutting method based on a distance function shape constraining |
CN107240114B (en) * | 2017-05-27 | 2019-11-12 | 南京林业大学 | A kind of semi-automatic medical image cutting method based on distance function shape constraining |
Also Published As
Publication number | Publication date |
---|---|
CN102163321B (en) | 2013-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102163321B (en) | Image segmentation method based on lattice Boltzman model | |
CN106685716B (en) | Network topology adaptive data visualization method and device | |
CN107220277A (en) | Image retrieval algorithm based on cartographical sketching | |
CN110992473B (en) | Tree branch modeling method and system based on vehicle-mounted laser scanning point cloud | |
CN101158985B (en) | Ultra-dimension fluvial dynamics self-adapting parallel monitoring method | |
CN105118090B (en) | A kind of point cloud filtering method of adaptive complicated landform structure | |
WO2015149302A1 (en) | Method for rebuilding tree model on the basis of point cloud and data driving | |
CN106022215A (en) | Face feature point positioning method and device | |
CN106651765A (en) | Method for automatically generating thumbnail by use of deep neutral network | |
CN113538689A (en) | Three-dimensional model mesh simplification method based on feature fusion of neural network | |
Doraiswamy et al. | Output-sensitive construction of Reeb graphs | |
CN101436303B (en) | Method for obtaining tetrahedral grid from object three-dimensional image | |
CN109740106A (en) | Large-scale network betweenness approximation method based on graph convolution neural network, storage device and storage medium | |
CN106023286A (en) | Fluid animation accelerated generation method based on data driving | |
CN104036550A (en) | Laser radar point-cloud interpretation and reconstruction method for building elevations on basis of shape semantics | |
CN108961385A (en) | A kind of SLAM patterning process and device | |
CN109948575A (en) | Eyeball dividing method in ultrasound image | |
CN112241676A (en) | Method for automatically identifying terrain sundries | |
CN107341449A (en) | A kind of GMS Calculation of precipitation method based on cloud mass changing features | |
Çavuşoğlu | The position of BIM tools in conceptual design phase: parametric design and energy modeling capabilities | |
CN117237542B (en) | Three-dimensional human body model generation method and device based on text | |
Kamatsuchi | Turbulent flow simulation around complex geometries with cartesian grid method | |
CN107862749A (en) | One kind point cloud geometric detail feature minimizing technology | |
CN109961424B (en) | Hand X-ray image data generation method | |
Röber et al. | In-Situ processing in climate science |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130828 Termination date: 20180615 |
|
CF01 | Termination of patent right due to non-payment of annual fee |