CN105894444B - A kind of method and device based on CBCT video generation panoramic dental image - Google Patents
A kind of method and device based on CBCT video generation panoramic dental image Download PDFInfo
- Publication number
- CN105894444B CN105894444B CN201610196374.4A CN201610196374A CN105894444B CN 105894444 B CN105894444 B CN 105894444B CN 201610196374 A CN201610196374 A CN 201610196374A CN 105894444 B CN105894444 B CN 105894444B
- Authority
- CN
- China
- Prior art keywords
- image
- pyramid
- layer
- panoramic dental
- pixel
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000007408 cone-beam computed tomography Methods 0.000 title claims abstract description 54
- 238000001914 filtration Methods 0.000 claims abstract description 55
- 238000009792 diffusion process Methods 0.000 claims abstract description 52
- 238000005070 sampling Methods 0.000 claims abstract description 19
- 238000010586 diagram Methods 0.000 claims description 76
- 230000002708 enhancing effect Effects 0.000 claims description 62
- 238000000354 decomposition reaction Methods 0.000 claims description 38
- 238000002156 mixing Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 14
- 230000004927 fusion Effects 0.000 claims description 13
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 8
- 239000010931 gold Substances 0.000 description 8
- 229910052737 gold Inorganic materials 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000037123 dental health Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005549 size reduction Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/06—Topological mapping of higher dimensional structures onto lower dimensional surfaces
- G06T3/067—Reshaping or unfolding 3D tree structures onto 2D planes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30036—Dental; Teeth
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a kind of method and devices based on CBCT video generation panoramic dental image.The described method includes: determining seed point in the current axial face picture of cone beam computed tomography equipment, the seed point is handled using curve fitting algorithm, to generate arch wire；Using the curved surface where the arch wire as center curved surface, the corresponding tissue density CT value of the pixel on the central surface is calculated according to the sampling thickness of setting, initial panoramic dental image is generated according to the tissue density CT value；The Graded factor for determining anisotropic diffusion filtering formula, using the initial panoramic dental image of the anisotropic diffusion filtering formula manipulation different scale, to generate final panoramic dental image.The present invention realizes based on CBCT three-dimensional data, and the purpose of the two-dimensional panoramic image of denture is obtained by image processing algorithm.
Description
Technical field
The present embodiments relate to technical field of image processing, more particularly to one kind to be based on CBCT video generation dental panoramic
The method and device of image.
Background technique
With the improvement of living standards, people increasingly pay attention to dental health.Dental imaging equipment is in odontopathy
Diagnosing and treating during play the role of it is vital.
Currently, passing through Digital tomography equipment (also known as cone beam computed tomography equipment, the English of rebuilding dentistry
Literary abbreviation CBCT) acquisition three-dimensional data, obtain the three-dimensional structure image of Buccodental, be conducive to the surgery planning of odontopathy
And treatment.But be only capable of obtaining the three-dimensional structure image of Buccodental by CBCT, the two-dimensional panoramic of denture cannot be obtained
Image is unfavorable for the diagnosis that doctor carries out odontopathy by CBCT image.
Summary of the invention
The present invention provides a kind of method and device based on CBCT video generation panoramic dental image, is passed through with providing one kind
The three-dimensional data of cone beam computed tomography equipment (CBCT) acquisition obtains the mode of the two-dimensional panoramic image of denture.
In a first aspect, the embodiment of the invention provides a kind of method based on CBCT video generation panoramic dental image, packet
It includes:
Seed point is determined in the current axial face picture of cone beam computed tomography equipment, utilizes curve fitting algorithm
The seed point is handled, to generate arch wire；
Using the curved surface where the arch wire as center curved surface, calculated on the central surface according to the sampling thickness of setting
The corresponding tissue density CT value of pixel, initial panoramic dental image is generated according to the tissue density CT value；
The Graded factor for determining anisotropic diffusion filtering formula, not using the anisotropic diffusion filtering formula manipulation
With the initial panoramic dental image of scale, to generate final panoramic dental image.
Second aspect, the embodiment of the invention also provides a kind of devices based on CBCT video generation panoramic dental image.
Described device includes:
Arch wire generation unit, for determining seed in the current axial face picture of cone beam computed tomography equipment
Point handles the seed point using curve fitting algorithm, to generate arch wire；
First image generation unit, for using the curved surface where the arch wire as center curved surface, according to the sampling of setting
The corresponding tissue density CT value of pixel on central surface described in THICKNESS CALCULATION, according to the tissue density CT value
Generate initial panoramic dental image；
Second image generation unit, for determining the Graded factor of anisotropic diffusion filtering formula, using it is described respectively to
Anisotropic diffusion Filtering Formula handles the initial panoramic dental image of different scale, to generate final panoramic dental image.
The embodiment of the present invention by determining seed point in the current axial face picture of cone beam computed tomography equipment,
The seed point is handled using curve fitting algorithm, to generate arch wire.It is bent centered on the curved surface where the arch wire again
Face calculates the corresponding tissue density CT value of the pixel on the central surface according to the sampling thickness of setting, according to institute
It states tissue density CT value and generates initial panoramic dental image.Then, it is determined that the gradient of anisotropic diffusion filtering formula because
Son, using the initial panoramic dental image of the anisotropic diffusion filtering formula manipulation different scale, to generate final dentistry
Panoramic picture.The embodiment of the present invention solves asking for the current panoramic picture that denture cannot be obtained by CBCT three-dimensional data
Topic is realized based on CBCT three-dimensional data, the purpose of the two-dimensional panoramic image of denture is obtained by image processing algorithm,
The three-dimensional structure image that Buccodental can be not only obtained by CBCT, but also the two-dimensional panoramic figure of available denture are reached
The effect of picture.
Detailed description of the invention
Fig. 1 a is the process of method of one of the embodiment of the present invention one based on CBCT video generation panoramic dental image
Figure；
Fig. 1 b is axial in method of one of the embodiment of the present invention one based on CBCT video generation panoramic dental image
The schematic diagram of arch wire in face picture；
Fig. 2 a is the stream of another method based on CBCT video generation panoramic dental image in the embodiment of the present invention two
Cheng Tu；
Fig. 2 b is more rulers in method of one of the embodiment of the present invention two based on CBCT video generation panoramic dental image
Spend the flow chart for the treatment of process；
Fig. 3 is that the structure of device of one of the embodiment of the present invention three based on CBCT video generation panoramic dental image is shown
It is intended to.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining the present invention rather than limiting the invention.It also should be noted that in order to just
Only the parts related to the present invention are shown in description, attached drawing rather than entire infrastructure.
Embodiment one
Fig. 1 a is a kind of stream for method based on CBCT video generation panoramic dental image that the embodiment of the present invention one provides
Cheng Tu, the present embodiment are applicable to based on the three-dimensional data of cone beam computed tomography equipment, generate denture
The case where two-dimensional panoramic image, this method can be by being executed based on the device of CBCT video generation panoramic dental image.It is described
Method based on CBCT video generation panoramic dental image specifically comprises the following steps:
Step 110 determines seed point in the current axial face picture of cone beam computed tomography equipment, utilizes curve
The fitting algorithm processing seed point, to generate arch wire.
Wherein, cone beam computed tomography equipment, i.e. Cone beam computed tomograhpy, english abbreviation
For CBCT, principle, which is x-ray generator, does digital the throwings photograph of annular according to body around throwing with lower quantity of X-rays X.Then, it will surround
Data obtained in " intersection " obtain 3-D image after recombinating in a computer after the multiple number throwing of transillumination body is shone.CBCT image
It is from three-dimensional angle, i.e., sagittal plain, Coronal and axle position come lesions showed tissue and normal organization.
Determined in the current axial face picture of cone beam computed tomography equipment seed point mode can by a variety of,
In order to enhance human-computer interaction, seed point can be determined in current axial face picture in such a way that user inputs.For example, can be with
It is repeatedly clicked on current axial face by user, using the corresponding pixel of each clicking operation as seed point.Also
It can be user and input a series of coordinates, using the corresponding pixel of the coordinate as seed point.
The current axial face picture of cone beam computed tomography equipment as shown in Figure 1 b, is processing with the seed point
Object, using curve fitting algorithm by the seed point it is smooth connect obtain a whole smooth curve, by the song
Line is as arch wire.Wherein, curve fitting algorithm includes conic section fitting algorithm, Bi_arc fitting algorithm and discrete data curve
Fitting algorithm (common list point, which describes method, bezier curve and B-spline curves).
Step 120, using the curved surface where the arch wire as center curved surface, according to setting sampling thickness calculate it is described in
It is complete to generate initial dentistry according to the tissue density CT value for the corresponding tissue density CT value of pixel on innermost being face
Scape image.
Wherein, the curved surface where the arch wire can be the curved surface perpendicular to current axial face picture.It will be currently axial
One be located in the pixel on central surface in face picture is used as current pixel point, determines and passes through institute on the central surface
State the normal of current pixel point.It is thick as sampling at a distance from the current pixel point that the pixel on the normal will be located at
Degree.
Pass through the 3-dimensional image that recombination obtains since CBCT image may be considered a series of axial face images.It can sequence
Determine that the width in the axial face image is current axial face picture, as shown in Figure 1 b, firstly, by vertical where arch wire
The curved surface centered on the curved surface of current axial face picture determines and is located on the central surface in the current axial face picture
Pixel.Then, one in the pixel is sequentially obtained as current pixel point, is determined and is passed through on the central surface
The normal of the current pixel point.Secondly, calculate meet on the normal setting condition the corresponding tissue of pixel it is close
Spend the average value of CT value.Calculation formula can be with are as follows:Wherein, Vmin≤Vp≤Vmax,
VpFor the threshold range [V of settingmin, Vmax] in the corresponding tissue density CT value of pixel, DpFor the threshold value of above-mentioned setting
Pixel in range at a distance from the current pixel point,For sampling thickness range.
Using the average value as the corresponding tissue density CT value of the current pixel point on the central surface,
In, described impose a condition includes sampling thickness range of the pixel at a distance from the current pixel point in setting (referring to Fig. 1 b
Middle shadow region is sampling thickness range) in, and the corresponding tissue density CT value of pixel is in the threshold range of setting.
Wherein, the threshold range of the setting can be the research range for the tissue density CT value that user selectes, and avoid because introducing
More pixel is differed with the tissue density CT value of current pixel point, and leads to the tissue for calculating current pixel point
There is large error in density CT value.For example, if by some CT for differing biggish pixel with current dental tissue density CT value
Value (may be gap) introduces mean value calculation formula, necessarily will affect calculated result, in turn, leads to mean value calculation result not
Accurately.Finally, according to the stacking order of the axial face image, by the corresponding tissue density of pixel on central surface
CT value is combined to obtain the initial panoramic dental image.
Step 103, the Graded factor for determining anisotropic diffusion filtering formula, it is public using the anisotropic diffusion filtering
Formula handles the initial panoramic dental image of different scale, to generate final panoramic dental image.
Wherein, the diffusion coefficient of the anisotropic diffusion filtering has the characteristics that directional distribution, realizes image
Non-linear filtering operation, so that image border can be kept well while eliminating noise.Utilize image pyramid method
The initial panoramic dental image is decomposed into multiple dimensioned decomposition pyramid diagram picture, it is one corresponding that each layer decomposes pyramid diagram picture
The decomposition image of the panoramic dental image of scale.For example, make the initial panoramic dental image by low-pass filter into
Row is smooth, and then, carrying out line sampling to smoothed out image, (general sampling proportion is all low one layer in the horizontal and vertical directions
Decompose the 1/2 of pyramid diagram picture), to obtain a series of image of diminutions.The image that each layer reduces is known as the layer (scale
Under) decomposition pyramid diagram picture.Then, use interpolation method (for the inverse process of decomposition, and using the low pass with above-mentioned decomposable process
The identical filtering core of filter) inversion process that pyramid diagram picture executes decomposition is decomposed to high one layer, it obtains under different scale
Extension pyramid diagram picture, wherein the extension pyramid diagram picture is identical as the low one layer decomposition dimension of picture of pyramid diagram picture.
Laplacian-pyramid image, and each layer of drawing are determined according to the difference of extension pyramid diagram picture and low one layer of decomposition pyramid diagram picture
This pyramid diagram picture of pula is that this layer decomposes pyramid diagram picture with its high one layer of decomposition pyramid diagram as obtaining after interpolated amplification
Extension pyramid diagram picture difference.The contrast of image is improved by enhancing the Laplacian-pyramid image of different layers,
Noise can be also enhanced while the contrast enhancing of image, in order to inhibit noise, utilize enhanced Laplce's gold word
When tower image synthesizes final panoramic dental image, by the pyramidal gradient of enhanced Laplce at as Graded factor,
Anisotropic diffusion filtering is carried out to the blending image of each layer.Until detecting that initial layers enhance Laplacian-pyramid image
(initial layers are minimum one layer of level), using the gradient of initial layers enhancing laplacian pyramid as anisotropy parameter
The Graded factor of Filtering Formula merges the initial layers enhancing Laplacian-pyramid image and obtains with high one layer of blending image
The blending image of initial layers.Final dentistry is obtained using the blending image of initial layers described in anisotropic diffusion filtering formula manipulation
Panoramic picture.Wherein, anisotropic diffusion filtering formula is as follows:
Wherein, i is the number of plies for decomposing pyramid diagram picture, TiThe blending image that (m+x, n+y) is i-th layer, G (m, n) are volume
Product core,For the gradient of i-th layer of Laplacian-pyramid image.
The technical solution of the present embodiment, by being determined in the current axial face picture of cone beam computed tomography equipment
Seed point handles the seed point using curve fitting algorithm, to generate arch wire.It is with the curved surface where the arch wire again
Central surface calculates the corresponding tissue density CT value of the pixel on the central surface according to the sampling thickness of setting,
Initial panoramic dental image is generated according to the tissue density CT value.Then, it is determined that anisotropic diffusion filtering formula
Graded factor, using the initial panoramic dental image of the anisotropic diffusion filtering formula manipulation different scale, to generate most
Whole panoramic dental image.The technical solution of the present embodiment solves and cannot obtain denture by CBCT three-dimensional data at present
The problem of panoramic picture, realizes that based on CBCT three-dimensional data, the two dimension for obtaining denture by image processing algorithm is complete
The purpose of scape image has reached the three-dimensional structure image that Buccodental can be not only obtained by CBCT, but also available dental pattern
The effect of the two-dimensional panoramic image of tooth.
Embodiment two
Fig. 2 a is the stream of another method based on CBCT video generation panoramic dental image in the embodiment of the present invention two
Cheng Tu.The technical solution of the present embodiment on the basis of the technical solution of above-described embodiment, to determine anisotropic diffusion filtering
The Graded factor of formula, using the initial panoramic dental image of the anisotropic diffusion filtering formula manipulation different scale, with
The step of generating final panoramic dental image carries out further description, specifically comprises the following steps:
The initial panoramic dental image is decomposed into multiple dimensioned decomposition gold using image pyramid method by step 210
Word tower image.
Wherein, described image pyramid is a kind of effective but simple structure of concept for carrying out interpretation of images with multiresolution.
The pyramid of piece image is a series of image collection that resolution ratio with Pyramid arrangement gradually reduce.Image pyramid
Change generally comprises two steps: 1, utilizing low-pass filter smoothed image；2, smoothed image is sampled, to obtain a system
The image of column size reduction.For example, can choose 5 × 5 Gaussian convolution core (in addition, according to actual needs it is also an option that 3 ×
3 Gaussian convolution core or 7 × 7 Gaussian convolution core etc.) current initial panoramic dental image is smoothed.Then,
Even number line sampling is carried out to the image obtained after smoothing processing, to obtain the decomposition gold word of desired different layers (different scale)
Tower image (G0、G1、G2、…、Gk-2、Gk-1、Gk).Formula in decomposable process is as follows:
Wherein, W (m, n) is 5 × 5 Gaussian convolution core, GLFor L layers of decomposition pyramid diagram picture, L value is 0,1,2,
3 ..., k, k are the top layer level number for decomposing pyramid diagram picture.
Step 220 decomposes the inversion process that pyramid diagram picture executes decomposition to high one layer using interpolation method, is expanded
Pyramid diagram picture.
Even number row and column insertion 0, then, uses what is used in above-mentioned decomposable process in high one layer of decomposition pyramid diagram picture
5 × 5 Gaussian convolution core is filtered, and obtains extension gold identical with the current layer decomposition feature size of pyramid diagram picture
Word tower image.Pyramid diagram picture is decomposed to each layer respectively according to above-mentioned interpolation method and executes the inversion process decomposed, obtains k-1 layers
Extend pyramid diagram picture (Gt0、Gt1、Gt2、…Gt(k-1)).Formula in expansion process is as follows:
Wherein, W (m, n) is 5 × 5 Gaussian convolution core, GtbFor b layers of extension pyramid diagram picture, G(b+1)It is b+1 layers
Pyramid diagram picture is decomposed, b value is 0,1,2,3 ..., and k-1, k are the top layer level number for decomposing pyramid diagram picture.
Step 230 determines that the current layer decomposes the pixel of pyramid diagram picture and corresponding extension pyramid diagram picture respectively
Difference determines the initial panoramic dental image in the high frequency section image array of different scale according to the pixel value difference.
As shown in Figure 2 b, G is calculated0Layer decomposes pyramid diagram picture and Gt0Layer extends the pixel difference of pyramid diagram picture, will include
The pattern matrix of pixel value difference is as L0Layer Laplacian-pyramid image.Similarly, G is calculated1Layer decompose pyramid diagram picture and
Gt1The pixel difference of layer extension pyramid diagram picture will include the pattern matrix of pixel value difference as L1Layer laplacian pyramid figure
Picture.L can be calculated using same method2Layer Laplacian-pyramid image, L3Layer Laplacian-pyramid image ...,
Lk-2Layer Laplacian-pyramid image and Lk-1Layer Laplacian-pyramid image.Laplacian-pyramid image be it is described most
High frequency section image array of the first panoramic dental image in different scale.
Step 240 carries out enhancing contrast processing to the high frequency section image array, obtains high frequency section image enhancement
Matrix, using the gradient of the high frequency section image enhancement matrix as the Graded factor of anisotropic diffusion filtering formula.
Firstly, executing enhancing processing respectively to the Laplacian-pyramid image of different scale, different scale is obtained
Enhancing Laplacian-pyramid image.By the modes such as linear transformation or curve transformation, make the La Pu of different scale respectively
Curve of the Lars pyramid multiplied by setting, the mesh that the high frequency section image array (image detail) is enhanced with realization
, thus, obtain the enhancing Laplacian-pyramid image under different scale.The formula of above-mentioned enhancing process are as follows: Lep=
Enhance(Lp), wherein LepEnhance Laplacian-pyramid image for pth layer, Enhance is enhancing operator, LpFor pth layer
Laplacian-pyramid image, the value of p are 0,1,2,3 ..., and k-1, k are the top layer level number for decomposing pyramid diagram picture.Such as Fig. 2 b
It is shown, L0The enhanced processing of layer Laplacian-pyramid image obtains Le0Layer enhancing Laplacian-pyramid image.Similarly, L1
The enhanced processing of layer Laplacian-pyramid image obtains Le1Layer enhancing Laplacian-pyramid image.It is obtained using same method
To Le2Layer enhancing Laplacian-pyramid image, Le3Layer enhancing Laplacian-pyramid image ..., Le(k-2)Layer enhancing La Pula
This pyramid diagram picture and Le(k-1)Layer enhancing Laplacian-pyramid image.
Secondly, the gradient of the enhancing Laplacian-pyramid image of each scale is determined, using the gradient as anisotropy
Graded factor in diffusing filter formula.
Enhance the enhancing Laplce's gold for each scale of calculated for pixel values that Laplacian-pyramid image includes according to each layer
The gradient of word tower image.For example, the point centered on any pixel point that current layer enhancing Laplacian-pyramid image includes, meter
The region average gradient in setting regions corresponding with the pixel is calculated, it is thus possible to obtain current layer enhancing Laplce's gold
The gradient of word tower image.The gradient of current layer enhancing Laplacian-pyramid image is each as carrying out to current layer blending image
The Graded factor of Anisotropic diffusion filtering processing.Wherein, the current layer blending image can indicate are as follows: Tk=Up (Dk+1)+
Lek, wherein TkFor the blending image of kth layer, Up represents the operator of the inverse process decomposed, Dk+1It is complete for+1 layer of kth of decomposition dentistry
Scape image, LekEnhance Laplacian-pyramid image for kth layer.
Step 250, using the initial panoramic dental image of the anisotropic diffusion filtering formula manipulation different scale, with
Generate final panoramic dental image.
Maximum layer extension pyramid diagram picture and maximum layer enhancing Laplacian-pyramid image are subjected to fusion treatment,
Obtain maximum layer blending image.As shown in Figure 2 b, to GkLayer decomposes pyramid diagram picture and is extended processing, will treated expansion
Open up pyramid diagram picture and Le(k-1)Layer enhancing laplacian pyramid is merged, and the blending image T for -1 layer of kth is obtainedk-1(m+
X, n+y).
Using the gradient of the corresponding maximum layer enhancing Laplacian-pyramid image as Graded factor, by described
Maximum layer blending image described in anisotropic diffusion filtering formula manipulation obtains the decomposition panoramic dental image of maximum layer.
Above-mentioned calculating process can indicate are as follows: Dk-1=Filter (Gt(k-1)+Le(k-1)), wherein Dk-1It is complete for -1 layer of kth of decomposition dentistry
Scape image, Gt(k-1)For -1 layer of kth of extension panoramic dental image, Le(k-1)For -1 layer of enhancing Laplacian-pyramid image of kth,
K is the top layer level number for decomposing pyramid diagram picture, and Filter represents anisotropic diffusion filtering process.
Successively using high one layer of decomposition panoramic dental image as present image, by low one layer of enhancing laplacian pyramid
Image and the present image carry out fusion treatment, with the gradient of corresponding low one layer of enhancing Laplacian-pyramid image
Low one layer of decomposition tooth is obtained by the fused image of anisotropic diffusion filtering formula manipulation as Graded factor
Section's panoramic picture, until the Graded factor is minimum one layer of enhancing Laplacian-pyramid image, corresponding low one layer of point
Solution panoramic dental image is final panoramic dental image.
As shown in Figure 2 b, with Dk-1Layer decomposes panoramic dental image as present image, by Le(k-2)Layer enhancing Laplce
Pyramid diagram picture and the present image carry out fusion treatment, obtain Tk-2Layer blending image.With Le(k-2)Layer enhancing Laplce
The gradient of pyramid diagram picture passes through the anisotropic diffusion filtering formula manipulation T as Graded factork-2Layer blending image, obtains
To Dk-2Layer decomposes panoramic dental image.And so on, with D2Layer decomposes panoramic dental image as present image, by Le1Layer increases
Strong Laplacian-pyramid image and the present image carry out fusion treatment, obtain T1Layer blending image.With Le1Layer enhancing is drawn
The gradient of this pyramid diagram picture of pula passes through the anisotropic diffusion filtering formula manipulation T as Graded factor1Layer fusion figure
Picture obtains D1Layer decomposes panoramic dental image.Finally, with D1Layer decomposes panoramic dental image as present image, by Le0Layer increases
Strong Laplacian-pyramid image and the present image carry out fusion treatment, obtain T0Layer blending image.With Le0Layer enhancing is drawn
The gradient of this pyramid diagram picture of pula passes through the anisotropic diffusion filtering formula manipulation T as Graded factor0Layer fusion figure
Picture obtains D0Layer decomposes panoramic dental image.Due to Le0Layer enhancing Laplacian-pyramid image is that minimum one layer of enhancing drawing is general
Lars pyramid diagram picture, corresponding D0It is final panoramic dental image that layer, which decomposes panoramic dental image,.
Above-mentioned calculating process can indicate are as follows: DN=Filter (Up (DN+1)+LeN), wherein DNFor the decomposition tooth of n-th layer
Section's panoramic picture, Up represent the operator of the inverse process decomposed, LeNEnhance Laplacian-pyramid image for n-th layer, the value of N is
0,1,2,3 ..., k-2, wherein k is the top layer level number for decomposing pyramid diagram picture, and Filter represents anisotropic diffusion filtering mistake
Journey.
The technical solution of the present embodiment, by using the gradient conduct of the enhancing Laplacian-pyramid image of different scale
The Graded factor of Filtering Formula carries out anisotropic diffusion filtering processing to fused image, realizes and is weakening noise
Meanwhile retaining the marginal information of image, enhance the picture quality of final panoramic dental image.
Embodiment three
Fig. 3 is that the structure of device of one of the embodiment of the present invention three based on CBCT video generation panoramic dental image is shown
It is intended to.Described device specifically includes:
Arch wire generation unit 310, for being determined in the current axial face picture of cone beam computed tomography equipment
Seed point handles the seed point using curve fitting algorithm, to generate arch wire；
First image generation unit 320, for using the curved surface where the arch wire as center curved surface, according to adopting for setting
The corresponding tissue density CT value of pixel on central surface described in sample THICKNESS CALCULATION, according to the tissue density CT
Value generates initial panoramic dental image；
Second image generation unit 330, for determining the Graded factor of anisotropic diffusion filtering formula, using described each
Anisotropic diffusion Filtering Formula handles the initial panoramic dental image of different scale, to generate final panoramic dental image.
The technical solution of the present embodiment passes through the working as in cone beam computed tomography equipment of arch wire generation unit 310
Front axle determines seed point into face picture, the seed point is handled using curve fitting algorithm, to generate arch wire.Again by first
Image generation unit 320 executes operations described below --- using the curved surface where the arch wire as center curved surface, according to adopting for setting
The corresponding tissue density CT value of pixel on central surface described in sample THICKNESS CALCULATION, according to the tissue density CT
Value generates initial panoramic dental image.Then, anisotropic diffusion filtering formula is determined by the second image generation unit 330
Graded factor, using the initial panoramic dental image of the anisotropic diffusion filtering formula manipulation different scale, to generate most
Whole panoramic dental image.The technical solution of the present embodiment solves and cannot obtain denture by CBCT three-dimensional data at present
The problem of panoramic picture, realizes that based on CBCT three-dimensional data, the two dimension for obtaining denture by image processing algorithm is complete
The purpose of scape image has reached the three-dimensional structure image that Buccodental can be not only obtained by CBCT, but also available dental pattern
The effect of the two-dimensional panoramic image of tooth.
Further, the first image generation unit 320 is specifically used for:
Sequence determines that the width in the axial face image is current axial face picture, by where arch wire perpendicular to working as
Front axle determines the picture being located on the central surface in the current axial face picture to curved surface centered on the curved surface of face picture
Vegetarian refreshments.
One in the pixel is sequentially obtained as current pixel point, determines and is worked as on the central surface by described
The normal of preceding pixel point.
The average value for meeting the corresponding tissue density CT value of pixel of setting condition on the normal is calculated, by institute
Average value is stated as the corresponding tissue density CT value of current pixel point on the central surface, wherein the setting item
Part include pixel at a distance from the current pixel point within the scope of the sampling thickness of setting, and the corresponding human body group of pixel
Density CT value is knitted in the threshold range of setting.
According to the stacking order of the axial face image, by the corresponding tissue density CT of pixel on central surface
Value is combined to obtain the initial panoramic dental image.
Further, the second image generation unit 330 includes:
Subelement is decomposed, it is multiple dimensioned for being decomposed into the initial panoramic dental image using image pyramid method
Decompose pyramid diagram picture.
Subelement is extended, executes the inversion process decomposed for decomposing pyramid diagram picture to high one layer using interpolation method,
Be expanded pyramid diagram picture, wherein the extension pyramid diagram picture decomposes the dimension of picture phase of pyramid diagram picture with low one layer
Together.
Image array determines subelement, for determining that the current layer decomposes pyramid diagram picture and corresponding extension gold respectively
The pixel value difference of word tower image determines the initial panoramic dental image in the radio-frequency head of different scale according to the pixel value difference
Partial image matrix.
Graded factor determines subelement, for carrying out enhancing contrast processing to the high frequency section image array, obtains
High frequency section image enhancement matrix, using the gradient of the high frequency section image enhancement matrix as anisotropic diffusion filtering formula
Graded factor.
Further, described image pyramid is gaussian pyramid, and the high frequency section image array is Laplce's gold
Word tower image；
And Graded factor determines that subelement is specifically used for: to the Laplacian-pyramid image point of different scale
Not Zhi Hang enhancing processing, obtain the enhancing Laplacian-pyramid image of different scale；Determine the enhancing Laplce of each scale
The gradient of pyramid diagram picture, using the gradient as the Graded factor in anisotropic diffusion filtering formula.
Further, second image generation unit 330 further include:
Maximum layer blending image determines subelement, for enhancing maximum layer extension pyramid diagram picture and maximum layer
Laplacian-pyramid image carries out fusion treatment, obtains maximum layer blending image.
The decomposition panoramic dental image of maximum layer determines subelement, general for being drawn with the corresponding maximum layer enhancing
The gradient of Lars pyramid diagram picture passes through maximum layer described in the anisotropic diffusion filtering formula manipulation as Graded factor
Blending image obtains the decomposition panoramic dental image of maximum layer.
Final panoramic dental image determines subelement, for successively using high one layer of decomposition panoramic dental image as currently
Low one layer of enhancing Laplacian-pyramid image and the present image are carried out fusion treatment, with corresponding described low by image
The gradient of one layer of enhancing Laplacian-pyramid image passes through the anisotropic diffusion filtering formula manipulation as Graded factor
Fused image obtains low one layer of decomposition panoramic dental image, until the Graded factor is that minimum one layer of enhancing drawing is general
Lars pyramid diagram picture, corresponding low one layer of decomposition panoramic dental image are final panoramic dental image.
The above-mentioned device based on CBCT video generation panoramic dental image can be performed provided by any embodiment of the invention
Based on the method for CBCT video generation panoramic dental image, have the corresponding functional module of execution method and beneficial effect.
Note that the above is only a better embodiment of the present invention and the applied technical principle.It will be appreciated by those skilled in the art that
The invention is not limited to the specific embodiments described herein, be able to carry out for a person skilled in the art it is various it is apparent variation,
It readjusts and substitutes without departing from protection scope of the present invention.Therefore, although being carried out by above embodiments to the present invention
It is described in further detail, but the present invention is not limited to the above embodiments only, without departing from the inventive concept, also
It may include more other equivalent embodiments, and the scope of the invention is determined by the scope of the appended claims.
Claims (8)
1. a kind of method based on CBCT video generation panoramic dental image characterized by comprising
Seed point is determined in the current axial face picture of cone beam computed tomography equipment, is handled using curve fitting algorithm
The seed point, to generate arch wire；
Using the curved surface where the arch wire as center curved surface, the picture on the central surface is calculated according to the sampling thickness of setting
The corresponding tissue density CT value of vegetarian refreshments generates initial panoramic dental image according to the tissue density CT value；
The Graded factor for determining anisotropic diffusion filtering formula, using the anisotropic diffusion filtering formula manipulation difference ruler
The initial panoramic dental image of degree, to generate final panoramic dental image；
Wherein, the corresponding tissue density CT value of the pixel on the central surface is calculated according to the sampling thickness of setting,
Initial panoramic dental image is generated according to the tissue density CT value, comprising:
Sequence determines that the width in the axial face image is current axial face picture, by where arch wire perpendicular to current axis
Curved surface centered on to the curved surface of face picture determines the pixel being located on the central surface in the current axial face picture
Point；
One in the pixel is sequentially obtained as current pixel point, is determined on the central surface through the current picture
The normal of vegetarian refreshments；
The average value for meeting the corresponding tissue density CT value of pixel of setting condition on the normal is calculated, it will be described flat
Mean value is as the corresponding tissue density CT value of current pixel point on the central surface, wherein the setting condition packet
Pixel is included at a distance from the current pixel point within the scope of the sampling thickness of setting, and the corresponding tissue of pixel is close
CT value is spent in the threshold range of setting；
According to the stacking order of the axial face image, by the corresponding tissue density CT value of pixel on central surface into
Row combination obtains the initial panoramic dental image.
2. the method according to claim 1, wherein determine anisotropic diffusion filtering formula Graded factor,
Include:
The initial panoramic dental image is decomposed into multiple dimensioned decomposition pyramid diagram picture using image pyramid method；
Pyramid diagram picture is decomposed to high one layer using interpolation method and executes the inversion process decomposed, be expanded pyramid diagram picture,
Wherein, the extension pyramid diagram picture is identical as the low one layer decomposition dimension of picture of pyramid diagram picture；
Determine that current layer decomposes the pixel value difference of pyramid diagram picture and corresponding extension pyramid diagram picture respectively, according to the pixel
Difference determines the initial panoramic dental image in the high frequency section image array of different scale；
Enhancing contrast processing is carried out to the high frequency section image array, obtains high frequency section image enhancement matrix, it will be described
Graded factor of the gradient of high frequency section image enhancement matrix as anisotropic diffusion filtering formula.
3. according to the method described in claim 2, it is characterized in that, described image pyramid is gaussian pyramid, the high frequency
Parts of images matrix is Laplacian-pyramid image；
And using the gradient of the high frequency section image enhancement matrix as the Graded factor of anisotropic diffusion filtering formula,
Include:
Enhancing processing is executed respectively to the Laplacian-pyramid image of different scale, the enhancing drawing for obtaining different scale is general
Lars pyramid diagram picture；
The gradient for determining the enhancing Laplacian-pyramid image of each scale, using the gradient as anisotropic diffusion filtering public affairs
Graded factor in formula.
4. according to the method described in claim 3, it is characterized in that, different using the anisotropic diffusion filtering formula manipulation
The initial panoramic dental image of scale, to generate final panoramic dental image, comprising:
Maximum layer extension pyramid diagram picture and maximum layer enhancing Laplacian-pyramid image are subjected to fusion treatment, obtained
Maximum layer blending image；
Using the gradient of corresponding maximum layer enhancing Laplacian-pyramid image as Graded factor, by it is described respectively to
The anisotropic diffusion Filtering Formula processing maximum layer blending image, obtains the decomposition panoramic dental image of maximum layer；
Successively using high one layer of decomposition panoramic dental image as present image, by low one layer of enhancing Laplacian-pyramid image
With the present image carry out fusion treatment, using it is corresponding it is described it is low one layer enhancing Laplacian-pyramid image gradient as
It is complete to obtain low one layer of decomposition dentistry by the fused image of anisotropic diffusion filtering formula manipulation for Graded factor
Scape image, until the Graded factor is minimum one layer of enhancing Laplacian-pyramid image, corresponding low one layer of decomposition tooth
Section's panoramic picture is final panoramic dental image.
5. a kind of device based on CBCT video generation panoramic dental image characterized by comprising
Arch wire generation unit, for determining seed point in the current axial face picture of cone beam computed tomography equipment,
The seed point is handled using curve fitting algorithm, to generate arch wire；
First image generation unit, for using the curved surface where the arch wire as center curved surface, according to the sampling thickness of setting
The corresponding tissue density CT value of pixel on the central surface is calculated, is generated according to the tissue density CT value
Initial panoramic dental image；
Second image generation unit, for determining the Graded factor of anisotropic diffusion filtering formula, using the anisotropy
The initial panoramic dental image of diffusing filter formula manipulation different scale, to generate final panoramic dental image；
Wherein, the first image generation unit is specifically used for:
Sequence determines that the width in the axial face image is current axial face picture, by where arch wire perpendicular to current axis
Curved surface centered on to the curved surface of face picture determines the pixel being located on the central surface in the current axial face picture
Point；
One in the pixel is sequentially obtained as current pixel point, is determined on the central surface through the current picture
The normal of vegetarian refreshments；
The average value for meeting the corresponding tissue density CT value of pixel of setting condition on the normal is calculated, it will be described flat
Mean value is as the corresponding tissue density CT value of current pixel point on the central surface, wherein the setting condition packet
Pixel is included at a distance from the current pixel point within the scope of the sampling thickness of setting, and the corresponding tissue of pixel is close
CT value is spent in the threshold range of setting；
According to the stacking order of the axial face image, by the corresponding tissue density CT value of pixel on central surface into
Row combination obtains the initial panoramic dental image.
6. device according to claim 5, which is characterized in that the second image generation unit includes:
Subelement is decomposed, for the initial panoramic dental image to be decomposed into multiple dimensioned decomposition using image pyramid method
Pyramid diagram picture；
Subelement is extended, the inversion process decomposed is executed for decomposing pyramid diagram picture to high one layer using interpolation method, obtains
Extend pyramid diagram picture, wherein the extension pyramid diagram picture is identical as the low one layer decomposition dimension of picture of pyramid diagram picture；
Image array determines subelement, for determining that current layer decomposes pyramid diagram picture and corresponding extension pyramid diagram picture respectively
Pixel value difference, determine the initial panoramic dental image in the high frequency section image moment of different scale according to the pixel value difference
Battle array；
Graded factor determines subelement, for carrying out enhancing contrast processing to the high frequency section image array, obtains high frequency
Parts of images enhances matrix, using the gradient of the high frequency section image enhancement matrix as the ladder of anisotropic diffusion filtering formula
Spend the factor.
7. device according to claim 6, which is characterized in that described image pyramid is gaussian pyramid, the high frequency
Parts of images matrix is Laplacian-pyramid image；
And Graded factor determines that subelement is specifically used for:
Enhancing processing is executed respectively to the Laplacian-pyramid image of different scale, the enhancing drawing for obtaining different scale is general
Lars pyramid diagram picture；
The gradient for determining the enhancing Laplacian-pyramid image of each scale, using the gradient as anisotropic diffusion filtering public affairs
Graded factor in formula.
8. device according to claim 7, which is characterized in that second image generation unit further include:
Maximum layer blending image determines subelement, general for drawing maximum layer extension pyramid diagram picture with maximum layer enhancing
Lars pyramid diagram picture carries out fusion treatment, obtains maximum layer blending image；
The decomposition panoramic dental image of maximum layer determines subelement, for enhancing Laplce with the corresponding maximum layer
The gradient of pyramid diagram picture is merged as Graded factor by maximum layer described in the anisotropic diffusion filtering formula manipulation
Image obtains the decomposition panoramic dental image of maximum layer；
Final panoramic dental image determines subelement, for successively using high one layer of decomposition panoramic dental image as current figure
Low one layer of enhancing Laplacian-pyramid image and the present image are carried out fusion treatment, with corresponding described low one by picture
The gradient of layer enhancing Laplacian-pyramid image is melted as Graded factor by the anisotropic diffusion filtering formula manipulation
Image after conjunction obtains low one layer of decomposition panoramic dental image, until the Graded factor is minimum one layer of enhancing La Pula
This pyramid diagram picture, corresponding low one layer of decomposition panoramic dental image are final panoramic dental image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610196374.4A CN105894444B (en) | 2016-03-31 | 2016-03-31 | A kind of method and device based on CBCT video generation panoramic dental image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610196374.4A CN105894444B (en) | 2016-03-31 | 2016-03-31 | A kind of method and device based on CBCT video generation panoramic dental image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105894444A CN105894444A (en) | 2016-08-24 |
| CN105894444B true CN105894444B (en) | 2018-12-21 |
Family
ID=57011722
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610196374.4A Active CN105894444B (en) | 2016-03-31 | 2016-03-31 | A kind of method and device based on CBCT video generation panoramic dental image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105894444B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107301622B (en) * | 2017-06-20 | 2018-04-06 | 广州华端科技有限公司 | The panoramic picture acquisition methods and system of oral cavity CBCT images |
| KR101970989B1 (en) * | 2017-08-03 | 2019-08-13 | 오스템임플란트 주식회사 | Method and system for applying filters for dental CT imaging |
| CN110335281A (en) * | 2018-03-28 | 2019-10-15 | 北京连心医疗科技有限公司 | A kind of tumor boundaries determine method, equipment and storage medium |
| CN109377534B (en) * | 2018-10-22 | 2022-04-26 | 南方医科大学 | Nonlinear oral cavity CT panoramic image synthesis method capable of automatically sampling thickness detection |
| CN110211200B (en) * | 2019-04-22 | 2023-04-18 | 深圳安科高技术股份有限公司 | Dental arch wire generating method and system based on neural network technology |
| CN112270643A (en) * | 2020-09-04 | 2021-01-26 | 深圳市菲森科技有限公司 | Three-dimensional imaging data splicing method and device, electronic equipment and storage medium |
| CN112102495B (en) * | 2020-09-15 | 2023-11-21 | 北京朗视仪器股份有限公司 | Dental arch surface generation method based on CBCT image |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102663387A (en) * | 2012-04-16 | 2012-09-12 | 南京大学 | Cortical bone width automatic calculating method on basis of dental panorama |
| CN104166974A (en) * | 2013-08-01 | 2014-11-26 | 上海联影医疗科技有限公司 | CT locating film image enhancing method and CT locating film image enhancing device |
| CN104408700A (en) * | 2014-11-21 | 2015-03-11 | 南京理工大学 | Morphology and PCA (principal component analysis) based contourlet fusion method for infrared and visible light images |
| CN104809713A (en) * | 2015-04-24 | 2015-07-29 | 上海理工大学 | CBCT panorama nonlinear sharpening enhancing method based on neighborhood information and Gaussian filter |
-
2016
- 2016-03-31 CN CN201610196374.4A patent/CN105894444B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102663387A (en) * | 2012-04-16 | 2012-09-12 | 南京大学 | Cortical bone width automatic calculating method on basis of dental panorama |
| CN104166974A (en) * | 2013-08-01 | 2014-11-26 | 上海联影医疗科技有限公司 | CT locating film image enhancing method and CT locating film image enhancing device |
| CN104408700A (en) * | 2014-11-21 | 2015-03-11 | 南京理工大学 | Morphology and PCA (principal component analysis) based contourlet fusion method for infrared and visible light images |
| CN104809713A (en) * | 2015-04-24 | 2015-07-29 | 上海理工大学 | CBCT panorama nonlinear sharpening enhancing method based on neighborhood information and Gaussian filter |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105894444A (en) | 2016-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105894444B (en) | A kind of method and device based on CBCT video generation panoramic dental image | |
| CN109409503B (en) | Neural network training method, image conversion method, device, equipment and medium | |
| Yang et al. | 4D‐CT motion estimation using deformable image registration and 5D respiratory motion modeling | |
| EP2442278B1 (en) | Medical image processing apparatus, method, and program | |
| JP5584006B2 (en) | Projection image generation apparatus, projection image generation program, and projection image generation method | |
| EP1522045B1 (en) | Motion artifact correction of tomographical images | |
| CN106725570A (en) | Imaging method and system | |
| CN106600609A (en) | Spine partition method and system in medical image | |
| WO2007023723A1 (en) | Image processing method, image processing program, and image processing device | |
| CN102024251B (en) | System and method for multi-image based virtual non-contrast image enhancement for dual source CT | |
| JP2007537811A (en) | How to detect nodules | |
| Rashed et al. | Statistical image reconstruction from limited projection data with intensity priors | |
| CN111598989B (en) | Image rendering parameter setting method and device, electronic equipment and storage medium | |
| Angoth et al. | A novel wavelet based image fusion for brain tumor detection | |
| CN112669401B (en) | CT image reconstruction method and system based on convolutional neural network | |
| CN111524200B (en) | Method, apparatus and medium for segmenting a metal object in a projection image | |
| CN108520542A (en) | A kind of matched method for reconstructing of PET/CT data phase | |
| CN106910179A (en) | Multimode medical image fusion method based on wavelet transformation | |
| JP2018011870A (en) | Image processing device, image processing system, image processing method, and program | |
| Uccheddu et al. | 3D printing of cardiac structures from medical images: an overview of methods and interactive tools | |
| CN112154342A (en) | Medical image conversion | |
| CN111612689B (en) | Medical image processing method, medical image processing device, computer equipment and readable storage medium | |
| EP2548172A1 (en) | Functional image data enhancement and/or enhancer | |
| Ibraheem | Reduction of artifacts in dental cone beam CT images to improve the three dimensional image reconstruction | |
| CN110473297A (en) | Image processing method, device, electronic equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |