WO2020159201A1 - Radiographic image processing method and radiography apparatus - Google Patents
Radiographic image processing method and radiography apparatus Download PDFInfo
- Publication number
- WO2020159201A1 WO2020159201A1 PCT/KR2020/001303 KR2020001303W WO2020159201A1 WO 2020159201 A1 WO2020159201 A1 WO 2020159201A1 KR 2020001303 W KR2020001303 W KR 2020001303W WO 2020159201 A1 WO2020159201 A1 WO 2020159201A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- pixel
- value
- target pixel
- pixel value
- abnormal
- Prior art date
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 20
- 238000002601 radiography Methods 0.000 title claims abstract 3
- 230000002159 abnormal effect Effects 0.000 claims abstract description 136
- 230000005855 radiation Effects 0.000 claims description 92
- 238000000034 method Methods 0.000 claims description 45
- 238000001514 detection method Methods 0.000 claims description 13
- 238000012795 verification Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- 230000005856 abnormality Effects 0.000 claims description 6
- 230000002093 peripheral effect Effects 0.000 description 23
- 238000003384 imaging method Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 8
- 229910021420 polycrystalline silicon Inorganic materials 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 2
- 238000007408 cone-beam computed tomography Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 229910044991 metal oxide Inorganic materials 0.000 description 2
- 150000004706 metal oxides Chemical class 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000019771 cognition Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- 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
-
- 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/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
-
- 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
Definitions
- the present invention relates to a radiographic image processing method and a radiographic imaging apparatus, and more particularly, to a radiographic image processing method and a radiographic imaging apparatus for detecting abnormal pixel values from a plurality of radiographic images obtained over time.
- a radiographic imaging apparatus includes a radiation generator that irradiates radiation to a subject and a radiation detector that detects radiation that has transmitted through the subject with a radiation detection panel.
- the radiation detection panel may be made of a material such as a complementary metal-oxide semiconductor (CMOS) or a crystal such as polycrystalline silicon (Poly-Si), and noise may occur according to the characteristics of the material.
- CMOS complementary metal-oxide semiconductor
- Poly-Si polycrystalline silicon
- the noise has a relatively high value
- the radiographic image may be distorted, thereby degrading the quality of the radiographic image.
- Patent Document 1 Korean Registered Patent Publication No. 10-0617517
- the present invention provides a radiographic image processing method and a radiographic imaging apparatus for detecting abnormal pixel values from a plurality of radiographic images acquired over time.
- a radiographic image processing method includes: (a) obtaining a plurality of radiographic images over time using a radiation detection panel having a plurality of pixels; (b) selecting a target pixel among the plurality of pixels from two or more selected radiation images among the plurality of radiation images; (c) calculating a pixel value change amount over time of the selected target pixel; (d) comparing the calculated amount of change of the pixel value with time of the target pixel with a preset reference value; And (e) determining a pixel value of the target pixel as an abnormal value in at least one of the selected two or more radiographic images when the amount of change in the pixel value over time of the target pixel is equal to or greater than the preset reference value. It can contain.
- the step of verifying the pixel value of the target pixel determined to be an abnormal value further comprising the step of verifying the pixel value of the target pixel, the neighboring pixels adjacent to the target pixel calculated from the selected two or more radiographic images Comparing the pixel value change with time and the preset reference value; And checking the pixel value of the target pixel determined as the abnormal value as an abnormal value when the number of surrounding pixels having a change in the pixel value over time is equal to or less than the preset determination number.
- the step of verifying the pixel value of the target pixel determined to be an abnormal value further comprising the step of verifying the pixel value of the target pixel, the neighboring pixels adjacent to the target pixel calculated from the selected two or more radiographic images Comparing a pixel value change amount with time and a pixel value change amount with time of the target pixel; And the pixel value of the target pixel determined as the abnormal value is abnormal when the number of pixel values variation over time and the number of surrounding pixels having a pixel value variation within a preset allowable deviation range is equal to or less than a preset determination number. It may include the process of confirming by value.
- the process of (b) to (e) is repeated for the whole of the plurality of radiographic images, and a process of storing the position of the target pixel in which the pixel value determined as the abnormal value occurs is further included.
- the pixel storing the position may be excluded from the selection of the target pixel, and the pixel value of the pixel storing the position may be determined as an abnormal value.
- the preset reference value may be generated even when radiation is not irradiated from each pixel as a whole to the plurality of pixels, or may be greater than the inherent noise value generated by the radiation source.
- the (b) process may be performed by selecting two or more radiographic images that are temporally continuous among the plurality of radiographic images.
- the step (b) two radiographic images that are temporally continuous among the plurality of radiographic images are selected, and in the step (e), two radiographic images in which the amount of change in a pixel value over time of the target pixel is greater than or equal to the preset reference value Among the pixels, the pixel value of the target pixel may be determined as an abnormal value in a radiographic image having a high pixel value of the target pixel.
- the target pixel may be selected from the pixel regions.
- the process may further include correcting a pixel value of a target pixel determined as an abnormal value by interpolating a pixel value of adjacent neighboring pixels on the same radiation image.
- a radiographic imaging apparatus includes a plurality of pixels, and a radiation detection panel that generates electrical signals for acquiring a plurality of radiographic images over time by irradiated radiation;
- a change amount calculating unit for calculating a pixel value change amount over time of a target pixel selected from the plurality of pixels from two or more radiation images selected from the plurality of radiation images and a pixel value change amount over time of neighboring pixels adjacent to the target pixel;
- a change amount comparison unit comparing the calculated change amount of the pixel value with time of the target pixel with a preset reference value;
- a pixel value abnormality determination unit determining a pixel value of the target pixel as an abnormal value in at least one radiation image among the selected two or more radiographic images when the amount of change of the pixel value over time of the target pixel is equal to or greater than the preset reference value. It may include.
- a first pixel value verification unit configured to determine a pixel value of the target pixel determined as the abnormal value as an abnormal value when the number of surrounding pixels having a change amount of the pixel value over time is equal to or less than a preset determination number; It may further include.
- the pixel value of the target pixel determined as the abnormal value is an abnormal value when the number of pixel values variation over time and the number of surrounding pixels having a pixel value variation within a preset allowable deviation range is equal to or less than a preset determination number. It may further include; a second pixel value verification unit determined by.
- It may further include a location storage unit for storing the location of the target pixel where the pixel value determined as the abnormal value occurs.
- It may further include a reference value setting unit for setting the preset reference value that is greater than the intrinsic noise value generated even when radiation is irradiated from each pixel as a whole, or is generated by the radiation source.
- an input signal over time may be reflected to determine whether a pixel value of each target pixel is an abnormal pixel value including random noise. Accordingly, by correcting by removing abnormal pixel values including random noise from a plurality of radiographic images, distortion of the radiographic image can be prevented, and quality of the radiographic image can be improved.
- the probability of random noise is high in a pixel where random noise has occurred once, so the location of the target pixel where the abnormal pixel value is generated is stored and the pixel value of the pixel is regarded as an abnormal pixel value in the subsequent radiation image.
- the number of pixels that need to be determined to be abnormal pixel values can be reduced, and the time for image processing to determine and correct whether there are abnormal pixel values for the entire plurality of radiographic images can be shortened.
- FIG. 1 is a flowchart illustrating a radiographic image processing method according to an embodiment of the present invention.
- FIG. 2 is a conceptual diagram illustrating a target pixel and surrounding pixels according to an embodiment of the present invention.
- FIG. 3 is a conceptual diagram illustrating verification of a pixel value of a target pixel through a preset reference value according to an embodiment of the present invention.
- FIG. 4 is a conceptual diagram illustrating verification of a pixel value of a target pixel through comparison of a pixel value change amount with time of a target pixel and a change amount of pixel value with time of a surrounding pixel according to an embodiment of the present invention.
- FIG. 5 is a block diagram showing the configuration of a radiographic apparatus according to an embodiment of the present invention.
- FIG. 1 is a flowchart illustrating a radiographic image processing method according to an embodiment of the present invention.
- a radiographic image processing method includes (a) obtaining a plurality of radiographic images over time using a radiation detection panel 110 having a plurality of pixels 111.
- Process S100; (b) selecting a target pixel 111a among the plurality of pixels 111 from two or more selected radiation images 10 and 20 among the plurality of radiation images (S200); (c) calculating the amount of change in the pixel value over time of the selected target pixel 111a (S300); (d) comparing the calculated amount of change of the pixel value over time of the target pixel 111a with a preset reference value R (S400); And (e) when the amount of change in pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R, the target pixel in at least one of the selected two or more radiographic images 10 and 20 ( And determining the pixel value of 111a) as an abnormal value (S500).
- a plurality of radiation images are acquired over time using the radiation detection panel 110 including a plurality of pixels 111 (S100).
- a radiation image can be generated, and charges generated by irradiation of radiation are stored in a plurality of pixels 111 of the radiation detection panel 110.
- the radiation image may be generated according to charges (amounts) stored in the plurality of pixels 111.
- the generated radiation image may be directly obtained, or the generated and stored radiation image may be obtained.
- CT computed tomography
- CBCT computed tomography
- a target pixel 111a is selected from the plurality of pixels 111 from two or more selected radiation images 10 and 20 among the plurality of radiation images (S200).
- the target pixel 111a may be one or a plurality of pixels.
- a pixel value change amount according to the time of the selected target pixel 111a is calculated (S300).
- the amount of change in pixel values over time of each target pixel 111a may be calculated, or the amount of change in pixel values over time of all target pixels 111a may be calculated at once. Through this, the amount of change in the pixel value over time of the entire target pixel 111a may be calculated.
- the calculated change in pixel value over time of the target pixel 111a is compared with a preset reference value R (S400).
- the calculated amount of change in pixel value over time of the target pixel 111a can be compared with a preset reference value R, and the amount of change in pixel value over time of each target pixel 111a can be compared with a preset reference value R, , It can be compared with a preset reference value R for all target pixels 111a.
- the predetermined reference value R may be used to identify an element (for example, random noise) generating.
- Elements commonly included in all the pixels 111 may be excluded by the preset reference value R, and when the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R, the target pixel ( Since a different element (or noise) is interposed more than other pixels in 111a), a pixel value of the target pixel 111a may be determined as an abnormal value in the selected two or more radiographic images 10 and 20.
- the pixel value of the target pixel 111a may be determined as an abnormal value in the selected two or more radiographic images 10 and 20 through a portion in which only the pixel value of the target pixel 111a is independently changed.
- the present invention by taking into account the amount of change in the pixel value over time of the target pixel 111a between the selected two or more radiographic images 10 and 20, the amount of change in pixel value over time of the target pixel 111a and a preset reference value (R) Through comparison, it may be determined whether the pixel value of the target pixel 111a is abnormal.
- the step (b) (S200) may be performed by selecting two or more consecutive radiographic images 10 and 20 among the multiple radiographic images. Among the plurality of radiographic images, two or more temporally continuous radiographic images 10 and 20 may be selected, and a temporal tendency may be maintained (or continued) between two or more temporally continuous radiographic images 10 and 20. Therefore, the amount of change in pixel value over time can be more effective in detecting abnormal pixel values.
- the probability of occurrence of a change in the input signal for each pixel 111 is high, and when a change in radiation intensity and/or a position of the subject gradually changes with time
- the amount of change in pixel value according to is inevitably increased.
- the amount of change in the pixel value over time mostly exceeds the preset reference value R, so that only the comparison of the amount of change in the pixel value over time with the preset reference value R of the target pixel 111a results in the pixel value of the target pixel 111a.
- the overall pixel value variation is large, so the ratio of random noise to the pixel value variation may be small, which may make it difficult to verify the pixel value of the target pixel 111a through the peripheral pixel 111b. .
- the random noise refers to noise generated randomly with respect to time (or along) for each pixel 111 independently, for example, a complementary metal-oxide semiconductor; CMOS) or polycrystalline silicon (Poly-Si), such as crystals, may be generated depending on the material characteristics of the radiation detection panel 110, and may also be referred to as random telegraph signal (RTS) noise (hereinafter RTS noise).
- RTS noise random telegraph signal
- the random noise has a relatively high value. Therefore, when a radiographic image is generated using a pixel value including random noise, the radiographic image may be distorted and the quality of the radiographic image may be deteriorated.
- RTS noise is a noise generated when one or several electrons or holes in a channel of a switching device are trapped or de-traded in an oxide, and discontinuously. It shows characteristics that repeatedly display multiple signal levels at very low frequencies. As a result, in the radiographic image, the signal of the pixel 111 leads to a defect that appears randomly at several levels.
- RTS noise is a phenomenon that occurs in almost all semiconductors, and when utilized in a sensor, RTS noise is generally removed by increasing the number of sampling and taking a median value. Although a method is used, it is impossible to perform sampling several times in an image sensor of a passive pixel sensor (PPS) type.
- PPS passive pixel sensor
- step (b), in step S200, two consecutive radiographic images 10 and 20 among the plurality of radiographic images may be selected, and in step (e), in step S500, at the time of the target pixel 111a.
- the pixel value of the target pixel 111a in the radiographic image 10 or 20 where the pixel value of the target pixel 111a is relatively high among the two radiographic images 10 and 20 in which the amount of change in the pixel value is greater than or equal to the preset reference value R Can be determined as an abnormal value.
- the pixel values of the target pixel 111a may be determined as both abnormal values in the two radiographic images 10 and 20, but any one of the two radiographic images 10 and 20 may be a radiographic image 10 or 20.
- the pixel value of the target pixel 111a may be a normal value
- the pixel value of the target pixel 111a which is the normal value among the two radiographic images 10 and 20 may be filtered to be continuously used in the radiographic image.
- the radiographic image can be imaged similarly to the optimal radiographic image generated only with charge generated by radiation without noise and/or correction. Since the random noise increases the pixel value, the pixel value including the random noise is forced to be higher, and the target pixel 111a in the radiographic image 10 or 20 where the pixel value of the target pixel 111a is relatively high.
- the pixel value of can be determined as an abnormal value.
- the preset reference value R may be generated even when radiation is not irradiated from each pixel 111 as a whole of the plurality of pixels 111 or may be greater (or more) than the inherent noise value generated by radiation.
- the preset reference value R may be set larger than the intrinsic noise value common to all the plurality of pixels 111.
- the intrinsic noise may have a different value, but refers to noise common to all of the plurality of pixels 111, and may be noise generated by a device offset. Since it occurs in common with the pixels 111 of, it is not a big problem if all the device offset values are removed from the pixel values of the plurality of pixels 111 constituting the radiation image.
- the intrinsic noise is generated without causing radiation due to a cause such as a thermal cause, poor insulation, etc., and causes dark current to accumulate charge in a plurality of pixels 111 even when radiation is not irradiated.
- photons and may include quantum noise generated by a radiation source such as a radiation generator.
- the preset reference value R is an average value of the intrinsic noise values of the plurality of pixels 111 or an intermediate value among the intrinsic noise values of the plurality of pixels 111. Alternatively, it can be set larger (or higher) than the maximum value. On the other hand, since the random noise value is much different from the intrinsic noise value, the preset reference value R may be set to about twice the average value of the intrinsic noise values of the plurality of pixels 111, or the plurality of pixels 111 It can be set to about twice the minimum, median or maximum value of the intrinsic noise value of
- the preset reference value R is set to be larger than the intrinsic noise value generated in each pixel 111 as a whole of the plurality of pixels 111, thereby thereby forming a plurality of pixels between the selected two or more radiographic images 10 and 20.
- the natural variation due to the intrinsic noise can be excluded from the variation in pixel values over time in (111), and the random noise can be detected from the pixel values of the plurality of pixels 111 constituting the radiation image more effectively. have.
- the preset reference value R may be set by using a difference value between the maximum and minimum values of the intrinsic noise that may occur in each pixel 111 and is greater than the difference between the maximum and minimum values of the intrinsic noise ( Above).
- the difference value between the maximum value and the minimum value of the intrinsic noise may be calculated to set a preset reference value R above a minimum value, a median value, a maximum value, or an average value of the calculated difference values.
- the preset reference value R may be set to be equal to or greater than a difference value between the maximum value of the intrinsic noise and the minimum value of the intrinsic noise in the entire plurality of pixels 111. In addition, it may be based on 1/2 of the difference between the maximum and minimum values of the inherent noise.
- FIG. 2 is a conceptual diagram for explaining a target pixel and a surrounding pixel according to an embodiment of the present invention
- FIG. 3 is a conceptual diagram for explaining a pixel value verification of a target pixel through a preset reference value according to an embodiment of the present invention
- 3(a) shows a case determined by an abnormal pixel value
- FIG. 3(b) shows a case determined by a normal pixel value.
- the radiographic image processing method according to the present invention may further include a process of verifying a pixel value of the target pixel 111a determined as an abnormal value (S600 ).
- the pixel value of the target pixel 111a determined to be an abnormal value may be verified (S600 ).
- S600 the pixel value of the target pixel 111a determined as the abnormal value
- the pixel value of the target pixel 111a determined as an abnormal value due to only a change in the input signal is determined as a normal value, thereby forming a radiation image.
- only the abnormal pixel value including (or attributable to) the random noise may be corrected to prevent distortion of the radiographic image, and quality of the radiographic image may be improved.
- the subject moves. It is possible to more accurately determine whether the pixel value of each target pixel 111a is an abnormal pixel value including the random noise by effectively reflecting the change of the input signal due to the like.
- the amount of change in pixel value over time of the neighboring pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20 and a preset reference value R may be compared (S610 ). Similar to the target pixel 111a, the amount of change in the pixel value over time of the adjacent pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10, 20 is also compared with a preset reference value R Can. Through this, it is possible to grasp the tendency of the amount of change in the pixel value over time in the peripheral area of the target pixel 111a to be similar to the amount of change in the pixel value over time of the target pixel 111a.
- the target pixel 111a may be represented by a in the drawing, and a plurality of peripheral pixels 111b may be represented by b1 to b8 as shown in FIGS. 2 and 3.
- the selected two or more radiographic images 10 and 20 may include a first radiographic image 10 and a second radiographic image 20, and the first radiographic image 10 is generated at a first time t1. It may be a radiographic image, and the second radiographic image 20 may be a radiographic image generated at a second time t2.
- the target pixel 111a may be represented as a1 as shown in FIG. 2B by reflecting the first time t1
- the peripheral pixels 111b may be represented as b11 to b81.
- the target pixel 111a may be represented by a2 as shown in FIG. 2(c) by reflecting the second time t2, and the peripheral pixels 111b may be expressed by b12 to b82. Can be represented.
- the pixel value of the target pixel 111a determined as the abnormal value is abnormal when the number of surrounding pixels 111b that is greater than or equal to the preset reference value R is less than or equal to the preset determination number n1. It can be confirmed (or confirmed) by the value (S620). As shown in (a) of FIG.
- the target pixel 111a and the target pixel are similar, and the reason that the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R can be considered to be attributed to the random noise.
- the pixel value of the target pixel 111a determined as the abnormal value may be determined (or determined) as the abnormal value. And, as shown in (b) of FIG.
- the target pixel when the number of surrounding pixels 111b having a change amount of the pixel value over the preset reference value R is greater than a preset determination number n1, the target pixel ( It can be considered that the amount of change in the pixel value over time in the peripheral region of 111a) and the target pixel 111a has a similar tendency, and it can be considered that there is a large change in the input signal rather than the random noise. Accordingly, the pixel value of the target pixel 111a, which is determined to be an abnormal value due to a change in the input signal, can be determined as a normal value through the overturning of the determination, thereby allowing the normal pixel value due to a change in the input signal to be an abnormal value. It can be excluded from judgment.
- the number of determinations (n1) is appropriately set to a number that can determine the tendency whether the amount of change in the pixel value over time of the target pixel 111a is similar to the amount of change in the pixel value over time in the peripheral area of the target pixel 111a. It can be set to three or more to reliably grasp similar trends.
- FIG. 4 is a conceptual diagram illustrating verification of a pixel value of a target pixel through comparison of a pixel value change amount with time of a target pixel and a change amount of pixel value with time of a surrounding pixel according to an embodiment of the present invention.
- a process (S600) of verifying a pixel value of the target pixel 111a is performed on a neighboring pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20. Comparing the amount of change in the pixel value over time and the amount of change in the pixel value over time of the target pixel (111a) (S630); And the number of surrounding pixels 111b having the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel within the preset tolerance (d) range is less than or equal to the preset number of determinations (n2). And checking the pixel value of the target pixel 111a determined as an abnormal value (S640 ).
- a pixel value change amount over time of a neighboring pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20 and a pixel value change amount over time of the target pixel 111a may be compared. (S630). Through this, it is possible to grasp the tendency of the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel value over time in the peripheral area of the target pixel 111a.
- the number of surrounding pixels 111b having the amount of change in the pixel value and the amount of change in the pixel within the preset tolerance (d) range ( ⁇ d) according to the time of the target pixel 111a is less than or equal to the preset determination number (n2)
- the pixel value of the target pixel 111a determined as the abnormal value may be checked (or confirmed) as the abnormal value (S640).
- the target pixel 111a when the number of pixel values variation over time of the target pixel 111a and the number of surrounding pixels 111b having a pixel value variation within a preset allowable deviation (d) range is less than or equal to a preset determination number (n2) ) And the amount of change in the pixel value over time in the peripheral area of the target pixel 111a cannot be confirmed, and the reason for the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R is caused by the random noise. It can be considered to be due to the abnormality, and the pixel value of the target pixel 111a determined as the abnormal value may be determined as the abnormal value.
- the number of surrounding pixels 111b having the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel within the preset tolerance (d) range is larger than the number of predetermined determinations (n2). It can be considered that the amount of change in the pixel value over time in the peripheral region of the target pixel 111a and the target pixel 111a has a similar tendency, and it can be considered that there is a large change in the input signal rather than the random noise.
- the pixel value of the target pixel 111a which is determined to be an abnormal value due to a change in the input signal, can be determined as a normal value through the overturning of the determination, thereby allowing the normal pixel value due to a change in the input signal to be an abnormal value. It can be excluded from judgment.
- the allowable deviation (d) may be appropriately set at a preset reference value (R), and may be appropriately determined within about 10% of a change in the pixel value over time of the target pixel (111a), which is not very large.
- the number of determinations n2 is suitably set to a number capable of determining the tendency of the amount of change in the pixel value over time of the target pixel 111a to be similar to the amount of change in the pixel value over time in the peripheral area of the target pixel 111a. It can be set to three or more to reliably grasp similar trends.
- the number of surrounding pixels 111b having the amount of change in pixel value over time of the target pixel 111a and the amount of change in pixel value within a preset tolerance (d) range is 8, and the pixel value of the target pixel 111a is set. It is decided by the normal value. Referring to FIG.
- the time of the peripheral pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20
- the amount of change with the pixel value with time of the target pixel 111a and the value of the pixel with time of the peripheral area of the target pixel 111a It is possible to grasp (or confirm) the tendency of the change amount more accurately.
- the present invention by selecting two or more radiographic images 10 and 20 among a plurality of radiographic images acquired over time, by considering the amount of change in pixel values between the selected two or more radiographic images 10 and 20, over time, In response to a change in the input signal, it may be determined whether the pixel value of each target pixel 111a is an abnormal pixel value including the random noise. Accordingly, by correcting by removing the abnormal pixel value including the random noise from a plurality of radiographic images, distortion of the radiographic image can be prevented, and the quality of the radiographic image can be improved.
- the (b) process (S200) to the (e) process (S500) may be repeated for the whole of the plurality of radiographic images, and the radiographic image processing method according to the present invention may generate a pixel value determined as an abnormal value.
- the process of storing the position of the pixel 111a (S700); may further include, and the pixel in which the position is stored in the (B) process (S200) after the process of storing the position of the target pixel 111a (S700)
- the pixel value of the pixel 111 in which the location is stored may be determined as an abnormal value by excluding (111) from the selection of the target pixel 111a.
- the (b) process (S200) to the (e) process (S500) may be repeated, and the (b) process (S200) to (e) process (S500) may be performed on the whole of the plurality of radiographic images. can do. Further, it may be determined whether the pixel value is an abnormal value for all the plurality of pixels 111 or all target pixels 111a for which inspection (or determination) is desired.
- the position of the target pixel 111a in which the pixel value determined (or determined) as the abnormal value has occurred may be stored (S700 ). Due to the nature of the random noise, since the probability of the random noise occurring again is high in the pixel 111 where the random noise has once occurred, the position of the target pixel 111a in which an abnormal pixel value has occurred is stored, and the position in the subsequent radiation image It may be regarded as an abnormal pixel value without performing the steps (b) (S200) to (e) (S500) of determining whether the stored pixel 111 is an abnormal pixel value.
- the pixel storing the position is excluded from the selection of the target pixel 111a, and the pixel storing the position ( The pixel value of 111) may be determined (or considered) as an abnormal value. That is, because of the nature of the random noise, since the probability of generating the random noise is high in the pixel 111 where the random noise has once occurred, the location of the target pixel 111a where the abnormal pixel value has occurred is stored in the subsequent radiation image. By considering the pixel value of the corresponding pixel 111 as an abnormal pixel value, the number of pixels 111 to determine whether the pixel value is abnormal can be reduced. In addition, the time of image processing for determining and correcting whether there are abnormal pixel values for all of the plurality of radiographic images may be shortened.
- the radiographic image processing method according to the present invention may further include a step (S150) of detecting a pixel region to which radiation is incident at an intensity within a predetermined range among a plurality of pixels 111 (S200). ), a target pixel 111a may be selected from the pixel areas.
- a pixel region in which radiation is incident at an intensity within a predetermined range among the plurality of pixels 111 may be detected (S150 ).
- a complex region in which the pixel value difference between adjacent pixels 111 is severe (or very large), such as an edge region of the subject and a boundary region between parts of the subject, and a central region of each part of the subject, such as A simple region in which the difference in pixel values between the pixels 111 is not large (or insignificant) may exist. Since the input signal to each pixel 111 is almost entirely different in the complex area, the tendency of the pixel value change over time of the target pixel 111a and the pixel value change over time of the peripheral area of the target pixel 111a can be grasped.
- the abnormal pixel value including the random noise it is not a big problem in the part where medical judgment is made from the radiographic image because it is unrecognizable to the eye in the radiographic image.
- the simple region can grasp the tendency of the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel value over time in the peripheral area of the target pixel 111a, the abnormal pixel value including the random noise
- the medical image is identified as the radiographic image containing the random noise, which is identified as an abnormal part in the radiographic image.
- a pixel region ie, the simple region in which radiation is incident at an intensity within a predetermined range among the plurality of pixels 111 is detected, and only the pixel(s) of the pixel region in which radiation is incident at an intensity within the predetermined range is pixels. It can be determined whether the value is an abnormal value.
- a target pixel 111a may be selected from the pixel areas.
- Target pixel 111a among the pixel areas to determine whether only the pixel(s) of the pixel area where radiation is incident at an intensity within the preset range capable of detecting the abnormal pixel value including the random noise is an abnormal value
- You can choose when the target pixel 111a is selected from the pixel areas, all of the plurality of pixels 111 do not determine whether the pixel values are abnormal values, but only the pixel(s) corresponding to some of the pixel areas have abnormal pixel values. By determining whether it is a value, the time for detecting the abnormal pixel value including the random noise in the selected radiation image can be shortened.
- the method of processing a radiographic image according to the present invention further includes a process of correcting a pixel value of a target pixel 111a determined as an abnormal value by interpolating a pixel value of adjacent peripheral pixels 111b on the same radiographic image (S800). Can.
- the pixel value of the target pixel 111a determined (or determined) as an abnormal value may be corrected by interpolating the pixel value of the adjacent pixels 111b on the adjacent (one) radiation image (S800).
- the pixel value of the target pixel 111a determined (or determined) as an abnormal value may be replaced with an average value or a median value of the pixel values of the adjacent peripheral pixels 111b.
- the adjacent peripheral pixel 111b includes a pixel 111 having a pixel value determined (or determined) as an abnormal value
- a peripheral pixel having a pixel value determined as an abnormal value in calculation of an average value or a median value (111b) can be excluded.
- the method of correcting the pixel value of the adjacent pixel 111b with the pixel value of the target pixel 111a determined to be an abnormal value may be performed in various ways in addition to this method.
- FIG. 5 is a block diagram showing the configuration of a radiographic imaging apparatus according to an embodiment of the present invention.
- the radiographic imaging apparatus according to another embodiment of the present invention will be described in more detail with reference to FIG. 5, and the overlapping with the above-described part in relation to the radiographic image processing method according to an embodiment of the present invention will be omitted.
- the radiographic imaging apparatus 100 includes a plurality of pixels 111, and a radiation detection panel for generating electrical signals for acquiring a plurality of radiographic images over time by irradiated radiation ( 110); The amount of change in pixel value over time of the selected target pixel 111a among the plurality of pixels 111 from two or more selected radiographic images 10 and 20 among the plurality of radiographic images and surrounding pixels adjacent to the target pixel 111a ( A change amount calculating unit 121 for calculating a pixel value change amount over time in 111b); A change amount comparison unit 122 comparing the calculated change amount of the pixel value with time of the target pixel 111a with a preset reference value R; And when the amount of change in the pixel value over time of the target pixel 111a is equal to or greater than the preset reference value R, the target pixel 111a is selected from at least one of the selected two or more radiographic images 10 and 20. It may include; a pixel value abnormality determination unit 123 for
- the radiation detection panel 110 may include a plurality of pixels 111 and may generate electrical signals for acquiring a plurality of radiation images over time by irradiated radiation.
- a radiographic image may be generated and obtained by the electric signal thus generated.
- the change amount calculating unit 121 may change the pixel value and the target pixel over time of the selected target pixel 111a among the plurality of pixels 111 from two or more selected radiation images 10 and 20 among the plurality of radiation images.
- the amount of change in the pixel value over time of the peripheral pixel 111b adjacent to 111a) may be calculated.
- the radiographic imaging apparatus 100 of the present invention is a pixel setting unit (not shown) that selects a target pixel 111a among a plurality of pixels 111 from two or more radiographic images 10 and 20 selected from the plurality of radiographic images. ); may be further included.
- the change amount comparison unit 122 may compare the calculated change amount of the pixel value according to the time of the target pixel 111a with a preset reference value R, and the pixel according to the time of the peripheral pixel 111b adjacent to the target pixel 111a The value change amount may be compared with a pixel value change amount and/or a preset reference value R over time of the target pixel 111a.
- the pixel value anomaly determining unit 123 may select the target pixel 111a from the selected two or more radiographic images 10 and 20 when the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to a preset reference value R.
- the pixel value of can be determined as an abnormal value.
- the radiographic imaging apparatus 100 of the present invention may further include an image processing unit 120 that performs image processing of the generated or obtained radiographic image, and the image processing unit 120 compares the change amount calculation unit 121 and the change amount
- the unit 122 and the pixel value abnormality determination unit 123 may be included.
- the radiographic imaging apparatus 100 is determined to be the abnormal value when the number of surrounding pixels 111b equal to or greater than the preset reference value R is less than or equal to the preset determination number n1. It may further include a first pixel value verification unit (not shown) that determines the pixel value of the target pixel 111a as an abnormal value.
- the first pixel value verification unit determines the abnormal value when the number of surrounding pixels 111b equal to or greater than the preset reference value R is less than or equal to the preset determination number n1.
- the pixel value of the target pixel 111a may be determined as an abnormal value to accurately verify whether the pixel value of the target pixel 111a is an abnormal value.
- the radiographic imaging apparatus 100 is the number of surrounding pixels 111b having a pixel value change amount over time and a predetermined value (d) within a range ( ⁇ d) of the target pixel 111a.
- a second pixel value verification unit (not shown) that determines the pixel value of the target pixel 111a determined as the abnormal value as an abnormal value when is less than or equal to the preset number of determinations n2 may be further included.
- the second pixel value verifying unit determines that the number of surrounding pixels 111b having the amount of change in pixel values over time of the target pixel 111a and the amount of change in pixel values within a preset tolerance (d) range is preset.
- the number n2 is less than or equal to the pixel value of the target pixel 111a determined as the abnormal value, the pixel value of the target pixel 111a can be accurately verified.
- the radiographic imaging apparatus 100 may further include a location storage unit (not shown) that stores the location of the target pixel 111a where the pixel value determined as the abnormal value has occurred.
- the location storage unit may store the location of the target pixel 111a where the pixel value determined as the abnormal value occurs. Due to the nature of random noise, since the probability of the random noise occurring again is high in the pixel 111 where the random noise has once occurred, the position of the target pixel 111a where the abnormal pixel value has occurred is stored, and the position in the subsequent radiation image is The stored pixel 111 may be regarded as an abnormal pixel value without a separate process of determining whether it is an abnormal pixel value.
- the radiographic imaging apparatus 100 is a reference value setting unit for setting the preset reference value greater than the intrinsic noise value generated even when radiation is irradiated from each pixel as a whole to the plurality of pixels or generated by a radiation source (Not shown); may further include.
- the reference value setting unit (not shown) is a preset reference value (or greater than or equal to) that is generated even when radiation is not irradiated from each pixel 111 as a whole, or greater than or equal to the intrinsic noise value generated by the radiation source.
- R can be set.
- the preset reference value R may be set to be larger (or higher) than the intrinsic noise value common to all the plurality of pixels 111.
- the intrinsic noise may have a different value, but refers to noise common to all the plurality of pixels 111, and may be noise generated by a device offset, and the noise may be any of the plurality of pixels 111 ), it is not a big problem if all of the device offset values are removed from the pixel values of the plurality of pixels 111 constituting the radiation image.
- the intrinsic noise is generated due to a cause of thermal (thermal), poor insulation, and the like due to radiation, and dark current and quantum of photons accumulating charge in the plurality of pixels 111 It is caused by nature, and may include quantum noise generated by a radiation generator.
- the present invention by selecting two or more radiographic images among a plurality of radiographic images acquired over time, and taking into account the amount of change in pixel values between the selected two or more radiographic images, the change in the input signal over time is reflected. It may be determined whether the pixel value of each target pixel is an abnormal pixel value including random noise. Accordingly, by correcting by removing abnormal pixel values including random noise from a plurality of radiographic images, distortion of the radiographic image can be prevented, and quality of the radiographic image can be improved.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Measurement Of Radiation (AREA)
Abstract
The present invention relates to a radiographic image processing method and a radiography apparatus, whereby abnormal pixel values are detected from a plurality of radiographic images acquired by time. At least two radiographic images, among the plurality of radiographic images acquired by time, are selected, and by considering pixel value variations between the selected at least two radiographic images, a determination may be made, by reflecting variations in an input signal by time, on whether the respective pixel values of target pixels are abnormal pixel values including random noise.
Description
본 발명은 방사선 영상 처리방법 및 방사선 촬영장치에 관한 것으로서, 보다 상세하게는 시간에 따라 획득된 복수의 방사선 이미지에서 비정상 픽셀값을 검출하는 방사선 영상 처리방법 및 방사선 촬영장치에 관한 것이다.The present invention relates to a radiographic image processing method and a radiographic imaging apparatus, and more particularly, to a radiographic image processing method and a radiographic imaging apparatus for detecting abnormal pixel values from a plurality of radiographic images obtained over time.
일반적으로, 방사선 촬영장치는 방사선을 피사체에 조사하는 방사선 발생장치(Radiation Generator)와 피사체를 투과한 방사선을 방사선 검출패널로 검출하는 방사선 디텍터(Radiation Detector)로 구성된다.In general, a radiographic imaging apparatus includes a radiation generator that irradiates radiation to a subject and a radiation detector that detects radiation that has transmitted through the subject with a radiation detection panel.
방사선 검출패널은 상보성 금속산화물 반도체(Complementary Metal-Oxide Semiconductor; CMOS) 또는 다결정 실리콘(Poly-Si)과 같은 결정체 등의 소재로 이루어질 수 있으며, 소재의 특성에 따른 노이즈가 발생할 수 있다.The radiation detection panel may be made of a material such as a complementary metal-oxide semiconductor (CMOS) or a crystal such as polycrystalline silicon (Poly-Si), and noise may occur according to the characteristics of the material.
이러한 노이즈는 그 값이 비교적 높게 나타나므로, 이러한 노이즈가 포함된 픽셀값을 사용하여 방사선 이미지를 생성할 경우에는 방사선 이미지를 왜곡시킬 수 있고, 이에 따라 방사선 이미지의 품질을 저하시킬 수 있다. Since the noise has a relatively high value, when a radiographic image is generated using a pixel value including such noise, the radiographic image may be distorted, thereby degrading the quality of the radiographic image.
이를 해결하고자 방사선 이미지에서 노이즈를 제거하기 위한 방법이 시도되고 있으나, 종래의 방법들은 단일 방사선 이미지에서 영역별 경향성(tendency)만을 판단하여 노이즈 픽셀을 검출하고 있으며, 시간에 따른 입력 신호(input signal)의 변화를 고려하지 않아 시간에 따라 획득된 복수의 방사선 이미지에서 시간에 대해 무작위(random)하게 발생하는 랜덤 노이즈의 경우에는 노이즈 픽셀을 검출하는 데에 한계가 있다.In order to solve this, a method for removing noise from a radiographic image has been attempted, but conventional methods detect noise pixels by determining only a tendency for each region in a single radiographic image, and an input signal over time. In the case of random noise generated randomly with respect to time in a plurality of radiographic images acquired over time because the change of the time is not taken into consideration, there is a limit in detecting a noise pixel.
(특허문헌 1) 한국등록특허공보 제10-0617517호(Patent Document 1) Korean Registered Patent Publication No. 10-0617517
본 발명은 시간에 따라 획득된 복수의 방사선 이미지에서 비정상 픽셀값을 검출하는 방사선 영상 처리방법 및 방사선 촬영장치를 제공한다.The present invention provides a radiographic image processing method and a radiographic imaging apparatus for detecting abnormal pixel values from a plurality of radiographic images acquired over time.
본 발명의 일실시예에 따른 방사선 영상 처리방법은 (a) 복수의 픽셀을 구비하는 방사선 검출패널을 이용하여 시간에 따라 복수의 방사선 이미지를 획득하는 과정; (b) 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지에서 상기 복수의 픽셀 중 대상 픽셀을 선택하는 과정; (c) 선택된 대상 픽셀의 시간에 따른 픽셀값 변화량을 산출하는 과정; (d) 산출된 상기 대상 픽셀의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값과 비교하는 과정; 및 (e) 상기 대상 픽셀의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지 중 적어도 하나의 방사선 이미지에서 상기 대상 픽셀의 픽셀값을 비정상값으로 판단하는 과정;을 포함할 수 있다.A radiographic image processing method according to an embodiment of the present invention includes: (a) obtaining a plurality of radiographic images over time using a radiation detection panel having a plurality of pixels; (b) selecting a target pixel among the plurality of pixels from two or more selected radiation images among the plurality of radiation images; (c) calculating a pixel value change amount over time of the selected target pixel; (d) comparing the calculated amount of change of the pixel value with time of the target pixel with a preset reference value; And (e) determining a pixel value of the target pixel as an abnormal value in at least one of the selected two or more radiographic images when the amount of change in the pixel value over time of the target pixel is equal to or greater than the preset reference value. It can contain.
비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 검증하는 과정;을 더 포함하고, 상기 대상 픽셀의 픽셀값을 검증하는 과정은, 상기 선택된 둘 이상의 방사선 이미지로부터 산출된 상기 대상 픽셀에 인접한 주변 픽셀의 시간에 따른 픽셀값 변화량과 상기 미리 설정된 기준값을 비교하는 과정; 및 상기 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 확인하는 과정을 포함할 수 있다.Further comprising the step of verifying the pixel value of the target pixel determined to be an abnormal value, and further comprising the step of verifying the pixel value of the target pixel, the neighboring pixels adjacent to the target pixel calculated from the selected two or more radiographic images Comparing the pixel value change with time and the preset reference value; And checking the pixel value of the target pixel determined as the abnormal value as an abnormal value when the number of surrounding pixels having a change in the pixel value over time is equal to or less than the preset determination number. .
비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 검증하는 과정;을 더 포함하고, 상기 대상 픽셀의 픽셀값을 검증하는 과정은, 상기 선택된 둘 이상의 방사선 이미지로부터 산출된 상기 대상 픽셀에 인접한 주변 픽셀의 시간에 따른 픽셀값 변화량과 상기 대상 픽셀의 시간에 따른 픽셀값 변화량을 비교하는 과정; 및 상기 대상 픽셀의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 확인하는 과정을 포함할 수 있다.Further comprising the step of verifying the pixel value of the target pixel determined to be an abnormal value, and further comprising the step of verifying the pixel value of the target pixel, the neighboring pixels adjacent to the target pixel calculated from the selected two or more radiographic images Comparing a pixel value change amount with time and a pixel value change amount with time of the target pixel; And the pixel value of the target pixel determined as the abnormal value is abnormal when the number of pixel values variation over time and the number of surrounding pixels having a pixel value variation within a preset allowable deviation range is equal to or less than a preset determination number. It may include the process of confirming by value.
상기 복수의 방사선 이미지 전체에 대해 상기 (b) 과정 내지 상기 (e) 과정을 반복하며, 비정상값으로 판단된 픽셀값이 발생한 대상 픽셀의 위치를 저장하는 과정;을 더 포함하고, 상기 대상 픽셀의 위치를 저장하는 과정 이후의 (b) 과정에서 위치가 저장된 픽셀을 상기 대상 픽셀의 선택에서 제외하고 상기 위치가 저장된 픽셀의 픽셀값을 비정상값으로 판단할 수 있다.The process of (b) to (e) is repeated for the whole of the plurality of radiographic images, and a process of storing the position of the target pixel in which the pixel value determined as the abnormal value occurs is further included. In the step (b) after the process of storing the position, the pixel storing the position may be excluded from the selection of the target pixel, and the pixel value of the pixel storing the position may be determined as an abnormal value.
상기 미리 설정된 기준값은 상기 복수의 픽셀 전체적으로 각각의 픽셀에서 방사선이 비조사되는 경우에도 발생되거나, 방사선원으로 인해 발생되는 고유 노이즈값보다 클 수 있다.The preset reference value may be generated even when radiation is not irradiated from each pixel as a whole to the plurality of pixels, or may be greater than the inherent noise value generated by the radiation source.
상기 (b) 과정은 상기 복수의 방사선 이미지 중 시간적으로 연속되는 둘 이상의 방사선 이미지를 선택하여 수행될 수 있다.The (b) process may be performed by selecting two or more radiographic images that are temporally continuous among the plurality of radiographic images.
상기 (b) 과정에서는 상기 복수의 방사선 이미지 중 시간적으로 연속되는 2개의 방사선 이미지를 선택하고, 상기 (e) 과정에서는 상기 대상 픽셀의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 2개의 방사선 이미지 중 상기 대상 픽셀의 픽셀값이 높은 방사선 이미지에서 상기 대상 픽셀의 픽셀값을 비정상값으로 판단할 수 있다.In the step (b), two radiographic images that are temporally continuous among the plurality of radiographic images are selected, and in the step (e), two radiographic images in which the amount of change in a pixel value over time of the target pixel is greater than or equal to the preset reference value Among the pixels, the pixel value of the target pixel may be determined as an abnormal value in a radiographic image having a high pixel value of the target pixel.
상기 복수의 픽셀 중 미리 설정된 범위 내의 세기로 방사선이 입사되는 픽셀 영역을 검출하는 과정;을 더 포함하고, 상기 (b) 과정에서는 상기 픽셀 영역 중 상기 대상 픽셀을 선택할 수 있다.And detecting a pixel region to which radiation is incident at an intensity within a predetermined range among the plurality of pixels. In the process (b), the target pixel may be selected from the pixel regions.
비정상값으로 판단된 대상 픽셀의 픽셀값을 동일한 방사선 이미지 상의 인접한 주변 픽셀의 픽셀값을 내삽하여 보정하는 과정;을 더 포함할 수 있다.The process may further include correcting a pixel value of a target pixel determined as an abnormal value by interpolating a pixel value of adjacent neighboring pixels on the same radiation image.
본 발명의 다른 실시예에 따른 방사선 촬영장치는 복수의 픽셀을 구비하며, 조사되는 방사선에 의해 시간에 따라 복수의 방사선 이미지를 획득하기 위한 전기신호를 생성하는 방사선 검출패널; 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지에서 상기 복수의 픽셀 중 선택된 대상 픽셀의 시간에 따른 픽셀값 변화량 및 상기 대상 픽셀에 인접한 주변 픽셀의 시간에 따른 픽셀값 변화량을 산출하는 변화량 산출부; 산출된 상기 대상 픽셀의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값과 비교하는 변화량 비교부; 및 상기 대상 픽셀의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지 중 적어도 하나의 방사선 이미지에서 상기 대상 픽셀의 픽셀값을 비정상값으로 판단하는 픽셀값 이상판단부;를 포함할 수 있다.A radiographic imaging apparatus according to another embodiment of the present invention includes a plurality of pixels, and a radiation detection panel that generates electrical signals for acquiring a plurality of radiographic images over time by irradiated radiation; A change amount calculating unit for calculating a pixel value change amount over time of a target pixel selected from the plurality of pixels from two or more radiation images selected from the plurality of radiation images and a pixel value change amount over time of neighboring pixels adjacent to the target pixel; A change amount comparison unit comparing the calculated change amount of the pixel value with time of the target pixel with a preset reference value; And a pixel value abnormality determination unit determining a pixel value of the target pixel as an abnormal value in at least one radiation image among the selected two or more radiographic images when the amount of change of the pixel value over time of the target pixel is equal to or greater than the preset reference value. It may include.
상기 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 결정하는 제1 픽셀값 검증부;를 더 포함할 수 있다.A first pixel value verification unit configured to determine a pixel value of the target pixel determined as the abnormal value as an abnormal value when the number of surrounding pixels having a change amount of the pixel value over time is equal to or less than a preset determination number; It may further include.
상기 대상 픽셀의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 결정하는 제2 픽셀값 검증부;를 더 포함할 수 있다.The pixel value of the target pixel determined as the abnormal value is an abnormal value when the number of pixel values variation over time and the number of surrounding pixels having a pixel value variation within a preset allowable deviation range is equal to or less than a preset determination number. It may further include; a second pixel value verification unit determined by.
비정상값으로 판단된 픽셀값이 발생한 대상 픽셀의 위치를 저장하는 위치저장부;를 더 포함할 수 있다.It may further include a location storage unit for storing the location of the target pixel where the pixel value determined as the abnormal value occurs.
상기 복수의 픽셀 전체적으로 각각의 픽셀에서 방사선이 비조사되는 경우에도 발생되거나, 방사선원으로 인해 발생되는 고유 노이즈값보다 큰 상기 미리 설정된 기준값을 설정하는 기준값 설정부;를 더 포함할 수 있다.It may further include a reference value setting unit for setting the preset reference value that is greater than the intrinsic noise value generated even when radiation is irradiated from each pixel as a whole, or is generated by the radiation source.
본 발명의 일실시예에 따른 방사선 영상 처리방법은 시간에 따라 획득된 복수의 방사선 이미지 중 둘 이상의 방사선 이미지를 선택하여, 선택된 둘 이상의 방사선 이미지 간의 픽셀값 변화량을 고려함으로써, 시간에 따른 입력 신호(input signal)의 변화를 반영하여 각 대상 픽셀의 픽셀값이 랜덤 노이즈를 포함하는 비정상 픽셀값인지를 판단할 수 있다. 이에 따라 복수의 방사선 이미지에서 랜덤 노이즈가 포함된 비정상 픽셀값을 제거하여 보정함으로써, 방사선 이미지의 왜곡을 방지할 수 있고, 방사선 이미지의 품질을 향상시킬 수 있다.In the radiographic image processing method according to an embodiment of the present invention, by selecting two or more radiographic images among a plurality of radiographic images acquired over time, an input signal over time ( Input signal) may be reflected to determine whether a pixel value of each target pixel is an abnormal pixel value including random noise. Accordingly, by correcting by removing abnormal pixel values including random noise from a plurality of radiographic images, distortion of the radiographic image can be prevented, and quality of the radiographic image can be improved.
또한, 대상 픽셀의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값을 비교하여 비정상 픽셀값을 판단한 후에 다시 주변 픽셀을 통해 비정상 픽셀값이 맞는지를 검증함으로써, 피사체의 움직임 등에 의한 입력 신호의 변화를 효과적으로 반영하여 각 대상 픽셀의 픽셀값이 랜덤 노이즈를 포함하는 비정상 픽셀값인지를 보다 정확하게 판단할 수 있다.In addition, by comparing the amount of change of the pixel value with time of the target pixel and a preset reference value, and determining the abnormal pixel value again, verifying that the abnormal pixel value is correct through the surrounding pixels, effectively reflecting the change in the input signal due to the movement of the subject. Thus, it can be more accurately determined whether the pixel value of each target pixel is an abnormal pixel value including random noise.
그리고 랜덤 노이즈의 특성상 한 번 랜덤 노이즈가 발생했던 픽셀에서는 랜덤 노이즈가 발생할 확률이 높기 때문에 비정상 픽셀값이 발생한 대상 픽셀의 위치를 저장하고 이후의 방사선 이미지에서 해당 픽셀의 픽셀값을 비정상 픽셀값으로 간주함으로써, 비정상 픽셀값인지를 판단해야 하는 픽셀의 개수가 줄어들 수 있으며, 복수의 방사선 이미지 전체에 대해 비정상 픽셀값이 있는지를 판단하고 보정하는 영상 처리의 시간이 단축될 수 있다.Also, because of the nature of random noise, the probability of random noise is high in a pixel where random noise has occurred once, so the location of the target pixel where the abnormal pixel value is generated is stored and the pixel value of the pixel is regarded as an abnormal pixel value in the subsequent radiation image. By doing so, the number of pixels that need to be determined to be abnormal pixel values can be reduced, and the time for image processing to determine and correct whether there are abnormal pixel values for the entire plurality of radiographic images can be shortened.
도 1은 본 발명의 일실시예에 따른 방사선 영상 처리방법을 나타낸 순서도.1 is a flowchart illustrating a radiographic image processing method according to an embodiment of the present invention.
도 2는 본 발명의 일실시예에 따른 대상 픽셀과 주변 픽셀을 설명하기 위한 개념도.2 is a conceptual diagram illustrating a target pixel and surrounding pixels according to an embodiment of the present invention.
도 3은 본 발명의 일실시예에 따른 미리 설정된 기준값을 통한 대상 픽셀의 픽셀값 검증을 설명하기 위한 개념도.3 is a conceptual diagram illustrating verification of a pixel value of a target pixel through a preset reference value according to an embodiment of the present invention.
도 4는 본 발명의 일실시예에 따른 대상 픽셀의 시간에 따른 픽셀값 변화량과 주변 픽셀의 시간에 따른 픽셀값 변화량의 비교를 통한 대상 픽셀의 픽셀값 검증을 설명하기 위한 개념도.4 is a conceptual diagram illustrating verification of a pixel value of a target pixel through comparison of a pixel value change amount with time of a target pixel and a change amount of pixel value with time of a surrounding pixel according to an embodiment of the present invention.
도 5는 본 발명의 일실시예에 따른 방사선 촬영장치의 구성을 나타낸 블록도.5 is a block diagram showing the configuration of a radiographic apparatus according to an embodiment of the present invention.
이하에서는 첨부된 도면을 참조하여 본 발명의 실시예를 더욱 상세히 설명하기로 한다. 그러나 본 발명은 이하에서 개시되는 실시예에 한정되는 것이 아니라 서로 다른 다양한 형태로 구현될 것이며, 단지 본 실시예들은 본 발명의 개시가 완전하도록 하며, 통상의 지식을 가진 자에게 발명의 범주를 완전하게 알려주기 위해 제공되는 것이다. 설명 중, 동일 구성에 대해서는 동일한 참조부호를 부여하도록 하고, 도면은 본 발명의 실시예를 정확히 설명하기 위하여 크기가 부분적으로 과장될 수 있으며, 도면상에서 동일 부호는 동일한 요소를 지칭한다.Hereinafter, embodiments of the present invention will be described in more detail with reference to the accompanying drawings. However, the present invention is not limited to the embodiments disclosed below, but will be implemented in various different forms, and only the present embodiments allow the disclosure of the present invention to be complete, and the scope of the invention to those skilled in the art is completely It is provided to inform you. In the description, the same reference numerals are assigned to the same components, and the drawings may be exaggerated in size in order to accurately describe the embodiments of the present invention, and the same reference numerals in the drawings refer to the same elements.
도 1은 본 발명의 일실시예에 따른 방사선 영상 처리방법을 나타낸 순서도이다.1 is a flowchart illustrating a radiographic image processing method according to an embodiment of the present invention.
도 1을 참조하면, 본 발명의 일실시예에 따른 방사선 영상 처리방법은 (a) 복수의 픽셀(111)을 구비하는 방사선 검출패널(110)을 이용하여 시간에 따라 복수의 방사선 이미지를 획득하는 과정(S100); (b) 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지(10,20)에서 상기 복수의 픽셀(111) 중 대상 픽셀(111a)을 선택하는 과정(S200); (c) 선택된 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 산출하는 과정(S300); (d) 산출된 상기 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값(R)과 비교하는 과정(S400); 및 (e) 상기 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값(R) 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지(10,20) 중 적어도 하나의 방사선 이미지에서 상기 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단하는 과정(S500);을 포함할 수 있다.Referring to FIG. 1, a radiographic image processing method according to an embodiment of the present invention includes (a) obtaining a plurality of radiographic images over time using a radiation detection panel 110 having a plurality of pixels 111. Process (S100); (b) selecting a target pixel 111a among the plurality of pixels 111 from two or more selected radiation images 10 and 20 among the plurality of radiation images (S200); (c) calculating the amount of change in the pixel value over time of the selected target pixel 111a (S300); (d) comparing the calculated amount of change of the pixel value over time of the target pixel 111a with a preset reference value R (S400); And (e) when the amount of change in pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R, the target pixel in at least one of the selected two or more radiographic images 10 and 20 ( And determining the pixel value of 111a) as an abnormal value (S500).
먼저, (a) 복수의 픽셀(111)을 구비하는 방사선 검출패널(110)을 이용하여 시간에 따라 복수의 방사선 이미지를 획득한다(S100). 방사선 발생장치를 이용하여 방사선 검출패널(110)에 방사선을 조사함으로써, 방사선 이미지를 생성할 수 있으며, 방사선의 조사에 의해 생성되는 전하가 방사선 검출패널(110)의 복수의 픽셀(111)에 저장될 수 있고, 복수의 픽셀(111)에 저장된 전하(량)에 따라 상기 방사선 이미지가 생성될 수 있다. 여기서, 생성된 방사선 이미지를 바로 획득할 수도 있고, 생성되어 저장된 방사선 이미지를 획득할 수도 있다. 이러한 상기 복수의 ㅂ방사선 이미지는 형광 투시법(Fluoroscopy), 전산화단층영상촬영기(Cone Beam Computed Tomography; CBCT) 등의 컴퓨터 단층촬영(Computed Tomography; CT)에서 사용될 수 있다.First, (a) a plurality of radiation images are acquired over time using the radiation detection panel 110 including a plurality of pixels 111 (S100). By irradiating radiation to the radiation detection panel 110 using a radiation generator, a radiation image can be generated, and charges generated by irradiation of radiation are stored in a plurality of pixels 111 of the radiation detection panel 110. The radiation image may be generated according to charges (amounts) stored in the plurality of pixels 111. Here, the generated radiation image may be directly obtained, or the generated and stored radiation image may be obtained. Such a plurality of X-ray images may be used in computed tomography (CT) such as fluorescence (Fluoroscopy), computed tomography (Cone Beam Computed Tomography; CBCT).
다음으로, (b) 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지(10,20)에서 상기 복수의 픽셀(111) 중 대상 픽셀(111a)을 선택한다(S200). 상기 복수의 방사선 이미지 중 비정상 픽셀값이 있는지를 확인할 둘 이상의 방사선 이미지(10,20)를 선택할 수 있으며, 선택된 둘 이상의 방사선 이미지(10,20)에서 각 픽셀(111)의 픽셀값에 대해 비정상값인지를 판단할 수 있다. 여기서, 대상 픽셀(111a)은 하나일 수도 있고, 복수일 수도 있다.Next, (b) a target pixel 111a is selected from the plurality of pixels 111 from two or more selected radiation images 10 and 20 among the plurality of radiation images (S200). Among the plurality of radiographic images, two or more radiographic images 10 and 20 to check whether there are abnormal pixel values may be selected, and an abnormality for a pixel value of each pixel 111 in the selected two or more radiographic images 10 and 20 You can judge cognition. Here, the target pixel 111a may be one or a plurality of pixels.
그 다음 (c) 선택된 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 산출한다(S300). 각 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 각각 산출할 수도 있고, 모든 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 한 번에 산출할 수도 있다. 이를 통해 대상 픽셀(111a) 전체의 시간에 따른 픽셀값 변화량을 산출할 수 있다.Next, (c) a pixel value change amount according to the time of the selected target pixel 111a is calculated (S300). The amount of change in pixel values over time of each target pixel 111a may be calculated, or the amount of change in pixel values over time of all target pixels 111a may be calculated at once. Through this, the amount of change in the pixel value over time of the entire target pixel 111a may be calculated.
그리고 (d) 산출된 상기 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값(R)과 비교한다(S400). 산출된 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값(R)과 비교할 수 있으며, 각 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값(R)과 비교할 수 있고, 모든 대상 픽셀(111a)에 대해 미리 설정된 기준값(R)과 비교할 수 있다. 여기서, 시간에 따라 기본적으로 모든 픽셀(111)에 개재될 수 있는 요소(또는 노이즈)의 값을 제거하여 실질적으로 다른 픽셀들(예를 들어, 주변 픽셀)의 시간에 따른 픽셀값 변화량과의 차이를 발생시키는 요소(예를 들어, 랜덤 노이즈)를 확인하기 위해 상기 미리 설정된 기준값(R)을 사용할 수 있다.Then, (d) the calculated change in pixel value over time of the target pixel 111a is compared with a preset reference value R (S400). The calculated amount of change in pixel value over time of the target pixel 111a can be compared with a preset reference value R, and the amount of change in pixel value over time of each target pixel 111a can be compared with a preset reference value R, , It can be compared with a preset reference value R for all target pixels 111a. Here, basically, by removing the value of the element (or noise) that may be interposed in all the pixels 111 basically with time, the difference from the pixel value change amount over time of substantially other pixels (eg, surrounding pixels) The predetermined reference value R may be used to identify an element (for example, random noise) generating.
그 다음 (e) 상기 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값(R) 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지(10,20)에서 상기 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단한다(S500). 미리 설정된 기준값(R)에 의해 모든 픽셀(111)에 공통적으로 포함되는 요소는 배제될 수 있으며, 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 경우에는 대상 픽셀(111a)에 다른 픽셀들보다 다른 요소(또는 노이즈)가 더 개재된 것이므로, 상기 선택된 둘 이상의 방사선 이미지(10,20)에서 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단할 수 있다. 즉, 대상 픽셀(111a)의 픽셀값이 주변 픽셀(111b)의 픽셀값과 공통적으로 변하는 부분은 미리 설정된 기준값(R)에 의해 제거되어 대상 픽셀(111a)의 픽셀값만이 변한 부분을 알 수 있고, 대상 픽셀(111a)의 픽셀값만이 독립적으로 변한 부분을 통해 상기 선택된 둘 이상의 방사선 이미지(10,20)에서 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단할 수 있다.Then (e) the pixel value of the target pixel 111a in the selected two or more radiographic images 10 and 20 when the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R It is determined as an abnormal value (S500). Elements commonly included in all the pixels 111 may be excluded by the preset reference value R, and when the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R, the target pixel ( Since a different element (or noise) is interposed more than other pixels in 111a), a pixel value of the target pixel 111a may be determined as an abnormal value in the selected two or more radiographic images 10 and 20. That is, the portion in which the pixel value of the target pixel 111a changes in common with the pixel value of the peripheral pixel 111b is removed by the preset reference value R, so that the portion where only the pixel value of the target pixel 111a is changed can be seen. In addition, the pixel value of the target pixel 111a may be determined as an abnormal value in the selected two or more radiographic images 10 and 20 through a portion in which only the pixel value of the target pixel 111a is independently changed.
본 발명에서는 상기 선택된 둘 이상의 방사선 이미지(10,20) 간의 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 고려함으로써, 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값(R)의 비교를 통해 간단하게 대상 픽셀(111a)의 픽셀값이 비정상값인지를 판단할 수 있다.In the present invention, by taking into account the amount of change in the pixel value over time of the target pixel 111a between the selected two or more radiographic images 10 and 20, the amount of change in pixel value over time of the target pixel 111a and a preset reference value (R) Through comparison, it may be determined whether the pixel value of the target pixel 111a is abnormal.
상기 (b) 과정(S200)은 상기 복수의 방사선 이미지 중 시간적으로 연속되는 둘 이상의 방사선 이미지(10,20)를 선택하여 수행될 수 있다. 상기 복수의 방사선 이미지 중에서 시간적으로 연속되는 둘 이상의 방사선 이미지(10,20)를 선택할 수 있으며, 시간적으로 연속되는 둘 이상의 방사선 이미지(10,20) 간에는 시간적 경향성(tendency)이 유지될(또는 이어질) 수 있어 시간에 따른 픽셀값 변화량이 비정상 픽셀값을 검출하는 데에 보다 효과적일 수 있다.The step (b) (S200) may be performed by selecting two or more consecutive radiographic images 10 and 20 among the multiple radiographic images. Among the plurality of radiographic images, two or more temporally continuous radiographic images 10 and 20 may be selected, and a temporal tendency may be maintained (or continued) between two or more temporally continuous radiographic images 10 and 20. Therefore, the amount of change in pixel value over time can be more effective in detecting abnormal pixel values.
예를 들어, 시간적으로 비연속적인 방사선 이미지 간에는 장치적 조작 및/또는 오류(error)에 의한 방사선의 조사 세기 변화, 고의 및/또는 실수에 의한 피사체의 위치 변화(예를 들어, 움직임 또는 이동) 등으로 인해 각 픽셀(111)마다 입력 신호(input signal)의 변화가 발생할 확률이 높을 수 밖에 없고, 시간에 따라 방사선의 조사 세기 변화 및/또는 피사체의 위치 변화가 점차적으로 발생하는 경우에는 시간에 따른 픽셀값 변화량이 커질 수 밖에 없다. 이러한 경우에는 시간에 따른 픽셀값 변화량이 대부분 미리 설정된 기준값(R)을 넘게 되어 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값(R)의 비교만으로는 대상 픽셀(111a)의 픽셀값이 비정상값인지를 판단하기에 어려움이 있다. 또한, 전체적인 픽셀값 변화량이 커서 픽셀값 변화량에서 랜덤 노이즈(random noise)가 차지하는 비율이 작아질 수 있으며, 이로 인해 주변 픽셀(111b)을 통한 대상 픽셀(111a)의 픽셀값 검증이 어려워질 수도 있다.For example, between temporally discontinuous radiation images, changes in the intensity of irradiation of radiation due to device manipulation and/or errors, and intentional and/or accidental changes in the position of a subject (eg, movement or movement) ), the probability of occurrence of a change in the input signal for each pixel 111 is high, and when a change in radiation intensity and/or a position of the subject gradually changes with time The amount of change in pixel value according to is inevitably increased. In this case, the amount of change in the pixel value over time mostly exceeds the preset reference value R, so that only the comparison of the amount of change in the pixel value over time with the preset reference value R of the target pixel 111a results in the pixel value of the target pixel 111a. It is difficult to determine whether this is an abnormal value. In addition, the overall pixel value variation is large, so the ratio of random noise to the pixel value variation may be small, which may make it difficult to verify the pixel value of the target pixel 111a through the peripheral pixel 111b. .
여기서, 상기 랜덤 노이즈는 각 픽셀(111)마다 독립적으로 시간에 대해(또는 따라) 무작위(random)하게 발생하는 노이즈(noise)를 말하며, 예를 들어, 상보성 금속산화물 반도체(Complementary Metal-Oxide Semiconductor; CMOS) 또는 다결정 실리콘(Poly-Si)과 같은 결정체 등 방사선 검출패널(110)의 소재 특성에 따라 발생될 수 있고, 이를 무작위 전신 신호(Random Telegraph Signal; RTS) 노이즈(이하 RTS 노이즈)라고 말할 수도 있다. 이러한 랜덤 노이즈는 그 값이 비교적 높게 나타나게 되고, 이로 인해 랜덤 노이즈가 포함된 픽셀값을 사용하여 방사선 이미지를 생성할 경우에는 방사선 이미지를 왜곡시킬 수 있고, 방사선 이미지의 품질을 저하시킬 수 있다. 여기서, RTS 노이즈는 스위칭(Switching) 소자의 채널(Channel) 내의 하나 또는 여러 개의 전자 또는 정공이 산화물(Oxide) 내에 트래핑(Trapping) 또는 디트래핑(De-trapping)되면서 발생하는 노이즈로서, 불연속적으로 여러 개의 신호 레벨(Signal Level)을 매우 낮은 주파수로 반복적으로 나타내는 특성을 보인다. 이로 인해 방사선 이미지에서는 픽셀(111)의 신호(Signal)가 여러 개의 레벨(Level)로 랜덤하게 나타나는 결함(Defect)으로 이어지게 된다. 또한, RTS 노이즈는 거의 모든 반도체에서 발생하는 현상이며, 검출소자(sensor)에서 활용하는 경우에는 일반적으로 샘플링(Sampling) 수를 늘려 중간값(Median value)을 취하는 등의 방법으로 RTS 노이즈를 제거하는 방식을 이용하지만, 간접 픽셀 센서(Passive Pixel Sensor; PPS) 방식(type)의 이미지 센서(Image Sensor)에서는 샘플링(Sampling)을 여러 번 수행하는 것이 불가능하다.Here, the random noise refers to noise generated randomly with respect to time (or along) for each pixel 111 independently, for example, a complementary metal-oxide semiconductor; CMOS) or polycrystalline silicon (Poly-Si), such as crystals, may be generated depending on the material characteristics of the radiation detection panel 110, and may also be referred to as random telegraph signal (RTS) noise (hereinafter RTS noise). have. The random noise has a relatively high value. Therefore, when a radiographic image is generated using a pixel value including random noise, the radiographic image may be distorted and the quality of the radiographic image may be deteriorated. Here, RTS noise is a noise generated when one or several electrons or holes in a channel of a switching device are trapped or de-traded in an oxide, and discontinuously. It shows characteristics that repeatedly display multiple signal levels at very low frequencies. As a result, in the radiographic image, the signal of the pixel 111 leads to a defect that appears randomly at several levels. In addition, RTS noise is a phenomenon that occurs in almost all semiconductors, and when utilized in a sensor, RTS noise is generally removed by increasing the number of sampling and taking a median value. Although a method is used, it is impossible to perform sampling several times in an image sensor of a passive pixel sensor (PPS) type.
상기 (b) 과정(S200)에서는 상기 복수의 방사선 이미지 중 시간적으로 연속되는 2개의 방사선 이미지(10,20)를 선택할 수 있고, 상기 (e) 과정(S500)에서는 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 2개의 방사선 이미지(10,20) 중 대상 픽셀(111a)의 픽셀값이 상대적으로 높은 방사선 이미지(10 or 20)에서 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단할 수 있다. 일반적으로 2개의 방사선 이미지(10,20)에서 대상 픽셀(111a)의 픽셀값을 모두 비정상값으로 판단할 수도 있으나, 2개의 방사선 이미지(10,20) 중 어느 하나의 방사선 이미지(10 or 20)에서 대상 픽셀(111a)의 픽셀값은 정상값일 수도 있으므로, 2개의 방사선 이미지(10,20) 중 정상값인 대상 픽셀(111a)의 픽셀값을 걸러내어 방사선 이미지에 계속 사용되도록 할 수 있다. 이를 통해 방사선 이미지를 노이즈 및/또는 보정 없이 방사선에 의해 발생된 전하로만 생성된 최적의 방사선 이미지와 유사하게 영상처리할 수 있다. 랜덤 노이즈는 픽셀값을 상승시키게 되므로, 상기 랜덤 노이즈가 포함된 픽셀값은 더 높아질 수 밖에 없고, 대상 픽셀(111a)의 픽셀값이 상대적으로 높은 방사선 이미지(10 or 20)에서 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단할 수 있다.In the step (b), in step S200, two consecutive radiographic images 10 and 20 among the plurality of radiographic images may be selected, and in step (e), in step S500, at the time of the target pixel 111a. The pixel value of the target pixel 111a in the radiographic image 10 or 20 where the pixel value of the target pixel 111a is relatively high among the two radiographic images 10 and 20 in which the amount of change in the pixel value is greater than or equal to the preset reference value R Can be determined as an abnormal value. In general, the pixel values of the target pixel 111a may be determined as both abnormal values in the two radiographic images 10 and 20, but any one of the two radiographic images 10 and 20 may be a radiographic image 10 or 20. Since the pixel value of the target pixel 111a may be a normal value, the pixel value of the target pixel 111a, which is the normal value among the two radiographic images 10 and 20, may be filtered to be continuously used in the radiographic image. Through this, the radiographic image can be imaged similarly to the optimal radiographic image generated only with charge generated by radiation without noise and/or correction. Since the random noise increases the pixel value, the pixel value including the random noise is forced to be higher, and the target pixel 111a in the radiographic image 10 or 20 where the pixel value of the target pixel 111a is relatively high. The pixel value of can be determined as an abnormal value.
미리 설정된 기준값(R)은 복수의 픽셀(111) 전체적으로 각각의 픽셀(111)에서 방사선이 비조사되는 경우에도 발생되거나, 방사선으로 인해 발생되는 고유 노이즈값보다 클(또는 이상일) 수 있다. 미리 설정된 기준값(R)은 모든 복수의 픽셀(111)에 공통적으로 발생하는 상기 고유 노이즈값보다 크게 설정될 수 있다. 상기 고유 노이즈(intrinsic noise)는 그 값은 다를 수 있으나 모든 복수의 픽셀(111)에 공통적으로 발생하는 노이즈를 말하며, 장치적 오프셋(offset)에 의해 발생하는 노이즈일 수 있고, 이러한 노이즈는 모든 복수의 픽셀(111)에 공통적으로 발생하기 때문에 방사선 이미지를 이루는 복수의 픽셀(111)의 픽셀값들에서 장치적 오프셋값을 모두 제거하게 되면 큰 문제가 되지 않는다. 예를 들어, 상기 고유 노이즈는 열적(thermal) 원인, 절연성 불량 등의 원인에 의해 방사선에 기인하지 않고 발생하여 방사선이 비조사되는 경우에도 복수의 픽셀(111)에 전하를 축적시키는 암전류(dark current) 및 광자가 가진 양자적 성질로 말미암아 생기며, 방사선 발생장치(Radiation Generator) 등의 방사선원(Radiation Source)으로 인해 발생하는 양자 잡음(quantum noise)을 포함할 수 있다.The preset reference value R may be generated even when radiation is not irradiated from each pixel 111 as a whole of the plurality of pixels 111 or may be greater (or more) than the inherent noise value generated by radiation. The preset reference value R may be set larger than the intrinsic noise value common to all the plurality of pixels 111. The intrinsic noise may have a different value, but refers to noise common to all of the plurality of pixels 111, and may be noise generated by a device offset. Since it occurs in common with the pixels 111 of, it is not a big problem if all the device offset values are removed from the pixel values of the plurality of pixels 111 constituting the radiation image. For example, the intrinsic noise is generated without causing radiation due to a cause such as a thermal cause, poor insulation, etc., and causes dark current to accumulate charge in a plurality of pixels 111 even when radiation is not irradiated. ) And photons, and may include quantum noise generated by a radiation source such as a radiation generator.
여기서, 상기 고유 노이즈값은 각각 픽셀(111)마다 다를 수 있으므로, 미리 설정된 기준값(R)은 복수의 픽셀(111)의 고유 노이즈값의 평균값 또는 복수의 픽셀(111)의 고유 노이즈값 중 중간값이나 최대값보다 크게(또는 이상으로) 설정될 수 있다. 한편, 상기 랜덤 노이즈값은 상기 고유 노이즈값과 차이가 많이 나므로, 미리 설정된 기준값(R)은 복수의 픽셀(111)의 고유 노이즈값의 평균값의 약 2배로 설정할 수도 있고, 복수의 픽셀(111)의 고유 노이즈값 중 최소값, 중간값 또는 최대값의 약 2배로 설정할 수도 있다.Here, since the intrinsic noise values may be different for each pixel 111, the preset reference value R is an average value of the intrinsic noise values of the plurality of pixels 111 or an intermediate value among the intrinsic noise values of the plurality of pixels 111. Alternatively, it can be set larger (or higher) than the maximum value. On the other hand, since the random noise value is much different from the intrinsic noise value, the preset reference value R may be set to about twice the average value of the intrinsic noise values of the plurality of pixels 111, or the plurality of pixels 111 It can be set to about twice the minimum, median or maximum value of the intrinsic noise value of
본 발명에서는 미리 설정된 기준값(R)을 복수의 픽셀(111) 전체적으로 각각의 픽셀(111)에서 발생하는 상기 고유 노이즈값보다 크게 설정함으로써, 상기 선택된 둘 이상의 방사선 이미지(10,20) 간의 복수의 픽셀(111)의 시간에 따른 픽셀값 변화량들에서 상기 고유 노이즈에 의한 자연적인 변화량을 배제할 수 있고, 보다 효과적으로 방사선 이미지를 이루는 복수의 픽셀(111)의 픽셀값들에서 상기 랜덤 노이즈를 검출할 수 있다.In the present invention, the preset reference value R is set to be larger than the intrinsic noise value generated in each pixel 111 as a whole of the plurality of pixels 111, thereby thereby forming a plurality of pixels between the selected two or more radiographic images 10 and 20. The natural variation due to the intrinsic noise can be excluded from the variation in pixel values over time in (111), and the random noise can be detected from the pixel values of the plurality of pixels 111 constituting the radiation image more effectively. have.
한편, 미리 설정된 기준값(R)은 각각의 픽셀(111)에서 발생 가능한 상기 고유 노이즈의 최대값과 최소값의 차이값을 이용하여 설정할 수도 있으며, 상기 고유 노이즈의 최대값과 최소값의 차이값보다 크게(이상으로) 설정할 수 있다. 이때, 각각의 픽셀(111)에 대해 상기 고유 노이즈의 최대값과 최소값의 차이값을 산출하여 산출된 차이값들의 최소값, 중간값, 최대값 또는 평균값 이상으로 미리 설정된 기준값(R)을 설정할 수도 있고, 복수의 픽셀(111) 전체에서 상기 고유 노이즈의 최대값과 상기 고유 노이즈의 최소값의 차이값 이상으로 미리 설정된 기준값(R)을 설정할 수도 있다. 또한, 상기 고유 노이즈의 최대값과 최소값의 차이값의 1/2을 기준으로 삼을 수도 있다.Meanwhile, the preset reference value R may be set by using a difference value between the maximum and minimum values of the intrinsic noise that may occur in each pixel 111 and is greater than the difference between the maximum and minimum values of the intrinsic noise ( Above). At this time, for each pixel 111, the difference value between the maximum value and the minimum value of the intrinsic noise may be calculated to set a preset reference value R above a minimum value, a median value, a maximum value, or an average value of the calculated difference values. The preset reference value R may be set to be equal to or greater than a difference value between the maximum value of the intrinsic noise and the minimum value of the intrinsic noise in the entire plurality of pixels 111. In addition, it may be based on 1/2 of the difference between the maximum and minimum values of the inherent noise.
도 2는 본 발명의 일실시예에 따른 대상 픽셀과 주변 픽셀을 설명하기 위한 개념도이며, 도 3은 본 발명의 일실시예에 따른 미리 설정된 기준값을 통한 대상 픽셀의 픽셀값 검증을 설명하기 위한 개념도로, 도 3의 (a)는 비정상 픽셀값으로 결정한 경우를 나타내고, 도 3의 (b)는 정상 픽셀값으로 결정한 경우를 나타낸다.2 is a conceptual diagram for explaining a target pixel and a surrounding pixel according to an embodiment of the present invention, and FIG. 3 is a conceptual diagram for explaining a pixel value verification of a target pixel through a preset reference value according to an embodiment of the present invention 3(a) shows a case determined by an abnormal pixel value, and FIG. 3(b) shows a case determined by a normal pixel value.
도 2 및 도 3을 참조하면, 본 발명에 따른 방사선 영상 처리방법은 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 검증하는 과정(S600);을 더 포함할 수 있다.2 and 3, the radiographic image processing method according to the present invention may further include a process of verifying a pixel value of the target pixel 111a determined as an abnormal value (S600 ).
그리고 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 검증할 수 있다(S600). 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 검증하여 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 원인이 입력 신호의 변화에 기인한 것인지, 상기 랜덤 노이즈에 기인한 것인지 아님 둘 다에 기인한 것인지를 확인할 수 있고, 입력 신호의 변화만이 기인한 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값은 정상값으로 결정하여 방사선 이미지를 이루는 픽셀값으로 그대로 사용할 수 있다. 이를 통해 상기 랜덤 노이즈가 포함(또는 기인)된 비정상 픽셀값만을 보정하여 방사선 이미지의 왜곡을 방지할 수 있고, 방사선 이미지의 품질을 향상시킬 수 있다.In addition, the pixel value of the target pixel 111a determined to be an abnormal value may be verified (S600 ). By verifying the pixel value of the target pixel 111a determined as the abnormal value, whether the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to a preset reference value R is due to a change in the input signal, or the random It is possible to check whether it is caused by noise or both, and the pixel value of the target pixel 111a determined as an abnormal value due to only a change in the input signal is determined as a normal value, thereby forming a radiation image. Can be used as is. Through this, only the abnormal pixel value including (or attributable to) the random noise may be corrected to prevent distortion of the radiographic image, and quality of the radiographic image may be improved.
즉, 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값(R)을 비교하여 비정상 픽셀값을 판단한 후에 다시 주변 픽셀(111b)을 통해 비정상 픽셀값이 맞는지를 검증함으로써, 피사체의 움직임 등에 의한 입력 신호의 변화를 효과적으로 반영하여 각 대상 픽셀(111a)의 픽셀값이 상기 랜덤 노이즈를 포함하는 비정상 픽셀값인지를 보다 정확하게 판단할 수 있다.That is, after comparing the amount of change of the pixel value over time of the target pixel 111a with a preset reference value R to determine the abnormal pixel value, and verifying that the abnormal pixel value is correct through the surrounding pixels 111b, the subject moves. It is possible to more accurately determine whether the pixel value of each target pixel 111a is an abnormal pixel value including the random noise by effectively reflecting the change of the input signal due to the like.
상기 대상 픽셀(111a)의 픽셀값을 검증하는 과정(S600)은 상기 선택된 둘 이상의 방사선 이미지(10,20)로부터 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값(R)을 비교하는 과정(S610); 및 상기 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n1) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 확인하는 과정(S620)을 포함할 수 있다.In the process of verifying the pixel value of the target pixel 111a (S600), the amount of change in the pixel value over time of the surrounding pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20 Comparing the preset reference value R with a step S610; And the pixel value of the target pixel 111a determined as the abnormal value is abnormal when the number of surrounding pixels 111b that is greater than or equal to the preset reference value R is less than or equal to the preset determination number n1. It may include a process of checking by value (S620).
상기 선택된 둘 이상의 방사선 이미지(10,20)로부터 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값(R)을 비교할 수 있다(S610). 대상 픽셀(111a)과 동일하게 상기 선택된 둘 이상의 방사선 이미지(10,20)로 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량도 미리 설정된 기준값(R)을 비교할 수 있다. 이를 통해 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량이 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 유사한지 경향성을 파악할 수 있다.The amount of change in pixel value over time of the neighboring pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20 and a preset reference value R may be compared (S610 ). Similar to the target pixel 111a, the amount of change in the pixel value over time of the adjacent pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10, 20 is also compared with a preset reference value R Can. Through this, it is possible to grasp the tendency of the amount of change in the pixel value over time in the peripheral area of the target pixel 111a to be similar to the amount of change in the pixel value over time of the target pixel 111a.
여기서, 대상 픽셀(111a)은 도면에서 a로 나타낼 수도 있으며, 주변 픽셀(111b)은 도 2 및 도 3과 같이 복수개일 수 있고, b1 ~ b8로 나타낼 수 있다. 그리고 상기 선택된 둘 이상의 방사선 이미지(10,20)는 제1 방사선 이미지(10)와 제2 방사선 이미지(20)를 포함할 수 있으며, 제1 방사선 이미지(10)는 제1 시간(t1)에 생성된 방사선 이미지일 수 있고, 제2 방사선 이미지(20)는 제2 시간(t2)에 생성된 방사선 이미지일 수 있다. 이때, 제1 방사선 이미지(10)에서는 제1 시간(t1)을 반영하여 도 2의 (b)와 같이 대상 픽셀(111a)을 a1으로 나타낼 수 있고, 주변 픽셀(111b)을 b11 ~ b81로 나타낼 수 있으며, 제2 방사선 이미지(20)에서는 제2 시간(t2)을 반영하여 도 2의 (c)와 같이 대상 픽셀(111a)을 a2으로 나타낼 수 있고, 주변 픽셀(111b)을 b12 ~ b82로 나타낼 수 있다.Here, the target pixel 111a may be represented by a in the drawing, and a plurality of peripheral pixels 111b may be represented by b1 to b8 as shown in FIGS. 2 and 3. Further, the selected two or more radiographic images 10 and 20 may include a first radiographic image 10 and a second radiographic image 20, and the first radiographic image 10 is generated at a first time t1. It may be a radiographic image, and the second radiographic image 20 may be a radiographic image generated at a second time t2. In this case, in the first radiographic image 10, the target pixel 111a may be represented as a1 as shown in FIG. 2B by reflecting the first time t1, and the peripheral pixels 111b may be represented as b11 to b81. In the second radiographic image 20, the target pixel 111a may be represented by a2 as shown in FIG. 2(c) by reflecting the second time t2, and the peripheral pixels 111b may be expressed by b12 to b82. Can be represented.
그리고 상기 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n1) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 확인(또는 확정)할 수 있다(S620). 도 3의 (a)와 같이, 상기 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n1) 이하인 경우에는 대상 픽셀(111a)과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량이 유사한 경향성을 확인할 수 없어 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 원인이 상기 랜덤 노이즈에 기인한다고 볼 수 있고, 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 결정(또는 확정)할 수 있다. 그리고 도 3의 (b)와 같이, 상기 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n1)보다 커서 개수가 많은 경우에는 대상 픽셀(111a)과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량이 유사한 경향성을 갖는다고 볼 수 있고, 상기 랜덤 노이즈보다는 입력 신호에 큰 변화가 있다고 간주할 수 있다. 이에 따라 입력 신호의 변화에 기인한 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 판단의 번복을 통해 정상값으로 결정할 수 있으며, 이를 통해 입력 신호의 변화에 기인한 정상 픽셀값을 비정상값으로의 판단에서 배제시킬 수 있다.The pixel value of the target pixel 111a determined as the abnormal value is abnormal when the number of surrounding pixels 111b that is greater than or equal to the preset reference value R is less than or equal to the preset determination number n1. It can be confirmed (or confirmed) by the value (S620). As shown in (a) of FIG. 3, when the number of surrounding pixels 111b that is greater than or equal to the preset reference value R is less than or equal to the preset determination number n1, the target pixel 111a and the target pixel The reason why the amount of change in the pixel value over time in the peripheral region of (111a) cannot be confirmed is similar, and the reason that the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R can be considered to be attributed to the random noise. In addition, the pixel value of the target pixel 111a determined as the abnormal value may be determined (or determined) as the abnormal value. And, as shown in (b) of FIG. 3, when the number of surrounding pixels 111b having a change amount of the pixel value over the preset reference value R is greater than a preset determination number n1, the target pixel ( It can be considered that the amount of change in the pixel value over time in the peripheral region of 111a) and the target pixel 111a has a similar tendency, and it can be considered that there is a large change in the input signal rather than the random noise. Accordingly, the pixel value of the target pixel 111a, which is determined to be an abnormal value due to a change in the input signal, can be determined as a normal value through the overturning of the determination, thereby allowing the normal pixel value due to a change in the input signal to be an abnormal value. It can be excluded from judgment.
여기서, 상기 판단 개수(n1)는 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량과 유사한지 경향성을 파악할 수 있는 개수로 알맞게 설정될 수 있으며, 유사한 경향성을 확실하게 파악할 수 있도록 3개 이상으로 설정될 수 있다.Here, the number of determinations (n1) is appropriately set to a number that can determine the tendency whether the amount of change in the pixel value over time of the target pixel 111a is similar to the amount of change in the pixel value over time in the peripheral area of the target pixel 111a. It can be set to three or more to reliably grasp similar trends.
도 4는 본 발명의 일실시예에 따른 대상 픽셀의 시간에 따른 픽셀값 변화량과 주변 픽셀의 시간에 따른 픽셀값 변화량의 비교를 통한 대상 픽셀의 픽셀값 검증을 설명하기 위한 개념도이다.4 is a conceptual diagram illustrating verification of a pixel value of a target pixel through comparison of a pixel value change amount with time of a target pixel and a change amount of pixel value with time of a surrounding pixel according to an embodiment of the present invention.
도 4를 참조하면, 상기 대상 픽셀(111a)의 픽셀값을 검증하는 과정(S600)은 상기 선택된 둘 이상의 방사선 이미지(10,20)로부터 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 비교하는 과정(S630); 및 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차(d) 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n2) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 확인하는 과정(S640)을 포함할 수 있다.Referring to FIG. 4, a process (S600) of verifying a pixel value of the target pixel 111a is performed on a neighboring pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20. Comparing the amount of change in the pixel value over time and the amount of change in the pixel value over time of the target pixel (111a) (S630); And the number of surrounding pixels 111b having the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel within the preset tolerance (d) range is less than or equal to the preset number of determinations (n2). And checking the pixel value of the target pixel 111a determined as an abnormal value (S640 ).
상기 선택된 둘 이상의 방사선 이미지(10,20)로부터 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 비교할 수 있다(S630). 이를 통해 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량이 유사한지 경향성을 파악할 수 있다.A pixel value change amount over time of a neighboring pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20 and a pixel value change amount over time of the target pixel 111a may be compared. (S630). Through this, it is possible to grasp the tendency of the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel value over time in the peripheral area of the target pixel 111a.
그리고 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차(d) 범위(± d) 이내인 픽셀값 변화량을 갖는 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n2) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 확인(또는 확정)할 수 있다(S640). 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차(d) 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n2) 이하인 경우에는 대상 픽셀(111a)과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량이 유사한 경향성을 확인할 수 없어 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 원인이 상기 랜덤 노이즈에 기인한다고 볼 수 있고, 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 결정할 수 있다. 그리고 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차(d) 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n2)보다 커서 개수가 많은 경우에는 대상 픽셀(111a)과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량이 유사한 경향성을 갖는다고 볼 수 있고, 상기 랜덤 노이즈보다는 입력 신호에 큰 변화가 있다고 간주할 수 있다. 이에 따라 입력 신호의 변화에 기인한 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 판단의 번복을 통해 정상값으로 결정할 수 있으며, 이를 통해 입력 신호의 변화에 기인한 정상 픽셀값을 비정상값으로의 판단에서 배제시킬 수 있다. 이때, 상기 허용편차(d)는 미리 설정된 기준값(R)에서 알맞게 설정될 수 있고, 그리 크지 않은 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량의 약 10 % 이내로 알맞게 정해질 수 있다.Also, when the number of surrounding pixels 111b having the amount of change in the pixel value and the amount of change in the pixel within the preset tolerance (d) range (± d) according to the time of the target pixel 111a is less than or equal to the preset determination number (n2) The pixel value of the target pixel 111a determined as the abnormal value may be checked (or confirmed) as the abnormal value (S640). The target pixel 111a when the number of pixel values variation over time of the target pixel 111a and the number of surrounding pixels 111b having a pixel value variation within a preset allowable deviation (d) range is less than or equal to a preset determination number (n2) ) And the amount of change in the pixel value over time in the peripheral area of the target pixel 111a cannot be confirmed, and the reason for the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to the preset reference value R is caused by the random noise. It can be considered to be due to the abnormality, and the pixel value of the target pixel 111a determined as the abnormal value may be determined as the abnormal value. Also, when the number of surrounding pixels 111b having the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel within the preset tolerance (d) range is larger than the number of predetermined determinations (n2). It can be considered that the amount of change in the pixel value over time in the peripheral region of the target pixel 111a and the target pixel 111a has a similar tendency, and it can be considered that there is a large change in the input signal rather than the random noise. Accordingly, the pixel value of the target pixel 111a, which is determined to be an abnormal value due to a change in the input signal, can be determined as a normal value through the overturning of the determination, thereby allowing the normal pixel value due to a change in the input signal to be an abnormal value. It can be excluded from judgment. At this time, the allowable deviation (d) may be appropriately set at a preset reference value (R), and may be appropriately determined within about 10% of a change in the pixel value over time of the target pixel (111a), which is not very large.
여기서, 상기 판단 개수(n2)는 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량과 유사한지 경향성을 파악할 수 있는 개수로 알맞게 설정될 수 있으며, 유사한 경향성을 확실하게 파악할 수 있도록 3개 이상으로 설정될 수 있다.Here, the number of determinations n2 is suitably set to a number capable of determining the tendency of the amount of change in the pixel value over time of the target pixel 111a to be similar to the amount of change in the pixel value over time in the peripheral area of the target pixel 111a. It can be set to three or more to reliably grasp similar trends.
또한, 도 4의 경우에는 상기 판단 개수(n1,n2)가 3개라고 할 때에 상기 선택된 둘 이상의 방사선 이미지(10,20)로 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값(R)을 비교하게 되면, 상기 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 주변 픽셀(111b)의 개수가 2개 밖에 되지 않아 대상 픽셀(111a)의 픽셀값을 비정상값으로 결정하게 된다. 하지만, 상기 선택된 둘 이상의 방사선 이미지(10,20)로 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 비교하게 되면, 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차(d) 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀(111b)의 개수가 8개로 대상 픽셀(111a)의 픽셀값을 정상값으로 결정하게 된다. 도 4를 보면, a와 b1~b8이 큰 편차 없이 유사한 경향성을 갖는 것으로 보이므로, 상기 선택된 둘 이상의 방사선 이미지(10,20)로 산출된 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 비교하는 경우에는 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량의 경향성을 보다 정확하게 파악(또는 확인)할 수 있다.In addition, in the case of FIG. 4, when the number of determinations (n1, n2) is 3, at the time of the surrounding pixels 111b adjacent to the target pixel 111a calculated as the selected two or more radiographic images 10 and 20, When the amount of change of the pixel value according to the preset reference value R is compared, the number of surrounding pixels 111b having the amount of change of the pixel value over the preset reference value R is only two, so that the target pixel 111a The pixel value is determined as an abnormal value. However, the amount of change in pixel value over time of the adjacent pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20 is compared with the amount of change in pixel value over time of the target pixel 111a. When doing so, the number of surrounding pixels 111b having the amount of change in pixel value over time of the target pixel 111a and the amount of change in pixel value within a preset tolerance (d) range is 8, and the pixel value of the target pixel 111a is set. It is decided by the normal value. Referring to FIG. 4, since a and b1 to b8 appear to have similar tendencies without large deviation, the time of the peripheral pixel 111b adjacent to the target pixel 111a calculated from the selected two or more radiographic images 10 and 20 When comparing the amount of change of the pixel value with time and the amount of change of the pixel value with time of the target pixel 111a, the amount of change with the pixel value with time of the target pixel 111a and the value of the pixel with time of the peripheral area of the target pixel 111a It is possible to grasp (or confirm) the tendency of the change amount more accurately.
이와 같이, 본 발명에서는 시간에 따라 획득된 복수의 방사선 이미지 중 둘 이상의 방사선 이미지(10,20)를 선택하여, 선택된 둘 이상의 방사선 이미지(10,20) 간의 픽셀값 변화량을 고려함으로써, 시간에 따른 입력 신호의 변화를 반영하여 각 대상 픽셀(111a)의 픽셀값이 상기 랜덤 노이즈를 포함하는 비정상 픽셀값인지를 판단할 수 있다. 이에 따라 복수의 방사선 이미지에서 상기 랜덤 노이즈가 포함된 비정상 픽셀값을 제거하여 보정함으로써, 방사선 이미지의 왜곡을 방지할 수 있고, 방사선 이미지의 품질을 향상시킬 수 있다.As described above, in the present invention, by selecting two or more radiographic images 10 and 20 among a plurality of radiographic images acquired over time, by considering the amount of change in pixel values between the selected two or more radiographic images 10 and 20, over time, In response to a change in the input signal, it may be determined whether the pixel value of each target pixel 111a is an abnormal pixel value including the random noise. Accordingly, by correcting by removing the abnormal pixel value including the random noise from a plurality of radiographic images, distortion of the radiographic image can be prevented, and the quality of the radiographic image can be improved.
상기 복수의 방사선 이미지 전체에 대해 상기 (b) 과정(S200) 내지 상기 (e) 과정(S500)을 반복할 수 있으며, 본 발명에 따른 방사선 영상 처리방법은 비정상값으로 판단된 픽셀값이 발생한 대상 픽셀(111a)의 위치를 저장하는 과정(S700);을 더 포함할 수 있고, 상기 대상 픽셀(111a)의 위치를 저장하는 과정(S700) 이후의 (b) 과정(S200)에서 위치가 저장된 픽셀(111)을 대상 픽셀(111a)의 선택에서 제외하고 상기 위치가 저장된 픽셀(111)의 픽셀값을 비정상값으로 판단할 수 있다.The (b) process (S200) to the (e) process (S500) may be repeated for the whole of the plurality of radiographic images, and the radiographic image processing method according to the present invention may generate a pixel value determined as an abnormal value. The process of storing the position of the pixel 111a (S700); may further include, and the pixel in which the position is stored in the (B) process (S200) after the process of storing the position of the target pixel 111a (S700) The pixel value of the pixel 111 in which the location is stored may be determined as an abnormal value by excluding (111) from the selection of the target pixel 111a.
상기 (b) 과정(S200) 내지 상기 (e) 과정(S500)을 반복할 수 있으며, 상기 복수의 방사선 이미지 전체에 대해 상기 (b) 과정(S200) 내지 상기 (e) 과정(S500)을 수행할 수 있다. 또한, 모든 복수의 픽셀(111) 또는 검사(또는 판단)을 원하는 모든 대상 픽셀(111a)에 대해 픽셀값이 비정상값인지를 판단할 수 있다.The (b) process (S200) to the (e) process (S500) may be repeated, and the (b) process (S200) to (e) process (S500) may be performed on the whole of the plurality of radiographic images. can do. Further, it may be determined whether the pixel value is an abnormal value for all the plurality of pixels 111 or all target pixels 111a for which inspection (or determination) is desired.
그리고 비정상값으로 판단(또는 결정)된 픽셀값이 발생한 대상 픽셀(111a)의 위치를 저장할 수 있다(S700). 상기 랜덤 노이즈의 특성상 한 번 상기 랜덤 노이즈가 발생했던 픽셀(111)에서는 상기 랜덤 노이즈가 다시 발생할 확률이 높기 때문에 비정상 픽셀값이 발생한 대상 픽셀(111a)의 위치를 저장하고 다음 이후의 방사선 이미지에서는 위치가 저장된 픽셀(111)에 대해 비정상 픽셀값인지 판단하는 상기 (b) 과정(S200) 내지 상기 (e) 과정(S500)의 수행 없이 비정상 픽셀값으로 간주할 수 있다.In addition, the position of the target pixel 111a in which the pixel value determined (or determined) as the abnormal value has occurred may be stored (S700 ). Due to the nature of the random noise, since the probability of the random noise occurring again is high in the pixel 111 where the random noise has once occurred, the position of the target pixel 111a in which an abnormal pixel value has occurred is stored, and the position in the subsequent radiation image It may be regarded as an abnormal pixel value without performing the steps (b) (S200) to (e) (S500) of determining whether the stored pixel 111 is an abnormal pixel value.
상기 대상 픽셀(111a)의 위치를 저장하는 과정(S700) 이후의 (b) 과정(S200)에서는 위치가 저장된 픽셀(111)을 대상 픽셀(111a)의 선택에서 제외하고, 상기 위치가 저장된 픽셀(111)의 픽셀값을 비정상값으로 판단(또는 간주)할 수 있다. 즉, 상기 랜덤 노이즈의 특성상 한 번 상기 랜덤 노이즈가 발생했던 픽셀(111)에서는 상기 랜덤 노이즈가 발생할 확률이 높기 때문에 비정상 픽셀값이 발생한 대상 픽셀(111a)의 위치를 저장하고, 이후의 방사선 이미지에서 해당 픽셀(111)의 픽셀값을 비정상 픽셀값으로 간주함으로써, 비정상 픽셀값인지를 판단해야 하는 픽셀(111)의 개수가 줄어들 수 있다. 또한, 상기 복수의 방사선 이미지 전체에 대해 비정상 픽셀값이 있는지를 판단하고 보정하는 영상 처리의 시간이 단축될 수 있다.In the step (B) after the process of storing the position of the target pixel 111a (S700) (S200), the pixel storing the position is excluded from the selection of the target pixel 111a, and the pixel storing the position ( The pixel value of 111) may be determined (or considered) as an abnormal value. That is, because of the nature of the random noise, since the probability of generating the random noise is high in the pixel 111 where the random noise has once occurred, the location of the target pixel 111a where the abnormal pixel value has occurred is stored in the subsequent radiation image. By considering the pixel value of the corresponding pixel 111 as an abnormal pixel value, the number of pixels 111 to determine whether the pixel value is abnormal can be reduced. In addition, the time of image processing for determining and correcting whether there are abnormal pixel values for all of the plurality of radiographic images may be shortened.
본 발명에 따른 방사선 영상 처리방법은 복수의 픽셀(111) 중 미리 설정된 범위 내의 세기로 방사선이 입사되는 픽셀 영역을 검출하는 과정(S150);을 더 포함할 수 있고, 상기 (b) 과정(S200)에서는 상기 픽셀 영역 중 대상 픽셀(111a)을 선택할 수 있다.The radiographic image processing method according to the present invention may further include a step (S150) of detecting a pixel region to which radiation is incident at an intensity within a predetermined range among a plurality of pixels 111 (S200). ), a target pixel 111a may be selected from the pixel areas.
그리고 복수의 픽셀(111) 중 미리 설정된 범위 내의 세기로 방사선이 입사되는 픽셀 영역을 검출할 수 있다(S150). 방사선 이미지에는 피사체의 경계(edge) 영역, 피사체의 각 부분 간의 경계 영역 등 인접 픽셀(111) 간의 픽셀값 차이가 심한(또는 매우 큰) 복잡(complicated) 영역과 피사체의 각 부분의 중앙 영역 등 인접 픽셀(111) 간의 픽셀값 차이가 크지 않은(또는 미미한) 단순(simple) 영역이 존재할 수 있다. 상기 복잡 영역은 각 픽셀(111)에 들어가는 입력 신호가 거의 모두 다르기 때문에 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량의 경향성을 파악할 수가 없어 상기 랜덤 노이즈가 포함된 비정상 픽셀값을 검출하기 어려우며, 상기 랜덤 노이즈가 포함된다 하더라도 방사선 이미지에서 눈으로 식별 불가한 정도이기 때문에 방사선 이미지로 의학적 판단을 하는 부분에 큰 문제가 되지 않는다. 하지만, 상기 단순 영역은 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 대상 픽셀(111a)의 주변 영역의 시간에 따른 픽셀값 변화량의 경향성 파악이 가능하기 때문에 상기 랜덤 노이즈가 포함된 비정상 픽셀값을 검출할 수 있고, 상기 랜덤 노이즈가 포함되게 되면, 방사선 이미지에서 이상 부분으로 식별되어 상기 랜덤 노이즈가 포함된 방사선 이미지로 의학적 판단을 하는 경우에 큰 문제가 될 수 있다. 이에 복수의 픽셀(111) 중 미리 설정된 범위 내의 세기로 방사선이 입사되는 픽셀 영역(즉, 상기 단순 영역)을 검출하여 상기 미리 설정된 범위 내의 세기로 방사선이 입사되는 픽셀 영역의 픽셀(들)만 픽셀값이 비정상값인지를 판단할 수 있다.In addition, a pixel region in which radiation is incident at an intensity within a predetermined range among the plurality of pixels 111 may be detected (S150 ). In the radiographic image, a complex region in which the pixel value difference between adjacent pixels 111 is severe (or very large), such as an edge region of the subject and a boundary region between parts of the subject, and a central region of each part of the subject, such as A simple region in which the difference in pixel values between the pixels 111 is not large (or insignificant) may exist. Since the input signal to each pixel 111 is almost entirely different in the complex area, the tendency of the pixel value change over time of the target pixel 111a and the pixel value change over time of the peripheral area of the target pixel 111a can be grasped. It is difficult to detect the abnormal pixel value including the random noise, and even if the random noise is included, it is not a big problem in the part where medical judgment is made from the radiographic image because it is unrecognizable to the eye in the radiographic image. However, since the simple region can grasp the tendency of the amount of change in the pixel value over time of the target pixel 111a and the amount of change in the pixel value over time in the peripheral area of the target pixel 111a, the abnormal pixel value including the random noise When it is possible to detect and the random noise is included, it may be a problem when the medical image is identified as the radiographic image containing the random noise, which is identified as an abnormal part in the radiographic image. Accordingly, a pixel region (ie, the simple region) in which radiation is incident at an intensity within a predetermined range among the plurality of pixels 111 is detected, and only the pixel(s) of the pixel region in which radiation is incident at an intensity within the predetermined range is pixels. It can be determined whether the value is an abnormal value.
상기 (b) 과정(S200)에서는 상기 픽셀 영역 중 대상 픽셀(111a)을 선택할 수 있다. 상기 랜덤 노이즈가 포함된 비정상 픽셀값을 검출 가능한 상기 미리 설정된 범위 내의 세기로 방사선이 입사되는 픽셀 영역의 픽셀(들)만 픽셀값이 비정상값인지를 판단하기 위해 상기 픽셀 영역 중 대상 픽셀(111a)을 선택할 수 있다. 또한, 상기 픽셀 영역 중 대상 픽셀(111a)을 선택하게 되면, 복수의 픽셀(111)을 모두 픽셀값이 비정상값인지 판단하지 않고 일부의 상기 픽셀 영역에 해당하는 픽셀(들)만 픽셀값이 비정상값인지 판단함으로써, 선택된 방사선 이미지에서 상기 랜덤 노이즈가 포함된 비정상 픽셀값을 검출하는 시간이 단축될 수 있다.In the step (b) (S200), a target pixel 111a may be selected from the pixel areas. Target pixel 111a among the pixel areas to determine whether only the pixel(s) of the pixel area where radiation is incident at an intensity within the preset range capable of detecting the abnormal pixel value including the random noise is an abnormal value You can choose In addition, when the target pixel 111a is selected from the pixel areas, all of the plurality of pixels 111 do not determine whether the pixel values are abnormal values, but only the pixel(s) corresponding to some of the pixel areas have abnormal pixel values. By determining whether it is a value, the time for detecting the abnormal pixel value including the random noise in the selected radiation image can be shortened.
본 발명에 따른 방사선 영상 처리방법은 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 동일한 방사선 이미지 상의 인접한 주변 픽셀(111b)의 픽셀값을 내삽하여 보정하는 과정(S800);을 더 포함할 수 있다.The method of processing a radiographic image according to the present invention further includes a process of correcting a pixel value of a target pixel 111a determined as an abnormal value by interpolating a pixel value of adjacent peripheral pixels 111b on the same radiographic image (S800). Can.
그리고 비정상값으로 판단(또는 결정)된 대상 픽셀(111a)의 픽셀값을 인접한 동일(한) 방사선 이미지 상의 주변 픽셀(111b)의 픽셀값을 내삽하여 보정할 수 있다(S800). 예를 들어, 비정상값으로 판단(또는 결정)된 대상 픽셀(111a)의 픽셀값을 인접한 주변 픽셀(111b)의 픽셀값의 평균값(mean value) 또는 중간값(median value)으로 치환할 수 있다. 이때, 인접한 주변 픽셀(111b)에 비정상값으로 판단(또는 결정)된 픽셀값을 갖는 픽셀(111)이 포함되어 있는 경우에는 평균값 또는 중간값의 계산에서 비정상값으로 판단된 픽셀값을 갖는 주변 픽셀(111b)을 제외할 수 있다. 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 인접한 주변 픽셀(111b)의 픽셀값을 보정하는 방법은 이러한 방법 외에도 다양하게 수행될 수 있다.Further, the pixel value of the target pixel 111a determined (or determined) as an abnormal value may be corrected by interpolating the pixel value of the adjacent pixels 111b on the adjacent (one) radiation image (S800). For example, the pixel value of the target pixel 111a determined (or determined) as an abnormal value may be replaced with an average value or a median value of the pixel values of the adjacent peripheral pixels 111b. At this time, when the adjacent peripheral pixel 111b includes a pixel 111 having a pixel value determined (or determined) as an abnormal value, a peripheral pixel having a pixel value determined as an abnormal value in calculation of an average value or a median value (111b) can be excluded. The method of correcting the pixel value of the adjacent pixel 111b with the pixel value of the target pixel 111a determined to be an abnormal value may be performed in various ways in addition to this method.
도 5는 본 발명의 일실시예에 따른 방사선 촬영장치의 구성을 나타낸 블록도이다.5 is a block diagram showing the configuration of a radiographic imaging apparatus according to an embodiment of the present invention.
도 5를 참조하여 본 발명의 다른 실시예에 따른 방사선 촬영장치를 보다 상세히 살펴보는데, 본 발명의 일실시예에 따른 방사선 영상 처리방법과 관련하여 앞서 설명된 부분과 중복되는 사항들은 생략하도록 한다.The radiographic imaging apparatus according to another embodiment of the present invention will be described in more detail with reference to FIG. 5, and the overlapping with the above-described part in relation to the radiographic image processing method according to an embodiment of the present invention will be omitted.
본 발명의 다른 실시예에 따른 방사선 촬영장치(100)는 복수의 픽셀(111)을 구비하며, 조사되는 방사선에 의해 시간에 따라 복수의 방사선 이미지를 획득하기 위한 전기신호를 생성하는 방사선 검출패널(110); 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지(10,20)에서 상기 복수의 픽셀(111) 중 선택된 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량 및 상기 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량을 산출하는 변화량 산출부(121); 산출된 상기 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값(R)과 비교하는 변화량 비교부(122); 및 상기 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값(R) 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지(10,20) 중 적어도 하나의 방사선 이미지에서 상기 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단하는 픽셀값 이상판단부(123);를 포함할 수 있다.The radiographic imaging apparatus 100 according to another embodiment of the present invention includes a plurality of pixels 111, and a radiation detection panel for generating electrical signals for acquiring a plurality of radiographic images over time by irradiated radiation ( 110); The amount of change in pixel value over time of the selected target pixel 111a among the plurality of pixels 111 from two or more selected radiographic images 10 and 20 among the plurality of radiographic images and surrounding pixels adjacent to the target pixel 111a ( A change amount calculating unit 121 for calculating a pixel value change amount over time in 111b); A change amount comparison unit 122 comparing the calculated change amount of the pixel value with time of the target pixel 111a with a preset reference value R; And when the amount of change in the pixel value over time of the target pixel 111a is equal to or greater than the preset reference value R, the target pixel 111a is selected from at least one of the selected two or more radiographic images 10 and 20. It may include; a pixel value abnormality determination unit 123 for determining the pixel value as an abnormal value.
방사선 검출패널(110)은 복수의 픽셀(111)을 구비할 수 있고, 조사되는 방사선에 의해 시간에 따라 복수의 방사선 이미지를 획득하기 위한 전기신호를 생성할 수 있다. 이렇게 생성된 전기신호에 의해 방사선 이미지가 생성되고 획득될 수 있다.The radiation detection panel 110 may include a plurality of pixels 111 and may generate electrical signals for acquiring a plurality of radiation images over time by irradiated radiation. A radiographic image may be generated and obtained by the electric signal thus generated.
변화량 산출부(121)는 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지(10,20)에서 상기 복수의 픽셀(111) 중 선택된 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량 및 상기 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량을 산출할 수 있다. 한편, 본 발명의 방사선 촬영장치(100)는 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지(10,20)에서 복수의 픽셀(111) 중 대상 픽셀(111a)을 선택하는 픽셀 설정부(미도시);를 더 포함할 수도 있다.The change amount calculating unit 121 may change the pixel value and the target pixel over time of the selected target pixel 111a among the plurality of pixels 111 from two or more selected radiation images 10 and 20 among the plurality of radiation images. The amount of change in the pixel value over time of the peripheral pixel 111b adjacent to 111a) may be calculated. Meanwhile, the radiographic imaging apparatus 100 of the present invention is a pixel setting unit (not shown) that selects a target pixel 111a among a plurality of pixels 111 from two or more radiographic images 10 and 20 selected from the plurality of radiographic images. ); may be further included.
변화량 비교부(122)는 산출된 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값(R)과 비교할 수 있으며, 대상 픽셀(111a)에 인접한 주변 픽셀(111b)의 시간에 따른 픽셀값 변화량을 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량 및/또는 미리 설정된 기준값(R)과 비교할 수도 있다.The change amount comparison unit 122 may compare the calculated change amount of the pixel value according to the time of the target pixel 111a with a preset reference value R, and the pixel according to the time of the peripheral pixel 111b adjacent to the target pixel 111a The value change amount may be compared with a pixel value change amount and/or a preset reference value R over time of the target pixel 111a.
픽셀값 이상판단부(123)는 상기 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지(10,20)에서 상기 대상 픽셀(111a)의 픽셀값을 비정상값으로 판단할 수 있다.The pixel value anomaly determining unit 123 may select the target pixel 111a from the selected two or more radiographic images 10 and 20 when the amount of change in the pixel value over time of the target pixel 111a is greater than or equal to a preset reference value R. The pixel value of can be determined as an abnormal value.
본 발명의 방사선 촬영장치(100)는 생성 또는 획득된 방사선 이미지의 영상 처리를 수행하는 영상처리부(120);를 더 포함할 수 있고, 영상처리부(120)는 변화량 산출부(121), 변화량 비교부(122) 및 픽셀값 이상판단부(123)를 포함할 수 있다. The radiographic imaging apparatus 100 of the present invention may further include an image processing unit 120 that performs image processing of the generated or obtained radiographic image, and the image processing unit 120 compares the change amount calculation unit 121 and the change amount The unit 122 and the pixel value abnormality determination unit 123 may be included.
본 발명에 따른 방사선 촬영장치(100)는 상기 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n1) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 결정하는 제1 픽셀값 검증부(미도시);를 더 포함할 수 있다.The radiographic imaging apparatus 100 according to the present invention is determined to be the abnormal value when the number of surrounding pixels 111b equal to or greater than the preset reference value R is less than or equal to the preset determination number n1. It may further include a first pixel value verification unit (not shown) that determines the pixel value of the target pixel 111a as an abnormal value.
제1 픽셀값 검증부(미도시)는 상기 시간에 따른 픽셀값 변화량이 미리 설정된 기준값(R) 이상인 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n1) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 결정하여 대상 픽셀(111a)의 픽셀값이 비정상값인지를 정확하게 검증할 수 있다.The first pixel value verification unit (not shown) determines the abnormal value when the number of surrounding pixels 111b equal to or greater than the preset reference value R is less than or equal to the preset determination number n1. The pixel value of the target pixel 111a may be determined as an abnormal value to accurately verify whether the pixel value of the target pixel 111a is an abnormal value.
본 발명에 따른 방사선 촬영장치(100)는 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차(d) 범위(± d) 이내인 픽셀값 변화량을 갖는 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n2) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 결정하는 제2 픽셀값 검증부(미도시);를 더 포함할 수 있다.The radiographic imaging apparatus 100 according to the present invention is the number of surrounding pixels 111b having a pixel value change amount over time and a predetermined value (d) within a range (± d) of the target pixel 111a. A second pixel value verification unit (not shown) that determines the pixel value of the target pixel 111a determined as the abnormal value as an abnormal value when is less than or equal to the preset number of determinations n2 may be further included.
제2 픽셀값 검증부(미도시)는 대상 픽셀(111a)의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차(d) 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀(111b)의 개수가 미리 설정된 판단 개수(n2) 이하인 경우에 상기 비정상값으로 판단된 대상 픽셀(111a)의 픽셀값을 비정상값으로 결정하여 대상 픽셀(111a)의 픽셀값이 비정상값인지를 정확하게 검증할 수 있다.The second pixel value verifying unit (not shown) determines that the number of surrounding pixels 111b having the amount of change in pixel values over time of the target pixel 111a and the amount of change in pixel values within a preset tolerance (d) range is preset. When the number n2 is less than or equal to the pixel value of the target pixel 111a determined as the abnormal value, the pixel value of the target pixel 111a can be accurately verified.
본 발명에 따른 방사선 촬영장치(100)는 비정상값으로 판단된 픽셀값이 발생한 대상 픽셀(111a)의 위치를 저장하는 위치저장부(미도시);를 더 포함할 수 있다.The radiographic imaging apparatus 100 according to the present invention may further include a location storage unit (not shown) that stores the location of the target pixel 111a where the pixel value determined as the abnormal value has occurred.
위치저장부(미도시)는 비정상값으로 판단된 픽셀값이 발생한 대상 픽셀(111a)의 위치를 저장할 수 있다. 랜덤 노이즈의 특성상 한 번 상기 랜덤 노이즈가 발생했던 픽셀(111)에서는 상기 랜덤 노이즈가 다시 발생할 확률이 높기 때문에 비정상 픽셀값이 발생한 대상 픽셀(111a)의 위치를 저장하고 다음 이후의 방사선 이미지에서는 위치가 저장된 픽셀(111)에 대해 비정상 픽셀값인지 판단하는 별도의 과정 없이 비정상 픽셀값으로 간주할 수 있다.The location storage unit (not shown) may store the location of the target pixel 111a where the pixel value determined as the abnormal value occurs. Due to the nature of random noise, since the probability of the random noise occurring again is high in the pixel 111 where the random noise has once occurred, the position of the target pixel 111a where the abnormal pixel value has occurred is stored, and the position in the subsequent radiation image is The stored pixel 111 may be regarded as an abnormal pixel value without a separate process of determining whether it is an abnormal pixel value.
본 발명에 따른 방사선 촬영장치(100)는 상기 복수의 픽셀 전체적으로 각각의 픽셀에서 방사선이 비조사되는 경우에도 발생되거나, 방사선원으로 인해 발생되는 고유 노이즈값보다 큰 상기 미리 설정된 기준값을 설정하는 기준값 설정부(미도시);를 더 포함할 수 있다.The radiographic imaging apparatus 100 according to the present invention is a reference value setting unit for setting the preset reference value greater than the intrinsic noise value generated even when radiation is irradiated from each pixel as a whole to the plurality of pixels or generated by a radiation source (Not shown); may further include.
기준값 설정부(미도시)는 복수의 픽셀(111) 전체적으로 각각의 픽셀(111)에서 방사선이 비조사되는 경우에도 발생되거나, 방사선원으로 인해 발생되는 고유 노이즈값보다 큰(또는 이상인) 미리 설정된 기준값(R)을 설정할 수 있다. 미리 설정된 기준값(R)은 모든 복수의 픽셀(111)에 공통적으로 발생하는 상기 고유 노이즈값보다 크게(또는 이상으로) 설정될 수 있다. 상기 고유 노이즈는 그 값은 다를 수 있으나 모든 복수의 픽셀(111)에 공통적으로 발생하는 노이즈를 말하며, 장치적 오프셋(offset)에 의해 발생하는 노이즈일 수 있고, 이러한 노이즈는 모든 복수의 픽셀(111)에 공통적으로 발생하기 때문에 방사선 이미지를 이루는 복수의 픽셀(111)의 픽셀값들에서 장치적 오프셋값을 모두 제거하게 되면 큰 문제가 되지 않는다. 예를 들어, 상기 고유 노이즈는 열적(thermal) 원인, 절연성 불량 등의 원인에 의해 방사선에 기인하지 않고 발생하여 복수의 픽셀(111)에 전하를 축적시키는 암전류(dark current) 및 광자가 가진 양자적 성질로 말미암아 생기며, 방사선 발생장치(Radiation Generator)로 인해 발생하는 양자 잡음(quantum noise)을 포함할 수 있다.The reference value setting unit (not shown) is a preset reference value (or greater than or equal to) that is generated even when radiation is not irradiated from each pixel 111 as a whole, or greater than or equal to the intrinsic noise value generated by the radiation source. R) can be set. The preset reference value R may be set to be larger (or higher) than the intrinsic noise value common to all the plurality of pixels 111. The intrinsic noise may have a different value, but refers to noise common to all the plurality of pixels 111, and may be noise generated by a device offset, and the noise may be any of the plurality of pixels 111 ), it is not a big problem if all of the device offset values are removed from the pixel values of the plurality of pixels 111 constituting the radiation image. For example, the intrinsic noise is generated due to a cause of thermal (thermal), poor insulation, and the like due to radiation, and dark current and quantum of photons accumulating charge in the plurality of pixels 111 It is caused by nature, and may include quantum noise generated by a radiation generator.
이처럼, 본 발명에서는 시간에 따라 획득된 복수의 방사선 이미지 중 둘 이상의 방사선 이미지를 선택하여, 선택된 둘 이상의 방사선 이미지 간의 픽셀값 변화량을 고려함으로써, 시간에 따른 입력 신호(input signal)의 변화를 반영하여 각 대상 픽셀의 픽셀값이 랜덤 노이즈를 포함하는 비정상 픽셀값인지를 판단할 수 있다. 이에 따라 복수의 방사선 이미지에서 랜덤 노이즈가 포함된 비정상 픽셀값을 제거하여 보정함으로써, 방사선 이미지의 왜곡을 방지할 수 있고, 방사선 이미지의 품질을 향상시킬 수 있다. 또한, 대상 픽셀의 시간에 따른 픽셀값 변화량과 미리 설정된 기준값을 비교하여 비정상 픽셀값을 판단한 후에 다시 주변 픽셀을 통해 비정상 픽셀값이 맞는지를 검증함으로써, 피사체의 움직임 등에 의한 입력 신호의 변화를 효과적으로 반영하여 각 대상 픽셀의 픽셀값이 랜덤 노이즈를 포함하는 비정상 픽셀값인지를 보다 정확하게 판단할 수 있다. 그리고 랜덤 노이즈의 특성상 한 번 랜덤 노이즈가 발생했던 픽셀에서는 랜덤 노이즈가 발생할 확률이 높기 때문에 비정상 픽셀값이 발생한 대상 픽셀의 위치를 저장하고 이후의 방사선 이미지에서 해당 픽셀의 픽셀값을 비정상 픽셀값으로 간주함으로써, 비정상 픽셀값인지를 판단해야 하는 픽셀의 개수가 줄어들 수 있으며, 복수의 방사선 이미지 전체에 대해 비정상 픽셀값이 있는지를 판단하고 보정하는 영상 처리의 시간이 단축될 수 있다.As described above, in the present invention, by selecting two or more radiographic images among a plurality of radiographic images acquired over time, and taking into account the amount of change in pixel values between the selected two or more radiographic images, the change in the input signal over time is reflected. It may be determined whether the pixel value of each target pixel is an abnormal pixel value including random noise. Accordingly, by correcting by removing abnormal pixel values including random noise from a plurality of radiographic images, distortion of the radiographic image can be prevented, and quality of the radiographic image can be improved. In addition, by comparing the amount of change of the pixel value with time of the target pixel and a preset reference value, and determining the abnormal pixel value again, verifying that the abnormal pixel value is correct through the surrounding pixels, effectively reflecting the change in the input signal due to the movement of the subject. Thus, it can be more accurately determined whether the pixel value of each target pixel is an abnormal pixel value including random noise. Also, because of the nature of random noise, the probability of random noise is high in a pixel where random noise has occurred once, so the location of the target pixel where the abnormal pixel value is generated is stored and the pixel value of the pixel is regarded as an abnormal pixel value in the subsequent radiation image. By doing so, the number of pixels that need to be determined to be abnormal pixel values can be reduced, and the time for image processing to determine and correct whether there are abnormal pixel values for the entire plurality of radiographic images can be shortened.
이상에서 본 발명의 바람직한 실시예에 대하여 도시하고 설명하였으나, 본 발명은 상기한 실시예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 본 발명이 속하는 분야에서 통상의 지식을 가진 자라면 이로부터 다양한 변형 및 균등한 타 실시예가 가능하다는 점을 이해할 것이다. 따라서, 본 발명의 기술적 보호범위는 아래의 특허청구범위에 의해서 정하여져야 할 것이다.Although the preferred embodiments of the present invention have been shown and described above, the present invention is not limited to the above-described embodiments, and general knowledge in the field to which the present invention pertains without departing from the gist of the present invention claimed in the claims Anyone who has a will understand that various modifications and other equivalent embodiments are possible. Therefore, the technical protection scope of the present invention should be defined by the following claims.
Claims (14)
- (a) 복수의 픽셀을 구비하는 방사선 검출패널을 이용하여 시간에 따라 복수의 방사선 이미지를 획득하는 과정;(a) obtaining a plurality of radiation images over time using a radiation detection panel having a plurality of pixels;(b) 상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지에서 상기 복수의 픽셀 중 대상 픽셀을 선택하는 과정;(b) selecting a target pixel among the plurality of pixels from two or more selected radiation images among the plurality of radiation images;(c) 선택된 대상 픽셀의 시간에 따른 픽셀값 변화량을 산출하는 과정;(c) calculating a pixel value change amount over time of the selected target pixel;(d) 산출된 상기 대상 픽셀의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값과 비교하는 과정; 및(d) comparing the calculated amount of change of the pixel value with time of the target pixel with a preset reference value; And(e) 상기 대상 픽셀의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지 중 적어도 하나의 방사선 이미지에서 상기 대상 픽셀의 픽셀값을 비정상값으로 판단하는 과정;을 포함하는 방사선 영상 처리방법.(e) determining a pixel value of the target pixel as an abnormal value in at least one of the selected two or more radiographic images when the amount of change in the pixel value over time of the target pixel is greater than or equal to the preset reference value; Radiographic image processing method.
- 청구항 1에 있어서,The method according to claim 1,비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 검증하는 과정;을 더 포함하고,Further comprising the step of verifying the pixel value of the target pixel determined to be an abnormal value;상기 대상 픽셀의 픽셀값을 검증하는 과정은,The process of verifying the pixel value of the target pixel,상기 선택된 둘 이상의 방사선 이미지로부터 산출된 상기 대상 픽셀에 인접한 주변 픽셀의 시간에 따른 픽셀값 변화량과 상기 미리 설정된 기준값을 비교하는 과정; 및Comparing the pixel value change amount with time of a neighboring pixel adjacent to the target pixel calculated from the selected two or more radiographic images and the preset reference value; And상기 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 확인하는 과정을 포함하는 방사선 영상 처리방법.And checking the pixel value of the target pixel determined as the abnormal value as an abnormal value when the number of surrounding pixels over the preset reference value is less than or equal to the preset determination number. Way.
- 청구항 1에 있어서,The method according to claim 1,비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 검증하는 과정;을 더 포함하고,Further comprising the step of verifying the pixel value of the target pixel determined to be an abnormal value;상기 대상 픽셀의 픽셀값을 검증하는 과정은,The process of verifying the pixel value of the target pixel,상기 선택된 둘 이상의 방사선 이미지로부터 산출된 상기 대상 픽셀에 인접한 주변 픽셀의 시간에 따른 픽셀값 변화량과 상기 대상 픽셀의 시간에 따른 픽셀값 변화량을 비교하는 과정; 및Comparing a pixel value change amount over time of a neighboring pixel adjacent to the target pixel calculated from the selected two or more radiographic images and a pixel value change amount over time of the target pixels; And상기 대상 픽셀의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 확인하는 과정을 포함하는 방사선 영상 처리방법.The pixel value of the target pixel determined as the abnormal value is an abnormal value when the number of pixel values variation over time and the number of surrounding pixels having a pixel value variation within a preset allowable deviation range is equal to or less than a preset determination number. A radiographic image processing method comprising the step of confirming.
- 청구항 1에 있어서,The method according to claim 1,상기 복수의 방사선 이미지 전체에 대해 상기 (b) 과정 내지 상기 (e) 과정을 반복하며,The above steps (b) to (e) are repeated for the whole of the plurality of radiographic images,비정상값으로 판단된 픽셀값이 발생한 대상 픽셀의 위치를 저장하는 과정;을 더 포함하고,Further comprising; storing the position of the target pixel where the pixel value determined as the abnormal value occurs;상기 대상 픽셀의 위치를 저장하는 과정 이후의 (b) 과정에서 위치가 저장된 픽셀을 상기 대상 픽셀의 선택에서 제외하고 상기 위치가 저장된 픽셀의 픽셀값을 비정상값으로 판단하는 방사선 영상 처리방법.A radiographic image processing method of determining a pixel value of a pixel in which the location is stored as an abnormal value by excluding a pixel in which the location is stored in the selection of the target pixel in the step (b) after the process of storing the location of the target pixel.
- 청구항 1에 있어서,The method according to claim 1,상기 미리 설정된 기준값은 상기 복수의 픽셀 전체적으로 각각의 픽셀에서 방사선이 비조사되는 경우에도 발생되거나, 방사선원으로 인해 발생되는 고유 노이즈값보다 큰 방사선 영상 처리방법.The preset reference value is generated even when radiation is not irradiated from each pixel as a whole to the plurality of pixels, or a method of processing a radiographic image larger than an inherent noise value generated by a radiation source.
- 청구항 1에 있어서,The method according to claim 1,상기 (b) 과정은 상기 복수의 방사선 이미지 중 시간적으로 연속되는 둘 이상의 방사선 이미지를 선택하여 수행되는 방사선 영상 처리방법.The step (b) is a radiographic image processing method performed by selecting two or more temporally continuous radiographic images among the plurality of radiographic images.
- 청구항 6에 있어서,The method according to claim 6,상기 (b) 과정에서는 상기 복수의 방사선 이미지 중 시간적으로 연속되는 2개의 방사선 이미지를 선택하고,In the step (b), two consecutive radiographic images are selected from the plurality of radiographic images,상기 (e) 과정에서는 상기 대상 픽셀의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 2개의 방사선 이미지 중 상기 대상 픽셀의 픽셀값이 높은 방사선 이미지에서 상기 대상 픽셀의 픽셀값을 비정상값으로 판단하는 방사선 영상 처리방법.In the step (e), a pixel value of the target pixel is determined as an abnormal value in a radiographic image having a high pixel value of the target pixel among two radiographic images in which the amount of change in the pixel value over time of the target pixel is greater than or equal to the preset reference value. Radiographic image processing method.
- 청구항 1에 있어서,The method according to claim 1,상기 복수의 픽셀 중 미리 설정된 범위 내의 세기로 방사선이 입사되는 픽셀 영역을 검출하는 과정;을 더 포함하고,Further comprising the step of detecting a pixel region to which the radiation is incident at an intensity within a predetermined range among the plurality of pixels;상기 (b) 과정에서는 상기 픽셀 영역 중 상기 대상 픽셀을 선택하는 방사선 영상 처리방법.In the step (b), a radiographic image processing method of selecting the target pixel among the pixel areas.
- 청구항 1에 있어서,The method according to claim 1,비정상값으로 판단된 대상 픽셀의 픽셀값을 동일한 방사선 이미지 상의 인접한 주변 픽셀의 픽셀값을 내삽하여 보정하는 과정;을 더 포함하는 방사선 영상 처리방법.And correcting the pixel value of the target pixel determined as the abnormal value by interpolating the pixel value of adjacent surrounding pixels on the same radiation image.
- 복수의 픽셀을 구비하며, 조사되는 방사선에 의해 시간에 따라 복수의 방사선 이미지를 획득하기 위한 전기신호를 생성하는 방사선 검출패널;A radiation detection panel having a plurality of pixels and generating electrical signals for acquiring a plurality of radiation images over time by the irradiated radiation;상기 복수의 방사선 이미지 중 선택된 둘 이상의 방사선 이미지에서 상기 복수의 픽셀 중 선택된 대상 픽셀의 시간에 따른 픽셀값 변화량 및 상기 대상 픽셀에 인접한 주변 픽셀의 시간에 따른 픽셀값 변화량을 산출하는 변화량 산출부;A change amount calculating unit for calculating a pixel value change amount over time of a target pixel selected from the plurality of pixels from two or more radiation images selected from the plurality of radiation images and a pixel value change amount over time of neighboring pixels adjacent to the target pixel;산출된 상기 대상 픽셀의 시간에 따른 픽셀값 변화량을 미리 설정된 기준값과 비교하는 변화량 비교부; 및A change amount comparison unit comparing the calculated change amount of the pixel value with time of the target pixel with a preset reference value; And상기 대상 픽셀의 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 경우에 상기 선택된 둘 이상의 방사선 이미지 중 적어도 하나의 방사선 이미지에서 상기 대상 픽셀의 픽셀값을 비정상값으로 판단하는 픽셀값 이상판단부;를 포함하는 방사선 촬영장치.A pixel value abnormality determination unit determining a pixel value of the target pixel as an abnormal value in at least one radiation image among the selected two or more radiographic images when the amount of change of the pixel value over time of the target pixel is equal to or greater than the preset reference value; Radiography apparatus comprising.
- 청구항 10에 있어서,The method according to claim 10,상기 시간에 따른 픽셀값 변화량이 상기 미리 설정된 기준값 이상인 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 결정하는 제1 픽셀값 검증부;를 더 포함하는 방사선 촬영장치.A first pixel value verification unit configured to determine a pixel value of the target pixel determined as the abnormal value as an abnormal value when the number of surrounding pixels having a change amount of the pixel value with respect to the time is equal to or less than a preset determination number; Radiographic apparatus further comprising.
- 청구항 10에 있어서,The method according to claim 10,상기 대상 픽셀의 시간에 따른 픽셀값 변화량과 미리 설정된 허용편차 범위 이내인 픽셀값 변화량을 갖는 주변 픽셀의 개수가 미리 설정된 판단 개수 이하인 경우에 상기 비정상값으로 판단된 상기 대상 픽셀의 픽셀값을 비정상값으로 결정하는 제2 픽셀값 검증부;를 더 포함하는 방사선 촬영장치.The pixel value of the target pixel determined as the abnormal value is an abnormal value when the number of pixel values change over time and the number of surrounding pixels having a pixel value variation within a preset allowable deviation range is equal to or less than a preset determination number. And a second pixel value verifying unit determined as.
- 청구항 10에 있어서,The method according to claim 10,비정상값으로 판단된 픽셀값이 발생한 대상 픽셀의 위치를 저장하는 위치저장부;를 더 포함하는 방사선 촬영장치.And a location storage unit which stores a location of a target pixel in which a pixel value determined as an abnormal value occurs.
- 청구항 10에 있어서,The method according to claim 10,상기 복수의 픽셀 전체적으로 각각의 픽셀에서 방사선이 비조사되는 경우에도 발생되거나, 방사선원으로 인해 발생되는 고유 노이즈값보다 큰 상기 미리 설정된 기준값을 설정하는 기준값 설정부;를 더 포함하는 방사선 촬영장치.And a reference value setting unit configured to set the preset reference value greater than the intrinsic noise value generated even when radiation is irradiated from each pixel as a whole, or greater than the inherent noise value generated by the radiation source.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020190011479A KR102206996B1 (en) | 2019-01-29 | 2019-01-29 | Method for processing radiation image and radiographic apparatus |
KR10-2019-0011479 | 2019-01-29 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2020159201A1 true WO2020159201A1 (en) | 2020-08-06 |
Family
ID=71842301
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/KR2020/001303 WO2020159201A1 (en) | 2019-01-29 | 2020-01-28 | Radiographic image processing method and radiography apparatus |
Country Status (2)
Country | Link |
---|---|
KR (1) | KR102206996B1 (en) |
WO (1) | WO2020159201A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113112495B (en) * | 2021-04-30 | 2024-02-23 | 浙江华感科技有限公司 | Abnormal image processing method and device, thermal imaging equipment and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050015720A (en) * | 2003-08-07 | 2005-02-21 | 삼성전자주식회사 | Error frame detecting appatrus of moving picture data and method thereof |
US20120328071A1 (en) * | 2009-08-05 | 2012-12-27 | Telesystems Co., Ltd. | Radiation imaging apparatus and imaging method using radiation |
KR20130096293A (en) * | 2010-11-12 | 2013-08-29 | 인디안 인스티튜트 오브 테크놀로지, 카라그푸르 | Method and apparatus for detecting the bad pixels in sensor array and concealing the error |
KR101367747B1 (en) * | 2009-05-12 | 2014-02-26 | 캐논 가부시끼가이샤 | Radiographic apparatus and control method for the same |
US20180341029A1 (en) * | 2015-10-30 | 2018-11-29 | Canon Kabushiki Kaisha | Radiation imaging system, signal processing apparatus, and, radiographic image signal processing method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1489435B1 (en) | 2003-06-18 | 2006-02-08 | Canon Kabushiki Kaisha | Radiographic apparatus and radiographic method |
-
2019
- 2019-01-29 KR KR1020190011479A patent/KR102206996B1/en active IP Right Grant
-
2020
- 2020-01-28 WO PCT/KR2020/001303 patent/WO2020159201A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050015720A (en) * | 2003-08-07 | 2005-02-21 | 삼성전자주식회사 | Error frame detecting appatrus of moving picture data and method thereof |
KR101367747B1 (en) * | 2009-05-12 | 2014-02-26 | 캐논 가부시끼가이샤 | Radiographic apparatus and control method for the same |
US20120328071A1 (en) * | 2009-08-05 | 2012-12-27 | Telesystems Co., Ltd. | Radiation imaging apparatus and imaging method using radiation |
KR20130096293A (en) * | 2010-11-12 | 2013-08-29 | 인디안 인스티튜트 오브 테크놀로지, 카라그푸르 | Method and apparatus for detecting the bad pixels in sensor array and concealing the error |
US20180341029A1 (en) * | 2015-10-30 | 2018-11-29 | Canon Kabushiki Kaisha | Radiation imaging system, signal processing apparatus, and, radiographic image signal processing method |
Also Published As
Publication number | Publication date |
---|---|
KR102206996B1 (en) | 2021-01-25 |
KR20200094019A (en) | 2020-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5828725A (en) | Processing images for removal of artifacts | |
US7623162B2 (en) | Signal processing device and method for reducing influence on column noise detection from defective pixels | |
EP1487193B1 (en) | Method and apparatus for correcting defect pixels in radiation imaging, computer program and computer-readable recording medium | |
US7567277B2 (en) | Image-pickup device and signal processing method for noise correction | |
KR101001431B1 (en) | Pixel defect detecting/correcting device and pixel defect detecting/correcting method | |
JP4995193B2 (en) | X-ray diagnostic imaging equipment | |
JP5269841B2 (en) | Image processing device | |
US11284852B2 (en) | Radiation imaging apparatus, radiation imaging system, method of controlling radiation imaging apparatus, and non-transitory computer-readable storage medium | |
JP2002281384A (en) | Method for activating image system of medical inspection device and medical inspection device | |
WO2017034075A1 (en) | Radiation image processing method and radiographic system | |
WO2020159201A1 (en) | Radiographic image processing method and radiography apparatus | |
US7471842B2 (en) | Method and apparatus for dynamically detecting pixel values | |
EP1178676B1 (en) | Method and device for detecting defective pixels of a solid-state colour image-pickup apparatus | |
JP7308609B2 (en) | RADIATION IMAGING APPARATUS, RADIATION IMAGING SYSTEM, CONTROL METHOD OF RADIATION IMAGING APPARATUS, AND PROGRAM | |
JP5262953B2 (en) | Image processing apparatus, image processing method, and program | |
JP2002034961A (en) | Radiographing apparatus and radiographing method | |
WO2015115714A1 (en) | Apparatus and method for processing scan data in the event of esd input | |
JP3257131B2 (en) | Automatic defect detection device for solid-state imaging device and solid-state imaging device | |
US7986354B2 (en) | Method for correcting pixel defect of image pickup device | |
US20230179878A1 (en) | Image sensor and camera | |
JP2017164374A (en) | Radiation detecting device and radiation inspecting apparatus | |
JP2003158744A (en) | Pixel defect detection/correction device and image pick- up device using the same | |
US20230309947A1 (en) | Radiation imaging system, control apparatus, radiation imaging apparatus, control methods therefor, and non-transitory computer-readable storage medium | |
CN113835115B (en) | CT detector state detection method | |
TWI528817B (en) | Image capture devices and image processing method thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20748586 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20748586 Country of ref document: EP Kind code of ref document: A1 |