CN106373165A - Tomography synthesis image reconstruction method and system - Google Patents
Tomography synthesis image reconstruction method and system Download PDFInfo
- Publication number
- CN106373165A CN106373165A CN201610796748.6A CN201610796748A CN106373165A CN 106373165 A CN106373165 A CN 106373165A CN 201610796748 A CN201610796748 A CN 201610796748A CN 106373165 A CN106373165 A CN 106373165A
- Authority
- CN
- China
- Prior art keywords
- voxel
- view
- value
- projection
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- 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
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/436—Limited angle
Abstract
The invention relates to a tomography synthesis image reconstruction method and system. The method comprises steps that S10, a non-updated region of an (n-1)th three-dimensional image is acquired, a voxel difference value of an mth reference image and the (n-1)th three-dimensional image is calculated, interpolation fitting is carried out at the non-updated region, and a compensation basic value of each voxel of the non-updated region is acquired; S20, an (n-1)th compensation image is acquired; S30, the projection data in an nth angle is reconstructed to acquire an nth three-dimensional image; S40, whether the n is equal to n0 is determined, if no, the n is updated into n+1, the process returns to the step S10, if yes, the process enters a step S50; S50, the nth three-dimensional image is taken as an mth reconstruction image result, whether the m is equal to the preset iteration times is determined, if no, the process enters a step S60, if yes, the mth reconstruction image result is taken as a final tomography synthesis image; and S60, weight processing on the mth reference image is carried out according to the mth reconstruction image result and the iteration times m, an (m+1) reference image is acquired, the m is updated into m+1, and the process returns to the step S10.
Description
Technical field
The present invention relates to tomography composograph processing technology field, more particularly to a kind of tomography composograph method for reconstructing
And system.
Background technology
Mammary gland x-ray camera chain is the important means of diagnosis of breast disease and examination, and it utilizes plate for forcing and gripper shoe
The measured targets such as breast are pressed, and is passed through the irradiated tissue to above-mentioned measured target for the soft x ray, by detector to wearing
Saturating x-ray is received and is processed, the two-dimensional projection image of measured target structure after being oppressed.However, above-mentioned mammary gland x penetrates
In line camera chain, the three dimensional structure of measured target cannot avoid tissue overlapping and mutual coverage on the two-dimensional projection image, and
And still cannot obtain the structural information of the longitudinal depth to measured targets such as breast.
For realizing the problems such as obtaining and solve tissue overlap therein of above-mentioned measured target three-dimensional information, grind at present
Send the mammary gland digital tomosynthesis system (digital based on digital tomography synthetic technology (digital tomosynthesis, dts)
Breast tomosynthesis, dbt), such system, on the basis of conventional mammary gland x-ray camera chain, sets up x ray tube
Independent rotational movement structure, keeps under rigid condition in compressor, measured target (breast etc.), detector, by bulb along circle
Arc orbital rotational motion, carries out multiple low dose exposure to the tissue in the measured target of pressing, so in limited angular range
Read several detector cells on detector afterwards and obtain a series of two-dimensional projection data, finally utilize image rebuilding method to generate
The measured target 3-D view of pressing.The 3-D view that dbt obtains can solve two-dimensional digital mammography image (digital
Mammogram, dm) tissue is overlapping and the problems such as false positive, and can with the dm image under equal oppression state realize good
Contrast and combination, improve the diagnosis efficiency of mastopathy to a certain extent.
If directly the analytic reconstruction method such as application filtered back projection fbp and fdk carries out dbt image reconstruction, the having of lack sampling
Limit sparse angular projection is more difficult to wave filter design, and the disappearance of projection information is difficult to obtain the tomography composite diagram being accurately satisfied with
Picture.And traditional backprojection reconstruction (bp) algorithm, though basic three dimensional structure information can quickly be obtained, reconstructed results lack essence
Thin details is it is impossible to directly apply to clinic.
Parallel algebraic reconstruction technique (simultaneous algebraic reconstruction technique,
Sart it is) that current dbt rebuilds accurate effective method relatively.However, in dbt imaging system, the x ray tube of circular motion and
The angle that detector normal is formed is bigger, and longitudinal depth information of acquisition is abundanter.But angle is crossed conference and is led to image objects
The subregion in region (field of vision, fov) exceeds the effective range of detector, thus causing projection information to block
And disappearance, form projection absent region, and the size of this projection absent region can be with the change of scanning angle and dimension of object
And change.When 3-D view voxel is updated under all angles in sart process of reconstruction with operation successively, corresponding angle
The voxel of lower projection absent region will be unable to obtain back projection's renewal, forms non-update area, and the projection of adjacent angular is blocked
Difference will lead to each non-update area adjacent edges pixel value in corresponding tomography composograph to occur repeatedly drastically being mutated, this
Kind phenomenon will ultimately result in the tomography composograph obtaining in reconstruction and a plurality of banding artifact along scanning direction lateral side regions, from
And affect quality and the effect that the tomography composograph of measured target is rebuild.
Content of the invention
Based on this it is necessary to affect the technical problem of measured target tomography composograph quality for traditional reconstruction scheme,
A kind of new tomography composograph method for reconstructing and system are provided.
A kind of tomography composograph method for reconstructing, comprises the steps:
S10, the non-update area of acquisition the (n-1)th 3-D view, calculate m reference picture and described (n-1)th 3-D view
Voxel difference, interpolation fitting is carried out to non-update area according to voxel difference, obtain described in each voxel of non-update area benefit
Repay base value;Wherein, n is more than 1 and to be less than or equal to n0Integer, n0For launching the child's hair twisted in a knot-childhood number of degrees of x-ray, m is current iteration time
Number, initial value is set as 1;
S20, according to described compensation base value the non-update area voxel value of the (n-1)th 3-D view is modified compensate, root
Obtain the (n-1)th compensation image according to the voxel value that correction-compensation obtains;
S30, using parallel algebraic reconstruction technique according to (n-1)th compensation image the data for projection of the n-th angle is rebuild,
Obtain the n-th 3-D view;
S40, judging that whether n is equal to n0, if it is not, then n is updated to n+1, return to step s10, if so, then entering step
s50;
S50, using the n-th 3-D view as m reconstruction image result, judge that whether m is equal to default iterationses, if
No, enter step s60, if so, then using m reconstruction image result as final tomography composograph;
S60, according to m reconstruction image result and iterationses m m reference picture is weighted process, obtain m+1
Reference picture, and m is updated to m+1, return to step s10.
A kind of tomography composograph reconstructing system, comprising:
Acquisition module, for obtaining the non-update area of the (n-1)th 3-D view, calculates m reference picture and described (n-1)th
The voxel difference of 3-D view, carries out interpolation fitting according to voxel difference to non-update area, obtain described in non-update area each
The compensation base value of voxel;Wherein, n is the integer more than 1 and less than or equal to n0, and n0 is the child's hair twisted in a knot-childhood number of degrees of transmitting x-ray, and m is
Current iteration number of times, initial value is set as 1;
Correcting module, for repairing to the non-update area voxel value of the (n-1)th 3-D view according to described compensation base value
Just compensate, the (n-1)th compensation image is obtained according to the voxel value that correction-compensation obtains;
Rebuild module, for compensating the data for projection to the n-th angle for the image using parallel algebraic reconstruction technique according to (n-1)th
Rebuild, obtained the n-th 3-D view;
Judge module, for judging whether n is equal to n0, if it is not, then n is updated to n+1, returns execution acquisition module, if
It is then to execute update module;
Update module, for as m reconstruction image result, judging the n-th 3-D view whether m is equal to default iteration
Number of times, if it is not, execution setup module, if so, then using m reconstruction image result as final tomography composograph;
Setup module, for m reference picture is weighted process with m reconstruction image result and iterationses m, obtains
To m+1 reference picture, and m is updated to m+1, returns execution acquisition module.
Above-mentioned tomography composograph method for reconstructing and system, can obtain the non-update area of the (n-1)th 3-D view, calculate
The voxel difference of m reference picture and the (n-1)th 3-D view, carries out interpolation fitting according to voxel difference to non-update area, obtains
To the compensation base value of each voxel of described non-update area, according to above-mentioned compensation base value, the voxel value of the (n-1)th 3-D view is carried out
Correction-compensation, is obtained the (n-1)th compensation image according to the voxel value that correction-compensation obtains, is carried out with the data for projection to the n-th angle
Rebuild, obtain the n-th 3-D view, and circulate said process, until obtaining the n-th 0 3-D views, using the n-th 0 3-D views as the
M reconstruction image result, and acquisition m+1 reference picture is processed with the weighting of this result, continue the data for projection of all angles is entered
Row iteration is rebuild after reaching preset times up to iterationses, using m reconstruction image result as final tomography composograph;
It can eliminate tomography composograph along a plurality of banding artifact of scanning direction lateral side regions, improves measured target tomography composite diagram
The reconstruction effect of picture.
Brief description
Fig. 1 is the tomography composograph method for reconstructing flow chart of an embodiment;
Fig. 2 is the dbt system schematic of an embodiment;
Fig. 3 is the x-ray directive measured target schematic diagram of an embodiment;
Fig. 4 is the digital breast model schematic diagram of an embodiment;
Fig. 5 is the first reference picture schematic diagram of an embodiment;
The dbt tomography composograph schematic diagram that Fig. 6 is rebuild by the present invention of an embodiment;
The dbt tomography composograph schematic diagram that Fig. 7 is rebuild by the traditional scheme of an embodiment;
Fig. 8 is the tomography composograph reconstructing system structural representation of an embodiment.
Specific embodiment
The specific embodiment of the tomography composograph method for reconstructing to the present invention and system is made in detail below in conjunction with the accompanying drawings
Description.
With reference to Fig. 1, Fig. 1 show the tomography composograph method for reconstructing flow chart of an embodiment, comprises the steps:
S10, the non-update area of acquisition the (n-1)th 3-D view, calculate m reference picture and described (n-1)th 3-D view
Voxel difference, interpolation fitting is carried out to non-update area according to voxel difference, obtain described in each voxel of non-update area benefit
Repay base value;Wherein, n be more than 1 and less than or equal to n0 integer, n0 be transmitting x-ray the child's hair twisted in a knot-childhood number of degrees, that is, dbt system sweep
The total projection number retouched;The initial value of n is 2, m is current iteration number of times, and initial value is set as 1;
Dbt system (mammary gland digital tomosynthesis system) can be as shown in Fig. 2 wherein plate for forcing, measured target and detector be protected
Hold and be relatively fixed, x-ray bulb launches x-ray in n0 angle set in advance to measured target successively along arc track,
To be scanned to corresponding measured target, gather several low dose exposure images and carry out tomography composograph reconstruction.As Fig. 2
Shown, when some angles (as larger angle angle) expose, measured target can be projected to outside detector for x-ray and detector
Portion, leads to the projection information in this region of measured target to block and lack, and forms corresponding projection absent region, causes corresponding reconstructed
The loss of learning of image.Using parallel algebraic reconstruction technique to carrying out the projection under the n-th angle and update revising the n-th three obtaining
Dimension image, there is corresponding absent region in this image.In the makeover process of data for projection under this angle of above-mentioned absent region
Less than the renewal correction of data for projection, form corresponding non-update area, because the n-th 3-D view is according to the (n-1)th compensation figure
Data for projection under picture and the n-th angle is updated revising gained, the non-update area of the therefore n-th 3-D view and the n-th angle
Under data for projection block that absent region is corresponding, that is, the non-update area of the n-th 3-D view is the projection under the n-th angle
Data block absent region.If the data for projection that detector is received directly carries out sart reconstruction, successively under all angles
When 3-D view voxel is updated with operation, the voxel of corresponding absent region will be unable to be updated, and the lacking of adjacent angular
Losing area differentiation will lead to absent region edge pixel values to occur repeatedly drastically being mutated, and cause in reconstruction image along scanning direction side
A plurality of banding artifact in face region.Above-mentioned measured target can include the tested organ of dbt system scanning.
Fig. 3 show the x-ray directive measured target schematic diagram under the (n-1)th angle, if as shown in figure 3, x-ray bulb from
, to the second end motion, the x-ray 101 obtaining near detector second end is flat with detector 103 place for the first end of detector
Face angulation 104, (first end including measured target 102 is away from spy for the position relationship of measured target 102 and detector 103
Survey the distance of device 103 first end, including the distance away from detector 103 second end for second end of measured target 102), can also obtain
Take the size (as width, height etc.) of measured target 102, relation is projected according to light and related objective geometric parameter just can calculate
The spatial information blocking absent region under (n-1)th angle, the absent region 105 including in figure (does not project to the area of detector
Domain) layer scope in measured target for the corresponding voxel, line range and row scope, thus can determine that in the (n-1)th angle reconstruction institute
The non-update area of the (n-1)th 3-D view obtaining.
Above-mentioned steps can obtain the calculation of the basic three dimensional structure of measured target using back projection (bp) method for reconstructing etc.
Method is rebuild to data for projection, to obtain the initial reference image of measured target, i.e. the first reference picture.Using back projection
(bp) method for reconstructing can be calculated using the global information of data for projection well, obtains preferably overall three-dimensional information and not
Banding artifact occurs, using above-mentioned back projection (bp) method for reconstructing, the data for projection under all angles is rebuild, permissible
Obtain the basic three dimensional structure of scanned object (measured target), mended with the correction carrying out non-update area when sart rebuilds
Repay.
S20, according to described compensation base value the non-update area voxel value of the (n-1)th 3-D view is modified compensate, root
Obtain the (n-1)th compensation image according to the voxel value that correction-compensation obtains;
Can be superimposed to compensating base value on the voxel value of the non-update area of the (n-1)th 3-D view in above-mentioned steps, with reality
The correction-compensation of existing (n-1)th 3-D view non-update area voxel.
S30, using parallel algebraic reconstruction technique according to (n-1)th compensation image the data for projection of the n-th angle is rebuild,
Obtain the n-th 3-D view;
In above-mentioned steps, angularly n carries out forward projection first can to compensate image to (n-1)th, obtains the n-th estimated projection number
According to, then the data for projection under above-mentioned n-th estimated projection data and the n-th angle is carried out weighted differences etc. process, with to front to throwing
After the result that shadow and weighted differences obtain carries out back projection, the (n-1)th compensation each voxel of image is weighted revising, is updated
The n-th 3-D view.When data for projection under above-mentioned n-th angle includes the x-ray irradiation measured target of n-th angle, detect
The data for projection that device receives.
S40, judging that whether n is equal to n0, if it is not, then n is updated to n+1, return to step s10, if so, then entering step
s50;
S50, using the n-th 3-D view as m reconstruction image result, judge that whether m is equal to default iterationses, if
No, enter step s60, if so, then using m reconstruction image result as final tomography composograph;
S60, according to m reconstruction image result and iterationses m m reference picture is weighted process, obtain m+1
Reference picture, and m is updated to m+1, return to step s10.
In the present embodiment, the projection angle that can initiate exposure from x-ray terminates the projection angle of exposure to x-ray, according to
Secondary data for projection corresponding to measured target under all angles carries out the reconstruction operation as described in step s10 to s30, when not up to
On stopping criterion for iteration (generally default iterationses, also can be using as other end conditions such as image gradient norms) repeats
State process of reconstruction, so that the tomography composograph finally giving, keep higher with the three-dimensional image information of corresponding measured target
Concordance.Above-mentioned iterationses can require be configured according to the reconstruction precision of tomography composograph, be such as set to 5 or
10 is equivalent.
In above-mentioned iterative process, the data for projection for first angle in first time iterative process, can be directly right
The data for projection of above-mentioned first angle is rebuild, and to obtain the first 3-D view, projects relation and scan geometry using light
The non-update area of parameter determination first 3-D view, calculates further according to above-mentioned first 3-D view and accordingly non-update area
Obtain the first compensation image, and carry out the weight of the second 3-D view according to the data for projection under this compensation image and second angle
Build.If in subsequently each iterative process (in addition to for the first time), can m reconstruction image obtained by above an iteration process
Result is as initial three-dimensional image, and is weighted by iterationses with m reconstruction image result and m reference picture,
Obtain m+1 reference picture, the data for projection according to above-mentioned initial three-dimensional image successively all angles is carried out as step s10 extremely
Reconstruction described in s30;I.e. the m reconstruction image result obtained by above an iteration process is as the first 3-D view, according to
Above-mentioned first 3-D view and the m+1 reference picture obtaining carry out the reconstruction to the second 3-D view.
The tomography composograph method for reconstructing that the present embodiment proposes, can calculate and obtain not updating of the (n-1)th 3-D view
Region, calculates the voxel difference of m reference picture and the (n-1)th 3-D view, according to voxel difference, non-update area is inserted
Value matching, obtain described in each voxel of non-update area compensation base value, according to above-mentioned compensation base value to the (n-1)th 3-D view not more
The each voxel value of new region is modified compensating, and obtains the (n-1)th compensation image according to the voxel value that correction-compensation obtains, according to n-th
The data for projection of angle is rebuild, and obtains the n-th 3-D view, and circulates said process, until obtaining the n-th 0 3-D views, will
The n-th 0 3-D views are as m reconstruction image result, and this time result weighting processes and obtains m+1 reference picture, continues to each
The data for projection of individual angle is iterated rebuilding after reaching default iterationses up to iteration, and m reconstruction image result is made
For final tomography composograph;It can eliminate tomography composograph along a plurality of banding artifact of scanning direction lateral side regions,
Improve the picture quality of measured target tomography composograph and rebuild effect.
In one embodiment, can also include before above-mentioned steps s10:
Multiple data for projection that detector is received carry out backprojection reconstruction, and the 3-D view that backprojection reconstruction is obtained sets
It is set to the first reference picture;
Using parallel algebraic reconstruction technique, the data for projection of first angle is rebuild, obtain the first 3-D view;
Obtain the non-update area of the first 3-D view, calculate the voxel of the first reference picture and described first 3-D view
Difference, carries out interpolation fitting according to voxel difference to non-update area, obtain described in each voxel of non-update area compensation base value;
Described compensation base value and the first reference picture are modified compensation to the voxel value of the first 3-D view, according to repairing
Just compensate the voxel value obtaining and obtain the first compensation image;
According to the first compensation image, the data for projection of second angle is rebuild using parallel algebraic reconstruction technique, obtain
Second 3-D view.
The non-update area of above-mentioned first 3-D view is the region that under first angle, detector does not receive data for projection,
I.e. under this angle, absent region is blocked in projection.The voxel difference of the first reference picture and described first 3-D view is above-mentioned the
The voxel value of each voxel difference to the voxel value of corresponding voxel in the first 3-D view respectively in one reference picture.
The present embodiment carries out regional compensation, and the projection to second angle according to the first 3-D view and the first reference picture
Carry out the reconstruction of the second 3-D view, the accuracy of the second 3-D view of reconstruction can be improved, and eliminate under first angle not
The voxel mutation of update area.
In one embodiment, above-mentioned multiple data for projection to detector reception carry out backprojection reconstruction, by back projection
The process that the 3-D view that reconstruction obtains is set to the first reference picture may include that
To detector receive multiple different angles data for projection according to detector position, ray source position, object chi
The three-dimensional relationship such as very little and light projection principle, carry out backprojection reconstruction to each voxel of target three-dimensional image, by back projection
Rebuild the 3-D view obtaining and be set to the first reference picture.
When above-mentioned dbt is to described measured target transmitting x-ray, detector can receive x-ray penetrator under different angles
Multiple data for projection of body.
As an embodiment, under the above-mentioned different angles to the detector described x-ray of reception, data for projection is counter is thrown
Shadow is rebuild, and the process that the 3-D view that backprojection reconstruction is obtained is set to the first reference picture may include that
The data for projection receiving described x-ray to detector carries out back projection (bp) method for reconstructing and processes, according in difference
The spatial informations such as the x-ray source under angle, scanning object (measured target), detector, determine that each voxel x under this angle penetrates
Line penetrates path, and obtains the pixel coordinate of the projected image under respective path, is sequentially overlapped the angled respective pixel coordinate of institute
The pixel value that is located after weighting process, you can obtain the basic three dimensional structure of measured target, described basic three dimensional structure is set
It is set to the first reference picture.
Above-mentioned back projection (bp) method for reconstructing can be calculated using the global information of data for projection well, obtains preferably
Overall three-dimensional information, be not in banding artifact, projection detector being received using above-mentioned back projection (bp) method for reconstructing
Data is rebuild, and can obtain the basic three dimensional structure of scanned object (measured target), as the first reference picture, then
Iterationses constantly weight renewal, to carry out the correction-compensation of the non-update area of sart reconstruction.
In one embodiment, above-mentioned steps s30 may include that
To the (n-1)th compensation image, angularly n carries out forward projection, obtains the n-th estimated projection data;
Carry out the weighted differences of pixel according to the data for projection under the n-th estimated projection data and the n-th angle;
After the result back projection that above-mentioned forward projection and weighted differences are obtained is processed, to the (n-1)th compensation each voxel of image
It is weighted revising, obtain the n-th 3-D view updating.
The present embodiment can angularly n carries out forward projection to the (n-1)th compensation image, and under simulation current angular n, bulb is sent out
Penetrate x-ray and penetrate scanning project objects to detector, calculate each pixel of projected image (detector unit) corresponding x-ray
Penetrate path, and spatially each voxel is weighed by respective path will to penetrate corresponding n-1 compensation image three-dimensional on path at this successively
It is weighted adding up again, obtain the n-th estimated projection data;Entered according to the data for projection under the n-th estimated projection data and the n-th angle
Row weighted differences, the pixel value of each pixel corresponding with acquired projections for estimated projection is carried out difference process, and presses this pixel
The total length that the x-ray that (detector unit) is located penetrates scanning object is weighted to this pixel difference processing;To forward projection and
The result that weighted differences obtain carries out back projection's weighting and revises, angularly corresponding bulb, detector, the space of scanning object under n
Position, obtains the projection corresponding x-ray of each coordinate of difference result, is pierced scanning object by the acquisition of this x-ray path corresponding
3 d space coordinate, above-mentioned projection difference result is pressed the path weight value of each voxel on the path of place, weighted superposition to (n-1)th
Compensate the pixel value in image corresponding three-dimensional coordinate pixel, obtain updating the n-th 3-D view revised.
The present embodiment, under the conditions of the dbt of detector finite size and wide-angle scanning rebuilds, can be mended by overall situation projection
The image repaying correction carries out sart reconstruction, overcomes the reconstruction image scanning direction lateral side regions banding that missing projection data causes pseudo-
Shadow, it is ensured that longitudinal depth information of the tomography composograph rebuild, improves the quality of image reconstruction simultaneously.
In one embodiment, above-mentioned steps s60 may include that
The voxel value of the voxel value to described m reconstruction image result and m reference picture substitutes into weighting more new formula, meter
Calculate the voxel value of m+1 reference picture;Wherein, described weighting more new formula is:
refm+1(i, j, k)=(refm(i,j,k)+(m-1)presm(i,j,k))/(1+(m-1)p),
Wherein, refm+1(i, j, k) represents the voxel value of the i-th row jth row kth layer in m+1 reference picture, refm(i,j,
K) voxel value of the i-th row jth row kth layer in m reference picture, res are representedm(i, j, k) represents the in m reconstruction image result
The voxel value of i row jth row kth layer, m represents current iteration number of times, and p is to update weight coefficient.Above-mentioned renewal weight coefficient p is permissible
It is set greater than 0 value.
The present embodiment in an iterative process, is carried out adding by iterationses to m reference picture according to m reconstruction image result
Power updates it is ensured that while the basic structure of initial reference image, stepping up the quality of reference picture according to iterationses, entering
One step improves the accuracy of sart process of reconstruction non-update area compensating approach.In one embodiment, calculate m reference picture
Include with the process of the voxel difference of described (n-1)th 3-D view:
The voxel value of the voxel value of described m reference picture and the (n-1)th 3-D view is substituted into voxel differences formula meter respectively
Calculate the voxel difference of the (n-1)th 3-D view;Wherein, described voxel differences formula is:
diffn-1(i, j, k)=imgn-1(i,j,k)-refm(i, j, k),
Wherein, imgn-1(i, j, k) represents the voxel value of the i-th row jth row kth layer in the (n-1)th 3-D view, refm(i,j,
K) voxel value of the i-th row jth row kth layer in m reference picture, diff are representedn-1(i, j, k) represents the in the (n-1)th 3-D view
The voxel difference of i row jth row kth layer, m represents current iteration number of times.
The present embodiment calculates the voxel difference of m reference picture and the (n-1)th 3-D view, to insert to non-update area
Value, the local voxel that can eliminate m reference picture with the n-th 3-D view of homolographic projection data reconstruction is mutated the mistake brought
Difference.
As an embodiment, above-mentioned interpolation fitting is carried out to non-update area according to voxel difference, obtain described in not more
The process of the compensation base value of each voxel of new region includes:
Described voxel difference is substituted into the neighborhood correction value that neighboring mean value formula calculates non-update area;
The neighborhood correction value of described non-update area is substituted into the compensation base value that interpolation fitting formula calculates non-update area;
Described neighboring mean value formula is:
Described interpolation fitting formula is:
Wherein, ndiffn-1(i, k) represents the neighborhood correction value of the i-th row kth layer in the (n-1)th 3-D view, and j represents (n-1)th
The border columns of voxel, jb in 3-D viewn-1(i, k) represents the border of non-update area in (n-1)th 3-D view the i-th row kth
Columns residing for layer, n represents the voxel columns of default neighborhood, comn-1(i, j, k) represents the i-th row jth in the (n-1)th 3-D view
The compensation base value of row kth layer.
Above-mentioned neighborhood can carry out pre- according to the voxel columns of the voxel columns of measured target and corresponding non-update area
If it is generally the case that when the data for projection that the (n-1)th 3-D view is carried out with the n-th angle updates, corresponding neighborhood could be arranged to
In (n-1)th 3-D view, from the upper Angles Projections non-update area (correspondence that i.e. data for projection of the (n-1)th angle cannot update
Region) boundary line start, to some row voxels of non-non- update area side;If above-mentioned neighborhood is set to never update area side
Boundary line starts, to two row voxels of non-non- update area side, then n=2.
The present embodiment, according to neighborhood columns n, neighborhood correction value ndiff (i, k), can successively be pressed scanning course bearing matching and work as
The compensating curve of each voxel of front non-update area.Above-mentioned fit approach may include linear fit, it is possible to use the side such as curve matching
Formula calculates the offset of each voxel of non-update area.
In one embodiment, above-mentioned according to described compensation base value the voxel value of the (n-1)th 3-D view is modified mend
The process repaid includes:
The voxel value of described compensation base value and the (n-1)th 3-D view is substituted into correction-compensation formula respectively and enters row operation;Institute
Stating correction-compensation formula is:
Wherein, fixn-1(i, j, k) represents the voxel value of (n-1)th compensation image the i-th row jth row kth layer, jbn-1(i, k) table
Show the border of the non-update area columns residing in (n-1)th 3-D view the i-th row kth layer, imgn-1(i, j, k) represents the (n-1)th three
The voxel value of the i-th row jth row kth layer voxel, com in dimension imagen-1(i, j, k) represents the compensation base of the i-th row jth row kth layer
Value, refm(i, j, k) represents the voxel value of m reference picture the i-th row jth row kth layer, and m represents current iteration number of times.
The present embodiment using superposition revise by the way of it is ensured that superposition after obtained non-update area voxel be worth to entirely
The compensation of office's projection information, it is to avoid under traditional scheme, data for projection blocks the voxel value mutation that disappearance causes.
As an embodiment, the above-mentioned voxel by described compensation base value superposition corresponding with the (n-1)th angle 3-D view
The process that value is overlapped can also adopt other stacked systems such as mean filter stacked system.Above-mentioned mean filter superposition can
For:
Wherein, fixn-1(i, j, k) is the voxel value that superposition gained (n-1)th compensates image the i-th row jth row kth layer, comn-1
(i, j, k) represents the compensation base value of (n-1)th 3-D view the i-th row jth row kth layer, refm(i, o, k) represents m reference picture
The voxel value of the i-th row o row kth layer, the length of l mean filter, m represents current iteration number of times.
Above-mentioned local interpolation method, in conjunction with the first reference picture using neighborhood matching and local interpolation method to not updating
Region carries out Interpolation compensation, obtains information using the global information of data for projection and more completely compensates 3-D view, above-mentioned benefit
Repay 3-D view and can be used for sart reconstruction projection renewal operation next time, can solve well to lack neighborhood border pixel values
It is mutated the problem of the gray value acute variation leading to, it is to avoid the generation of banding artifact.
The tomography composograph method for reconstructing that the present invention provides, by predetermined angle range scans object to be reconstructed, difference
Obtain corresponding data for projection under each angle, backprojection reconstruction method process is carried out to the data for projection obtaining, obtain scanned
The basic three dimensional structure of object, and as initial reference image (the first reference picture);According to mechanical scanning geological information and light
Line projects relation, respectively data for projection back projection under each angle is processed, and accurately calculates dividing of corresponding non-update area
Cloth;Sart reconstruction is carried out to data for projection, combines the first reference picture and to last time angle reconstruction image not under current angular
Update area carries out local interpolation, carries out data correction to the voxel of non-update area;Carried out according to revised 3-D view
Forward projection weighted differences and back projection's weighting are revised, and complete the renewal of 3-D view under Current projection, then complete all angles
Data for projection under degree updates, and obtains when time iteration result;According to iterationses and iteration result, reference picture is updated, continue
Iterative approximation obtains final dbt tomography composograph.Its ensure detector size finite sum wide-angle scanning on the premise of,
The reconstruction image that missing projection data causes can be overcome along scanning direction lateral side regions banding artifact it is ensured that dbt image is longitudinal
Depth information, improves the quality of image reconstruction simultaneously.
With reference to a concrete application scene, the present invention above-mentioned tomography composograph method for reconstructing is illustrated.
(model matrix size is 400*600*50, and voxel is big as measured target to select digital breast model as shown in Figure 4
Little for 0.5*0.5*1mm).Dbt rebuilds (reconstruction of tomography composograph) step and may include that
1) breast model is carried out with dbt scanning, x-ray bulb range of movement, by -30 ° to+30 °, is exposed every 3 °
And projection data acquisitions, gather 21 data for projection proj altogethern, n ∈ [1,21] it is assumed that detector matrix size be 600*800,
Pixel size is 0.4*0.4mm, and the data for projection matrix that correspondence obtains is 600*800*21;
2) 21 data for projection proj simultaneously to collection for back projection (bp) method for reconstructing are utilizednCarry out backprojection reconstruction,
Obtain reference picture as shown in Figure 5, what the reconstruction of back projection (bp) method for reconstructing obtained obtains basic three dimensional structure, above-mentioned in order to make
First reference picture ref of measured target1;
3) preset initial three-dimensional image img0For 0, according to the data for projection proj of the 1st angle1Carry out the forward projection of sart
After weighted differences, back projection obtains the 1st reconstruction image img after revising to reconstruction domain weighting1;
4), from the beginning of the 2nd angle, that is, during n > 1, according to the perspective geometry setting and light projection principle, calculate (n-1)th
3-D view imgn-1Non- update area, calculate m reference picture refm-1With described (n-1)th 3-D view imgn-1Voxel
Difference, according to voxel difference to non-update area unn-1Carry out interpolation fitting, obtain described in each voxel of non-update area compensation
Base value comn-1;Wherein, n be more than 1 and less than or equal to n0 integer, n0 be transmitting x-ray the child's hair twisted in a knot-childhood number of degrees, n0=herein
21, m is current iteration number of times, and initial value is set as 1;
Described neighboring mean value formula is:
Described interpolation fitting formula is:
Wherein, ndiffn-1(i, k) represents the neighborhood correction value of (n-1)th 3-D view the i-th row kth layer, and j represents the (n-1)th three
The border columns of voxel, jb in dimension imagen-1(i, k) represents the border of non-update area in (n-1)th 3-D view the i-th row kth layer
Residing columns, n represents the voxel columns of default neighborhood, n=3 herein, comn-1(i, j, k) represents the i-th row jth row kth layer
Compensate base value, the matching compensating base value is fitted using linear mode.
5) utilize described correction-compensation formula and compensate base value comn-1To the (n-1)th 3-D view imgn-1Non- update area
Voxel value is modified compensating, and obtains the (n-1)th compensation image fix according to the voxel value that correction-compensation obtainsn-1;Described correction is mended
Repaying formula is:
Wherein, fixn-1(i, j, k) represents the voxel value of (n-1)th compensation image the i-th row jth row kth layer, jbn-1(i, k) table
Show the border of the non-update area columns residing in (n-1)th 3-D view the i-th row kth layer, imgn-1(i, j, k) represents the (n-1)th three
The voxel value of the i-th row jth row kth layer voxel, com in dimension imagen-1(i, j, k) represents that (n-1)th 3-D view the i-th row jth arranges the
The compensation base value of k layer, refm(i, j, k) represents the voxel value of m reference picture the i-th row jth row kth layer, and m represents current iteration
Number of times.
6) utilize parallel algebraic reconstruction sart technology according to the (n-1)th compensation image fixn-1Data for projection to the n-th angle
projnRebuild, that is, obtained the n-th 3-D view imgn;
To the (n-1)th compensation image fixn-1Angularly n carries out forward projection, bulb transmitting x-ray under simulation current angular n
And penetrate scanning project objects to detector, calculate each pixel of projected image (detector unit) corresponding x-ray penetrates road
Footpath, and corresponding n-1 compensation image fix on path will be penetrated at this successivelyn-1On three dimensions, each voxel presses respective path weight
It is weighted adding up, obtain the n-th estimated projection data aprojn;
According to the n-th estimated projection data aprojnWith the data for projection proj under the n-th anglenCarry out weighted differences, will estimate
The pixel value projecting each pixel corresponding with acquired projections carries out difference process, and the x being located by this pixel (detector unit)
The total length that ray penetrates scanning object is weighted to this pixel difference processing;
Result dproj that forward projection and weighted differences are obtainednCarry out back projection's weighting to revise, angularly corresponding under n
Bulb, detector, the locus of scanning object, obtain projection difference result dprojnThe corresponding x-ray of each coordinate, by this x
Ray path obtains and is pierced the scanning corresponding 3 d space coordinate of object, above-mentioned projection difference result is pressed each on the path of place
The path weight value of individual voxel, weighted superposition to the (n-1)th compensation image fixn-1Pixel value in corresponding three-dimensional coordinate pixel, obtains
Update the n-th 3-D view img revisingn.
7) judge whether n is equal to 21, if it is not, then n is updated to n+1, return to step 4) carry out under each angle not successively
Update area calculates, 3-D view compensates, data for projection is rebuild, and the renewal if so, then completing all data for projection is rebuild, and enters
Step next step 8);
8) by the n-th 3-D view imgnAs m reconstruction image result resm, judge whether m is equal to 5 (default iteration time
Number), if so, then using m reconstruction image result as final tomography composograph, if it is not, then to m reconstruction image result
resmVoxel value and m reference picture refmVoxel value substitute into weighting more new formula, calculate m+1 reference picture refm+1;
Wherein, described weighting more new formula is:
refm+1(i, j, k)=(refm(i,j,k)+(m-1)presm(i,j,k))/(1+(m-1)p),
Wherein, refm+1(i, j, k) represents the voxel value of the i-th row jth row kth layer in m+1 reference picture, refm(i,j,
K) voxel value of the i-th row jth row kth layer in m reference picture, res are representedm(i, j, k) represents the in m reconstruction image result
The voxel value of i row jth row kth layer, m represents current iteration number of times, and p is the renewal weight coefficient more than 0, and this application scenarios is arranged
For 1.
Fig. 6 show the dbt tomography composograph that application the present invention program is rebuild, and Fig. 7 show application traditional scheme
The dbt tomography composograph rebuild, as shown in fig. 7, using traditional sart method for reconstructing, due to limited by detector size and
Angle is excessive, leads to data for projection disappearance under some angles to be blocked, leads to the imaging fov of this angle back projection disappearance area
Domain, when each projection updates iterative operation, producing absent region boundary voxel mutation, forming many banding artifact, thus affecting
The effect that corresponding tomography composograph is rebuild.The tomography composograph method for reconstructing that the application present invention is provided is rebuild to be obtained
Tomography composograph (as Fig. 6) is compared with the dbt tomography composograph shown in Fig. 7, it can be shown that disconnected using the present invention
Lamination becomes image rebuilding method, the structural information basic of the reconstruction image of scanned object (measured target) and preferable die body
Cause;The tomography composograph (as Fig. 7) of many banding artifact, tomography composograph above-mentioned by the present invention occur with respect to border
It is clear that relatively sharp accurate, the above-mentioned wide-angle of effectively solving scans and detector tomography composograph obtained by method for reconstructing
The size-constrained banding artifact problem leading to, and peripheral components need not be increased and realize, execution efficiency is high, and stability is strong.
With reference to Fig. 8, Fig. 8 show the tomography composograph reconstructing system structural representation of an embodiment, comprising:
Acquisition module 10, for obtaining the non-update area of the (n-1)th 3-D view, calculates m reference picture and described the
The voxel difference of n-1 3-D view, carries out interpolation fitting according to voxel difference to non-update area, obtain described in non-update area
The compensation base value of each voxel;Wherein, n be more than 1 and less than or equal to n0 integer, n0 be transmitting x-ray the child's hair twisted in a knot-childhood number of degrees, m
For current iteration number of times, initial value is set as 1;
Correcting module 20, for carrying out to the non-update area voxel value of the (n-1)th 3-D view according to described compensation base value
Correction-compensation, obtains the (n-1)th compensation image according to the voxel value that correction-compensation obtains;
Rebuild module 30, for compensating the projection number to the n-th angle for the image using parallel algebraic reconstruction technique according to (n-1)th
According to being rebuild, obtain the n-th 3-D view;
Judge module 40, for judging whether n is equal to n0, if it is not, then n is updated to n+1, returns execution acquisition module,
If so, then execute update module;
Update module 50, for as m reconstruction image result, judging the n-th 3-D view whether m is equal to default changing
Generation number, if it is not, execution setup module, if so, then using m reconstruction image result as final tomography composograph;
Setup module 60, for m reconstruction image result and iterationses m are weighted to m reference picture processing,
Obtain m+1 reference picture, and m is updated to m+1, return execution acquisition module.
In one embodiment, above-mentioned reconstruction module is further used for:
To the (n-1)th compensation image, angularly n carries out forward projection, obtains the n-th estimated projection data;
Carry out the weighted differences of pixel according to the data for projection under the n-th estimated projection data and the n-th angle;
After the result back projection that above-mentioned forward projection and weighted differences are obtained is processed, to the (n-1)th compensation each voxel of image
It is weighted revising, obtain the n-th 3-D view updating.
In one embodiment, above-mentioned setup module is further used for:
The voxel value of the voxel value to described m reconstruction image result and m reference picture substitutes into weighting more new formula, meter
Calculate the voxel value of m+1 reference picture;Wherein, described weighting more new formula is:
refm+1(i, j, k)=(refm(i,j,k)+(m-1)presm(i,j,k))/(1+(m-1)p),
Wherein, refm+1(i, j, k) represents the voxel value of the i-th row jth row kth layer in m+1 reference picture, refm(i,j,
K) voxel value of the i-th row jth row kth layer in m reference picture, res are representedm(i, j, k) represents the in m reconstruction image result
The voxel value of i row jth row kth layer, m represents current iteration number of times.
The tomography composograph method for reconstructing one that the tomography composograph reconstructing system that the present invention provides is provided with the present invention
One is corresponding, and the technical characteristic illustrating in the embodiment of described tomography composograph method for reconstructing and its advantage are all applied to disconnected
Lamination becomes in the embodiment of image re-construction system, hereby give notice that.
Each technical characteristic of embodiment described above can arbitrarily be combined, for making description succinct, not to above-mentioned reality
The all possible combination of each technical characteristic applied in example is all described, as long as however, the combination of these technical characteristics is not deposited
In contradiction, all it is considered to be the scope of this specification record.
Embodiment described above only have expressed the several embodiments of the present invention, and its description is more concrete and detailed, but simultaneously
Can not therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art
Say, without departing from the inventive concept of the premise, some deformation can also be made and improve, these broadly fall into the protection of the present invention
Scope.Therefore, the protection domain of patent of the present invention should be defined by claims.
Claims (10)
1. a kind of tomography composograph method for reconstructing is it is characterised in that comprise the steps:
S10, the non-update area of acquisition the (n-1)th 3-D view, calculate the body of m reference picture and described (n-1)th 3-D view
Plain difference, carries out interpolation fitting according to voxel difference to non-update area, obtain described in each voxel of non-update area compensation base
Value;Wherein, n be more than 1 and less than or equal to n0 integer, n0 be transmitting x-ray the child's hair twisted in a knot-childhood number of degrees, m be current iteration number of times,
Initial value is set as 1;
S20, the non-update area voxel value of the (n-1)th 3-D view is modified compensate according to described compensation base value, according to repairing
Just compensate the voxel value obtaining and obtain the (n-1)th compensation image;
S30, using parallel algebraic reconstruction technique according to (n-1)th compensation image the data for projection of the n-th angle is rebuild, obtain
N-th 3-D view;
S40, judging that whether n is equal to n0, if it is not, then n is updated to n+1, return to step s10, if so, then entering step s50;
S50, using the n-th 3-D view as m reconstruction image result, judge that whether m is equal to default iterationses, if it is not, entering
Enter step s60, if so, then using m reconstruction image result as final tomography composograph;
S60, according to m reconstruction image result and iterationses m m reference picture is weighted process, obtain m+1 reference
Image, and m is updated to m+1, return to step s10.
2. tomography composograph method for reconstructing according to claim 1 is it is characterised in that also wrap before described step s10
Include:
Multiple data for projection that detector is received carry out backprojection reconstruction, and the 3-D view that backprojection reconstruction is obtained is set to
First reference picture;
Using parallel algebraic reconstruction technique, the data for projection of first angle is rebuild, obtain the first 3-D view;
Obtain the non-update area of the first 3-D view, calculate the voxel differences of the first reference picture and described first 3-D view
Value, carries out interpolation fitting according to voxel difference to non-update area, obtain described in each voxel of non-update area compensation base value;
Described compensation base value and the first reference picture are modified to the voxel value of the first 3-D view compensating, mend according to revising
Repay the voxel value obtaining and obtain the first compensation image;
According to the first compensation image, the data for projection of second angle is rebuild using parallel algebraic reconstruction technique, obtain second
3-D view.
3. tomography composograph method for reconstructing according to claim 1 is it is characterised in that described step s30 includes:
To the (n-1)th compensation image, angularly n carries out forward projection, obtains the n-th estimated projection data;
Carry out the weighted differences of pixel according to the data for projection under the n-th estimated projection data and the n-th angle;
After the result back projection that above-mentioned forward projection and weighted differences are obtained is processed, the (n-1)th compensation each voxel of image is carried out
Weighting is revised, and obtains the n-th 3-D view updating.
4. tomography composograph method for reconstructing according to claim 1 is it is characterised in that described step s60 includes:
The voxel value of the voxel value to described m reconstruction image result and m reference picture substitutes into weighting more new formula, calculates the
The voxel value of m+1 reference picture;Wherein, described weighting more new formula is:
refm+1(i, j, k)=(refm(i,j,k)+(m-1)presm(i,j,k))/(1+(m-1)p),
Wherein, refm+1(i, j, k) represents the voxel value of the i-th row jth row kth layer in m+1 reference picture, refm(i, j, k) table
Show the voxel value of the i-th row jth row kth layer in m reference picture, resm(i, j, k) represents the i-th row in m reconstruction image result
The voxel value of jth row kth layer, m represents current iteration number of times, and p is to update weight coefficient.
5. tomography composograph method for reconstructing according to claim 1 is it is characterised in that described calculating m reference picture
Include with the process of the voxel difference of described (n-1)th 3-D view:
The voxel value of the voxel value of described m reference picture and the (n-1)th 3-D view is substituted into voxel differences formula respectively and calculates the
The voxel difference of n-1 3-D view;Wherein, described voxel differences formula is:
diffn-1(i, j, k)=imgn-1(i,j,k)-refm(i,j,k);
Wherein, imgn-1(i, j, k) represents the voxel value of the i-th row jth row kth layer in the (n-1)th 3-D view, refm(i, j, k) table
Show the voxel value of the i-th row jth row kth layer in m reference picture, diffn-1(i, j, k) represents the i-th row in the (n-1)th 3-D view
The voxel difference of jth row kth layer, m represents current iteration number of times.
6. tomography composograph method for reconstructing according to claim 5 it is characterised in that described according to voxel difference to not
Update area carries out interpolation fitting, obtain described in the process of compensation base value of each voxel of non-update area include:
Described voxel difference is substituted into the neighborhood correction value that neighboring mean value formula calculates non-update area;
The neighborhood correction value of described non-update area is substituted into the compensation base value that interpolation fitting formula calculates non-update area;
Described neighboring mean value formula is:
Described interpolation fitting formula is:
Wherein, ndiffn-1(i, k) represents the neighborhood correction value of the i-th row kth layer in the (n-1)th 3-D view, and j represents that (n-1)th is three-dimensional
The border columns of voxel, jb in imagen-1(i, k) represents the border of non-update area in (n-1)th 3-D view the i-th row kth layer institute
The columns at place, n represents the voxel columns of default neighborhood, comn-1(i, j, k) represents (n-1)th 3-D view the i-th row jth row kth layer
Compensation base value.
7. tomography composograph method for reconstructing according to claim 1 it is characterised in that described according to compensating base value to the
The process that the voxel value of n-1 3-D view is modified compensating includes:
The voxel value of described compensation base value and the (n-1)th 3-D view is substituted into correction-compensation formula respectively and enters row operation;Described repair
Positive compensation formula is:
Wherein, fixn-1(i, j, k) represents the voxel value of (n-1)th compensation image the i-th row jth row kth layer, jbn-1(i, k) represents not
Columns residing in (n-1)th 3-D view the i-th row kth layer for the border of update area, imgn-1(i, j, k) represents the (n-1)th graphics
The voxel value of the i-th row jth row kth layer voxel, com in picturen-1(i, j, k) represents (n-1)th 3-D view the i-th row jth row kth layer
Compensation base value, refm(i, j, k) represents the voxel value of m reference picture the i-th row jth row kth layer.
8. a kind of tomography composograph reconstructing system is it is characterised in that include:
Acquisition module, for obtaining the non-update area of the (n-1)th 3-D view, calculates m reference picture three-dimensional with described (n-1)th
The voxel difference of image, carries out interpolation fitting according to voxel difference to non-update area, obtain described in each voxel of non-update area
Compensation base value;Wherein, n be more than 1 and less than or equal to n0 integer, n0 be transmitting x-ray the child's hair twisted in a knot-childhood number of degrees, m is current
Iterationses, initial value is set as 1;
Correcting module, for being modified to the non-update area voxel value of the (n-1)th 3-D view mending according to described compensation base value
Repay, the (n-1)th compensation image is obtained according to the voxel value that correction-compensation obtains;
Rebuild module, for carrying out to the data for projection of the n-th angle according to the (n-1)th compensation image using parallel algebraic reconstruction technique
Rebuild, obtain the n-th 3-D view;
Judge module, for judging whether n is equal to n0, if it is not, then n is updated to n+1, returns execution acquisition module, if so, then
Execution update module;
Update module, for as m reconstruction image result, judging the n-th 3-D view whether m is equal to default iteration time
Number, if it is not, execution setup module, if so, then using m reconstruction image result as final tomography composograph;
Setup module, for m reference picture is weighted process with m reconstruction image result and iterationses m, obtains the
M+1 reference picture, and m is updated to m+1, return execution acquisition module.
9. tomography composograph reconstructing system according to claim 8 is it is characterised in that described reconstruction module is used further
In:
To the (n-1)th compensation image, angularly n carries out forward projection, obtains the n-th estimated projection data;
Carry out the weighted differences of pixel according to the data for projection under the n-th estimated projection data and the n-th angle;
After the result back projection that above-mentioned forward projection and weighted differences are obtained is processed, the (n-1)th compensation each voxel of image is carried out
Weighting is revised, and obtains the n-th 3-D view updating.
10. tomography composograph reconstructing system according to claim 8 is it is characterised in that described setup module is further
For:
The voxel value of the voxel value to described m reconstruction image result and m reference picture substitutes into weighting more new formula, calculates the
The voxel value of m+1 reference picture;Wherein, described weighting more new formula is:
refm+1(i, j, k)=(refm(i,j,k)+(m-1)presm(i,j,k))/1+(m-1)p,
Wherein, refm+1(i, j, k) represents the voxel value of the i-th row jth row kth layer in m+1 reference picture, refm(i, j, k) table
Show the voxel value of the i-th row jth row kth layer in m reference picture, resm(i, j, k) represents the i-th row in m reconstruction image result
The voxel value of jth row kth layer, m represents current iteration number of times.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610796748.6A CN106373165B (en) | 2016-08-31 | 2016-08-31 | Tomography composograph method for reconstructing and system |
PCT/CN2016/098738 WO2018040126A1 (en) | 2016-08-31 | 2016-09-12 | Tomosynthesis image reconstruction method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610796748.6A CN106373165B (en) | 2016-08-31 | 2016-08-31 | Tomography composograph method for reconstructing and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106373165A true CN106373165A (en) | 2017-02-01 |
CN106373165B CN106373165B (en) | 2017-11-28 |
Family
ID=57899881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610796748.6A Active CN106373165B (en) | 2016-08-31 | 2016-08-31 | Tomography composograph method for reconstructing and system |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN106373165B (en) |
WO (1) | WO2018040126A1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108010096A (en) * | 2017-11-30 | 2018-05-08 | 上海联影医疗科技有限公司 | CBCT image rebuilding methods, device and CBCT equipment |
CN109658465A (en) * | 2018-12-07 | 2019-04-19 | 广州华端科技有限公司 | Data processing, image rebuilding method and device in image reconstruction process |
CN109978988A (en) * | 2017-12-21 | 2019-07-05 | 西门子保健有限责任公司 | For rebuilding method, the biplane X-ray device of three-dimensional image data sets |
CN110335204A (en) * | 2019-05-07 | 2019-10-15 | 中国人民解放军陆军工程大学 | A kind of graphic images Enhancement Method |
CN110706338A (en) * | 2019-09-30 | 2020-01-17 | 东软医疗系统股份有限公司 | Image reconstruction method and device, CT (computed tomography) equipment and CT system |
CN113168721A (en) * | 2018-11-28 | 2021-07-23 | 皇家飞利浦有限公司 | System for reconstructing an image of an object |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019047545A1 (en) * | 2018-05-04 | 2019-03-14 | 西安大医集团有限公司 | Low-dose imaging method and apparatus |
CN109147035B (en) * | 2018-08-03 | 2023-08-08 | 上海电气集团股份有限公司 | Three-dimensional model display method and system |
CN112258627B (en) * | 2020-09-18 | 2023-09-15 | 中国科学院计算技术研究所 | Local fault three-dimensional reconstruction system |
CN115147378B (en) * | 2022-07-05 | 2023-07-25 | 哈尔滨医科大学 | CT image analysis and extraction method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006038142A1 (en) * | 2004-10-08 | 2006-04-13 | Koninklijke Philips Electronics N.V. | Image reconstruction with voxel dependent interpolation |
CN102419866A (en) * | 2010-09-27 | 2012-04-18 | 北京农业智能装备技术研究中心 | Construction method of filter function for CT (computerized tomography) image tomography reconstruction and tomography reconstruction method |
CN103971349A (en) * | 2013-01-30 | 2014-08-06 | 上海西门子医疗器械有限公司 | Computer tomography image reconstruction method and computer tomography device |
US20150221124A1 (en) * | 2014-02-05 | 2015-08-06 | University Of Utah | Iterative reconstruction of image data in ct |
CN105184835A (en) * | 2015-09-15 | 2015-12-23 | 上海联影医疗科技有限公司 | Mammary gland tomographic image reconstruction method and device |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103757B (en) * | 2010-12-27 | 2012-09-19 | 中国科学院深圳先进技术研究院 | Cone beam image rebuilding method and device |
CN103455989B (en) * | 2013-09-24 | 2016-06-15 | 南京大学 | A kind of method improving limited angle CT image quality in conjunction with ultrasonoscopy |
CN105488823B (en) * | 2014-09-16 | 2019-10-18 | 株式会社日立制作所 | CT image rebuilding method, CT equipment for reconstructing image and CT system |
-
2016
- 2016-08-31 CN CN201610796748.6A patent/CN106373165B/en active Active
- 2016-09-12 WO PCT/CN2016/098738 patent/WO2018040126A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006038142A1 (en) * | 2004-10-08 | 2006-04-13 | Koninklijke Philips Electronics N.V. | Image reconstruction with voxel dependent interpolation |
CN102419866A (en) * | 2010-09-27 | 2012-04-18 | 北京农业智能装备技术研究中心 | Construction method of filter function for CT (computerized tomography) image tomography reconstruction and tomography reconstruction method |
CN103971349A (en) * | 2013-01-30 | 2014-08-06 | 上海西门子医疗器械有限公司 | Computer tomography image reconstruction method and computer tomography device |
US20150221124A1 (en) * | 2014-02-05 | 2015-08-06 | University Of Utah | Iterative reconstruction of image data in ct |
CN105184835A (en) * | 2015-09-15 | 2015-12-23 | 上海联影医疗科技有限公司 | Mammary gland tomographic image reconstruction method and device |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108010096A (en) * | 2017-11-30 | 2018-05-08 | 上海联影医疗科技有限公司 | CBCT image rebuilding methods, device and CBCT equipment |
US11734862B2 (en) | 2017-11-30 | 2023-08-22 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image reconstruction |
CN109978988A (en) * | 2017-12-21 | 2019-07-05 | 西门子保健有限责任公司 | For rebuilding method, the biplane X-ray device of three-dimensional image data sets |
CN109978988B (en) * | 2017-12-21 | 2023-03-24 | 西门子保健有限责任公司 | Method for reconstructing a three-dimensional image data set, biplane X-ray device |
CN113168721A (en) * | 2018-11-28 | 2021-07-23 | 皇家飞利浦有限公司 | System for reconstructing an image of an object |
CN109658465A (en) * | 2018-12-07 | 2019-04-19 | 广州华端科技有限公司 | Data processing, image rebuilding method and device in image reconstruction process |
CN110335204A (en) * | 2019-05-07 | 2019-10-15 | 中国人民解放军陆军工程大学 | A kind of graphic images Enhancement Method |
CN110706338A (en) * | 2019-09-30 | 2020-01-17 | 东软医疗系统股份有限公司 | Image reconstruction method and device, CT (computed tomography) equipment and CT system |
CN110706338B (en) * | 2019-09-30 | 2023-05-02 | 东软医疗系统股份有限公司 | Image reconstruction method, device, CT equipment and CT system |
Also Published As
Publication number | Publication date |
---|---|
WO2018040126A1 (en) | 2018-03-08 |
CN106373165B (en) | 2017-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106373165B (en) | Tomography composograph method for reconstructing and system | |
KR101728046B1 (en) | Tomography apparatus and method for reconstructing a tomography image thereof | |
US9155514B2 (en) | Reconstruction with partially known attenuation information in time of flight positron emission tomography | |
US5552605A (en) | Motion correction based on reprojection data | |
RU2510080C2 (en) | Image processing device, image processing method and long-term information storage medium | |
US20150213633A1 (en) | System, method and computer-accessible medium for providing a panoramic cone beam computed tomography (cbct) | |
US20100172472A1 (en) | Collecting images for image stitching with rotating a radiation detector | |
US20160140749A1 (en) | Methods for generation of edge=preserving synthetic mammograms from tomosynthesis data | |
ES2780152T3 (en) | System and procedure to generate an improved density distribution in three-dimensional bone models | |
US11179128B2 (en) | Methods and systems for motion detection in positron emission tomography | |
US9848844B2 (en) | Iterative reconstruction process | |
CN101529475A (en) | Visualization of 3D images in combination with 2D projection images | |
CN102016923A (en) | Image artifact reduction | |
US20140212018A1 (en) | System optics in at least in one of backprojection and forward projection for model-based iterative reconstruction | |
CN117425433A (en) | Artificial intelligence training using multiple motion pulse X-ray source tomosynthesis imaging systems | |
US9965875B2 (en) | Virtual projection image method | |
US10062168B2 (en) | 5D cone beam CT using deformable registration | |
WO2008063817A2 (en) | A method and system for grouping images in a tomosynthesis imaging system | |
US10089758B1 (en) | Volume image reconstruction using projection decomposition | |
EP3518182B1 (en) | Tilted slices in dbt | |
US9237873B2 (en) | Methods and systems for CT projection domain extrapolation | |
RU2551791C2 (en) | Multi-section alignment of imaging data | |
CN112204607B (en) | Scattering correction for X-ray imaging | |
US10966670B2 (en) | Imaging system and method for dual-energy and computed tomography | |
US20180336708A1 (en) | Poly-energetic reconstruction method for metal artifacts reduction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |