WO2015052890A1 - X線画像処理装置および方法、並びにx線撮像装置、プログラム - Google Patents
X線画像処理装置および方法、並びにx線撮像装置、プログラム Download PDFInfo
- Publication number
- WO2015052890A1 WO2015052890A1 PCT/JP2014/004962 JP2014004962W WO2015052890A1 WO 2015052890 A1 WO2015052890 A1 WO 2015052890A1 JP 2014004962 W JP2014004962 W JP 2014004962W WO 2015052890 A1 WO2015052890 A1 WO 2015052890A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- analysis value
- value
- analysis
- image processing
- statistic
- Prior art date
Links
- 238000012545 processing Methods 0.000 title claims abstract description 90
- 238000003384 imaging method Methods 0.000 title claims description 5
- 238000000034 method Methods 0.000 title abstract description 36
- 238000004458 analytical method Methods 0.000 claims abstract description 298
- 238000012937 correction Methods 0.000 claims description 45
- 238000003672 processing method Methods 0.000 claims description 9
- 238000012790 confirmation Methods 0.000 claims description 4
- 239000006185 dispersion Substances 0.000 description 13
- 230000008569 process Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- RTZKZFJDLAIYFH-UHFFFAOYSA-N Diethyl ether Chemical compound CCOCC RTZKZFJDLAIYFH-UHFFFAOYSA-N 0.000 description 6
- 238000010191 image analysis Methods 0.000 description 6
- 238000010187 selection method Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 4
- 238000002594 fluoroscopy Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 238000004321 preservation Methods 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 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/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
- A61B6/487—Diagnostic techniques involving generating temporal series of image data involving fluoroscopy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Definitions
- the present invention relates to an X-ray image processing apparatus and method for performing gradation processing on an X-ray image, an X-ray imaging apparatus, and a program.
- the X-ray fluoroscopy system monitors the X-ray dose reaching the sensor and controls the incident dose so that the value becomes constant. To do. If the arrival dose to the sensor is small, a command to increase the incident dose is issued, and if the arrival dose to the sensor is large, a command to decrease the incidence dose is issued to control the X-ray irradiation amount. In this specification, such dose control is referred to as automatic dose control (hereinafter referred to as ADC).
- ADC automatic dose control
- the X-ray irradiation amount is optimized by the ADC, and then the subject is continuously seen at a fixed dose.
- ADC is stable.
- the ADC is stable, the gradation of the X-ray image observed by the user does not change, and an image without luminance fluctuation is displayed.
- the output generally has a certain degree of fluctuation.
- the brightness fluctuation may occur in the X-ray image observed by the user.
- an X-ray image is optimized by further performing image processing on an image obtained by optimization by an ADC, and brightness fluctuation is suppressed.
- the input image is analyzed, a gradation conversion curve is created so that the region of interest always has a constant luminance value, and the result of processing the image using the created gradation conversion curve is the final image. Is output.
- a subject region is extracted by extracting a collimation region or an unexposed region (region where X-rays are directly incident) in an X-ray image, and the extracted subject region has an optimum gradation.
- the analysis value for becoming is calculated. Further, by improving the accuracy of such analysis, a more optimal final image can be output.
- Patent Document 1 As a method for correcting the influence of such flicker, there is a method of reducing the influence of the failed analysis value by smoothing the analysis value calculated by the image analysis in the time axis direction (Patent Document 1). As another method, there is a method of calculating an analysis value by analyzing an image, and adding a difference value between the analysis value and an analysis value smoothed in the time axis direction to the analysis value (Patent Document 2). A value for stabilization from the smoothing in the time direction is calculated, and flicker is removed by matching with the value.
- the influence of flicker can be reduced by smoothing, but if the amount of failure is large, the reduction effect is weakened and flicker occurs.
- the analysis value fluctuates due to the fluctuation of the X-ray generator.
- an average value is used as a correction value for a value determined to have failed image analysis, the fluctuation of the X-ray generator cannot be taken into account, and the luminance fluctuates.
- An X-ray image processing apparatus has the following configuration. That is, An X-ray image processing apparatus that performs image processing on frames continuously obtained by an X-ray sensor, Analyzing means for analyzing each of the continuously obtained frames to obtain a first analysis value and a second analysis value; A calculating means for calculating a time-direction statistic of the first analysis value and the second analysis value of the continuously obtained frames as a first statistic and a second statistic; Determination means for determining whether to correct the first analysis value based on the first statistic and the second statistic; And a correction unit that corrects the first analysis value based on the first statistic when it is determined to be corrected by the determination unit.
- the block diagram which shows the function structural example of the X-ray image processing apparatus by embodiment.
- the block diagram which shows the hardware structural example of an X-ray image processing apparatus.
- the flowchart which shows the acquisition process of the analysis value for a gradation process.
- the block diagram which shows the function structural example of the X-ray image processing apparatus by 4th embodiment.
- the flowchart explaining the selection procedure of a solution The flowchart explaining the selection procedure of a solution.
- the X-ray image processing method according to the present embodiment includes an analysis unit 101, a dose state determination unit 102, a statistic calculation unit 103, an analysis value comparison unit 104, an analysis value correction unit 105, and a gradation processing unit 106. Realized by the device.
- the analysis unit 101 receives an image obtained by performing predetermined preprocessing on the X-ray image acquired by the X-ray sensor (hereinafter referred to as preprocessed image), and outputs an analysis value for use in gradation processing. .
- Pre-processing is processing for correcting sensor characteristics such as offset correction, gain correction, and defect correction, and is processing for maintaining a correlation with surrounding pixels.
- the dose state determination unit 102 receives the analysis value output from the analysis unit 101 and outputs a determination result as to whether or not the ADC is stable.
- the statistic calculation unit 103 receives the analysis value from the analysis unit 101 and the determination result from the dose state determination unit 102 and calculates a statistic in the time direction of the analysis value.
- the analysis value comparison unit 104 receives the statistics in the time direction calculated by the statistic calculation unit 103 as input, and determines whether or not the analysis value is an appropriate value. When the analysis value comparison unit 104 determines that the analysis value is not an appropriate value, the analysis value correction unit 105 outputs a result of correcting the analysis value to an appropriate value.
- the gradation processing unit 106 performs gradation processing on the preprocessed image using either the analysis value output from the analysis unit 101 or the corrected analysis value output from the analysis value correction unit 105. Note that the gradation processing unit 106 determines whether to use the corrected analysis value based on the determination result of the analysis value comparison unit 104.
- FIG. 1 When the configuration of FIG. 1 is realized using an information processing apparatus such as a computer (hereinafter referred to as a control PC), an example of the hardware configuration is as shown in FIG.
- the control PC 201 is connected to the X-ray sensor 202 as an X-ray detection panel and the X-ray generator 210 by a gigabit ether 222.
- the signal line is not limited to gigabit ether, and may be a CAN (Controller Area Network) or an optical fiber.
- the X-ray sensor 202 outputs an electrical signal corresponding to the X-ray dose, and the output electrical signal is stored in the external storage unit 211 as an X-ray image, for example.
- An external display unit 209, an external storage unit 211, and a network interface unit 212 are connected to the gigabit ether 222.
- a CPU (Central Processing Unit) 203, a RAM (Random Access Memory) 204, a ROM (Read Only Memory) 205, an input unit 206, a display unit 207, and a storage unit 208 are connected to the control PC 201, for example.
- the A command is sent to the X-ray sensor 202 and the external display unit 209 via the control PC 201.
- the processing content for each shooting mode is stored as a software module in the storage unit 208, and the software module is read into the RAM 204 and executed by the CPU 203.
- each configuration shown in FIG. 1 may be mounted as a dedicated image processing board, and an optimal mounting may be performed according to the purpose. Details of the X-ray image processing method according to the first embodiment will be described below.
- an X-ray fluoroscopy system is exemplified as the X-ray imaging apparatus.
- the X-ray image processing in the X-ray fluoroscopy system according to the first embodiment will be described with reference to the configuration diagram of FIG. 1 and the flowchart of FIG.
- frames (pre-processed X-ray images) obtained continuously from the X-ray sensor 202 are subjected to image processing (gradation processing in the present embodiment).
- the analysis unit 101 analyzes each of the continuously obtained frames and outputs a plurality of analysis values including an analysis value 1 and an analysis value 2 described later.
- the preprocessed image is stored in, for example, the external storage unit 211 and is read by the analysis unit 101 or the gradation processing unit 106.
- the analysis unit 101 calculates the analysis value 1 and the analysis value 2 by performing analysis processing on different regions of one frame. First, the analysis unit 101 outputs a value representing the pixel value of the subject area from the preprocessed image as the analysis value 1 (step 301). There are various methods for extracting the subject area, and any known method may be used. For example, a subject area can be extracted by extracting a collimation area and a blank area from an image.
- Patent Document 3 There is, for example, a method described in Patent Document 3 for extracting a collimation region.
- the X axis and the Y axis are set along two adjacent sides of the outline of the rectangular irradiation field, and the image data is added and summed in the set X axis direction and Y axis direction.
- the addition total data in the irradiation field is higher than that in the irradiation field where X-rays are hardly irradiated.
- a position on the Y axis where the sum total data in the X axis direction is equal to or greater than the predetermined threshold TH and a position on the X axis where the sum total data in the Y axis direction is equal to or greater than the predetermined threshold TH are calculated.
- a rectangular area surrounded by the straight line in the X-axis direction at the calculated position on the Y-axis and the straight line in the Y-axis direction at the position on the X-axis is defined as a collimation area.
- the blank region is a region where the irradiated X-rays reach the sensor directly without passing through the subject.
- the background missing area is concentrated in the high pixel value area. For this reason, by removing several percent from the high pixel value side of the histogram range, it is possible to extract (remove) the missing region.
- region is not limited to these.
- the remaining area obtained by removing the collimation area and the blank area obtained as described above from the X-ray image is the subject area.
- a value representing the subject area for example, an average value or an intermediate value of the subject area can be used.
- the calculated analysis value 1 is stored in the external storage unit 211, for example.
- the image to be processed by the analysis unit 101 is a preprocessed image. However, in addition to the preprocessing, an image after processing such as log conversion, enhancement processing, and noise suppression processing may be input.
- the analysis unit 101 calculates an analysis value 2 used for the ADC from the preprocessed image (step 302).
- the analysis value used for the ADC is calculated by a simpler method than the analysis value 1 in order to suppress the risk of analysis failure. For example, a representative value of a region of a predetermined size at a predetermined position of the preprocessed image is used. In the present embodiment, the average value of an area of several inches square in the center of the preprocessed image is calculated as the analysis value 2.
- the calculated analysis value 1 and analysis value 2 are stored in the external storage unit 211 in association with the frame number, for example.
- the dose state determination unit 102 calculates a statistic in the time direction for the analysis value 2 of frames obtained continuously, and determines whether or not the X-ray dose is stable. More specifically, the dose state determination unit 102 inputs the analysis value 2 of the current frame (hereinafter, the current frame) and the analysis value 2 calculated in the past frame from the external storage unit 211, and these analysis values Statistical processing. Whether the ADC is stable or unstable is determined based on the statistics obtained by the statistical processing (steps 303 and 304). Here, as the statistic, the variance or the standard deviation of the analysis value 2 can be used. When the ADC is activated, the analysis value 2 fluctuates. Therefore, the fluctuation amount of the analysis value 2 is calculated and threshold processing is performed.
- the value of the analysis value 2 is R 2 (t), t is the frame number, and the variance V 2 (t) is calculated as the following expression 1 as an index representing the fluctuation amount. Note that n 2 represents the number of samplings in the time direction.
- the gradation processing unit 106 performs gradation processing of the preprocessed image using the uncorrected analysis value 1 for gradation processing (step 305).
- the statistic calculation unit 103 and the analysis value comparison unit 104 determine whether or not to execute the correction process of the analysis value 1. For this determination, the statistic calculation unit 103 calculates a statistic in the time direction of the analysis value 1 of frames obtained continuously. That is, the statistic calculation unit 103 reads the analysis value 1 in the current frame and the result of the analysis value 1 calculated in the past frame from the external storage unit 211, and whether the analysis value 1 in the current frame has succeeded or failed. A statistic for determining whether or not it is determined is calculated (step 306). For this statistic, the variance and standard deviation can be used as in the case of the analysis value 2.
- the variance V 1 (t) in the time direction of the analysis value 1 uses Equation 2 as the statistic. Is calculated.
- the value of the analysis value 1 is R 1 (t)
- t is the frame number
- n 1 is the number of samplings in the time direction.
- the analysis value comparison section 104 as input and variance V 2 of variance V 1 (t) and the analysis value 2 of the analysis value 1 (t), determines whether the analysis value 1 has failed (step 307, 308).
- the analysis value comparison unit 104 calculates a dispersion ratio obtained by dividing V 1 (t) by V 2 (t) for this determination, and compares the calculated ratio with a threshold value. If the ADC is stable and the analysis value 1 does not fail, the variances V 1 (t) and V 2 (t) should be approximately the same value. Therefore, when the variance ratio is less than or equal to the threshold value, the analysis value 1 is determined to be correct. If it is determined that the analysis value 1 is correct, the gradation processing unit 106 performs gradation processing on the preprocessed image using the analysis value 1 without correction (step 305).
- step 308 when it is determined that analysis value 1 is unsuccessful (acquisition of analysis value 1 is unsuccessful), correction of analysis value 1 in step 310 may be executed.
- the analysis value 1 since it is not known which analysis value among the sampled analysis values 1 is a failure analysis value, if the analysis value 1 is simply corrected according to the result of step 308, the analysis value that does not need to be corrected 1 may be corrected. Therefore, in this embodiment, it is confirmed whether or not acquisition of the analysis value 1 of the current frame has failed by referring to a predetermined number of frames in the past (step 309). By this confirmation processing, it is possible to more accurately determine whether or not correction for the analysis value 1 of the current frame is necessary. If the failure of analysis value 1 of the current frame cannot be confirmed in step 309, correction of analysis value 1 is prohibited and the analysis value 1 is used as it is.
- the confirmation process of this embodiment will be described.
- the analysis value used for calculating the dispersion V 1 (t) includes the analysis value that has failed.
- the variance ratio up to the previous frame is equal to or less than the threshold value (step 309). If the variance ratio is less than or equal to the threshold until the previous frame, but the newly added analysis value 1 in the current frame exceeds the variance ratio threshold, the failed analysis is an analysis in the current frame. Become.
- the gradation processing unit 106 performs gradation processing on the preprocessed image using the analysis value 1 without correction (step 305).
- analyzing value correcting unit 105 analyzes value 1 and the variance V 1 (t), and inputs a dispersion V 2 of the analysis value 2 (t), the analysis value 1
- the analysis value 3 corrected by is output (step 310).
- the analysis value correction unit 105 corrects the analysis value 1 using Equation 3 so that the variance ratio including the analysis value 1 of the current frame is within the threshold thF. If the analysis value 3 is R 3 (t), the quadratic equation for R 3 (t) is solved.
- Equation 3 when the solution is a real number, two results R 31 (t) and R 32 (t) are obtained.
- a solution having the smallest difference from the analysis value 1 before correction is selected as the analysis value after correction.
- This process is shown in steps 6011, 6012, and 6013 in FIG. 6A. That is, the absolute value of the difference between each of R 31 (t) and R 32 (t) and R 1 (t) is calculated (step 6011), and is calculated from R 31 (t) and R 32 (t). The smaller difference is selected (Equation 4, Step 6012, S6013).
- the analysis value 1 is corrected so that the variance ratio between the analysis value 1 and the analysis value 2 obtained from two different regions is maintained. Therefore, even if the fluctuation of the analysis value includes the fluctuation of the X-ray generator, the fluctuation component can be corrected so as to be retained. As a result, flicker can be effectively suppressed.
- the correction of the analysis value 1 is necessary, so that there is some fluctuation in the X-ray dose. Even if it occurs, it becomes possible to effectively suppress flicker.
- the gradation processing unit 106 performs gradation processing on the preprocessed image using the analysis value 3 (step 311).
- the image to be processed by the gradation processing unit 106 is a preprocessed image, it may be an image after processing such as enhancement processing and noise suppression processing in addition to preprocessing.
- V 1 (t) the ratio of V 1 (t) to V 2 (t) is used.
- V 2 (t) is equal to or lower than the threshold value
- V The analysis value 1 may be corrected so that 1 (t) becomes a predetermined value. That is, when the dose state determination unit 102 determines that the X-ray dose is stable based on the statistic of the analysis value 2, it is determined whether the correction of the analysis value 1 is necessary or the correction value of the analysis value 1. May be executed based on the statistic of analysis value 1.
- the X-ray image processing apparatus performs image processing on frames continuously obtained by the X-ray sensor, and analyzes each of the continuously obtained frames to obtain the first analysis value and the second analysis value.
- Analysis means for example, the analysis unit 101 for obtaining the analysis value of the first analysis value and the second analysis value in the time direction of the first analysis value and the second analysis value of the frames obtained in succession.
- a calculation unit for example, a statistic calculation unit 103 that calculates as a statistic, and a determination unit that determines whether to correct the first analysis value based on the first statistic and the second statistic (for example, a dose state determination unit 102), a correction unit (for example, an analysis value correction unit 105) that corrects the first analysis value based on the first statistic when it is determined to be corrected by the determination unit, Is provided. Therefore, flicker can be effectively suppressed even if the X-ray irradiation dose fluctuates.
- the X-ray image processing apparatus performs image processing on frames continuously obtained by the X-ray sensor, and analyzes means for analyzing each of the continuously obtained frames to obtain an analysis value (for example, the analysis unit 101). And a calculation means (for example, a statistic calculation unit 103) for calculating a statistic in the time direction of analysis values of frames obtained continuously, and a correction means (for example, an analysis value) for correcting the analysis value based on the statistic. This can be realized even with a configuration including the correction unit 105).
- ⁇ Second embodiment> In the first embodiment, the processing when the analysis value 1 fails in the current frame has been described. In the second embodiment, processing in the next frame when it is determined that analysis value 1 has failed in the previous frame will be described with reference to the configuration diagram in FIG. 1 and the flowchart in FIG. 3.
- the external storage unit 211 is used as a holding unit for holding a predetermined number of past frames. Steps 301 to 306 are the same as the processing up to the previous frame, and are the same as those in the first embodiment, and thus the description thereof is omitted.
- Analysis value comparison section 104 as input and variance V 2 analysis value 2 and variance V 1 of the analysis value 1 (t) (t), determines whether the analysis value 1 has failed (step 307, 308).
- the analysis value comparison unit 104 calculates a dispersion ratio obtained by dividing V 2 (t) by V 1 (t) and compares it with a threshold value. When the dispersion ratio exceeds the threshold value in a state where the ADC is stable, the analysis value for which the variance V 1 (t) is calculated includes the failed analysis value. In order to clarify that the analysis in the target frame has failed, it is confirmed that the variance value up to the previous frame held in the holding unit is equal to or less than the threshold value (step 309).
- the variance value when calculating the variance value up to the previous frame in step 309, the variance value is calculated without using the analysis value 1 of the frame in which the analysis value 1 is determined to have failed. That is, in step 309, the variance value is calculated using the analysis value 1 that has been successfully analyzed.
- FIGS. 4A to 4C A method of selecting the analysis value 1 calculated in the past frame is shown in FIGS. 4A to 4C.
- 4A to 4C are diagrams showing the storage status in the holding unit of the external storage unit 211 when the number of samples is four.
- the holding unit stores the value R 1 (t) of the analysis value 1 at the frame number t, and holds the analysis value 1 of the past four frames.
- the oldest value of R 1 (1) is discarded and a new value of R 1 (5) is stored. This is because the oldest one is discarded in order to compare with the latest distributed state.
- R 1 (6) is determined from these four values.
- the value of R 1 (5) is not stored in the holding unit. This is because using the failed analysis value 1 makes it impossible to clarify that the analysis in the target frame has failed.
- the analysis value 1 continues to fail, the old variance value is maintained in the holding unit indefinitely. Therefore, when the variance value is gradually changed due to a trend component or the like, correct determination cannot be performed.
- the analysis of the current frame is forcibly performed.
- the content held in the holding unit is updated with value 1.
- FIG. 4C is a diagram for explaining a storage state in the holding unit when the threshold value (the predetermined number of frames described above) is 9 frames.
- the threshold value the predetermined number of frames described above
- R 1 (10) was an analysis value that failed, since the oldest R 1 (1) exceeded the threshold value (9 frames), it is forcibly discarded and the value of R 1 (10) is stored. At this time, since the dispersion value up to the previous frame becomes equal to or higher than the threshold, the correction of the analysis value R 1 is not performed. In this way, when the analysis value 1 of the past frame has failed, the determination is performed without using that value.
- the old analysis value 1 is discarded. However, if there is room in the memory, it may be retained without being discarded.
- a new selection method may be used.
- the processing from step 310 to step 311 is the same as in the first embodiment.
- the failed analysis value 1 is not excluded, but as in step 309, the failed analysis value 1 may be excluded and the dispersion ratio may be calculated.
- the analysis value 1 that has been determined to have failed in the past is not used. It is not limited to this.
- the statistic is calculated by excluding the analysis value determined to be corrected (analysis value 1 determined as failure) from the analysis value 1 of the past frame. You may make it perform.
- the method shown in FIG. 6A is shown as a method for selecting one solution when the solution of the quadratic equation related to the analysis value 3 is a real number.
- the method is not limited to this.
- a modification of such a solution selection method is shown.
- a solution is selected by the following process. (1) Among the plurality of solutions, a solution that matches the sign of the differential value in the time direction of the analysis value 1 with the sign of the differential value in the time direction of the analysis value 2 is selected as the analysis value 1 after correction.
- step 310 the analysis value correction unit 105, the analysis value of 1 and the variance V 1 (t), and inputs the variance V 2 of the analysis value 2 (t), and failure determination has been analyzed value 1 is corrected
- An analysis value 3 is generated.
- Equation 3 is used as shown in the first embodiment so that the dispersion ratio falls within the threshold thF.
- the solution becomes a real number, two solutions (R 31 (t) and R 32 (t)) are obtained, so it is necessary to select one of the two.
- a solution is selected by the procedure shown in FIG. 6B.
- the sign of the differential value in the time direction of analysis value 1 (R 31 (t) -R 1 (t-1) or R 32 (t) -R 1 ( t-1)) and the differential value of the analysis value 2 in the time direction (R 2 (t) -R 2 (t-1)) are equal in sign. Therefore, out of the two solutions, the one in which the analysis value 2 and the sign of such a differential value are equal is selected (steps 6021 and 6022).
- the signs of the differential values of the two solutions are equal, the selection is made by paying attention to the magnitude of the differential value. That is, the user recognizes the flicker when the luminance difference and contrast difference from the previous frame are large.
- Equation 6 the one with the smaller differential value is selected (steps 6023, 6024, 6025). If the analysis value used in the gradation process of the t-1 frame is not corrected, R 3 (t-1) in Equation 6 is equal to R 1 (t-1).
- ⁇ Fourth embodiment> Another example of a method of selecting one solution when the solution of the quadratic equation related to the analysis value 3 is a real number will be described with reference to the configuration diagram of FIG. 5 and the flowchart of FIG.
- a prediction value (prediction analysis value) of analysis value 1 of the current frame is acquired from analysis values 1 obtained from a plurality of past frames, and a solution is selected based on the prediction analysis value. Therefore, in the functional configuration for realizing the X-ray image processing method shown in FIG. 5, an analysis value predicting unit 507 is added to the configuration of FIG.
- the analysis value prediction unit 507 receives the analysis value 1 in the past frame as an input, and outputs a prediction analysis value obtained as a result of predicting the analysis value 1 in the current frame.
- Steps 701 to 709 are the same as steps 301 to 309, and will not be described.
- Analyzing value correcting unit 105 the analysis value of 1 and the variance V 1 (t), and inputs the variance V 2 of the analysis value 2 (t), failure result outputs an analysis value 3 that has been corrected (step 710) .
- the analysis value 1 is corrected using Equation 3 so that the dispersion ratio falls within the threshold thF.
- two solutions R 31 (t) and R 32 (t) are obtained.
- the analysis value prediction unit 507 reads the analysis value 1 and the value of the analysis value 1 calculated in the past frame from the external storage unit 211, and outputs R 3 (t) ′ as the prediction analysis value (step 711). .
- a method of calculating the predicted analysis value for example, a method of extrapolating the analysis value 1 in the time direction using spline interpolation or the like can be cited.
- the analysis value correction unit 105 selects one of the two analysis values obtained in Step 710 using Equation 7 using the predicted analysis value R 3 (t) ′ (Step 712).
- the distance between each of R 31 (t) and R 32 (t) and R 3 (t) ′ is calculated (step 6031), and the shorter of the calculated distances Are determined to be more optimal solutions (steps 6032 and 6033).
- the analysis value 1 calculated as a value representing the subject area is used.
- the present invention is not limited to this.
- the above-described processing may be performed by defining the result of performing flicker countermeasure processing or the like on the calculated value representing the subject area as the analysis value 1.
- the processing of each of the above embodiments can be operated as a suppression method for interpolating flicker that cannot be suppressed by existing flicker countermeasure processing.
- Embodiments The present invention is also realized by executing the following processing. That is, software (program) that realizes the functions of the above-described embodiments is supplied to a system or apparatus via a network or various storage media, and a computer (or CPU, MPU, or the like) of the system or apparatus reads the program. It is a process to be executed.
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Geometry (AREA)
Abstract
X線センサにより連続して得られるフレームを画像処理するX線画像処理装置は、連続して得られるフレームの各々を解析して第1の解析値と第2の解析値を取得し、連続して得られるフレームの第1および第2の解析値の時間方向の統計量を第1および第2の統計量として算出する。X線画像処理装置は、算出された第1および第2の統計量に基づいて、第1の解析値を補正するか否かを判定し、補正すると判定された場合に、第1の統計量に基づいて第1の解析値を補正する。
Description
本発明は、X線画像を階調処理をするX線画像処理装置および方法、並びにX線撮像装置、プログラムに関する。
一般に、X線透視システムは、センサから出力されるX線画像の一定画質を保持するため、センサに到達しているX線の線量を監視し、その値が一定になるように入射線量を制御する。センサへの到達線量が少なければ入射線量を多くする命令が、センサへの到達線量が多ければ入射線量を少なくする命令が発行され、X線照射量が制御される。本明細書ではこのような線量の制御を自動線量制御(以下、ADC)と称する。
X線画像の連続撮影によって同じ被写体を透視し続ける場合、ADCによってX線照射量を最適化した後は、一定の線量で透視しつづけることになる。このような状態を「ADCが安定した」という。ADCが安定している状態ではユーザが観察するX線画像の階調は変化することなく、輝度のふらつきのない画像が表示されることになる。
しかしながら、X線発生装置は完全に一定の線量を出力することはできないため、ある程度のふらつきをもった出力になるのが一般的である。その結果、同じ被写体を透視し続けているのにも関わらず、ユーザが観察するX線画像に輝度のふらつきが発生してしまうことがある。そのためADCによって最適化されて得られた画像に対して、さらに画像処理を施すことによってX線画像を最適化し、輝度のふらつきを抑制することが行われている。たとえば、入力された画像を解析して、関心領域が常に一定の輝度値になるように階調変換カーブを作成し、作成された階調変換カーブを用いて画像を処理した結果が最終画像として出力される。
ところで、X線発生装置のふらつきを抑えた最適な階調変換カーブを生成するには、X線画像がどのような状態であるか画像解析をする必要がある。このような画像解析では、X線画像内のコリメーション領域や、素抜け領域(X線が直接入射される領域)を抽出することにより被写体領域を抽出し、抽出した被写体領域が最適な階調になるための解析値を算出する。そして、このような解析の精度を高めることで、より最適な最終画像を出力することができる。
しかしながら、解析の精度を高めれば高めるほど処理が困難になるため、画像解析が失敗する確率が高くなる。画像解析が失敗すると、最適な階調変換カーブを得ることができないため、輝度が不適切な画像が出力される。同じ被写体を透視し続けている場合、画像解析に失敗したフレームは輝度が不適切な状態になるため、ユーザにはフリッカが発生したように感じられてしまう。
そのようなフリッカの影響を補正する方法として、画像解析によって算出された解析値を時間軸方向で平滑化することで失敗した解析値の影響を低減する方法が挙げられる(特許文献1)。また別の方法として、画像を解析して解析値を算出し、この解析値と時間軸方向で平滑化した解析値との差分値を解析値に加算する方法がある(特許文献2)。時間方向の平滑化から安定するための値を算出し、その値に合わせ込むことでフリッカを取り除いている。
特許文献1の方法では、平滑化によりフリッカの影響を低減することはできるが、失敗量が大きい場合は低減効果が弱くなり、フリッカが発生してしまう。特許文献2の方法においても、X線発生装置のふらつきがあるため解析値もふらつくことになる。また、画像解析が失敗と判断された値に対し平均値を補正値として使用すると、X線発生装置のふらつきが考慮できなくなり、逆に輝度がふらついてしまう。
本発明の一態様によるX線画像処理装置は以下の構成を備える。すなわち、
X線センサにより連続して得られるフレームを画像処理するX線画像処理装置であって、
前記連続して得られるフレームの各々を解析して第1の解析値と第2の解析値を得る解析手段と、
前記連続して得られるフレームの前記第1の解析値と第2の解析値の時間方向の統計量を第1の統計量と第2の統計量として算出する算出手段と、
前記第1の統計量と第2の統計量に基づいて、前記第1の解析値を補正するか否かを判定する判定手段と、
前記判定手段により補正すると判定された場合に、前記第1の統計量に基づいて前記第1の解析値を補正する補正手段と、を備える。
X線センサにより連続して得られるフレームを画像処理するX線画像処理装置であって、
前記連続して得られるフレームの各々を解析して第1の解析値と第2の解析値を得る解析手段と、
前記連続して得られるフレームの前記第1の解析値と第2の解析値の時間方向の統計量を第1の統計量と第2の統計量として算出する算出手段と、
前記第1の統計量と第2の統計量に基づいて、前記第1の解析値を補正するか否かを判定する判定手段と、
前記判定手段により補正すると判定された場合に、前記第1の統計量に基づいて前記第1の解析値を補正する補正手段と、を備える。
本発明によれば、X線の照射量にふらつきが生じても効果的にフリッカを抑制することが可能になる。
本発明のその他の特徴及び利点は、添付図面を参照とした以下の説明により明らかになるであろう。なお、添付図面においては、同じ若しくは同様の構成には、同じ参照番号を付す。
添付図面は明細書に含まれ、その一部を構成し、本発明の実施の形態を示し、その記述と共に本発明の原理を説明するために用いられる。
実施形態によるX線画像処理装置の機能構成例を示すブロック図。
X線画像処理装置のハードウェア構成例を示すブロック図。
階調処理のための解析値の取得処理を示すフローチャート。
解析値保存を説明する図。
解析値保存を説明する図。
解析値保存を説明する図。
第四実施形態によるX線画像処理装置の機能構成例を示すブロック図。
解の選択手順を説明するフローチャート。
解の選択手順を説明するフローチャート。
解の選択手順を説明するフローチャート。
第四実施形態の、階調処理のための解析値の取得処理を示すフローチャート。
以下、添付の図面を参照して本発明の好適な実施形態の例を説明する。
<第一実施形態>
第一実施形態によるX線撮像装置におけるX線画像処理方法を実現するための機能構成例を図1に示す。本実施形態によれば、X線の照射量にふらつきが生じても効果的にフリッカを抑制することが可能になる。本実施形態のX線画像処理方法は、解析部101、線量状態判定部102、統計量算出部103、解析値比較部104、解析値補正部105、階調処理部106を有するX線画像処理装置により実現される。
第一実施形態によるX線撮像装置におけるX線画像処理方法を実現するための機能構成例を図1に示す。本実施形態によれば、X線の照射量にふらつきが生じても効果的にフリッカを抑制することが可能になる。本実施形態のX線画像処理方法は、解析部101、線量状態判定部102、統計量算出部103、解析値比較部104、解析値補正部105、階調処理部106を有するX線画像処理装置により実現される。
解析部101は、X線センサにて取得されたX線画像に所定の前処理を施した画像(以下、前処理済み画像)を入力とし、階調処理で利用するための解析値を出力する。前処理とは、オフセット補正、ゲイン補正、欠陥補正などのセンサの特性を補正する処理であり、周辺画素との相関関係が保たれている状態にする処理である。線量状態判定部102は、解析部101から出力された解析値を入力とし、ADCが安定しているか否かの判定結果を出力する。統計量算出部103は、解析部101からの解析値と、線量状態判定部102からの判定結果を入力とし、解析値の時間方向の統計量を計算する。解析値比較部104は、統計量算出部103で算出された時間方向の統計量を入力とし、解析値が適切な値であるか否かを判定する。解析値補正部105は、解析値比較部104により解析値が適切な値ではないと判定された場合に、解析値を適切な値に補正した結果を出力する。階調処理部106は、解析部101から出力された解析値または解析値補正部105から出力された補正された解析値のいずれかを用いて前処理済み画像を階調処理する。なお、階調処理部106は、補正後の解析値を用いるか否かを解析値比較部104の判定結果に基づいて決定する。
図1の構成を、コンピュータ等の情報処理装置(以下、コントロールPCという)を使って実現する場合、そのハードウェア構成例は図2のようになる。コントロールPC201とX線検出パネルであるX線センサ202、X線発生装置210がギガビットイーサ222でつながっている。なお、信号線はギガビットイーサに限られるものではなく、CAN(Controller Area Network)や光ファイバーなどであってもよい。X線センサ202はX線量に応じた電気信号を出力し、出力された電気信号は、たとえばX線画像として外部記憶部211に格納される。
ギガビットイーサ222には、外部表示部209、外部記憶部211、ネットワークインタフェース部212が接続されている。コントロールPC201には例えば、バス221に対して、CPU(中央演算装置)203、RAM(Random Access Memory)204、ROM(Read Only Memory)205、入力部206、表示部207、記憶部208が接続される。このコントロールPC201を介して、X線センサ202や外部表示部209などにコマンドが送られる。コントロールPC201では、撮影モードごとの処理内容がソフトウェアモジュールとして記憶部208に格納され、そのソフトウェアモジュールがRAM204に読み込まれ、CPU203により実行される。図1に示した各構成は、ソフトウェアモジュールとして記憶部208に格納されている。もちろん図1に示した各構成を専用の画像処理ボードとして実装してもよく、目的に応じて最適な実装を行うようにすればよい。以下、第一実施形態によるX線画像処理方法の詳細について説明する。
本実施形態では、X線撮像装置としてX線透視システムを例示する。図1の構成図と図3のフローチャートを参照して、第一実施形態によるX線透視システムにおけるX線画像処理を説明する。本実施形態のX線画像処理では、X線センサ202より連続して得られるフレーム(前処理済みのX線画像)を画像処理(本実施形態では階調処理)する。解析部101は連続して得られるフレームの各々を解析して、後述の解析値1や解析値2を含む複数の解析値を出力する。なお、前処理済みの画像は、たとえば、外部記憶部211に格納されており、解析部101や階調処理部106により読み出されるものとする。
本実施形態の解析部101は、1つのフレームの異なる領域について解析処理を行うことにより解析値1と解析値2を算出する。まず、解析部101は、前処理済み画像から被写体領域の画素値を代表する値を解析値1として出力する(ステップ301)。被写体領域の抽出には種々の方法があり、公知のいかなる手法が用いられてもよい。たとえば、画像からコリメーション領域と素抜け領域を抽出することにより被写体領域を抽出することができる。
コリメーション領域の抽出には、例えば特許文献3に記載された方法がある。この方法では、矩形照射野の輪郭の隣り合う2辺に沿ってX軸とY軸を設定し、画像データを設定したX軸方向およびY軸方向に加算集計する。ここで、照射野内の加算集計データはX線が殆ど照射されない照射野外に比べ高くなる。よって、X軸方向の加算集計データが所定の閾値TH以上となるY軸上の位置と、Y軸方向の加算集計データが所定の閾値TH以上となるX軸上の位置を算出する。そして、算出したY軸上の位置にあるX軸方向の直線と、X軸上の位置にあるY軸方向の直線で囲まれる矩形領域をコリメーション領域とする。
また、素抜け領域とは照射されたX線が被写体を通過せず、直接センサに到達した領域のことである。例えば画像ヒストグラムを算出すると、素抜け領域は高画素値領域に集中する。そのためヒストグラムレンジの高画素値側から数%を除去することで、素抜け領域を抽出〈除去〉することができる。なお、コリメーション領域や素抜け領域の抽出方法はこれらに限定されるものではない。
以上のようにして得られたコリメーション領域と素抜け領域をX線画像から取り除いた残りの領域が被写体領域となる。被写体領域を代表する値としては、例えば被写体領域の平均値や中間値等を利用することができる。算出された解析値1はたとえば外部記憶部211に保存される。なお、解析部101の処理対象となる画像を前処理済み画像としているが、前処理の他にLog変換、強調処理、ノイズ抑制処理などの処理を施した後の画像を入力としてもよい。
次に、解析部101は、前処理済み画像からADCに使用される解析値2を算出する(ステップ302)。ADCに利用される解析値は、解析失敗のリスクを抑えるために、解析値1に比べて単純な方法で算出される。例えば、前処理済み画像の所定位置における所定サイズの領域の代表値を用いる。本実施形態では、前処理済み画像の中心における数インチ四方の領域の平均値を解析値2として算出する。算出された解析値1や解析値2は、たとえばフレーム番号と対応付けられて外部記憶部211に保存される。
次に、線量状態判定部102は、連続して得られるフレームの解析値2について時間方向の統計量を算出し、X線の線量が安定しているか否かを判断する。より具体的には、線量状態判定部102は、現在のフレーム(以下、現フレーム)の解析値2と、過去のフレームで算出された解析値2を外部記憶部211から入力し、これら解析値を統計処理する。この統計処理により得られる統計量に基づいてADCが安定であるか不安定であるかが判定される(ステップ303、304)。ここで統計量としては、解析値2の分散または標準偏差を用いることができる。ADCが働くときは、解析値2が変動するタイミングとなるため、解析値2の変動量を計算し、閾値処理を行う。解析値2の値をR2(t)とし、tはフレーム番号とし、変動量を表す指標として分散V2(t)を以下の数1のように算出する。なお、n2は時間方向のサンプリング数を表す。
分散V2(t)の値が閾値を超えない場合、X線制御が行われていないと判断し、ADCが安定であると判断する。ADCが不安定であると判定された場合は、解析値1が成功しているか失敗しているかを判定するのは困難であるため、解析値1の補正をしない。すなわち、階調処理部106は、補正されていない解析値1を階調処理に用いて前処理済み画像の階調処理を行う(ステップ305)。
他方、ADCが安定であると判定された場合は、統計量算出部103と解析値比較部104により解析値1の補正処理を実行するか否かを判定する。この判定のために、統計量算出部103は、連続して得られるフレームの解析値1の時間方向の統計量を算出する。すなわち、統計量算出部103は現フレームでの解析値1と、過去のフレームで算出された解析値1の結果を外部記憶部211から読み出し、現フレームでの解析値1が成功しているか失敗しているかを判定するための統計量を算出する(ステップ306)。この統計量には、解析値2の場合と同様に分散や標準偏差を用いることができ、本実施形態では、統計量として解析値1の時間方向の分散V1(t)が数2を用いて算出される。なお、数2において、解析値1の値をR1(t)とし、tはフレーム番号、n1は時間方向のサンプリング数を表す。
次に、解析値比較部104が解析値1の分散V1(t)と解析値2の分散V2(t)とを入力とし、解析値1が失敗しているかを判断する(ステップ307、308)。解析値比較部104は、この判断のためにV1(t)をV2(t)で割った分散比を算出し、算出された比を閾値と比較する。ADCが安定している状態でかつ解析値1が失敗していなければ、分散V1(t)とV2(t)はほぼ同じ値になるはずである。したがって分散比が閾値以下であるときは、解析値1は正しいと判断される。解析値1が正しいと判断されると、階調処理部106は解析値1を補正せずに用いて、前処理済み画像の階調処理を行う(ステップ305)。
ステップ308の比較の結果、解析値1が失敗である(解析値1の取得に失敗している)と判断された場合に、ステップ310による解析値1の補正を実行することも考えられる。しかしながら、サンプリングされた解析値1の中でどの解析値が失敗した解析値であるかは判らないため、単純にステップ308の結果に応じて解析値1を補正すると、補正する必要の無い解析値1を補正してしまう可能性がある。そこで本実施形態では、過去の所定数のフレームを参照することにより現フレームの解析値1の取得に失敗したかどうかを確認する(ステップ309)。この確認処理により、現フレームの解析値1に対する補正の実行要否をより正確に判定することができる。なお、ステップ309で現フレームの解析値1の失敗を確認できなかった場合は、解析値1の補正は禁止され、当該解析値1がそのまま用いられる。以下、本実施形態の確認処理について説明する。
ADCが安定している状態で分散比が閾値を超えた場合、分散V1(t)の算出に用いられた解析値に失敗した解析値が含まれていることになる。本実施形態では、対象フレームでの解析が失敗したことを明確にするために、前フレームまでの分散比が閾値以下になっていることを確認する(ステップ309)。前フレームまで分散比が閾値以下であったものが、新たに加えた現フレームでの解析値1によって分散比の閾値を超えることになれば、失敗した解析は現フレームでの解析であることになる。ここで、前フレームまでの分散比が閾値を超えている場合は、どのフレームでの解析が失敗しているのか断定できない。そのため、階調処理部106は、解析値1を補正せずに用いて前処理済み画像の階調処理を行う(ステップ305)。
前フレームまで分散比が閾値以下であった場合、解析値補正部105は、解析値1とその分散V1(t)、解析値2の分散V2(t)とを入力とし、解析値1が補正された解析値3を出力する(ステップ310)。本実施形態では、解析値補正部105は、現フレームの解析値1を含めた分散比が閾値thF以内に収まるように、数3を用いて解析値1を補正する。解析値3をR3(t)とすると、R3(t)に関する二次方程式を解くことになる。
数3によれば、解が実数のときは2つの結果R31(t)、R32(t)が得られることになる。2つ得られた結果からうち1つを選択する方法としては、たとえば、補正前の解析値1との差が最も小さくなる解を補正後の解析値として選択することが挙げられる。この処理を、図6Aのステップ6011、6012、6013に示す。すなわち、R31(t)、R32(t)のそれぞれとR1(t)との差の絶対値を算出し(ステップ6011)、R31(t)、R32(t)のうち算出された差が小さい方が選択される(数4、ステップ6012、S6013)。こうして、R31(t)、R32(t)のうちR1(t)に近い値が選択されることにより、R1(t)の値の過補正が防止される。
一方、解が虚数になる場合は、分散比を閾値に最も近づける値、すなわち分散比を最小にする値を算出し直す。分散比を最小にするには、二次方程式の最小値を求めればよい。数3に示した数式のうち3番目の数式は下に凸の形状になるため、その微分値を計算すればよいことになる。
以上のように、第一実施形態によれば、2つの異なる領域から得られた解析値1と解析値2の分散比を保持するように解析値1が補正される。そのため、解析値の変動にX線発生装置のふらつきが入っていたとしても、そのふらつき成分を保持するように補正することができる。その結果、フリッカを効果的に抑制することができる。また、フレームごとに算出される解析値2に基づいてX線量の状態の安定が確認された場合に、解析値1の補正の要否を判断するので、X線の照射量に多少のふらつきが生じても効果的にフリッカを抑制することが可能になる。その後、階調処理部106が解析値3を用いて前処理済み画像に階調処理を行う(ステップ311)。なお、階調処理部106の処理対象となる画像を前処理済み画像としているが、前処理の他に強調処理やノイズ抑制処理などの処理をした後の画像であってもよい。
なお、上記実施形態では、V1(t)とV2(t)の比を用いたが、X線量の安定(V2(t)が閾値以下であること)が確認されているので、V1(t)が所定値となるように解析値1を補正するようにしてもよい。すなわち、線量状態判定部102により解析値2の統計量に基づいてX線の線量が安定していると判定された場合に、解析値1の補正の要否の判定や解析値1の補正値の算出を解析値1の統計量に基づいて実行するようにしてもよい。
本実施形態によれば、X線センサにより連続して得られるフレームを画像処理するX線画像処理装置であって、連続して得られるフレームの各々を解析して第1の解析値と第2の解析値を得る解析手段(たとえば、解析部101)と、連続して得られるフレームの第1の解析値と第2の解析値の時間方向の統計量を第1の統計量と第2の統計量として算出する算出手段(たとえば、統計量算出部103)と、第1の統計量と第2の統計量に基づいて、第1の解析値を補正するか否かを判定する判定手段(たとえば、線量状態判定部102)と、判定手段により補正すると判定された場合に、第1の統計量に基づいて第1の解析値を補正する補正手段(たとえば、解析値補正部105)と、を備える。よって、X線の照射量にふらつきが生じても効果的にフリッカを抑制することが可能になる。
なお、X線センサにより連続して得られるフレームを画像処理するX線画像処理装置であって、連続して得られるフレームの各々を解析して解析値を得る解析手段(たとえば、解析部101)と、連続して得られるフレームの解析値の時間方向の統計量を算出する算出手段(たとえば、統計量算出部103)と、統計量に基づいて解析値を補正する補正手段(たとえば、解析値補正部105)を備える形態であっても実現することができる。
<第二実施形態>
上記第一実施形態では、現フレームで解析値1が失敗した場合の処理を説明した。第二実施形態では、前フレームで解析値1が失敗したと判断された場合の次フレームでの処理を、図1の構成図と図3のフローチャートを参照して説明する。なお、本実施形態では、外部記憶部211を、過去の所定数のフレームを保持するための保持部として用いる。また、ステップ301から306までは前フレームまでの処理と同じであり、第一実施形態と同様であるため、説明を省略する。
上記第一実施形態では、現フレームで解析値1が失敗した場合の処理を説明した。第二実施形態では、前フレームで解析値1が失敗したと判断された場合の次フレームでの処理を、図1の構成図と図3のフローチャートを参照して説明する。なお、本実施形態では、外部記憶部211を、過去の所定数のフレームを保持するための保持部として用いる。また、ステップ301から306までは前フレームまでの処理と同じであり、第一実施形態と同様であるため、説明を省略する。
解析値比較部104は解析値1の分散V1(t)と解析値2の分散V2(t)とを入力とし、解析値1が失敗しているかを判断する(ステップ307、308)。この判断をするために、解析値比較部104はV2(t)をV1(t)で割った分散比を算出し閾値と比較する。ADCが安定している状態で分散比が閾値を超えた場合、分散V1(t)を算出した解析値に失敗した解析値が含まれていることになる。対象フレームでの解析が失敗したことを明確にするために、保持部に保持されている前フレームまでの分散値が閾値以下になっていることを確認する(ステップ309)。第二実施形態では、ステップ309において前フレームまでの分散値を算出する際、解析値1が失敗したと判定されたフレームの解析値1を使用せずに分散値を算出する。すなわち、ステップ309において、解析に成功した解析値1を用いて分散値が算出される。
過去フレームで算出した解析値1の選択方法を図4A~図4Cに示す。図4A~図4Cはサンプル数が4の場合の外部記憶部211の保持部における格納状況を示す図である。保持部には、フレーム番号tにおける解析値1の値R1(t)が格納されており、過去の4つのフレームの解析値1が保持されている。図4Aはt=5のときに、解析値1が成功であると判断された場合の保持部における解析値1の保持動作を説明する図である。この場合、最も古いR1(1)の値が破棄され、新たにR1(5)の値が格納される。これはより最新の分散状態と比較したいため、古いものから破棄するようにしているためである。t=6のときはこれらの4つの値からR1(6)の判定を行うことになる。
一方、図4Bはt=5のときに、解析値1が失敗であると判断された場合の保持部における格納状況を示す図である。この場合、R1(5)の値は保持部に格納されない。失敗した解析値1を使うと、対象フレームでの解析が失敗したことを明確にすることができなくなるためである。しかし、解析値1が失敗し続けるといつまでも古い分散値が保持部に維持されてしまうため、トレンド成分などで少しずつ分散値が変わっている場合等に、正しい判定ができなくなってしまう。したがって、所定フレーム数にわたって保持部の保持内容が更新されなかった場合(最も古い解析値1のフレームから現フレームまでのフレーム数が所定フレーム数を超えた場合)は、強制的に現フレームの解析値1により保持部の保持内容を更新する。
図4Cはその閾値(上述した所定フレーム数)を9フレームとした場合の保持部における格納状況を説明する図である。R1(10)が失敗した解析値であったが、最も古いR1(1)がその閾値(9フレーム)を超えたため、強制的に破棄されR1(10)の値が格納される。このとき、前フレームまでの分散値が閾値以上になってしまうため、解析値R1の補正は行われない。このように過去フレームの解析値1が失敗している場合は、その値を使わないようにして判定を行う。
なお、本実施形態では古い解析値1を破棄していったが、メモリに余裕がある場合は破棄をせずに保持していてもかまわない。選択方法として新しいものから選択する方法でもよい。ステップ310からステップ311の処理は第一実施形態と同様である。また、ステップ307の分散比の算出では、失敗した解析値1の除外を行っていないが、ステップ309と同様に、失敗した解析値1を除外して分散比を算出するようにしてもよい。
また、上記第二実施形態では、現フレームの解析値1が失敗であるか否かを判定する処理(ステップ309)において、過去に失敗と判定された解析値1を用いないようにしたが、これに限られるものではない。たとえば、ステップ306における統計量算出部103において、過去のフレームの解析値1のうち、補正が必要と判定された解析値(失敗と判定された解析値1)を除外して統計量の算出を行うようにしてもよい。
<第三実施形態>
第一実施形態では、解析値3に関する二次方程式の解が実数のときに、解を1つに選択する方法として図6Aに示す方法を示したが、これに限られるものではない。第三実施形態では、そのような解の選択方法の変形例を示す。第三実施形態では以下のような処理により解が選択される。
(1)複数の解のうち、解析値1の時間方向の微分値の符号を解析値2の時間方向の微分値の符号と一致させる解を補正後の解析値1として選択する。
(2)解析値1の時間方向の微分値が同じである複数解のうち、解析値1の時間方向の微分値と解析値2の時間方向の微分値との差が最も小さくなる解を補正後の解析値1として選択する。以下、詳述する。
第一実施形態では、解析値3に関する二次方程式の解が実数のときに、解を1つに選択する方法として図6Aに示す方法を示したが、これに限られるものではない。第三実施形態では、そのような解の選択方法の変形例を示す。第三実施形態では以下のような処理により解が選択される。
(1)複数の解のうち、解析値1の時間方向の微分値の符号を解析値2の時間方向の微分値の符号と一致させる解を補正後の解析値1として選択する。
(2)解析値1の時間方向の微分値が同じである複数解のうち、解析値1の時間方向の微分値と解析値2の時間方向の微分値との差が最も小さくなる解を補正後の解析値1として選択する。以下、詳述する。
ステップ310において、解析値補正部105が、解析値1とその分散V1(t)、解析値2の分散V2(t)とを入力とし、失敗と判定された解析値1が補正された解析値3が生成される。解析値1の補正方法では、分散比が閾値thF以内に収まるように、第一実施形態で示したように、数3が用いられる。解が実数になるときは2つの解(R31(t)、R32(t))が得られるため、2つのうち1つを選択する必要がある。本実施形態では、図6Bに示される手順で解が選択される。
分散の主原因がX線発生装置のばらつきにあるため、解析値1の時間方向の微分値の符号(R31(t)-R1(t-1)あるいはR32(t)-R1(t-1))と解析値2の時間方向の微分値(R2(t)-R2(t-1))の符号とは等しくなる。そこで、2つの解のうち、解析値2とそのような微分値の符号が等しくなる方が選択される(ステップ6021、6022)。2つの解の微分値の符号が等しい場合には、微分値の大きさに着目して選択が行われる。すなわち、ユーザがフリッカとして認識するのは前フレームとの輝度差、コントラスト差が大きい時である。したがって、数6に示されるように微分値の小さい方が選択される(ステップ6023、6024、6025)。
なお、t-1フレームの階調処理で使われた解析値が補正されたものでなければ、数6のR3(t-1)はR1(t-1)に等しくなる。
<第四実施形態>
解析値3に関する二次方程式の解が実数のときに、解を1つに選択する方法のさらに別の例を、図5の構成図と図7のフローチャートを参照して説明する。第四実施形態では、過去の複数フレームから得られた解析値1から、現フレームの解析値1の予測値(予測解析値)を取得し、この予測解析値に基づいて解の選択を行う。したがって、図5に示されるX線画像処理方法を実現する機能構成では、図1の構成に対して解析値予測部507が追加されている。解析値予測部507は過去フレームにおける解析値1を入力とし、現フレームにおける解析値1を予測した結果出る予測解析値を出力とする。ステップ701から709まではステップ301から309と同じになるため、省略する。
解析値3に関する二次方程式の解が実数のときに、解を1つに選択する方法のさらに別の例を、図5の構成図と図7のフローチャートを参照して説明する。第四実施形態では、過去の複数フレームから得られた解析値1から、現フレームの解析値1の予測値(予測解析値)を取得し、この予測解析値に基づいて解の選択を行う。したがって、図5に示されるX線画像処理方法を実現する機能構成では、図1の構成に対して解析値予測部507が追加されている。解析値予測部507は過去フレームにおける解析値1を入力とし、現フレームにおける解析値1を予測した結果出る予測解析値を出力とする。ステップ701から709まではステップ301から309と同じになるため、省略する。
解析値補正部105は、解析値1とその分散V1(t)、解析値2の分散V2(t)とを入力とし、失敗結果が補正された解析値3を出力する(ステップ710)。ここで、解析値1は分散比が閾値thF以内に収まるように数3を用いて補正される。そして、解が実数になるときは2つの解(R31(t)、R32(t))が得られる。
解析値予測部507は、解析値1と、過去のフレームで算出された解析値1の値を外部記憶部211から読み出して、予測解析値としてR3(t)’を出力する(ステップ711)。予測解析値を算出する方法としては、例えば、スプライン補間などを用いて、解析値1を時間方向に外挿する方法が挙げられる。本実施形態では、解析値補正部105はステップ710で得られた2つの解析値のうちの1つを、予測解析値R3(t)’を用いて数7により選択する(ステップ712)。
すなわち、図6Cのフローチャートに示されるように、R31(t)とR32(t)のそれぞれとR3(t)’との距離が計算され(ステップ6031)、計算された距離の近い方がより最適な解であると判定される(ステップ6032、6033)。
<第五実施形態>
なお、上記第一実施形態から第四実施形態では、被写体領域を代表する値として計算された解析値1を利用したが、これに限られるものではない。たとえば、計算された被写体領域を代表する値に対してフリッカ対策処理などを施した結果を解析値1と定義して、上述した処理を行うようにしてもよい。この場合は、既存のフリッカ対策処理で抑制しきれなかったフリッカを補間する抑制方法として、上記各実施形態の処理を動作をさせることが可能である。
なお、上記第一実施形態から第四実施形態では、被写体領域を代表する値として計算された解析値1を利用したが、これに限られるものではない。たとえば、計算された被写体領域を代表する値に対してフリッカ対策処理などを施した結果を解析値1と定義して、上述した処理を行うようにしてもよい。この場合は、既存のフリッカ対策処理で抑制しきれなかったフリッカを補間する抑制方法として、上記各実施形態の処理を動作をさせることが可能である。
他の実施形態
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ(またはCPUやMPU等)がプログラムを読み出して実行する処理である。
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ(またはCPUやMPU等)がプログラムを読み出して実行する処理である。
本発明は上記実施の形態に制限されるものではなく、本発明の精神及び範囲から離脱することなく、様々な変更及び変形が可能である。従って、本発明の範囲を公にするために、以下の請求項を添付する。
本願は、2013年10月10日提出の日本国特許出願特願2013-213224を基礎として優先権を主張するものであり、その記載内容の全てを、ここに援用する。
Claims (21)
- X線センサにより連続して得られるフレームを画像処理するX線画像処理装置であって、
前記連続して得られるフレームの各々を解析して第1の解析値と第2の解析値を得る解析手段と、
前記連続して得られるフレームの前記第1の解析値と第2の解析値の時間方向の統計量を第1の統計量と第2の統計量として算出する算出手段と、
前記第1の統計量と第2の統計量に基づいて、前記第1の解析値を補正するか否かを判定する判定手段と、
前記判定手段により補正すると判定された場合に、前記第1の統計量に基づいて前記第1の解析値を補正する補正手段と、を備えることを特徴とするX線画像処理装置。 - 前記判定手段は、前記第2の統計量に基づいてX線の線量が安定していると判定された場合に、前記第1の統計量と第2の統計量に基づいて前記第1の解析値を補正するか否かを判定することを特徴とする請求項1に記載のX線画像処理装置。
- 前記解析手段は、1つのフレームの異なる領域について解析を行うことにより前記第1の解析値と第2の解析値を得ることを特徴とする請求項1または2に記載のX線画像処理装置。
- 前記解析手段は、前記フレームの被写体領域を解析して前記第1の解析値を取得し、前記フレームの所定位置における所定サイズの領域を解析して前記第2の解析値を得ることを特徴とする請求項3に記載のX線画像処理装置。
- 前記判定手段は、前記第1の統計量と前記第2の統計量の比が閾値を超える場合に、前記第1の解析値の補正が必要であると判定することを特徴とする請求項2に記載のX線画像処理装置。
- 前記補正手段は、前記比が前記閾値と等しくなるように前記第1の解析値を補正することを特徴とする請求項5に記載のX線画像処理装置。
- 前記補正手段は、前記比が前記閾値に最も近づくように前記第1の解析値を補正することを特徴とする請求項5に記載のX線画像処理装置。
- 前記算出手段は、前記第1の解析値と第2の解析値の分散または標準偏差を前記統計量として算出することを特徴とする請求項1乃至7のいずれか1項に記載のX線画像処理装置。
- 前記算出手段は、過去のフレームの解析値のうち、前記判定手段により補正が必要と判定された解析値を除外して前記統計量の算出を行うことを特徴とする請求項1乃至8のいずれか1項に記載のX線画像処理装置。
- 前記補正手段は、前記比を前記閾値と等しくするように得られた複数の解析値のうち、補正前の前記第1の解析値との差が最も小さい解析値を補正後の前記第1の解析値として出力することを特徴とする請求項6に記載のX線画像処理装置。
- 過去の複数のフレームから得られた前記第1の解析値から、予測解析値を取得する予測手段を更に備え、
前記補正手段は、前記比を前記閾値と等しくするように得られた複数の解析値のうち、前記予測手段が予測した解析値との差が最も小さい解析値を補正後の前記第1の解析値として出力することを特徴とする請求項5に記載のX線画像処理装置。 - 前記補正手段は、前記比を前記閾値と等しくするように得られた複数の解析値のうち、前記第1の解析値の時間方向の微分値の符号と前記第2の解析値の時間方向の微分値の符号とを一致させる解析値を補正後の前記第1の解析値として選択することを特徴とする請求項5に記載のX線画像処理装置。
- 前記補正手段は、前記比を前記閾値と等しくするように得られた複数の解析値のうち、前記第1の解析値の時間方向の微分値と前記第2の解析値の時間方向の微分値との差が最も小さくなる解析値を補正後の前記第1の解析値として選択することを特徴とする請求項12に記載のX線画像処理装置。
- 過去の所定数のフレームにおける前記第1の解析値を保持する保持手段と、
前記所定数のフレームの前記第1の解析値の時間方向の統計量に基づいて、現在のフレームで前記第1の解析値の取得に失敗したことを確認する確認手段と、を更に備え、
前記補正手段は、前記確認手段により現在のフレームで前記第1の解析値の取得に失敗したことを確認できない場合に、前記第1の解析値の補正を禁止することを特徴とする請求項1乃至13のいずれか1項に記載のX線画像処理装置。 - 前記保持手段は、
保持されている最も古いフレームに対応する第1の解析値を破棄して新たに計算された前記第1の解析値を保持し、
前記補正手段による補正が必要とされた第1の解析値を保持しないことを特徴とする請求項14に記載のX線画像処理装置。 - 前記保持手段は、最も古い第1の解析値に対応するフレームと現在のフレームとの間のフレーム数が所定フレーム数を超える場合は、該現在のフレームの第1の解析値を、補正が必要であったか否かに関わらず保持することを特徴とする請求項15に記載のX線画像処理装置。
- X線センサにより連続して得られるフレームを画像処理するX線画像処理装置であって、
前記連続して得られるフレームの各々を解析して解析値を得る解析手段と、
前記連続して得られるフレームの前記解析値の時間方向の統計量を算出する算出手段と、
前記統計量に基づいて前記解析値を補正する補正手段と、を備えることを特徴とするX線画像処理装置。 - X線センサにより連続して得られるフレームを画像処理するX線画像処理方法であって、
前記連続して得られるフレームの各々を解析して第1の解析値と第2の解析値を得る解析工程と、
前記連続して得られるフレームの前記第1の解析値と第2の解析値の時間方向の統計量を第1の統計量と第2の統計量として算出する算出工程と、
前記第1の統計量と第2の統計量に基づいて、前記第1の解析値を補正するか否かを判定する判定工程と、
前記判定工程で補正すると判定された場合に、前記第1の統計量に基づいて前記第1の解析値を補正する補正工程と、を有することを特徴とするX線画像処理方法。 - X線センサにより連続して得られるフレームを画像処理するX線画像処理方法であって、
前記連続して得られるフレームの各々を解析して解析値を得る解析工程と、
前記連続して得られるフレームの前記解析値の時間方向の統計量を算出する算出工程と、
前記統計量に基づいて前記解析値を補正する補正工程と、を有することを特徴とするX線画像処理方法。 - コンピュータに、請求項18または19に記載の画像処理方法の各工程を実行させるためのプログラム。
- 請求項1乃至17のいずれか1項に記載されたX線画像処理装置と、
X線量に応じた電気信号を出力する前記X線センサと、
X線を発生するX線発生装置と、を備えることを特徴とするX線撮像装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/025,450 US10016176B2 (en) | 2013-10-10 | 2014-09-29 | X-ray image processing apparatus and method, and X-ray imaging apparatus |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013-213224 | 2013-10-10 | ||
JP2013213224A JP6254813B2 (ja) | 2013-10-10 | 2013-10-10 | X線画像処理装置および方法、並びにx線撮像装置、プログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2015052890A1 true WO2015052890A1 (ja) | 2015-04-16 |
Family
ID=52812726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2014/004962 WO2015052890A1 (ja) | 2013-10-10 | 2014-09-29 | X線画像処理装置および方法、並びにx線撮像装置、プログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US10016176B2 (ja) |
JP (1) | JP6254813B2 (ja) |
WO (1) | WO2015052890A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6957170B2 (ja) | 2017-03-08 | 2021-11-02 | キヤノン株式会社 | 放射線撮影装置、放射線撮影システム、放射線撮影方法、及びプログラム |
JP7212530B2 (ja) | 2019-01-24 | 2023-01-25 | キヤノン株式会社 | 放射線撮影システム、制御装置および制御方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6392336A (ja) * | 1986-10-08 | 1988-04-22 | 株式会社東芝 | X線診断装置 |
JPH04139587A (ja) * | 1990-09-30 | 1992-05-13 | Shimadzu Corp | 画像処理装置 |
JPH06268915A (ja) * | 1993-03-12 | 1994-09-22 | Hitachi Medical Corp | ディジタル画像処理装置 |
JPH07240886A (ja) * | 1994-02-28 | 1995-09-12 | Shimadzu Corp | X線テレビジョン装置 |
JPH0877329A (ja) * | 1994-09-02 | 1996-03-22 | Konica Corp | 時系列処理画像の表示装置 |
JP2009033293A (ja) * | 2007-07-25 | 2009-02-12 | Rohm Co Ltd | 画像処理回路、半導体装置、画像処理装置 |
JP2011125462A (ja) * | 2009-12-16 | 2011-06-30 | Canon Inc | X線画像処理装置、x線画像処理方法、及びコンピュータプログラム |
JP2013127773A (ja) * | 2011-10-07 | 2013-06-27 | Zakrytoe Akcionernoe Obshchestvo (Impul's) | デジタルx線フレームシリーズにおけるノイズ低減方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63133760A (ja) | 1986-11-25 | 1988-06-06 | Fuji Photo Film Co Ltd | 照射野認識方法および画像処理条件決定方法 |
JPH01212983A (ja) | 1988-02-19 | 1989-08-25 | Shimadzu Corp | ディジタルx線tv装置 |
US5675624A (en) * | 1995-12-14 | 1997-10-07 | General Electric Company | Adaptive x-ray brightness and display control for a medical imaging system |
NO2628146T3 (ja) * | 2010-07-01 | 2018-04-07 | ||
US8995614B2 (en) * | 2010-09-29 | 2015-03-31 | Konica Minolta Medical & Graphic, Inc. | Method for displaying medical images and medical image display system |
EP2716222A4 (en) * | 2011-05-30 | 2014-11-26 | Fujifilm Corp | METHOD AND DEVICE FOR DETECTING A BEAM EXPOSURE LEVEL AND X-RAY IMAGING SYSTEM |
-
2013
- 2013-10-10 JP JP2013213224A patent/JP6254813B2/ja active Active
-
2014
- 2014-09-29 WO PCT/JP2014/004962 patent/WO2015052890A1/ja active Application Filing
- 2014-09-29 US US15/025,450 patent/US10016176B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6392336A (ja) * | 1986-10-08 | 1988-04-22 | 株式会社東芝 | X線診断装置 |
JPH04139587A (ja) * | 1990-09-30 | 1992-05-13 | Shimadzu Corp | 画像処理装置 |
JPH06268915A (ja) * | 1993-03-12 | 1994-09-22 | Hitachi Medical Corp | ディジタル画像処理装置 |
JPH07240886A (ja) * | 1994-02-28 | 1995-09-12 | Shimadzu Corp | X線テレビジョン装置 |
JPH0877329A (ja) * | 1994-09-02 | 1996-03-22 | Konica Corp | 時系列処理画像の表示装置 |
JP2009033293A (ja) * | 2007-07-25 | 2009-02-12 | Rohm Co Ltd | 画像処理回路、半導体装置、画像処理装置 |
JP2011125462A (ja) * | 2009-12-16 | 2011-06-30 | Canon Inc | X線画像処理装置、x線画像処理方法、及びコンピュータプログラム |
JP2013127773A (ja) * | 2011-10-07 | 2013-06-27 | Zakrytoe Akcionernoe Obshchestvo (Impul's) | デジタルx線フレームシリーズにおけるノイズ低減方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6254813B2 (ja) | 2017-12-27 |
US10016176B2 (en) | 2018-07-10 |
JP2015073800A (ja) | 2015-04-20 |
US20160270752A1 (en) | 2016-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5954063B2 (ja) | 画像の霧除去方法及びシステム | |
JP5889013B2 (ja) | 画像処理装置及び画像処理方法 | |
JP6390847B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
US9025053B2 (en) | Imaging apparatus and image restoration method | |
US9922409B2 (en) | Edge emphasis in processing images based on radiation images | |
WO2011132416A1 (ja) | 撮像装置及び画像復元方法 | |
JP5091994B2 (ja) | 動きベクトル検出装置 | |
US10319083B2 (en) | Image artifact detection and correction in scenes obtained from multiple visual images | |
JP6254813B2 (ja) | X線画像処理装置および方法、並びにx線撮像装置、プログラム | |
JP6253633B2 (ja) | 画像処理装置、画像処理方法、プログラムおよびコンピュータ記録媒体 | |
JP7409606B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
US10755389B2 (en) | Image processing apparatus, image processing method, and medium | |
US20100260386A1 (en) | Image processing apparatus and control method of image processing apparatus | |
CN111429369B (zh) | 高动态范围图像生成方法、装置、电子设备及存储介质 | |
US10339637B2 (en) | Image processing apparatus, image pickup apparatus, image processing method, and non-transitory computer-readable storage medium for correcting deterioration of image | |
JP5674659B2 (ja) | 画像処理装置、x線画像診断装置 | |
JP6121013B2 (ja) | 画像処理装置及び画像処理方法 | |
JP7409607B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
JP6250988B2 (ja) | 画像処理方法および装置 | |
US20230070520A1 (en) | Image processing apparatus, radiation imaging system, image processing method, and computer-readable medium | |
CN112381836B (zh) | 图像处理方法和装置、计算机可读存储介质、电子设备 | |
EP3719740B1 (en) | Image processing device, image processing method, and image processing program | |
KR101796389B1 (ko) | 문서 캡처 영상의 영상 보정 방법 | |
JP2013195402A (ja) | 信号分析装置、方法及びプログラム | |
JP2021128497A (ja) | 画像処理装置、画像処理方法およびプログラム |
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: 14852447 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 15025450 Country of ref document: US |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 14852447 Country of ref document: EP Kind code of ref document: A1 |