CN102841220B - Clay scanning electron micrograph segmentation method based on porosity - Google Patents
Clay scanning electron micrograph segmentation method based on porosity Download PDFInfo
- Publication number
- CN102841220B CN102841220B CN201210244143.8A CN201210244143A CN102841220B CN 102841220 B CN102841220 B CN 102841220B CN 201210244143 A CN201210244143 A CN 201210244143A CN 102841220 B CN102841220 B CN 102841220B
- Authority
- CN
- China
- Prior art keywords
- sample
- clay
- threshold value
- dry density
- microstructure
- 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.)
- Expired - Fee Related
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention provides a clay scanning electron micrograph segmentation method based on porosity. The clay scanning electron micrograph segmentation method comprises the following steps: adopting a conventional experimental method in soil mechanics to measure dry density of a clay sample and figuring out two thirds power of the dry density; preparing a microstructural sample of the clay sample in a laboratory; adopting a ZD-A3 freeze-dryer to freeze-drying the prepared sample, vacuumizing the sample, scanning the electron microscope when vacuum degree reaches 10-6Pa, observing a microstructure of the sample to obtain a scanning photograph; adopting image analysis software to extract binaryzation data corresponding to two thresholds; comparing parameter features of granule and porosity extracted under the two thresholds and figuring out area ratio of pore bodies; and stopping approaching when the area ratio equals to the two thirds power of the dry density and considering the threshold as picture segmentation threshold corresponding to the dry density. The beneficial effects of the clay scanning electron micrograph segmentation method provided by the invention are as follows: scientific change from three-dimensional parameters to two-dimensional parameters is realized, interference by man-made subjective factors is eliminated, and study for the microstructure of a soil body is more accurate.
Description
Technical field
The invention belongs to civil engineering work subject ground microstructure study field, be specifically related to a kind of clay microstructure parameters based on scanning electron microscope and determine method.
Technical background
The particle of ground and Agglomerate and hole and pore-body are called as the microstructured bodies of ground, the micromechanism that is totally referred to as ground or the microstructure of its size, shape, orientation, arrangement and each other contact and connection relation, it is the fundamental that forms soil.The engineering characteristic of material and residing state are corresponding with the relation between its microstructure.For ground class material, microstructure parameters has not only reflected its forming process, physical property and mechanical property, and is again the reflection of its current state.The microstructure of clay not only forms relevantly with its Material Source, particle, also depends on its stress state and stress history, and determines its mechanics and hydraulic properties.Say in a sense, between native microstructure and its macromechanics characteristic, exist and must be related to.On the one hand, different mechanical process produces material impact to its micromechanism; On the other hand, micromechanism has determined again its downstream mechanical process.In addition, the hydraulic characteristic such as water suction/dehydration feature, moisture holding capacity, infiltration coefficient of soil also depends on its micromechanism.Therefore, the Study on Microstructure of soil and quantification thereof are significant to disclosing different native hydraulics processes.Therefore, be familiar with native microstructure features extremely important to research soil mechanics internal action mechanism, therefore the interrelationship study of ground microstructure quantification and its mechanical property is also subject to generally attention both domestic and external, and has obtained lot of experimental data and theoretical result.
The Microstructure observation means of soil are common mercury experiment, SEM method, the experiment of aperture ratio surface area, scanning electron microscope method, the CT technology of pressing.Due to technology maturation, the reason such as cost is low, experimentation is simple, SEM method becomes that native microstructure appearance is observed and one of main method of research gradually.SEM image can reflect size, form and the distribution characteristics of particle and hole, and research SEM image contributes to understand the soil body origin cause of formation and related physical mechanical property.
But the SEM image of soil is two dimensional gray photo, in data pick-up process, is to distinguish particle and aperture ranges with the relative size of gray scale and threshold value, and further reflect between soil particle and hole be connected and structural.The gray scale of SEM image is determined by factors such as electron beam distance, particulate mineral composition, accelerating potential, coating film thickness, sample size and fixing meanss, even for same observation station, different operating personnel, different scanning electron microscope, operation time period not even together, often obtain different gray-scale values.And the selection of threshold value often shows as subjectivity, randomness, uncertainty and nonuniqueness are stronger.For same analysis domain, because threshold value is different, often can obtain different image segmentation results, and finally obtain different object micro-parameters.So determining of threshold value just becomes vital link in image information leaching process.
In fact, the micromechanism of a certain soil sample is objectively, not different with " the apparent microstructure difference " that division result was reflected of different threshold values.Therefore the measurement that, adopts scientific and reasonable Threshold definite threshold and carry out microstructure parameters seems very meaningful.
Except " with sensation " setting threshold of routine, digital elevation model and the Tobin's mean variance model based on Otsu algorithm based on generalized information system in addition of nearest Introduction of Literatures, but these methods still need artificially to arrange a certain parameter, therefore still do not break away from the impact of subjective factor.
The subject matter existing is at present: in the microstructure evaluation process of ground SEM photo, threshold value determine the impact still do not break away from subjective factor, therefore follow-up image is cut apart, the work such as microstructure parameters extractions, parameter quantitative research and thus obtained all results and conclusion be all subjective and not objective.Therefore, that can not exaggerate says, science, rationally, objectively Threshold is prerequisite and the basis of carrying out native microstructure study, to from microstructure angle assessment ground aspect terms of mechanics, hydraulics and the otherwise otherness such as thermodynamics, electromagnetics have very important significance.
Summary of the invention
For solving the problem existing in above-mentioned technology, the object of the invention is to propose a kind of clay microstructure parameters based on scanning electron microscope and determine method, reach more objective, reasonable to native Study on Microstructure, the object that reliability is strong.
For achieving the above object, the technical solution used in the present invention is to provide a kind of clay stereoscan photograph image partition method based on porosity, and the method comprises the following steps:
(1) adopt soil mechanics normal experiment method to record the dry density of clay sample, three/second power of calculating dry density is standby;
(2) the microstructure sample of clay sample is prepared in laboratory;
(3) adopt ZD-A3 freeze-drying instrument by the sample freeze-drying of step (2) making and vacuumize, vacuum tightness reaches to 10
-6pa;
(4) sample of processing through step (3) is carried out to scanning electron microscope, observe the microstructure of sample and obtain its SEM photo;
(5) corresponding binaryzation data while adopting Leica QWin image analysis software to extract threshold value and be respectively a and b SEM photo, a+10=b wherein, 0≤a≤255, a is artificial definite primary election threshold value arbitrarily, b is for selecting threshold value according to another definite time of " a+10=b " principle;
(6) the binaryzation data that obtain according to (5) are calculated under a, two threshold conditions of b soil body pore-body Area Ratio ξ in SEM picture
awith ξ
b, ξ wherein
athe area of the SEM picture mesoporosity while being a for threshold value and the Area Ratio of particle, ξ
bthe area of the SEM picture mesoporosity while being b for threshold value and the Area Ratio of particle;
(7) compare ξ
a, ξ
bwith the size of three/second power of the described dry density of step (1), determine that follow-up threshold value adjustment is to raise or lower, and raises with ξ
bfor basic point, 5 gray-scale values of each increase, lower with ξ later
afor basic point, reduce by 5 gray-scale values later at every turn, so constantly change threshold value value, until pore-body Area Ratio is while equaling three/second power big or small of the described dry density of step (1), stop approaching, threshold value is now the image segmentation threshold corresponding with surveyed dry density, and the method that obtains this threshold value is the clay stereoscan photograph image partition method based on porosity.
Effect of the present invention is: 1. adopt this clay microstructure parameters to determine method, by three-dimensional parameter two dimensionization, present stage, this similar conversion method belonged to the first in ground field; 2. determining of image segmentation threshold objectifies more, got rid of the interference of artificial subjective factor, and the research that makes Soil microstructure aspect is accuracy more, and the result precision that utilizes the method to analyze is higher than other method, and result approaches reality more.
Accompanying drawing explanation
Fig. 1 is the corresponding dry density 1.39g/cm of the inventive method
3the effect of 50kPa consolidation pressure complete after SEM figure;
Fig. 2 is the binary picture of the corresponding threshold value of the inventive method Fig. 1 that is 112;
Fig. 3 is the binary picture of the corresponding threshold value of the inventive method Fig. 1 that is 92.
Embodiment
By reference to the accompanying drawings a kind of clay stereoscan photograph image partition method implementation procedure based on porosity of the present invention is illustrated.
The mechanism of action of the present invention is: the dry density of soil is a kind of the measuring of its granule and pore-body distribution proportion, and in the expression soil body, soil particle occupies the ratio in whole analysis domain space.Soil stereoscan photograph come from real particle and distribution of pores state, after image is cut apart, stereoscan photograph in whole analysis domain with regard to discrete granule and pore-body two parts of turning to.The ratio of pore-body area and whole analysis domain area is defined as pore-body Area Ratio, and it represents the density degree of plane of scanning motion granule.Therefore, theoretically, pore-body Area Ratio should be able to reflect the density degree of the original soil body, and has positive connection with indexs such as conventional void ratio, porositys.The concepts such as conventional void ratio, porosity are based on three-dimensional, and stereoscan photograph derives from two-dimensional scan plane.Therefore, the pore-body Area Ratio based on stereoscan photograph should equal three/second power of conventional porosity, and now the critical gray scale of pore-body and granule is the threshold value based on conventional void ratio.
A kind of clay stereoscan photograph image partition method based on porosity of the present invention, the clay stereoscan photograph image segmentation threshold based on porosity is determined method, the method be take the porosity of sample and is found the threshold value that can reflect soil body compaction rate as benchmark, take this threshold value carries out cutting apart of granule and pore-body to two dimensional gray photo as critical point, binary image is carried out to information extraction, adopt on this basis multidimensional analysis method to assess the microstructure parameters of granule and pore-body, so far just completed the whole assessment of native microstructure.
A kind of clay stereoscan photograph image partition method based on porosity of the present invention includes following steps:
(1) adopt soil mechanics normal experiment method to record the dry density of clay sample, three/second power of calculating dry density is standby;
(2) the microstructure sample of clay sample is prepared in laboratory;
(3) adopt ZD-A3 freeze-drying instrument by the sample freeze-drying of step (2) making and vacuumize, vacuum tightness reaches to 10
-6pa;
(4) sample of processing through step (3) is carried out to scanning electron microscope, observe the microstructure of sample and obtain its SEM photo;
(5) corresponding binaryzation data while adopting Leica QWin image analysis software to extract threshold value and be respectively a and b SEM photo, a+10=b wherein, 0≤a≤255, a is primary election threshold value, artificial definite arbitrarily, b is for selecting threshold value according to another definite time of " a+10=b " principle;
(6) the binaryzation data that obtain according to (5) are calculated under a, two threshold conditions of b soil body pore-body Area Ratio ξ in SEM picture
awith ξ
b, ξ wherein
athe area of the SEM picture mesoporosity while being a for threshold value and the Area Ratio of particle, ξ
bthe area of the SEM picture mesoporosity while being b for threshold value and the Area Ratio of particle;
(7) compare ξ
a, ξ
bwith the size of three/second power of the described dry density of step (1), determine that follow-up threshold value adjustment is to raise or lower.Raise with ξ
bfor basic point, 5 gray-scale values of each increase, lower with ξ later
afor basic point, reduce by 5 gray-scale values later at every turn, so constantly change threshold value value, until pore-body Area Ratio is while equaling three/second power big or small of the described dry density of step (1), stop approaching, threshold value is now the image segmentation threshold corresponding with surveyed dry density, and the method that obtains this threshold value is the clay stereoscan photograph image partition method based on porosity.
Below in conjunction with accompanying drawing, a kind of clay stereoscan photograph image partition method implementation procedure based on porosity of the present invention is described in detail.
The original SEM Fig. 1 of take is example with its corresponding binary picture 2 and Fig. 3:
(1) adopting soil mechanics normal experiment method to record clay sample is that the dry density of a1 is 1.39g/cm
3, three/second power of calculating dry density is 1.25g/cm
3;
(2) the microstructure sample of clay sample is prepared in laboratory;
(3) adopt ZD-A3 freeze-drying instrument by the sample freeze-drying of step (2) making and vacuumize, vacuum tightness reaches to 10
-6pa;
(4) sample of processing through step (3) is carried out to scanning electron microscope, observe the microstructure of sample and obtain the SEM photo that clay sample is Fig. 1;
(5) adopt Leica QWin image analysis software to extract threshold value to SEM photo and be respectively 82 and 92 o'clock corresponding binaryzation data.Threshold value is 82 o'clock, and its binary image is Fig. 2, and in Fig. 2, the pixel number that gray scale is 1 is 37619, and the pixel number that gray scale is 0 is for it than being respectively 35362; Threshold value is 92 o'clock, and its binary image is Fig. 3, and in Fig. 3, the number of the pixel that gray scale is 1 is 33207, and the pixel number that gray scale is 0 is 39774; 0.94,0.97;
(6) the binaryzation data that obtain according to (5) are calculated under 82,92 two threshold conditions soil body pore-body Area Ratio 0.94 and 0.97 in SEM picture;
(7) size of the three/second power 1.25 of comparison 0.94,0.97 and the described dry density of step (1), determines that follow-up threshold value raises.Raise that to take the threshold value 92 of 0.97 correspondence be basic point, increase by 5 gray-scale values, so constantly change threshold value value, when threshold value adjusts to 112, pore-body Area Ratio equals 1.25, and the equal and opposite in direction with the three/second power of (1) described dry density, stops approaching, threshold value 112 is now the image segmentation threshold corresponding with surveyed dry density, and the method that obtains this threshold value is the clay stereoscan photograph image partition method based on porosity.
A kind of clay stereoscan photograph image partition method based on porosity of the present invention, by 3 D stereo parameter of soil body porosity, electron-microscope scanning image is two dimensional surface parameter, by utilizing the operation rule of three/second power, three-dimensional amount is rationally transformed into two dimension amount, has realized the reasonable transformation of three-dimensional parameter research with two-dimensional parameter research.Thereby the acquisition that makes the porosity of clay microstructure is the experimental data of objective, utilize above-mentioned operation rule to change choosing of rear guide image threshold value the experimental data of this objective, there is objective rationality, got rid of the interference of artificial subjective factor.
A kind of clay stereoscan photograph image partition method based on porosity of the present invention, gets up the segmentation threshold organic connections of native dry density and its microstructure image, makes the granule that obtains after cutting apart and the pore-body objectivity that more becomes.
Claims (1)
1. the clay stereoscan photograph image partition method based on porosity, the method adopts soil mechanics normal experiment method to record the dry density of clay sample, and three/second power of calculating dry density is standby; The microstructure sample of preparing clay sample through laboratory; Adopt ZD-A3 freeze-drying instrument by the sample freeze-drying of making and vacuumize, vacuum tightness reaches to 10
-6pa; To prepare the sample that clay sample pretreating crosses through laboratory, carry out scanning electron microscope, the microstructure of observing sample also obtains its SEM photo; Also include following steps simultaneously:
(1) adopt Leica QWin image analysis software to extract threshold value to SEM photo, corresponding binaryzation data when threshold value is respectively a and b, a+10=b wherein, 0≤a≤255, a is artificial definite primary election threshold value arbitrarily, and b is for selecting threshold value according to another definite time of " a+10=b " principle;
(2) the binaryzation data that obtain according to step (1) are calculated under a, two threshold conditions of b soil body pore-body Area Ratio ξ in SEM picture
awith ξ
b, ξ wherein
athe area of the SEM picture mesoporosity while being a for threshold value and the Area Ratio of particle, ξ
bthe area of the SEM picture mesoporosity while being b for threshold value and the Area Ratio of particle;
(3) more described ξ
a, ξ
bwith the size of three/second power of described dry density, determine that follow-up threshold value adjustment is to raise or lower, and raises with ξ
bfor basic point, increase by 5 gray-scale values later at every turn; Lower with ξ
afor basic point, reduce by 5 gray-scale values later at every turn, so constantly change threshold value value, until pore-body Area Ratio is while equaling three/second power big or small of described dry density, stop approaching, threshold value is now the image segmentation threshold corresponding with surveyed dry density, and the method that obtains this threshold value is the clay stereoscan photograph image partition method based on porosity.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210244143.8A CN102841220B (en) | 2012-07-16 | 2012-07-16 | Clay scanning electron micrograph segmentation method based on porosity |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210244143.8A CN102841220B (en) | 2012-07-16 | 2012-07-16 | Clay scanning electron micrograph segmentation method based on porosity |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102841220A CN102841220A (en) | 2012-12-26 |
CN102841220B true CN102841220B (en) | 2014-10-01 |
Family
ID=47368768
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210244143.8A Expired - Fee Related CN102841220B (en) | 2012-07-16 | 2012-07-16 | Clay scanning electron micrograph segmentation method based on porosity |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102841220B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103471974B (en) * | 2013-09-05 | 2015-08-12 | 东华大学 | A kind of image method measures the method for fabric theoretic porosity |
CN106053316A (en) * | 2016-06-08 | 2016-10-26 | 燕山大学 | Program and calculation method for pore distribution and fractal of microstructure of soil mass |
CN107462744A (en) * | 2017-08-28 | 2017-12-12 | 中国矿业大学(北京) | A kind of pore space anisotropy evaluation method based on electron-microscope scanning |
CN108061697B (en) * | 2017-12-06 | 2020-09-25 | 石家庄铁道大学 | Method for calculating three-dimensional porosity of soil body |
CN108225989A (en) * | 2017-12-20 | 2018-06-29 | 中国电建集团河南省电力勘测设计院有限公司 | A kind of great soil group quickly names method and storage medium |
CN110458816B (en) * | 2019-08-06 | 2021-04-27 | 北京工商大学 | Fiber material porosity analysis method based on threshold value regression |
CN114034621B (en) * | 2021-11-09 | 2023-07-07 | 中国电建集团成都勘测设计研究院有限公司 | Soil body permeability coefficient acquisition method and device and model construction method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002168756A (en) * | 2000-12-01 | 2002-06-14 | Jeol Ltd | Surface irregularity measuring device |
JP3435684B2 (en) * | 1997-09-30 | 2003-08-11 | 株式会社アドバンテスト | Image information processing device |
WO2011068508A1 (en) * | 2009-12-02 | 2011-06-09 | Hewlett-Packard Development Company, Lp | System and method of foreground-background segmentation of digitized images |
CN102435628A (en) * | 2011-11-02 | 2012-05-02 | 电子科技大学 | Characterization method for structural morphology of porous electrode material |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7171047B2 (en) * | 2002-12-20 | 2007-01-30 | Lsi Logic Corporation | Adaptive Sem edge recognition algorithm |
-
2012
- 2012-07-16 CN CN201210244143.8A patent/CN102841220B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3435684B2 (en) * | 1997-09-30 | 2003-08-11 | 株式会社アドバンテスト | Image information processing device |
JP2002168756A (en) * | 2000-12-01 | 2002-06-14 | Jeol Ltd | Surface irregularity measuring device |
WO2011068508A1 (en) * | 2009-12-02 | 2011-06-09 | Hewlett-Packard Development Company, Lp | System and method of foreground-background segmentation of digitized images |
CN102435628A (en) * | 2011-11-02 | 2012-05-02 | 电子科技大学 | Characterization method for structural morphology of porous electrode material |
Non-Patent Citations (3)
Title |
---|
冻土力学性质影响因素的显著性和交互作用研究;李顺群等;《岩土力学》;20120430;第33卷(第4期);1173-1177 * |
李顺群等.冻土力学性质影响因素的显著性和交互作用研究.《岩土力学》.2012,第33卷(第4期),1173-1177. |
李顺群等.阈值选取对黏土微结构参数的影响.《解放军理工大学学报(自然科学版)》.2011,第12卷(第4期),354-360. * |
Also Published As
Publication number | Publication date |
---|---|
CN102841220A (en) | 2012-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102841220B (en) | Clay scanning electron micrograph segmentation method based on porosity | |
Zou et al. | Two-dimensional SEM image-based analysis of coal porosity and its pore structure | |
Xu et al. | An investigation into the relationship between saturated permeability and microstructure of remolded loess: A case study from Chinese Loess Plateau | |
Lanaro et al. | 3-D-laser measurements and representation of roughness of rock fractures | |
Reingruber et al. | Quantitative characterization of microfiltration membranes by 3D reconstruction | |
CN110021030A (en) | A kind of segmentation threshold of material of rock and soil digital picture determines method | |
CN104501735A (en) | Method for observing three-dimensional deformation of side slope by utilizing circular marking points | |
Wang et al. | Measurement of the deformation of sand in a plane strain compression experiment using incremental digital image correlation | |
Wang et al. | Morphological approach to quantifying soil cracks: Application to dynamic crack patterns during wetting‐drying cycles | |
Zhou et al. | A novel 3D geometrical reconstruction model for porous rocks | |
CN108061697B (en) | Method for calculating three-dimensional porosity of soil body | |
Fomin et al. | Soil pore structure dynamics under steady-state wetting-drying cycle | |
Levitz et al. | Topology of evolving pore networks | |
Liu et al. | Magnification effects on the interpretation of SEM images of expansive soils | |
Yin et al. | Stochastic reconstruction of Gosford sandstone from surface image | |
Li et al. | Microstructure characteristics of nonwoven geotextiles using SEM and CT methods | |
CN110599477B (en) | Method and system for detecting effective porosity of grass planting concrete | |
CN112730102A (en) | Method for measuring evolution of particle material shear band based on standard deviation ellipse | |
CN108875260B (en) | Thermal spraying coating performance parameter calculation method | |
Jian et al. | Surface roughness measure based on average texture cycle | |
Gong et al. | Characterization of micro-scale pore structure and permeability simulation of peat soil based on 2D/3D X-ray computed tomography images | |
Shin et al. | Porosity and pore-size distribution of geomaterials from X-ray CT scans | |
CN102967694B (en) | Microstructure parameter-based clay consolidation coefficient prediction method | |
CN114689820A (en) | Expansive soil swelling and shrinking characteristic testing and evaluating method based on three-dimensional scanning | |
CN114049471A (en) | Hole filling method based on ordered layered simplification scattered point cloud |
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: 20141001 Termination date: 20150716 |
|
EXPY | Termination of patent right or utility model |