CN107657647A - A kind of taper beam X-ray oral cavity fault image algorithm - Google Patents
A kind of taper beam X-ray oral cavity fault image algorithm Download PDFInfo
- Publication number
- CN107657647A CN107657647A CN201710863873.9A CN201710863873A CN107657647A CN 107657647 A CN107657647 A CN 107657647A CN 201710863873 A CN201710863873 A CN 201710863873A CN 107657647 A CN107657647 A CN 107657647A
- Authority
- CN
- China
- Prior art keywords
- image
- projection
- oral cavity
- value
- ray
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- 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/10116—X-ray image
-
- 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)
Abstract
The invention discloses a kind of taper beam X-ray oral cavity fault image algorithm, comprise the following steps：（1）Sliding-model control is carried out to the image that needs are rebuild using grid, and assigns an initial value；（2）In initial valueOn the basis of to image carry out all angles under projection, its projection simulation numerical value is calculated, and establish projection coefficient matrix W；（3）In the mutual restriction range of POCS method for reconstructing, to step（2）In the projection simulation numerical value that calculates difference is drawn compared with its real projection values under equal angular, and amendment is iterated to the pixel value of each grid of image according to this difference；（4）Repeat step（3）Terminate to iteration, obtain reconstruction image.The present invention realizes Enhanced Imaging effect and three-dimensional reconstruction effect by being improved to existing iterative algorithm.
Description
Technical field
The present invention relates to medical imaging field, and in particular to a kind of taper beam X-ray oral cavity fault image algorithm.
Background technology
X-ray tomography imaging is the attenuation characteristic using ray to material, nondestructively obtains internal structure of body information
A kind of technological means.Oral cavity CT has many superiority as a kind of detection technique of hurtless measure, and the technology is in medical diagnosis
It is widely applied with industrial nondestructive testing.In recent years, cone-beam obtains more and more extensive application in medical science and industrially,
Have the characteristics that high accuracy, roentgen dose X are small, sweep speed is fast, turn into one of future technical advances trend.With dentistry
Development, the special cone-beam shape of dentistry is applied widely.Oral cavity can be from different angles to organizing reflect, Neng Goufa
Now trickleer lesion, bone tissue situation, lower song articular conditions are accurately evaluated, for assisting a physician, progress is more accurate
Clinical diagnosis played important function.
The core of taper beam X-ray oral cavity fault image is exactly algorithm for reconstructing.Imaging algorithm master for oral cavity CT equipment
It is divided into two classes：Analytic reconstruction algorithm and algebraic reconstruction algorithm.The speed that analytic method amount of calculation during reconstruction is smaller, rebuilds
Degree is also very fast, and the reconstruction image of high quality can be obtained in the case where obtaining complete and uniform data for projection.But
Not exclusively or in the case of strong noise, the algorithm is no longer applicable data for projection.And during practical application, due in the presence of one
The reason for a little objective, tend not to enough meet its requirement, the data for projection of collection is often uneven or imperfect.Iterative algorithm
There is certain advantage in terms of for handling incomplete data for projection, preferable figure can be reconstructed under less data for projection
As quality, and also there is certain suppression to noise.
The content of the invention
The goal of the invention of the present invention is to provide a kind of taper beam X-ray oral cavity fault image algorithm, by existing iteration
Algorithm is improved, and realizes Enhanced Imaging effect and three-dimensional reconstruction effect.
To achieve the above object of the invention, the technical solution adopted by the present invention is：A kind of taper beam X-ray oral cavity fault image
Algorithm, comprise the following steps：
（1）Sliding-model control is carried out to the image that needs are rebuild using grid, and assigns an initial value；
（2）In initial valueOn the basis of to image carry out all angles under projection, its projection simulation numerical value is calculated, and build
Vertical projection coefficient matrix W；
（3）In the mutual restriction range of POCS method for reconstructing, to step（2）In the projection simulation numerical value that calculates with it in phase
It is compared with the real projection values under angle and draws difference, and the pixel value of each grid of image is entered according to this difference
Row iteration amendment；
（4）Repeat step（3）Terminate to iteration, obtain reconstruction image.
Above, when in view of image actual pixel value, makeTake 255 when more than 255, less than 0 when
Time takes 0.
Preferably, the step（1）In, initial valueTake null vector.
Preferably, the step（3）In the formula of iterated revision be,
Wherein,To rebuild image vector component,For+1 iteration result of kth,For kth time iteration result,For reality
The projection value of border measurement,The line segment length intercepted by i bars ray in n pixel,To be intercepted in n pixel
Line segment,For weighted factor.
Because above-mentioned technical proposal is used, the present invention has following advantages compared with prior art：
1. the present invention is rebuild using iterative reconstruction algorithm to image, suitable for the image reconstruction of incomplete data for projection, and
The picture quality reconstructed is high, and has certain inhibitory action to noise；
2. the present invention can effectively solve the problem that two dimensional image Problems of Reconstruction, and can extend to three dimensions.
Brief description of the drawings
Fig. 1 is the schematic diagram of the embodiment of the present invention one.
Embodiment
Below in conjunction with the accompanying drawings and embodiment the invention will be further described：
Embodiment one：
It is shown in Figure 1, a kind of taper beam X-ray oral cavity fault image algorithm, comprise the following steps：
（1）Sliding-model control is carried out to the image that needs are rebuild using grid, and assigns an initial value；
（2）In initial valueOn the basis of to image carry out all angles under projection, its projection simulation numerical value is calculated, and build
Vertical projection coefficient matrix W；
（3）In the mutual restriction range of POCS method for reconstructing, to step（2）In the projection simulation numerical value that calculates with it in phase
It is compared with the real projection values under angle and draws difference, and the pixel value of each grid of image is entered according to this difference
Row iteration amendment；
（4）Repeat step（3）Terminate to iteration, obtain reconstruction image.
Above, when in view of image actual pixel value, makeTake 255 when more than 255, less than 0 when
Time takes 0.
In the present embodiment, the step（1）In, initial valueTake null vector.
The step（3）In the formula of iterated revision be,
Wherein,To rebuild image vector component,For+1 iteration result of kth,For kth time iteration result,For reality
The projection value of border measurement,The line segment length intercepted by i bars ray in n pixel,To be intercepted in n pixel
Line segment,For weighted factor.
The present invention iteration theorem be：Arbitrarily one group of initial value of selection, projection fortune is carried out on the basis of the initial value
Calculate, after the amendment compared with real projection value, be met the value of the ray equation in a direction；Then willBased on be worth, then carry out the value that above-mentioned calculating process is met the ray equation in next direction；Constantly repeat
The process stated, after the ray equation on all directions all is calculated into completion, obtain faultage image estimate.Such one
Secondary iteration just completes, and then again willAs initial value, to carry out the interative computation of next round, until meeting iteration stopping
Untill condition, the image for reconstructing and is finally obtained.
The foregoing description of the disclosed embodiments, professional and technical personnel in the field are enable to realize or using the present invention.
A variety of modifications to above-described embodiment will be apparent for those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, it is of the invention
Above-described embodiment shown in this article is not intended to be limited to, and is to fit to and principles disclosed herein and features of novelty phase one
The most wide scope caused.
Claims (3)
1. a kind of taper beam X-ray oral cavity fault image algorithm, it is characterised in that comprise the following steps：
（1）Sliding-model control is carried out to the image that needs are rebuild using grid, and assigns an initial value；
（2）In initial valueOn the basis of to image carry out all angles under projection, its projection simulation numerical value is calculated, and build
Vertical projection coefficient matrix W；
（3）In the mutual restriction range of POCS method for reconstructing, to step（2）In the projection simulation numerical value that calculates with it in phase
It is compared with the real projection values under angle and draws difference, and the pixel value of each grid of image is entered according to this difference
Row iteration amendment；
（4）Repeat step（3）Terminate to iteration, obtain reconstruction image.
2. taper beam X-ray oral cavity fault image algorithm according to claim 1, it is characterised in that：The step（1）
In, initial valueTake null vector.
3. taper beam X-ray oral cavity fault image algorithm according to claim 1, it is characterised in that：The step（3）In
The formula of iterated revision be,
Wherein,To rebuild image vector component,For+1 iteration result of kth,For kth time iteration result,For reality
The projection value of border measurement,The line segment length intercepted by i bars ray in n pixel,To be intercepted in n pixel
Line segment,For weighted factor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710863873.9A CN107657647A (en) | 2017-09-22 | 2017-09-22 | A kind of taper beam X-ray oral cavity fault image algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710863873.9A CN107657647A (en) | 2017-09-22 | 2017-09-22 | A kind of taper beam X-ray oral cavity fault image algorithm |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107657647A true CN107657647A (en) | 2018-02-02 |
Family
ID=61130850
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710863873.9A Withdrawn CN107657647A (en) | 2017-09-22 | 2017-09-22 | A kind of taper beam X-ray oral cavity fault image algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107657647A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109389575A (en) * | 2018-10-09 | 2019-02-26 | 山东理工大学 | A kind of quick partial reconstruction method of image based on Algebraic Iterative Method |
| CN112602119A (en) * | 2018-09-03 | 2021-04-02 | 爱克发有限公司 | Method for generating enhanced tomographic image of object |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102590248A (en) * | 2012-03-13 | 2012-07-18 | 重庆大学 | Method for online detecting electronic element by translational type micro focus CT (Computerized Tomography) detection device |
| CN104021582A (en) * | 2014-05-28 | 2014-09-03 | 山东大学 | CT (Computed Tomography) iterative image reconstruction method |
-
2017
- 2017-09-22 CN CN201710863873.9A patent/CN107657647A/en not_active Withdrawn
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102590248A (en) * | 2012-03-13 | 2012-07-18 | 重庆大学 | Method for online detecting electronic element by translational type micro focus CT (Computerized Tomography) detection device |
| CN104021582A (en) * | 2014-05-28 | 2014-09-03 | 山东大学 | CT (Computed Tomography) iterative image reconstruction method |
Non-Patent Citations (1)
| Title |
|---|
| 曹涵: "口腔CT成像的迭代重建算法研究", 《医药卫生科技辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112602119A (en) * | 2018-09-03 | 2021-04-02 | 爱克发有限公司 | Method for generating enhanced tomographic image of object |
| CN109389575A (en) * | 2018-10-09 | 2019-02-26 | 山东理工大学 | A kind of quick partial reconstruction method of image based on Algebraic Iterative Method |
| CN109389575B (en) * | 2018-10-09 | 2023-12-01 | 山东理工大学 | Algebraic iteration method-based rapid local image reconstruction method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7711087B2 (en) | Patient setup using tomosynthesis techniques | |
| Leng et al. | High temporal resolution and streak-free four-dimensional cone-beam computed tomography | |
| JP6130811B2 (en) | CT apparatus and method based on motion compensation | |
| McClelland et al. | Inter-fraction variations in respiratory motion models | |
| Fleming et al. | Description of pulmonary deposition of radiolabeled aerosol by airway generation using a conceptual three dimensional model of lung morphology | |
| CN110268446A (en) | System and method for real time imagery | |
| US20020127523A1 (en) | Mathematical model and a method and apparatus for utilizing the model | |
| Maier et al. | Fast simulation of x-ray projections of spline-based surfaces using an append buffer | |
| Badea et al. | 4D micro-CT for cardiac and perfusion applications with view under sampling | |
| Segars et al. | Extension of the 4D NCAT phantom to dynamic x-ray CT simulation | |
| CN111260748A (en) | Digital synthesis X-ray tomography method based on neural network | |
| Schultze et al. | Performance of hull-detection algorithms for proton computed tomography reconstruction | |
| CN105488826A (en) | Energy spectrum CT iterative imaging method and system based on EBP | |
| Sano et al. | Super-resolution method and its application to medical image processing | |
| Do et al. | A decomposition-based CT reconstruction formulation for reducing blooming artifacts | |
| TW201514919A (en) | Computed tomography image reconstruction method | |
| CN107657647A (en) | A kind of taper beam X-ray oral cavity fault image algorithm | |
| JP2007198866A (en) | General saddle cone beam ct system, and three-dimensional reconstitution method | |
| WO2003105017A2 (en) | Analysis of a multi-dimensional structure | |
| WO2021059165A1 (en) | Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device | |
| EP3867871A1 (en) | System and method for automatic delineation of scanned images | |
| Yuan et al. | A modified ray tracing method for ultrasound computed tomography in breast imaging | |
| Xiu et al. | An Innovative Beam Hardening Correction Method for Computed Tomography Systems. | |
| Pengpan et al. | Cone Beam CT using motion-compensated algebraic reconstruction methods with limited data | |
| Manescu et al. | Motion compensation for PET image reconstruction using deformable tetrahedral meshes |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180202 |