CN106667512A - Geometric correction method of X-ray imaging equipment and breast tomography equipment - Google Patents
Geometric correction method of X-ray imaging equipment and breast tomography equipment Download PDFInfo
- Publication number
- CN106667512A CN106667512A CN201611248790.0A CN201611248790A CN106667512A CN 106667512 A CN106667512 A CN 106667512A CN 201611248790 A CN201611248790 A CN 201611248790A CN 106667512 A CN106667512 A CN 106667512A
- Authority
- CN
- China
- Prior art keywords
- label
- projected image
- projected
- image
- edge
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012937 correction Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000003384 imaging method Methods 0.000 title claims abstract description 49
- 210000000481 breast Anatomy 0.000 title abstract description 3
- 238000003325 tomography Methods 0.000 title abstract 2
- 238000001514 detection method Methods 0.000 claims description 7
- 210000005075 mammary gland Anatomy 0.000 claims description 7
- 238000003745 diagnosis Methods 0.000 abstract description 3
- 239000003550 marker Substances 0.000 abstract 5
- 229910000831 Steel Inorganic materials 0.000 description 121
- 239000010959 steel Substances 0.000 description 121
- 239000011159 matrix material Substances 0.000 description 20
- 238000007689 inspection Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
- A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
The invention relates to a geometric correction method of X-ray imaging equipment, an image reconstruction method and breast tomography equipment. The geometric correction method comprises the following steps: acquiring space position of a marker; generating a projection image of the marker, and acquiring a projection position of the marker in the projection image according to the projection image; and acquiring a geometric correction relation according to the space position of the marker and the projection position of the marker. According to the technical scheme, the geometric correction method is simple, the speed is high, the degree of accuracy of an acquired imaging geometric correction relation is high, after projection coordinates are transformed into space coordinates on the basis of the imaging geometric correction relation in a reconstruction process, the quality of a reconstructed image meets actual clinic requirements, and missed diagnosis and misdiagnosis can be avoided.
Description
Technical field
The present invention relates to technical field of image processing, more particularly to a kind of geometric correction method of x-ray imaging equipment, figure
As method for reconstructing, mammary gland Laminographic device.
Background technology
For medical imaging devices, the accuracy of radiographic source and detector space geometry position relationship is to cross sectional reconstruction
The degree of accuracy of result plays vital effect.However, the real space position of radiographic source and detector is always due to installing
Precision, feedback accuracy, certainty of measurement etc. introduce error, in addition, with the increase of medical imaging devices service life, radiographic source and
The degree of accuracy of detector locus can also decline, above-mentioned reason may cause the noise of the final image rebuild and obtain compared with
Greatly, edge blurry or there is artifact, if image that this reconstruction is obtained is used to diagnose, it is likely that cause to fail to pinpoint a disease in diagnosis or mistaken diagnosis phenomenon
Generation.
In order to the picture quality after rebuilding is met actual clinical demand, it is necessary to radiographic source and detector it
Between geometry site carry out geometric correction.
Therefore, how the geometry site between radiographic source and detector is corrected, obtains geometric correction relation,
And then rebuild using the geometric correction relation in process of reconstruction to improve the quality of the image after rebuilding, to cause it
Meet actual clinical demand, as one of current problem demanding prompt solution.
The content of the invention
The problem to be solved in the present invention is to provide a kind of geometric correction method of x-ray imaging equipment, with simple accurate
Imaging geometry correction relationship is obtained, and then when causing for the imaging geometry correction relationship to be used for image reconstruction, the image after reconstruction
Quality meet actual clinical demand.
To solve the above problems, technical solution of the present invention provides a kind of geometric correction method of x-ray imaging equipment, bag
Include：
Obtain the locus of label；
The projected image of the label is generated, and the label is obtained in the perspective view according to the projected image
Projected position as in；
The projected position of locus and the label according to the label obtains imaging geometry correction relationship；
Wherein, projection coordinate table of the projected position of the label with the label center in the projected image
Show.
Optionally, obtaining projection coordinate of the label center in the projected image includes：
Obtain edge of the label in the projected image；
The projection coordinate at the label center is obtained according to the edge.
Optionally, the edge for obtaining the label in the projected image includes：
Detect the edge in the projected image；
Initial projections coordinate of the label center in the projected image is determined based on the edge for detecting；
Chosen in the projected image based on initial projections coordinate of the label center in the projected image
Area-of-interest corresponding with the label；And
Edge of the label in the projected image is obtained based on area-of-interest corresponding with the label.
Optionally, the initial projections coordinate based on the label center in the projected image is in the projection
Area-of-interest corresponding with the label is chosen in image to be included：It is the area-of-interest with the initial projections coordinate
Center, and physical dimension based on the label determines the scope of area-of-interest corresponding with the label.
Optionally, it is described that the label is obtained in the perspective view based on area-of-interest corresponding with the label
Edge as in includes：Edge in pair area-of-interest corresponding with the label is fitted to obtain the label
Edge in the projected image.
Optionally, when the label is by row and/or by rows multiple, the geometric correction method also includes：
Judge whether to detect projected position of each label in the projected image；
When projected position of each label in the projected image is not detected by, based on the multiple label
The projected position of locus and the label being had detected that in the projected image is obtained and not examined in the projected image
The projected position of the label for measuring.
To solve the above problems, technical solution of the present invention also provides a kind of image rebuilding method, including：
Imaging geometry correction relationship is obtained using above-mentioned geometric correction method；
Rebuild based on the imaging geometry correction relationship and analytic reconstruction method or iterative reconstruction.
To solve the above problems, technical solution of the present invention also provides a kind of mammary gland Laminographic device, including：Can be located at mark
The X-ray bulb of note thing side, can be located at the detector that label opposite side and the X-ray bulb are oppositely arranged, and place
The label and the support of the mark label locus, the mammary gland Laminographic device also include：
Image generation unit, the X-ray information for being received according to the detector generates the perspective view of the label
Picture；
Detection unit, for detecting projected position of the label in the projected image；
Acquiring unit, the projected position acquisition imaging for the locus according to the label and the label is several
What correction relationship.
Compared with prior art, technical solution of the present invention has advantages below：
The locus of label is obtained, the projected image of the label is generated, according to the space bit of the label
Put the projection coordinate with the label in the projected image to obtain imaging geometry correction relationship, method is simple and speed
Hurry up.
Further, the projection coordinate with the label center in the projected image is the label described
Projected position in projected image, and then locus and the label center according to the label center is in the throwing
The relation between projected position in shadow image obtains imaging geometry correction relationship, and amount of calculation is small, to a certain extent can be with
Improve the speed for obtaining imaging geometry correction relationship.
Further, when projection coordinate of the label center in the projected image is obtained, according to the label
Edge in the projected image is simple to obtain projection coordinate's method at the label center, and in the label for obtaining
The degree of accuracy of the projection coordinate of the heart is high.
Further, when obtaining the edge of the label in the projected image, in first detecting the projected image
Edge, and initial projections coordinate of the label center in the projected image is determined based on the edge for detecting, profit
Area-of-interest corresponding with label is chosen with the initial projections coordinate, area-of-interest corresponding with label is recycled
To obtain edge of the label in the projected image, label can accurately be obtained in projection using this kind of mode
Edge in image, and then projection coordinate of the label center in projected image is determined according to accurate edge, it is final to obtain
Projection coordinate's degree of accuracy it is high.Therefore, in the imaging obtained according to the locus of degree of accuracy projection coordinate high and label
The degree of accuracy of geometric correction relation is high.
Further, imaging geometry correction relationship is obtained with above-mentioned geometric correction method, and based on the imaging geometry school
Positive relationship is rebuild with parsing reconstruction method or iterative reconstruction, due to the space geometry position between radiographic source and detector
Putting has carried out accurate correction, thus rebuild after picture quality meet actual clinical demand.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of the geometric correction method of the x-ray imaging equipment of embodiment of the present invention；
Fig. 2 is the structural representation of the die body of the embodiment of the present invention；
Fig. 3 is projected image of the die body of the embodiment of the present invention under a certain projection angle；
Fig. 4 is edge of a certain steel ball that detects of the embodiment of the present invention in projected image；
Fig. 5 is the schematic diagram of the area-of-interest of the projecting edge comprising a certain steel ball of the embodiment of the present invention；
Fig. 6 is schematic diagram of the steel ball center for detecting of the embodiment of the present invention in the projected image of die body.
Specific embodiment
It is understandable to enable the above objects, features and advantages of the present invention to become apparent, below in conjunction with the accompanying drawings to the present invention
Specific embodiment be described in detail.Elaborate detail in order to fully understand the present invention in the following description.But
It is that the present invention can be implemented with various different from other manner described here, those skilled in the art can be without prejudice to originally
Various changes are done in the case of invention intension.Therefore the present invention is not limited by following public specific embodiment.
As described in prior art, the geometry site between radiographic source and detector is to the image after reconstruction
Quality plays vital effect, (is penetrated to realize being corrected the geometry site between radiographic source and detector
There is theoretical geometry site and actual geometry site, the reason between radiographic source and detector between line source and detector
Refer to the geometry site that should be met between radiographic source and detector according to the design requirement of equipment by geometry site,
In general, theoretical geometric position is to set up coordinate system according to the theoretical position of flat board, and the position of radiographic source is then in the seat
Position under mark system, this coordinate also characterizes position of the radiographic source in theoretical coordinate system.After radiographic source and detector are installed
Geometry site be referred to as actual geometry site, being typically due to the reasons such as installation accuracy can cause radiographic source and detector
Between actual geometry site and theoretical geometry site deviation), the present invention using label locus
And relation between projected position of the label in projected image obtains imaging geometry correction relationship.
Referring to Fig. 1, Fig. 1 is the schematic flow sheet of the geometric correction method of the x-ray imaging equipment of embodiment of the present invention,
As shown in figure 1, the geometric correction method includes：
S101：Obtain the locus of label；
S102：The projected image of the label is generated, and the label is obtained described according to the projected image
Projected position in projected image；
S103：The projected position of locus and the label according to the label obtains imaging geometry correction and closes
System；Wherein, projection coordinate of the projected position of the label with the label center in the projected image represents.
The geometric correction method of embodiment of the present invention can be used to be related to radiographic source and detector and the projection to collecting
The x-ray imaging equipment that image is rebuild, such as：Digital galactophore fault imaging (DBT, Digital Breast
Tomosynthesis) equipment etc., illustrates in the present embodiment so that the x-ray imaging equipment is DBT equipment as an example.
The geometric correction method of embodiment of the present invention is described in detail below in conjunction with specific embodiment, this reality
Apply in example, be arranged in the die body of test with label, having in die body in other words as a example by the label that can position is carried out
Explanation.Specifically, the die body includes matrix and the label being arranged in described matrix, and the label can be to be embedded in institute
State the steel ball in matrix, or filament being embedded in described matrix etc..Described matrix is usually uniform material, and the mark
Remember that the material of thing is different from the absorption coefficient of the material of described matrix, so that can when being imaged to die body
So that readily the label to be detected from projected image.Additionally, the laying of the label in described matrix need to be caused
Under same projection angle, position of each label in the projected image of the die body is misaligned, the mark in described matrix
Remember thing quantity can according to the actual requirements depending on, the quantity of label is conducive to improving the imaging geometry schools of final acquisition
The degree of accuracy of positive relationship.In the present embodiment by the label be in being embedded in described matrix with certain size steel ball as a example by
Illustrate.The label can not also be embedded in the base in other embodiments, but directly by multiple steel balls or thin
Silk composition.Referring to Fig. 2, Fig. 2 is the structural representation of the die body of the embodiment of the present invention, as shown in Fig. 2 the matrix of die body is a length
Cube, the upper surface of cuboid is embedded with 5 row × 4 row totally 20 steel balls, and the lower surface of cuboid is embedded with 5 row × 5 row totally 25
The projected image and 5 rows of the steel ball of 5 row × 4 row in individual steel ball, and the projected image of the die body for the die body shoot acquisition
The projected image of the steel ball of × 5 row is misaligned.
S101 is performed, the locus of label is obtained, it is enterprising the die body is placed on into detector in the present embodiment
During row imaging, the die body carries to be marked on support, and the support to be had, and will be placed on detector with standoff die body
On a precalculated position, now can according on support mark scale know position of the die body in space, additionally, due to
When die body is designed, label position wherein is also what is known, therefore, it can be known in the position in space according to die body
Label is in locus.
Perform S102, generate the projected image of the label, in this implementations due to be embedded in the label described in
Illustrated as a example by die body, therefore, it can that first the die body is placed on the compressing platform of DBT equipment, then at certain
The die body is exposed under one projection angle, to generate the projected image of the projection angle lower mold body, and then is also created
The projected image of label.
Fig. 3 is projected image of the die body of the embodiment of the present invention under a certain projection angle, the present embodiment as described above
The upper surface of the matrix of middle die body embedded in 5 row × 4 row totally 20 steel balls, and lower surface embedded in 5 row × 5 row totally 25 steel balls,
Therefore in the projected image of the die body under a certain projection angle, as shown in Figure 3, top the first row in the projected image of die body
4 density bullet points be 4 steel balls for being embedded in the first row in the steel ball of 5 row × 4 of the body upper surface of die body row, and position
In 5 density bullet points under it then for be embedded in 5 row × 5 of the lower surface of the matrix of die body row steel ball in the first row 5
Individual steel ball.Locus in the base is embedded according to steel ball and the throwing that can know each steel ball in generation is distributed accordingly
Approximate location on shadow image.
After the projected image for obtaining label, can be according to the quality and the demand of reality of the projected image of generation
The projected image is pre-processed accordingly, the pretreatment can include：Smothing filtering, grey scale mapping etc..Smooth filter
Ripple can be smoothed by setting filtering core to the projected image, and noise in projected image is reduced with this.Grey scale mapping
The projected image of the label of generation can be then mapped to the tonal range specified according to actual demand.Projection to generating
Image carries out the degree of accuracy that corresponding pretreatment is conducive to improving projected position of the detection label in the projected image.
Further, since when exposure is performed to the die body, the position that die body is put on the compressing platform is known
, therefore the locus of label is also known, so label can also be first estimated by way of theoretical calculation existing
Position in projected image, can directly to the region where the label estimated in the projected image in actual treatment
Processed, in can so avoiding projected image, the interference of region beyond label to it, and then can improve to mark
The accuracy in detection of thing, it is also possible to improve the speed of detection label.
After the projected image for obtaining label, projected position of the label in the projected image is obtained.This reality
Apply in example, projected position of the label in projected image is represented with projection coordinate of the label center in projected image.
Set up in the coordinate in space by the projection coordinate and the steel ball at each the steel ball center obtained in the projected image
Equation, obtains projection coordinate to the geometric correction relation of imaging space coordinate.
Specifically, projection coordinate of the steel ball center in the projected image is obtained in the present embodiment in the following way.
Edge of the steel ball in the projected image is obtained first, specifically, is obtained in the following way in the present embodiment
Edge of the steel ball in the projected image：
Detect generation steel ball projected image in edge, specifically can be by edge detection algorithm, such as：By ladder
Sobel operators, canny operators, prewitt operators in degree operator etc. detect the edge in the projected image, this implementation
Sobel operators are used to detect the edge in the projected image in example.As shown in figure 4, the side in projected image is detected
After edge, size following according to steel ball and the edge detected in the projected image determine that each steel ball center exists
Initial projections coordinate in projected image, can detect each steel ball center in perspective view by Hough transformation in the present embodiment
Initial projections coordinate as in, specifically, the real radius and magnifying power in the present embodiment according to steel ball are set for Hough
The ballot radius of conversion, determines steel ball edge in projected image in theory further according to the real radius and magnifying power of steel ball
The number of pixels setting that should be accounted for votes threshold value for detecting initial projections coordinate of the steel ball center in projected image.Ballot threshold value
Setting can be according to actual situation depending on, can also set multiple ballot threshold values.In other embodiments can also be using ash
The method of center of gravity is spent to detect initial projections coordinate of each steel ball center in projected image.Due in the projected image
May include other reasonses generation in the edge for detecting is not the edge at steel ball edge, or perhaps some are invalid
Edge, is likely to without the complete edge including steel ball completely, therefore determined by the above method in the edge for detecting in addition
Projection coordinate of the steel ball center in projected image, namely initial projections coordinate is not accurate coordinates, but the seat of estimation
Mark using the projection coordinate of estimation, it is necessary to accurately determine projection coordinate of the steel ball center in projected image.
Specifically, first with the initial projections coordinate of each steel ball of above-mentioned determination chosen in the projected image with often
The corresponding area-of-interest of individual steel ball, should include the edge of correspondence steel ball in each area-of-interest, the purpose of this step is used for smart
Really calculate the projection coordinate of the steel ball center in the projected image.In the present embodiment, can according to the physical dimension of steel ball,
And generation steel ball projected image when magnifying power determine the size of area-of-interest, area-of-interest can be with general
It is circle, rectangle, square etc., corresponding parameter can be determined in practical application according to the shape of area-of-interest, such as feels emerging
When interesting region is for circle, then need to determine circular radius, when being square, then need to determine the foursquare length of side.For example：If
The radius of steel ball is 10 pixels, and magnifying power when generating the projected image of steel ball is 1.1 times, if area-of-interest is pros
Shape, then the foursquare length of side can be 2 × 10 × 1.1+E, and wherein E is extension surplus, can be 4 in the present embodiment, in reality
In E can according to the actual requirements depending on, as long as the edge in being able to ensure that square area at least includes whole steel
Pearl edge.Referring to Fig. 5, Fig. 5 is showing for the area-of-interest of the projecting edge comprising a certain steel ball of the embodiment of the present invention
It is intended to.As shown in figure 5, the region that white dashed line is included in figure is area-of-interest corresponding with a certain steel ball, in white portion
The edge of white is the steel ball edge detected in the projected image.After the size that area-of-interest is determined, with
Initial projections coordinate of the steel ball center in the projected image determined in said process is center interested determines sense
The scope in interest region.
It is determined that after the area-of-interest for containing real edge of the steel ball in projected image, based on corresponding with steel ball
Area-of-interest to obtain edge of the steel ball in the projected image, be exactly pair corresponding with each steel ball specifically
Steel ball edge in area-of-interest is fitted, using least square method come emerging to the corresponding sense of each steel ball in the present embodiment
Steel ball edge in interesting region is fitted, in other embodiments can also be using Newton iteration method come in area-of-interest
Steel ball edge be fitted.
So far edge of the steel ball in projected image is obtained through the above way, is next being thrown according to the steel ball for obtaining
Edge in shadow image obtains the projection coordinate at steel ball center, specifically, according to the shape at the edge after fitting in the present embodiment
To determine projection coordinate of the steel ball center in the projected image.In the present embodiment, the edge after fitting is circle, the ginseng of fitting
Number is set as the radius of round central coordinate of circle and circle, and the central coordinate of circle as steel ball center that fitting is obtained is in the projected image
Projection coordinate.
It should be noted that illustrating that beam perpendicular to the ray in other words of the plane where detector in the present embodiment
Schematic diagram of the steel ball in projected image when source is 0 °, in practical application, the angle of beam can for -10 °, -5 ° ,+
5 ° ,+10 ° etc., now edge of the steel ball in the projected image be approximately oval, it is determined that region of interest corresponding with steel ball
During domain, the shape of the area-of-interest can be rectangle, when being fitted to steel ball edge in the region of interest, can be with
Use the method for ellipse fitting to be fitted to steel ball edge to obtain edge of the steel ball in the projected image, now, intend
The parameter of conjunction can be oval centre coordinate, oval major axis short axle, the oval centre coordinate as steel ball that fitting is obtained
Projection coordinate of the center in the projected image.
Additionally, being by row and multiple by rows for label in the die body shown in Fig. 2 during actually detected
When (label is likely to only by rows in other embodiments, or only by row arrangement), in detection steel ball in the projection
During projected position in image, it is possible that the situation of missing inspection or flase drop, actually detected each steel ball is in projection
Position relationship on image should have certain corresponding relation with each steel ball between the position relationship in space.Still with ray
As a example by the projected image of the die body that source generates when being 0 °, if the right-on words of testing result, the projected image for eventually detecting
In position relationship between each steel ball should as shown in Figure 3 between each steel ball position relationship.In the present embodiment,
In order to obtain accurate geometric correction relation, after the steel ball in the projected image is detected, can first sentence
It is disconnected that whether whether all of steel ball (missing inspection) is detected in the projected image, then judge currently detected steel ball
Projected position it is whether correct (whether flase drop), it is also possible to first judge whether the projected position of currently detected steel ball correct,
Judge whether to detect all of steel ball in the projected image again.In the present embodiment, first to judge whether
Judge whether to be illustrated as a example by flase drop after missing inspection.In the present embodiment, judge especially by the mode of packet and sequence
Currently whether there is missing inspection and flase drop.Physical location due to steel ball in space is known, be therefore, it can first to detecting
Projection coordinate of the steel ball center in projected image is grouped, and Fig. 6 is that the steel ball center for detecting of the embodiment of the present invention exists
Schematic diagram in the projected image of die body, as shown in fig. 6, according to steel ball space laying situation, can be according to X-coordinate in figure
Size the steel ball center shown in figure is divided into 9 groups and is sorted accordingly, as shown in Figure 6：Asterisk marks (the 2nd, 4,
6th, 8 groups) corresponded to embedded body upper surface shown in Fig. 25 row × 4 row steel ball center in the projected image of die body
Projection coordinate, the steel ball center that (the 1st, 3,5,7,9 groups) of circle sign have corresponded to 5 row × 5 row of embedded matrix lower surface exists
Projection coordinate in the projected image of die body, the steel ball not missing inspection by after packet, could be aware that embedded matrix lower surface is embedding
Entering the steel ball of body upper surface has missing inspection, and the number of missing inspection is 6.Now can be according to the geometric position between steel ball in die body
Relation calculates projection coordinate of the steel ball center of missing inspection in the projected image of die body, namely the steel ball center for obtaining missing inspection
Projection coordinate in steel ball projected image.Specifically, due to the steel ball in the present embodiment die body, the upper and lower are respectively
Away from arrangement, that is to say, that steel ball is in the horizontal direction equidistant laying, and steel ball is also vertically equidistant laying,
Undetected steel ball center now can be determined according to the projection coordinate for having been detected as steel ball center in the horizontal direction
With the projection coordinate of vertical direction, the projection coordinate at the actually steel ball center on matrix upper strata can constitute the grid of 5 row × 4 row,
And the projection coordinate at the steel ball center of matrix lower floor can constitute the grid of 5 row × 5 row, the projection coordinate at the steel ball center of missing inspection is
It is the intersecting point coordinate of grid.
Next judge whether projection coordinate of the steel ball center for detecting in the projected image of die body be correct, specifically
Ground, still can judge steel ball center in the projected image of die body by the geometry site between steel ball in die body
Whether projection coordinate is correct, specifically, from the above can be according to the size of X-coordinate in Fig. 6 by the steel ball shown in figure
The heart is divided into 9 groups and is sorted accordingly, and in grouping process, according to the projection coordinate at the steel ball center having detected that, and
Position relationship and magnifying power between adjacent steel ball can generate the grid constituted with the projection coordinate at steel ball center, each in grid
The coordinate of intersection point is the projection coordinate at steel ball center, therefore when the projection coordinate at steel ball center is judged either with or without flase drop, can
Directly to be judged accordingly according to the lattice relationship for being formed, if detecting projection coordinate and the grid at steel ball center
The ratio of intersecting point coordinate range difference is larger, then judge the projection coordinate's mistake at steel ball center for detecting.In the steel that judgement is detected
After projection coordinate's mistake at pearl center, the projection coordinate will not participate in follow-up calculating.
It should be noted that when the geometric correction relation of projection coordinate and aerial image coordinate is actually calculated, and differ
Surely the projection coordinate at all of steel ball center is needed to participate in calculating, as long as the correct steel ball center for detecting is in the perspective view
Projection coordinate as in is enough, can be without judging whether missing inspection again.
So far projection coordinate of the steel ball center in the projected image detected by the above method, next performs
S103, the space coordinates according to each steel ball center obtains imaging geometry correction with the projection coordinate at each steel ball center and closes
System, imaging geometry correction relationship is obtained in the present embodiment by setting up equation therebetween.In particular by equation below
Obtain：
[u,v,1]T=P [x, y, z, 1]T (1)
Wherein (u, v) is the projection coordinate at the steel ball center detected in the projected image, and (x, y, z) is in steel ball
The space coordinates of the heart；
P is transformation matrix, P=K [R | t], K for 3 × 3 the intrinsic matrix of system, R is 3 × 3 spin matrix, and t is 3 × 1
Translation matrix.Bring the projection coordinate of the space coordinates at all steel ball centers and all steel ball centers into equation (1), foundation side
Journey group, solves each coefficient in equation group, then the imaging geometry correction that can be obtained between projection coordinate and space coordinates is closed
System.
Through the above way, projection coordinate to the geometric correction relation of imaging space coordinate is obtained, next can be by
The imaging geometry correction relationship for obtaining through the above way is used in image reconstruction, to cause that the image after rebuilding meets actual
Clinical demand.
The present embodiment also provides a kind of image rebuilding method, first passes through above-mentioned geometric correction method and obtains imaging geometry school
Positive relationship, namely obtain the three dimensions how projection coordinate in projected image being correctly mapped under theoretical coordinate system
Coordinate, next can use analytic reconstruction method or iterative reconstruction to be rebuild with the image after being rebuild.This implementation
In example, by projected image between projection coordinate and space coordinates the correction of relation realize to radiographic source and detector it
Between geometrical relationship correction, and then rebuild based on the space coordinates after correction, the quality of the image after reconstruction meets reality
The clinical demand on border, avoids missing inspection or the generation of flase drop to a certain extent.
The embodiment of the present invention also provides a kind of mammary gland Laminographic device, including：Can be located at the X-ray ball of label side
Pipe, can be located at the detector that label opposite side and the X-ray bulb are oppositely arranged, and place the label and mark
The support of the label locus, the mammary gland Laminographic device also includes：
Image generation unit, the X-ray information for being received according to the detector generates the perspective view of the label
Picture；
Detection unit, for detecting projected position of the label in the projected image；
Acquiring unit, the projected position acquisition imaging for the locus according to the label and the label is several
What correction relationship.
In sum, geometric correction method, the image reconstruction side of the x-ray imaging equipment that embodiment of the present invention is provided
Method, at least has the advantages that：
The locus of label is obtained, the projected image of the label is generated, according to the space bit of the label
Put the projection coordinate with the label in the projected image to obtain imaging geometry correction relationship, method is simple and speed
Hurry up.
Further, the projection coordinate with the label center in the projected image is the label described
Projected position in projected image, and then locus and the label center according to the label center is in the throwing
The relation between projected position in shadow image obtains imaging geometry correction relationship, and amount of calculation is small, to a certain extent can be with
Improve the speed for obtaining imaging geometry correction relationship.
Further, when projection coordinate of the label center in the projected image is obtained, according to the label
Edge in the projected image is simple to obtain projection coordinate's method at the label center, and in the label for obtaining
The degree of accuracy of the projection coordinate of the heart is high.
Further, when obtaining the edge of the label in the projected image, in first detecting the projected image
Edge, and initial projections coordinate of the label center in the projected image is determined based on the edge for detecting, profit
Area-of-interest corresponding with label is chosen with the initial projections coordinate, area-of-interest corresponding with label is recycled
To obtain edge of the label in the projected image, label can accurately be obtained in projection using this kind of mode
Edge in image, and then when determining projection coordinate of the label center in projected image according to accurate edge, finally obtain
The projection coordinate's degree of accuracy for obtaining is high.Therefore, obtained according to the locus of degree of accuracy projection coordinate high and label center
Imaging geometry correction relationship the degree of accuracy it is high.
Further, imaging geometry correction relationship is obtained with above-mentioned geometric correction method, and based on the imaging geometry school
Positive relationship is rebuild with parsing reconstruction method or iterative reconstruction, due to the space geometry position between radiographic source and detector
Putting has carried out accurate correction, thus rebuild after picture quality meet actual clinical demand.
Although the present invention is disclosed as above with preferred embodiment, it is not for limiting the present invention, any this area
Technical staff without departing from the spirit and scope of the present invention, may be by the methods and techniques content of the disclosure above to this hair
Bright technical scheme makes possible variation and modification, therefore, every content without departing from technical solution of the present invention, according to the present invention
Any simple modification, equivalent variation and modification for being made to above example of technical spirit, belong to technical solution of the present invention
Protection domain.
Claims (9)
1. a kind of geometric correction method of x-ray imaging equipment, it is characterised in that including：
Obtain the locus of label；
The projected image of the label is generated, and according to the projected image acquisition label in the projected image
Projected position；
The projected position of locus and the label according to the label obtains imaging geometry correction relationship.
2. geometric correction method according to claim 1, it is characterised in that the projected position of the label is the mark
Projection coordinate of the Ji Wu centers in the projected image, the locus of the label is the space at the label center
Position.
3. geometric correction method according to claim 1, it is characterised in that obtain the label center in the projection
Projection coordinate in image includes：
Obtain edge of the label in the projected image；
The projection coordinate at the label center is obtained according to the edge.
4. geometric correction method according to claim 3, it is characterised in that the acquisition label is in the projection
Edge in image includes：
Detect the edge in the projected image；
Initial projections coordinate of the label center in the projected image is determined based on the edge for detecting；
Chosen and institute in the projected image based on initial projections coordinate of the label center in the projected image
State the corresponding area-of-interest of label；And
Edge of the label in the projected image is obtained based on area-of-interest corresponding with the label.
5. geometric correction method according to claim 4, it is characterised in that it is described based on the label center described
Initial projections coordinate in projected image chooses area-of-interest corresponding with the label in the projected image to be included：
With the center that the initial projections coordinate is the area-of-interest, and physical dimension based on the label determine with it is described
The scope of the corresponding area-of-interest of label.
6. geometric correction method according to claim 5, it is characterised in that described based on sense corresponding with the label
Interest region obtains edge of the label in the projected image to be included：Pair region of interest corresponding with the label
Edge in domain is fitted to obtain edge of the label in the projected image.
7. geometric correction method as claimed in claim 1, it is characterised in that the label is by row and/or by rows
It is multiple when, the geometric correction method also includes：
Judge whether to detect projected position of each label in the projected image；
When projected position of each label in the projected image is not detected by, the space based on the multiple label
The projected position of position and the label being had detected that in the projected image is obtained and is not detected by the projected image
Label projected position.
8. a kind of image rebuilding method, it is characterised in that including：
Imaging geometry correction relationship is obtained using the geometric correction method described in any one of claim 1~6；
Rebuild based on the imaging geometry correction relationship and analytic reconstruction method or iterative reconstruction.
9. a kind of mammary gland Laminographic device, including the X-ray bulb of label side is can be located at, can be located at label another
The detector that side and the X-ray bulb are oppositely arranged, and place the label and mark the label locus
Support, it is characterised in that also include：
Image generation unit, the X-ray information for being received according to the detector generates the projected image of the label；
Detection unit, for detecting projected position of the label in the projected image；
Acquiring unit, the projected position for the locus according to the label and the label obtains imaging geometry school
Positive relationship；Wherein, projection coordinate table of the projected position of the label with the label center in the projected image
Show.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611248790.0A CN106667512B (en) | 2016-12-29 | 2016-12-29 | Geometric correction method of X-ray imaging device and breast tomography device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611248790.0A CN106667512B (en) | 2016-12-29 | 2016-12-29 | Geometric correction method of X-ray imaging device and breast tomography device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106667512A true CN106667512A (en) | 2017-05-17 |
| CN106667512B CN106667512B (en) | 2020-10-09 |
Family
ID=58872535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611248790.0A Active CN106667512B (en) | 2016-12-29 | 2016-12-29 | Geometric correction method of X-ray imaging device and breast tomography device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106667512B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109452950A (en) * | 2018-12-29 | 2019-03-12 | 上海联影医疗科技有限公司 | A kind of geometric markers method of die body and imaging device |
| CN109658465A (en) * | 2018-12-07 | 2019-04-19 | 广州华端科技有限公司 | Data processing, image rebuilding method and device in image reconstruction process |
| CN110648739A (en) * | 2019-11-19 | 2020-01-03 | 河北医科大学第一医院 | Appendix CT image reporting method based on spatial anatomical position |
| CN110766629A (en) * | 2019-10-17 | 2020-02-07 | 广州华端科技有限公司 | CBCT system geometric correction method, device, computer equipment and storage medium |
| CN111544020A (en) * | 2020-04-17 | 2020-08-18 | 北京东软医疗设备有限公司 | Geometric correction method and device for X-ray imaging equipment |
| CN112294349A (en) * | 2020-10-29 | 2021-02-02 | 上海联影医疗科技股份有限公司 | X-ray equipment correction method, imaging method and imaging system |
| CN112529820A (en) * | 2020-12-29 | 2021-03-19 | 广州华端科技有限公司 | Half-fan-beam geometric correction method, device, equipment and medium for CBCT system |
| CN112634353A (en) * | 2020-12-17 | 2021-04-09 | 广州华端科技有限公司 | Self-calibration method, device and medium for geometric calibration die body of CBCT system |
| WO2021104153A1 (en) * | 2019-11-28 | 2021-06-03 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for automated calibration |
| CN113926176A (en) * | 2020-07-13 | 2022-01-14 | Sg游戏公司 | Gaming environment tracking system calibration |
| CN115018720A (en) * | 2022-05-18 | 2022-09-06 | 深圳市安健科技股份有限公司 | Geometric correction method and device for multilayer flat panel detector |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0623884A2 (en) * | 1993-05-05 | 1994-11-09 | Philips Patentverwaltung GmbH | Process and device for determining the amount of distortion in radiographic images |
| US5442674A (en) * | 1993-01-27 | 1995-08-15 | Ge Medical Systems | Device and automatic method for the geometrical calibration of an X-ray imaging system |
| CN101750021A (en) * | 2009-12-04 | 2010-06-23 | 深圳先进技术研究院 | Calibration method, device and calibration phantomof geometric parameters in CT (computer tomography) system |
| CN103006251A (en) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | Calibration phantom, calibration device and calibration method for calibrating geometric parameters in CT (Computed Tomography) system |
| CN103489181A (en) * | 2013-09-03 | 2014-01-01 | 上海逸动医学科技有限公司 | Perspective X-ray image correcting device and correcting method thereof |
| CN103729827A (en) * | 2012-10-11 | 2014-04-16 | 南京普爱射线影像设备有限公司 | Method for correcting three-dimensional reconstruction overlapped artifact of image intensifier C-shaped arm X-ray system |
| CN104783824A (en) * | 2014-01-20 | 2015-07-22 | 上海联影医疗科技有限公司 | Calibration method for X-ray imaging system |
| EP2952227A1 (en) * | 2014-06-03 | 2015-12-09 | Qualiformed | Object-test for correction of imager spurious motion output fitted to an apparatus for treatment by external radiotherapy when this is in motion |
| CN105960204A (en) * | 2013-11-29 | 2016-09-21 | 普兰梅德有限公司 | Mammography imaging arrangement for tomosynthesis |
-
2016
- 2016-12-29 CN CN201611248790.0A patent/CN106667512B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5442674A (en) * | 1993-01-27 | 1995-08-15 | Ge Medical Systems | Device and automatic method for the geometrical calibration of an X-ray imaging system |
| EP0623884A2 (en) * | 1993-05-05 | 1994-11-09 | Philips Patentverwaltung GmbH | Process and device for determining the amount of distortion in radiographic images |
| CN101750021A (en) * | 2009-12-04 | 2010-06-23 | 深圳先进技术研究院 | Calibration method, device and calibration phantomof geometric parameters in CT (computer tomography) system |
| CN103729827A (en) * | 2012-10-11 | 2014-04-16 | 南京普爱射线影像设备有限公司 | Method for correcting three-dimensional reconstruction overlapped artifact of image intensifier C-shaped arm X-ray system |
| CN103006251A (en) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | Calibration phantom, calibration device and calibration method for calibrating geometric parameters in CT (Computed Tomography) system |
| CN103489181A (en) * | 2013-09-03 | 2014-01-01 | 上海逸动医学科技有限公司 | Perspective X-ray image correcting device and correcting method thereof |
| CN105960204A (en) * | 2013-11-29 | 2016-09-21 | 普兰梅德有限公司 | Mammography imaging arrangement for tomosynthesis |
| CN104783824A (en) * | 2014-01-20 | 2015-07-22 | 上海联影医疗科技有限公司 | Calibration method for X-ray imaging system |
| EP2952227A1 (en) * | 2014-06-03 | 2015-12-09 | Qualiformed | Object-test for correction of imager spurious motion output fitted to an apparatus for treatment by external radiotherapy when this is in motion |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109658465A (en) * | 2018-12-07 | 2019-04-19 | 广州华端科技有限公司 | Data processing, image rebuilding method and device in image reconstruction process |
| CN109452950A (en) * | 2018-12-29 | 2019-03-12 | 上海联影医疗科技有限公司 | A kind of geometric markers method of die body and imaging device |
| CN110766629B (en) * | 2019-10-17 | 2022-03-01 | 广州华端科技有限公司 | CBCT system geometric correction method, device, computer equipment and storage medium |
| CN110766629A (en) * | 2019-10-17 | 2020-02-07 | 广州华端科技有限公司 | CBCT system geometric correction method, device, computer equipment and storage medium |
| CN110648739A (en) * | 2019-11-19 | 2020-01-03 | 河北医科大学第一医院 | Appendix CT image reporting method based on spatial anatomical position |
| CN110648739B (en) * | 2019-11-19 | 2023-02-28 | 河北医科大学第一医院 | Appendix CT image reporting method based on spatial anatomical position |
| US11461929B2 (en) | 2019-11-28 | 2022-10-04 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for automated calibration |
| WO2021104153A1 (en) * | 2019-11-28 | 2021-06-03 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for automated calibration |
| US11676305B2 (en) | 2019-11-28 | 2023-06-13 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for automated calibration |
| CN111544020A (en) * | 2020-04-17 | 2020-08-18 | 北京东软医疗设备有限公司 | Geometric correction method and device for X-ray imaging equipment |
| CN111544020B (en) * | 2020-04-17 | 2023-08-01 | 北京东软医疗设备有限公司 | Geometric correction method and device for X-ray imaging equipment |
| CN113926176A (en) * | 2020-07-13 | 2022-01-14 | Sg游戏公司 | Gaming environment tracking system calibration |
| CN112294349A (en) * | 2020-10-29 | 2021-02-02 | 上海联影医疗科技股份有限公司 | X-ray equipment correction method, imaging method and imaging system |
| CN112634353A (en) * | 2020-12-17 | 2021-04-09 | 广州华端科技有限公司 | Self-calibration method, device and medium for geometric calibration die body of CBCT system |
| CN112634353B (en) * | 2020-12-17 | 2024-03-26 | 广州华端科技有限公司 | Self-calibration method, device and medium for geometric calibration die body of CBCT (computed tomography) system |
| CN112529820A (en) * | 2020-12-29 | 2021-03-19 | 广州华端科技有限公司 | Half-fan-beam geometric correction method, device, equipment and medium for CBCT system |
| CN115018720A (en) * | 2022-05-18 | 2022-09-06 | 深圳市安健科技股份有限公司 | Geometric correction method and device for multilayer flat panel detector |
| CN115018720B (en) * | 2022-05-18 | 2024-09-10 | 深圳市安健科技股份有限公司 | Geometric correction method and device for multi-layer flat panel detector |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106667512B (en) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106667512A (en) | Geometric correction method of X-ray imaging equipment and breast tomography equipment | |
| CN100528078C (en) | Method and system for airway measurement | |
| Mutaf et al. | The impact of temporal inaccuracies on 4DCT image quality | |
| US11935229B2 (en) | Automated scan quality monitoring system | |
| CN102265307B (en) | Mr segmentation using transmission data in hybrid nuclear/mr imaging | |
| CN103181775B (en) | For detecting the method and system of patient body's cursor position | |
| US8000435B2 (en) | Method and system for error compensation | |
| US20090202127A1 (en) | Method And System For Error Compensation | |
| CN106688014B (en) | Apparatus and method for the picture quality for determining radiographic image | |
| CN105931202A (en) | Geometric correction die body correction method and system | |
| CN106548492A (en) | Determine method and device, the image acquiring method of matching double points | |
| CN109419526A (en) | Method and system for locomotion evaluation and correction in the synthesis of digital breast tomography | |
| US20220313176A1 (en) | Artificial Intelligence Training with Multiple Pulsed X-ray Source-in-motion Tomosynthesis Imaging System | |
| Hamill et al. | Respiratory‐gated CT as a tool for the simulation of breathing artifacts in PET and PET/CT | |
| CN109671128A (en) | Data processing, image rebuilding method and device in image reconstruction process | |
| George et al. | Robust automatic rigid registration of MRI and X‐ray using external fiducial markers for XFM‐guided interventional procedures | |
| Shafer et al. | The quantitative potential for breast tomosynthesis imaging | |
| CN109146857A (en) | The Lung neoplasm self-sensing method of chest CT images | |
| CN102081051A (en) | Computed tomography method, computer software, computing device and computed tomography system for determining a volumetric representation of a sample | |
| CN108399354A (en) | The method and apparatus of Computer Vision Recognition tumour | |
| CN106373164A (en) | Method and application for eliminating microscopic CT image geometric artifact | |
| Ulrich et al. | Automated model‐based quantitative analysis of phantoms with spherical inserts in FDG PET scans | |
| CN114494076A (en) | Geometric correction method, image reconstruction method and device | |
| Robert et al. | The geometric calibration of cone–beam systems with arbitrary geometry | |
| CN110151211B (en) | Method for improving scattering estimation accuracy of scanning visual field boundary part in imaging system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| 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 Jiading District Industrial Zone Jiading Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |