CN105184835B - Mammary gland tomographic image reconstructing process and device - Google Patents
Mammary gland tomographic image reconstructing process and device Download PDFInfo
- Publication number
- CN105184835B CN105184835B CN201510583397.6A CN201510583397A CN105184835B CN 105184835 B CN105184835 B CN 105184835B CN 201510583397 A CN201510583397 A CN 201510583397A CN 105184835 B CN105184835 B CN 105184835B
- Authority
- CN
- China
- Prior art keywords
- image
- mammary gland
- artifact
- edge
- intersection point
- 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
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention provides a kind of mammary gland tomographic image reconstructing process and devices.This approach includes the following steps：The image parameter of the mechanical coordinate of mammary gland computer tomography device and projected image is initialized；The projected image is pre-processed；To being filtered back projection by pretreated projected image, mammary gland three-dimension disclocation reconstruction image is generated；Remove the detector edge artifact in the mammary gland three-dimension disclocation reconstruction image；Remove the mammary gland edge artifact in the mammary gland three-dimension disclocation reconstruction image；And remove sawtooth artifact in the mammary gland three-dimension disclocation reconstruction image.The present invention can adapt to the mammary gland computer tomography device in multiray source.
Description
Technical field
The invention mainly relates to mammary gland tomography more particularly to a kind of mammary gland tomographic image reconstructing process and dress
It sets.
Background technology
Breast cancer is all the important diseases for seriously threatening women's health in the world.Breast X-ray photography quilt at present
It is known as the preferred test mode of breast cancer.In recent years, as image documentation equipment is constantly updated, digitlization mammary gland tomography synthesizes skill
The appearance of art, also known as digital galactophore tomography (Digital Breast Tomosynthesis, DBT), makes the morning of breast cancer
Phase detects and diagnosis has further raising.
Digital galactophore tomography is a 3 Dimension Image Technique, it is obtained by a series of shootings from different perspectives
A series of mammary gland three-dimension disclocation of the synthesis comprising high-resolution faultage images is rebuild after low dosage mammary gland projected image is reconstructed
Image.These fault images individually show or are carried out in the form of continuously playing Dynamic Announce.Each faultage image shows mammary gland
Each tomography structure.
Radiographic source used by traditional mammary gland computer tomography device is all single radiographic source.Therefore, traditional mammary gland is disconnected
Tomographic image reconstruction and corresponding artifact correction algorithm and flow are also designed single light source.With x-ray light source skill
The development of art has developed the New Type X-Ray Tube in multiray source.But traditional reconstruction flow cannot be satisfied this
The demand of the novel multiray source x-ray bulb of kind.
Invention content
The technical problem to be solved in the present invention is to provide a kind of mammary gland tomographic image reconstructing sides can adapt to multiray source
Method and device.
In order to solve the above technical problems, the present invention provides a kind of mammary gland tomographic image reconstructing process, include the following steps：
Projected image is pre-processed；To being filtered back projection by pretreated projected image, mammary gland three-dimension disclocation weight is generated
Build image；Remove the detector edge artifact in the mammary gland three-dimension disclocation reconstruction image；Remove the mammary gland three-dimension disclocation reconstruction figure
Mammary gland edge artifact as in；And remove sawtooth artifact in the mammary gland three-dimension disclocation reconstruction image.
In one embodiment of this invention, carrying out pretreated step to the projected image includes：By mammary region from this
It is split in projected image；Negative film operation is carried out to the projected image after segmentation；And to the geometric error of mechanical equipment into
Row correction.
In one embodiment of this invention, to being filtered back projection by pretreated projected image, mammary gland three is generated
Tie up tomographic reconstruction images the step of include：This is filtered by pretreated projected image；To filtered perspective view
As being filtered correction to remove highlighted artifact；X-ray attenuation artifact correction is carried out to the projected image of filtered correction；With
And to carrying out backprojection reconstruction by the projected image of X-ray attenuation artifact correction, generate mammary gland three-dimension disclocation reconstruction image.
In one embodiment of this invention, the step of removing the sawtooth artifact in the mammary gland three-dimension disclocation reconstruction image is wrapped
It includes：It calculates in each Angles Projections image in corresponding faultage image of the detector edge in the mammary gland three-dimension disclocation reconstruction image
Mapping position；It determines in the faultage image, the intersection point and its location information of galactophore image edges of regions and each mapping position；Base
In each intersection point and dislocation information corresponding with each intersection point, detector side in each Angles Projections image is removed in galactophore image region
Sawtooth artifact of the edge at the mapping position in the faultage image, wherein dislocation information corresponding with intersection point refers to and the intersection point pair
The tip for the sawtooth artifact answered is relative to intersection point offset in the horizontal direction.
In one embodiment of this invention, detector edge is calculated in each Angles Projections image in the mammary gland three-dimension disclocation
The step of mapping position in corresponding faultage image in reconstruction image includes：Based on the geometry position between radiographic source, detector
The geometry site between relationship, each Angles Projections image and corresponding faultage image is set, it will be each using image interpolation algorithm
The coordinate of pixel is mapped in corresponding faultage image in Angles Projections image；It is determined based on effective imaging region of detector
Mapping position of the detector edge in corresponding faultage image in each Angles Projections image.
In one embodiment of this invention, which includes：Bilinear interpolation and nearest-neighbor interpolation method.
In one embodiment of this invention, each intersection point and dislocation information corresponding with each intersection point are based in galactophore image region
Sawtooth artifact of the detector edge at the mapping position in the faultage image includes in middle each Angles Projections image of removal：It creates
0-1 mammary gland templates corresponding with the faultage image；In the 0-1 mammary gland templates, it is pseudo- that each sawtooth is removed successively from mediad both sides
Shadow, including：The edge of the corresponding sawtooth artifact of the intersection point of one mapping position is moved inward one and is equal to the corresponding dislocation of the intersection point
Then galactophore image region part between edge after movement and the edge before movement is set to direct exposure by the distance of information
Region part, and redefine the intersection point of galactophore image edges of regions and next mapping position；It steps be repeated alternatively until
Except each sawtooth artifact；And the faultage image is multiplied with the 0-1 mammary gland templates, the faultage image after being corrected.
In one embodiment of this invention, carrying out pretreated step before to projected image further includes：By mammary gland tomography
The mechanical coordinate of photographic equipment and the image parameter of projected image are initialized.
The present invention also proposes a kind of mammary gland tomographic image reconstructing device, including：Preprocessing module, for projected image into
Row pretreatment；Filtered back projection's module connects the preprocessing module, for being filtered instead by pretreated projected image
Projection generates mammary gland three-dimension disclocation reconstruction image；Detector edge artifact cancellation module, connection filtered back projection module, is used for
Remove the detector edge artifact in the mammary gland three-dimension disclocation reconstruction image；Mammary gland edge artifact cancellation module, connects the detection
Device edge artifact cancellation module, for removing the mammary gland edge artifact in the mammary gland three-dimension disclocation reconstruction image；And sawtooth is pseudo-
Shadow cancellation module connects the mammary gland edge artifact cancellation module, for removing the sawtooth in the mammary gland three-dimension disclocation reconstruction image
Artifact.
In one embodiment of this invention, which executes following operation：By mammary region from projected image
It splits；Negative film operation is carried out to the projected image after segmentation；And the geometric error of mechanical equipment is corrected.
In one embodiment of this invention, which executes following operation：Pretreated throwing is passed through to this
Shadow image is filtered；Correction is filtered to remove highlighted artifact to filtered projected image；To filtered correction
Projected image carry out X-ray attenuation artifact correction；And it is thrown the projected image by X-ray attenuation artifact correction is counter
Shadow is rebuild, and mammary gland three-dimension disclocation reconstruction image is generated.
In one embodiment of this invention, which includes：First submodule, for calculating each angle
Spend mapping position of the detector edge in corresponding faultage image in projected image；The second submodule, for determining the tomograph
As in, the intersection point and its location information of galactophore image edges of regions and each mapping position；Third submodule, for being based on each intersection point
And dislocation information corresponding with each intersection point, it is disconnected at this that detector edge in each Angles Projections image is removed in galactophore image region
Sawtooth artifact at mapping position in tomographic image, wherein and the corresponding dislocation information of intersection point refers to sawtooth corresponding with the intersection point
The tip of artifact is relative to intersection point offset in the horizontal direction.
In one embodiment of this invention, which executes following operation to calculate in each Angles Projections image
The mapping position of corresponding faultage image of the detector edge in mammary gland three-dimension disclocation reconstruction image：Based on radiographic source, detector
Between geometry site, the geometry site between each Angles Projections image and corresponding faultage image, use image
The coordinate of pixel in each Angles Projections image is mapped in corresponding faultage image by interpolation algorithm；Based on the effective of detector
Imaging region determines mapping position of the detector edge in corresponding faultage image in each Angles Projections image.
In one embodiment of this invention, which includes：Bilinear interpolation and nearest-neighbor interpolation method.
In one embodiment of this invention, which executes following operation：It creates corresponding with the faultage image
0-1 mammary gland templates；In the 0-1 mammary gland templates, each sawtooth artifact is removed successively from mediad both sides, including：By a mapped bits
The edge for the corresponding sawtooth artifact of intersection point set moves inward a distance for being equal to the corresponding dislocation information of the intersection point, then will move
The galactophore image region part between the edge before edge and movement after dynamic is set to direct exposure area part, and again really
Determine the intersection point of galactophore image edges of regions and next mapping position；It steps be repeated alternatively until each sawtooth artifact of removal；And it will
The faultage image is multiplied with the 0-1 mammary gland templates, the faultage image after being corrected.
In one embodiment of this invention, above-mentioned apparatus further includes：Parameter initialization module connects the preprocessing module,
For the image parameter of the mechanical coordinate of mammary gland computer tomography device and projected image to be initialized.
Compared with prior art, the present invention passes through to the edge artifact and sawtooth artifact in mammary gland three-dimension disclocation reconstruction image
Correction, can adapt to the mammary gland computer tomography device in multiray source.
Description of the drawings
Fig. 1 is the mammary gland tomographic image reconstructing process flow chart of one embodiment of the invention.
Fig. 2 is the projected image pretreatment process figure of method shown in Fig. 1.
Fig. 3 is filtered back projection's flow chart of method shown in Fig. 1.
Fig. 4 is that the sawtooth artifact of method shown in Fig. 1 eliminates flow chart.
Fig. 5 is the single sawtooth artifact correction flow of method shown in Fig. 4.
Fig. 6 is the reconstruction image schematic diagram before sawtooth artifact correction.
Fig. 7 is the sawtooth artifact correction schematic diagram of one embodiment of the invention.
Fig. 8 A-8D are the sawtooth artifact correction schematic diagrames of one embodiment of the invention.
Fig. 9 is the sawtooth artifact correction result of one embodiment of the invention.
Figure 10 is the artifact correction result of existing mammary gland tomographic reconstruction images.
Figure 11 is detector edge artifact schematic diagram.
Figure 12 is the structure chart of the mammary gland computer tomography device of one embodiment of the invention.
Figure 13 is the structure chart of the mammary gland tomographic image reconstructing device of one embodiment of the invention.
Specific implementation mode
For the above objects, features and advantages of the present invention can be clearer and more comprehensible, below in conjunction with attached drawing to the tool of the present invention
Body embodiment elaborates.
Many details are elaborated in the following description to facilitate a thorough understanding of the present invention, still the present invention can be with
Implemented different from other manner described here using other, therefore the present invention do not limited by following public specific embodiment
System.
Figure 12 shows the structure chart of the mammary gland computer tomography device of one embodiment of the invention.With reference to shown in figure 12, this implementation
The mammary gland computer tomography device 1200 of example includes multiple radiographic sources 1210, detector 1220 and plate for forcing 1230.Detector 1220
Breast is clamped from both direction with plate for forcing 1230.Multiple radiographic sources 1210 are worn from angular illumination breast different at a distance, ray
Enter detector 1220 after crossing breast, to leave projected image on detector 1220.
The characteristics of this mammary gland computer tomography device 1200, is, multiple radiographic sources 1210 has been used to irradiate breast, so as to
The breast area for enough covering bigger, radiographic source need not be moved again.
Mammary gland computer tomography device 1200 can obtain the projected image that the irradiation of radiographic source 1210 of each different angle is formed,
It is then transferred to computer support cooperative work equipment 1200, carries out image reconstruction described below.It can be obtained after image reconstruction
Three-dimensional mammary gland tomographic reconstruction images, it includes many faultage images (or sectioning images).Each faultage image can show mammary gland
A tomography.
It will appear various artifacts in the reconstruction process of mammary gland three-dimensional tomographic image.Artifact refers to not having reflection in image
The change of any density of the real difference in attenuation of body.The presence of these artifacts not only influences the quality of image, it is also possible to lead to vacation
It is positive as a result, covering true lesion and causing to fail to pinpoint a disease in diagnosis or mistaken diagnosis.Therefore it is projected image process to eliminate these artifacts
In necessary work.
Fig. 1 is the mammary gland tomographic image reconstructing process flow chart of one embodiment of the invention.Refering to what is shown in Fig. 1, the present embodiment
Mammary gland tomographic image reconstructing process flow 100 include the following steps：
Step 101, parameter initialization is carried out.Specifically, mammary gland computer tomography device 800 to be transmitted through to the mechanical coordinate come
And the image parameter of projected image 111 is initialized.
In step 102, projected image pretreatment is carried out.This step input projected image 112 is pre-processed, and is exported pre-
Projected image that treated.
In step 103, it is filtered back projection.This step inputs pretreated projected image, after being filtered back projection,
Export mammary gland three-dimension disclocation reconstruction image 113 (hereinafter reconstruction image).
In step 104, removal detector edge artifact.
Since detector is too small (with reference to figure 7), as shown in figure 11, can be generated at the edge of each faultage image of reconstruction image
Detector edge artifact 1100.This step 104 inputs reconstruction image 113, removes detector edge puppet movie queen, exports new reconstruction
Image.
In step 105, removal mammary gland edge artifact.
In step 106, sawtooth artifact is removed.
Zigzag artifact can be left at detector edge artifact, which is that scanning angle is small to be made since detector is too small
At, therefore the artifact is removed.The reconstruction image of this step input removal mammary gland edge puppet movie queen, removal sawtooth are pseudo-
Movie queen exports new reconstruction image.
In step 104-106, the step of removing various artifacts, can be directed to entire reconstruction image and carry out, and can also distinguish needle
Each faultage image is carried out.
Fig. 2 shows the projected image pre-treatment step flow charts of method shown in Fig. 1.Refering to what is shown in Fig. 2, this step is further
Include the following steps：
In step 201, mammary region is split from projected image.
In step 202, negative film operation is carried out to the projected image after segmentation.
In step 203, geometric error existing for this equipment is corrected.
Fig. 3 shows filtered back projection's flow chart of steps of method shown in Fig. 1.Refering to what is shown in Fig. 3, this step includes following step
Suddenly：
In step 301, projected image filtering is carried out.This step generates corresponding filter function, to passing through pretreated throwing
Shadow image carries out projected image filtering.
In step 302, projected image filtering and calibration is carried out.A circle can be generated after being filtered due to projected image in organization edge
High bright limb so that high bright limb occurs in corresponding reconstructed result, therefore is filtered correction in step 302, removes and highlights artifact.
In step 303, X-ray attenuation artifact correction is carried out.Due to X-ray attenuation, decaying artifact can be generated on the image,
Therefore decaying artifact correction is carried out in step 303.
In step 304, backprojection reconstruction is carried out to image, generates reconstruction image.
Figure 10 is the artifact correction result of existing mammary gland tomographic reconstruction images.The puppet of existing mammary gland tomographic reconstruction images
Zigzag artifact can occur in region A in shadow correction result, influence to diagnose.The embodiment of the present invention introduces removal sawtooth artifact
The step of eliminate this zigzag artifact.
Fig. 6 is the reconstruction image schematic diagram before sawtooth artifact correction, refering to what is shown in Fig. 6, at the upper and lower edge of reconstruction image
Sawtooth artifact is left, the position being directed toward such as arrow A1, A2.Fig. 4 is the sawtooth artifact correction schematic diagram of one embodiment of the invention.
As shown in figure 4, the method includes the following steps：
In step 401, calculate detector edge in each Angles Projections image reconstruction image corresponding faultage image
Mapping position.
Fig. 7 is the sawtooth artifact correction schematic diagram of one embodiment of the invention, and Fig. 7 top halfs are the top of mammary gland, lower half
Part is the lower part of mammary gland.Refering to what is shown in Fig. 7, mapping position of the detector edge in corresponding faultage image is the level in figure
Dotted line L1-L4.The purpose of this step is to calculate the position of these horizontal dotted lines L1-L4.It these positions can be in subsequent step
Position intersection point.
In step 402, it determines in the faultage image, the intersection point and intersection point of galactophore image edges of regions and each mapping position
Location information.
Shown in Fig. 7 top halfs, in faultage image, grey parts are the corresponding image-regions of breast tissue,
That is galactophore image region；Black portions are the image-region not comprising breast tissue, i.e., direct exposure area.Galactophore image region
The serrate line of demarcation of a band is formed in edge and direct exposure area.Galactophore image edges of regions (at line of demarcation) with it is each
The intersection point of a mapping position (horizontal dotted line) is the position that each arrow is directed toward in figure.It can be seen that these intersection points P1-P4
(being pointed out with arrow in figure) is generally in the same horizontal position with the tips each sawtooth S1-S4.
In step 403, it is based on each intersection point and dislocation information corresponding with each intersection point removes respectively in galactophore image region
Sawtooth artifact of the detector edge at the mapping position in corresponding faultage image in Angles Projections image.
Here, and the corresponding dislocation information of intersection point refer to sawtooth artifact corresponding with the intersection point with the intersection point in the horizontal direction
Offset.As shown in fig. 7, in each point of intersection, due to the presence of sawtooth, galactophore image region is in the boundary with direct exposure area
Place misplaces, therefore once obtains dislocation information, you can removes sawtooth according to this dislocation information.For example, in level
At dotted line L1, sawtooth S1 results in the offset that amplitude is D1.Therefore by deviate the dislocation information that D1 is indicated, you can from mammary gland
Sawtooth S1 is removed in image-region.
In above-mentioned step 401, detector edge in each Angles Projections image can be calculated by following manner
Mapping position in corresponding faultage image：First based on geometry site, each Angles Projections figure between radiographic source, detector
As and the geometry site between corresponding faultage image, using image interpolation algorithm by pixel in each Angles Projections image
Coordinate is mapped in corresponding faultage image；The effective imaging region for being then based on detector is determined and is visited in each Angles Projections image
Survey mapping position of the device edge in corresponding faultage image.Such as assume that projected image resolution ratio is 1000*1000 (i.e. length and width models
Enclose all is 0~1000) size, the pixel in projected image be mapped to after faultage image might have coordinate be negative value and
There is (the corresponding pixel of these coordinates is all the pixel except detector) in pixel of the coordinate more than 1000, pass through coordinate
0 and 1000 the two critical values just can determine that detector edge, the horizontal dotted line in as Fig. 7.It can make in result of calculation
The point coordinates (x, y) above each dotted line is indicated with a two-dimensional array.
It in one embodiment, can be by bilinear interpolation or nearest-neighbor interpolation method by the pixel of projected image
Coordinate be mapped to rebuild each layer in space faultage image in.
Fig. 5 is the single sawtooth artifact correction flow of method shown in Fig. 4.Refering to what is shown in Fig. 5, artifact correction flow includes such as
Lower step：
Step 501,0-1 mammary gland templates corresponding with faultage image are created.
0-1 mammary gland templates be the gray value of the pixel in the galactophore image region in faultage image is set to 1 (white), and
The gray value for the pixel that direct exposure area is set is set to 0 (black).This helps to reduce the resource needed for image procossing.Figure
8A illustrates the part of 0-1 mammary gland templates, including 4 sawtooth S1-S4.
In step 202, in 0-1 mammary gland templates, each sawtooth artifact is removed successively from mediad both sides.
Shown in Fig. 8 A, arrow it is downward be mammary gland center, arrow it is upward be mammary gland one side edge.Such as Fig. 8 A
Shown, by the edge E1 of the corresponding sawtooth artifact S1 of the intersection point P1 of a mapping position L1, inwardly (for the left in figure), movement one is equal to
The distance of the corresponding dislocation information D1 of the point of intersection S 1, the edge E1 ' (as shown in Figure 8 B) after being moved.Then such as Fig. 8 C institutes
Show, it is (black that the gray value of the galactophore image region part between the edge E1 ' after movement and the edge E1 before movement is set to 0
Color).Finally, the intersection point of galactophore image edges of regions and next mapping position is redefined.Compare the position of new intersection point P2 ' in Fig. 8 D
It sets with the position of old acquaintance's point P2 in Fig. 8 A as it can be seen that this intersection point is inwardly moved a distance (for the left in figure), the distance is big
About dislocation information D1.
Sawtooth artifact S2, S3 and S4 are removed successively according to the flow of Fig. 8 A-8D, you can remove all sawtooth in this region
Artifact.
It returns to shown in Fig. 7, the sawtooth artifact of faultage image lower half portion is also similarly removed, to be free of
There are the 0-1 mammary gland templates of sawtooth artifact.
In step 503, faultage image is multiplied with 0-1 mammary gland templates, the faultage image after being corrected.
Fig. 9 is the sawtooth artifact correction of one embodiment of the invention as a result, left side is the part before the removal of sawtooth artifact in Fig. 9
Faultage image, right side are the local fault image after the removal of sawtooth artifact.Refering to what is shown in Fig. 9, the region B1 of left-side images has saw
Tooth artifact has had been removed sawtooth artifact, to obtain after sawtooth artifact correction above-mentioned at the region B2 of image right
The faultage image of pseudo- movie queen must be removed.
Figure 13 is the structure chart of the mammary gland tomographic image reconstructing device of one embodiment of the invention.With reference to shown in figure 13, mammary gland
Tomographic image reconstructing device can be implemented in computer support cooperative work equipment 900 comprising parameter initialization module
1301, preprocessing module 1302, filtered back projection's module 1303, detector edge artifact cancellation module 1304, mammary gland edge are pseudo-
Shadow cancellation module 1305 and sawtooth artifact cancellation module 1306.Parameter initialization module 1301 is for setting mammary gland tomography
The image parameter of standby mechanical coordinate and projected image is initialized.Preprocessing module 1302 connects the parameter initialization mould
Block 1301, for being pre-processed to 111 projected images.Filtered back projection's module 1303 connects preprocessing module 1302, is used for
To being filtered back projection by pretreated projected image 113, reconstruction image is generated.Detector edge artifact cancellation module
1304 connection filtered back projection modules 1303, for removing the detector edge artifact in reconstruction image 113.Mammary gland edge artifact
Cancellation module 1305 connects detector edge artifact cancellation module 1304, for removing the puppet of the mammary gland edge in reconstruction image 113
Shadow.Sawtooth artifact cancellation module 1306 connects mammary gland edge artifact cancellation module 1305, for removing the sawtooth in reconstruction image
Artifact.
The mammary gland tomographic image reconstructing and artifact correction method of the above embodiment of the present invention can such as computer software,
It is implemented in the computer-readable medium of the combination of hardware or computer software and hardware.For hardware implementation, this
Embodiment described in invention can be in one or more application-specific integrated circuits (ASIC), digital signal processor (DSP), number
Signal processor (DAPD), programmable logic device (PLD), field programmable gate array (FPGA), processor, controller,
The selection of microcontroller, microprocessor, other electronic devices for executing above-mentioned function or above-mentioned apparatus combines to be subject to reality
It applies.In some circumstances, this kind of embodiment can be implemented by controller.
For software implementation, embodiment described in the present invention can by such as program module (procedures) and
The independent software modules such as function module (functions) are implemented, wherein each module execute it is one or more this
The function and operation of described in the text.Software code can be implemented by the application software write in properly programmed language,
It can store in memory, be executed by controller or processor.
The mammary gland tomographic image reconstructing and artifact correction device of the above embodiment of the present invention can be in such as computer hardwares
In be implemented.For example, embodiment described in the present invention can be in one or more application-specific integrated circuits (ASIC), number letter
Number processor (DSP), digital signal processing device (DAPD), programmable logic device (PLD), field programmable gate array
(FPGA), processor, controller, microcontroller, microprocessor, other electronic devices for executing above-mentioned function or above-mentioned dress
The selection set combines to be implemented.In some circumstances, this kind of embodiment can be implemented by controller.Namely
It says, one or more modules of mammary gland tomographic image reconstructing and artifact correction device can be implemented with hardware circuit.
Mammary gland tomographic image reconstructing and artifact correction method and device, can be adapted for described in the embodiment of the present invention
The mammary gland computer tomography device of the radiographic source of various quantity.The pattern of these radiographic sources can be linear, arc line shaped or its
Its kenel.Also, the number of these radiographic sources can be one or more.
The mammary gland tomographic image reconstructing process and device of the present invention, passes through the edge artifact and sawtooth artifact to reconstruction image
Correction, can adapt to the mammary gland computer tomography device in multiray source.
Although the present invention is described with reference to current specific embodiment, those of ordinary skill in the art
It should be appreciated that above embodiment is intended merely to illustrate the present invention, can also make in the case of no disengaging spirit of that invention
Go out various equivalent change or replacement, therefore, as long as to the variation of above-described embodiment, change in the spirit of the present invention
Type will all be fallen in the range of following claims.
Claims (16)
1. a kind of mammary gland tomographic image reconstructing process, includes the following steps：
Projected image is pre-processed；
To being filtered back projection by pretreated projected image, mammary gland three-dimension disclocation reconstruction image is generated；
Remove the detector edge artifact in the mammary gland three-dimension disclocation reconstruction image；
Remove the mammary gland edge artifact in the mammary gland three-dimension disclocation reconstruction image；And
Remove the sawtooth artifact in the mammary gland three-dimension disclocation reconstruction image.
2. according to the method described in claim 1, it is characterized in that, including to the pretreated step of projected image progress：
Mammary region is split from the projected image；
Negative film operation is carried out to the projected image after segmentation；And
The geometric error of mechanical equipment is corrected.
3. according to the method described in claim 1, it is characterized in that, to being filtered anti-throwing by pretreated projected image
Shadow, generate mammary gland three-dimension disclocation reconstruction image the step of include：
This is filtered by pretreated projected image；
Correction is filtered to remove highlighted artifact to filtered projected image；
X-ray attenuation artifact correction is carried out to the projected image of filtered correction；And
To carrying out backprojection reconstruction by the projected image of X-ray attenuation artifact correction, mammary gland three-dimension disclocation reconstruction image is generated.
4. according to the method described in claim 1, it is characterized in that, the sawtooth removed in the mammary gland three-dimension disclocation reconstruction image is pseudo-
The step of shadow includes：
Calculate corresponding faultage image of the detector edge in the mammary gland three-dimension disclocation reconstruction image in each Angles Projections image
In mapping position；
It determines in the faultage image, the intersection point and its location information of galactophore image edges of regions and each mapping position；
Based on each intersection point and dislocation information corresponding with each intersection point, removes in galactophore image region and visited in each Angles Projections image
Sawtooth artifact of the device edge at the mapping position in the faultage image is surveyed, wherein dislocation information corresponding with intersection point refers to and is somebody's turn to do
The tip of the corresponding sawtooth artifact of intersection point is relative to intersection point offset in the horizontal direction.
5. according to the method described in claim 4, it is characterized in that, calculating in each Angles Projections image detector edge at this
The step of mapping position in corresponding faultage image in mammary gland three-dimension disclocation reconstruction image includes：
Based between geometry site, each Angles Projections image and the corresponding faultage image between radiographic source, detector
The coordinate of pixel in each Angles Projections image is mapped to corresponding tomograph by geometry site using image interpolation algorithm
As in；
Determine that detector edge is in corresponding faultage image in each Angles Projections image based on effective imaging region of detector
Mapping position.
6. according to the method described in claim 5, it is characterized in that, the image interpolation algorithm includes：Bilinear interpolation and most
Neighbour domain interpolation method.
7. according to the method described in claim 4, it is characterized in that, being existed based on each intersection point and dislocation information corresponding with each intersection point
Sawtooth of the detector edge at the mapping position in the faultage image in each Angles Projections image is removed in galactophore image region
Artifact includes：
Create 0-1 mammary gland templates corresponding with the faultage image；
In the 0-1 mammary gland templates, each sawtooth artifact is removed successively from mediad both sides, including：By the intersection point of a mapping position
The edge of corresponding sawtooth artifact moves inward a distance for being equal to the corresponding dislocation information of the intersection point, then by the side after movement
The galactophore image region part between edge before edge and movement is set to direct exposure area part, and redefines mammary gland figure
As the intersection point of edges of regions and next mapping position；It steps be repeated alternatively until each sawtooth artifact of removal；And
The faultage image is multiplied with the 0-1 mammary gland templates, the faultage image after being corrected.
8. according to the method described in claim 1, it is characterized in that, to projected image also wrap before pretreated step
It includes：
The image parameter of the mechanical coordinate of mammary gland computer tomography device and projected image is initialized.
9. a kind of mammary gland tomographic image reconstructing device, including：
Preprocessing module, for being pre-processed to projected image；
Filtered back projection's module connects the preprocessing module, for being filtered back projection to passing through pretreated projected image,
Generate mammary gland three-dimension disclocation reconstruction image；
Detector edge artifact cancellation module, connection filtered back projection module, for removing the mammary gland three-dimension disclocation reconstruction image
In detector edge artifact；
Mammary gland edge artifact cancellation module connects the detector edge artifact cancellation module, for removing the mammary gland three-dimension disclocation
Mammary gland edge artifact in reconstruction image；And
Sawtooth artifact cancellation module connects the mammary gland edge artifact cancellation module, for removing the mammary gland three-dimension disclocation reconstruction figure
Sawtooth artifact as in.
10. device according to claim 9, which is characterized in that the preprocessing module executes following operation：By mammary region
It is split from projected image；Negative film operation is carried out to the projected image after segmentation；And the geometric error to mechanical equipment
It is corrected.
11. device according to claim 9, which is characterized in that filtered back projection's module executes following operation：To the warp
Pretreated projected image is crossed to be filtered；Correction is filtered to remove highlighted artifact to filtered projected image；It is right
The projected image of filtered correction carries out X-ray attenuation artifact correction；And to the projection by X-ray attenuation artifact correction
Image carries out backprojection reconstruction, generates mammary gland three-dimension disclocation reconstruction image.
12. device according to claim 9, which is characterized in that the sawtooth artifact cancellation module includes：
First submodule, for calculating mapped bits of the detector edge in corresponding faultage image in each Angles Projections image
It sets；
The second submodule, for determining in the faultage image, intersection point and its position of the galactophore image edges of regions with each mapping position
Confidence ceases；
Third submodule removes each for being based on each intersection point and dislocation information corresponding with each intersection point in galactophore image region
Sawtooth artifact of the detector edge at the mapping position in the faultage image in Angles Projections image, wherein corresponding with intersection point
The information that misplaces refers to the tip of sawtooth artifact corresponding with the intersection point relative to intersection point offset in the horizontal direction.
13. device according to claim 12, which is characterized in that it is each to calculate that first submodule executes following operation
The mapping position of corresponding faultage image of the detector edge in mammary gland three-dimension disclocation reconstruction image in Angles Projections image：It is based on
Geometry site, each Angles Projections image between radiographic source, detector and the geometric position between corresponding faultage image
The coordinate of pixel in each Angles Projections image is mapped in corresponding faultage image by relationship using image interpolation algorithm；Base
Mapped bits of the detector edge in corresponding faultage image in each Angles Projections image are determined in effective imaging region of detector
It sets.
14. device according to claim 13, which is characterized in that the image interpolation algorithm includes：Bilinear interpolation and
Nearest-neighbor interpolation method.
15. device according to claim 12, which is characterized in that the third submodule executes following operation：It creates and is somebody's turn to do
The corresponding 0-1 mammary gland template of faultage image；In the 0-1 mammary gland templates, each sawtooth artifact is removed successively from mediad both sides, is wrapped
It includes：The edge of the corresponding sawtooth artifact of the intersection point of one mapping position is moved inward one and is equal to the corresponding dislocation information of the intersection point
Then galactophore image region part between edge after movement and the edge before movement is set to direct exposure area portion by distance
Point, and redefine the intersection point of galactophore image edges of regions and next mapping position；It steps be repeated alternatively until each saw of removal
Tooth artifact；And the faultage image is multiplied with the 0-1 mammary gland templates, the faultage image after being corrected.
16. device according to claim 9, which is characterized in that further include：
Parameter initialization module connects the preprocessing module, is used for the mechanical coordinate of mammary gland computer tomography device and projection
The image parameter of image is initialized.
Priority Applications (10)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510583397.6A CN105184835B (en) | 2015-09-15 | 2015-09-15 | Mammary gland tomographic image reconstructing process and device |
| CN201680053776.7A CN108352078B (en) | 2015-09-15 | 2016-09-14 | Image reconstruction system and method |
| US15/317,382 US10140735B2 (en) | 2015-09-15 | 2016-09-14 | Image reconstruction system and method |
| PCT/CN2016/099061 WO2017045618A1 (en) | 2015-09-15 | 2016-09-14 | Image reconstruction system and method |
| GB1706273.8A GB2547360B (en) | 2015-09-15 | 2016-09-14 | Image reconstruction system and method |
| US15/460,187 US9697623B1 (en) | 2015-09-15 | 2017-03-15 | Image reconstruction system and method |
| US15/608,935 US9875558B2 (en) | 2015-09-15 | 2017-05-30 | Image reconstruction system and method |
| US16/199,016 US10600214B2 (en) | 2015-09-15 | 2018-11-23 | Image reconstruction system and method |
| US16/199,025 US10586355B2 (en) | 2015-09-15 | 2018-11-23 | Image reconstruction system and method |
| US16/826,717 US11335041B2 (en) | 2015-09-15 | 2020-03-23 | Image reconstruction system and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510583397.6A CN105184835B (en) | 2015-09-15 | 2015-09-15 | Mammary gland tomographic image reconstructing process and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105184835A CN105184835A (en) | 2015-12-23 |
| CN105184835B true CN105184835B (en) | 2018-11-06 |
Family
ID=54906887
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510583397.6A Active CN105184835B (en) | 2015-09-15 | 2015-09-15 | Mammary gland tomographic image reconstructing process and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105184835B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2547360B (en) | 2015-09-15 | 2018-02-14 | Shanghai United Imaging Healthcare Co Ltd | Image reconstruction system and method |
| US9697623B1 (en) | 2015-09-15 | 2017-07-04 | Shanghai United Imaging Healthcare Co., Ltd. | Image reconstruction system and method |
| CN106251378B (en) * | 2016-07-22 | 2019-10-08 | 上海联影医疗科技有限公司 | The minimizing technology and device of edge artifact in medical image |
| CN106373165B (en) * | 2016-08-31 | 2017-11-28 | 广州华端科技有限公司 | Tomography composograph method for reconstructing and system |
| US10993689B2 (en) * | 2017-08-31 | 2021-05-04 | General Electric Company | Method and system for motion assessment and correction in digital breast tomosynthesis |
| CN112150426B (en) * | 2020-09-17 | 2022-08-30 | 南京邮电大学 | Digital mammary gland tomosynthesis photography reconstruction method based on nonparametric nuclear density estimation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102768759A (en) * | 2012-07-04 | 2012-11-07 | 深圳安科高技术股份有限公司 | Intraoperative CT (Computed Tomography) image beam hardening artifact correction method and device |
| CN102800073A (en) * | 2012-06-28 | 2012-11-28 | 西北工业大学 | Automatic judgment and correction method of cone beam CT annulus artifact |
| CN104135931A (en) * | 2012-02-22 | 2014-11-05 | 皇家飞利浦有限公司 | Method and system for reducing localized artifacts in imaging data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8712134B2 (en) * | 2011-10-18 | 2014-04-29 | Kabushiki Kaisha Toshiba | Method and system for expanding axial coverage in iterative reconstruction in computer tomography (CT) |
-
2015
- 2015-09-15 CN CN201510583397.6A patent/CN105184835B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104135931A (en) * | 2012-02-22 | 2014-11-05 | 皇家飞利浦有限公司 | Method and system for reducing localized artifacts in imaging data |
| CN102800073A (en) * | 2012-06-28 | 2012-11-28 | 西北工业大学 | Automatic judgment and correction method of cone beam CT annulus artifact |
| CN102768759A (en) * | 2012-07-04 | 2012-11-07 | 深圳安科高技术股份有限公司 | Intraoperative CT (Computed Tomography) image beam hardening artifact correction method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105184835A (en) | 2015-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105184835B (en) | Mammary gland tomographic image reconstructing process and device | |
| CN104545969A (en) | Determining value of recording parameter by use of anatomic landmark | |
| CN105078492B (en) | The bearing calibration of sawtooth artifact and device in digital galactophore cross sectional reconstruction | |
| AU2018301580B2 (en) | Three-dimensional ultrasound image display method | |
| CN109313816B (en) | Depth enhanced tomosynthesis reconstruction | |
| CN105326524B (en) | The medical imaging procedure and device of the artifact in image can be reduced | |
| CN115209808A (en) | Learning model creation method, image generation method, and image processing device | |
| JP7459101B2 (en) | Method and device for generating a panoramic tomographic image of an object to be recorded - Patents.com | |
| Ibraheem | Reduction of artifacts in dental cone beam CT images to improve the three dimensional image reconstruction | |
| US12115017B2 (en) | Method of calibrating x-ray projection geometry in x-ray cone beam computed tomography | |
| KR102662878B1 (en) | Method for reconstructing 2D image from multiple X-ray images | |
| JPWO2020067475A1 (en) | Tomographic image generator, method and program | |
| JP5487172B2 (en) | Medical image diagnostic apparatus and medical image processing apparatus | |
| JPH0759763A (en) | Three-dimensional object measuring system | |
| JP7439075B2 (en) | Device and method for editing panoramic radiographic images | |
| CN104545962A (en) | Medical imaging method and system capable of reducing artifacts in images | |
| KR20230165022A (en) | Cycle-GAN-based System for learning for image quality improvement | |
| JP7532379B2 (en) | Method and apparatus for generating panoramic layered images - Patents.com | |
| US11786193B2 (en) | Metal artifacts reduction in cone beam reconstruction | |
| JP2024500828A (en) | peripheral perfusion measurement | |
| US11317887B2 (en) | Computed tomography reconstruction of moving bodies | |
| CN107495976B (en) | Method and device for acquiring maximum value and gray value image in image reconstruction | |
| JP7326807B2 (en) | Image processing device, display device, tomographic image display system, image processing method, and program | |
| JP6707643B2 (en) | Computation and display of composite mammogram during scan acquisition | |
| CN118177852A (en) | System and method for projection enhancement for 2D composite image generation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |