CN1723855A - The tooth proximal surface dental caries are decreased the processing method of X line image with computer - Google Patents
The tooth proximal surface dental caries are decreased the processing method of X line image with computer Download PDFInfo
- Publication number
- CN1723855A CN1723855A CN 200510012014 CN200510012014A CN1723855A CN 1723855 A CN1723855 A CN 1723855A CN 200510012014 CN200510012014 CN 200510012014 CN 200510012014 A CN200510012014 A CN 200510012014A CN 1723855 A CN1723855 A CN 1723855A
- Authority
- CN
- China
- Prior art keywords
- tooth
- dental caries
- bad
- line
- 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.)
- Pending
Links
- 208000002925 dental caries Diseases 0.000 title claims abstract description 82
- 238000003672 processing method Methods 0.000 title claims abstract description 8
- 230000003247 decreasing effect Effects 0.000 title abstract description 6
- 230000007423 decrease Effects 0.000 claims abstract description 13
- 239000000284 extract Substances 0.000 claims abstract description 5
- 210000004268 dentin Anatomy 0.000 claims description 23
- 239000006185 dispersion Substances 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000007373 indentation Methods 0.000 claims description 3
- 230000007170 pathology Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 12
- 238000003745 diagnosis Methods 0.000 abstract description 8
- 210000003141 lower extremity Anatomy 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 210000003298 dental enamel Anatomy 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
The present invention relates to a kind of processing method of the tooth proximal surface dental caries being decreased the X line image with computer; this method comprises the X line image that extracts pending tooth; in the X line image input computer with the pending tooth that obtains; mark and contain the edge that dental caries decrease tooth in the input picture; find out and in original image, mark and comprise that tooth is wide; bad number of dental caries; bad position of dental caries; bad width of dental caries, the bad degree of depth of dental caries, point data such as bad density of dental caries; it is of the present invention that method is simple; and can accurately mark the tooth dental caries through the X line original image that this method is handled and decrease the position, dental caries decrease density; the scope and the degree of depth have increased substantially diagnosis rate.
Description
Technical field
The present invention relates to a kind of processing method of the tooth proximal surface dental caries being decreased the X line image with computer.
Background technology
The X line original image that the early stage dental caries of traditional tooth decreases at interior proximal surface tooth dental caries as shown in Figure 1, can only be made a definite diagnosis with doctors experience for doctor's finding diagnosis reference again.Can not mark the tooth dental caries and decrease the position in X line original image, the tooth dental caries decrease density, scope and the degree of depth.The doctor can cause the mistaken diagnosis of the early stage dental caries of a large amount of teeth according to the diagnosis of X line original image and fail to pinpoint a disease in diagnosis.
Summary of the invention
The purpose of this invention is to provide a kind of processing method of the tooth proximal surface dental caries being decreased the X line image with computer, method is simple for this, and can accurately mark the tooth dental caries through the X line original image that this method is handled and decrease the position, dental caries decrease density, scope and the degree of depth, to increase substantially diagnosis rate.
Processing method of the present invention, it comprises the steps:
(1) extracts the X line image of pending tooth, in the X line image input computer with the pending tooth that obtains;
(2) mark the edge that contains dental caries damage tooth in the input picture;
The computing formula of grey scale curve such as (3) analysis of image grey scale change is inwardly done at tooth edge, the left and right sides in the image that marks, and obtains the grey scale curve such as level and smooth of tooth in the image, and is level and smooth is:
(4) with the grey scale curve such as level and smooth of obtaining tooth in the image, compare with the tooth dental caries damage data that the micro-X ray examination of pathology tooth slice has obtained, dental caries decrease the position and show as along the localized indentation zone in the tooth edge certain thickness scope on the grey scale curve such as level and smooth of tooth in obtaining image, according to comparison result, determine the minimum sinking degree and the scope in recessed zone, find all recessed zones rower of going forward side by side to annotate, greater than zero criterion as recessed zone, the second dervative computing formula is with second dervative:
Δ＝f(y-5)-2×f(y)+f(y+5)，
Calculate each recessed provincial characteristics parameter, width, the degree of depth, curvature and dispersion, the relevant parameter at the bad position of dental caries that itself and the micro-X ray examination of tooth slice are obtained is compared, and dispersion＜0.2 is satisfied one of following condition simultaneously and is recessed zone, bad position of dental caries:
(a) width＞tooth is wide by 9～11%, and the degree of depth＞tooth is wide by 4～6%;
(b) wide 1～3%＜width of tooth＜tooth is wide by 9～11%, and the degree of depth＞tooth is wide by 7～9%, and second dervative maximum＞2 in the zone;
Determining of the bad degree of depth of dental caries, width and glaze dentin line:
The determining of the bad degree of depth of dental caries scans the recessed zone of confirming, bad position of dental caries, finds out in this recessed zone in all level and smooth equi intensity curves apart from the edge concave point farthest, and this distance of putting the edge is defined as the bad degree of depth of dental caries;
Bad width of dental caries determine to find out in this recessed zone concave point the highest in all level and smooth equi intensity curves and minimum concave point, the vertical distance of this point to point is confirmed as bad width of dental caries;
The definite position of determining glaze dentin line is to extract left side glaze dentin line earlier, calculate tooth body left side Grad distribution, take out the point of gradient greater than this scope intermediate value, these some great majority all belong to the point on the glaze dentin line, remove spaced point, the remaining point that belongs to glaze dentin line is carried out certain left side glaze dentin line that smoothly just obtained.
The gradient calculation formula is as follows: and (x, y)=I (x-3, y)-I (x+3, y)
Tooth body right-hand part is repeated the left side identical operations, obtain right side glaze dentin line;
(5) computer marks the tooth decay zone with point, and gets the report of result data according to above data in original image, and described data comprise that tooth is wide, bad number of dental caries, bad position of dental caries, bad width of dental caries, the bad degree of depth of dental caries, bad density of dental caries.
Advantage of the present invention is simple and easy to do, and can accurately mark tooth dental caries damage position through the X line original image that this method is handled, and dental caries decrease density, scope and the degree of depth, have increased substantially diagnosis rate.
Description of drawings
Fig. 1 is an X line original image
Fig. 2 for the present invention along pending tooth edge the image when tooth body inside is made grey scale change and analyzed
Fig. 3 is grey scale curve figure (top is the curve chart of tooth body left side proximal surface, and the bottom is the curve chart of tooth body right side proximal surface) such as level and smooth of the present invention
Fig. 4 is through method processed images of the present invention
The specific embodiment
As shown in Figures 1 to 4, a kind of processing method of the tooth proximal surface dental caries being decreased the X line image with computer, it comprises the steps:
(1) extracts the X line image of pending tooth, in the X line image input computer with the pending tooth that obtains;
At first obtain the X line original image of bad tooth of dental caries, X line original image comprise by traditional X-ray sheet film (is good with the parallel projection bite wing film) scanning back obtain X-ray sheet beginning image and directly with digital X line original image indirectly; The former is memory image behind scanner, and the latter then directly deposits computer in pictorial form.
(2) mark the edge that contains dental caries damage tooth in the input picture; At first, from pending tooth X line (a tooth X line comprises many teeth), roughly scale off dental imaging to be detected (, seeing accompanying drawing 1) hereinafter to be referred as tooth to be measured.Determining to have comprised in this image whole tooth to be measured with naked eyes gets final product.
Tooth picture material to be measured is divided brighter prospect and darker background.Prospect promptly comprises a complete tooth to be measured that is positioned at picture centre, and the part of adjacent two teeth in its left and right sides, and root of the tooth is downward.The purpose of edge extracting is with tooth to be measured complete separating from image.
Edge extracting adopts rectangular histogram and fractal bonded method:
Separate prospect and background with histogram method, obtain prospect profile figure.Find the tooth (hereinafter to be referred as prospect) in the tooth sheet and the boundary gray value of tooth place oral environment background (hereinafter to be referred as background) with rectangular histogram.With this gray value of demarcating is that threshold value is separated prospect and background, obtains prospect profile figure.Contain the part that complete pending tooth reaches two teeth adjacent with it in the prospect profile of this moment.Therefore, prospect profile edge that extract this moment and fict tooth body edge have only provided the approximate location at tooth body edge, and tooth to be measured and its left and right sides be the marginal approximate location of adjacent teeth and the rough profile of tooth to be measured mutually.Demarcation line, tooth to be measured left side approximate location determine to get prospect profile figure left-half, tooth to be measured with its on the left of mutually the validity score boundary line of adjacent teeth be included in this part.Determine prospect top edge and lower limb at a distance of nearest point, and connect upward lower limb with straight line herein.This connecting line has provided marginal approximate location.The determining of demarcation line, tooth to be measured right side approximate location got prospect profile figure right half part, does the same treatment of determining with demarcation line, left side approximate location, obtains pending tooth and its right side marginal approximate location of adjacent teeth mutually.Determining of the rough profile of tooth to be measured is the common rough contour line that constitutes tooth to be measured of lower limb on demarcation line, the above-mentioned left and right sides and the prospect.The picture centre zone that they surround jointly is exactly the approximate range of tooth to be measured.
Determine to comprise the steps: the definite contour edge of tooth to be measured with fractal method
(a) put a little, on the rough contour line of tooth to be measured, put a little, and be that the heart is done the square block of 17 * 17 (or 33 * 33), select the piece size with these points with 8 (or 16) pel spacing, make these pieces cover true tooth border fully, the spacing of each piece central point is half of the piece length of side.
(b) block operations is appointed and is got one as original block.It is dwindled half (17 * 17 → 9 * 9, or 33 * 33 → 17 * 17) constitute the template piece.In original block, find out the fritter the most similar to template.Because comprise the tooth border in the original block, and the tooth border is obvious characteristics in the template, so still comprise real border in the fritter that so obtains in interior (can template to fritter poor quadratic sum or difference absolute value and minimum conduct " similar " criterion).As mentioned above, find out two non-intersect fritters in original block, these two fritters still can cover the border.
(c) traversal carries out operating among the b to all pieces that find among a, obtains 2 times of fritter and central points thereof to the original block number.As described in b, all fritters still cover the border.
(d) iteration, with among the c as a result fritter be new original block, repeat b, c step.Repeat all to make the piece yardstick to be reduced into original half, and all fritters cover borders at every turn.Repeat b, c step repeatedly, be reduced into 1 * 1 until the piece yardstick.This moment, all central points all were the real border points.With the expansion etch boundary point that finds is connected into continuous slick curve, obtain the tooth to be measured edge profile of cutting edge really.
(3) inwardly do the analysis of image grey scale change along tooth edge, the left and right sides in the image that marks, obtain the grey scale curve such as level and smooth of tooth in the image, along the tooth edge to grey scale curve such as the inner works of tooth body, gray scale with the layer of structure that shows X line image upper tooth body normal anomaly zone is distributed and situation of change, sees accompanying drawing 2.Grey scale curve such as level and smooth, computing formula is:
Wherein y is the upright position coordinate of putting on the equi intensity curve, and f (y) is a horizontal coordinate, and the i representative is counted.
(4) with the grey scale curve such as level and smooth of obtaining tooth in the image, the tooth dental caries that obtained with (human body of being made by the doctor) micro-X ray examination of pathology tooth slice decrease data (importing in the computer) and compare, in obtaining image on the grey scale curve such as level and smooth of tooth, dental caries decrease the position and show as along the localized indentation zone in the tooth edge certain thickness scope, according to comparison result, determine the minimum sinking degree and the scope in recessed zone, find all recessed zones rower of going forward side by side to annotate
Tooth to be measured left and right sides proximal surface is carried out recessed area marking.Therefore along tooth edge, the left and right sides etc. grey scale curve, in the certain thickness scope, search for, find all recessed zones rower of going forward side by side to annotate.Greater than zero criterion as " recessed ", second dervative computing formula used herein is with second dervative:
Δ＝f(y-5)-2×f(y)+f(y+5)，
Δ is the second dervative value in the formula, and wherein y is the upright position coordinate of putting on the equi intensity curve, and f (y) is a horizontal coordinate, and f (y) is a horizontal coordinate.
Remove and the bad irrelevant recessed zone of dental caries.
Calculate each recessed provincial characteristics parameter.Comprise that the width, the degree of depth perpendicular to edge direction, the curvature of second dervative, the dispersion that are parallel to edge direction are concave point in recessed region area (width x depth)/recessed zone (second dervative＞0) number; the relevant parameter at the bad position of dental caries that recessed provincial characteristics parameter and the micro-X ray examination of tooth slice are obtained is compared, and removes the recessed zone that does not conform to bad positional parameter of dental caries.Removing the recessed zone of being left behind the irrelevant person is exactly final dental caries bad (doubtful) district of assert.Dispersion＜0.2 is satisfied one of following condition person simultaneously and is regarded as dental caries bad (doubtful) position:
(a) width＞tooth is wide by 10%, and the degree of depth＞tooth is wide by 5%;
(b) the wide 2%＜width of tooth＜tooth is wide by 10%, and the degree of depth＞tooth is wide by 8%, and second dervative maximum＞2 in the zone.
Determining of the bad degree of depth of dental caries: bad doubtful district is scanned to the dental caries of confirming, finds out in the zone in all equi intensity curves apart from the edge concave point farthest, and this distance of putting the edge is defined as the bad degree of depth of dental caries.
Determining of bad width of dental caries: find out concave point the highest in interior all equi intensity curves in zone and minimum concave point, the vertical distance of this point to point is confirmed as bad width of dental caries
Dento enamel junction is that the higher part of tooth external body density is omited the marginal definite of lower part tooth body structure with internal density.As the tolerance of the bad order of severity of dental caries, the relative position of needs report dental caries evil idea and glaze dentin line.Therefore be necessary to determine the definite position of glaze dentin line.
Left side glaze dentin line drawing is because glaze dentin line is the demarcation line that higher part of tooth body marginal density and internal density omit lower part tooth body structure, so the gradient of image and gray scale of glaze dentin line position is bigger than other position gradient of tooth body.Calculate tooth body left side Grad distribution, take out the point of gradient greater than this scope intermediate value.These somes overwhelming majority belongs to glaze dentin line, has only indivedual spaced points.Remove spaced point, the remaining point that belongs to glaze dentin line is carried out certain left side glaze dentin line that smoothly just obtained.
The gradient calculation formula is as follows: and (x, y)=I (x-3, y)-I (x+3, y).Wherein (x, y) and f (x, y) be respectively image point (x, Grad y) and gray value, x wherein, y is a pixel coordinates, (x y) is x to I, the gray value of y point pixel.
Tooth body right-hand part is repeated the left side identical operations, obtain right side glaze dentin line.
(5) computer marks the tooth decay zone with point, as shown in Figure 4 according to above data in original image.And the report of getting result data, described data comprise that tooth is wide, bad number of dental caries, bad position of dental caries, bad width of dental caries, the bad degree of depth of dental caries, bad density of dental caries.
In conjunction with the data that obtain, in X line original image, mark the image in band tooth decay district, and, finish through method processed images of the present invention with its output print.
In a word, the inventive method is simple and easy to do, improves sensitivity, the accuracy of X radiodiagnosis x, reaches and exceed expert's naked eyes diagnostic level of film, and with guiding clinical treatment, other can be used for the Basic Experiment Study of early stage dental caries disease, can promote the use of.
Claims (1)
1, a kind of with the processing method of computer to tooth proximal surface dental caries damage X line image, it is characterized in that it comprises the steps:
(1) extracts the X line image of pending tooth, in the X line image input computer with the pending tooth that obtains;
(2) mark the edge that contains dental caries damage tooth in the input picture;
The computing formula of grey scale curve such as (3) analysis of image grey scale change is inwardly done at tooth edge, the left and right sides in the image that marks, and obtains the grey scale curve such as level and smooth of tooth in the image, and is level and smooth is:
(4) with the grey scale curve such as level and smooth of obtaining tooth in the image, compare with the tooth dental caries damage data that the micro-X ray examination of pathology tooth slice has obtained, dental caries decrease the position and show as along the localized indentation zone in the tooth edge certain thickness scope on the grey scale curve such as level and smooth of tooth in obtaining image, according to comparison result, determine the minimum sinking degree and the scope in recessed zone, find all recessed zones rower of going forward side by side to annotate, greater than zero criterion as recessed zone, the second dervative computing formula is with second dervative:
Δ＝f(y-5)-2×f(y)+f(y+5)，
Calculate each recessed provincial characteristics parameter, width, the degree of depth, curvature and dispersion, the relevant parameter at the bad position of dental caries that itself and the micro-X ray examination of tooth slice are obtained is compared, and dispersion＜0.2 is satisfied one of following condition simultaneously and is recessed zone, bad position of dental caries:
(a) width＞tooth is wide by 9～11%, and the degree of depth＞tooth is wide by 4～6%;
(b) wide 1～3%＜width of tooth＜tooth is wide by 9～11%, and the degree of depth＞tooth is wide by 7～9%, and second dervative maximum＞2 in the zone;
Determining of the bad degree of depth of dental caries, width and glaze dentin line:
The determining of the bad degree of depth of dental caries scans the recessed zone of confirming, bad position of dental caries, finds out in this recessed zone in all level and smooth equi intensity curves apart from the edge concave point farthest, and this distance of putting the edge is defined as the bad degree of depth of dental caries;
Bad width of dental caries determine to find out in this recessed zone concave point the highest in all level and smooth equi intensity curves and minimum concave point, the vertical distance of this point to point is confirmed as bad width of dental caries;
The definite position of determining glaze dentin line is to extract left side glaze dentin line earlier, calculate tooth body left side Grad distribution, take out the point of gradient greater than this scope intermediate value, these some great majority all belong to the point on the glaze dentin line, remove spaced point, the remaining point that belongs to glaze dentin line is carried out certain left side glaze dentin line that smoothly just obtained
The gradient calculation formula is as follows: and (x, y)=I (x-3, y)-I (x+3, y)
Tooth body right-hand part is repeated the left side identical operations, obtain right side glaze dentin line;
(5) computer marks the tooth decay zone with point, and gets the report of result data according to above data in original image, and described data comprise that tooth is wide, bad number of dental caries, bad position of dental caries, bad width of dental caries, the bad degree of depth of dental caries, bad density of dental caries.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200510012014 CN1723855A (en) | 2005-06-27 | 2005-06-27 | The tooth proximal surface dental caries are decreased the processing method of X line image with computer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200510012014 CN1723855A (en) | 2005-06-27 | 2005-06-27 | The tooth proximal surface dental caries are decreased the processing method of X line image with computer |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN1723855A true CN1723855A (en) | 2006-01-25 |
Family
ID=35923573
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200510012014 Pending CN1723855A (en) | 2005-06-27 | 2005-06-27 | The tooth proximal surface dental caries are decreased the processing method of X line image with computer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1723855A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102564923A (en) * | 2011-12-29 | 2012-07-11 | 上海交通大学 | Device for automatically identifying fiber distribution density of bone marrow biopsy |
| US8737706B2 (en) | 2009-06-16 | 2014-05-27 | The University Of Manchester | Image analysis method |
| CN104657519A (en) * | 2013-11-18 | 2015-05-27 | 同方威视技术股份有限公司 | Method for building enamel-dentinal junction statistical averaging model |
| CN111528846A (en) * | 2020-04-30 | 2020-08-14 | 赤峰学院附属医院 | Oral craniomaxillofacial scanning device and scanning method and electronic device |
-
2005
- 2005-06-27 CN CN 200510012014 patent/CN1723855A/en active Pending
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8737706B2 (en) | 2009-06-16 | 2014-05-27 | The University Of Manchester | Image analysis method |
| CN102564923A (en) * | 2011-12-29 | 2012-07-11 | 上海交通大学 | Device for automatically identifying fiber distribution density of bone marrow biopsy |
| CN102564923B (en) * | 2011-12-29 | 2013-07-17 | 上海交通大学 | Device for automatically identifying fiber distribution density of bone marrow biopsy |
| CN104657519A (en) * | 2013-11-18 | 2015-05-27 | 同方威视技术股份有限公司 | Method for building enamel-dentinal junction statistical averaging model |
| CN111528846A (en) * | 2020-04-30 | 2020-08-14 | 赤峰学院附属医院 | Oral craniomaxillofacial scanning device and scanning method and electronic device |
| CN111528846B (en) * | 2020-04-30 | 2021-05-14 | 赤峰学院附属医院 | Oral craniomaxillofacial scanning device and scanning method and electronic device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9439610B2 (en) | Method for teeth segmentation and alignment detection in CBCT volume | |
| Skinner et al. | Enamel-dentine junction (EDJ) morphology distinguishes the lower molars of Australopithecus africanus and Paranthropus robustus | |
| Kumar et al. | Descriptive analysis of dental X-ray images using various practical methods: A review | |
| CN110889850B (en) | CBCT tooth image segmentation method based on central point detection | |
| JPH11510419A (en) | Quantitative dental caries inspection system and method | |
| US7460700B2 (en) | Image processing method and computer-readable recording medium containing image processing program | |
| CN108648825B (en) | A kind of leucoderma hickie appraisal procedure based on image recognition | |
| Meintjes et al. | A stereo-photogrammetric method to measure the facial dysmorphology of children in the diagnosis of fetal alcohol syndrome | |
| CN113223010A (en) | Method and system for fully automatically segmenting multiple tissues of oral cavity image | |
| Fornai et al. | Variability of Australopithecus second maxillary molars from Sterkfontein Member 4 | |
| CN1723855A (en) | The tooth proximal surface dental caries are decreased the processing method of X line image with computer | |
| CN1505490A (en) | Method and apparatus for the non-invasive imaging of anatomic tissue structures | |
| CN111325754B (en) | Automatic lumbar vertebra positioning method based on CT sequence image | |
| CN106780492B (en) | Method for extracting key frame of CT pelvic image | |
| CN1764930A (en) | Imaging inside structure | |
| CN105678711B (en) | A kind of attenuation correction method based on image segmentation | |
| CN113786262A (en) | Preparation method of dental crown lengthening operation guide plate | |
| Lazzari et al. | Quantitative three-dimensional topography in taxonomy applied to the dental morphology of catarrhines | |
| Hopp et al. | Breast tissue characterization by sound speed: correlation with mammograms using a 2D/3D image registration | |
| Cao et al. | Liver fibrosis identification based on ultrasound images | |
| KR101334029B1 (en) | Method for extracting muscle of head and neck and method of ＲＯＩ analysis | |
| CN106340014A (en) | Binary curve processing method | |
| Schwartz | Taxonomic and functional aspects of enamel cap structure in South Africa Plio-Pleistocene hominids: A high-resolution computed tomographic study | |
| Dhar et al. | Automatic tracing of mandibular canal pathways using deep learning | |
| Nirmala et al. | Measurement of nuchal translucency thickness in first trimester ultrasound fetal images for detection of chromosomal abnormalities |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |