CN105184835B - Mammary gland tomographic image reconstructing process and device - Google Patents

Mammary gland tomographic image reconstructing process and device Download PDF

Info

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
Application number
CN201510583397.6A
Other languages
Chinese (zh)
Other versions
CN105184835A (en
Inventor
张娜
杨乐
陈皓
牛杰
韩业成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201510583397.6A priority Critical patent/CN105184835B/en
Publication of CN105184835A publication Critical patent/CN105184835A/en
Priority to US15/317,382 priority patent/US10140735B2/en
Priority to CN201680053776.7A priority patent/CN108352078B/en
Priority to PCT/CN2016/099061 priority patent/WO2017045618A1/en
Priority to GB1706273.8A priority patent/GB2547360B/en
Priority to US15/460,187 priority patent/US9697623B1/en
Priority to US15/608,935 priority patent/US9875558B2/en
Application granted granted Critical
Publication of CN105184835B publication Critical patent/CN105184835B/en
Priority to US16/199,016 priority patent/US10600214B2/en
Priority to US16/199,025 priority patent/US10586355B2/en
Priority to US16/826,717 priority patent/US11335041B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Mammary gland tomographic image reconstructing process and device
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.
CN201510583397.6A 2015-09-15 2015-09-15 Mammary gland tomographic image reconstructing process and device Active CN105184835B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

Patent Citations (3)

* Cited by examiner, † Cited by third party
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.