CN107192726B - The quick high-resolution 3 D cone-beam computer tomography method of plate shell object and device - Google Patents
The quick high-resolution 3 D cone-beam computer tomography method of plate shell object and device Download PDFInfo
- Publication number
- CN107192726B CN107192726B CN201710312113.9A CN201710312113A CN107192726B CN 107192726 B CN107192726 B CN 107192726B CN 201710312113 A CN201710312113 A CN 201710312113A CN 107192726 B CN107192726 B CN 107192726B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- image
- plate shell
- image sequence
- indicates
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/30—Accessories, mechanical or electrical features
- G01N2223/306—Accessories, mechanical or electrical features computer control
Abstract
The invention discloses a kind of quick high-resolution 3 D cone-beam computer tomography method of plate shell object and devices, this method comprises: obtaining the collected two-dimensional projection's sequence of detector;Logarithm demodulation is carried out to the two-dimensional projection image sequence and obtains two-dimensional line integral image sequence;Object dimensional CT sectioning image is obtained using image reconstruction is carried out to two-dimensional line integral image sequence with azimuthal novel three-dimensional filter back-projection reconstruction algorithm.The embodiment of the present invention can effectively reduce reconstruction volume compared to existing plate shell CT scan technology, to improve reconstruction efficiency;The CT faultage image reconstructed, without passing through the operation such as image rotation and space interpolation, the exact matching of coordinate system Yu plate shell object three-dimensional orthogonal anatomical structure can be realized, reach smooth vertical and horizontal visual effect, image processing time is not only saved, reconstructed image quality is improved, and enhances the readability of image.
Description
Technical field
The present embodiments relate to X ray computer chromatography imaging technique field more particularly to a kind of plate shell object are quick
High-resolution 3 D cone-beam computer tomography method and device.
Background technique
In X ray computer tomography (Computed Tomography, abbreviation CT) system, x-ray source issues X
Ray, passes through a certain region of object to be detected from different perspectives, and the detector for being placed in radiographic source opposite connects in respective angles
It receives.Then, the decaying different degrees of according to each angle ray carries out operation using certain algorithm for reconstructing and computer, rebuilds
Object is scanned the ray line attenuation coefficient distribution map image in region out, to realize reconstructing image from projection, nondestructively again
The features such as existing Media density, ingredient and the structural form of object in the area.
Rebuilding efficiency, reconstruction precision and three-dimensional visualization complexity is the imaging indicators that current CT technology is paid close attention to.
It is restricted by image-forming principle, solar energy sailboard, circuit board etc. have large area-thickness ratio plate shell object, traditional CT technology
It is difficult to obtain good reconstructed results.
In recent years, Helfen L etc., On the implementation of computed laminography using
Synchrotron radiation, Review of Scientific Instruments, vol.82, No.6, pp:129-
151,2011, realize a kind of dip sweeping CT imaging technique (Computed specifically for plate and shell structure
Laminography, abbreviation CL), it can preferably solve the above problems.But there is also following limitations for the technology: when inclination turntable is flat
When face is uneven or plate shell object surface projecting feature causes principal plane that can not be bonded with turntable plane, CL slice is difficult to and plate shell
Surface is parallel, leads to that reconstruction time is long, reconstruction precision is low and three-dimensional visualization is difficult.
Currently, it is not yet found that the plate shell X ray computer chromatography imaging technique of the above problem can be overcome.
Summary of the invention
The embodiment of the present invention provides a kind of quick high-resolution 3 D cone-beam computer tomography method of plate shell object and dress
It sets, rebuilds volume to reduce plate shell object CT, improve imaging efficiency and reconstruction precision, enhance image readability.
The one aspect of the embodiment of the present invention is to provide a kind of quick high-resolution 3 D cone-beam Computerized chromatographic of plate shell object
Imaging method includes the following steps:
Step 1 obtains the collected two-dimensional projection image sequence of detector, and the two-dimensional projection image sequence is object edge
In axial rotation process, multiple two dimensions for being obtained after the ray projection of the excessively described object of detector periodicity acquisition of transmission
Projected image, and the two-dimensional projection image of an a sampling period corresponding object;
Step 2 carries out logarithm demodulation acquisition two-dimensional line integral image sequence to the two-dimensional projection image sequence;Utilize band
Azimuthal novel three-dimensional filter back-projection reconstruction algorithm carries out image reconstruction to two-dimensional line integral image sequence, obtains the object
The corresponding three dimensional CT sectioning image of body.
Further, step 1 further include: plate shell object be placed in inclination turntable on, during the scanning process plate shell object with
Turntable rotation and rotate;
The plate shell object is along 360 degree of rotation center uniform rotation of rotary shaft, and imaging area quilt during rotation
Cone-beam covering.
Further, described that logarithm demodulation acquisition two-dimensional line integral image sequence is carried out to the two-dimensional projection image sequence
Column, comprising:
Logarithm demodulation is carried out to the two-dimensional projection image sequence I (i, m, n) according to formula (1) and obtains two-dimensional line integrogram
As sequence p (i, m, n):
Wherein, i indicates scanning angle, and (m, n) indicates two coordinates of the detection channels in the vertical plane of the detector,
Ln indicates natural logrithm operation, and mean indicates that two-dimentional mean operation, 1:10 indicate that getting 10,1:K from 1 indicates to get K, I from 1
(i, 1:10,1:K) indicates that width is 10, is highly the area of K in i-th of two-dimensional projection image in the two-dimensional projection image sequence
Domain.
Further, described using corresponding to the plate shell object with azimuthal novel filter back-projection reconstruction algorithm
Two-dimensional line integral image sequence carries out image reconstruction and obtains the corresponding three dimensional CT sectioning image of the plate shell object, comprising:
Image reconstruction acquisition is carried out to the corresponding two-dimensional line integral image sequence of the plate shell object according to formula (2)-(6)
The three dimensional CT sectioning image f of the plate shell object:
X=(x " cos α+(and y " sin α+z " cos α) sin γ) cos β-(y " cos α-z " sin α) sin β (4)
Wherein, f (x ", y ", z ") indicates that reconstructed three dimensional panel shell object, (x ", y ", z ") indicate that three-dimension object particle exists
Coordinate position in object coordinates system, (x, y, z) indicate three-dimension object particle in the case where some rotates angle in system coordinate system
Coordinate position, D indicates that radiographic source arrives the distance of detector, and p (s, v, β) indicates that the two-dimensional line that obtains accumulates after carrying out logarithm demodulation
Partial image sequence, (s, v) indicate that coordinate position of certain pixel in two dimensional image coordinate system, corresponding three-dimensional object particle are mapped to
Projection address on two-dimensional line integral image, β indicate rotation angle, and h (s) is one-dimensional filtering device, α,γ distinguishes characterization board shell
Angle between three main shafts of object and three-dimensional cartesian rectangular coordinate system.
The other side of the embodiment of the present invention is to provide a kind of quick high-resolution 3 D cone-beam computer layer of plate shell object
Analyse imaging device, comprising:
Module is obtained, for obtaining the collected two-dimensional projection image sequence of detector, the two-dimensional projection image sequence
It is object axially in rotation process, is obtained after the ray projection of the excessively described object of detector periodicity acquisition of transmission
Multiple two-dimensional projection images, and the two-dimensional projection image of an a sampling period corresponding object;
Computing module obtains two-dimensional line integral image sequence for carrying out logarithm demodulation to the two-dimensional projection image sequence
Column;Image reconstruction is carried out to two-dimensional line integral image sequence using with azimuthal novel three-dimensional filter back-projection reconstruction algorithm,
Obtain the corresponding three dimensional CT sectioning image of the object.
Further, module is obtained further include: plate shell object is placed on inclination turntable, during the scanning process plate shell object
It is rotated as turntable rotates;
The plate shell object is along 360 degree of rotation center uniform rotation of rotary shaft, and imaging area quilt during rotation
Cone-beam covering.
Further, the computing module be specifically used for according to formula (1) to the two-dimensional projection image sequence I (i, m,
N) it carries out logarithm demodulation and obtains two-dimensional line integral image sequence p (i, m, n):
Wherein, i indicates scanning angle, and (m, n) indicates two coordinates of the detection channels in the vertical plane of the detector,
Ln indicates natural logrithm operation, and mean indicates that two-dimentional mean operation, 1:10 indicate that getting 10,1:K from 1 indicates to get K, I from 1
(i, 1:10,1:K) indicates that width is 10, is highly the area of K in i-th of two-dimensional projection image in the two-dimensional projection image sequence
Domain.
Further, the computing module is specifically used for the two-dimensional line integrogram according to formula (2)-(6) to the object
As sequence carries out the three dimensional CT sectioning image f that image reconstruction obtains the object:
X=(x " cos α+(and y " sin α+z " cos α) sin γ) cos β-(y " cos α-z " sin α) sin β (4)
Wherein, f (x ", y ", z ") indicates that reconstructed three dimensional panel shell object, (x ", y ", z ") indicate that three-dimension object particle exists
Coordinate position in object coordinates system, (x, y, z) indicate three-dimension object particle in the case where some rotates angle in system coordinate system
Coordinate position, D indicates that radiographic source arrives the distance of detector, and p (s, v, β) indicates that the two-dimensional line that obtains accumulates after carrying out logarithm demodulation
Partial image sequence, (s, v) indicate that coordinate position of certain pixel in two dimensional image coordinate system, corresponding three-dimensional object particle are mapped to
Projection address on two-dimensional line integral image, β indicate rotation angle, and h (s) is one-dimensional filtering device, α,γ distinguishes characterization board shell
Angle between three main shafts of object and three-dimensional cartesian rectangular coordinate system.
The advantages of the present invention over the prior art are that: the present invention is able to solve plate shell object in three-dimensional orthogonal Descartes
Azimuth problem under coordinate system, quick high accuracy directly reconstruct out the sectional slice image for being parallel to plate shell object surface;Step
It is rapid simple, it reduces and rebuilds volume to improve reconstruction efficiency;Object space position is corrected during reconstruction, rather than
Image rotation and space interpolation are carried out to data again after rebuilding, without interpolation error, thus improve reconstructed image quality;
It reconstructs the sectioning image come and shows with plate shell object parallel, legibility is strong, three-dimensional visualization precision height.
Detailed description of the invention
Fig. 1 is the quick high-resolution 3 D cone-beam computer tomography method stream of plate shell object provided in an embodiment of the present invention
Cheng Tu;
Fig. 2 is that the quick high-resolution 3 D cone-beam computerized tomography system of plate shell object provided in an embodiment of the present invention is former
Reason figure;
Fig. 3 is that plate shell object novel three-dimensional filter back-projection reconstruction algorithm provided in an embodiment of the present invention rebuilds coordinate system original
Reason figure;
Fig. 4 is that the actual object obtained using the quick high-resolution 3 D cone-beam computer tomography method of plate shell object is existed
Two-dimensional line integral image under some scanning angles;
Fig. 5 is the three-dimensional after the actual object obtained using plate shell object novel three-dimensional filter back-projection reconstruction algorithm is rebuild
CT image;
Fig. 6 is the three-dimensional CT image that the actual object obtained is rebuild using traditional reconstruction algorithm;
Fig. 7 is the quick high-resolution 3 D cone-beam computer laminated imaging device of plate shell object provided in an embodiment of the present invention
Structure chart.
Specific embodiment
With reference to the accompanying drawing and specific embodiment further illustrates the present invention.
Fig. 1 is the quick high-resolution 3 D cone-beam computer tomography method stream of plate shell object provided in an embodiment of the present invention
Cheng Tu;The embodiment of the present invention for existing plate shell object CT be imaged, sectioning image be difficult to it is parallel with plate shell surface, cause rebuild when
Between it is long, reconstruction precision is low and three-dimensional visualization is difficult the problems such as, provide the plate based on novel filter back-projection reconstruction algorithm
The quick high-resolution 3 D cone-beam computer tomography method of shell object, specific step is as follows for this method:
Step S101, the collected two-dimensional projection image sequence of detector is obtained, the two-dimensional projection image sequence is object
Body axially in rotation process, obtains multiple after the ray projection of the excessively described object of detector periodicity acquisition of transmission
Two-dimensional projection image, and the two-dimensional projection image of an a sampling period corresponding object.
Fig. 2 is that the quick high-resolution 3 D cone-beam computerized tomography system of plate shell object provided in an embodiment of the present invention is former
Reason figure;As shown in Fig. 2, plate and shell structure 3 D cone-beam computerized tomography system includes: x-ray source 20, X-ray beam 21, inclines
Oblique turntable 22, board under test shell object 23, detector 24, computer 25.Wherein, board under test shell object 23 is placed on inclination turntable 22
On, detector 24, inclination turntable 22 and x-ray source 20 are connect with computer 25 respectively.Computer 25 is for controlling x-ray source 20
Intensity and the time, control inclination 22 velocity of rotation of turntable of cone-beam, and the sampling period of control detector 24 occurs.X-ray
21 measuring targets 23 of cone-beam that source 20 is formed are transmitted, and 23 imaging area of object under test is covered by cone-beam 21, and turntable 22 is at the uniform velocity
Axially rotating, the periodical acquisition of transmission of detector 24 crosses the ray projection of object under test 23, after turntable 22 rotates 360 °,
Computer 25 controls the stopping of detector 24 sampling, control x-ray source 20 stops that cone-beam occurs while controlling turntable 22 to stop turning
Dynamic, entire 3 D cone-beam computerized tomography system completes primary single circular orbit cone-beam scan.Detector 24 is by each sampling
Period, collected ray projection formed a two-dimensional projection image, and two-dimensional projection image is transferred to computer 25, then turned
After platform 22 rotates 360 °, detector 24 forms a series of two-dimensional projection images i.e. two-dimensional projection image sequence, and by two-dimensional projection
Image sequence is transferred to computer 25.Fig. 3 is plate shell object novel three-dimensional filtered back projection provided in an embodiment of the present invention reconstruction
Algorithm rebuilds coordinate system schematic diagram.
Preferably, the object axially rotates, comprising: the object along rotary shaft rotation center uniform rotation
360 degree, and the image objects area is covered by cone-beam during rotation;The sampling period of the detector can flexibly set
It is fixed, it is preferred that the sampling period is the object axially 1 degree of uniform rotation of time.
In addition, the embodiment of the present invention does not limit the sampling period of detector 24.
Step S102, logarithm demodulation is carried out to the two-dimensional projection image sequence and obtains two-dimensional line integral image sequence.
The collected two-dimensional projection image sequence of the detector that computer 25 obtains it carries out logarithm demodulation and obtains two dimension
Line integral image sequence specially carries out logarithm demodulation to the two-dimensional projection image sequence I (i, m, n) according to formula (1) and obtains
Obtain two-dimensional line integral image sequence p (i, m, n):
Wherein, i indicates scanning angle, is that the multiple object axially at the uniform velocity turns simultaneously when the sampling period of detector
When dynamic 1 degree of the time, the value of i is from 1 to 360, and (m, n) indicates two of the detection channels in the vertical plane of the detector
Coordinate, ln indicate natural logrithm operation, and mean indicates that two-dimentional mean operation, 1:10 indicate that getting 10,1:K from 1 indicates to get from 1
K, I (i, 1:10,1:K) indicate that width is 10, is highly K in i-th of two-dimensional projection image in the two-dimensional projection image sequence
Region.
Step S103, using with azimuthal novel filter back-projection reconstruction algorithm to the corresponding two dimension of the plate shell object
Line integral image sequence carries out image reconstruction and obtains the corresponding three dimensional CT sectioning image of the plate shell object.Specially according to formula
(2)-(6) two dimension line integral subimage sequence p corresponding to the plate shell object carries out image reconstruction and obtains the plate shell object
Corresponding three dimensional CT sectioning image f:
X=(x " cos α+(and y " sin α+z " cos α) sin γ) cos β-(y " cos α-z " sin α) sin β (4)
Wherein, f (x ", y ", z ") indicates that reconstructed three dimensional panel shell object, (x ", y ", z ") indicate that three-dimension object particle exists
Coordinate position in object coordinates system, (x, y, z) indicate three-dimension object particle in the case where some rotates angle in system coordinate system
Coordinate position, D indicates that radiographic source arrives the distance of detector, and p (s, v, β) indicates that the two-dimensional line that obtains accumulates after carrying out logarithm demodulation
Partial image sequence, (s, v) indicate that coordinate position of certain pixel in two dimensional image coordinate system, corresponding three-dimensional object particle are mapped to
Projection address on two-dimensional line integral image, β indicate rotation angle, and h (s) is one-dimensional filtering device.α,γ is mentioned in the present invention
Three important angle parameters out, the angle between three main shafts of characterization board shell object and three-dimensional cartesian rectangular coordinate system.
For the embodiment of the present invention compared to existing Cone-Beam CT, CL technology, step is simple, reduces and rebuilds volume to improve
Reconstruction efficiency;Object space position is corrected during reconstruction, rather than image rotation is carried out to data again after rebuilding
And space interpolation, without interpolation error, thus improve reconstructed image quality;Reconstruct the sectioning image and plate shell object come
Show that legibility is strong in parallel, three-dimensional visualization precision is high.
In order to prove that the effect of above-described embodiment, the embodiment of the present invention have carried out following experiment, experimental procedure is as follows:
(1) experiment condition is set.The plate shell object that this experiment uses is two layers PCB printed circuit board.The printing
Circuit board has the characteristic feature of plate shell object, and has many characteristics, such as many cavities, solder joint, route on each layer of the circuit board.X
Radiographic source specifically uses X-ray tube, and the tube voltage of the X-ray tube is 110KV, tube current 1.5mA, x-ray source to detector
Distance be 1180mm, the angular samples number of detector is 360, and the size of two-dimensional projection image is 960*768.
(2) according to setup parameter, computer control inclination turntable rotates, detector acquired projections data, and according to projection
Data generate two-dimensional projection image sequence.
(3) computer carries out logarithm demodulation to two-dimensional projection image sequence according to above-mentioned formula (1) and obtains two-dimentional line integral
Image sequence.
(4) computer carries out the corresponding two-dimensional line integral image sequence of the plate shell object according to above-mentioned formula (2)-(6)
Image reconstruction obtains the corresponding three dimensional CT sectioning image of the object.
Fig. 4 is that the actual object obtained using the quick high-resolution 3 D cone-beam computer tomography method of plate shell object is existed
Two-dimensional line integral image under some scanning angles;Fig. 5 is to be obtained using plate shell object novel three-dimensional filter back-projection reconstruction algorithm
Three-dimensional CT image after the actual object reconstruction obtained.By Figure 4 and 5 it is found that CT image correctly discloses the PCB printed circuit board
Structure shows the validity of present invention method.Meanwhile it can be observed that with present invention method reconstruct come
Circuit board tomography " smooth vertical and horizontal ", it is parallel with plate shell object surface, it is very intuitively, readable strong.
Fig. 6 is the three-dimensional CT image that the actual object obtained is rebuild using traditional reconstruction algorithm.It can significantly be seen by Fig. 6
It observes, since traditional algorithm does not account for existing folder between plate shell object and each main shaft of three-dimensional cartesian rectangular coordinate system
Angle, therefore the object tomography reconstructed is not parallel with circuit board surface, there is obvious angle, prevent reconstruction image from it is intuitive, have
The information of the validity response circuit board.
By Fig. 4,5 and 6 it is found that energy quick high accuracy of the embodiment of the present invention realize the three-dimensional tomography of plate shell object,
The subject image for directly reconstructing out " smooth vertical and horizontal " does not need the process of subsequent image rotation interpolation, rebuilds high-efficient, reconstruction essence
Degree is high, and three-dimensional visualization effect is good.Algorithm for reconstructing realizes that process is simply easily realized.
Fig. 7 is the quick high-resolution 3 D cone-beam computer laminated imaging device of plate shell object provided in an embodiment of the present invention
Structure chart.The quick high-resolution 3 D cone-beam computer laminated imaging device of plate shell object provided in an embodiment of the present invention can execute
The process flow that the 3 D cone-beam computer tomography method embodiment of plate shell object provides, as shown in fig. 7, being directed to plate shell object
Body 3 D cone-beam computer laminated imaging device 70 includes obtaining module 71 and computing module 72, wherein obtains module 71 and is used for
The collected two-dimensional projection image sequence of detector is obtained, the two-dimensional projection image sequence is that the plate shell object turns along inclination
In platform axial direction rotation process, obtained after the ray projection of the excessively described plate shell object of detector periodicity acquisition of transmission multiple
Two-dimensional projection image, and the two-dimensional projection image of an a sampling period corresponding plate shell object;Computing module 72 is used for
Logarithm demodulation is carried out to the two-dimensional projection image sequence and obtains two-dimensional line integral image sequence, is proposed using the embodiment of the present invention
With azimuthal novel filter back-projection reconstruction algorithm to the corresponding two-dimensional line integral image sequence of the plate shell object carry out
Image reconstruction obtains the corresponding three dimensional CT sectioning image of the plate shell object.
For the embodiment of the present invention compared to existing Cone-Beam CT, CL technology, step is simple, reduces and rebuilds volume to improve
Reconstruction efficiency;Object space position is corrected during reconstruction, rather than image rotation is carried out to data again after rebuilding
And space interpolation, without interpolation error, thus improve reconstructed image quality;Reconstruct the sectioning image and plate shell object come
Show that legibility is strong in parallel, three-dimensional visualization precision is high.
On the basis of the above embodiments, the plate shell object is along inclination 360 degree of rotation center uniform rotation of turntable,
And the image objects area is covered by cone-beam during rotation;The sampling period of the detector can flexibly set, preferably
, the sampling period is the plate shell object axially 1 degree of uniform rotation of time.
Computing module 72 is specifically used for carrying out logarithm solution to the two-dimensional projection image sequence I (i, m, n) according to formula (1)
It adjusts and obtains two-dimensional line integral image sequence p (i, m, n):
Wherein, i indicates scanning angle, and (m, n) indicates two coordinates of the detection channels in the vertical plane of the detector,
Ln indicates natural logrithm operation, and mean indicates that two-dimentional mean operation, 1:10 indicate that getting 10,1:K from 1 indicates to get K, I from 1
(i, 1:10,1:K) indicates that width is 10, is highly the area of K in i-th of two-dimensional projection image in the two-dimensional projection image sequence
Domain.
Computing module 72 is specifically used for according to above-mentioned formula (2)-(6) to the corresponding two-dimensional line integral image of the plate shell object
Sequence carries out image reconstruction and obtains the corresponding three dimensional CT sectioning image of the object.
X=(x " cos α+(and y " sin α+z " cos α) sin γ) cos β-(y " cos α-z " sin α) sin β (4)
Wherein, f (x ", y ", z ") indicates that reconstructed three dimensional panel shell object, (x ", y ", z ") indicate that three-dimension object particle exists
Coordinate position in object coordinates system, (x, y, z) indicate three-dimension object particle in the case where some rotates angle in system coordinate system
Coordinate position, D indicates that radiographic source arrives the distance of detector, and p (s, v, β) indicates that the two-dimensional line that obtains accumulates after carrying out logarithm demodulation
Partial image sequence, (s, v) indicate that coordinate position of certain pixel in two dimensional image coordinate system, corresponding three-dimensional object particle are mapped to
Projection address on two-dimensional line integral image, β indicate rotation angle, and h (s) is one-dimensional filtering device.α,γ is mentioned in the present invention
Three important angle parameters out, the angle between three main shafts of characterization board shell object and three-dimensional cartesian rectangular coordinate system.
The quick high-resolution 3 D cone-beam computer laminated imaging device of plate shell object provided in an embodiment of the present invention can have
Body is for executing embodiment of the method provided by above-mentioned Fig. 1, and details are not described herein again for concrete function.
The three-dimensional tomography for realizing to energy quick high accuracy of the embodiment of the present invention plate shell object directly reconstructs out " horizontal flat
Subject image vertically " does not need the process of subsequent image rotation interpolation, rebuilds efficiency and precision is high, three-dimensional visualization effect
Good, legibility is strong, realizes that process is simple.
In conclusion the embodiment of the present invention, compared to existing imaging technique, reconstruction image volume substantially reduces, save
Conventional method is used to rebuild the reality of nonsensical volume of air, improves reconstruction efficiency;Subsequent image rotation is not needed to entangle
Positive object space, does not generate interpolation error, to improve CT picture quality, simplifies three-dimensional visualization process.
In several embodiments provided by the present invention, it should be understood that disclosed device and method can pass through it
Its mode is realized.For example, the apparatus embodiments described above are merely exemplary, for example, the division of the unit, only
Only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or components can be tied
Another system is closed or is desirably integrated into, or some features can be ignored or not executed.Another point, it is shown or discussed
Mutual coupling, direct-coupling or communication connection can be through some interfaces, the INDIRECT COUPLING or logical of device or unit
Letter connection can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list
Member both can take the form of hardware realization, can also realize in the form of hardware adds SFU software functional unit.
The above-mentioned integrated unit being realized in the form of SFU software functional unit can store and computer-readable deposit at one
In storage media.Above-mentioned SFU software functional unit is stored in a storage medium, including some instructions are used so that a computer
It is each that equipment (can be personal computer, server or the network equipment etc.) or processor (processor) execute the present invention
The part steps of embodiment the method.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory (Read-
Only Memory, ROM), random access memory (Random Access Memory, RAM), magnetic or disk etc. it is various
It can store the medium of program code.
Those skilled in the art can be understood that, for convenience and simplicity of description, only with above-mentioned each functional module
Division progress for example, in practical application, can according to need and above-mentioned function distribution is complete by different functional modules
At the internal structure of device being divided into different functional modules, to complete all or part of the functions described above.On
The specific work process for stating the device of description, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations;To the greatest extent
Pipe present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its according to
So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into
Row equivalent replacement;And these are modified or replaceed, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution
The range of scheme.
Claims (4)
1. a kind of quick high-resolution 3 D cone-beam computer tomography method of plate shell object, which is characterized in that including walking as follows
It is rapid:
Step 1 obtains the collected two-dimensional projection image sequence of detector, and the two-dimensional projection image sequence is object along axis
Into rotation process, multiple two-dimensional projections for being obtained after the ray projection of the excessively described object of detector periodicity acquisition of transmission
Image, and the two-dimensional projection image of an a sampling period corresponding object;
Step 2 carries out logarithm demodulation acquisition two-dimensional line integral image sequence to the two-dimensional projection image sequence;Utilize band orientation
The novel three-dimensional filter back-projection reconstruction algorithm at angle carries out image reconstruction to two-dimensional line integral image sequence, obtains the object pair
The three dimensional CT sectioning image answered;
It is described that logarithm demodulation acquisition two-dimensional line integral image sequence is carried out to the two-dimensional projection image sequence, comprising:
Logarithm demodulation is carried out to the two-dimensional projection image sequence I (i, m, n) according to formula (1) and obtains two-dimensional line integral image sequence
It arranges p (i, m, n):
Wherein, i indicates scanning angle, and (m, n) indicates the two-dimensional coordinate of the detection channels in the vertical plane of the detector, ln table
Showing natural logrithm operation, mean indicates two-dimentional mean operation, and 1:10 indicates that getting 10,1:K from 1 indicates to get K from 1, I (i, 1:
10,1:K) indicate that width is 10, is highly the region of K in i-th of two-dimensional projection image in the two-dimensional projection image sequence;
It is described to use with azimuthal novel filter back-projection reconstruction algorithm to the corresponding two-dimensional line integrogram of the plate shell object
The corresponding three dimensional CT sectioning image of the plate shell object is obtained as sequence carries out image reconstruction, comprising:
The corresponding two-dimensional line integral image sequence of the plate shell object is carried out described in image reconstruction acquisition according to formula (2)-(6)
The three dimensional CT sectioning image f of plate shell object:
X=(x " cos α+(and y " sin α+z " cos α) sin γ) cos β-(y " cos α-z " sin α) sin β (4)
Wherein, f (x ", y ", z ") indicates that reconstructed three dimensional panel shell object, (x ", y ", z ") indicate three-dimension object particle in object
Coordinate position in coordinate system, (x, y, z) indicate seat of the three-dimension object particle in the case where some rotates angle in system coordinate system
Cursor position, the distance of D expression radiographic source to detector, p (s, v, β) indicate the two-dimensional line integrogram obtain after logarithm demodulation
As sequence, (s, v) indicates that coordinate position of certain pixel in two dimensional image coordinate system, corresponding three-dimensional object particle are mapped to two dimension
Projection address on line integral image, β indicate rotation angle, and h (s) is one-dimensional filtering device, α,γ distinguishes characterization board shell object
With the angle between three main shafts of three-dimensional cartesian rectangular coordinate system.
2. the quick high-resolution 3 D cone-beam computer tomography method of plate shell object according to claim 1, feature
Be, step 1 further include: plate shell object be placed in inclination turntable on, during the scanning process plate shell object with turntable rotate and
Rotation;
The plate shell object along 360 degree of rotation center uniform rotation of rotary shaft, and during rotation imaging area by cone-beam
Covering.
3. a kind of quick high-resolution 3 D cone-beam computer laminated imaging device of plate shell object characterized by comprising
Module is obtained, for obtaining the collected two-dimensional projection image sequence of detector, the two-dimensional projection image sequence is object
Body axially in rotation process, obtains multiple after the ray projection of the excessively described object of detector periodicity acquisition of transmission
Two-dimensional projection image, and the two-dimensional projection image of an a sampling period corresponding object;
Computing module obtains two-dimensional line integral image sequence for carrying out logarithm demodulation to the two-dimensional projection image sequence;Benefit
Image reconstruction is carried out to two-dimensional line integral image sequence with azimuthal novel three-dimensional filter back-projection reconstruction algorithm, obtains institute
State the corresponding three dimensional CT sectioning image of object;
The computing module is specifically used for carrying out logarithm demodulation to the two-dimensional projection image sequence I (i, m, n) according to formula (1)
It obtains two-dimensional line integral image sequence p (i, m, n):
Wherein, i indicates scanning angle, and (m, n) indicates the two-dimensional coordinate of the detection channels in the vertical plane of the detector, ln table
Showing natural logrithm operation, mean indicates two-dimentional mean operation, and 1:10 indicates that getting 10,1:K from 1 indicates to get K from 1, I (i, 1:
10,1:K) indicate that width is 10, is highly the region of K in i-th of two-dimensional projection image in the two-dimensional projection image sequence;
The computing module is specifically used for carrying out image according to two-dimensional line integral image sequence of formula (2)-(6) to the object
Rebuild the three dimensional CT sectioning image f for obtaining the object:
X=(x " cos α+(and y " sin α+z " cos α) sin γ) cos β-(y " cos α-z " sin α) sin β (4)
Wherein, f (x ", y ", z ") indicates that reconstructed three dimensional panel shell object, (x ", y ", z ") indicate three-dimension object particle in object
Coordinate position in coordinate system, (x, y, z) indicate seat of the three-dimension object particle in the case where some rotates angle in system coordinate system
Cursor position, the distance of D expression radiographic source to detector, p (s, v, β) indicate the two-dimensional line integrogram obtain after logarithm demodulation
As sequence, (s, v) indicates that coordinate position of certain pixel in two dimensional image coordinate system, corresponding three-dimensional object particle are mapped to two dimension
Projection address on line integral image, β indicate rotation angle, and h (s) is one-dimensional filtering device, α,γ distinguishes characterization board shell object
With the angle between three main shafts of three-dimensional cartesian rectangular coordinate system.
4. the quick high-resolution 3 D cone-beam computer laminated imaging device of plate shell object according to claim 3, feature
It is, obtains module further include: plate shell object is placed on inclination turntable, and plate shell object is rotated with turntable during the scanning process
And it rotates;
The plate shell object along 360 degree of rotation center uniform rotation of rotary shaft, and during rotation imaging area by cone-beam
Covering.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710312113.9A CN107192726B (en) | 2017-05-05 | 2017-05-05 | The quick high-resolution 3 D cone-beam computer tomography method of plate shell object and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710312113.9A CN107192726B (en) | 2017-05-05 | 2017-05-05 | The quick high-resolution 3 D cone-beam computer tomography method of plate shell object and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107192726A CN107192726A (en) | 2017-09-22 |
CN107192726B true CN107192726B (en) | 2019-11-12 |
Family
ID=59874079
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710312113.9A Active CN107192726B (en) | 2017-05-05 | 2017-05-05 | The quick high-resolution 3 D cone-beam computer tomography method of plate shell object and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107192726B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107796835B (en) * | 2017-10-20 | 2021-05-25 | 北京航空航天大学 | X-ray cylindrical surface three-dimensional cone beam computer tomography method and device |
CN108680589B (en) * | 2018-05-31 | 2021-07-16 | 北京航空航天大学 | Three-dimensional cone beam computer tomography method and device based on transverse dislocation grating |
CN110660123B (en) * | 2018-06-29 | 2021-08-31 | 清华大学 | Three-dimensional CT image reconstruction method and device based on neural network and storage medium |
CN109325428B (en) * | 2018-09-05 | 2020-11-27 | 周军 | Human activity posture recognition method based on multilayer end-to-end neural network |
CN111462266B (en) * | 2020-03-20 | 2023-08-04 | 北京东软医疗设备有限公司 | Image reconstruction method, device, CT equipment and CT system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101718719A (en) * | 2009-12-01 | 2010-06-02 | 北京航空航天大学 | Method for determining angle increment of 3D cone-beam industry CT through continuous scan |
CN102004111A (en) * | 2010-09-28 | 2011-04-06 | 北京航空航天大学 | CT imaging method using tilted multi-cone-beam linear track |
CN103489206A (en) * | 2013-09-25 | 2014-01-01 | 华南理工大学 | Hybrid filtering method based on fan-beam X-ray CT filtering back-projection resetting |
CN103698401A (en) * | 2013-12-23 | 2014-04-02 | 电子科技大学 | Tomography method for lamb waves within any convex region |
CN105717145A (en) * | 2016-02-03 | 2016-06-29 | 北京航空航天大学 | Multi-mounted three-dimensional cone beam computer tomography method and device |
-
2017
- 2017-05-05 CN CN201710312113.9A patent/CN107192726B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101718719A (en) * | 2009-12-01 | 2010-06-02 | 北京航空航天大学 | Method for determining angle increment of 3D cone-beam industry CT through continuous scan |
CN102004111A (en) * | 2010-09-28 | 2011-04-06 | 北京航空航天大学 | CT imaging method using tilted multi-cone-beam linear track |
CN103489206A (en) * | 2013-09-25 | 2014-01-01 | 华南理工大学 | Hybrid filtering method based on fan-beam X-ray CT filtering back-projection resetting |
CN103698401A (en) * | 2013-12-23 | 2014-04-02 | 电子科技大学 | Tomography method for lamb waves within any convex region |
CN105717145A (en) * | 2016-02-03 | 2016-06-29 | 北京航空航天大学 | Multi-mounted three-dimensional cone beam computer tomography method and device |
Non-Patent Citations (1)
Title |
---|
Large field of view computed laminography with the asymmetric rotational scanning geometry;FU Jian 等;《Science china technological sciences》;20100831;第53卷(第8期);第2261-2271页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107192726A (en) | 2017-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107192726B (en) | The quick high-resolution 3 D cone-beam computer tomography method of plate shell object and device | |
US7950849B2 (en) | Method and device for geometry analysis and calibration of volumetric imaging systems | |
CN105717145B (en) | Multi-joint dress 3 D cone-beam computer tomography method and device | |
US8724889B2 (en) | Method and apparatus for CT image reconstruction | |
US8805037B2 (en) | Method and system for reconstruction of tomographic images | |
WO2006078085A1 (en) | Method for reconstructing a local high resolution x-ray ct image and apparatus for reconstructing a local high resolution x-ray ct image | |
CN102973291B (en) | C-arm tomography imaging method using semi-accurate filtered back-projection | |
JP2007319674A (en) | Image reconstruction method for its target and device for performing its method | |
WO2009043224A1 (en) | Image reconstructing method using x-ray volume photography | |
CN107796835B (en) | X-ray cylindrical surface three-dimensional cone beam computer tomography method and device | |
CN110057847B (en) | TR (transmitter-receiver) tomography projection rearrangement method and device | |
CN102631208A (en) | Device and method used for correcting geometric error of digital tomosynthsis imaging system | |
CN110133014A (en) | A kind of chip interior defect inspection method and system | |
CN105319225B (en) | A kind of scan method for realizing plaques high-resolution large-viewing open country CL imaging | |
US8861829B2 (en) | Method and system for reconstruction of tomographic images | |
JP4129572B2 (en) | Reconstruction method of tilted three-dimensional X-ray CT image | |
CN105973917B (en) | X ray CT turntable unilateral side helical scanning single-layer rearrangement method for reconstructing twice | |
Sun et al. | A reconstruction method for cone-beam computed laminography based on projection transformation | |
JP2000316843A (en) | Method and device for three-dimensional computerized tomographic imaging | |
Schrapp et al. | Artifact reduction in non-destructive testing by means of complementary data fusion of x-ray computed tomography and ultrasonic pulse-echo testing | |
JP2003329616A (en) | Eucentric inclined three-dimensional x-ray ct and imaging method of three-dimensional image by using the same | |
Guo et al. | C-arm CT imaging with the extended line-ellipse-line trajectory: first implementation on a state-of-the-art robotic angiography system | |
Meng et al. | Cone beam CT imaging with limited angle of projections and prior knowledge for volumetric verification of non-coplanar beam radiation therapy: a proof of concept study | |
Kong et al. | Analytic reconstruction approach for parallel translational computed tomography | |
Kang et al. | A projection method for reconstructing X-ray images of arbitrary cross-section |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |