CN103700102A - Rock core target extracting method based on CT (Computed Tomography) images - Google Patents
Rock core target extracting method based on CT (Computed Tomography) images Download PDFInfo
- Publication number
- CN103700102A CN103700102A CN201310706021.0A CN201310706021A CN103700102A CN 103700102 A CN103700102 A CN 103700102A CN 201310706021 A CN201310706021 A CN 201310706021A CN 103700102 A CN103700102 A CN 103700102A
- Authority
- CN
- China
- Prior art keywords
- image
- pixel
- rock core
- gray
- scale value
- 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
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a rock core target extracting method based on a CT (Computed Tomography) image. The method mainly comprises the following steps: reading the CT image of a rock core; performing threshold conversion on the CT image of the rock core; extracting an outline of the CT image of the rock core; acquiring a connected region with the maxima pixel number of the CT image of the rock core; acquiring the outline of a rock core target; acquiring a rock core target image. The rock core target extracting method based on the CT image is used for extracting the internal information and outline of the rock core target completely based on the natural forms of the rock core target, so that the rock core target extracting method based on the CT image can be used for accurately acquiring the rock core target.
Description
Technical field
The invention belongs to technical field of image processing, be specifically related to a kind of rock core target extraction method based on CT image.
Background technology
Computed tomography (Computed Tomography, CT) can under harmless condition, to complex object, detect, realization is to the qualitative analysis of the information such as testee inner structure, Density Distribution and size and quantitatively detection, and not affected by measurand external and internal compositions complexity, therefore be widely applied to the fields such as machinery, ground and agriculture and forestry.On ground field, utilize CT to detect, its object is exactly that the fine structure of rock interior is explored and analyzed.Compare as thin slice authenticate technology, scanning electron microscope etc. with conventional art, the advantage of CT technology is under the form that does not change rock interior, structural condition, obtains rapidly and accurately testing result.But the current core analysis technology based on CT image is not utilized whole CT view data of rock core sample conventionally completely, but adopts the Local C T view data in self defined area, thereby has certain error.When particularly verifying with traditional physical method, can cause the inconsistent problem of sample.
Summary of the invention
The object of the invention is to overcome the problems referred to above of the prior art, a kind of rock core target extraction method based on CT image that can Obtaining Accurate rock core target is provided.
For solving the problems of the technologies described above, the present invention by the following technical solutions:
A rock core target extraction method based on CT image, comprises the following steps:
S1. read the CT image of rock core, and gained CT image is designated as to I
0: f (x, y)=g, wherein, (x, y) is image pixel coordinate, g is the gray-scale value that image pixel coordinate (x, y) is corresponding, g ∈ [0,1];
S2. the threshold transformation of rock core CT image: specifically comprise the following steps:
The threshold value of S21. establishing gray-scale value conversion is 0;
S22. travel through the pixel of whole Zhang Yanxin CT image, if the gray-scale value of pixel is 0, keep its gray-scale value constant, if the gray-scale value of pixel is greater than 0, its gray-scale value is set to 1, and corresponding expression formula is
The result of gained bianry image is designated as I
1;
S3. to image I
1carry out profile extraction;
S4. each profile of S3 gained is filled, the gray-scale value that is about to the pixel of each profile inside is made as 1, and obtaining image is I
2: f
2(x, y);
S5. each connected region of S4 gained is carried out to number of pixels statistics, record the number of pixels that each connected region comprises, and statistics is compared, obtain the connected region with maximum pixel number;
S6. by being all set to zero except thering is the gray-scale value of the pixel of the connected region the connected region of maximum pixel number, obtain image I
3: f
3(x, y);
S7. to image I
3carry out profile extraction, obtain rock core objective contour and be labeled as I
4: f
4(x, y);
S8. by image I
3gray-scale value and the original image I of each pixel
0the gray-scale value of the pixel of middle relevant position multiplies each other, and obtains the rock core target image based on CT image, is designated as I
5: f
5(x, y)=f
0(x, y) * f
3(x, y);
S9. by target image I
5with contour images I
4add computing, obtain rock core target and rock core objective contour based on CT image.
Further technical scheme is, described profile extracts and specifically comprises the following steps:
S31. travel through I
1each pixel, the pixel that is 1 for gray-scale value, if the gray-scale value of eight neighborhood points of this pixel is all 1, the gray-scale value of this pixel is set to zero, can extract I
1in all profiles;
S32. according to from left to right, from top to bottom along successively S31 gained profile being carried out to eight neighborhood point trackings, obtain uniformizable profile.
Further technical scheme is, described eight neighborhood point tracking processes are specific as follows: finding first gray-scale value is 1 pixel, be assumed to be the frontier point of lower left, this point is designated as to the first frontier point, in the right side of this pixel, upper right, upper, four of upper lefts neighborhood, have at least one to be frontier point, be designated as the second boundary point, from the second boundary point, by the right side, upper right, upper, upper left, a left side, lower-left, order lower and bottom right, look for the frontier point in consecutive point, be designated as the 3rd frontier point, if the 3rd frontier point is exactly the first frontier point, show that closed contour followed the trail of; Otherwise, from the 3rd frontier point, continue, until find the first frontier point.
Compared with prior art, the invention has the beneficial effects as follows: the rock core target extraction method based on CT image of the present invention can not cause any interference to the irregular contour of rock core target internal information and rock core target, natural form based on rock core target is extracted the internal information of rock core target and profile completely, therefore, adopt the rock core target extraction method based on CT image of the present invention can Obtaining Accurate rock core target.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of the rock core target extraction method based on CT image of the present invention;
Fig. 2 is the rock core CT image in the embodiment of the present invention;
Fig. 3 is the bianry image of the rock core CT image that extracts in the embodiment of the present invention;
Fig. 4 is the position view of 8 neighborhood points in the embodiment of the present invention;
Fig. 5 is that the profile of the bianry image of Fig. 3 extracts result;
Fig. 6 is for filling the image of gained to the profile of Fig. 5;
Fig. 7 is the largest connected region in Fig. 6;
Fig. 8 is the objective contour of the rock core CT image in the embodiment of the present invention;
Fig. 9 is the target image of the rock core CT image in the embodiment of the present invention;
Figure 10 is the objective contour of the rock core CT image in the embodiment of the present invention and the general image of target image.
Embodiment
In order to make object of the present invention, technical scheme and advantage clearer, below in conjunction with drawings and Examples, the present invention is further elaborated.Should be appreciated that specific embodiment described herein, only in order to explain the present invention, is not intended to limit the present invention.
As described in Figure 1, the rock core target extraction method based on CT image in the present embodiment, mainly comprises the following steps:
S1. read the CT image of rock core, gained rock core CT image as shown in Figure 2, is designated as I by this CT image
0: f (x, y)=g, wherein, (x, y) is image pixel coordinate, g is the gray-scale value that image pixel coordinate (x, y) is corresponding, g ∈ [0,1];
S2. the threshold transformation of rock core CT image pixel: specifically comprise the following steps:
The threshold value of S21. establishing gray-scale value conversion is 0;
S22. travel through the pixel of whole Zhang Yanxin CT image, if the gray-scale value of pixel is 0, keep its gray-scale value constant, if the gray-scale value of pixel is greater than 0, be set to 1, corresponding expression formula is
The result of gained bianry image is designated as I
1, as shown in Figure 3;
S3. to image I
1carry out profile extraction, specifically comprise the following steps:
S31. travel through image I
1each pixel, the pixel that is 1 for gray-scale value, if the gray-scale value of eight neighborhood points (position view of 8 field points as shown in Figure 4) of this pixel is all 1, the gray-scale value of this pixel is set to zero, can extract image I
1in all profiles;
S32. according to from left to right, from top to bottom along successively S31 gained profile being carried out to eight neighborhood point trackings, obtain uniformizable profile, as shown in Figure 5; Wherein, the detailed process of following the trail of is as follows: finding first gray-scale value is 1 pixel, be assumed to be the frontier point of lower left, this point is designated as to the first frontier point, in the right side of this pixel, upper right, upper, four of upper lefts neighborhood, have at least one to be frontier point, be designated as the second boundary point, from the second boundary point, by the right side, upper right, upper, upper left, a left side, lower-left, order lower and bottom right, look for the frontier point in consecutive point, be designated as the 3rd frontier point, if the 3rd frontier point is exactly the first frontier point, show that closed contour followed the trail of; Otherwise, from the 3rd frontier point, continue, until find the first frontier point.
S4. each profile of S3 gained is filled, the grey scale pixel value that is about to each profile inside is made as 1, and obtaining image is I
2: f
2(x, y), as shown in Figure 6;
S5. each connected region of S4 gained is carried out to number of pixels statistics, record the number of pixels that each connected region comprises, and statistics is compared, obtain the connected region with maximum pixel number;
S6. by being all set to zero except thering is the gray-scale value of the pixel of the connected region the connected region of maximum pixel number, obtain image I
3: f
3(x, y), as shown in Figure 7;
S7. to image I
3carry out profile extraction, obtain objective contour and be labeled as I
4: f
4(x, y), as shown in Figure 8;
S8. by image I
3gray-scale value and the original image I of each pixel
0the gray-scale value of the pixel of middle relevant position multiplies each other, and obtains the rock core target image based on CT image, is designated as I
5: f
5(x, y)=f
0(x, y) * f
3(x, y), as shown in Figure 9;
S9. by target image I
5with contour images I
4add computing, obtain rock core target and rock core objective contour based on CT image, as shown in figure 10.
Because the rock core target extraction method based on CT image of the present embodiment is the CT image based on rock core target, irregular contour to rock core target internal information and rock core target can not cause any interference, the complete natural form based on rock core target, therefore, the rock core target extraction method based on CT image of the present embodiment can Obtaining Accurate rock core target.
Those of ordinary skill in the art will appreciate that, embodiment described here is in order to help reader understanding's principle of the present invention, should be understood to that protection scope of the present invention is not limited to such special statement and embodiment.Those of ordinary skill in the art can make various other various concrete distortion and combinations that do not depart from essence of the present invention according to these technology enlightenments disclosed by the invention, and these distortion and combination are still in protection scope of the present invention.
Claims (3)
1. the rock core target extraction method based on CT image, is characterized in that: comprise the following steps:
S1. read the CT image of rock core, and gained CT image is designated as to I
0: f (x, y)=g, wherein, (x, y) is image pixel coordinate, g is the gray-scale value that image pixel coordinate (x, y) is corresponding, g ∈ [0,1];
S2. the threshold transformation of rock core CT image: specifically comprise the following steps:
The threshold value of S21. establishing gray-scale value conversion is 0;
S22. travel through the pixel of whole Zhang Yanxin CT image, if the gray-scale value of pixel is 0, keep its gray-scale value constant, if the gray-scale value of pixel is greater than 0, its gray-scale value is set to 1, and corresponding expression formula is
The result of gained bianry image is designated as I
1;
S3. to image I
1carry out profile extraction;
S4. each profile of S3 gained is filled, the gray-scale value that is about to the pixel of each profile inside is made as 1, and obtaining image is I
2: f
2(x, y);
S5. each connected region of S4 gained is carried out to number of pixels statistics, record the number of pixels that each connected region comprises, and statistics is compared, obtain the connected region with maximum pixel number;
S6. by being all set to zero except thering is the gray-scale value of the pixel of the connected region the connected region of maximum pixel number, obtain image I
3: f
3(x, y);
S7. to image I
3carry out profile extraction, obtain objective contour and be labeled as I
4: f
4(x, y);
S8. by image I
3gray-scale value and the original image I of each pixel
0the gray-scale value of the pixel of middle relevant position multiplies each other, and obtains the rock core target image based on CT image, is designated as I
5: f
5(x, y)=f
0(x, y) * f
3(x, y);
S9. by target image I
5with contour images I
4add computing, obtain rock core target and rock core objective contour based on CT image.
2. the rock core target extraction method based on CT image according to claim 1, is characterized in that: the profile described in S3 extracts and specifically comprises the following steps:
S31. travel through I
1each pixel, the pixel that is 1 for gray-scale value, if the gray-scale value of eight neighborhood points of this pixel is all 1, the gray-scale value of this pixel is set to zero, can extract I
1in all profiles;
S32. according to from left to right, from top to bottom along successively S31 gained profile being carried out to eight neighborhood point trackings, obtain uniformizable profile.
3. the rock core target extraction method based on CT image according to claim 2, it is characterized in that: described eight neighborhood point trackings are specific as follows: finding first gray-scale value is 1 pixel, be assumed to be the frontier point of lower left, this point is designated as to the first frontier point, the right side of this pixel, upper right, on, in the neighborhood of four of upper lefts, have at least one to be frontier point, be designated as the second boundary point, from the second boundary point, by right, upper right, on, upper left, left, lower-left, order lower and bottom right is looked for the frontier point in consecutive point, be designated as the 3rd frontier point, if the 3rd frontier point is exactly the first frontier point, show that closed contour followed the trail of, otherwise, from the 3rd frontier point, continue, until find the first frontier point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310706021.0A CN103700102A (en) | 2013-12-20 | 2013-12-20 | Rock core target extracting method based on CT (Computed Tomography) images |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310706021.0A CN103700102A (en) | 2013-12-20 | 2013-12-20 | Rock core target extracting method based on CT (Computed Tomography) images |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103700102A true CN103700102A (en) | 2014-04-02 |
Family
ID=50361621
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310706021.0A Pending CN103700102A (en) | 2013-12-20 | 2013-12-20 | Rock core target extracting method based on CT (Computed Tomography) images |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103700102A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105654525A (en) * | 2014-11-10 | 2016-06-08 | 中国石油化工股份有限公司 | Two-dimensional model construction method for reflecting real rock core connectivity |
| CN107194304A (en) * | 2016-03-15 | 2017-09-22 | 中国石油化工股份有限公司 | The connected region recognition methods of 3-dimensional digital rock core and device |
| CN107204003A (en) * | 2016-03-15 | 2017-09-26 | 中国石油化工股份有限公司 | The connected region recognition methods of two-dimensional digital rock core and device |
| CN110490890A (en) * | 2019-08-21 | 2019-11-22 | 青岛科技大学 | A method of quality digital core image processing is carried out with neural network |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090316219A1 (en) * | 2008-06-18 | 2009-12-24 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method and computer-readable storage medium |
| CN101867790A (en) * | 2010-04-23 | 2010-10-20 | 刘文萍 | Millimeter-wave image analysis method, fire monitoring method and system |
| CN101950352A (en) * | 2010-05-31 | 2011-01-19 | 北京智安邦科技有限公司 | Target detection method capable of removing illumination influence and device thereof |
| CN103235939A (en) * | 2013-05-08 | 2013-08-07 | 哈尔滨工业大学 | Datum point positioning method based on machine vision |
| CN103400134A (en) * | 2013-06-28 | 2013-11-20 | 哈尔滨工程大学 | Non-contact method for extracting region of interest of finger vein sample |
-
2013
- 2013-12-20 CN CN201310706021.0A patent/CN103700102A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090316219A1 (en) * | 2008-06-18 | 2009-12-24 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method and computer-readable storage medium |
| CN101867790A (en) * | 2010-04-23 | 2010-10-20 | 刘文萍 | Millimeter-wave image analysis method, fire monitoring method and system |
| CN101950352A (en) * | 2010-05-31 | 2011-01-19 | 北京智安邦科技有限公司 | Target detection method capable of removing illumination influence and device thereof |
| CN103235939A (en) * | 2013-05-08 | 2013-08-07 | 哈尔滨工业大学 | Datum point positioning method based on machine vision |
| CN103400134A (en) * | 2013-06-28 | 2013-11-20 | 哈尔滨工程大学 | Non-contact method for extracting region of interest of finger vein sample |
Non-Patent Citations (2)
| Title |
|---|
| 李朝友: "指纹图像增强", 《中国优秀硕士学位论文全文数据库 信息科技辑》, 15 January 2011 (2011-01-15) * |
| 苏统华: "图像的边缘追踪及其在VC环境中的实现", 《电脑编程技巧与维护》, no. 6, 30 June 2002 (2002-06-30) * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105654525A (en) * | 2014-11-10 | 2016-06-08 | 中国石油化工股份有限公司 | Two-dimensional model construction method for reflecting real rock core connectivity |
| CN105654525B (en) * | 2014-11-10 | 2018-06-29 | 中国石油化工股份有限公司 | Reflect the 2D model building methods of true core connectivity |
| CN107194304A (en) * | 2016-03-15 | 2017-09-22 | 中国石油化工股份有限公司 | The connected region recognition methods of 3-dimensional digital rock core and device |
| CN107204003A (en) * | 2016-03-15 | 2017-09-26 | 中国石油化工股份有限公司 | The connected region recognition methods of two-dimensional digital rock core and device |
| CN107204003B (en) * | 2016-03-15 | 2020-04-17 | 中国石油化工股份有限公司 | Method and device for identifying connected area of two-dimensional digital core |
| CN107194304B (en) * | 2016-03-15 | 2020-09-04 | 中国石油化工股份有限公司 | Method and device for identifying connected area of three-dimensional digital core |
| CN110490890A (en) * | 2019-08-21 | 2019-11-22 | 青岛科技大学 | A method of quality digital core image processing is carried out with neural network |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105469094B (en) | A kind of edge vectors line drawing method of road surface bianry image | |
| US9625611B2 (en) | Method for detecting zonal underground target in mountain land based on ridge heat radiation correction | |
| CN102043958B (en) | High-definition remote sensing image multi-class target detection and identification method | |
| CN103198480B (en) | Based on the method for detecting change of remote sensing image of region and Kmeans cluster | |
| CN103871062B (en) | A kind of lunar surface rock detection method described based on super-pixel | |
| CN108109163A (en) | A kind of moving target detecting method for video of taking photo by plane | |
| CN103455813A (en) | Method for facula center positioning of CCD image measurement system | |
| CN104268853A (en) | Infrared image and visible image registering method | |
| CN104197933A (en) | Method for enhancing and extracting high stars and other sliding fixed stars in view field of telescope | |
| CN104408711A (en) | Multi-scale region fusion-based salient region detection method | |
| CN104680538A (en) | SAR image CFAR target detection method on basis of super pixels | |
| CN103700102A (en) | Rock core target extracting method based on CT (Computed Tomography) images | |
| CN103914829B (en) | Method for detecting edge of noisy image | |
| CN104376564A (en) | Method for extracting rough image edge based on anisotropism Gaussian directional derivative filter | |
| CN105631872A (en) | Remote sensing image registration method based on multiple feature points | |
| CN103822932B (en) | Based on the X ray realtime graphic weld defect method for detecting of multi-scale filtering operator | |
| CN104200458A (en) | MeanShift based high-resolution remote sensing image segmentation distance measurement optimization method | |
| Sundaresan et al. | Robustness of change detection algorithms in the presence of registration errors | |
| CN105913425A (en) | Self-adaptive oval blocking and wavelet transformation-based multi-pig contour extraction method | |
| CN104637060A (en) | Image partition method based on neighbor-hood PCA (Principal Component Analysis)-Laplace | |
| CN108764343B (en) | Method for positioning tracking target frame in tracking algorithm | |
| CN104063682A (en) | Pedestrian detection method based on edge grading and CENTRIST characteristic | |
| CN103065296B (en) | High-resolution remote sensing image residential area extraction method based on edge feature | |
| CN103955936B (en) | Significant object detection method based on stack-typed denoising self-coding machine | |
| CN103903258B (en) | Method for detecting change of remote sensing image based on order statistic spectral clustering |
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 |
Application publication date: 20140402 |
|
| WD01 | Invention patent application deemed withdrawn after publication |