CN111445397B - Flat panel detector ghost correction method and device, storage medium and medical equipment - Google Patents
Flat panel detector ghost correction method and device, storage medium and medical equipment Download PDFInfo
- Publication number
- CN111445397B CN111445397B CN202010158713.6A CN202010158713A CN111445397B CN 111445397 B CN111445397 B CN 111445397B CN 202010158713 A CN202010158713 A CN 202010158713A CN 111445397 B CN111445397 B CN 111445397B
- Authority
- CN
- China
- Prior art keywords
- ghost
- exposure
- unsaturated
- image
- saturated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012937 correction Methods 0.000 title claims abstract description 25
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 181
- 238000003702 image correction Methods 0.000 claims abstract description 8
- 238000004590 computer program Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 206010047571 Visual impairment Diseases 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000003705 background correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 101100074187 Caenorhabditis elegans lag-1 gene Proteins 0.000 description 1
- 101100510615 Caenorhabditis elegans lag-2 gene Proteins 0.000 description 1
- 101100134058 Caenorhabditis elegans nth-1 gene Proteins 0.000 description 1
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 1
- 229910021417 amorphous silicon Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- 239000011669 selenium Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10124—Digitally reconstructed radiograph [DRR]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P70/00—Climate change mitigation technologies in the production process for final industrial or consumer products
- Y02P70/50—Manufacturing or production processes characterised by the final manufactured product
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
The application provides a flat panel detector ghost correction method, a device, a storage medium and medical equipment, which are used for rapidly and accurately correcting the flat panel detector ghost. In the method, for any group of DR images which are subjected to saturation exposure before a DR image of a current frame, the residual image value of the bright field image of the DR image of the current frame is determined according to the determined saturation exposure residual image compensation coefficient, the bright field image in the group of DR images and the determined final continuous saturated residual image attenuation curve; for any group of DR images subjected to unsaturated exposure before, determining the residual image value of the bright field image of the DR images in the current frame according to the determined unsaturated exposure residual image compensation coefficient, dark field images in the group of DR images and the determined final continuous unsaturated residual image attenuation curve; and then subtracting the sum of residual image values of the bright field images in each group of DR images in the current frame DR image from the current frame DR image to obtain the DR image after residual image correction.
Description
Technical Field
The application relates to the technical field of medical images, in particular to a method and a device for correcting residual shadows of a flat panel detector, a storage medium and medical equipment.
Background
With the continuous development of medical technology, digital X-ray imaging systems (DR imaging systems) are widely used, and flat panel detectors are key devices in DR imaging systems. The flat panel detector mainly comprises an amorphous silicon flat panel detector and an amorphous selenium flat panel detector according to the materials, but has a common problem of ghost shadow, namely, due to the physical characteristics of the flat panel detector, the signal of the X-ray received by the flat panel detector after exposure is not immediately disappeared under a certain dose exposure, but gradually decays and disappears along with the time, if the exposure is carried out again before the residual signal is not disappeared, a 'shadow' of the last exposure exists on an exposure image, namely, the ghost shadow can cause misdiagnosis of clinical diagnosis.
In order to solve the problem of the ghost shadow of the flat panel detector, two methods of hardware and software can be used at present, the ghost shadow can be removed through hardware methods such as hardware design or improvement process, but the cost can be greatly improved, the effect is not particularly ideal, the software method for removing the ghost shadow not only saves cost and manpower, but also has strong flexibility, and can be applicable to all products of the same material, so the method is generally favored, but the software method in the related technology has high algorithm complexity, much time consumption, the removal effect is not particularly ideal, and excessive correction or insufficient correction is easily caused.
Disclosure of Invention
In view of the foregoing, the present application provides a method, an apparatus, a storage medium and a medical device for correcting a flat panel detector ghost.
In a first aspect, an embodiment of the present application provides a method for correcting a residual image of a flat panel detector, where the method includes:
obtaining a current frame DR image obtained by scanning a detected body and a plurality of groups of previous continuous DR images, wherein each group of DR images comprises a bright field image and a dark field image;
judging whether exposure corresponding to the DR images belongs to saturated exposure or unsaturated exposure for any one of the continuous DR images;
if the image is saturated exposure, determining a residual image value of the bright field image in the DR image of the current frame according to the determined saturated exposure residual image compensation coefficient, the bright field image in the DR image of the group and the determined final continuous saturated residual image attenuation curve;
if the image is unsaturated exposure, determining a residual image value of a bright field image in the DR image of the group of DR images in the current frame according to the determined unsaturated exposure residual image compensation coefficient, a dark field image in the DR image of the group of DR images and the determined final continuous unsaturated residual image attenuation curve;
Subtracting the sum of residual image values of the bright field images in each group of DR images in the current frame DR image from the current frame DR image to obtain a DR image after residual image correction.
In one possible implementation, the determination of the final continuous saturated/unsaturated ghost decay curve includes the steps of:
acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load saturated exposure/unsaturated exposure, and collecting a dark field image at each set time interval in a set time period after exposure;
generating a discrete saturated/unsaturated residual image attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load saturated exposure/unsaturated exposure and each dark field image after exposure;
fitting the discrete saturated/unsaturated ghost attenuation curve to obtain the final continuous saturated/unsaturated ghost attenuation curve.
In one possible implementation manner, the fitting the discrete saturated/unsaturated ghost attenuation curve to obtain the final continuous saturated/unsaturated ghost attenuation curve includes:
fitting the discrete saturated/unsaturated ghost attenuation curve to obtain an original continuous saturated/unsaturated ghost attenuation curve;
Determining a saturated/unsaturated decay time threshold corresponding to the original continuous saturated/unsaturated ghost decay curve according to the original continuous saturated/unsaturated ghost decay curve and a set threshold;
and correcting the original continuous saturated/unsaturated ghost attenuation curve according to the saturated/unsaturated attenuation time threshold to obtain the final continuous saturated/unsaturated ghost attenuation curve.
In one possible implementation manner, the correcting the original continuous saturated/unsaturated ghost attenuation curve according to the saturated/unsaturated attenuation time threshold includes:
maintaining the original continuous saturated/unsaturated ghost decay curve when the ghost decay time is less than or equal to the saturated/unsaturated decay time threshold;
and when the ghost decay time is greater than the saturated/unsaturated decay time threshold, replacing the ghost values of points on the original continuous saturated/unsaturated ghost decay curve which are greater than the decay time threshold with the ghost values corresponding to the saturated/unsaturated decay time threshold.
In one possible implementation manner, the determination of the saturated exposure residual image compensation coefficient includes the following steps:
Collecting a first bright field image under the condition of no-load saturated exposure;
collecting a bright field image under the condition of loaded saturated exposure;
collecting a second bright field image under the condition of no-load saturated exposure again;
determining a ghost value in the second bright-field image according to a pending saturated exposure ghost compensation coefficient, the bright-field image under the condition of loaded saturated exposure and the continuous saturated ghost attenuation curve;
and subtracting the residual image value in the second bright field image from the second bright field image to approximate the first bright field image, and reversely calculating the saturated exposure residual image compensation coefficient.
In one possible implementation manner, the determining of the unsaturated exposure residual compensation coefficient includes the following steps:
collecting a third bright field image under the condition of no-load unsaturated exposure;
collecting a bright field image under the condition of loaded unsaturated exposure and a dark field image after exposure;
acquiring a fourth bright field image under the condition of no-load unsaturated exposure again;
determining a ghost value in the fourth bright-field image according to a to-be-determined unsaturated exposure ghost compensation coefficient, the dark-field image subjected to loaded unsaturated exposure and the continuous unsaturated ghost attenuation curve;
And subtracting the residual image value in the fourth bright field image from the fourth bright field image to approximate the third bright field image, and reversely calculating the unsaturated exposure residual image compensation coefficient.
In a second aspect, embodiments of the present application further provide a flat panel detector ghost correction apparatus, including a module for performing the flat panel detector ghost correction method in the first aspect or any possible implementation manner of the first aspect.
In a third aspect, embodiments of the present application further provide a storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the flat panel detector ghost correction method of the first aspect or any possible implementation manner of the first aspect.
In a fourth aspect, embodiments of the present application further provide a medical device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the flat panel detector ghost correction method in the first aspect or any possible implementation manner of the first aspect when the program is executed.
The technical scheme provided by the embodiment of the application at least brings the following beneficial effects:
according to the technical scheme provided by the application, for the image subjected to saturation exposure before the current frame DR image, the residual image value of the saturated exposure image remained in the current frame DR image is determined according to the determined saturated exposure residual image compensation coefficient, the bright field image in the group of DR images and the determined final continuous saturated residual image attenuation curve, and for the image subjected to unsaturated exposure before the current frame DR image, the residual image value of the unsaturated exposure image remained in the current frame DR image is determined according to the determined unsaturated exposure residual image compensation coefficient, the dark field image in the group of DR images and the determined final continuous unsaturated residual image attenuation curve, that is, according to different scanning conditions, different models and parameters are respectively applied to carry out residual image removal, so that the residual image of the flat panel detector can be quickly and accurately corrected.
Drawings
Fig. 1 is a schematic flow chart of a method for correcting a residual image of a flat panel detector according to an embodiment of the present application;
fig. 2 is a schematic structural diagram of a flat panel detector ghost correction device according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a second structure of a device for correcting residual shadows of a flat panel detector according to an embodiment of the present disclosure;
fig. 4 is a schematic diagram of a third structure of the flat panel detector ghost correction device according to the embodiment of the present application;
fig. 5 is a schematic structural diagram of a medical device according to an embodiment of the present application.
Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples are not representative of all implementations consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with some aspects of the present application as detailed in the accompanying claims.
The terminology used in the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the present application. As used in this application and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any or all possible combinations of one or more of the associated listed items.
It should be understood that although the terms first, second, third, etc. may be used herein to describe various information, these information should not be limited by these terms. These terms are only used to distinguish one type of information from another. For example, a first message may also be referred to as a second message, and similarly, a second message may also be referred to as a first message, without departing from the scope of the present application. The word "if" as used herein may be interpreted as "at … …" or "at … …" or "responsive to a determination", depending on the context.
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application.
Referring to fig. 1, an embodiment of the present application provides a method for correcting a residual image of a flat panel detector, which may be used in a DR imaging system, and the method may include the following steps:
s101, obtaining a current frame DR image obtained by scanning a detected body and a plurality of groups of previous continuous DR images, wherein each group of DR images comprises a bright field image and a dark field image;
in this embodiment of the present application, the multiple sets of continuous DR images refer to images acquired before the current frame DR image, for example, the current frame DR image acquired by the present exposure is a bright field image acquired by the nth exposure (or an nth frame bright field image), and the multiple sets of continuous DR images may be a set of DR images acquired by the 1 st exposure (including the 1 st frame bright field image and the 1 st frame dark field image) to a set of DR images acquired by the nth-1 st exposure (including the n-1 st frame bright field image and the n-1 st frame dark field image).
S102, judging whether exposure corresponding to any one of the plurality of groups of continuous DR images belongs to saturated exposure or unsaturated exposure;
in this embodiment of the present application, when DR images are acquired, it may be determined, according to the received X-ray signal, whether the current exposure is saturated exposure or unsaturated exposure, and the mark is performed by using a corresponding identifier, and when performing the residual image correction, it may be determined whether the exposure corresponding to the set of DR images belongs to saturated exposure or unsaturated exposure by querying the identifier corresponding to the set of DR images.
S103, if the image is saturated exposure, determining a residual image value of the bright field image in the DR image of the current frame according to the determined saturated exposure residual image compensation coefficient, the bright field image in the DR image of the group and the determined final continuous saturated residual image attenuation curve;
s104, if the image is in unsaturated exposure, determining a residual image value of a bright field image in the DR image of the group of DR images in the current frame according to the determined unsaturated exposure residual image compensation coefficient, dark field images in the DR image of the group of DR images and the determined final continuous unsaturated residual image attenuation curve;
s105, subtracting the sum of residual image values of the bright field images in each group of DR images in the current frame DR image from the current frame DR image to obtain a DR image after residual image correction.
In the embodiment of the present application, the following equation (1) may be used to calculate the residual image value of the previously exposed image in the DR image of the current frame, but is not limited to the following equation (1).
Wherein Ghost (i) is the residual image value of the i-th frame Bright field image remained in the current frame DR image, kb is an unsaturated exposure compensation coefficient, kc is a saturated exposure compensation coefficient, dark (i) is the image mean value of the i-th frame Dark field image, bright (i) is the image mean value of the i-th frame Bright field image, lag1 (t) is a continuous unsaturated residual image attenuation curve, lag2 (t) is a continuous saturated residual image attenuation curve, t ref For the reference frame acquisition time, i.e. the time interval from the ith bright field image to the ith dark field image, t i For the time interval from the ith Bright field image to the current frame DR image, T is a saturated gray threshold, bright (i) < T indicates that the ith exposure is unsaturated exposure, and right (i) > T indicates that the ith exposure is saturated exposure. Therefore, if the i-th exposure is an unsaturated exposure, the ghost value can be calculated using the formula of the upper part in formula (1), and if the i-th exposure is a saturated exposure, the ghost value can be calculated using the formula of the lower part in formula (1).
In the embodiment of the present application, after determining the residual image value of the bright field image in each set of DR images that remains in the DR image of the current frame, the DR image after the residual image correction may be obtained using the following formula (2).
Wherein CaliImg is a DR image after the ghost correction, colImg is a DR image of a current frame (or DR image before correction), N is a bright-field image of N frames before the DR image of the current frame, and there is a ghost value remaining in the DR image of the current frame.
In this embodiment of the present application, the continuous saturated/unsaturated ghost attenuation curve and the saturated/unsaturated exposure ghost compensation coefficient may be determined before delivery, or may be determined and stored after user testing, which is not limited in this application.
In some embodiments, the determining of the final continuous saturated ghost attenuation curve includes the following steps:
acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load saturated exposure, and collecting one dark field image every set time (delta t) in a set time period (for example, 90 s) after exposure;
generating a discrete saturated residual shadow attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load saturated exposure and each dark field image after exposure;
fitting the discrete saturated ghost attenuation curve to obtain the final continuous saturated ghost attenuation curve.
In this embodiment of the present application, from the dark field image before exposure, the bright field image in the case of no-load saturated exposure, and each of the dark field images after exposure, a discrete saturated ghost attenuation curve Lag (n) may be generated by the following formula (3):
Wherein, dark (0) is the image mean of Dark field image before exposure, dark (N) is the image mean of Dark field image collected N time after exposure, bright is the image mean of Bright field image under no-load saturated exposure, n=1, 2, …, N.
In some embodiments, the discrete saturated ghost attenuation curve may be fitted to obtain an original continuous saturated ghost attenuation curve, and then used as a final continuous saturated ghost attenuation curve.
In this embodiment, according to the geometric rule of the ghost attenuation curve, the discrete saturated ghost attenuation curve may be fitted by using a double-exponential model to obtain a final continuous saturated ghost attenuation curve, where the final continuous saturated ghost attenuation curve may be represented by, for example, lang (t) =a×exp (b×t) +c×exp (d×t), and the parameters a, b, c and d are determined according to the ghost values in the discrete saturated ghost attenuation curve.
In other embodiments, to improve the accuracy of the ghost correction, the fitting the discrete saturated ghost attenuation curve to obtain the final continuous saturated ghost attenuation curve includes:
fitting the discrete saturated ghost attenuation curve to obtain an original continuous saturated ghost attenuation curve;
Determining a saturation attenuation time threshold corresponding to the original continuous saturated ghost attenuation curve according to the original continuous saturated ghost attenuation curve and a set threshold;
and correcting the original continuous saturated ghost attenuation curve according to the saturated attenuation time threshold value to obtain the final continuous saturated ghost attenuation curve.
In this embodiment of the present application, according to the original continuous saturated ghost attenuation curve and the set threshold, a saturated attenuation time threshold corresponding to the original continuous saturated ghost attenuation curve is determined, for example, when after a certain point on the original continuous saturated ghost attenuation curve passes, the difference between the ghost values of any two adjacent points is smaller than the set threshold, and then the point can be determined as the saturated attenuation time threshold corresponding to the original continuous saturated ghost attenuation curve.
Generally, after a certain period of time after exposure, the image sticking value of the flat panel detector is basically stabilized at a constant value, and the image sticking attenuation curve fitted on the flat panel detector is a function of attenuation until the image sticking attenuation curve is close to 0, if the image sticking value is calculated according to the image sticking attenuation curve continuously when the image sticking attenuation time is greater than the saturation attenuation time threshold, the image sticking is not removed completely, so that in order to improve the accuracy of image sticking correction, the original continuous saturated image sticking attenuation curve can be corrected.
In some embodiments, the correcting the original continuous saturated ghost attenuation curve according to the saturation attenuation time threshold includes:
when the ghost attenuation time is smaller than or equal to the saturation attenuation time threshold, maintaining the original continuous saturation ghost attenuation curve;
and when the ghost decay time is larger than the saturation decay time threshold, replacing the ghost values of points larger than the decay time threshold on the original continuous saturated ghost decay curve with the ghost values corresponding to the saturation decay time threshold.
In this embodiment, the original continuous saturated ghost attenuation curve may be corrected by the following formula (4):
wherein t is c Is the saturation decay time threshold.
In this embodiment of the present application, the method for determining the final continuous unsaturated ghost attenuation curve is similar to the method for determining the final continuous saturated ghost attenuation curve, and the final continuous unsaturated ghost attenuation curve can be determined by referring to the method for determining the final continuous saturated ghost attenuation curve.
In some embodiments, the determining of the final continuous unsaturated ghost attenuation curve includes the following steps:
Acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load unsaturated exposure, and collecting a dark field image at intervals in a set time period after exposure;
generating a discrete unsaturated residual image attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load unsaturated exposure and each dark field image after exposure;
fitting the discrete unsaturated ghost attenuation curve to obtain the final continuous unsaturated ghost attenuation curve.
In some embodiments, to improve the accuracy of the ghost correction, the fitting the discrete unsaturated ghost attenuation curve to obtain the final continuous unsaturated ghost attenuation curve includes:
fitting the discrete unsaturated ghost attenuation curve to obtain an original continuous unsaturated ghost attenuation curve;
determining an unsaturated attenuation time threshold corresponding to the original continuous unsaturated ghost attenuation curve according to the original continuous unsaturated ghost attenuation curve and a set threshold;
and correcting the original continuous unsaturated ghost attenuation curve according to the unsaturated attenuation time threshold to obtain the final continuous unsaturated ghost attenuation curve.
In some embodiments, the correcting the original continuous unsaturated ghost attenuation curve according to the unsaturated attenuation time threshold includes:
when the ghost decay time is less than or equal to the unsaturated decay time threshold, maintaining the original continuous unsaturated ghost decay curve;
and when the ghost decay time is greater than the unsaturated decay time threshold, replacing the ghost values of points, which are greater than the decay time threshold, on the original continuous unsaturated ghost decay curve with the ghost values corresponding to the unsaturated decay time threshold.
In some embodiments, the determining the saturated exposure residual compensation coefficient includes the following steps:
collecting a first bright field image under the condition of no-load saturated exposure;
collecting a bright field image under the condition of loaded saturated exposure;
collecting a second bright field image under the condition of no-load saturated exposure again;
determining a ghost value in the second bright-field image according to a pending saturated exposure ghost compensation coefficient, the bright-field image under the condition of loaded saturated exposure and the continuous saturated ghost attenuation curve;
and subtracting the residual image value in the second bright field image from the second bright field image to approximate the first bright field image, and reversely calculating the saturated exposure residual image compensation coefficient (Kc).
In some embodiments, the determining the unsaturated exposure residual compensation coefficient includes the following steps:
collecting a third bright field image under the condition of no-load unsaturated exposure;
collecting a bright field image under the condition of loaded unsaturated exposure and a dark field image after exposure;
acquiring a fourth bright field image under the condition of no-load unsaturated exposure again;
determining a ghost value in the fourth bright-field image according to a to-be-determined unsaturated exposure ghost compensation coefficient, the dark-field image subjected to loaded unsaturated exposure and the continuous unsaturated ghost attenuation curve;
and subtracting the residual image value in the fourth bright field image from the fourth bright field image to approximate the third bright field image, and reversely calculating the unsaturated exposure residual image compensation coefficient (Kb).
It should be noted that, before processing, the collected image in the present application generally performs background correction, gain correction and dead pixel correction on the collected image, and then performs subsequent processing on the image after the background correction, gain correction and dead pixel correction.
Based on the same inventive concept, as shown in fig. 2, an embodiment of the present application further provides a device for correcting a residual image of a flat panel detector, where the device includes: the image acquisition module 11, the judgment module 12, the first ghost value determination module 13, the second ghost value determination module 14 and the ghost correction module 15.
An image acquisition module 11 configured to acquire a current frame DR image scanned from a subject and a plurality of consecutive DR images before, each of the DR images including a bright-field image and a dark-field image;
a judging module 12 configured to judge, for any one of the plurality of sets of continuous DR images, whether exposure corresponding to the set of DR images belongs to saturated exposure or unsaturated exposure;
a first ghost value determining module 13 configured to determine, if the first image is saturated, a ghost value of a bright field image in the set of DR images remaining in the current frame DR image according to the determined saturation exposure ghost compensation coefficient, the bright field image in the set of DR images, and the determined final continuous saturated ghost attenuation curve;
a second ghost value determining module 14 configured to determine, if the image is an unsaturated exposure, a ghost value of a bright field image in the set of DR images remaining in the current frame DR image according to the determined unsaturated exposure ghost compensation coefficient, the dark field images in the set of DR images, and the determined final continuous unsaturated ghost decay curve;
the afterimage correction module 15 is configured to subtract the sum of afterimage values of bright field images in each group of DR images remaining in the current frame DR image from the current frame DR image to obtain an afterimage corrected DR image.
In some embodiments, as shown in fig. 3, the apparatus for correcting a residual image of a flat panel detector further includes: the ghost attenuation curve determination module 16.
The ghost attenuation curve determination module 16 is configured to:
acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load saturated exposure/unsaturated exposure, and collecting a dark field image at each set time interval in a set time period after exposure;
generating a discrete saturated/unsaturated residual image attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load saturated exposure/unsaturated exposure and each dark field image after exposure;
fitting the discrete saturated/unsaturated ghost attenuation curve to obtain the final continuous saturated/unsaturated ghost attenuation curve.
In some embodiments, the ghost decay curve determination module 16 is configured to:
fitting the discrete saturated/unsaturated ghost attenuation curve to obtain an original continuous saturated/unsaturated ghost attenuation curve;
determining a saturated/unsaturated decay time threshold corresponding to the original continuous saturated/unsaturated ghost decay curve according to the original continuous saturated/unsaturated ghost decay curve and a set threshold;
And correcting the original continuous saturated/unsaturated ghost attenuation curve according to the saturated/unsaturated attenuation time threshold to obtain the final continuous saturated/unsaturated ghost attenuation curve.
In some embodiments, the ghost decay curve determination module 16 is configured to: maintaining the original continuous saturated/unsaturated ghost decay curve when the ghost decay time is less than or equal to the saturated/unsaturated decay time threshold;
and when the ghost decay time is greater than the saturated/unsaturated decay time threshold, replacing the ghost values of points on the original continuous saturated/unsaturated ghost decay curve which are greater than the decay time threshold with the ghost values corresponding to the saturated/unsaturated decay time threshold.
In some embodiments, as shown in fig. 4, the apparatus for correcting a residual image of a flat panel detector further includes: the first compensation coefficient determination module 17.
The first compensation coefficient determination module 17 is configured to:
collecting a first bright field image under the condition of no-load saturated exposure;
collecting a bright field image under the condition of loaded saturated exposure;
collecting a second bright field image under the condition of no-load saturated exposure again;
determining a ghost value in the second bright-field image according to a pending saturated exposure ghost compensation coefficient, the bright-field image under the condition of loaded saturated exposure and the continuous saturated ghost attenuation curve;
And subtracting the residual image value in the second bright field image from the second bright field image to approximate the first bright field image, and reversely calculating the saturated exposure residual image compensation coefficient.
In some embodiments, as shown in fig. 4, the apparatus for correcting a residual image of a flat panel detector further includes: the second compensation coefficient determination module 18.
The second compensation coefficient determination module 18 is configured to:
collecting a third bright field image under the condition of no-load unsaturated exposure;
collecting a bright field image under the condition of loaded unsaturated exposure and a dark field image after exposure;
acquiring a fourth bright field image under the condition of no-load unsaturated exposure again;
determining a ghost value in the fourth bright-field image according to a to-be-determined unsaturated exposure ghost compensation coefficient, the dark-field image subjected to loaded unsaturated exposure and the continuous unsaturated ghost attenuation curve;
and subtracting the residual image value in the fourth bright field image from the fourth bright field image to approximate the third bright field image, and reversely calculating the unsaturated exposure residual image compensation coefficient.
The implementation process of the functions and roles of each unit in the above device is specifically shown in the implementation process of the corresponding steps in the above method, and will not be described herein again.
For the device embodiments, reference is made to the description of the method embodiments for the relevant points, since they essentially correspond to the method embodiments. The apparatus embodiments described above are merely illustrative, wherein the elements illustrated as separate elements may or may not be physically separate, and the elements shown as elements may or may not be physical elements, may be located in one place, or may be distributed over a plurality of network elements. Some or all of the modules may be selected according to actual needs to achieve the purposes of the present application. Those of ordinary skill in the art will understand and implement the present invention without undue burden.
Based on the same inventive concept, the embodiments of the present application also provide a storage medium having stored thereon a computer program, which when executed by a processor, implements the steps of the flat panel detector ghost correction method in any of the possible implementations described above.
Alternatively, the storage medium may be a non-transitory computer readable storage medium, which may be, for example, ROM, random Access Memory (RAM), CD-ROM, magnetic tape, floppy disk, optical data storage device, and the like.
Based on the same inventive concept, referring to fig. 5, the embodiment of the present application further provides a medical device, including a memory 71 (e.g. a non-volatile memory), a processor 72, and a computer program stored on the memory 71 and executable on the processor 72, wherein the steps of the flat panel detector ghost correction method in any of the possible implementations described above are implemented when the processor 72 executes the program. The medical device may be, for example, a DR imaging system.
As shown in fig. 5, the medical device may generally further include: memory 73, network interface 74, and internal bus 75. In addition to these components, other hardware may be included, which is not described in detail.
It should be noted that the above-mentioned device for correcting the residual image of the flat panel detector may be implemented by software, and is a device in a logic sense, and is formed by reading, by the processor 72 of the medical device in which the device is located, the computer program instructions stored in the nonvolatile memory into the memory 73 for operation.
Embodiments of the subject matter and the functional operations described in this specification can be implemented in: digital electronic circuitry, tangibly embodied computer software or firmware, computer hardware including the structures disclosed in this specification and structural equivalents thereof, or a combination of one or more of them. Embodiments of the subject matter described in this specification can be implemented as one or more computer programs, i.e., one or more modules of computer program instructions encoded on a tangible, non-transitory program carrier for execution by, or to control the operation of, data processing apparatus. Alternatively or additionally, the program instructions may be encoded on a manually-generated propagated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal, that is generated to encode and transmit information to suitable receiver apparatus for execution by data processing apparatus. The computer storage medium may be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or a combination of one or more of them.
The processes and logic flows described in this specification can be performed by one or more programmable computers executing one or more computer programs to perform corresponding functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application-specific integrated circuit).
Computers suitable for executing computer programs include, for example, general purpose and/or special purpose microprocessors, or any other type of central processing unit. Typically, the central processing unit will receive instructions and data from a read only memory and/or a random access memory. The essential elements of a computer include a central processing unit for carrying out or executing instructions and one or more memory devices for storing instructions and data. Typically, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto-optical disks, or optical disks, etc. However, a computer does not have to have such a device. Furthermore, the computer may be embedded in another device, such as a mobile phone, a Personal Digital Assistant (PDA), a mobile audio or video player, a game console, a Global Positioning System (GPS) receiver, or a portable storage device such as a Universal Serial Bus (USB) flash drive, to name a few.
Computer readable media suitable for storing computer program instructions and data include all forms of non-volatile memory, media and memory devices including, for example, semiconductor memory devices (e.g., EPROM, EEPROM, and flash memory devices), magnetic disks (e.g., internal hard disk or removable disks), magneto-optical disks, and CD-ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.
While this specification contains many specific implementation details, these should not be construed as limitations on the scope of any invention or of what may be claimed, but rather as descriptions of features of specific embodiments of particular inventions. Certain features that are described in this specification in the context of separate embodiments can also be implemented in combination in a single embodiment. On the other hand, the various features described in the individual embodiments may also be implemented separately in the various embodiments or in any suitable subcombination. Furthermore, although features may be acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a subcombination or variation of a subcombination.
Similarly, although operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results. In some cases, multitasking and parallel processing may be advantageous. Moreover, the separation of various system modules and components in the embodiments described above should not be understood as requiring such separation in all embodiments, and it should be understood that the described program components and systems can generally be integrated together in a single software product or packaged into multiple software products.
Thus, particular embodiments of the subject matter have been described. Other embodiments are within the scope of the following claims. In some cases, the actions recited in the claims can be performed in a different order and still achieve desirable results. Furthermore, the processes depicted in the accompanying drawings are not necessarily required to be in the particular order shown, or sequential order, to achieve desirable results. In some implementations, multitasking and parallel processing may be advantageous.
The foregoing description of the preferred embodiments of the present invention is not intended to limit the invention to the precise form disclosed, and any modifications, equivalents, improvements and alternatives falling within the spirit and principles of the present invention are intended to be included within the scope of the present invention.
Claims (8)
1. A method for correcting a residual shadow of a flat panel detector, the method comprising:
obtaining a current frame DR image obtained by scanning a detected body and a plurality of groups of previous continuous DR images, wherein each group of DR images comprises a bright field image and a dark field image;
judging whether exposure corresponding to the DR images belongs to saturated exposure or unsaturated exposure for any one of the continuous DR images; the unsaturated exposure means that the incident dose of X-rays does not exceed a saturation threshold value during exposure, and the image mean value of the DR image does not exceed a saturation gray level threshold value; the saturation exposure means that the incident dose of X-rays exceeds a saturation threshold value during exposure, and the image mean value of the DR image exceeds a saturation gray level threshold value;
if the image is saturated exposure, determining a residual image value of the bright field image in the DR image of the current frame according to the determined saturated exposure residual image compensation coefficient, the bright field image in the DR image of the group and the determined final continuous saturated residual image attenuation curve;
if the image is unsaturated exposure, determining a residual image value of a bright field image in the DR image of the group of DR images in the current frame according to the determined unsaturated exposure residual image compensation coefficient, a dark field image in the DR image of the group of DR images and the determined final continuous unsaturated residual image attenuation curve;
Subtracting the sum of residual image values of bright field images in each group of DR images in the current frame DR image from the current frame DR image to obtain a DR image after residual image correction;
the determination of the final continuous saturated ghost attenuation curve comprises the following steps:
acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load saturated exposure, and collecting a dark field image at each set time interval in a set time period after exposure;
generating a discrete saturated residual shadow attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load saturated exposure and each dark field image after exposure;
fitting the discrete saturated ghost attenuation curve to obtain the final continuous saturated ghost attenuation curve; the determination of the final continuous unsaturated ghost attenuation curve comprises the following steps:
acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load unsaturated exposure, and collecting a dark field image at intervals in a set time period after exposure;
generating a discrete unsaturated residual image attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load unsaturated exposure and each dark field image after exposure;
Fitting the discrete unsaturated ghost attenuation curve to obtain the final continuous unsaturated ghost attenuation curve;
the determination of the saturated exposure residual shadow compensation coefficient comprises the following steps:
collecting a first bright field image under the condition of no-load saturated exposure;
collecting a bright field image under the condition of loaded saturated exposure;
collecting a second bright field image under the condition of no-load saturated exposure again;
determining a ghost value in the second bright-field image according to a pending saturated exposure ghost compensation coefficient, the bright-field image under the condition of loaded saturated exposure and the continuous saturated ghost attenuation curve;
the difference value obtained by subtracting the first bright field image from the second bright field image is compared with the residual image value in the second bright field image, and the saturated exposure residual image compensation coefficient is calculated;
the determination of the unsaturated exposure residual image compensation coefficient comprises the following steps:
collecting a third bright field image under the condition of no-load unsaturated exposure;
collecting a bright field image under the condition of loaded unsaturated exposure and a dark field image after exposure;
acquiring a fourth bright field image under the condition of no-load unsaturated exposure again;
Determining a ghost value in the fourth bright-field image according to a to-be-determined unsaturated exposure ghost compensation coefficient, the dark-field image subjected to loaded unsaturated exposure and the continuous unsaturated ghost attenuation curve;
and comparing the difference value obtained by subtracting the third bright field image from the fourth bright field image with the residual image value in the fourth bright field image, and calculating the unsaturated exposure residual image compensation coefficient.
2. The method of claim 1, wherein said fitting the discrete saturated ghost decay curve to obtain the final continuous saturated ghost decay curve comprises:
fitting the discrete saturated ghost attenuation curve to obtain an original continuous saturated ghost attenuation curve;
determining a saturation attenuation time threshold corresponding to the original continuous saturated ghost attenuation curve according to the original continuous saturated ghost attenuation curve and a set threshold;
correcting the original continuous saturated ghost attenuation curve according to the saturated attenuation time threshold to obtain the final continuous saturated ghost attenuation curve;
fitting the discrete unsaturated ghost attenuation curve to obtain the final continuous unsaturated ghost attenuation curve, wherein the fitting comprises the following steps:
Fitting the discrete unsaturated ghost attenuation curve to obtain an original continuous unsaturated ghost attenuation curve;
determining an unsaturated attenuation time threshold corresponding to the original continuous unsaturated ghost attenuation curve according to the original continuous unsaturated ghost attenuation curve and a set threshold;
and correcting the original continuous unsaturated ghost attenuation curve according to the unsaturated attenuation time threshold to obtain the final continuous unsaturated ghost attenuation curve.
3. The method of claim 2, wherein modifying the original continuous saturated ghost decay curve based on the saturation decay time threshold comprises:
when the ghost attenuation time is smaller than or equal to the saturation attenuation time threshold, maintaining the original continuous saturation ghost attenuation curve;
when the ghost attenuation time is greater than the saturation attenuation time threshold, replacing the ghost values of points on the original continuous saturated ghost attenuation curve which are greater than the attenuation time threshold with the ghost values corresponding to the saturation attenuation time threshold;
the correcting the original continuous unsaturated ghost attenuation curve according to the unsaturated attenuation time threshold includes:
When the ghost decay time is less than or equal to the unsaturated decay time threshold, maintaining the original continuous unsaturated ghost decay curve;
and when the ghost decay time is greater than the unsaturated decay time threshold, replacing the ghost values of points, which are greater than the decay time threshold, on the original continuous unsaturated ghost decay curve with the ghost values corresponding to the unsaturated decay time threshold.
4. A flat panel detector ghost correction device, the device comprising:
an image acquisition module configured to acquire a current frame DR image obtained by scanning a subject and a plurality of previous continuous DR images, each of the DR images including a bright-field image and a dark-field image;
a judging module configured to judge, for any one of the plurality of sets of continuous DR images, whether exposure corresponding to the set of DR images belongs to saturated exposure or unsaturated exposure; the unsaturated exposure means that the incident dose of X-rays does not exceed a saturation threshold value during exposure, and the image mean value of the DR image does not exceed a saturation gray level threshold value; the saturation exposure means that the incident dose of X-rays exceeds a saturation threshold value during exposure, and the image mean value of the DR image exceeds a saturation gray level threshold value;
A first ghost value determining module configured to determine, if the first image is saturated, a ghost value of a bright field image in the set of DR images remaining in the current frame DR image according to the determined saturated exposure ghost compensation coefficient, the bright field image in the set of DR images, and the determined final continuous saturated ghost attenuation curve;
a second ghost value determining module configured to determine, if the image is unsaturated, a ghost value of a bright field image in the set of DR images remaining in the current frame DR image according to the determined unsaturated exposure ghost compensation coefficient, a dark field image in the set of DR images, and the determined final continuous unsaturated ghost attenuation curve;
a residual image correction module configured to subtract a sum of residual image values of bright field images in each group of DR images from the current frame DR image to obtain a DR image after residual image correction;
the ghost attenuation curve determining module;
if the exposure is saturated, the ghost attenuation curve determining module is configured to:
acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load saturated exposure, and collecting a dark field image at each set time interval in a set time period after exposure;
Generating a discrete saturated residual shadow attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load saturated exposure and each dark field image after exposure;
fitting the discrete saturated ghost attenuation curve to obtain the final continuous saturated ghost attenuation curve;
if the exposure is unsaturated, the ghost decay curve determination module is configured to:
acquiring a dark field image before exposure;
collecting bright field images under the condition of no-load unsaturated exposure, and collecting a dark field image at intervals in a set time period after exposure;
generating a discrete unsaturated residual image attenuation curve according to the dark field image before exposure, the bright field image under the condition of no-load unsaturated exposure and each dark field image after exposure;
fitting the discrete unsaturated ghost attenuation curve to obtain the final continuous unsaturated ghost attenuation curve;
a first compensation coefficient determination module;
the first compensation coefficient determination module is configured to:
collecting a first bright field image under the condition of no-load saturated exposure;
collecting a bright field image under the condition of loaded saturated exposure;
Collecting a second bright field image under the condition of no-load saturated exposure again;
determining a ghost value in the second bright-field image according to a pending saturated exposure ghost compensation coefficient, the bright-field image under the condition of loaded saturated exposure and the continuous saturated ghost attenuation curve;
the difference value obtained by subtracting the first bright field image from the second bright field image is compared with the residual image value in the second bright field image, and the saturated exposure residual image compensation coefficient is calculated;
a second compensation coefficient determination module;
the second compensation coefficient determination module is configured to:
collecting a third bright field image under the condition of no-load unsaturated exposure;
collecting a bright field image under the condition of loaded unsaturated exposure and a dark field image after exposure;
acquiring a fourth bright field image under the condition of no-load unsaturated exposure again;
determining a ghost value in the fourth bright-field image according to a to-be-determined unsaturated exposure ghost compensation coefficient, the dark-field image subjected to loaded unsaturated exposure and the continuous unsaturated ghost attenuation curve;
and comparing the difference value obtained by subtracting the third bright field image from the fourth bright field image with the residual image value in the fourth bright field image, and calculating the unsaturated exposure residual image compensation coefficient.
5. The apparatus of claim 4, wherein the device comprises a plurality of sensors,
if the exposure is saturated, the ghost attenuation curve determining module is configured to:
fitting the discrete saturated ghost attenuation curve to obtain an original continuous saturated ghost attenuation curve;
determining a saturation attenuation time threshold corresponding to the original continuous saturated ghost attenuation curve according to the original continuous saturated ghost attenuation curve and a set threshold;
correcting the original continuous saturated ghost attenuation curve according to the saturated attenuation time threshold to obtain the final continuous saturated ghost attenuation curve;
if the exposure is unsaturated, the ghost decay curve determination module is configured to:
fitting the discrete unsaturated ghost attenuation curve to obtain an original continuous unsaturated ghost attenuation curve;
determining an unsaturated attenuation time threshold corresponding to the original continuous unsaturated ghost attenuation curve according to the original continuous unsaturated ghost attenuation curve and a set threshold;
and correcting the original continuous unsaturated ghost attenuation curve according to the unsaturated attenuation time threshold to obtain the final continuous unsaturated ghost attenuation curve.
6. The apparatus of claim 5, wherein the device comprises a plurality of sensors,
if the exposure is saturated, the ghost attenuation curve determining module is configured to:
when the ghost attenuation time is smaller than or equal to the saturation attenuation time threshold, maintaining the original continuous saturation ghost attenuation curve;
when the ghost attenuation time is greater than the saturation attenuation time threshold, replacing the ghost values of points on the original continuous saturated ghost attenuation curve which are greater than the attenuation time threshold with the ghost values corresponding to the saturation attenuation time threshold;
if the exposure is unsaturated, the ghost decay curve determination module is configured to:
when the ghost decay time is less than or equal to the unsaturated decay time threshold, maintaining the original continuous unsaturated ghost decay curve;
and when the ghost decay time is greater than the unsaturated decay time threshold, replacing the ghost values of points, which are greater than the decay time threshold, on the original continuous unsaturated ghost decay curve with the ghost values corresponding to the unsaturated decay time threshold.
7. A storage medium having stored thereon a computer program, which when executed by a processor performs the steps of the method of any of claims 1-3.
8. A medical device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of the method of any of claims 1-3 when the program is executed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010158713.6A CN111445397B (en) | 2020-03-09 | 2020-03-09 | Flat panel detector ghost correction method and device, storage medium and medical equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010158713.6A CN111445397B (en) | 2020-03-09 | 2020-03-09 | Flat panel detector ghost correction method and device, storage medium and medical equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111445397A CN111445397A (en) | 2020-07-24 |
| CN111445397B true CN111445397B (en) | 2023-07-04 |
Family
ID=71655770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010158713.6A Active CN111445397B (en) | 2020-03-09 | 2020-03-09 | Flat panel detector ghost correction method and device, storage medium and medical equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111445397B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118154699A (en) * | 2024-04-30 | 2024-06-07 | 合肥埃科光电科技股份有限公司 | Calibration method, image afterglow correction method, device and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003185752A (en) * | 2001-12-17 | 2003-07-03 | Shimadzu Corp | Method and device for correcting image of radiation detector, and radiation imaging device using the same |
| CN106097282A (en) * | 2016-07-01 | 2016-11-09 | 上海奕瑞光电子科技有限公司 | Eliminate method and the flat panel detector of flat panel detector image ghost |
| CN106296595A (en) * | 2015-06-08 | 2017-01-04 | 上海奕瑞光电子科技有限公司 | A kind of flat panel detector and the method for reduction flat panel detector image ghost |
| CN108171765A (en) * | 2017-12-08 | 2018-06-15 | 上海奕瑞光电子科技股份有限公司 | The bearing calibration of flat panel detector ghost and its means for correcting |
| CN108172659A (en) * | 2017-12-20 | 2018-06-15 | 上海奕瑞光电子科技股份有限公司 | The generation method of flat panel detector and its ghost tables of data, ghost compensation correction method |
-
2020
- 2020-03-09 CN CN202010158713.6A patent/CN111445397B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003185752A (en) * | 2001-12-17 | 2003-07-03 | Shimadzu Corp | Method and device for correcting image of radiation detector, and radiation imaging device using the same |
| CN106296595A (en) * | 2015-06-08 | 2017-01-04 | 上海奕瑞光电子科技有限公司 | A kind of flat panel detector and the method for reduction flat panel detector image ghost |
| CN106097282A (en) * | 2016-07-01 | 2016-11-09 | 上海奕瑞光电子科技有限公司 | Eliminate method and the flat panel detector of flat panel detector image ghost |
| CN108171765A (en) * | 2017-12-08 | 2018-06-15 | 上海奕瑞光电子科技股份有限公司 | The bearing calibration of flat panel detector ghost and its means for correcting |
| CN108172659A (en) * | 2017-12-20 | 2018-06-15 | 上海奕瑞光电子科技股份有限公司 | The generation method of flat panel detector and its ghost tables of data, ghost compensation correction method |
Non-Patent Citations (1)
| Title |
|---|
| 医用诊断X射线设备技术资料审评要点探讨;王博等;《中国医疗器械信息》;20111125(第11期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111445397A (en) | 2020-07-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111667448B (en) | Image processing method, device and equipment | |
| US6920198B2 (en) | Methods and apparatus for processing a fluoroscopic image | |
| US10045746B2 (en) | Radiation image processing apparatus, method, and medium | |
| EP2293239B1 (en) | Line noise reduction apparatus and method | |
| CN106920245B (en) | Boundary detection method and device | |
| CN102096917A (en) | Automatic eliminating method for redundant image data of capsule endoscope | |
| US20080012967A1 (en) | Defective-area correction apparatus, method and program and radiation detection apparatus | |
| US8355064B2 (en) | Noise reduction for machine vision systems | |
| JP2017000565A (en) | Image processing apparatus, radiographic apparatus, image processing method, program, and storage medium | |
| JP5561112B2 (en) | Image processing apparatus, imaging apparatus, and program | |
| CN111445397B (en) | Flat panel detector ghost correction method and device, storage medium and medical equipment | |
| WO2007110798A1 (en) | Temperature artifact correction | |
| CN105141853A (en) | Image processing method and electronic device | |
| US20140086390A1 (en) | Radiation exposure dose obtaining method and apparatus, and radiation image capturing system | |
| JP6491545B2 (en) | Image processing apparatus, radiation imaging apparatus, image processing method, program, and storage medium | |
| CN113034389A (en) | Image processing method, image processing device, computer equipment and storage medium | |
| CN109903245B (en) | Non-uniform correction method for infrared image | |
| WO2011119478A2 (en) | Transient pixel defect detection and correction | |
| CN115266779A (en) | Flat panel detector pixel response nonlinear correction method for CT scanning | |
| CN111612773A (en) | Thermal infrared imager and real-time automatic blind pixel detection processing method | |
| US8090215B2 (en) | Method and system for noise level detection in image data | |
| CN110428375B (en) | DR image processing method and device | |
| US20130039561A1 (en) | Image processing apparatus, image processing method and storage medium | |
| CN117314771A (en) | Belt artifact processing method and device, electronic equipment and storage medium | |
| CN108171765B (en) | Method and device for correcting residual image of flat panel detector |
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 |