CN105894444A - Method and device for generating dental panoramic image on the basis of CBCT image - Google Patents
Method and device for generating dental panoramic image on the basis of CBCT image Download PDFInfo
- Publication number
- CN105894444A CN105894444A CN201610196374.4A CN201610196374A CN105894444A CN 105894444 A CN105894444 A CN 105894444A CN 201610196374 A CN201610196374 A CN 201610196374A CN 105894444 A CN105894444 A CN 105894444A
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 85
- 238000007408 cone-beam computed tomography Methods 0.000 title claims abstract description 54
- 238000009792 diffusion process Methods 0.000 claims abstract description 46
- 238000005070 sampling Methods 0.000 claims abstract description 21
- 238000010586 diagram Methods 0.000 claims description 71
- 238000001914 filtration Methods 0.000 claims description 49
- 230000008569 process Effects 0.000 claims description 47
- 238000000354 decomposition reaction Methods 0.000 claims description 32
- 230000002708 enhancing effect Effects 0.000 claims description 22
- 230000004927 fusion Effects 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000005728 strengthening Methods 0.000 claims description 10
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 4
- 239000010931 gold Substances 0.000 claims description 4
- 229910052737 gold Inorganic materials 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 210000002455 dental arch Anatomy 0.000 abstract 2
- 210000001519 tissue Anatomy 0.000 abstract 2
- 230000000694 effects Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 208000010641 Tooth disease Diseases 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000005056 cell body Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000037123 dental health Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 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
- 230000000717 retained effect Effects 0.000 description 1
- 238000005549 size reduction Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000002560 therapeutic 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 method and device for generating a dental panoramic image on the basis of a CBCT image. The method comprises steps of: determining a seed point in a current axial surface image of a CBCT apparatus, processing the seed point by using a curve fitting algorithm to generate a dental arch line; using a curve surface where the dental arch line is located as a central curve surface, computing body tissue density CT values corresponding to the pixels in the central curve surface according to a set sampling thickness, and generating initial dental panoramic images according to the body tissue density CT values; and determining the gradient factor of an anisotropic diffusion filter formula and processing the initial dental panoramic images with different scales by using the anisotropic diffusion filter formula to generate an final dental panoramic image. The method and device achieve a purpose of obtaining the two-dimensional panoramic image of the entire teeth by using the image processing algorithm on the basis of CBCT three-dimensional data.
Description
Technical field
The present embodiments relate to technical field of image processing, particularly relate to a kind of based on CBCT video generation
The method and device of panoramic dental image.
Background technology
Along with growth in the living standard, people increasingly pay attention to for dental health.Dental imaging equipment is at tooth
The diagnosis of tooth disease and therapeutic process serve vital effect.
At present, (it is also called cone beam computed tomography to set by the Digital tomography equipment of reconstruction dentistry
Standby, English abbreviation CBCT) three-dimensional data that gathers, obtain the three-dimensional structure image of Buccodental, be conducive to
The surgery planning of odontopathy and treatment.But, it is only capable of obtaining the three-dimensional structure of Buccodental by CBCT
Image, it is impossible to obtain the two-dimensional panoramic image of denture, be unfavorable for that doctor carries out tooth by CBCT image
The diagnosis of tooth disease.
Summary of the invention
The present invention provides a kind of method and device based on CBCT video generation panoramic dental image, to provide
A kind of three-dimensional data gathered by cone beam computed tomography equipment (CBCT) obtains the two of denture
The mode of dimension panoramic picture.
First aspect, embodiments provides a kind of based on CBCT video generation panoramic dental image
Method, including:
In the current axially face picture of cone beam computed tomography equipment, determine Seed Points, utilize curve to intend
Hop algorithm processes described Seed Points, to generate arch wire；
Curved surface centered by the curved surface at described arch wire place, calculates described center according to the sampling thickness set
The tissue density CT value that pixel on curved surface is corresponding, generates according to described tissue density CT value
Initially panoramic dental image；
Determine the Graded factor of anisotropic diffusion filtering formula, use described anisotropic diffusion filtering formula
Process the initial panoramic dental image of different scale, to generate final panoramic dental image.
Second aspect, the embodiment of the present invention additionally provides a kind of based on CBCT video generation panoramic dental image
Device.Described device includes:
Arch wire signal generating unit, for true in the current axially face picture of cone beam computed tomography equipment
Determine Seed Points, utilize curve fitting algorithm to process described Seed Points, to generate arch wire；
First image generation unit, for curved surface centered by the curved surface at described arch wire place, according to setting
Sampling thickness calculate the tissue density CT value that the pixel on described central surface is corresponding, according to institute
State tissue density CT value and generate initial panoramic dental image；
Second image generation unit, for determining the Graded factor of anisotropic diffusion filtering formula, uses institute
State the initial panoramic dental image of anisotropic diffusion filtering formula manipulation different scale, to generate final dentistry
Panoramic picture.
The embodiment of the present invention is by determining kind in the current axially face picture of cone beam computed tomography equipment
Sub-, utilize curve fitting algorithm to process described Seed Points, to generate arch wire.Again with described arch wire institute
Curved surface centered by curved surface, calculate the pixel on described central surface according to the sampling thickness set corresponding
Tissue density CT value, according to described tissue density CT value generate initial panoramic dental image.
Then, it is determined that the Graded factor of anisotropic diffusion filtering formula, use described anisotropic diffusion filtering public
Formula processes the initial panoramic dental image of different scale, to generate final panoramic dental image.The present invention implements
The problem illustrating the panoramic picture that can not be obtained denture the most at present by CBCT three-dimensional data, it is achieved with
Based on CBCT three-dimensional data, obtained the mesh of the two-dimensional panoramic image of denture by image processing algorithm
, reach both to be obtained by CBCT the three-dimensional structure image of Buccodental, full mouth can have been obtained again
The effect of the two-dimensional panoramic image of tooth.
Accompanying drawing explanation
Fig. 1 a is a kind of based on CBCT video generation panoramic dental image the method in the embodiment of the present invention one
Flow chart；
Fig. 1 b is a kind of based on CBCT video generation panoramic dental image the method in the embodiment of the present invention one
In the axial schematic diagram of arch wire in face picture；
Fig. 2 a is the side based on CBCT video generation panoramic dental image of the another kind in the embodiment of the present invention two
The flow chart of method；
Fig. 2 b is a kind of based on CBCT video generation panoramic dental image the method in the embodiment of the present invention two
In the flow chart of multiple dimensioned processing procedure；
Fig. 3 is a kind of based on CBCT video generation panoramic dental image the device in the embodiment of the present invention three
Structural representation.
Detailed description of the invention
The present invention is described in further detail with embodiment below in conjunction with the accompanying drawings.It is understood that this
Specific embodiment described by place is used only for explaining the present invention, rather than limitation of the invention.The most also need
It is noted that for the ease of describing, accompanying drawing illustrate only part related to the present invention and not all knot
Structure.
Embodiment one
A kind of based on CBCT video generation panoramic dental image the side that Fig. 1 a provides for the embodiment of the present invention one
The flow chart of method, the present embodiment is applicable to based on the three-dimensional data of cone beam computed tomography equipment,
Generating the situation of the two-dimensional panoramic image of denture, the method can be by based on CBCT video generation dentistry
The device of panoramic picture performs.Described method based on CBCT video generation panoramic dental image is specifically wrapped
Include following steps:
Step 110, in the current axially face picture of cone beam computed tomography equipment, determine Seed Points, profit
Described Seed Points is processed, to generate arch wire with curve fitting algorithm.
Wherein, cone beam computed tomography equipment, i.e. Cone beam computed tomograhpy, English
CBCT be abbreviated as in literary composition, and its principle is that x-ray generator does annular numeral around throwing according to body with relatively low quantity of X-rays X
Formula is thrown and is shone.Then, will be around transillumination body repeatedly numeral to throw according to the data obtained in rear " common factor " in calculating
Machine obtains 3-D view after restructuring.CBCT image is from three-dimensional angle, i.e. sagittal plain, Coronal and
Lesions showed tissue and normal organization are come in axle position.
In the current axially face picture of cone beam computed tomography equipment, determine that the mode of Seed Points can be by
Multiple, in order to strengthen man-machine interaction, can determine in current axially face picture by the way of user inputs
Seed Points.For example, it is possible to repeatedly clicked on current axially face by user, behaviour will be clicked on every time
Make corresponding pixel as Seed Points.Can also is that user inputs a series of coordinate, corresponding with described coordinate
Pixel as Seed Points.
The current axially face picture of cone beam computed tomography equipment as shown in Figure 1 b, with described Seed Points
For processing object, use curve fitting algorithm that smooth the coupling together of described Seed Points is obtained a whole piece and smooth
Curve, using described curve as arch wire.Wherein, curve fitting algorithm include conic section fitting algorithm,
(conventional list point describes method bezier curve for Bi_arc fitting algorithm and discrete data curve fitting algorithm
And B-spline curves).
Step 120, centered by the curved surface at described arch wire place curved surface, according to set sampling thickness meter
Calculate the tissue density CT value that the pixel on described central surface is corresponding, close according to described tissue
Degree CT value generates initial panoramic dental image.
Wherein, the curved surface at described arch wire place can be perpendicular to the curved surface of current axially face picture.Ought
Front axle is positioned in the pixel on central surface as current pixel point in face picture, determines described
By the normal of described current pixel point on central surface.The pixel that will be located on described normal is worked as with described
The distance of preceding pixel point is as sampling thickness.
Owing to CBCT image may be considered the 3-dimensional image that a series of axial face image obtains through restructuring.
Can sequentially determine that the width in described axial face image is current axially face picture, as shown in Figure 1 b, first,
Using the curved surface being perpendicular to current axially face picture at arch wire place as central surface, determine described current axis
The pixel on described central surface it is positioned in face picture.Subsequently, order obtains in described pixel
Individual as current pixel point, determine the normal by described current pixel point on described central surface.Secondly,
Calculate the mean value meeting tissue density CT value corresponding to the pixel imposed a condition on described normal.
Computing formula can be:Wherein, Vmin≤Vp≤Vmax,VpFor
Threshold range [the V setmin, VmaxThe tissue density CT value that pixel in] is corresponding, DpFor above-mentioned
Pixel in the threshold range set and the distance of described current pixel point,For sampling thickness model
Enclose.
Using described mean value as tissue density CT corresponding to the current pixel point on described central surface
Value, wherein, described in the distance including pixel and described current pixel point that imposes a condition thick in the sampling set
In degree scope (seeing shadow region in Fig. 1 b is sampling thickness scope), and the tissue that pixel is corresponding
Density CT value is in the threshold range set.Wherein, the threshold range of described setting can be that user selectes
The research range of tissue density CT value, it is to avoid because introducing the tissue density with current pixel point
The pixel that CT value difference is more, and cause the tissue density CT value calculating current pixel point to occur relatively
Big error.Such as, if some to be differed the CT value of bigger pixel with current dental tissue density CT value
(being probably space) introduces mean value calculation formula, will necessarily affect result of calculation, and then, cause average
Value result of calculation is inaccurate.Finally, according to the stacking order of described axial face image, by central surface
The tissue density CT value that pixel is corresponding is combined obtaining described initial panoramic dental image.
Step 103, determine the Graded factor of anisotropic diffusion filtering formula, use described anisotropy to expand
Dissipate Filtering Formula and process the initial panoramic dental image of different scale, to generate final panoramic dental image.
Wherein, the diffusion coefficient of described anisotropic diffusion filtering has the feature of directional distribution, it is achieved that
The non-linear filtering operation of image so that while eliminating noise, it is possible to keep image border well.
Utilize image pyramid method that described initial panoramic dental image is decomposed into multiple dimensioned decomposition pyramid diagram picture,
Each layer decomposes the exploded view picture of the panoramic dental image of the corresponding yardstick of pyramid diagram picture.Such as, described in making
Initially panoramic dental image smooths through a low pass filter, then, carries out the image after smooth
Line sampling (general sampling proportion is the most all low one layer and decomposes the 1/2 of pyramid diagram picture),
Thus obtain a series of image reduced.The image that each layer reduces is referred to as the decomposition of this layer (under this yardstick)
Pyramid diagram picture.Subsequently, interpolation method is used (to be the inverse process decomposed, and use and above-mentioned decomposable process
The filtering core that low pass filter is identical) decompose, to high one layer, the inversion process that pyramid diagram picture performs to decompose,
Obtain the extension pyramid diagram picture under different scale, wherein, described extension pyramid diagram picture and low one layer of decomposition
The dimension of picture of pyramid diagram picture is identical.According to extension pyramid diagram picture with low one layer decompose pyramid diagram picture it
Difference determines Laplacian-pyramid image, and each layer of Laplacian-pyramid image is this layer of decomposition pyramid
High with it one layer of image decomposes pyramid diagram as the difference extending pyramid diagram picture obtained after interpolated amplification.Logical
Cross strengthen different layers Laplacian-pyramid image to improve the contrast of image, image contrast increase
While Qiang, noise also can be enhanced, and in order to suppress noise, is utilizing enhanced laplacian pyramid figure
During as synthesizing final panoramic dental image, the pyramidal gradient of enhanced Laplce is become as gradient
The factor, carries out Anisotropic diffusion filtering to the fused images of each layer.Until, detect that initial layers strengthens and draw
This pyramid diagram picture of pula (initial layers be one layer that level is minimum), strengthens Laplce with described initial layers
Pyramidal gradient, as the Graded factor of anisotropic diffusion filtering formula, merges the enhancing of described initial layers and draws
This pyramid diagram picture of the pula fused images with high a layer obtains the fused images of initial layers.Use anisotropy
Described in diffusing filter formula manipulation, the fused images of initial layers obtains final panoramic dental image.Wherein, respectively to
Anisotropic diffusion Filtering Formula is as follows:
Wherein, i is the number of plies decomposing pyramid diagram picture, Ti(m+x, n+y) is the fused images of i-th layer,
G (m, n) is convolution kernel,It it is the gradient of the Laplacian-pyramid image of i-th layer.
The technical scheme of the present embodiment, by the current axially face picture at cone beam computed tomography equipment
In determine Seed Points, utilize curve fitting algorithm process described Seed Points, to generate arch wire.Again with described
Curved surface centered by the curved surface at arch wire place, calculates the picture on described central surface according to the sampling thickness set
The tissue density CT value that vegetarian refreshments is corresponding, generates initial dentistry according to described tissue density CT value complete
Scape image.Then, it is determined that the Graded factor of anisotropic diffusion filtering formula, described anisotropy is used to expand
Dissipate Filtering Formula and process the initial panoramic dental image of different scale, to generate final panoramic dental image.This
The technical scheme of embodiment solves the panorama sketch that can not be obtained denture at present by CBCT three-dimensional data
The problem of picture, it is achieved based on CBCT three-dimensional data, obtain denture by image processing algorithm
The purpose of two-dimensional panoramic image, has reached both to be obtained the three-dimensional structure image of Buccodental by CBCT,
The effect of the two-dimensional panoramic image of denture can be obtained again.
Embodiment two
Fig. 2 a is the side based on CBCT video generation panoramic dental image of the another kind in the embodiment of the present invention two
The flow chart of method.The technical scheme of the present embodiment is on the basis of the technical scheme of above-described embodiment, to determining
The Graded factor of anisotropic diffusion filtering formula, uses described anisotropic diffusion filtering formula manipulation different
The initial panoramic dental image of yardstick, is carried out further generating the step of final panoramic dental image
Bright, specifically include following steps:
Step 210, image pyramid method is utilized to be decomposed into multiple dimensioned by described initial panoramic dental image
Decompose pyramid diagram picture.
Wherein, described image pyramid be with multiresolution come interpretation of images a kind of effectively but concept is simple
Structure.The pyramid of piece image is a series of images gradually reduced with the resolution ratio of Pyramid arrangement
Set.Image pyramid generally comprises two steps: 1, utilize low pass filter smoothed image；2, right
Smoothed image is sampled, thus obtains the image of a series of size reduction.For example, it is possible to select 5 × 5
Gaussian convolution core is (it addition, can also select Gaussian convolution core or the Gauss of 7 × 7 of 3 × 3 according to actual needs
Convolution kernel etc.) current initial panoramic dental image is smoothed.Then, obtain after smoothing processing
To image carry out even number line sampling, thus obtain the decomposition pyramid of the different layers (different scale) wanted
Image (G0、G1、G2、…、Gk-2、Gk-1、Gk).Formula in decomposable process is as follows:
Wherein, W (m, n) be 5 × 5 Gaussian convolution core, GLIt is that L layer decomposes pyramid diagram picture, L value
It is 0,1,2,3 ..., k, k are the top layer level number decomposing pyramid diagram picture.
Step 220, employing interpolation method decompose, to high one layer, the inversion process that pyramid diagram picture performs to decompose,
Be expanded pyramid diagram picture.
Decompose even number row and column in pyramid diagram picture at high one layer and insert 0, then, use above-mentioned decomposable process
The Gaussian convolution core of the 5 × 5 of middle employing is filtered processing, and obtains and the figure of current layer decomposition pyramid diagram picture
The extension pyramid diagram picture that shape size is identical.Respectively each layer is decomposed pyramid diagram picture according to above-mentioned interpolation method to hold
The inversion process that row decomposes, obtains k-1 layer extension pyramid diagram picture (Gt0、Gt1、Gt2、…Gt(k-1))。
Formula in expansion process is as follows:
Wherein, W (m, n) be 5 × 5 Gaussian convolution core, GtbIt is that b layer extends pyramid diagram picture, G(b+1)For
B+1 layer decomposes pyramid diagram picture, and b value is 0, and 1,2,3 ..., k-1, k are the top layer decomposing pyramid diagram picture
Level number.
Step 230, determine that described current layer decomposes pyramid diagram picture and the extension pyramid diagram picture of correspondence respectively
Pixel value difference, determine the described initial panoramic dental image high frequency at different scale according to described pixel value difference
Parts of images matrix.
As shown in Figure 2 b, G is calculated0Layer decomposes pyramid diagram picture and Gt0The pixel of layer extension pyramid diagram picture is poor,
Using the pattern matrix that includes pixel value difference as L0Layer Laplacian-pyramid image.Similarly, G is calculated1
Layer decomposes pyramid diagram picture and Gt1The pixel of layer extension pyramid diagram picture is poor, will include the figure square of pixel value difference
Battle array is as L1Layer Laplacian-pyramid image.Use same method can be calculated L2Layer Laplce
Pyramid diagram picture, L3Layer Laplacian-pyramid image ..., Lk-2Layer Laplacian-pyramid image and Lk-1
Layer Laplacian-pyramid image.Laplacian-pyramid image is i.e. that described initial panoramic dental image is not
HFS image array with yardstick.
Step 240, described HFS image array is carried out strengthen contrast process, obtain HFS
Image enhaucament matrix, the gradient of described HFS image enhaucament matrix is public as anisotropic diffusion filtering
The Graded factor of formula.
First, the described Laplacian-pyramid image of different scale is performed enhancing process respectively, obtains not
Enhancing Laplacian-pyramid image with yardstick.By modes such as linear transformation or curve conversion, respectively
The laplacian pyramid making different scale is multiplied by the curve of setting, to realize described HFS image moment
Battle array (image detail) carries out the purpose strengthened, thus, obtain the enhancing Laplce's gold word under different scale
Tower image.The formula of above-mentioned enhancing process is: Lep=Enhance (Lp), wherein, LepIncrease for pth layer
Strong Laplacian-pyramid image, Enhance is for strengthening operator, LpFor pth layer laplacian pyramid figure
Picture, the value of p is 0,1,2,3 ..., k-1, k are the top layer level number decomposing pyramid diagram picture.As shown in Figure 2 b,
L0The layer enhanced process of Laplacian-pyramid image obtains Le0Layer strengthens Laplacian-pyramid image.With
Sample ground, L1The layer enhanced process of Laplacian-pyramid image obtains Le1Layer strengthens laplacian pyramid figure
Picture.Same method is used to obtain Le2Layer strengthens Laplacian-pyramid image, Le3Layer strengthens Laplce
Pyramid diagram picture ..., Le(k-2)Layer strengthens Laplacian-pyramid image and Le(k-1)Layer strengthens Laplce
Pyramid diagram picture.
Secondly, the gradient strengthening Laplacian-pyramid image of each yardstick is determined, using described gradient as respectively
Graded factor in Anisotropic diffusion Filtering Formula.
The enhancing La Pula of each yardstick of calculated for pixel values that Laplacian-pyramid image includes is strengthened according to each layer
The gradient of this pyramid diagram picture.Such as, strengthen, with current layer, arbitrary picture that Laplacian-pyramid image includes
Point centered by vegetarian refreshments, calculates the zone leveling gradient in the setting regions corresponding with this pixel, thus, can
To obtain the gradient of current layer enhancing Laplacian-pyramid image.Current layer is strengthened laplacian pyramid
The gradient of image is as the Graded factor that current layer fused images carries out anisotropic diffusion filtering process.Its
In, described current layer fused images can be expressed as: Tk=Up (Dk+1)+Lek, wherein, TkFor kth
The fused images of layer, Up represents the operator of the inverse process decomposed, Dk+1Decomposition dental panoramic for kth+1 layer
Image, LekLaplacian-pyramid image is strengthened for kth layer.
Step 250, use the initial dental panoramic of described anisotropic diffusion filtering formula manipulation different scale
Image, to generate final panoramic dental image.
Maximum layer extends pyramid diagram picture merge with maximum layer enhancing Laplacian-pyramid image
Process, obtain maximum layer fused images.As shown in Figure 2 b, to GkLayer decomposes pyramid diagram picture and expands
Exhibition processes, the extension pyramid diagram picture after processing and Le(k-1)Layer strengthens laplacian pyramid and merges,
Obtain fused images T for kth-1 layerk-1(m+x, n+y).
The gradient of Laplacian-pyramid image is strengthened as Graded factor using corresponding described maximum layer, logical
Cross maximum layer fused images described in described anisotropic diffusion filtering formula manipulation, obtain dividing of maximum layer
Solve panoramic dental image.Above-mentioned calculating process can be expressed as: Dk-1=Filter (Gt(k-1)+Le(k-1)),
Wherein, Dk-1For the decomposition panoramic dental image of kth-1 layer, Gt(k-1)Extension dental panoramic for kth-1 layer
Image, Le(k-1)Strengthening Laplacian-pyramid image for kth-1 layer, k is the top layer decomposing pyramid diagram picture
Level number, Filter represents anisotropic diffusion filtering process.
Successively using the decomposition panoramic dental image of high a layer as present image, strengthen Laplce by low one layer
Pyramid diagram picture and described present image carry out fusion treatment, strengthen Laplce with corresponding described low one layer
The gradient of pyramid diagram picture is as Graded factor, after being merged by described anisotropic diffusion filtering formula manipulation
Image, obtain the decomposition panoramic dental image of low a layer, until described Graded factor is minimum one layer of enhancing
Laplacian-pyramid image, the decomposition panoramic dental image of corresponding low one layer is final panoramic dental image.
As shown in Figure 2 b, with Dk-1Layer decomposition panoramic dental image is as present image, by Le(k-2)Layer strengthens
Laplacian-pyramid image and described present image carry out fusion treatment, obtain Tk-2Layer fused images.With
Le(k-2)The gradient of layer enhancing Laplacian-pyramid image, as Graded factor, is expanded by described anisotropy
Dissipate Filtering Formula and process Tk-2Layer fused images, obtains Dk-2Layer decomposes panoramic dental image.By that analogy,
With D2Layer decomposition panoramic dental image is as present image, by Le1Layer strengthen Laplacian-pyramid image with
Described present image carries out fusion treatment, obtains T1Layer fused images.With Le1Layer strengthens laplacian pyramid
The gradient of image is as Graded factor, by described anisotropic diffusion filtering formula manipulation T1Layer fused images,
Obtain D1Layer decomposes panoramic dental image.Finally, with D1Layer decompose panoramic dental image as present image,
By Le0Layer strengthens Laplacian-pyramid image and carries out fusion treatment with described present image, obtains T0Layer merges
Image.With Le0The gradient of layer enhancing Laplacian-pyramid image is as Graded factor, by described each to different
Property diffusing filter formula manipulation T0Layer fused images, obtains D0Layer decomposes panoramic dental image.Due to Le0Layer
Strengthening Laplacian-pyramid image is minimum one layer of enhancing Laplacian-pyramid image, corresponding D0Layer point
Solving panoramic dental image is final panoramic dental image.
Above-mentioned calculating process can be expressed as: DN=Filter (Up (DN+1)+LeN), wherein, DNIt is
The decomposition panoramic dental image of N layer, Up represents the operator of the inverse process decomposed, LeNStrengthen for n-th layer and draw
This pyramid diagram picture of pula, the value of N is 0,1,2,3 ..., k-2, wherein, k is the top decomposing pyramid diagram picture
Layer by layer number, Filter represents anisotropic diffusion filtering process.
The technical scheme of the present embodiment, by using the ladder strengthening Laplacian-pyramid image of different scale
Spend the Graded factor as Filtering Formula, the image after merging is carried out anisotropic diffusion filtering process, real
Show while weakening noise, retained the marginal information of image, enhance the figure of final panoramic dental image
Picture element amount.
Embodiment three
Fig. 3 is a kind of based on CBCT video generation panoramic dental image the device in the embodiment of the present invention three
Structural representation.Described device specifically includes:
Arch wire signal generating unit 310, for the current axially face picture at cone beam computed tomography equipment
In determine Seed Points, utilize curve fitting algorithm process described Seed Points, to generate arch wire；
First image generation unit 320, is used for curved surface centered by the curved surface at described arch wire place, according to
The sampling thickness set calculates the tissue density CT value that the pixel on described central surface is corresponding, root
Initial panoramic dental image is generated according to described tissue density CT value；
Second image generation unit 330, for determining the Graded factor of anisotropic diffusion filtering formula, adopts
With the initial panoramic dental image of described anisotropic diffusion filtering formula manipulation different scale, final to generate
Panoramic dental image.
The technical scheme of the present embodiment, is set at cone beam computed tomography by arch wire signal generating unit 310
Standby current axially face picture determines Seed Points, utilizes curve fitting algorithm to process described Seed Points, with life
Become arch wire.Operations described below is performed with described arch wire place again by the first image generation unit 320
Curved surface centered by curved surface, calculates, according to the sampling thickness set, the people that the pixel on described central surface is corresponding
Soma's density CT value, generates initial panoramic dental image according to described tissue density CT value.Then,
Determined the Graded factor of anisotropic diffusion filtering formula by the second image generation unit 330, use described
The initial panoramic dental image of anisotropic diffusion filtering formula manipulation different scale is complete to generate final dentistry
Scape image.The technical scheme of the present embodiment solves and can not obtain dental pattern by CBCT three-dimensional data at present
The problem of the panoramic picture of tooth, it is achieved based on CBCT three-dimensional data, obtained by image processing algorithm
The purpose of the two-dimensional panoramic image of denture, has reached both to be obtained the three of Buccodental by CBCT
Dimension structural images, can obtain again the effect of the two-dimensional panoramic image of denture.
Further, the first image generation unit 320 specifically for:
Order determines that the width in described axial face image is current axially face picture, hanging down arch wire place
The straight curved surface in current axially face picture, as central surface, determines in described current axially face picture and is positioned at institute
State the pixel on central surface.
Order obtains in described pixel as current pixel point, determines and passes through on described central surface
The normal of described current pixel point.
Calculate and on described normal, meet the average of tissue density CT value corresponding to the pixel imposed a condition
Value, using described mean value as tissue density CT corresponding to the current pixel point on described central surface
Value, wherein, described in the distance including pixel and described current pixel point that imposes a condition thick in the sampling set
In the range of degree, and tissue density CT value corresponding to pixel is in the threshold range set.
According to the stacking order of described axial face image, by tissue corresponding for the pixel on central surface
Density CT value is combined obtaining described initial panoramic dental image.
Further, the second image generation unit 330 includes:
Decompose subelement, be used for utilizing image pyramid method to be decomposed into many by described initial panoramic dental image
The decomposition pyramid diagram picture of yardstick.
Extension subelement, for using interpolation method to decompose, to high one layer, the inverse transformation that pyramid diagram picture performs to decompose
Processing, be expanded pyramid diagram picture, wherein, and described extension pyramid diagram picture and low one layer of decomposition pyramid
The dimension of picture of image is identical.
Image array determines subelement, for determine respectively described current layer decompose pyramid diagram as and correspondence
According to described pixel value difference, the pixel value difference of extension pyramid diagram picture, determines that described initial panoramic dental image exists
The HFS image array of different scale.
Graded factor determines subelement, processes for described HFS image array carries out enhancing contrast,
Obtain HFS image enhaucament matrix, using the gradient of described HFS image enhaucament matrix as each to different
The Graded factor of property diffusing filter formula.
Further, described image pyramid is gaussian pyramid, and described HFS image array is general for drawing
Lars pyramid diagram picture；
And, Graded factor determines that subelement is specifically for described laplacian pyramid to different scale
Image performs enhancing respectively and processes, and obtains the enhancing Laplacian-pyramid image of different scale；Determine each chi
The gradient strengthening Laplacian-pyramid image of degree, using described gradient as anisotropic diffusion filtering formula
In Graded factor.
Further, described second image generation unit 330 also includes:
Maximum layer fused images determines subelement, for maximum layer extending pyramid diagram picture with the highest by one
Layer strengthens Laplacian-pyramid image and carries out fusion treatment, obtains maximum layer fused images.
The decomposition panoramic dental image of maximum layer determines subelement, for increasing with corresponding described maximum layer
The gradient of strong Laplacian-pyramid image is as Graded factor, by described anisotropic diffusion filtering formula
Process described maximum layer fused images, obtain the decomposition panoramic dental image of maximum layer.
Final panoramic dental image determines subelement, for successively with the decomposition panoramic dental image work of high a layer
For present image, strengthen Laplacian-pyramid image by low one layer and carry out fusion treatment with described present image,
Using corresponding described low one layer of gradient strengthening Laplacian-pyramid image as Graded factor, by described
Image after the fusion of anisotropic diffusion filtering formula manipulation, obtains the decomposition panoramic dental image of low a layer,
Until described Graded factor is minimum one layer strengthens Laplacian-pyramid image, the decomposition of corresponding low one layer
Panoramic dental image is final panoramic dental image.
Above-mentioned device based on CBCT video generation panoramic dental image can perform any embodiment of the present invention institute
The method based on CBCT video generation panoramic dental image provided, possesses execution method corresponding function mould
Block and beneficial effect.
Note, above are only presently preferred embodiments of the present invention and institute's application technology principle.Those skilled in the art
It will be appreciated that the invention is not restricted to specific embodiment described here, can enter for a person skilled in the art
Row various obvious changes, readjust and substitute without departing from protection scope of the present invention.Therefore, though
So by above example, the present invention is described in further detail, but the present invention be not limited only to
Upper embodiment, without departing from the inventive concept, it is also possible to include other Equivalent embodiments more,
And the scope of the present invention is determined by scope of the appended claims.
Claims (10)
1. a method based on CBCT video generation panoramic dental image, it is characterised in that including:
In the current axially face picture of cone beam computed tomography equipment, determine Seed Points, utilize curve to intend
Hop algorithm processes described Seed Points, to generate arch wire；
Curved surface centered by the curved surface at described arch wire place, calculates described center according to the sampling thickness set
The tissue density CT value that pixel on curved surface is corresponding, generates according to described tissue density CT value
Initially panoramic dental image；
Determine the Graded factor of anisotropic diffusion filtering formula, use described anisotropic diffusion filtering formula
Process the initial panoramic dental image of different scale, to generate final panoramic dental image.
Method the most according to claim 1, it is characterised in that calculate institute according to the sampling thickness set
State the tissue density CT value that the pixel on central surface is corresponding, according to described tissue density CT
Value generates initial panoramic dental image, including:
Order determines that the width in described axial face image is current axially face picture, hanging down arch wire place
The straight curved surface in current axially face picture, as central surface, determines in described current axially face picture and is positioned at institute
State the pixel on central surface；
Order obtains in described pixel as current pixel point, determines and passes through on described central surface
The normal of described current pixel point；
Calculate and on described normal, meet the average of tissue density CT value corresponding to the pixel imposed a condition
Value, using described mean value as tissue density CT corresponding to the current pixel point on described central surface
Value, wherein, described in the distance including pixel and described current pixel point that imposes a condition thick in the sampling set
In the range of degree, and tissue density CT value corresponding to pixel is in the threshold range set；
According to the stacking order of described axial face image, by tissue corresponding for the pixel on central surface
Density CT value is combined obtaining described initial panoramic dental image.
Method the most according to claim 1, it is characterised in that determine anisotropic diffusion filtering formula
Graded factor, including:
Utilize image pyramid method that described initial panoramic dental image is decomposed into multiple dimensioned decomposition pyramid
Image；
Using interpolation method to decompose, to high one layer, the inversion process that pyramid diagram picture performs to decompose, be expanded gold
Word tower image, wherein, described extension pyramid diagram picture with low one layer decompose pyramid diagram picture dimension of picture phase
With；
Determine that described current layer decomposes pyramid diagram picture and the pixel value difference of corresponding extension pyramid diagram picture respectively,
The described initial panoramic dental image HFS image array at different scale is determined according to described pixel value difference；
Described HFS image array carries out strengthen contrast process, obtain HFS image enhaucament square
Battle array, using the gradient of described HFS image enhaucament matrix as anisotropic diffusion filtering formula gradient because of
Son.
Method the most according to claim 3, it is characterised in that described image pyramid is Gauss gold word
Tower, described HFS image array is Laplacian-pyramid image；
And, using the gradient of described HFS image enhaucament matrix as anisotropic diffusion filtering formula
Graded factor, including:
The described Laplacian-pyramid image of different scale is performed enhancing process respectively, obtains different scale
Enhancing Laplacian-pyramid image；
Determine the gradient strengthening Laplacian-pyramid image of each yardstick, using described gradient as anisotropy
Graded factor in diffusing filter formula.
Method the most according to claim 4, it is characterised in that use described anisotropic diffusion filtering
The initial panoramic dental image of formula manipulation different scale, to generate final panoramic dental image, including:
Maximum layer extends pyramid diagram picture merge with maximum layer enhancing Laplacian-pyramid image
Process, obtain maximum layer fused images；
The gradient of Laplacian-pyramid image is strengthened as Graded factor using corresponding described maximum layer, logical
Cross maximum layer fused images described in described anisotropic diffusion filtering formula manipulation, obtain dividing of maximum layer
Solve panoramic dental image；
Successively using the decomposition panoramic dental image of high a layer as present image, strengthen Laplce by low one layer
Pyramid diagram picture and described present image carry out fusion treatment, strengthen Laplce with corresponding described low one layer
The gradient of pyramid diagram picture is as Graded factor, after being merged by described anisotropic diffusion filtering formula manipulation
Image, obtain the decomposition panoramic dental image of low a layer, until described Graded factor is minimum one layer of enhancing
Laplacian-pyramid image, the decomposition panoramic dental image of corresponding low one layer is final panoramic dental image.
6. a device based on CBCT video generation panoramic dental image, it is characterised in that including:
Arch wire signal generating unit, for true in the current axially face picture of cone beam computed tomography equipment
Determine Seed Points, utilize curve fitting algorithm to process described Seed Points, to generate arch wire；
First image generation unit, for curved surface centered by the curved surface at described arch wire place, according to setting
Sampling thickness calculate the tissue density CT value that the pixel on described central surface is corresponding, according to institute
State tissue density CT value and generate initial panoramic dental image；
Second image generation unit, for determining the Graded factor of anisotropic diffusion filtering formula, uses institute
State the initial panoramic dental image of anisotropic diffusion filtering formula manipulation different scale, to generate final dentistry
Panoramic picture.
Device the most according to claim 6, it is characterised in that the first image generation unit specifically for:
Order determines that the width in described axial face image is current axially face picture, hanging down arch wire place
The straight curved surface in current axially face picture, as central surface, determines in described current axially face picture and is positioned at institute
State the pixel on central surface；
Order obtains in described pixel as current pixel point, determines and passes through on described central surface
The normal of described current pixel point；
Calculate and on described normal, meet the average of tissue density CT value corresponding to the pixel imposed a condition
Value, using described mean value as tissue density CT corresponding to the current pixel point on described central surface
Value, wherein, described in the distance including pixel and described current pixel point that imposes a condition thick in the sampling set
In the range of degree, and tissue density CT value corresponding to pixel is in the threshold range set；
According to the stacking order of described axial face image, by tissue corresponding for the pixel on central surface
Density CT value is combined obtaining described initial panoramic dental image.
Device the most according to claim 6, it is characterised in that the second image generation unit includes:
Decompose subelement, be used for utilizing image pyramid method to be decomposed into many by described initial panoramic dental image
The decomposition pyramid diagram picture of yardstick；
Extension subelement, for using interpolation method to decompose, to high one layer, the inverse transformation that pyramid diagram picture performs to decompose
Processing, be expanded pyramid diagram picture, wherein, and described extension pyramid diagram picture and low one layer of decomposition pyramid
The dimension of picture of image is identical；
Image array determines subelement, for determine respectively described current layer decompose pyramid diagram as and correspondence
According to described pixel value difference, the pixel value difference of extension pyramid diagram picture, determines that described initial panoramic dental image exists
The HFS image array of different scale；
Graded factor determines subelement, processes for described HFS image array carries out enhancing contrast,
Obtain HFS image enhaucament matrix, using the gradient of described HFS image enhaucament matrix as each to different
The Graded factor of property diffusing filter formula.
Device the most according to claim 8, it is characterised in that described image pyramid is Gauss gold word
Tower, described HFS image array is Laplacian-pyramid image；
And, Graded factor determine subelement specifically for:
The described Laplacian-pyramid image of different scale is performed enhancing process respectively, obtains different scale
Enhancing Laplacian-pyramid image；
Determine the gradient strengthening Laplacian-pyramid image of each yardstick, using described gradient as anisotropy
Graded factor in diffusing filter formula.
Device the most according to claim 9, it is characterised in that described second image generation unit is also
Including:
Maximum layer fused images determines subelement, for maximum layer extending pyramid diagram picture with the highest by one
Layer strengthens Laplacian-pyramid image and carries out fusion treatment, obtains maximum layer fused images；
The decomposition panoramic dental image of maximum layer determines subelement, for increasing with corresponding described maximum layer
The gradient of strong Laplacian-pyramid image is as Graded factor, by described anisotropic diffusion filtering formula
Process described maximum layer fused images, obtain the decomposition panoramic dental image of maximum layer；
Final panoramic dental image determines subelement, for successively with the decomposition panoramic dental image work of high a layer
For present image, strengthen Laplacian-pyramid image by low one layer and carry out fusion treatment with described present image,
Using corresponding described low one layer of gradient strengthening Laplacian-pyramid image as Graded factor, by described
Image after the fusion of anisotropic diffusion filtering formula manipulation, obtains the decomposition panoramic dental image of low a layer,
Until described Graded factor is minimum one layer strengthens Laplacian-pyramid image, the decomposition of corresponding low one layer
Panoramic dental image is 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 true CN105894444A (en) | 2016-08-24 |
| CN105894444B 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) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107301622A (en) * | 2017-06-20 | 2017-10-27 | 广州华端科技有限公司 | The panoramic picture acquisition methods and system of oral cavity CBCT images |
| CN109377534A (en) * | 2018-10-22 | 2019-02-22 | 南方医科大学 | The non-linear oral cavity CT panoramic image synthesis method of automatic sampling thickness detection |
| CN109389652A (en) * | 2017-08-03 | 2019-02-26 | 奥齿泰因普兰特株式会社 | Filter usability methods and system for dental CT scan image |
| CN110211200A (en) * | 2019-04-22 | 2019-09-06 | 深圳安科高技术股份有限公司 | A kind of arch wire generation method and its system based on nerual network technique |
| CN110335281A (en) * | 2018-03-28 | 2019-10-15 | 北京连心医疗科技有限公司 | A kind of tumor boundaries determine method, equipment and storage medium |
| CN112102495A (en) * | 2020-09-15 | 2020-12-18 | 北京朗视仪器有限公司 | Dental arch curved surface generation method based on CBCT image |
| CN112270643A (en) * | 2020-09-04 | 2021-01-26 | 深圳市菲森科技有限公司 | Three-dimensional imaging data splicing method and device, electronic equipment and storage medium |
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 |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107301622A (en) * | 2017-06-20 | 2017-10-27 | 广州华端科技有限公司 | The panoramic picture acquisition methods and system of oral cavity CBCT images |
| CN107301622B (en) * | 2017-06-20 | 2018-04-06 | 广州华端科技有限公司 | The panoramic picture acquisition methods and system of oral cavity CBCT images |
| CN109389652A (en) * | 2017-08-03 | 2019-02-26 | 奥齿泰因普兰特株式会社 | Filter usability methods and system for dental CT scan image |
| CN109389652B (en) * | 2017-08-03 | 2023-10-13 | 奥齿泰因普兰特株式会社 | Filter application method and system for dental computerized tomography image |
| CN110335281A (en) * | 2018-03-28 | 2019-10-15 | 北京连心医疗科技有限公司 | A kind of tumor boundaries determine method, equipment and storage medium |
| CN109377534A (en) * | 2018-10-22 | 2019-02-22 | 南方医科大学 | The non-linear oral cavity CT panoramic image synthesis method of automatic sampling thickness detection |
| CN109377534B (en) * | 2018-10-22 | 2022-04-26 | 南方医科大学 | Nonlinear oral cavity CT panoramic image synthesis method capable of automatically sampling thickness detection |
| CN110211200A (en) * | 2019-04-22 | 2019-09-06 | 深圳安科高技术股份有限公司 | A kind of arch wire generation method and its system based on nerual network technique |
| CN112270643A (en) * | 2020-09-04 | 2021-01-26 | 深圳市菲森科技有限公司 | Three-dimensional imaging data splicing method and device, electronic equipment and storage medium |
| CN112102495A (en) * | 2020-09-15 | 2020-12-18 | 北京朗视仪器有限公司 | Dental arch curved surface generation method based on CBCT image |
| CN112102495B (en) * | 2020-09-15 | 2023-11-21 | 北京朗视仪器股份有限公司 | Dental arch surface generation method based on CBCT image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105894444B (en) | 2018-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105894444A (en) | Method and device for generating dental panoramic image on the basis of CBCT image | |
| Segars et al. | Modeling respiratory mechanics in the MCAT and spline-based MCAT phantoms | |
| CN104318536B (en) | The bearing calibration of CT images and device | |
| CN103971403B (en) | The method and system being imaged for Curve Reconstruction | |
| EP1522045B1 (en) | Motion artifact correction of tomographical images | |
| CN102715906B (en) | Method and system for 3D cardiac motion estimation from single scan of c-arm angiography | |
| CN106725570A (en) | Imaging method and system | |
| CN106600609A (en) | Spine partition method and system in medical image | |
| CN106683144A (en) | Image iteration reconstruction method and device | |
| CN106296615B (en) | The method and system of metal artifacts is corrected in CT images | |
| CN106934821A (en) | A kind of conical beam CT and CT method for registering images based on ICP algorithm and B-spline | |
| CN105118039B (en) | Realize the method and system that pyramidal CT image is rebuild | |
| Zhou et al. | DuDoUFNet: dual-domain under-to-fully-complete progressive restoration network for simultaneous metal artifact reduction and low-dose CT reconstruction | |
| CN102024251A (en) | System and method for multi-image based virtual non-contrast image enhancement for dual source CT | |
| CN112669401B (en) | CT image reconstruction method and system based on convolutional neural network | |
| CN111612689B (en) | Medical image processing method, medical image processing device, computer equipment and readable storage medium | |
| Uccheddu et al. | 3D printing of cardiac structures from medical images: an overview of methods and interactive tools | |
| Sano et al. | Super-resolution method and its application to medical image processing | |
| Bayaraa et al. | A two-stage approach for beam hardening artifact reduction in low-dose dental CBCT | |
| CN107240140A (en) | A kind of method and device of image reconstruction | |
| Rohkohl et al. | Interventional 4-D motion estimation and reconstruction of cardiac vasculature without motion periodicity assumption | |
| CN110473297A (en) | Image processing method, device, electronic equipment and storage medium | |
| DE102007009019B3 (en) | A method of providing advanced capabilities in the imaging of a patient's heart | |
| US8571171B2 (en) | Reconstructing a tomographic image with reduced artifacts | |
| KR20210069389A (en) | Apparatus and method for removing metal artifact of computer tomography image based on artificail intelligence |
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 |