WO2015029695A1 - 画像処理装置、画像処理方法及びプログラム - Google Patents

画像処理装置、画像処理方法及びプログラム Download PDF

Info

Publication number
WO2015029695A1
WO2015029695A1 PCT/JP2014/070338 JP2014070338W WO2015029695A1 WO 2015029695 A1 WO2015029695 A1 WO 2015029695A1 JP 2014070338 W JP2014070338 W JP 2014070338W WO 2015029695 A1 WO2015029695 A1 WO 2015029695A1
Authority
WO
WIPO (PCT)
Prior art keywords
data loss
data
image data
loss rate
fitting
Prior art date
Application number
PCT/JP2014/070338
Other languages
English (en)
French (fr)
Inventor
野田　剛司
Original Assignee
キヤノン株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by キヤノン株式会社 filed Critical キヤノン株式会社
Priority to US14/913,667 priority Critical patent/US10311549B2/en
Publication of WO2015029695A1 publication Critical patent/WO2015029695A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/70—Denoising; Smoothing
    • 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
    • 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/24—Measuring radiation intensity with semiconductor detectors
    • G01T1/247—Detector read-out circuitry
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/60—Noise processing, e.g. detecting, correcting, reducing or removing noise
    • H04N25/63—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to dark current
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00—Details of television systems
    • H04N5/30—Transforming light or analogous information into electric information
    • H04N5/32—Transforming X-rays
    • 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4233—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
    • 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 technique for correcting image data captured by detecting radiation.
  • FPDs flat-panel detectors
  • Such an X-ray imaging apparatus is integrally formed with a conventional modality and is configured as a dedicated apparatus.
  • an imaging medium such as an existing modality film or imaging plate with an FPD. is there.
  • Patent Document 1 discloses a technique in which an interface is not provided between the X-ray generator and the FPD, and X-ray irradiation is detected on the FPD side to start a charge accumulation operation.
  • FPD converts photons into electrons with a photodiode
  • dark charges are actually generated by thermal excitation or leakage current even when there are no photons.
  • this dark charge is accumulated, it is added to the actual X-ray charge, resulting in unevenness in the data and artifacts.
  • the amount of X-rays that can be accumulated decreases due to the accumulation of dark charges. For this reason, in the FPD, a reset operation for continuously turning on the TFT and discharging dark charges is performed before X-ray imaging.
  • Patent Document 2 discloses a technique for performing a reset operation by performing a reset operation at the time of X-ray detection for each line and switching an even line and an odd line for each frame.
  • the artifacts can be corrected by interpolation processing.
  • the grid is one of basic tools in X-ray photography, and has a function of blocking X-rays scattered inside the subject and suppressing blurring of the X-ray image data.
  • imaging using a grid in the technique disclosed in Patent Document 2, if the direction of the grid line and the direction of the scanning line coincide with each other, the grid stripe and an artifact generated at the time of X-ray detection interfere with each other.
  • this X-ray image data is corrected by interpolation processing, new artifacts may occur or grid stripes may be partially erased.
  • Patent Document 3 discloses a technique for removing grid stripes by model fitting. However, in the technique disclosed in Patent Document 3, it is not assumed that there is an artifact that occurs during X-ray detection, and if the grid stripe is partially erased, the grid stripe is appropriately removed. There are things that cannot be done.
  • an object of the present invention is to appropriately correct data deficiencies in image data generated during radiation detection even when radiation imaging is performed using a grid.
  • An image processing apparatus includes an acquisition unit that acquires image data captured by detecting radiation, a calculation unit that calculates a data loss rate in the image data, and a model function that represents data loss. Fitting means for fitting to a rate, and correction means for correcting data loss in the image data based on the model function fitted by the fitting means.
  • FIG. 1 is a diagram illustrating a configuration of an X-ray imaging apparatus according to an embodiment of the present invention. It is a flowchart which shows the process of the X-ray imaging apparatus which concerns on embodiment of this invention. It is a figure which shows the circuit structure of FPD. 3 is a timing chart showing driving of the FPD in the first embodiment. It is a figure which shows the example of the X-ray image data acquired in 1st Embodiment. It is a figure for demonstrating an X-ray waveform model. It is a figure which shows the data loss rate calculated from the X-ray-image data image
  • FIG. 6 is a timing chart illustrating driving of an FPD according to a second embodiment. It is a figure which shows the example of the X-ray image data acquired in 2nd Embodiment. It is a figure which shows the data loss rate calculated in the data loss rate calculation part in 2nd Embodiment, and the model function fitted with the said data loss rate.
  • a sensor using a PIN photodiode is applied, but a sensor using a MIS photodiode may be applied.
  • a direct FPD that directly converts X-rays into electrons may also be applied.
  • FIG. 1 is a diagram showing a configuration of an X-ray imaging apparatus 100 according to the first embodiment of the present invention.
  • an X-ray tube 101 irradiates a subject 103 with X-rays.
  • the X-ray generator 104 applies a high voltage pulse to the X-ray tube 101 when the exposure switch is pressed to generate X-rays.
  • the FPD 102 Based on the control of the FPD control unit 105, the FPD 102 converts X-rays that have passed through the subject 103 into visible light using a phosphor, and detects an electrical signal using a photodiode.
  • the detected electrical signal is AD converted into a digital signal and output to the FPD control unit 105 as X-ray image data.
  • the grid 114 is disposed in front of the FPD 102, blocks X-rays scattered inside the subject, and suppresses blurring of the X-ray image data.
  • the FPD control unit 105 is connected to the image processing unit 109, the image storage unit 108, the monitor 106, and the operation unit 107.
  • the FPD control unit 105 includes one or a plurality of computers.
  • the computer includes, for example, a main control unit such as a CPU, and a storage unit such as a ROM (Read Only Memory) and a RAM (Random Access Memory).
  • the computer also includes a graphic control unit such as a GPU (Graphics Processing Unit), a communication unit such as a network card, and an input / output unit such as a keyboard, a display, or a touch panel.
  • a bus or the like is controlled by the main control unit executing a program stored in the storage unit.
  • the monitor 106 displays the X-ray image data input from the FPD control unit 105.
  • the operation unit 107 inputs instructions for the image processing unit 109 and the FPD 102.
  • the image storage unit 108 stores the X-ray image data output from the FPD control unit 105 and the X-ray image data processed by the image processing unit 109.
  • the image processing unit 109 includes a function calculation unit 110, a data loss rate calculation unit 111, a fitting unit 112, and a correction processing unit 113, and processes the acquired X-ray image data in accordance with an instruction from the FPD control unit 105.
  • the function calculation unit 110 calculates a model function representing an artifact shape due to data loss caused by the reset operation until the FPD 102 detects the X-ray irradiation and starts the charge accumulation operation.
  • the data loss rate calculation unit 111 calculates the data loss rate by calculating the ratio of the number of pixels in which data loss has occurred and the number of pixels in which no data loss has occurred.
  • the fitting unit 112 fits the model function calculated by the function calculating unit 110 to the data missing rate calculated by the data missing rate calculating unit 111.
  • the correction processing unit 113 performs correction processing on the X-ray image data using the model function fitted by the fitting unit 112.
  • the FPD control unit 105 acquires X-ray image data. This is performed by applying a high voltage pulse to the X-ray tube 101 by the X-ray generator 104 and irradiating the subject with X-rays.
  • a high voltage pulse to the X-ray tube 101 by the X-ray generator 104 and irradiating the subject with X-rays.
  • the distance between the FPD 102 and the X-ray tube 101 is generally set to about 100 cm to 150 cm. In this embodiment, no interface is provided between the X-ray generator 104 and the FPD 102, and the FPD 102 detects the X-ray irradiation and starts the charge accumulation operation.
  • the circuit configuration of the FPD 102 will be described with reference to FIG.
  • the FPD 102 one pixel is formed by combining a TFT 302 and a photodiode 303 on a glass substrate 301.
  • the FPD 102 is composed of one million pixels of 1000 vertical pixels ⁇ 1000 horizontal pixels. In FIG. 3, only 16 pixels out of 1 million pixels are shown for convenience of explanation. Note that the FPD 102 is an example of a radiation detection apparatus.
  • the scanning line control circuit 306 sequentially applies an ON signal to G1 to G1000 (G5 and later are not shown) to turn on the TFT 302. Thereby, charge can be read from the photodiode 303 line by line.
  • S1 to S1000 (not shown after S5) are signal lines.
  • the charges read from the photodiode 303 are transmitted through the signal lines and read by the signal detection circuit 305.
  • the signal detection circuit 305 performs holding, amplification, AD conversion, and the like of the read charge, and outputs the digital signal to the FPD control unit 105.
  • the power supply 304 supplies an operating voltage to each photodiode 303.
  • FIG. 4 is a timing chart showing driving of the FPD 102.
  • the FPD 102 stands by in a state called “X-ray detection drive” when X-rays are not irradiated.
  • the scanning line control circuit 306 has a scanning line spacing of one line such as G1, G3, G5,. After that, only the odd-numbered lines are sequentially driven, and the dark charges of the pixels connected to the odd-numbered scanning lines are read and reset.
  • the scanning line control circuit 306 sequentially scans even-numbered lines with a line interval of one line (line interval of one line or more) such as G2, G4, and G6. Drive, read and reset pixels connected to even-numbered scan lines.
  • the odd lines are reset in the odd frames
  • the even lines are reset in the even frames.
  • the signal detection circuit 305 detects X-ray irradiation when the output of the signal line exceeds a predetermined threshold.
  • the scanning line control circuit 306 turns off all TFTs 302. As a result, the charge accumulation operation starts.
  • X-rays can be detected more quickly as the threshold value is smaller.
  • appropriate values are set according to detection sensitivity and noise level.
  • step S102 the data loss rate calculation unit 111 calculates the data loss rate from the X-ray image data acquired in step S101.
  • FIG. 5 is a diagram illustrating an example of the X-ray image data acquired in step S101.
  • one rectangle in FIG. 5 is considered as one pixel.
  • the even lines of G156 to G190 are deficient in charges that should be accumulated by X-ray irradiation by the reset operation.
  • the pixel values of adjacent lines are close to each other. Therefore, in the present embodiment, as shown in the following Expression 1, by utilizing the fact that the adjacent line (the line adjacent in the y coordinate direction) is normal, the pixel of the line where the data loss has occurred and the data loss have occurred.
  • the data loss rate A T (x, y) is calculated by calculating the ratio with the non-existing pixel.
  • P (x, y) represents a pixel value at the position coordinates (x, y).
  • step S103 the fitting unit 112 fits the model function calculated by the function calculating unit 110 to the data loss rate calculated in step S102.
  • the true data loss rate is calculated by removing the high frequency components of the edges and grids from the data loss rate calculated in step S202.
  • Data loss from the start of X-ray irradiation to X-ray detection is caused by the loss of charge accumulated from the start of X-ray irradiation to the resetting of a predetermined line. Therefore, data loss is obtained by time integration of the X-ray waveform from the start of X-ray irradiation to the resetting of a predetermined line.
  • the X-ray waveform may be considered to be equal to the voltage pulse shape of the X-ray generator 104.
  • a general voltage pulse rises with a round due to electrical impedance. This waveform model is shown in Equation 2.
  • the rise of the voltage pulse can be expressed by a time constant ⁇ , which is the product of the resistance R and the capacitance C shown in FIG. 6A, and has a waveform as shown in FIG. It becomes. Therefore, the data loss can be expressed by a model function of the following formula 3 that integrates from the X-ray irradiation start time t0 to the reset time t of a predetermined line.
  • Equation 3 is unknowns determined by design values and generation conditions of the X-ray generator 104, and are obtained by fitting. Specifically, ⁇ , V, and t 0 are obtained so as to minimize the following Expression 4. As a minimization algorithm, a generally known Newton method or conjugate gradient method may be used.
  • FIG. 7 is a diagram showing a data loss rate ( ⁇ ) calculated from X-ray image data taken in a state where no subject actually exists and a model function (solid line) fitted with the data loss rate. As shown in FIG. 7, it can be seen that the model function after fitting well matches the data loss rate calculated from the X-ray image data.
  • FIG. 8 is a diagram showing a data loss rate (broken line and ⁇ ) calculated from X-ray image data photographed by inserting a grid 114 having a periodic structure close to the Nyquist frequency of the FPD 102.
  • X-ray image data is captured in a state where the grid 114 is inserted so that the grid lines are parallel to the scanning lines.
  • the data loss rate calculated from the X-ray image data captured with the grid 114 inserted so that the grid line is parallel to the scanning line is the data loss rate in FIG. It can be seen that the vibration accompanied by “beat” is added.
  • the model function after fitting has the shape of the data loss rate based on the X-ray waveform and the driving of the FPD 102. Therefore, the model function after the fitting is not affected by the grid stripes or the “beat” that have entered the data loss rate calculated from the X-ray image data, and the true model function caused by the driving of the X-ray waveform and the FPD 102. It can be seen that the data loss rate can be extracted.
  • step S104 the correction processing unit 113 corrects the X-ray image data using the model function after fitting.
  • the correction processing unit 113 uses the model function after fitting to correct data deficiency that occurs between the start of X-ray irradiation and the start of the charge accumulation operation for X-ray image data.
  • This correction processing is performed according to the following Expression 5, where the pixel value after correction is P C (x, y) and the pixel value before correction is P (x, y).
  • the correction processing unit 113 can appropriately correct data deficiency in X-ray image data that occurs during radiation detection.
  • the FPD control unit 105 acquires X-ray image data as in step S101 in FIG.
  • FIG. 9 is a timing chart showing driving of the FPD 102 in the second embodiment.
  • the reset operation in the “X-ray detection driving” state is sequentially performed as G1, G2, G3, and G4 to reset all the scanning lines in one frame. Therefore, in the second embodiment, data loss occurs continuously as shown in FIG.
  • the data loss rate calculation unit 111 calculates the data loss rate AT from the X-ray image data acquired by the FPD control unit 105 according to the following formula 6.
  • y detect is a line where X-rays are detected. That is, one on the line y detect, since the accumulation operation of the normal charge is started from the y detect +1, data loss ratio is calculated that line as a reference.
  • FIG. 11 is a diagram illustrating a data loss rate A T ( ⁇ ) calculated by the data loss rate calculation unit 111 according to the second embodiment and a model function (solid line) fitted with the data loss rate.
  • the data loss rate AT varies with the influence of the subject.
  • the fitting unit 112 fits the model function of Expression 3 with the data loss rate AT .
  • the model function shown in Expression 3 has a data loss rate shape based on the X-ray waveform and driving of the FPD 102. Therefore, the model function after fitting is not affected by the subject that has entered the data loss rate calculated from the X-ray image data, and the true data loss rate due to the driving of the X-ray waveform and the FPD 102 can be extracted. I understand that.
  • the correction processing unit 113 corrects the X-ray image data using the model function of Equation 3 (model coefficient after fitting) in which the coefficient is determined by fitting.
  • the correction processing unit 113 uses the model function of Formula 3 (model coefficient after fitting) in which the coefficient is determined by fitting, and the data loss that has occurred from the start of X-ray irradiation to the start of charge accumulation operation regarding the X-ray image data. Correct.
  • This correction is performed by the following Expression 7, where the pixel value after correction is P C (x, y) and the pixel value before correction is P (x, y).
  • the correction processing unit 113 can appropriately correct data deficiency in X-ray image data that occurs during radiation detection.
  • the present invention can also be 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)
  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Surgery (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

　グリッドを用いて放射線撮影が行われた場合であっても、放射線検出の際に生じたデータ欠損を適切に補正する。ＦＰＤ制御部１０５は、放射線を検知することで撮影された画像データを取得する。データ欠損率算出部１１１は、上記画像データにおけるデータ欠損率を算出する。フィッティング部１１２は、データ欠損を表すモデル関数を上記データ欠損率にフィッティングさせる。補正処理部１１３は、フィッティングされた上記モデル関数に基づいて、上記画像データを補正する。

Description

画像処理装置、画像処理方法及びプログラム
　本発明は、放射線を検知することで撮影された画像データを補正する技術に関するものである。
　近年、ＴＦＴアクティブマトリクス基板上に蛍光体を配置し、Ｘ線を電荷として蓄積し、それをデジタル信号に変換してＸ線画像データを提供するＦＰＤ（flat panel detector）が実用化され、盛んに用いられている。このようなＸ線画像撮影装置は、従来のモダリティと一体的に形成され、専用装置として構成されていたが、既存のモダリティのフィルムやイメージングプレート等の撮影媒体をＦＰＤに置き換えたいという要求も多くある。
　既存のモダリティの撮影媒体をＦＰＤに置き換える場合、Ｘ線発生装置とＦＰＤとの間のインタフェースの構築が困難なことがある。これに対して、特許文献１には、Ｘ線発生装置とＦＰＤとの間にインタフェースを設けず、ＦＰＤ側でＸ線照射を検出して電荷の蓄積動作を開始する技術が開示されている。
　ところで、ＦＰＤは、フォトダイオードで光子を電子に変換するが、実際は光子がない場合でも熱励起やリーク電流により暗電荷が発生する。この暗電荷が蓄積されると、実際のＸ線による電荷に加算されてデータにムラが生じ、アーチファクトが発生する。また、暗電荷が蓄積されることでＸ線を蓄積できる量が低下してしまう。このため、ＦＰＤでは、Ｘ線撮影前に絶えずＴＦＴをオンにして暗電荷を排出するリセット動作を行っている。
　しかしながら、上述したようにＸ線照射を検出して電荷の蓄積動作を開始する場合、電荷の蓄積動作を開始するまでは、ある程度のＸ線量が必要であり、その間はリセット動作が行われてしまう。このため、電荷の蓄積動作を開始するまでに照射されたＸ線による電荷は排出されてしまい、データ欠損が生じ、アーチファクトが発生する。
　これに対して、特許文献２には、Ｘ線検出時のリセット動作を１ライン毎に行い、フレーム毎に偶数ラインと奇数ラインとを入れ替えて、リセット動作を行う技術が開示されている。特許文献２に開示された技術では、Ｘ線照射から電荷の蓄積動作の開始までに生じるデータ欠損によるアーチファクトは１ライン置きになるため、補間処理によりアーチファクトを補正することができる。
　ところで、Ｘ線撮影においては、グリッドを用いた撮影が頻繁に行われる。グリッドはＸ線撮影における基本用具の一つであり、被写体内部で散乱されたＸ線を遮断し、Ｘ線画像データがボケてしまうのを抑制する働きがある。グリッドを用いた撮影において、特許文献２に開示された技術では、グリッドラインの方向と走査線の方向とが一致した場合、グリッド縞とＸ線検出時に発生したアーチファクトとが干渉してしまう。このＸ線画像データに補間処理で補正を行うと、新たなアーチファクトが発生する場合や、グリッド縞が部分的に消去されてしまうことがある。
　特許文献３には、モデルフィッティングにより、グリッド縞を除去する技術が開示されている。しかしながら、特許文献３に開示された技術では、Ｘ線検出時に発生するアーチファクトがある場合は想定されておらず、グリッド縞が部分的に消去されてしまったりする場合は、グリッド縞を適切に除去できないことがある。
特開２０１２－０７５０７７号公報 特開２０１１－２４９８９１号公報 特開２００５－５２５５３号公報
　しかしながら、従来技術では、グリッドがセットされる方向に応じてアーチファクトが発生する可能性がある。これにより、再撮影が必要になると、患者は無駄な被曝を受けることになる。
　そこで、本発明の目的は、グリッドを用いて放射線撮影が行われた場合であっても、放射線検出の際に生じた画像データにおけるデータ欠損を適切に補正することにある。
　本発明の画像処理装置は、放射線を検知することで撮影された画像データを取得する取得手段と、前記画像データにおけるデータ欠損率を算出する算出手段と、データ欠損を表すモデル関数を前記データ欠損率にフィッティングさせるフィッティング手段と、前記フィッティング手段によりフィッティングされた前記モデル関数に基づいて、前記画像データにおけるデータ欠損を補正する補正手段とを有する。
　本発明によれば、グリッドを用いて放射線撮影が行われた場合であっても、放射線検出の際に生じた画像データにおけるデータ欠損を適切に補正することが可能となる。
本発明の実施形態に係るＸ線画像撮影装置の構成を示す図である。 本発明の実施形態に係るＸ線画像撮影装置の処理を示すフローチャートである。 ＦＰＤの回路構成を示す図である。 第１の実施形態におけるＦＰＤの駆動を示すタイミングチャートである。 第１の実施形態で取得されるＸ線画像データの例を示す図である。 Ｘ線波形モデルを説明するための図である。 実際に被写体が存在しない状態で撮影されたＸ線画像データから算出されたデータ欠損率と、当該データ欠損率でフィッティングしたモデル関数（実線）とを示す図である。 ＦＰＤのナイキスト周波数に近い周期構造を持ったグリッドを挿入して撮影されたＸ線画像データから算出されたデータ欠損率を示す図である。 第２の実施形態におけるＦＰＤの駆動を示すタイミングチャートである。 第２の実施形態で取得されるＸ線画像データの例を示す図である。 第２の実施形態におけるデータ欠損率算出部で算出されたデータ欠損率と、当該データ欠損率でフィッティングしたモデル関数とを示す図である。
　以下、本発明を適用した好適な実施形態を、添付図面を参照しながら詳細に説明する。
　先ず、本発明の第１の実施形態について説明する。なお、本実施形態に係る放射線撮影装置としてのＸ線画像撮影装置では、ＰＩＮ型フォトダイオードを用いたセンサを適用しているが、ＭＩＳ型フォトダイオードを用いたセンサを適用してもよい。また、Ｘ線を直接電子に変換する直接型ＦＰＤを適用してもよい。
　図１は、本発明の第１の実施形態に係るＸ線画像撮影装置１００の構成を示す図である。図１において、Ｘ線管１０１は、被検体１０３に対してＸ線を照射する。Ｘ線発生装置１０４は、曝射スイッチの押下でＸ線管１０１に高電圧パルスを与え、Ｘ線を発生させる。ＦＰＤ１０２は、ＦＰＤ制御部１０５の制御に基づき、被検体１０３を通過したＸ線を蛍光体により可視光に変換し、フォトダイオードで電気信号を検出する。検出された電気信号は、デジタル信号にＡＤ変換され、Ｘ線画像データとしてＦＰＤ制御部１０５に出力される。グリッド１１４は、ＦＰＤ１０２の前面に配置され、被写体内部で散乱されたＸ線を遮断し、Ｘ線画像データがボケてしまうのを抑制する。
　ＦＰＤ制御部１０５は、画像処理部１０９、画像保存部１０８、モニタ１０６及び操作部１０７と接続される。また、ＦＰＤ制御部１０５には、一又は複数のコンピュータが内蔵される。コンピュータは、例えば、ＣＰＵ等の主制御部、並びに、ＲＯＭ（Read Only Memory）及びＲＡＭ（Random Access Memory）等の記憶部を備える。また、コンピュータは、ＧＰＵ（Graphics Processing Unit）等のグラフィック制御部、ネットワークカード等の通信部、及び、キーボード、ディスプレイ又はタッチパネル等の入出力部等を備える。なお、これら各構成は、バス等により接続され、主制御部が記憶部に記憶されたプログラムを実行することで制御される。
　モニタ１０６は、ＦＰＤ制御部１０５から入力したＸ線画像データを表示する。操作部１０７は、画像処理部１０９やＦＰＤ１０２に対する指示を入力する。画像保存部１０８は、ＦＰＤ制御部１０５から出力されたＸ線画像データや画像処理部１０９で処理されたＸ線画像データを保存する。
　画像処理部１０９は、関数算出部１１０、データ欠損率算出部１１１、フィッティング部１１２及び補正処理部１１３を備え、ＦＰＤ制御部１０５からの指示に従って、取得したＸ線画像データを処理する。
　関数算出部１１０は、ＦＰＤ１０２がＸ線照射を検出して電荷の蓄積動作を開始するまでに、リセット動作で生じるデータ欠損によるアーチファクト形状を表すモデル関数を算出する。データ欠損率算出部１１１は、データ欠損が生じた画素数とデータ欠損が生じていない画素数との比を算出することで、データ欠損率を算出する。フィッティング部１１２は、データ欠損率算出部１１１で算出されたデータ欠損率に関数算出部１１０で算出されたモデル関数をフィッティングする。補正処理部１１３は、フィッティング部１１２でフィッティングされたモデル関数を用いて、Ｘ線画像データに対して補正処理を行う。
　次に、図２を参照しながら、本実施形態に係るＸ線画像撮影装置１００の処理について説明する。ステップＳ１０１において、ＦＰＤ制御部１０５は、Ｘ線画像データを取得する。これは、Ｘ線発生装置１０４により高電圧パルスをＸ線管１０１に印加し、Ｘ線を被検体に照射することで行う。Ｘ線画像データの撮影条件としては任意の値を設定することが可能であるが、胸部等の撮影では、１００ｋＶ及び２ｍＡｓ程度に設定すればよい。また、ＦＰＤ１０２とＸ線管１０１との間の距離は、一般的に１００ｃｍ～１５０ｃｍ程度に設定される。本実施形態では、Ｘ線発生装置１０４とＦＰＤ１０２との間にインタフェースを設けておらず、ＦＰＤ１０２は、Ｘ線照射を検知して電荷の蓄積動作を開始する。
　ここで、図３を参照しながら、ＦＰＤ１０２の回路構成について説明する。ＦＰＤ１０２は、ガラス基板３０１上においてＴＦＴ３０２及びフォトダイオード３０３の組み合わせで一つの画素を形成している。ここでは、ＦＰＤ１０２は、縦１０００個×横１０００個の百万個の画素から構成されているものとする。なお、図３では、説明の便宜上、百万個の画素のうちの１６個の画素をのみを示している。なお、ＦＰＤ１０２は、放射線検知装置の例となる構成である。
　走査線制御回路３０６は、Ｇ１～Ｇ１０００（Ｇ５以降は図示せず）に順次オン信号を印加してＴＦＴ３０２のスイッチをオンにする。これにより、１ラインずつフォトダイオード３０３から電荷を読み出すことができる。Ｓ１～Ｓ１０００（Ｓ５以降は図示せず）は信号線であり、フォトダイオード３０３から読み出された電荷は信号線を伝達して、信号検出回路３０５で読み取られる。信号検出回路３０５は、読み取った電荷の保持、増幅及びＡＤ変換等を実行して、デジタル信号としてＦＰＤ制御部１０５に出力する。電源３０４は、各フォトダイオード３０３に動作電圧を供給する。
　図４は、ＦＰＤ１０２の駆動を示すタイミングチャートである。図４において、ＦＰＤ１０２は、Ｘ線が照射されていないとき、「Ｘ線検出駆動」と呼ばれる状態で待機する。「Ｘ線検出駆動」状態では、走査線制御回路３０６は、第１のフレームにおいて、走査線をＧ１、Ｇ３、Ｇ５、・・・というように１ラインのライン間隔（１ライン以上のライン間隔）を空けて奇数番目のラインだけを順次駆動し、奇数番目の走査線に接続された画素の暗電荷を読み取ってリセットしていく。次に、走査線制御回路３０６は、第２のフレームにおいて、走査線をＧ２、Ｇ４、Ｇ６というように１ラインのライン間隔（１ライン以上のライン間隔）を空けて偶数番目のラインだけを順次駆動し、偶数番目の走査線に接続された画素を読み取ってリセットしていく。このように、本実施形態では、奇数フレームでは奇数ラインをリセットし、偶数フレームでは偶数ラインをリセットする。
　Ｘ線が照射されると、フォトダイオード３０３により電荷が生成されるため、信号線の出力が上昇する。信号検出回路３０５は、信号線の出力が予め定められた閾値を超えたとき、Ｘ線照射を検知する。信号検出回路３０５によりＸ線照射が検知されると、走査線制御回路３０６は、全ＴＦＴ３０２をオフにする。これにより、電荷の蓄積動作が開始する。ここで、閾値が小さいほどＸ線を素早く検知することができるが、ノイズ等の誤動作が生じるため、検知感度及びノイズレベルに応じた適切な値が設定される。
　上述したＸ線検知動作では、Ｘ線照射開始とＸ線検知との間に必ずタイムラグが生じる。そして、電荷の蓄積動作に入るまでに照射されたＸ線に対応する電荷はリセット動作により排出されるため、Ｘ線画像データに寄与することができず、データ欠損によるアーチファクトが生じることになる。
　ステップＳ１０２において、データ欠損率算出部１１１は、ステップＳ１０１で取得されたＸ線画像データからデータ欠損率を算出する。図５は、ステップＳ１０１で取得されたＸ線画像データの例を示す図である。ここでは、図５における１つの矩形を１つの画素と考える。例えば、偶数フレームのＧ１５６ラインでＸ線照射が開始され、Ｇ１９０ラインでＸ線が検知されたとする。この場合、奇数ラインはリセットされないため、正常な画素値の電荷が蓄積される。他方、Ｇ１５６～Ｇ１９０の偶数ラインはリセット動作によりＸ線照射で蓄積されるはずの電荷が欠損している。
　一般的に、自然画像データは連続的に画素値が変化するため、被写体にエッジやグリッド等の高周波成分が含まれない場合、隣接ライン同士の画素値は近い値となる。従って、本実施形態では、次の式１に示すように、隣接ライン（ｙ座標方向で隣接するライン）が正常なことを利用して、データ欠損が生じたラインの画素とデータ欠損が生じていない画素との比を算出することでデータ欠損率ＡT（ｘ，ｙ）を算出する。ここで、Ｐ（ｘ，ｙ）は、位置座標（ｘ，ｙ）における画素値を表している。
Figure JPOXMLDOC01-appb-M000001
　ステップＳ１０３において、フィッティング部１１２は、ステップＳ１０２で算出されたデータ欠損率に関数算出部１１０で算出されたモデル関数をフィッティングする。これにより、ステップＳ２０２で算出されたデータ欠損率からエッジやグリッドの高周波成分を除いた真のデータ欠損率が算出される。
　ここで、モデル関数のフィッティング処理について詳細に説明する。Ｘ線照射開始からＸ線検知までのデータ欠損は、Ｘ線照射開始から所定のラインのリセットまでに蓄積された電荷が欠損することで生じる。従って、Ｘ線照射開始から所定のラインのリセットまでＸ線波形を時間積分したものがデータ欠損となる。
　Ｘ線波形は、Ｘ線発生装置１０４の電圧パルスの形状に等しいと考えてよい。一般的な電圧パルスは、電気的なインピーダンスによるなまりを持って立ち上がる。この波形モデルを式２に示す。
Figure JPOXMLDOC01-appb-M000002
　式２に示すように、電圧パルスの立ち上がりは、図６（ａ）に示す抵抗ＲとキャパシタンスＣとの積である時定数τで表現することができ、図６（ｂ）に示すような波形となる。従って、データ欠損は、これをＸ線照射開始時間ｔ０から所定のラインのリセット時間ｔまでを積分した次の式３のモデル関数で表すことができる。
Figure JPOXMLDOC01-appb-M000003
　なお、式３におけるτ、Ｖ及びｔ0は、Ｘ線発生装置１０４の設計値や発生条件により決まる未知数であり、フィッティングにより求められる。具体的には、次の式４を最小化するようにτ、Ｖ、ｔ0が求められる。最小化アルゴリズムは、一般に知られているニュートン法や共役勾配法等を用いればよい。
Figure JPOXMLDOC01-appb-M000004
　図７は、実際に被写体が存在しない状態で撮影されたＸ線画像データから算出されたデータ欠損率（□）と、当該データ欠損率でフィッティングしたモデル関数（実線）とを示す図である。図７に示すように、フィッティング後のモデル関数がＸ線画像データから算出されたデータ欠損率によく合致していることがわかる。
　図８は、ＦＰＤ１０２のナイキスト周波数に近い周期構造を持ったグリッド１１４を挿入して撮影されたＸ線画像データから算出されたデータ欠損率（破線及び□）を示す図である。この例では、グリッドラインが走査線と並行となるようにグリッド１１４を挿入した状態でＸ線画像データが撮影されている。図８に示すように、グリッドラインが走査線と並行となるようにグリッド１１４を挿入して状態で撮影されたＸ線画像データから算出されたデータ欠損率は、図７のデータ欠損率に「うなり」を伴った振動が加わっていることがわかる。
　フィッティング後のモデル関数は、Ｘ線波形及びＦＰＤ１０２の駆動に基づいたデータ欠損率の形状を有する。従って、フィッティング後のモデル関数は、Ｘ線画像データから算出されたデータ欠損率に入り込んだグリッドの縞や上記「うなり」からの影響を受けず、Ｘ線波形及びＦＰＤ１０２の駆動に起因した真のデータ欠損率を抽出できていることがわかる。
　ステップＳ１０４において、補正処理部１１３は、フィッティング後のモデル関数を用いて、Ｘ線画像データを補正する。補正処理部１１３は、フィッティング後のモデル関数を用いて、Ｘ線画像データに関し、Ｘ線照射開始から電荷の蓄積動作開始までに生じたデータ欠損を補正する。この補正処理は、補正後の画素値をＰC（ｘ，ｙ）、補正前の画素値をＰ（ｘ，ｙ）とすると、次の式５で行われる。
Figure JPOXMLDOC01-appb-M000005
　その結果、グリッドラインがＦＰＤ１０２の走査線と平行な方向に入っている場合でもＸ線照射開始から電荷の蓄積動作開始までに生じたデータ欠損を適切に補正し、Ｘ線画像データを復元することができる。つまり、グリッドを用いて放射線撮影が行われた場合であっても、補正処理部１１３は、放射線検出の際に生じたＸ線画像データにおけるデータ欠損を適切に補正することができる。
　次に、本発明の第２の実施形態について説明する。なお、第２の実施形態に係るＸ線画像撮影装置は、図１に示した構成と同様であるため、以下の説明においても、図１の符号を用いるものとする。
　第１の実施形態では、「Ｘ線検出駆動」中に、フレーム毎に偶数ラインと奇数ラインとのリセットを切り替える駆動を行うものである。しかしながら、このような特殊な駆動をしなくても、Ｘ線照射開始からＸ線検知までに生じるデータ欠損を補正することが可能である。
　ＦＰＤ制御部１０５は、図２のステップＳ１０１と同様に、Ｘ線画像データを取得する。図９は、第２の実施形態におけるＦＰＤ１０２の駆動を示すタイミングチャートである。第２の実施形態におけるＦＰＤ１０２の駆動においては、「Ｘ線検出駆動」状態におけるリセット動作をＧ１、Ｇ２、Ｇ３、Ｇ４と順次行い、１フレームで全ての走査線をリセットする。従って、第２の実施形態において、データ欠損は、図１０に示すように連続的に生じることになる。
　データ欠損率算出部１１１は、次の式６に従って、ＦＰＤ制御部１０５で取得されたＸ線画像データからデータ欠損率ＡTを算出する。
Figure JPOXMLDOC01-appb-M000006
　式６において、ｙdetectはＸ線が検知されたラインである。即ち、ｙdetectの１つ上のライン、ｙdetect＋１からは正常に電荷の蓄積動作が開始されているため、そのラインを基準としてデータ欠損率が算出される。
　図１１は、第２の実施形態におけるデータ欠損率算出部１１１で算出されたデータ欠損率ＡT（□）と、当該データ欠損率でフィッティングしたモデル関数（実線）とを示す図である。データ欠損率ＡTは、被写体の影響を受けて変動している。フィッティング部１１２は、上記データ欠損率ＡTで式３のモデル関数をフィッティングする。式３に示すモデル関数は、Ｘ線波形及びＦＰＤ１０２の駆動に基づいたデータ欠損率の形状を有する。従って、フィッティング後のモデル関数は、Ｘ線画像データから算出されたデータ欠損率に入り込んだ被写体の影響を受けず、Ｘ線波形及びＦＰＤ１０２の駆動に起因した真のデータ欠損率を抽出できていることがわかる。
　補正処理部１１３は、フィッティングで係数が定まった式３のモデル関数（フィッティング後のモデル係数）を用いてＸ線画像データを補正する。補正処理部１１３は、フィッティングで係数が定まった式３のモデル関数（フィッティング後のモデル係数）を用いて、Ｘ線画像データに関し、Ｘ線照射開始から電荷の蓄積動作開始までに生じたデータ欠損を補正する。この補正は、補正後の画素値をＰC（ｘ，ｙ）、補正前の画素値をＰ（ｘ，ｙ）とすると、次の式７で行われる。
Figure JPOXMLDOC01-appb-M000007
　その結果、ＦＰＤの駆動で連続的にラインがリセットされる場合でも、Ｘ線照射開始から電荷の蓄積動作開始までに生じたデータ欠損を適切に補正し、Ｘ線画像データを復元することができる。つまり、グリッドを用いて放射線撮影が行われた場合であっても、補正処理部１１３は、放射線検出の際に生じたＸ線画像データにおけるデータ欠損を適切に補正することができる。
　また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
　本発明は上記実施の形態に制限されるものではなく、本発明の精神及び範囲から離脱することなく、様々な変更及び変形が可能である。従って、本発明の範囲を公にするために以下の請求項を添付する。
　本願は、2013年08月27日提出の日本国特許出願特願2013-175939を基礎として優先権を主張するものであり、その記載内容の全てをここに援用する。
１００　Ｘ線画像撮影装置
１０１　Ｘ線管
１０２　ＦＰＤ
１０３　被検体
１０４　Ｘ線発生装置
１０５　ＦＰＤ制御部
１０６　モニタ
１０７　操作部
１０８　画像保存部
１０９　画像処理部
１１０　関数算出部
１１１　データ欠損率算出部
１１２　フィッティング部
１１３　補正処理部
１１４　グリッド

Claims (6)

  1. 　放射線を検知することで撮影された画像データを取得する取得手段と、
    前記画像データにおけるデータ欠損率を算出する算出手段と、
    　データ欠損を表すモデル関数を前記データ欠損率にフィッティングさせるフィッティング手段と、
    　前記フィッティング手段によりフィッティングされた前記モデル関数に基づいて、前記画像データにおけるデータ欠損を補正する補正手段とを有する画像処理装置。
  2. 　前記モデル関数は、所定のＸ線波形を時間積分した関数であることを特徴とする請求項１に記載の画像処理装置。
  3. 　前記算出手段は、前記画像データにおいてデータ欠損が生じた画素とデータ欠損が生じていない画素との比に基づいて、前記データ欠損率を算出することを特徴とする請求項１又は２に記載の画像処理装置。
  4. 　前記画像データは、前記放射線を電荷に変換する複数のラインの画素を備える放射線検知装置において１ライン以上の間隔を空けて又はライン間隔を空けずに読み出された電荷に基づく画像データであることを特徴とする請求項１乃至３の何れか１項に記載の画像処理装置。
  5. 　画像処理装置によって実行される画像処理方法であって、
    　放射線を検知することで撮影された画像データを取得する取得ステップと、
    　前記画像データにおけるデータ欠損率を算出する算出ステップと、
    　データ欠損を表すモデル関数を前記データ欠損率にフィッティングさせるフィッティングステップと、
    　前記フィッティングステップによりフィッティングされた前記モデル関数に基づいて、前記画像データにおけるデータ欠損を補正する補正ステップとを有する画像処理方法。
  6. 　放射線を検知することで撮影された画像データを取得する取得ステップと、
    　前記画像データにおけるデータ欠損率を算出する算出ステップと、
    　データ欠損を表すモデル関数を前記データ欠損率にフィッティングさせるフィッティングステップと、
    　前記フィッティングステップによりフィッティングされた前記モデル関数に基づいて、前記画像データにおけるデータ欠損を補正する補正ステップとをコンピュータに実行させるためのプログラム。
PCT/JP2014/070338 2013-08-27 2014-08-01 画像処理装置、画像処理方法及びプログラム WO2015029695A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US14/913,667 US10311549B2 (en) 2013-08-27 2014-08-01 Image processing apparatus, image processing method, and program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013-175939 2013-08-27
JP2013175939A JP6188488B2 (ja) 2013-08-27 2013-08-27 画像処理装置、画像処理方法及びプログラム

Publications (1)

Publication Number Publication Date
WO2015029695A1 true WO2015029695A1 (ja) 2015-03-05

Family

ID=52586269

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/070338 WO2015029695A1 (ja) 2013-08-27 2014-08-01 画像処理装置、画像処理方法及びプログラム

Country Status (3)

Country Link
US (1) US10311549B2 (ja)
JP (1) JP6188488B2 (ja)
WO (1) WO2015029695A1 (ja)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013085632A (ja) * 2011-10-17 2013-05-13 Konica Minolta Medical & Graphic Inc 放射線画像撮影システム、画像処理装置および放射線画像撮影装置
JP2013138280A (ja) * 2011-12-28 2013-07-11 Konica Minolta Inc 放射線画像撮影装置および放射線画像撮影システム

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5440647A (en) * 1993-04-22 1995-08-08 Duke University X-ray procedure for removing scattered radiation and enhancing signal-to-noise ratio (SNR)
US6314160B1 (en) * 1999-12-17 2001-11-06 General Electric Company Method and apparatus for performing fluoroscopic noise reduction
US7142705B2 (en) * 2001-05-01 2006-11-28 Canon Kabushiki Kaisha Radiation image processing apparatus, image processing system, radiation image processing method, storage medium, and program
US7054474B1 (en) * 2001-07-25 2006-05-30 3D Sharp, Inc. Image noise reduction
US6404853B1 (en) * 2001-11-02 2002-06-11 Ge Medical Systems Global Technology Company, Llc Method for identifying and correcting pixels with excess pixel lag in a solid state x-ray detector
JP4522044B2 (ja) * 2002-11-15 2010-08-11 キヤノン株式会社 放射線撮影装置
JP3903027B2 (ja) 2003-08-07 2007-04-11 キヤノン株式会社 放射線画像処理方法及び装置並びにグリッドの選別方法及び装置
DE102004027163B4 (de) * 2004-06-03 2008-04-10 Siemens Ag Verfahren zur Kompensation von Bildstörungen bei einer Strahlungsbildaufnahme
US7796792B2 (en) * 2005-06-29 2010-09-14 Agfa Healthcare, N.V. Method of identifying disturbing frequencies originating from the presence of an anti-scatter grid during acquisition of a radiation image
WO2007141689A1 (en) * 2006-06-02 2007-12-13 Philips Intellectual Property & Standards Gmbh X-ray image apparatus and device for and method of calibrating an x-ray image apparatus
JP4350738B2 (ja) * 2006-10-27 2009-10-21 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線断層撮影装置およびアーチファクトの低減方法
US8270564B2 (en) * 2008-12-01 2012-09-18 Teratech Corporation Digital integration with detector correction
US8718348B2 (en) * 2009-08-07 2014-05-06 Carestream Health, Inc. Grid suppression in imaging
JP5418599B2 (ja) * 2009-10-21 2014-02-19 株式会社島津製作所 放射線撮像装置
JP5459066B2 (ja) 2010-05-24 2014-04-02 コニカミノルタ株式会社 放射線画像撮影装置
JP5653823B2 (ja) 2010-06-30 2015-01-14 富士フイルム株式会社 放射線検出素子、及び放射線画像撮影装置
US8680473B2 (en) * 2011-10-31 2014-03-25 Bruker Axs, Inc. Multiply-sampled CMOS sensor for X-ray diffraction measurements with corrections for non-ideal sensor behavior
WO2013065515A1 (ja) * 2011-11-01 2013-05-10 富士フイルム株式会社 放射線撮影装置、放射線撮影システム及び放射線撮影方法
EP2664280A3 (en) * 2012-05-14 2013-12-04 Samsung Electronics Co., Ltd X-ray imaging apparatus and control method therefor
JP5860003B2 (ja) * 2012-07-27 2016-02-16 富士フイルム株式会社 放射線画像検出装置、並びに放射線撮影システムおよびその作動方法
KR101999266B1 (ko) * 2012-12-31 2019-07-12 삼성전자주식회사 엑스선 영상 장치 및 그 제어 방법
US10058301B2 (en) * 2014-09-24 2018-08-28 Fujifilm Corporation Image analysis device, image analysis method, and program

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013085632A (ja) * 2011-10-17 2013-05-13 Konica Minolta Medical & Graphic Inc 放射線画像撮影システム、画像処理装置および放射線画像撮影装置
JP2013138280A (ja) * 2011-12-28 2013-07-11 Konica Minolta Inc 放射線画像撮影装置および放射線画像撮影システム

Also Published As

Publication number Publication date
JP2015043834A (ja) 2015-03-12
US20160203588A1 (en) 2016-07-14
US10311549B2 (en) 2019-06-04
JP6188488B2 (ja) 2017-08-30

Similar Documents

Publication Publication Date Title
JP5010637B2 (ja) 画像補正方法および画像補正装置
JP5747913B2 (ja) 放射線画像処理装置
JP2014183475A (ja) 撮像システム
US20130068955A1 (en) Control device for radiation imaging apparatus and control method therefor
JP6491545B2 (ja) 画像処理装置、放射線撮影装置、画像処理方法、プログラム、および記憶媒体
JP2011193306A (ja) 放射線画像撮影装置および放射線画像撮影システム
JP2014183542A5 (ja) 放射線撮像システム、コンピュータ及びプログラム
JP6202855B2 (ja) 画像処理装置、画像処理装置の制御方法、及びプログラム
JP2009201586A (ja) 放射線画像撮影装置
US9036783B2 (en) Gain calibration technique for digital imaging systems
JP6917782B2 (ja) 放射線撮影装置、放射線撮影方法およびプログラム
JP6188488B2 (ja) 画像処理装置、画像処理方法及びプログラム
JP2009201552A (ja) 放射線画像撮影装置
JP6305008B2 (ja) 放射線撮影装置およびその制御方法、放射線画像処理装置および方法、並びに、プログラムおよびコンピュータ可読記憶媒体
JP2016201634A (ja) 撮像装置、撮像方法及びプログラム
JP5186425B2 (ja) 画像処理方法および放射線画像撮影装置
WO2015059886A1 (ja) 放射線撮影装置およびその制御方法、放射線画像処理装置および方法、並びに、プログラムおよびコンピュータ可読記憶媒体
JP2013255040A (ja) 撮像装置、撮像システム、撮像装置の制御方法およびプログラム
JP7383673B2 (ja) 画像処理装置、放射線撮影システム、画像処理方法及びプログラム
JP2018153447A (ja) 放射線撮影装置及びその制御方法
JP2015144629A (ja) 情報処理装置、放射線撮影システム、情報処理方法、及びプログラム
JP2004020300A (ja) 放射線診断装置
US11782171B2 (en) Radiation imaging apparatus, method for controlling same, and storage medium
JP6501513B2 (ja) 放射線画像撮影装置及びその制御方法
JP2013138280A (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: 14840732

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 14913667

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14840732

Country of ref document: EP

Kind code of ref document: A1