CN106053316A - Program and calculation method for pore distribution and fractal of microstructure of soil mass - Google Patents
Program and calculation method for pore distribution and fractal of microstructure of soil mass Download PDFInfo
- Publication number
- CN106053316A CN106053316A CN201610402716.3A CN201610402716A CN106053316A CN 106053316 A CN106053316 A CN 106053316A CN 201610402716 A CN201610402716 A CN 201610402716A CN 106053316 A CN106053316 A CN 106053316A
- Authority
- CN
- China
- Prior art keywords
- image
- hole
- calculation
- soil
- fractal
- 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.)
- Pending
Links
- 239000002689 soil Substances 0.000 title claims abstract description 47
- 239000011148 porous material Substances 0.000 title claims abstract description 16
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 13
- 239000002245 particle Substances 0.000 claims description 12
- 238000000205 computational method Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000003086 colorant Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
The invention discloses a program and calculation method for pore distribution and fractal of a microstructure of soil mass. The method comprises the following steps: images subjected to cutting, gray processing, contrast control, image binaryzation with appropriate selected threshold value as well as noise reduction processing are copied to a current folder before calculation; a well edited fig file is opened, the size of epsilon is selected in the fig file, then a system pops up a dialog box automatically to prompt selection of an image needing processing, after the image is selected, the system calculates the porosity of the image and the number of meshes occupied by pores instantly, and a calculation result is shown in a pie chart form and a dialog box form. The operation steps are simple and convenient, pixel is taken as a unit, the calculation speed and the calculation precision are high, the calculation efficiency can be increased greatly, and important parameters such as the pore distribution, fractal dimension and the like of the microstructure of the soil mass can be obtained.
Description
Art
The present invention relates to one and be applicable to Soil Microstructure distribution of pores and fractal program and computational methods.
Background technology
Soil Microstructure is the basis that its macroscopic view exists, and the research to Soil Microstructure is broadly divided into qualitative and quantitative
In two stages, the qualitative investigation to Soil Microstructure at present has had the biggest progress, and quantitative study becomes current soil
The research emphasis of body microstructure.Distribution of pores is the critical nature of the soil body, is also the important indicator of reflection soil body microscopic characteristics,
It is thus desirable to a kind of accurate method calculates the distribution of pores of the soil body, on the other hand, the development of fractal theory is advanced by leaps and bounds, soil
Body is to be made up of soil particle, and the character of himself has certain self-similarity, which dictates that the soil body is applicable to fractal theory, profit
Combine with fractal theory by image processing techniques, it is possible to obtain numerous information of Soil Microstructure quantitative study, propose one
It is also abnormal the most important that the approach of kind calculates fractal dimension.
Current technology has the disadvantage that complex operation step；Computational accuracy can not be guaranteed；The quantitative letter obtained
Breath is few, is not suitable for Soil Microstructure analysis.It is thus desirable to design a kind of program being applicable to the soil body and computational methods are counted
The important indicators such as the distribution of pores and the fractal dimension that calculate Soil Microstructure.
Summary of the invention
It is an object of the invention to provide a kind of simplicity easily realize and can quickly calculate be applicable to Soil Microstructure hole
Gap distribution and fractal program and computational methods, can calculate many important parameters such as porosity, fractal dimension by the method.
The technical scheme is that a kind of Soil Microstructure distribution of pores and fractal program and computational methods, the party
Method content is as follows: before calculating, and the soil sample image of the bmp form of scanned Electronic Speculum shooting need to be through cutting out, at gray scale
Reason, setting contrast, the suitable threshold value of selection carry out image binaryzation, image noise reduction early stage processes, and the image carrying out calculating divides
For two kinds of colors of black and white, wherein black be hole, white be soil particle；In units of a pixel, it is assumed that required calculating image
Size be m × n-pixel, the square with the length of side as ε divides the image into the grid into (m/ ε) × (n/ ε), by with MATLAB
It is N (ε) that the program write calculates containing the hole i.e. grid sum of black part, or calculates the grid containing aperture sections
Sum is N (ε), changes the size of ε, then obtains the number of grid shared by corresponding hole, and concrete operation step is: at needs
The binary image of reason copies to, under current file folder, open the fig file editted, select the value of ε, click on this ε value corresponding
Button, eject current bmp FileDialog, click on and need image to be processed, system will quickly calculate result, from meter
Calculate in result it can be seen directly that proportion shared by hole, and the meshes number of image area, hole and soil particle and each
Shared proportion；
And the data of gained are plotted in log-log coordinate system, described log-log coordinate system is with lg ε as abscissa,
With lgN (ε) as vertical coordinate；If the figure in coordinate system has linear character, then curve is carried out linear fit, after matching
Slope is k, then opposite number-the k of slope is fractal dimension；On the other hand, program also can calculate shared by black hole automatically
Number of pixels, thus draw the porosity of the soil body, result of calculation can present with the form of pie chart, can find out hole intuitively
Distribution situation.
Owing to using technique scheme, compared with prior art the present invention has such beneficial effect:
1, the image that pixel is higher can be calculated, applied widely；
2, operating procedure is simple, calculates speed fast, can be greatly improved computational efficiency；
3, in units of pixel, computational accuracy is high；
4, the quantitative information of result of calculation gained is many, it is possible to obtain distribution of pores and the fractal dimension etc. of Soil Microstructure
Important parameter；Fractal dimension can be calculated, it is also possible to from pie chart, find out distribution of pores situation intuitively.
Accompanying drawing explanation
Fig. 1 is fig file schematic diagram；
Fig. 2 by calculating picture black (hole) and white (soil particle) proportion size pie chart；
Fig. 3 by the data knot of number of pixels and proportion shared by calculating picture black (hole) and white (soil particle)
Really schematic diagram；
Fig. 4 is for calculating Soil Microstructure porosity and fractal dimension flow chart.
Detailed description of the invention
A kind of Soil Microstructure distribution of pores that the present invention proposes and fractal program and computational methods, the method can be fast
Speed calculates the multiple quantitative information such as Soil Microstructure porosity and fractal dimension, and the method is simple and convenient, calculates
Speed is fast, it is thus achieved that quantitative information many, computational efficiency can be greatly improved, below in conjunction with the accompanying drawings and example is described in detail:
A kind of Soil Microstructure distribution of pores and fractal program and computational methods, the method concrete operation step is as follows:
Before calculating, the soil sample image of the bmp form of scanned Electronic Speculum shooting need to through cutting out, gray proces,
Setting contrast, select suitable threshold value carry out image binaryzation, image noise reduction early stage process, the image carrying out calculating is divided into black
White two kinds of colors, wherein black is that hole, white are for soil particle；In units of a pixel, it is assumed that required calculating image big
Little for m × n-pixel, the square with the length of side as ε divides the image into the grid into (m/ ε) × (n/ ε), by writing with MATLAB
Program to calculate containing the hole i.e. grid sum of black part be N (ε), or calculate the grid sum containing aperture sections
For N (ε), changing the size of ε, then obtain the number of grid shared by corresponding hole, concrete operation step is: will need to be processed
Binary image copies to, under current file folder, open the fig file editted, as it is shown in figure 1, select the value of ε, click on this ε
Value corresponding button, eject current bmp FileDialog, click on and need image to be processed, its result of calculation by: calculated
The pie chart of picture black (hole) and white (soil particle) proportion size is as in figure 2 it is shown, calculated picture black (hole)
With white (soil particle) shared by number of pixels and proportion data result schematic diagram as shown in Figure 3.Can from result of calculation
Directly to find out the proportion shared by hole, and the meshes number of image area, hole and soil particle and each shared proportion；
When ε=1, i.e. with the square grid division that the length of side is 1 pixel, write corresponding M file, calculate when ε=1
Time hole shared by meshes number N (ε1) and porosity, this porosity calculated can reflect the distribution of pores of the soil body really；
When ε=2, i.e. with the square grid division that the length of side is 2 pixels, write corresponding M file, calculate when ε=2
Time hole shared by meshes number be N (ε2)；
When ε=3, i.e. with the square grid division that the length of side is 3 pixels, write corresponding M file, calculate when ε=3
Time hole shared by meshes number be N (ε3)；
When ε=4, i.e. with the square grid division that the length of side is 4 pixels, write corresponding M file, calculate when ε=4
Time hole shared by meshes number be N (ε4)；
When ε=5, i.e. with the square grid division that the length of side is 5 pixels, write corresponding M file, calculate when ε=5
Time hole shared by number of pixels be N (ε5)；
When ε=6, i.e. with the square grid division that the length of side is 6 pixels, write corresponding M file, when calculating ε=6
Meshes number shared by hole is N (ε6)。
2) editor's fig file as shown in Figure 1, then edit corresponding M file and fig file is reached as follows require: when
When clicking on ε=1 button, calculate N (ε1) and porosity；When clicking on ε=2, calculate N (ε2)；When clicking on ε=3, calculate N
(ε3)；When clicking on ε=4, calculate N (ε4)；When clicking on ε=5, calculate N (ε5)；When clicking on ε=6, calculate N (ε6)。
3) selecting the value of ε, fig file will call corresponding M file and take different ε values to calculate selected image
Grid number N (ε shared by lower black holei), porosity, the shared grid number of white and white ratio shared by soil particle area
Weight, result of calculation can present with the form of pie chart and dialog box, as shown in Figure 2 and Figure 3.
Finally, and the data of gained being plotted in log-log coordinate system, described log-log coordinate system is with lg ε as horizontal stroke
Coordinate, with logN (εi) it is vertical coordinate, if the figure in coordinate system is linear, then curve is carried out linear fit and obtain
The slope k of straight line, then-k is fractal dimension, and calculating Soil Microstructure porosity and fractal dimension flow chart are as shown in Figure 4.
On the other hand, program also can calculate the number of pixels shared by black hole automatically, thus draws the porosity of the soil body, calculates knot
Fruit can present with the form of pie chart, can find out the distribution situation of hole intuitively.The present invention is simply effective, and result of calculation is permissible
From pie chart, find out the distribution situation of soil body hole intuitively, the calculating time can be greatly saved, and accuracy is higher.
Claims (1)
1. a Soil Microstructure distribution of pores and fractal program and computational methods, it is characterised in that: the method content is as follows:
Before calculating, the soil sample image of the bmp form of scanned Electronic Speculum shooting need to through cutting out, gray proces, contrast adjust
Whole, select suitable threshold value carry out image binaryzation, image noise reduction early stage process, the image carrying out calculating is divided into two kinds of face of black and white
Color, wherein black is that hole, white are for soil particle；In units of a pixel, it is assumed that the size of required calculating image is m × n
Pixel, the square with the length of side as ε divides the image into the grid into (m/ ε) × (n/ ε), by the program meter write with MATLAB
Calculating containing the hole i.e. grid sum of black part is N (ε), or to calculate the grid sum containing aperture sections be N (ε),
Changing the size of ε, then obtain the number of grid shared by corresponding hole, concrete operation step is: will need binaryzation to be processed
Under copying image is pressed from both sides to current file, open the fig file editted, select the value of ε, click on the button that this ε value is corresponding, eject
Current bmp FileDialog, clicks on and needs image to be processed, and system will quickly calculate result, can from result of calculation
Directly to find out the proportion shared by hole, and the meshes number of image area, hole and soil particle and each shared proportion；
And the data of gained are plotted in log-log coordinate system, described log-log coordinate system is with lg ε as abscissa, with
LgN (ε) is vertical coordinate；If the figure in coordinate system has linear character, then curve is carried out linear fit, oblique after matching
Rate is k, then opposite number-the k of slope is fractal dimension；On the other hand, program also can calculate shared by black hole automatically
Number of pixels, thus draw the porosity of the soil body, result of calculation can present with the form of pie chart, can find out hole intuitively
Distribution situation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610402716.3A CN106053316A (en) | 2016-06-08 | 2016-06-08 | Program and calculation method for pore distribution and fractal of microstructure of soil mass |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610402716.3A CN106053316A (en) | 2016-06-08 | 2016-06-08 | Program and calculation method for pore distribution and fractal of microstructure of soil mass |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106053316A true CN106053316A (en) | 2016-10-26 |
Family
ID=57169796
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610402716.3A Pending CN106053316A (en) | 2016-06-08 | 2016-06-08 | Program and calculation method for pore distribution and fractal of microstructure of soil mass |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106053316A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106645244A (en) * | 2016-12-26 | 2017-05-10 | 首钢总公司 | Evaluation method for Ti element dispersity in Ti-containing furnace burden |
| CN107240132A (en) * | 2017-05-09 | 2017-10-10 | 中国人民解放军装甲兵工程学院 | A kind of method that utilization gray level co-occurrence matrixes calculate fractal dimension |
| CN108061697A (en) * | 2017-12-06 | 2018-05-22 | 石家庄铁道大学 | Soil body three-dimensional pore space rate computational methods |
| CN108876778A (en) * | 2018-06-15 | 2018-11-23 | 西安建筑科技大学 | A kind of quantitative analysis method on loam wall surface and section crack |
| CN108961246A (en) * | 2018-07-10 | 2018-12-07 | 吉林大学 | A kind of scanning electron microscope image hole recognition methods based on artificial intelligence |
| CN109444015A (en) * | 2018-10-31 | 2019-03-08 | 成都理工大学 | More kens, it is multiple dimensioned under shale reservoir microcellular system identification method |
| CN110411911A (en) * | 2019-08-28 | 2019-11-05 | 四川大学 | A kind of bulk solid geometry composition characteristic evaluation method based on fractal dimension |
| CN113609736A (en) * | 2021-08-09 | 2021-11-05 | 安徽理工大学 | Numerical calculation model construction method based on hole crack structure digital image |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005104002A1 (en) * | 2004-03-31 | 2005-11-03 | Exxonmobil Upstream Research Company | Method for simulating and estimating sandstone properties |
| US20100154514A1 (en) * | 2008-12-18 | 2010-06-24 | Algive Lionnel | Method of determining the evolution of petrophysical properties of a rock during diagenesis |
| CN102841220A (en) * | 2012-07-16 | 2012-12-26 | 天津城市建设学院 | Clay scanning electron micrograph segmentation method based on porosity |
| CN102901697A (en) * | 2012-09-26 | 2013-01-30 | 北京农业信息技术研究中心 | Porosity detection method for soil |
| CN103033456A (en) * | 2012-12-13 | 2013-04-10 | 北京农业信息技术研究中心 | Soil porosity detection method based on SFS (Shape from Shading) algorithm |
| CN104268830A (en) * | 2014-09-11 | 2015-01-07 | 河海大学 | Digital-image-based method for determining permeability coefficient of nonuniform geotechnical material |
| CN104931400A (en) * | 2015-05-29 | 2015-09-23 | 宁夏大学 | Quantitative test and image analysis method of particle material pore fabric |
-
2016
- 2016-06-08 CN CN201610402716.3A patent/CN106053316A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005104002A1 (en) * | 2004-03-31 | 2005-11-03 | Exxonmobil Upstream Research Company | Method for simulating and estimating sandstone properties |
| US20100154514A1 (en) * | 2008-12-18 | 2010-06-24 | Algive Lionnel | Method of determining the evolution of petrophysical properties of a rock during diagenesis |
| CN102841220A (en) * | 2012-07-16 | 2012-12-26 | 天津城市建设学院 | Clay scanning electron micrograph segmentation method based on porosity |
| CN102901697A (en) * | 2012-09-26 | 2013-01-30 | 北京农业信息技术研究中心 | Porosity detection method for soil |
| CN103033456A (en) * | 2012-12-13 | 2013-04-10 | 北京农业信息技术研究中心 | Soil porosity detection method based on SFS (Shape from Shading) algorithm |
| CN104268830A (en) * | 2014-09-11 | 2015-01-07 | 河海大学 | Digital-image-based method for determining permeability coefficient of nonuniform geotechnical material |
| CN104931400A (en) * | 2015-05-29 | 2015-09-23 | 宁夏大学 | Quantitative test and image analysis method of particle material pore fabric |
Non-Patent Citations (1)
| Title |
|---|
| 刘熙媛: "基于分形理论的土体微观结构研究", 《建筑科学》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106645244B (en) * | 2016-12-26 | 2019-10-08 | 首钢集团有限公司 | The evaluation method of Ti element dispersion degree in a kind of titaniferous furnace charge |
| CN106645244A (en) * | 2016-12-26 | 2017-05-10 | 首钢总公司 | Evaluation method for Ti element dispersity in Ti-containing furnace burden |
| CN107240132A (en) * | 2017-05-09 | 2017-10-10 | 中国人民解放军装甲兵工程学院 | A kind of method that utilization gray level co-occurrence matrixes calculate fractal dimension |
| CN107240132B (en) * | 2017-05-09 | 2019-11-12 | 中国人民解放军装甲兵工程学院 | A method of fractal dimension is calculated using gray level co-occurrence matrixes |
| CN108061697B (en) * | 2017-12-06 | 2020-09-25 | 石家庄铁道大学 | Method for calculating three-dimensional porosity of soil body |
| CN108061697A (en) * | 2017-12-06 | 2018-05-22 | 石家庄铁道大学 | Soil body three-dimensional pore space rate computational methods |
| CN108876778A (en) * | 2018-06-15 | 2018-11-23 | 西安建筑科技大学 | A kind of quantitative analysis method on loam wall surface and section crack |
| CN108961246A (en) * | 2018-07-10 | 2018-12-07 | 吉林大学 | A kind of scanning electron microscope image hole recognition methods based on artificial intelligence |
| CN108961246B (en) * | 2018-07-10 | 2021-05-07 | 吉林大学 | Scanning electron microscope image pore identification method based on artificial intelligence |
| CN109444015A (en) * | 2018-10-31 | 2019-03-08 | 成都理工大学 | More kens, it is multiple dimensioned under shale reservoir microcellular system identification method |
| CN110411911A (en) * | 2019-08-28 | 2019-11-05 | 四川大学 | A kind of bulk solid geometry composition characteristic evaluation method based on fractal dimension |
| CN113609736A (en) * | 2021-08-09 | 2021-11-05 | 安徽理工大学 | Numerical calculation model construction method based on hole crack structure digital image |
| CN113609736B (en) * | 2021-08-09 | 2024-05-17 | 安徽理工大学 | Numerical calculation model construction method based on hole fracture structure digital image |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106053316A (en) | Program and calculation method for pore distribution and fractal of microstructure of soil mass | |
| Cai et al. | Measurement and characterization of porosity in aluminium selective laser melting parts using X-ray CT | |
| CN108961303B (en) | Image processing method and device, electronic equipment and computer readable medium | |
| US9811933B2 (en) | Image editing using selective editing tools | |
| US10055884B2 (en) | Three-dimensional fluid micromodels | |
| Cocks et al. | A guide to analysis and reconstruction of serial block face scanning electron microscopy data | |
| CN108682008B (en) | A kind of Leukocyte Image clarity evaluation method and device | |
| CN111488921A (en) | Panoramic digital pathological image intelligent analysis system and method | |
| CN110148088B (en) | Image processing method, image rain removing method, device, terminal and medium | |
| WO2007148284A3 (en) | A method, a system and a computer program for determining a threshold in an image comprising image values | |
| CN111062854B (en) | Method, device, terminal and storage medium for detecting watermark | |
| CN104484410A (en) | Data fusion method and system applied to big data system | |
| WO2009015058A2 (en) | High dynamic range image hallucination | |
| CN111597845A (en) | Two-dimensional code detection method, device and equipment and readable storage medium | |
| EP0204593B1 (en) | Method for regulating a grey scale in digital image reproduction, and apparatus for carrying out this method | |
| CN111724331A (en) | Porous medium image reconstruction method based on generation network | |
| CN110728675A (en) | Pulmonary nodule analysis device, model training method, device and analysis equipment | |
| CN110334693A (en) | A kind of Remote Sensing Target sample generating method towards deep learning | |
| Perring et al. | Using automated digital image analysis to provide quantitative petrographic data on olivine–phyric basalts | |
| CN110738658B (en) | Image quality evaluation method | |
| CN111079528A (en) | Primitive drawing checking method and system based on deep learning | |
| CN109584224B (en) | Method for analyzing and displaying X-ray image of casting | |
| CN110659724A (en) | Target detection convolutional neural network construction method based on target scale range | |
| US10002210B2 (en) | Simulation model generation method for filler mixed material | |
| CN111429440B (en) | Method, system, equipment, device and medium for detecting sufficiency of microscopic pathology image cells |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161026 |