WO2011145377A1 - 放射線画像処理装置 - Google Patents
放射線画像処理装置 Download PDFInfo
- Publication number
- WO2011145377A1 WO2011145377A1 PCT/JP2011/054280 JP2011054280W WO2011145377A1 WO 2011145377 A1 WO2011145377 A1 WO 2011145377A1 JP 2011054280 W JP2011054280 W JP 2011054280W WO 2011145377 A1 WO2011145377 A1 WO 2011145377A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image data
- data
- radiation
- value
- average value
- Prior art date
Links
- 238000012545 processing Methods 0.000 title claims abstract description 153
- 238000012937 correction Methods 0.000 claims abstract description 169
- 238000007906 compression Methods 0.000 claims abstract description 53
- 230000006835 compression Effects 0.000 claims abstract description 52
- 230000003044 adaptive effect Effects 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 230000005855 radiation Effects 0.000 claims description 215
- 238000001514 detection method Methods 0.000 claims description 94
- 238000000034 method Methods 0.000 claims description 83
- 238000012935 Averaging Methods 0.000 claims description 10
- 238000003384 imaging method Methods 0.000 abstract description 88
- 238000005192 partition Methods 0.000 abstract description 2
- 238000000638 solvent extraction Methods 0.000 abstract 1
- 230000000875 corresponding effect Effects 0.000 description 27
- 239000000758 substrate Substances 0.000 description 20
- 210000000988 bone and bone Anatomy 0.000 description 16
- 210000000056 organ Anatomy 0.000 description 16
- 238000006243 chemical reaction Methods 0.000 description 7
- 230000001678 irradiating effect Effects 0.000 description 7
- 230000002159 abnormal effect Effects 0.000 description 5
- 230000007547 defect Effects 0.000 description 5
- 230000002950 deficient Effects 0.000 description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 239000002313 adhesive film Substances 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- 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
- 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/70—Denoising; Smoothing
-
- 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/67—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response
- H04N25/671—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction
- H04N25/673—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction by using reference sources
- H04N25/674—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction by using reference sources based on the scene itself, e.g. defocusing
-
- 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/67—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response
- H04N25/671—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction
- H04N25/677—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction for reducing the column or line fixed pattern noise
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4283—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by a detector unit being housed in a cassette
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4405—Constructional features of apparatus for radiation diagnosis the apparatus being movable or portable, e.g. handheld or mounted on a trolley
-
- 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/5252—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data removing objects from field of view, e.g. removing patient table from a CT 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/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/20004—Adaptive image processing
- G06T2207/20012—Locally adaptive
-
- 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/20192—Edge enhancement; Edge preservation
-
- 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/20212—Image combination
- G06T2207/20224—Image subtraction
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/409—Edge or detail enhancement; Noise or error suppression
Definitions
- the present invention relates to a radiographic image processing apparatus, and more particularly to an image processing apparatus that removes streaky irregularities from a radiographic image captured by a radiographic image capturing apparatus.
- This type of radiographic imaging device is known as an FPD (Flat Panel Detector) and has been developed as a so-called dedicated machine integrally formed with a support stand or the like (see, for example, Patent Documents 1 and 2). .
- FPD Full Panel Detector
- a portable radiographic image capturing apparatus in which a radiation detection element or the like is housed in a housing and can be carried has been developed and put into practical use (see, for example, Patent Documents 3 and 4).
- a so-called direct type radiographic imaging device that generates a charge by a detection element in accordance with a dose of irradiated radiation such as X-rays and converts it into an electrical signal, or a scintillator that radiates irradiated radiation
- a scintillator that radiates irradiated radiation
- indirect radiation image that is converted into an electric signal by generating electric charge with a photoelectric conversion element such as a photodiode according to the energy of the converted electromagnetic wave after being converted into electromagnetic wave of other wavelengths such as visible light
- the detection element in the direct type radiographic imaging apparatus and the photoelectric conversion element in the indirect type radiographic imaging apparatus are collectively referred to as a radiation detection element.
- Patent Documents 5 and 6 a radiation image processing apparatus for removing a non-uniformity in a radiographic image from a radiographic image obtained by a conventional CR apparatus and generating a non-uniform radiographic image. Etc. have been proposed.
- a nonuniformity radiographic image is generated by extracting a nonuniformity component from a captured radiographic image and removing the extracted nonuniformity component from the original radiographic image.
- radiation images are taken by irradiating a plate having a photostimulable phosphor layer with radiation through a subject, and the plates are transported to the radiation image reading device. A radiation image is read out.
- the laser light emitted from the laser diode 101 is converted into parallel light by the collimator lens 102, and a part thereof is reflected by the beam splitter 103 and sent to the detector 104.
- the result of monitoring by the detector 104 is fed back to the laser driver circuit 105, and the laser output from the laser diode 101 is adjusted.
- the imaging lens 106 formed of a cylindrical lens, reflected on the mirror surface of the polygon mirror 107, and the reflected light is irradiated to the image plane of the plate P through the f ⁇ lens 108 and the cylindrical mirror 109.
- the polygon mirror 109 rotates, the laser light incident on the image plane of the plate P moves in the main scanning direction X on the image plane as excitation light, and scans the image plane along the reading line Z.
- image data for each pixel on the plate P is acquired.
- the direction perpendicular to the main scanning direction X is the sub-scanning direction Y, and the plate P is scanned with the laser beam that is the excitation light as described above while being moved little by little in the sub-scanning direction Y.
- the image data of each pixel on the plate P is read in a two-dimensional manner, and the image data is read.
- the read image data are arranged two-dimensionally, and each strip-shaped area extending in the sub-scanning direction Y, for example, is arranged.
- the image data distribution of each pixel arranged in the sub-scanning direction Y is obtained for each strip-shaped area Sa.
- a high-pass filter process or a low-pass filter process is performed on the distribution of the image data to extract a non-uniform component from the original radiation image p, and the extracted non-uniform component is extracted from the original radio image p. It has been proposed to generate a radiation image p having no unevenness by removing it.
- Unevenness that occurs in the radiographic image p photographed by the CR device mainly causes vibration when the polygon mirror 107 of the radiographic image reading device 100 rotates or the cylindrical mirror 109 vibrates. It is generated by physical vibration.
- the distribution of the image data appears on the lower frequency side (that is, on the long cycle side). For this reason, for example, by applying a high-pass filter process to the distribution of the image data, an uneven component remains, and an image data component obtained by photographing an organ, a bone, or the like in the human body is cut off. Therefore, in this case, only the uneven component can be extracted from the radiation image p by the high-pass filter process.
- each apparatus or member constituting the radiographic imaging apparatus is adjusted so that the period of vibration thereof is generated on the high frequency side, for example. It is practically difficult to take such measures.
- a radiographic imaging apparatus normally, as shown in FIG. 5 and FIG. 7 described later, for example, a plurality of scanning lines 5 and a plurality of signal lines 6 are provided on the substrate 4 so as to intersect each other.
- a plurality of radiation detecting elements 7 are two-dimensionally arranged in each section r partitioned by 5 and the signal line 6.
- each radiation detection element 7 connected to the same line L of the scanning line 5 in each image data d read from each radiation detection element 7.
- on-voltages are applied to the lines L1 to Lx of the scanning line 5 at different timings, and the magnitude of noise generated in the power supply circuit 15a, the bias power supply 14 and the like at each timing. Is different.
- the present invention has been made in view of the above points, and provides a radiation image processing apparatus capable of accurately removing streaky unevenness superimposed on image data captured by an FPD type radiation image capturing apparatus.
- the purpose is to do.
- the radiographic image processing apparatus of the present invention includes: Each radiation detection element is provided in each section partitioned by a plurality of scanning lines and a plurality of signal lines arranged to cross each other, and the radiation detection elements are arranged in a two-dimensional manner.
- a radiographic image processing apparatus that performs image processing on each image data of each radiation detection element captured by a radiographic image capturing apparatus, When each image data is arranged two-dimensionally so as to correspond to each of the plurality of radiation detection elements arranged two-dimensionally, an area formed by arranging each image data two-dimensionally Dividing means for dividing the region into divided regions by dividing the image data into a predetermined number of the image data at least in the extending direction of the scanning line on the detection unit; A process of calculating an average value of the image data arranged in the extending direction of the same scanning line among the image data in the divided areas is performed for each scanning line, and for each divided area.
- Means for calculating average value In the profile in the extending direction of the signal line of the average value of the image data for each scanning line calculated for each segmented area, each image generated at the boundary between the photographed object and its surroundings
- Edge compression means for creating the profile obtained by compressing the difference between the average values of the data for each of the divided regions;
- Filter processing means for performing processing for applying an adaptive filter to the profile in which the difference is compressed for each of the divided regions; Based on each data on the profile to which the adaptive filter has been applied, each data at other positions in the region is restored, and correction data creating means for creating each correction data arranged two-dimensionally; With The correction data at a corresponding position is subtracted from the image data arranged in a two-dimensional manner to generate the image data from which noise has been removed for each radiation detection element. To do.
- the division means the average value calculation means, the edge compression means, the filter processing means, the correction data for the image data photographed by the FPD type radiation image photographing apparatus.
- Each process of the creation means is performed, and only the streaky unevenness component (see FIG. 29) extending in the extending direction of the scanning line peculiar to the FPD type radiographic imaging apparatus is accurately extracted from the image data.
- the streaky unevenness component see FIG. 29
- FIG. 3 is a sectional view taken along line XX in FIG. 2. It is a top view which shows the structure of the board
- FIG. 22 is a photograph showing an enlarged view of a radiation image generated based on original image data that is captured by the radiation image capturing apparatus illustrated in FIG.
- radiographic image capturing system 50 for capturing a radiographic image processed by the radiographic image processing apparatus.
- the radiographic image capturing system 50 is provided in a hospital, a doctor's office, or the like, and the case where the patient's body or the like is captured using the radiographic image capturing apparatus 1 is assumed.
- the radiographic imaging apparatus 1 includes a scintillator and the like, so-called indirect type radiographic imaging apparatus that converts irradiated radiation into electromagnetic waves of other wavelengths such as visible light and obtains an electrical signal with a photoelectric conversion element.
- a direct type radiographic imaging apparatus may be used.
- the radiographic imaging device is portable will be described, but it may be a radiographic imaging device formed integrally with a support base or the like.
- the radiation image processing apparatus 70 is provided in the front room Rb in a one-to-one correspondence with the imaging room Ra, and the radiation image processing apparatus 70 is provided in the imaging room Ra.
- a computer so-called console
- a console that is associated with the imaging room Ra and controls each device such as the imaging room Ra is provided separately, and the radiographic image processing apparatus 70 is installed in a place different from the imaging room Ra and the front room Rb.
- the radiation image processing apparatus 70 and the console may be connected via a network such as a LAN (Local Area Network).
- LAN Local Area Network
- the radiation image processing device 70 and the plurality of imaging rooms Ra, or the plurality of radiation image processing devices 70 and the plurality of imaging rooms Ra can be connected via a network.
- the imaging system 50 is appropriately configured so that the radiographic image p can be appropriately processed by the radiographic image processing apparatus 70.
- the radiographing room Ra is a room for taking radiographic images by irradiating a subject that is a part of the patient's body (that is, a radiographed region of the patient), and generating radiation of a radiation irradiating apparatus for irradiating the subject with radiation.
- a radiation source 52 of the apparatus 57 is disposed.
- a bucky device 51 into which the radiographic imaging device 1 can be loaded is provided in the radiographing room Ra.
- FIG. 2 is an external perspective view of the radiographic image capturing apparatus according to the present embodiment
- FIG. 3 is an external perspective view of the radiographic image capturing apparatus as viewed from the opposite side.
- 4 is a cross-sectional view taken along line XX in FIG.
- the radiographic imaging apparatus 1 houses a sensor panel SP including a scintillator 3, a substrate 4, and the like in a housing 2.
- a hollow rectangular tube-shaped housing body 2A having a radiation incident surface R in the housing 2 is made of a material such as a carbon plate or plastic that transmits radiation.
- the housing 2 is formed by closing the openings on both sides of the housing body 2A with lid members 2B and 2C.
- a so-called monocoque type for example, a so-called lunch box type formed of a frame plate and a back plate can be used.
- the lid member 2B on one side of the housing 2 is composed of a power switch 37, a selection switch 38, a connector 39, an LED for displaying a battery state, an operating state of the radiographic imaging apparatus 1, and the like.
- the indicator 40 and the like are arranged.
- an antenna device 41 that is a communication unit for wirelessly transferring image data or the like to the radiation image processing device 70 is embedded in the lid member 2 ⁇ / b> C on the opposite side of the housing 2.
- the antenna device 41 is provided, the arrangement location and the number of antenna devices 41 on the housing 2 are appropriately determined.
- a base 31 is disposed below the substrate 4 of the sensor panel SP, and a PCB substrate 33 on which electronic components 32 and the like are disposed on the base 31.
- a buffer member 34 and the like are attached.
- a glass substrate 35 for protecting the substrate 4 and the scintillator 3 on the radiation incident surface R side is disposed, and between the sensor panel SP and the side surface of the housing 2, A cushioning material 36 is provided to prevent them from colliding with each other.
- the scintillator 3 is affixed to a detection part P (described later) of the substrate 4.
- the scintillator 3 is composed of a phosphor as a main component and, when receiving radiation, converts it into an electromagnetic wave having a wavelength of 300 to 800 nm, that is, an electromagnetic wave centered on visible light and outputs it. .
- the substrate 4 is formed of a glass substrate. As shown in FIG. 5, a plurality of scanning lines 5 and a plurality of signal lines are provided on a surface 4 a of the substrate 4 facing the scintillator 3. 6 are arranged so as to cross each other. A radiation detection element 7 is provided in each section r partitioned by a plurality of scanning lines 5 and a plurality of signal lines 6 on the surface 4 a of the substrate 4.
- the entire section r provided with a plurality of radiation detection elements 7 arranged two-dimensionally in each section r partitioned by the scanning line 5 and the signal line 6, that is, shown by a one-dot chain line in FIG.
- the region is a detection unit P.
- a photodiode is used as the radiation detection element 7, but other than this, for example, a phototransistor or the like can also be used.
- Each radiation detection element 7 is connected to a source electrode 8s of a TFT 8 serving as a switch means, as shown in FIG.
- the drain electrode 8 d of the TFT 8 is connected to the signal line 6.
- the TFT 8 is turned on when a turn-on voltage is applied to the gate electrode 8g via the scanning line 5 from the scanning driving means 15 described later, and is accumulated in the radiation detection element 7 via the source electrode 8s and the drain electrode 8d. The charged electric charge is discharged to the signal line 6.
- the TFT 8 is turned off when an off voltage is applied to the gate electrode 8g via the connected scanning line 5, and the emission of the charge from the radiation detecting element 7 to the signal line 6 is stopped, and the radiation detecting element The electric charge is held in 7.
- the bias lines 9 are respectively connected to the plurality of radiation detection elements 7 arranged in a row, and as shown in FIG. 5, each bias line 9 is located outside the detection portion P of the substrate 4. It is bound to one connection 10 at a position.
- connection lines 10 of the scanning lines 5, the signal lines 6, and the bias lines 9 are connected to input / output terminals (also referred to as pads) 11 provided near the edge of the substrate 4.
- each input / output terminal 11 has a COF (Chip On Film) 12 in which a chip such as an IC 12 a is incorporated on a film, an anisotropic conductive adhesive film (Anisotropic Conductive Film) or an anisotropic conductive film. They are connected via an anisotropic conductive adhesive material 13 such as paste (Anisotropic Conductive Paste).
- the COF 12 is routed to the back surface 4b side of the substrate 4 and connected to the PCB substrate 33 described above on the back surface 4b side.
- substrate 4 part of sensor panel SP of the radiographic imaging apparatus 1 is formed.
- illustration of the electronic component 32 and the like is omitted.
- a bias line 9 is connected to one electrode of each radiation detection element 7, and each bias line 9 is bound to a connection 10 and connected to a bias power supply 14.
- the bias power supply 14 applies a bias voltage (reverse bias voltage in this embodiment) to the electrode of each radiation detection element 7 via the connection 10 and each bias line 9.
- each radiation detecting element 7 is connected to the source electrode 8s (denoted as S in FIG. 7) of the TFT 8, and the gate electrode 8g (denoted as G in FIG. 7) of each TFT 8. Are connected to the lines L1 to Lx of the scanning line 5 extending from the gate driver 15b of the scanning driving means 15, respectively. Further, the drain electrode 8 d (denoted as D in FIG. 7) of each TFT 8 is connected to each signal line 6.
- the scanning drive unit 15 includes a power supply circuit 15a that supplies an on voltage and an off voltage to the gate driver 15b, and a gate driver 15b that switches a voltage applied to each of the lines L1 to Lx of the scanning line 5 between the on voltage and the off voltage. It has. As described above, the gate driver 15b switches the voltage applied to the gate electrode 8g of the TFT 8 via the lines L1 to Lx of the scanning line 5 between the on-voltage and the off-voltage, It is designed to control the off state.
- the gate driver 15b of the scanning drive unit 15 scans the scanning line as shown in FIG. 28, for example. 5 are sequentially applied to the respective lines L1 to Lx, and the charges accumulated in the respective radiation detecting elements 7 by the irradiation of radiation are respectively emitted to the respective signal lines 6.
- Each signal line 6 is connected to each readout circuit 17 formed in the readout IC 16.
- the readout circuit 17 includes an amplifier circuit 18, a correlated double sampling (Correlated ⁇ ⁇ ⁇ ⁇ Sampling) circuit 19, an analog multiplexer 21, and an A / D converter 20.
- the scintillator 3 converts the radiation into electromagnetic waves of other wavelengths, and the converted electromagnetic waves are directly below the radiation detecting element 7. Is irradiated. Then, charges are generated in the radiation detection element 7 in accordance with the dose of irradiated radiation (that is, the amount of electromagnetic wave).
- each TFT 8 In the reading process of the image data D from each radiation detection element 7, when the ON voltage is applied to the predetermined line Ln of the scanning line 5 from the gate driver 15 b of the scanning driving unit 15 as described above, the scanning line 5. A turn-on voltage is applied to the gate electrode 8g of each TFT 8 connected thereto via the line Ln, and each TFT 8 is turned on. From the radiation detecting element 7 connected to each TFT 8 that is turned on, each of the TFTs 8 is turned on. Electric charges are discharged to the signal line 6 through the TFT 8.
- a voltage value is output from the amplifier circuit 18 in accordance with the amount of charge emitted from the radiation detection element 7, and is correlated double-sampled by the correlated double sampling circuit 19, and the analog value image data D is sent to the multiplexer 21. Is output.
- the image data D sequentially output from the multiplexer 21 is sequentially converted to digital image data D by the A / D converter 20, output to the storage means 23 and sequentially stored.
- the control means 22 is composed of a CPU (Central Processing Unit), ROM (Read Only Memory), RAM (Random Access Memory), a computer in which an input / output interface and the like are connected to a bus, an FPGA (Field Programmable Gate Array), and the like. ing. It may be configured by a dedicated control circuit.
- CPU Central Processing Unit
- ROM Read Only Memory
- RAM Random Access Memory
- FPGA Field Programmable Gate Array
- the control unit 22 controls the operation of each functional unit such as the scanning drive unit 15 and the readout circuit 17 of the radiographic image capturing apparatus 1.
- the control means 22 is connected to a storage means 23 constituted by a DRAM (Dynamic RAM) or the like, and a battery 24 that supplies power to each functional unit of the radiographic image capturing apparatus 1. Further, the above-described antenna device 41 is connected to the control means 22.
- the control means 22 When the radiographic imaging is performed, the control means 22 performs a readout process of reading out charges, that is, image data D from each radiation detection element 7 as described above. Further, in the present embodiment, the control means 22 thins out the image data D at a predetermined ratio based on the image data D read from each radiation detection element 7 and stored in the storage means 23, and the thinned data for preview. Dt is created.
- the control unit 22 when the radiographic image capturing is completed and the image data D is read from each radiation detection element 7 and stored in the storage unit 23, the control unit 22 quickly creates the thinning data Dt, and the thinning data Dt is described later. It transmits to the radiographic image processing apparatus 70 to perform.
- the control unit 22 transmits the thinned data Dt, the image data D, and the like via the bucky device 51.
- the control unit 22 transmits the thinned data Dt and the like via the antenna device 41. It has become.
- the thinned data Dt is, for example, image data D for one pixel every 3 ⁇ 3 pixels or 4 ⁇ 4 pixels when each image data D is arranged corresponding to each radiation detecting element 7 arranged in a two-dimensional form. Or the image data D from each radiation detection element 7 connected to each line L1, L4, L7,... Of the scanning line 5, respectively. The image data D from each radiation detection element 7 connected to each line Ln at predetermined intervals may be extracted and created.
- control unit 22 transmits the thinned data Dt to the radiation image processing apparatus 70, and then automatically sends the corresponding image data D on which the thinned data Dt is created to the radiation image processing apparatus 70. To be sent automatically.
- control means 22 applies image data D obtained by radiographic image capture every time one radiographic image capture is completed or at the time when a series of radiographic image captures are completed.
- image data D obtained by radiographic image capture every time one radiographic image capture is completed or at the time when a series of radiographic image captures are completed.
- a so-called dark reading process is automatically performed. It is also possible to perform a dark reading process before single or a series of radiographic imaging starts.
- each TFT 8 of the radiographic image capturing apparatus 1 is turned off, and left for a predetermined time in a state where the radiation image capturing apparatus 1 is not irradiated with radiation. Therefore, the dark charges and the like accumulated in them are read out as a so-called dark read value d and stored in the storage means 23.
- control unit 22 sets the read dark reading value d for each radiation detection element 7 as an offset correction value O, or performs dark reading processing a plurality of times to obtain a plurality of values obtained for each radiation detection element 7.
- the offset correction value O is calculated by, for example, averaging the dark read values d of the images, and transmitted to the radiation image processing apparatus 70.
- the bucky device 51 can be used by loading the portable radiographic imaging device 1 into a cassette holding portion (also referred to as a cassette holder) 51a. .
- the shooting room Ra is provided with the bucky devices 51A and 51B for standing-up shooting and lying-down shooting as the bucky device 51, but only one of them is provided. It may be done.
- the radiographic image capturing apparatus 1 when the radiographic image capturing apparatus 1 is used by being loaded into the bucky device 51, before the radiographic image capturing device 1 is loaded into the bucky device 51, as shown in FIG. A connector 51b provided at the tip of a cable extending from the connector is connected to the connector 39 of the radiation image capturing apparatus 1.
- the connector 51b is provided inside the cassette holding part 51a of the bucky device 51 and the radiographic imaging apparatus 1 is loaded into the cassette holding part 51a, the connector 39 and the connector 51b of the radiographic imaging apparatus 1 are automatically connected. It can also be configured to be connected.
- the imaging room Ra is provided with at least one radiation source 52 for irradiating a subject with radiation.
- one radiation source 52A among the radiation sources 52 is, for example,
- the radiographic imaging device 1 is mounted on the standing-up imaging bucky device 51A or the standing-up imaging bucky device 51B by changing the irradiation direction or the like by being suspended from the ceiling of the imaging room Ra. Can be irradiated with radiation.
- a portable radiation generator 52B is also provided.
- the portable radiation generator 52B can be carried anywhere in the imaging room Ra and can irradiate radiation in an arbitrary direction. ing. Then, the radiographic imaging device 1 in a single state that is not loaded in the bucky device 51 is applied to the body part of the patient as the subject, or the bucky device 51B for lying down photography or a bed (not shown) and the patient's body. It is possible to irradiate radiation from the portable radiation generating device 52B while being inserted between them.
- a base station 54 including a wireless antenna 53 that relays communication between the radiographic image capturing apparatus 1 and the radiographic image processing apparatus 70 is provided.
- the bucky device 51 and the base station 54 are connected by a cable or the like.
- the radiographic imaging device 1 is loaded into the bucky device 51 and used, it is output from the radiographic imaging device 1.
- the thinned data Dt and the like are transmitted to the radiation image processing device 70 via the Bucky device 51, the base station 54, and the like.
- a cradle 55 is connected to the base station 54.
- the cradle 55 may be installed in either the imaging room Ra or the front room Rb.
- the operation of the radiation generator provided with an exposure switch 56 for instructing the radiation source 52 to start radiation irradiation or the like in the front chamber (also referred to as an operation room) Rb.
- a table 57 is provided.
- An operator such as a radiologist irradiates the radiation imaging apparatus 1 with radiation from the radiation source 52 by operating the exposure switch 56.
- the radiation image processing apparatus 70 is provided in the front chamber Rb, and is configured by a computer or the like in which a CPU, a ROM, a RAM, an input / output interface, etc. (not shown) are connected to the bus. Has been. A predetermined program is stored in the ROM, and the radiation image processing apparatus 70 reads out a necessary program, expands it in a work area of the RAM, and executes various processes according to the program.
- the radiation image processing apparatus 70 is provided with a display unit 70a composed of an LCD (Liquid Crystal Display) or the like, and is connected to input means (not shown) such as a keyboard and a mouse.
- the radiographic image processing apparatus 70 is connected to a storage means 70b constituted by a hard disk or the like.
- image data output from another computer or the radiographic image processing apparatus 70 is included.
- An external device such as an imager that records and outputs a radiographic image on an image recording medium such as a film is connected via a LAN or the like.
- the radiographic image processing apparatus 70 displays a preview image on the display unit 70b based on the thinning data Dt. It has become.
- the display of the preview image is performed so that an operator such as a radiologist can check whether the radiographic imaging has been performed normally by looking at the preview image and determine whether or not re-imaging is necessary. If the operator determines that re-imaging is necessary, the radiographic imaging is performed again.
- the radiation image processing apparatus 70 Image processing is performed on the image data D obtained by image capturing.
- the radiation image processing apparatus 70 performs image processing on the image data D according to the flowchart shown in FIG.
- the radiation image processing apparatus 70 is functionally provided with at least a dividing means 71, an average value calculating means 72, an edge compressing means 73, a filter processing means 74, and a correction data creating means 75. Has been.
- step S3 in FIG. 9 the above-mentioned streaks are obtained from the image data D (that is, original image data Dori described later) pre-processed by these means 71 to 75 as described later.
- Each image from which streak-like unevenness has been removed is created by creating correction data Ds consisting only of the components of the shape unevenness (see FIG. 29) and subtracting the corresponding correction data Ds from each original image data Dori. Data D * is generated.
- the processing characteristic of the present invention is the unevenness removal processing (step S3) as described above, and the other processing is a known processing performed by normal image processing on the image data D obtained by radiographic imaging. Since this is a process, it will be briefly described or omitted.
- the radiation image processing apparatus 70 receives the image data D (that is, so-called raw data) and the offset correction value O for each radiation detection element 7 from the radiation image capturing apparatus 1, each of the data is transmitted. Is first stored in the storage means 70b, and first, gain / offset correction processing (step S1) is performed.
- the radiographic image processing apparatus 70 reads the gain correction value G for each radiation detection element 7 of the radiographic image capturing apparatus 1 stored in the storage unit 70b. 1) Based on the image data D, the original image data Dori is calculated for each radiation detection element 7 according to the equation (1).
- Dori G ⁇ (DO) (1)
- the original image data Dori means image data that is a basis for creating the correction data Ds in the following unevenness removal process (step S3), and the created correction data is subtracted to produce streaky unevenness.
- the radiation image processing apparatus 70 refers to the defective pixel map stored in the storage unit 70b and corrects the defective pixel.
- the defective pixel means the radiation detection element 7 from which the abnormal image data D is read out or the abnormal image data D read out from the radiation detection element 7.
- the radiation image processing apparatus 70 refers to the defective pixel map and discards the original image data Dori of the abnormal radiation detection element 7 indicated on the map as a defective pixel as invalid.
- the original image data Dori of the abnormal radiation detection element 7 is obtained by linear interpolation using the original image data Dori of the normal radiation detection element 7 adjacent to the abnormal radiation detection element 7 on the detection unit P.
- the radiation image processing apparatus 70 performs the unevenness removal processing (step S3) characteristic of the present invention.
- the unevenness removal process (step S3) is preferably performed after the gain / offset correction process (step S1) and the defect correction process (step S2) as described above.
- step S4 when the logarithmic conversion process or the like (step S4) described later is performed, the original image data Dori having a relatively small value among the original image data Dori calculated according to the above equation (1) is different in value.
- the original image data Dori having a relatively large value is compressed and the difference in value is compressed, so that it becomes difficult to correct streaky irregularities. Therefore, this unevenness removal process (step S3) is preferably performed before the logarithmic conversion process or the like (step S4).
- step S3 is composed of only the above-described streaky unevenness component (see FIG. 29) from the original image data Dori as described above by each means from the dividing means 71 to the correction data creating means 75.
- Correction data Ds is created. Then, by subtracting the corresponding correction data Ds from each original image data Dori, each image data D * from which streaky unevenness has been removed is generated.
- each of the original image data Dori is a plurality of radiations arranged two-dimensionally on the detection unit P (see FIG. 5 and FIG. 7).
- the detection unit P see FIG. 5 and FIG. 7.
- the original image data Dori is represented as Dori (m, n).
- the original image data Dori (m, n) is simply represented by (m, n).
- the extending direction of the scanning line 5 on the detection unit P (hereinafter referred to as the scanning line direction) is shown in the figure. It is represented by an arrow A, and the extending direction of the signal line 6 (hereinafter referred to as a signal line direction) is represented by an arrow B in the drawing.
- the dividing means 71 of the radiographic image processing apparatus 70 forms a region Q formed by arranging the original image data Dori (m, n) in a two-dimensional manner as shown in FIG. At least in the scanning line direction A, the image is divided into a predetermined number of original image data Dori, and the region Q is divided into divided regions Sk.
- the predetermined number of original image data Dori in the scanning line direction A of the segmented region Sk is set to an appropriate number, for example, 20 or the like.
- the signal line direction B of the segmented region Sk can be all the original image data Dori in the signal line direction B of the region Q, and the segmented region Sk is defined as the signal line direction B.
- a predetermined number of original image data Dori may be divided and the divided region Sk may be divided in the signal line direction B.
- the average value calculating means 72 of the radiation image processing apparatus 70 includes the original image data Dori (m, n) in each divided area Sk created by dividing the area Q by the dividing means 71. Then, the average value Da of the original image data Dori (m, n) arranged in the extending direction of the line Ln of the same scanning line 5 is calculated. This processing is performed for each of the lines L1 to Lx (see FIG. 7) of the scanning line 5.
- each average value Da in the scanning line direction A of each original image data Dori (m, n) becomes a value for each line L1 to Lx of the scanning line 5, and the line number n Variable. Therefore, hereinafter, the average value Da for the line Ln of the scanning line 5 is expressed as Da (n). Further, when those average values Da (n) are arranged in the signal line direction B, a profile of each average value Da (n) in the signal line direction B is obtained as shown in FIG. Each average value Da (n) is calculated for each segmented region Sk.
- streaky irregularities extending in the scanning line direction A appear characteristically.
- the average value Db (m) in the signal line direction B is calculated, the streaky unevenness components extending in the scanning line direction A cancel each other when the average value Db (m) is calculated. Therefore, the streaky unevenness component cannot be extracted from the original image data Dori.
- the original image data Dori that has been subjected to the gain / offset processing (step S1 in FIG. 9) and the defect correction processing (step S2) is superimposed not only with streaky unevenness components but also with various noise components.
- the streaky uneven component remains, but other various types Noise components are reduced because they cancel each other.
- the average value Da of the original image data Dori (m, n) in the scanning line direction A the average value Da becomes data including mainly streaky unevenness components as noise components.
- each average value Da (n) in the signal line direction B is obtained in this way, for example, the average value is obtained at a portion where an organ, a bone, or the like in the patient's human body is photographed and at other portions. Da (n) becomes a different value.
- FIG. 14 shows a case where the average value Da (n) of a portion where an organ, a bone, or the like is imaged is larger than the average value Da (n) of the surrounding portion.
- the edge compression unit 73 of the radiographic image processing apparatus 70 uses the average value Da (n) calculated for each segmented region Sk in the profile in the signal line direction B, so A profile obtained by compressing the difference of the average value Da (n) generated at the boundary portion between the object such as an organ or bone and its surroundings is created for each segmented region Sk.
- the edge compression unit 73 looks at the profile in the signal line direction B of the average value Da (n) calculated for each segmented region Sk, for example, in order from the smallest n.
- the absolute value of the difference Da (n) ⁇ Da (n ⁇ 1) calculated by subtracting the previous average value Da (n ⁇ 1) from a certain average value Da (n).
- -Da (n-1) is greater than or equal to a preset threshold value ⁇ Dath
- the difference is subtracted from Da (n + 2),...
- the average value after compressing the difference is expressed as an average value DA (n) in order to distinguish it from the average value Da (n).
- the edge compression means 73 calculates the difference Da (a) ⁇ Da (a ⁇ 1) from the average value Da (a) and the subsequent average values Da (a + 1), Da (a + 2),. ) Is subtracted to calculate the average value DA (a) after compression.
- the edge compression means 73 By performing the processing in this way in the edge compression means 73, for example, from the profile of the average value Da (n) of the segmented region Sk shown in FIG. 14, the difference as shown in FIG. The profile of n) is created.
- the edge compression means 73 performs this process for each segmented region Sk.
- ⁇ Dath a difference between values of adjacent average values Da (n) in an edge portion such as an organ or a bone and a value of adjacent average values Da (n) in a streaky uneven component are obtained.
- An appropriate value is set so that the difference can be accurately distinguished.
- the streaky uneven component remains in the profile of the average value DA (n) after compression without being subjected to the above-described edge compression processing, and edge portions such as organs and bones With respect to, it is possible to remove from the profile of the average value DA (n) after compression by performing the above-described edge compression processing.
- edge portions such as organs and bones are included in correction data Ds to be described later created based on the profile of the average value DA (n) after compression. It becomes possible not to remain, and it becomes possible to make the correction data Ds data only of the streaky unevenness component.
- edge compression processing in the edge compression means 73 of the radiation image processing apparatus 70 is, as described above, the difference Da calculated by subtracting the average value Da (n ⁇ 1) immediately before it from a certain average value Da (n).
- the difference Da calculated by subtracting the average value Da (n ⁇ 1) immediately before it from a certain average value Da (n).
- an average value Da (n) in a predetermined range such as n-5 to n + 5 including n of interest is calculated, and the average value Da (n) in the n and the average value in the predetermined range are calculated. It is also possible to calculate the difference with the average value of Da (n) and compress the edge as described above when the absolute value of the difference exceeds a preset threshold value.
- the edge compression processing in the edge compression means 73 is performed as described above in the signal line direction B profile of the average value Da (n) calculated for each segmented region Sk. Any method can be used as long as it can appropriately create a profile obtained by compressing the difference between the average values Da (n) generated at the boundary between an object such as an organ or bone and its surroundings, and is not limited to a specific method. .
- the filter processing means 74 of the radiation image processing apparatus 70 performs processing for applying an adaptive filter to the profile of the average value DA (n) after compression generated as described above for each divided region Sk. To do.
- a Wiener filter is used as the adaptive filter.
- the Wiener filter is an adaptive filter that minimizes a square error between true data on which noise is not superimposed and data on which noise is superimposed, and is applied to the case of this embodiment.
- the Wiener filter has a form represented by the following formula (5), for example.
- Df (n) represents the value of the data after the filtering process
- DAave represents the average value of the average value DA (n) after compression for each segmented region Sk
- ⁇ v 2 represents a preset variance
- ⁇ 2 (n) represents the variance of DA (n). That is, the variance ⁇ 2 (n) of DA (n) is the sum ⁇ (DA (n) ⁇ DAave) of (DA (n) ⁇ DAave) 2 in a minute range including n (for example, a range of n ⁇ 5 to n + 5). 2 ).
- the above equation (5) is Df (n) ⁇ DAave (7)
- the filtered data Df (n) is a value close to the average value DAave of the average value DA (n) after compression.
- the filtered data Df (n) calculated according to the above equation (5) is true data in which noise such as streaky irregularities in each segmented region Sk is not superimposed as described above.
- the profile is as shown.
- step S3 what is desired to be obtained by this unevenness removal processing (step S3) is, on the contrary, correction data Ds consisting only of the streaky unevenness component (see FIG. 29).
- the value of the variance ⁇ v 2 set in advance in the above equation (5) is set appropriately. Is required. Then, if the value of the variance ⁇ v 2 is set in advance to a value comparable to the value of the variance corresponding to the streak-like unevenness, this object can be achieved.
- the variance ⁇ of the average value DA (n) after compression is When 2 (n) is very large compared to the variance ⁇ v 2 (that is, the variance corresponding to streaky unevenness) (that is, when ⁇ 2 (n) >> ⁇ v 2 ), the above equation (6) shows As described above, Df (n) ⁇ DA (n), and when this is substituted into the above equation (8), the correction data ds (n) for each segmented region Sk becomes almost zero.
- the variance ⁇ 2 (n) of DA (n) calculated for a minute range including n, for example, n ⁇ 5 to n + 5, is compared with the variance ⁇ v 2, that is, the variance corresponding to the streak-like unevenness.
- the very large case represents, for example, a case where the average value DA (n) after compression in the minute range is generally larger or smaller than the average value DAave of the average value DA (n) after compression. ing.
- a long-period component derived from, for example, an organ or a bone in the patient's human body is more dominant than a muscular uneven component.
- the variance ⁇ 2 of the average value DA (n) after compression is given.
- (n) is substantially equal to the variance ⁇ v 2 (that is, the variance corresponding to the streaky unevenness) (that is, when ⁇ 2 (n) ⁇ v 2 )
- Df ( n) ⁇ DAave
- the correction data ds (n) for each segmented region Sk is substantially DA (n) ⁇ DAave.
- the correction data ds (n) for each segmented region Sk is not 0, but is almost equal to the streaky unevenness component included in the compressed average value DA (n).
- the value of the variance ⁇ v 2 is set to a value similar to the value of the variance corresponding to the streak-like unevenness, if the above-described minute range reflects the streak-like unevenness, the divided region Sk
- Each correction data ds (n) has a value substantially equal to the streaky unevenness component.
- the streaky unevenness is reflected in the portion of the profile in which the streaky unevenness component is reflected.
- the value of the variance ⁇ v 2 is approximately the same as the value of the variance corresponding to the streak-like unevenness. If the value is set, the streaky unevenness component can be accurately left in the correction data ds (n) for each segmented region Sk in the portion of the profile in which the streaky unevenness component is reflected. Therefore, also in this respect, it is desirable to set the value of the variance ⁇ v 2 to a value approximately equal to the value of the variance corresponding to the streaky unevenness.
- the value of the variance ⁇ v 2 has the above purpose for each radiographic imaging apparatus 1. It is set to an appropriate value that can be achieved.
- the Wiener filter As described above, by using the Wiener filter, the high-pass filter process and the low-pass filter process are not performed on the profile of the average value DA (n) after compression as in the conventional case described above. As shown in FIG. 5, it becomes possible to accurately extract the correction data ds (n) for each segmented region Sk.
- the average value DA (n) after compression subjected to the low-pass filter process is substituted for DA (n) in the above formula (5), and DAave is the average value after compression subjected to the low-pass filter process. This is the average value of DA (n).
- the correction data creation means 75 of the radiation image processing apparatus 70 corrects the correction data Ds (m, n) based on the profile of the correction data ds (n) for each segmented region Sk obtained for each segmented region Sk. Is supposed to create.
- the correction data creating unit 75 sets the correction data ds (n) for each segmented region Sk at, for example, the center position in the scanning line direction A in each segmented region Sk.
- the correction data Ds (m, n) at the position between the corresponding correction data ds (n) of the adjacent segment areas Sk arranged in the signal line direction B is created by, for example, linear interpolation, and arranged in a two-dimensional manner.
- Each corrected data Ds (m, n) is created.
- each process from the dividing means 71 to the correction data creating means 75 of the radiation image processing apparatus 70 is composed of only the streaky unevenness component (see FIG. 29) from the original image data Dori (m, n).
- Correction data Ds (m, n) is extracted.
- the radiation image processing apparatus 70 subtracts the corresponding correction data Ds (m, n) from the original image data Dori (m, n), thereby performing the original image data Dori. Striped unevenness is removed from (m, n), and each image data D * from which the striped unevenness has been removed is generated.
- the radiation image processing apparatus 70 When the radiation image processing apparatus 70 generates each image data D * from which the streak-like unevenness has been removed, subsequently, normalization processing, logarithmic conversion processing, gradation processing, and the like are performed on the generated image data D * . Then, the generated image data D * is converted into final image data (that is, diagnostic image data) (step S4 in FIG. 9).
- imaging may be performed by attaching a grid to the radiation incident surface R (see FIG. 2 and the like) of the radiographic imaging apparatus 1.
- unevenness hereinafter referred to as “grid moire” in order to distinguish from the above-mentioned streaky unevenness
- this grid moire is compared with the above streaky unevenness. Since the image appears on the image with a long period, the grid moire cannot be completely removed or removed by the unevenness removal process (step S3).
- the radiation image processing apparatus 70 performs image processing for removing known moire on the final image data generated from each of the generated image data D * to generate an image.
- a process of removing grid moire from the data is performed (step S5).
- the radiation image processing apparatus 70 stores the final image data generated as described above in the storage unit 70b (step S6). Further, if necessary, a radiographic image based on the final image data is recorded and output on an image recording medium such as a film by an imager, or transmitted to a main computer in a hospital or clinic to be stored.
- an image recording medium such as a film by an imager
- each of the dividing unit 71, the average value calculating unit 72, the edge compressing unit 73, the filter processing unit 74, and the correction data creating unit 75 Only the streak-like unevenness component is accurately extracted from the original image data Dori obtained by pre-processing the data D, and correction data Ds consisting only of the streak-like unevenness component is created.
- radiographic images are obtained based on the original image data Dori that has been subjected to normal preprocessing such as gain / offset correction processing and defect correction processing. Is generated, a streak-like unevenness (see FIG. 29) extending in the extending direction of the scanning line 5 is generated on the radiographic image, which is peculiar to the FPD type radiographic imaging apparatus.
- the average value calculation means 72 calculates the average value of the original image data Dori arranged in the extending direction of the same scanning line 5 in each divided region Sk, so that the streaky unevenness components cancel each other out.
- a streak-like uneven component can be extracted without matching.
- by applying a Wiener filter as shown in the above equation (5) in the filter processing unit 74 it is possible to accurately extract data in which streaky unevenness components appear strongly.
- the average value calculation unit 72 of the radiation image processing apparatus 70 calculates the average value Da (n) in the scanning line direction A of each original image data Dori (m, n) for each divided region Sk. By calculating, noise components other than the streaky unevenness component are canceled and reduced, and the average value Da (n) can be used as data having the streaky unevenness component as the main component of the noise component. .
- the dividing means 71 forms at least the region Q formed by arranging the original image data Dori (m, n) in a two-dimensional manner at least in the scanning line direction.
- A a predetermined number of original image data Dori is divided, and the region Q is not divided into divided regions Sk, and the original image data Dori (m, n) is scanned in the scanning line direction A over the entire region Q. It may be considered that the average value Da (n) may be calculated.
- the average value calculating means 72 uses the original image data Dori (as described above). When the average value Da (n) of the scanning line direction A of m, n) is calculated, the degree to which noise components other than the streaky unevenness component cancel each other out is reduced, and the average value Da (n) A large amount of noise components other than streaky unevenness components are contained therein.
- the predetermined number when the dividing unit 71 divides the region Q into a predetermined number of original image data Dori in the scanning line direction A is set to an appropriate number.
- the average value calculating means 72 performs the above average for each divided region Sk.
- the number of original image data Dori to be averaged when calculating the value Da (n) is the same, and the processing configuration of the arithmetic processing is simplified.
- the calculated average values Da (n) can be handled equally.
- the radiographic image becomes an image in which the boundary between a portion where an organ, a bone, or the like in the human body is imaged and the other portion is blurred, Image quality deteriorates.
- the edge compression unit 73 needs to appropriately detect the edge portion and accurately compress the edge portion.
- the edge compression unit 73 can appropriately detect the edge portion even if the edge compression is included, and each process in the edge compression unit 73 and the filter processing unit 74 is performed.
- a re-averaging unit 76 is provided between the filter processing unit 74 and the correction data generating unit 75, and each original image data Dori of each divided area Sk finely divided as described above is stored.
- the re-averaging means 76 After performing edge compression processing and filter processing on the profile of the average value Da (n) in the scanning line direction A, the re-averaging means 76 performs a predetermined number of adjacent divided regions Sk for each divided region Sk.
- the average value ds (n) ave in the scanning line direction A of the correction data ds (n) can be calculated.
- the correction data creating unit 75 Based on each average value ds (n) ave calculated by the re-averaging unit 76, the correction data creating unit 75 converts each average value ds (n) ave as described above, for example, as shown in FIG.
- the data at other positions in the region Q are restored by arranging them in the signal line direction B at the center position in the scanning line direction A of each of the predetermined number of segmented regions Sk adjacent to each other and performing linear interpolation, etc.
- the correction data Ds arranged in a dimension is created.
- the edge compression means 73 can accurately detect and compress the edge portions.
- a filtering process is performed on the profile of the average value DA (n) in the scanning line direction A of each original image data Dori of the segmented region Sk not including the edge portion or the segmented region Sk in which the edge portion is accurately compressed.
- the re-averaging means 76 calculates the average value ds (n) ave in the scanning line direction A of the correction data ds (n) for each divided region Sk for a predetermined number of adjacent divided regions Sk.
- the calculated average value ds (n) ave becomes data having a streaky uneven component as a main component of a noise component.
- the edge compression unit 73 compresses the edge part in the profile of the average value Da (n) in the scanning line direction A of each original image data Dori of each divided area Sk, the edge part is assumed to have a streak. Even if there is an uneven component of the shape, by compressing the average value Da (n), the striped uneven component of that portion is erased.
- the average value DA (n) of the edge portion is not necessarily reliable data.
- each segmented area Sk is stored in the storage unit 70b.
- each of the divided regions Sk is determined for each divided region Sk based on the profile of the correction data ds (n) for each divided region Sk obtained for each divided region Sk.
- a profile of the average value ds (n) all in the scanning line direction A of the correction data ds (n) for each segmented region Sk is created.
- the correction data ds (n) for each segmented region Sk at the edge of each segmented region Sk the data ds (n) calculated for each segmented region Sk is discarded and the average value ds for all the segmented regions Sk.
- each scanning line 5 is arranged over the entire area of the detection portion P (that is, in the case of FIG. 5, all the areas from the left side to the right side). Therefore, each original image data Dori that has been subjected to gain / offset correction processing or the like is two-dimensionally corresponding to the plurality of radiation detection elements 7 arranged two-dimensionally on the detector P. Arranged and each processing was performed by each means after the dividing means 71.
- each scanning line 5 is divided into the scanning line 5a and the scanning line 5b in the scanning line direction A on the detection unit P, respectively.
- radiographic imaging device 1 * configured as described above.
- the radiation image p1 is generated based on the original image data Dori created by pre-processing the image data D captured by the radiation image capturing apparatus 1 * , the image is enlarged in FIG. As shown, streaky irregularities extending in the scanning line direction A appear in different states in the regions corresponding to the regions Pa and Pb of the detection unit P, respectively.
- streaky irregularities extending in the signal line direction B in the radiation image are derived from the read characteristics of the read circuits 17 connected to the signal lines 6. This is different from the streaky unevenness extending in the scanning line direction A in the present invention.
- the striped unevenness extending in the signal line direction B is appropriately removed separately.
- the original image data Dori is two-dimensionally corresponding to the plurality of radiation detection elements 7 arranged two-dimensionally in the entire area of the detection unit P.
- each processing is performed by each means after the dividing means 71, as shown in an enlarged view in FIG. 23, each area is displayed on the radiation image p2, particularly in the vicinity of the boundary Bo between the areas Pa and Pb.
- the image data D photographed by the radiation image photographing apparatus 1 * configured such that at least each scanning line 5 is divided in the scanning line direction A on the detection unit P in this way is divided.
- the image processing described in the first embodiment may be performed on each of the image data D obtained in each of the areas Pa and Pb of the detection unit P where the scanning lines 5a and 5b respectively exist. desirable.
- each correction data Ds will be extracted and produced separately from each original image data Dori based on each image data D obtained by each area
- the correction data Ds extracted and created from each original image data Dori based on each image data D obtained in the region Pa of the detection unit P is subtracted and obtained in the region Pb of the detection unit P. Since the correction data Ds extracted and created from each original image data Dori based on each image data D is subtracted, each of the areas Pa and Pb of the detection unit P in which streaky unevenness is removed Image data D * can be generated accurately.
- pre-processing such as gain / offset correction processing is performed on each image data D read out by the radiation image capturing apparatus 1 to create each original image data Dori, and each original image
- Each process from the dividing means 71 to the correction data creation process 75 is performed on the data Dori to create the correction data Ds, and the corresponding correction data Ds is subtracted from each original image data Dori, thereby causing streak unevenness.
- a case has been described in which each image data D * from which is removed is generated.
- the radiation image capturing apparatus 1 is irradiated with radiation as described above. Except for not irradiating, the same processing as in radiographic imaging is performed and read out. Therefore, the stripe-shaped unevenness extending in the scanning line direction A similar to the image data D is also superimposed on the offset correction value O. ing.
- the dark reading process cannot be performed simultaneously with the radiographic image capturing and is performed before or after the radiographic image capturing, streaky unevenness superimposed on each offset correction value O read by the dark reading process, The streaky unevenness superimposed on each image data D read out by the reading process in radiographic imaging is in a completely different state.
- each original image data Dori is generated by multiplying the value DO, on which the streaky unevenness is superimposed, by the gain correction value G, and image processing is performed on each original image data Dori. It is also possible to do. This corresponds to the case described in the first embodiment or the second embodiment.
- the first and second offset correction values O The same processing as that for each original image data Dori shown in the embodiment is performed to generate each offset correction value O from which streak-like unevenness is removed.
- each offset correction value O from which the stripe-like unevenness has been removed is subtracted from the image data D on which the stripe-like unevenness is superimposed. Is possible.
- the offset correction value O is acquired by the dark reading process performed without irradiating the radiation image capturing apparatus 1 with radiation. Therefore, unlike the case of each original image data Dori, the offset correction value O does not include information on organs, bones, and the like in the human body, and there is no edge portion as described above.
- each offset correction value O corresponds to each of the plurality of radiation detection elements 7 arranged in a two-dimensional manner on the detection unit P (see FIGS. 5 and 7).
- the region Q in the case of the two-dimensional arrangement is not divided by the dividing unit 71, and each process in each unit from the average value calculating unit 72 to the correction data creating unit 75 is performed on the entire region Q. It is also possible to configure.
- the radiation image capturing apparatus 1 performs a dark reading process a plurality of times in advance to obtain a plurality of offset correction values O for each radiation detection element 7, and an average value of each of the radiation detection elements 7.
- the offset correction value O is calculated and stored in the storage unit 70b of the radiation image processing apparatus 70.
- an offset correction value Oc serving as a reference.
- the dark reading process for obtaining the reference offset correction value Oc is performed a plurality of times, and the state of streaky unevenness superimposed on the offset correction value O for each dark reading process is determined. 29. That is, as shown in FIG. 29, since the line L of the scanning line 5 where the streak-like unevenness appears strongly and the line L of the scanning line 5 where it appears weakly differ for each dark reading process, each dark reading is performed. By averaging the offset correction values O acquired in the process, the streaky irregularities superimposed on the offset correction values O cancel each other. Therefore, the reference offset correction value Oc is in a state in which no streaky unevenness is superimposed.
- the offset correction value Oc serving as a reference.
- the radiation image processing device 70 subtracts the reference offset correction value Oc for each radiation detection element 7 from the actual offset correction value O acquired in each dark reading process, for each radiation detection element 7.
- Each difference ⁇ O is calculated (step S11). As described above, since the streaky unevenness is not superimposed on the reference offset correction value Oc, the streaky unevenness superimposed on the actual offset correction value O is superimposed on each difference ⁇ O. It becomes a state.
- each process in the correction data creating means 75 from the dividing means 71 described above is performed on each calculated difference ⁇ O (step S12). Also in this case, each difference ⁇ O does not include information on organs and bones in the human body, and there is no edge portion as described above. Absent.
- the correction data creating means 75 creates the correction data ⁇ Os related to each difference ⁇ O
- the corresponding correction data ⁇ Os is subtracted from the original difference ⁇ O to obtain each difference ⁇ O * from which streaky unevenness has been removed. Generate.
- the radiation image processing apparatus 70 adds the generated differences ⁇ O * from which the streaky unevenness has been removed to the reference offset correction values Oc (step S13). As described above, since the streaky unevenness is not superimposed on the reference offset correction value Oc and the streaky unevenness is removed from each difference ⁇ O * , these added values are also streaked. The unevenness is removed.
- the radiological image processing apparatus 70 stores the added value of each difference ⁇ O * from which streak-like unevenness is removed and each offset correction value Oc as a reference as an offset correction value O * from which streak-like unevenness is removed. Save in means 70b. In this way, the radiation image processing apparatus 70 creates the offset correction value O * from which the streaky unevenness is removed from the actual offset correction value O (step S14).
- the original image data Dori is created by substituting the offset correction value O * from which the streak-like unevenness created in this way is removed into O in the above equation (1). Then, each process in the correction data creating means 75 is performed again from the dividing means 71, and the created correction data Ds is subtracted from the original image data Dori, and the respective image data D from which the streaky unevenness is removed. It is needless to explain that * is generated.
- the created offset correction value O * from which the streak-like unevenness has been removed is subtracted from the image data D. It is possible to appropriately perform each process for removing streak-like unevenness on the original image data Dori, and to accurately generate each image data D * from which the streak-like unevenness has been removed.
- It may be used in the field of radiographic imaging (especially in the medical field).
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Mathematical Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
FPD型の放射線画像撮影装置で撮影された画像データに重畳されている筋状のムラを的確に除去可能な放射線画像処理装置を提供する。 放射線画像処理装置70は、放射線画像撮影装置1で撮影された画像データDが二次元状に配列された領域Qを各区分領域Skに分割する分割手段71と、各区分領域Sk内の同一の走査線5の延在方向に並ぶ画像データDの平均値Da(n)を算出する平均値算出手段72と、撮影された対象物とその周囲との境界部分に生じる平均値Da(n)の差分を圧縮するエッジ圧縮手段73と、差分が圧縮されたプロファイルに対して適応フィルタを施すフィルタ処理手段74と、適応フィルタが施されたプロファイル上の各データds(n)に基づいて補正データDsを作成する補正データ作成手段75とを備え、画像データDから対応する補正データDsを減算処理してノイズが除去された画像データD*を生成する。
Description
本発明は、放射線画像処理装置に係り、特に、放射線画像撮影装置で撮影された放射線画像から筋状のムラを除去する画像処理装置に関する。
病気診断等を目的として、X線画像に代表される放射線を用いて撮影された放射線画像が広く用いられている。こうした医療用の放射線画像は、従来からスクリーンフィルムを用いて撮影されていたが、放射線画像のデジタル化を図るために輝尽性蛍光体シートを用いたCR(Computed Radiography)装置が開発され、最近では、照射された放射線をフォトダイオード(光電変換素子)等の放射線検出素子で検出してデジタル画像データとして取得する放射線画像撮影装置が開発されている。
このタイプの放射線画像撮影装置はFPD(Flat Panel Detector)として知られており、従来は、支持台等と一体的に形成された、いわゆる専用機として開発された(例えば特許文献1、2参照)。また、近年、放射線検出素子等をハウジングに収納して持ち運び可能とした可搬型の放射線画像撮影装置が開発され、実用化されている(例えば特許文献3、4参照)。
放射線画像撮影装置としては、照射されたX線等の放射線の線量に応じて検出素子で電荷を発生させて電気信号に変換するいわゆる直接型の放射線画像撮影装置や、照射された放射線をシンチレータ等で可視光等の他の波長の電磁波に変換した後、変換され照射された電磁波のエネルギに応じてフォトダイオード等の光電変換素子で電荷を発生させて電気信号に変換するいわゆる間接型の放射線画像撮影装置が知られている。なお、本発明では、直接型の放射線画像撮影装置における検出素子や、間接型の放射線画像撮影装置における光電変換素子を、あわせて放射線検出素子という。
ところで、特許文献5、6では、従来のCR装置で得られた放射線画像に対して、放射線画像中に生じているムラを除去して、ムラがない放射線画像を生成するための放射線画像処理装置等が提案されている。この放射線画像処理装置等では、撮影された放射線画像の中からムラの成分を抽出し、抽出したムラの成分を元の放射線画像から除去することにより、ムラのない放射線画像を生成する。
CR装置の場合、輝尽性蛍光体層を有するプレート等に被写体を介して放射線を照射して放射線画像撮影を行い、プレート等が放射線画像読取装置に搬送されて、放射線画像読取装置でプレートから放射線画像が読み出される。
その際、例えば図26に示すように、放射線画像読取装置100では、レーザダイオード101から出射したレーザ光をコリメートレンズ102で平行光とし、その一部をビームスプリッタ103で反射してディテクタ104に送り、ディテクタ104によるモニタの結果がレーザドライバ回路105にフィードバックされてレーザダイオード101からのレーザ出力が調整される。
また、シリンドリカルレンズからなる結像レンズ106で平行光の大半は屈折し、ポリゴンミラー107の鏡面上で反射し、反射光をfθレンズ108、シリンドリカルミラー109を経てプレートPの像面に照射する。その際、ポリゴンミラー109が回転することにより、プレートPの像面に入射するレーザ光は励起光として像面上を主走査方向Xに移動し、読取線Zに沿って像面を走査する。
プレートPの像面に励起光であるレーザ光が入射すると、プレートPの輝尽性蛍光体層のレーザ光が入射された位置に蓄積されている放射線エネルギが蛍光として放射される。そして、放射された蛍光は集光器110の導光管110aを通って集光器110で集光され、集光された蛍光がフォトマルやフォトダイオード111に導かれる。フォトダイオード111等は、集光された蛍光に感応して出力電流を出力する。
このようにして、プレートP上の各画素ごとの画像データが取得される。なお、主走査方向Xに直交する方向が副走査方向Yであり、プレートPは副走査方向Yに少しずつ移動されながら、上記のように励起光であるレーザ光の走査を受ける。このようにして、プレートP上の各画素の画像データが二次元状に読み取られて、画像データが読み取られる。
そして、例えば特許文献6に記載された放射線画像処理装置等では、図27に示すように、読み取られた画像データを二次元状に並べ、例えば副走査方向Yに延在する短冊状の各領域Saに分割し、短冊状の各領域Saごとに、副走査方向Yに並ぶ各画素の画像データの分布を得る。そして、この画像データの分布に対してハイパスフィルタ処理やローパスフィルタ処理を施す等して、もともとの放射線画像pの中からムラの成分を抽出し、抽出したムラの成分を元の放射線画像pから除去することで、ムラのない放射線画像pを生成することが提案されている。
CR装置で撮影された放射線画像pに生じるムラは、主に、放射線画像読取装置100のポリゴンミラー107が回転する際に振動したりシリンドリカルミラー109が振動するなど、放射線画像読取装置100の各部材の物理的な振動によって発生するものである。
そこで、上記の特許文献6に記載の発明を適用する場合、放射線画像読取装置100の各部材を調整して、それらの振動の周期が例えば高周波側で発生するように構成する等の措置がとられることが多い。放射線画像読取装置100をこのように構成すると、放射線画像pに生じているムラの成分の周期を、例えば高周波側に寄せることができる。
そして、人体内の臓器や骨等が撮影された画像データでは、画像データの分布はより低周波側(すなわち長周期側)に現れる。そのため、上記の画像データの分布に対して例えばハイパスフィルタ処理を施すことで、ムラの成分は残り、人体内の臓器や骨等が撮影された画像データの成分はカットオフされる。そのため、この場合は、ハイパスフィルタ処理により放射線画像pからムラの成分のみを抽出することができる。
しかしながら、放射線画像撮影を上記のFPD型の放射線画像撮影装置で行う場合、放射線画像撮影装置を構成する各装置や部材を調整して、それらの振動の周期を例えば高周波側で発生するように構成する等の措置をとることは、事実上、困難である。
また、放射線画像撮影装置では、上記の放射線画像読取装置の場合とは異なる原因で、撮影された放射線画像p上にムラが生じることが分かっている。
すなわち、放射線画像撮影装置では、通常、例えば後述する図5や図7に示すように、基板4上に複数の走査線5と複数の信号線6とが互いに交差するように設けられ、走査線5と信号線6により区画された各区画rに複数の放射線検出素子7が二次元状に配列されて構成される。
そして、被写体を介して放射線画像撮影装置に放射線が照射されて放射線画像撮影が行われた後、走査駆動手段15のゲートドライバ15bから、例えば図28に示すように、走査線5の各ラインL1~Lxにオン電圧が順次印加され、放射線の照射により各放射線検出素子7に蓄積された電荷がそれぞれ各信号線6を介して各読み出し回路17で読み出されることで、各放射線検出素子7から画像データが順次読み出される。
その際、走査駆動手段15の電源回路15a(図7参照)やバイアス電源14等ではそれぞれ各走査線5やバイアス線9に印加する電圧にノイズが生じている。そして、走査線5のあるラインLnにオン電圧が印加されて、走査線5の当該ラインLnに接続されている各放射線検出素子7から読み出される各電荷には、読み出し処理の瞬間に電源回路15aやバイアス電源14等で生じているノイズが電荷ノイズとして重畳される。
そのため、走査線5の当該ラインLnに接続されている各放射線検出素子7から読み出される各電荷には、同じ電荷ノイズが重畳される。また、走査線5の次のラインLn+1にオン電圧が印加されて、走査線5の当該ラインLn+1に接続されている各放射線検出素子7から読み出される各電荷にも、同様にして、同じ電荷ノイズが重畳される。
このように、各放射線検出素子7から読み出される各画像データdには、走査線5の同一のラインLに接続されている各放射線検出素子7については、同じ電荷ノイズが重畳される。しかし、図28に示したように、走査線5の各ラインL1~Lxにはそれぞれ異なるタイミングでオン電圧が印加され、各タイミングで、電源回路15aやバイアス電源14等で生じているノイズの大きさが異なる。
そのため、各画像データそのものではなく、各画像データに重畳されるノイズ成分だけを見ると、同じ走査線5に接続されている各放射線検出素子7では同じノイズが重畳されているが、走査線5が異なると重畳されるノイズの大きさが異なる値になり、図29に拡大して示すように、放射線画像p上に、走査線5の延在方向に延びる筋状のムラが現れる。
このように、FPD型の放射線画像撮影装置で画像データを取得すると、それらの画像データに基づいて生成された放射線画像p上に、走査線5の延在方向に延びる筋状のムラが現れるという特徴がある。なお、図29では、走査線5の延在方向が矢印Aで、それに直交する信号線6の延在方向が矢印Bで、それぞれ表されている。
そこで、FPD型の放射線画像撮影装置で取得された画像データに対する画像処理においては、このように特徴的に現れる筋状のムラを的確に除去することが望まれる。
本発明は、上記の点を鑑みてなされたものであり、FPD型の放射線画像撮影装置で撮影された画像データに重畳されている筋状のムラを的確に除去可能な放射線画像処理装置を提供することを目的とする。
前記の問題を解決するために、本発明の放射線画像処理装置は、
互いに交差するように配設された複数の走査線および複数の信号線により区画された各区画に各放射線検出素子が設けられ、前記各放射線検出素子が二次元状に配列された検出部を備える放射線画像撮影装置により撮影された前記各放射線検出素子ごとの各画像データに対して画像処理を行う放射線画像処理装置であって、
前記各画像データを前記二次元状に配列された複数の放射線検出素子にそれぞれ対応するように二次元状に配列した場合に、前記各画像データが二次元状に配列されて形成された領域を、少なくとも前記検出部上での前記走査線の延在方向において所定個数の前記画像データごとに区分して、前記領域を各区分領域に分割する分割手段と、
前記各区分領域内の前記各画像データのうち、同一の前記走査線の延在方向に並ぶ前記各画像データの平均値を算出する処理を前記各走査線ごとに行い、前記各区分領域ごとに行う平均値算出手段と、
前記各区分領域ごとに算出した前記各走査線ごとの前記各画像データの平均値の前記信号線の延在方向のプロファイルにおいて、撮影された対象物とその周囲との境界部分に生じる前記各画像データの平均値の差分を圧縮した前記プロファイルを前記各区分領域ごとに作成するエッジ圧縮手段と、
前記差分が圧縮された前記プロファイルに対して適応フィルタを施す処理を前記各区分領域ごとに行うフィルタ処理手段と、
前記適応フィルタが施された前記プロファイル上の各データに基づいて、前記領域における他の位置の各データを復元して、二次元状に配列された各補正データを作成する補正データ作成手段と、
を備え、
前記二次元状に配列された前記各画像データから、対応する位置の前記補正データを減算処理して、前記各放射線検出素子ごとにノイズが除去された前記各画像データを生成することを特徴とする。
互いに交差するように配設された複数の走査線および複数の信号線により区画された各区画に各放射線検出素子が設けられ、前記各放射線検出素子が二次元状に配列された検出部を備える放射線画像撮影装置により撮影された前記各放射線検出素子ごとの各画像データに対して画像処理を行う放射線画像処理装置であって、
前記各画像データを前記二次元状に配列された複数の放射線検出素子にそれぞれ対応するように二次元状に配列した場合に、前記各画像データが二次元状に配列されて形成された領域を、少なくとも前記検出部上での前記走査線の延在方向において所定個数の前記画像データごとに区分して、前記領域を各区分領域に分割する分割手段と、
前記各区分領域内の前記各画像データのうち、同一の前記走査線の延在方向に並ぶ前記各画像データの平均値を算出する処理を前記各走査線ごとに行い、前記各区分領域ごとに行う平均値算出手段と、
前記各区分領域ごとに算出した前記各走査線ごとの前記各画像データの平均値の前記信号線の延在方向のプロファイルにおいて、撮影された対象物とその周囲との境界部分に生じる前記各画像データの平均値の差分を圧縮した前記プロファイルを前記各区分領域ごとに作成するエッジ圧縮手段と、
前記差分が圧縮された前記プロファイルに対して適応フィルタを施す処理を前記各区分領域ごとに行うフィルタ処理手段と、
前記適応フィルタが施された前記プロファイル上の各データに基づいて、前記領域における他の位置の各データを復元して、二次元状に配列された各補正データを作成する補正データ作成手段と、
を備え、
前記二次元状に配列された前記各画像データから、対応する位置の前記補正データを減算処理して、前記各放射線検出素子ごとにノイズが除去された前記各画像データを生成することを特徴とする。
本発明のような方式の放射線画像処理装置によれば、FPD型の放射線画像撮影装置で撮影された画像データに対して、分割手段や平均値算出手段、エッジ圧縮手段、フィルタ処理手段、補正データ作成手段の各手段での各処理を行い、画像データから、FPD型の放射線画像撮影装置特有の走査線の延在方向に延びる筋状のムラ(図29参照)の成分のみを的確に抽出して、筋状のムラの成分のみからなる補正データを作成することが可能となる。
そのため、画像データからそれぞれ対応する補正データを減算処理することによって、画像データから、それに重畳されている筋状のムラを的確に除去することが可能となり、筋状のムラが除去された画像データを的確に生成することが可能となる。そして、筋状のムラが除去された画像データに基づいて放射線画像を生成することで、放射線画像上から筋状のムラが的確に除去することが可能となり、筋状のムラのない放射線画像を生成することが可能となる。
以下、本発明に係る放射線画像処理装置の実施の形態について、図面を参照して説明する。ただし、本発明は以下の図示例のものに限定されるものではない。
以下では、放射線画像処理装置について説明する前に、放射線画像処理装置が処理する放射線画像を撮影するための放射線画像撮影システムについて説明する。本実施形態では、放射線画像撮影システム50は、病院や医院等に設けられることが想定されており、放射線画像撮影装置1を用いて、患者の身体等を撮影する場合が想定されている。
なお、以下では、放射線画像撮影装置1が、シンチレータ等を備え、照射された放射線を可視光等の他の波長の電磁波に変換し光電変換素子で電気信号を得るいわゆる間接型の放射線画像撮影装置である場合について説明するが、直接型の放射線画像撮影装置であってもよい。また、以下では、放射線画像撮影装置が可搬型である場合について説明するが、支持台等と一体的に形成された放射線画像撮影装置であってもよい。
また、以下では、後述するように、放射線画像処理装置70が、撮影室Raと1対1に対応付けられて前室Rb内に設けられており、放射線画像処理装置70が撮影室Ra内等の各装置を制御するコンピュータ(いわゆるコンソール)の機能を兼用する場合について説明する。
しかし、例えば、撮影室Raに対応付けられ、撮影室Ra内等の各装置を制御するコンソールは別に設けておき、放射線画像処理装置70を撮影室Raや前室Rbとは別の場所に設置し、放射線画像処理装置70とコンソールとをLAN(Local Area Network)等のネットワークを介して接続するように構成することも可能である。
さらに、例えば、放射線画像処理装置70と複数の撮影室Raや、複数の放射線画像処理装置70と複数の撮影室Raとをネットワークを介して接続するように構成することも可能であり、放射線画像撮影システム50は、放射線画像処理装置70で放射線画像pを適切に画像処理できるように適宜構成される。
撮影室Raは、患者の身体の一部である被写体(すなわち患者の撮影部位)に放射線を照射して放射線画像撮影を行う部屋であり、被写体に放射線を照射するための放射線照射装置の放射線発生装置57の放射線源52等が配置されている。また、撮影室Raには、放射線画像撮影装置1を装填可能なブッキー装置51が設けられている。
ここで、放射線画像撮影装置1について説明する。図2は、本実施形態に係る放射線画像撮影装置の外観斜視図であり、図3は、放射線画像撮影装置を反対側から見た外観斜視図である。また、図4は、図2のX-X線に沿う断面図である。放射線画像撮影装置1は、図2~図4に示すように、筐体状のハウジング2内にシンチレータ3や基板4等で構成されるセンサパネルSPが収納されている。
図2や図3に示すように、本実施形態では、筐体2のうち、放射線入射面Rを有する中空の角筒状のハウジング本体部2Aは、放射線を透過するカーボン板やプラスチック等の材料で形成されており、ハウジング本体部2Aの両側の開口部を蓋部材2B、2Cで閉塞することで筐体2が形成されている。なお、筐体2をこのようないわゆるモノコック型として形成する代わりに、例えば、フレーム板とバック板とで形成された、いわゆる弁当箱型とすることも可能である。
図2に示すように、筐体2の一方側の蓋部材2Bには、電源スイッチ37や選択スイッチ38、コネクタ39、バッテリ状態や放射線画像撮影装置1の稼働状態等を表示するLED等で構成されたインジケータ40等が配置されている。
また、図3に示すように、筐体2の反対側の蓋部材2Cには、画像データ等を放射線画像処理装置70に無線で転送するための通信手段であるアンテナ装置41が埋め込まれている。なお、画像データ等を放射線画像処理装置70に有線方式で転送するように構成することも可能であり、その場合、例えば、前述したコネクタ39にケーブル等を接続して送受信するように構成される。また、アンテナ装置41を設ける場合には、アンテナ装置41の筐体2上の配置場所や配置する個数は適宜決められる。
筐体2の内部には、図4に示すように、センサパネルSPの基板4の下方側に基台31が配置され、基台31には、電子部品32等が配設されたPCB基板33や緩衝部材34等が取り付けられている。本実施形態では、基板4やシンチレータ3の放射線入射面R側には、それらを保護するためのガラス基板35が配設されており、センサパネルSPと筐体2の側面との間には、それらがぶつかり合うことを防止するための緩衝材36が設けられている。
シンチレータ3は、基板4の後述する検出部Pに貼り合わされるようになっている。本実施形態では、シンチレータ3は、例えば、蛍光体を主成分とし、放射線の入射を受けると300~800nmの波長の電磁波、すなわち可視光を中心とした電磁波に変換して出力するものが用いられる。
基板4は、本実施形態では、ガラス基板で構成されており、図5に示すように、基板4のシンチレータ3に対向する側の面4a上には、複数の走査線5と複数の信号線6とが互いに交差するように配設されている。基板4の面4a上の複数の走査線5と複数の信号線6により区画された各区画rには、放射線検出素子7がそれぞれ設けられている。
このように、走査線5と信号線6で区画された各区画rに二次元状に配列された複数の放射線検出素子7が設けられた区画rの全体、すなわち図5に一点鎖線で示される領域が検出部Pとされている。
本実施形態では、放射線検出素子7としてフォトダイオードが用いられているが、この他にも例えばフォトトランジスタ等を用いることも可能である。各放射線検出素子7は、後述する図7に示すように、スイッチ手段であるTFT8のソース電極8sに接続されている。また、TFT8のドレイン電極8dは信号線6に接続されている。
そして、TFT8は、後述する走査駆動手段15から走査線5を介してゲート電極8gにオン電圧が印加されるとオン状態となり、ソース電極8sやドレイン電極8dを介して放射線検出素子7内に蓄積されている電荷を信号線6に放出させるようになっている。また、TFT8は、接続された走査線5を介してゲート電極8gにオフ電圧が印加されるとオフ状態となり、放射線検出素子7から信号線6への電荷の放出を停止して、放射線検出素子7内に電荷を保持するようになっている。
本実施形態では、列状に配置された複数の放射線検出素子7にそれぞれバイアス線9が接続されており、図5に示すように、各バイアス線9は、基板4の検出部Pの外側の位置で1本の結線10に結束されている。
また、各走査線5や各信号線6、バイアス線9の結線10は、それぞれ基板4の端縁部付近に設けられた入出力端子(パッドともいう)11に接続されている。各入出力端子11には、図6に示すように、IC12a等のチップがフィルム上に組み込まれたCOF(Chip On Film)12が異方性導電接着フィルム(Anisotropic Conductive Film)や異方性導電ペースト(Anisotropic Conductive Paste)等の異方性導電性接着材料13を介して接続されている。
また、COF12は、基板4の裏面4b側に引き回され、裏面4b側で前述したPCB基板33に接続されるようになっている。このようにして、放射線画像撮影装置1のセンサパネルSPの基板4部分が形成されている。なお、図6では、電子部品32等の図示が省略されている。
ここで、図7を用いて放射線画像撮影装置1の回路構成について説明する。各放射線検出素子7の一方の電極にはそれぞれバイアス線9が接続されており、各バイアス線9は結線10に結束されてバイアス電源14に接続されている。バイアス電源14は、結線10および各バイアス線9を介して各放射線検出素子7の電極にそれぞれバイアス電圧(本実施形態では逆バイアス電圧)を印加するようになっている。
また、各放射線検出素子7の他方の電極はTFT8のソース電極8s(図7中ではSと表記されている。)に接続されており、各TFT8のゲート電極8g(図7中ではGと表記されている。)は、走査駆動手段15のゲートドライバ15bから延びる走査線5の各ラインL1~Lxにそれぞれ接続されている。また、各TFT8のドレイン電極8d(図7中ではDと表記されている。)は各信号線6にそれぞれ接続されている。
走査駆動手段15は、ゲートドライバ15bにオン電圧やオフ電圧を供給する電源回路15aと、走査線5の各ラインL1~Lxに印加する電圧をオン電圧とオフ電圧の間で切り替えるゲートドライバ15bとを備えている。ゲートドライバ15bは、前述したように、走査線5の各ラインL1~Lxを介してTFT8のゲート電極8gに印加する電圧をオン電圧とオフ電圧との間で切り替えて、各TFT8のオン状態とオフ状態とを制御するようになっている。
そして、被写体を介して放射線画像撮影装置1に放射線が照射されて放射線画像撮影が行われた後の読み出し処理では、走査駆動手段15のゲートドライバ15bは、例えば図28に示したように走査線5の各ラインL1~Lxにオン電圧を順次印加して、放射線の照射により各放射線検出素子7内に蓄積された電荷をそれぞれ各信号線6に放出させるようになっている。
各信号線6は、読み出しIC16内に形成された各読み出し回路17にそれぞれ接続されている。読み出し回路17は、増幅回路18と、相関二重サンプリング(Correlated Double Sampling)回路19と、アナログマルチプレクサ21と、A/D変換器20とで構成されている。
放射線画像撮影の際に、被写体を介して放射線画像撮影装置1に放射線が照射されると、シンチレータ3で放射線が他の波長の電磁波に変換され、変換された電磁波がその直下の放射線検出素子7に照射される。そして、照射された放射線の線量(すなわち電磁波の光量)に応じて放射線検出素子7内で電荷が発生する。
そして、各放射線検出素子7からの画像データDの読み出し処理において、上記のように走査駆動手段15のゲートドライバ15bから走査線5の所定のラインLnにオン電圧が印加されると、走査線5の当該ラインLnを介してそれに接続されている各TFT8のゲート電極8gにオン電圧が印加されて各TFT8がオン状態となり、オン状態となった各TFT8と接続されている放射線検出素子7から各TFT8を介して信号線6に電荷が放出される。
そして、放射線検出素子7から放出された電荷量に応じて増幅回路18から電圧値が出力され、それを相関二重サンプリング回路19で相関二重サンプリングしてアナログ値の画像データDがマルチプレクサ21に出力される。マルチプレクサ21から順次出力された画像データDは、A/D変換器20で順次デジタル値の画像データDに変換され、記憶手段23に出力されて順次保存されるようになっている。
制御手段22は、CPU(Central Processing Unit)やROM(Read Only Memory)、RAM(Random Access Memory)、入出力インターフェース等がバスに接続されたコンピュータや、FPGA(Field Programmable Gate Array)等により構成されている。専用の制御回路で構成されていてもよい。
制御手段22は、放射線画像撮影装置1の走査駆動手段15や読み出し回路17等の各機能部の動作等を制御するようになっている。制御手段22には、DRAM(Dynamic RAM)等で構成される記憶手段23や、放射線画像撮影装置1の各機能部に電力を供給するバッテリ24が接続されている。また、制御手段22には、前述したアンテナ装置41が接続されている。
制御手段22は、放射線画像撮影が行われると、上記のようにして、各放射線検出素子7から電荷すなわち画像データDを読み出す読み出し処理を行うようになっている。また、本実施形態では、制御手段22は、各放射線検出素子7から読み出して記憶手段23に保存した各画像データDに基づいて、各画像データDを所定の割合で間引いてプレビュー用の間引きデータDtを作成するようになっている。
本実施形態では、制御手段22は、放射線画像撮影が終了して各放射線検出素子7から画像データDが読み出して記憶手段23に保存すると、速やかに間引きデータDtを作成し、間引きデータDtを後述する放射線画像処理装置70に送信するようになっている。
なお、本実施形態では、後述するように、放射線画像撮影装置1がブッキー装置51に装填されている場合は、制御手段22はブッキー装置51を介して間引きデータDtや画像データD等を放射線画像処理装置70に送信し、放射線画像撮影装置1がブッキー装置51に装填されておらず単独で用いられている場合は、制御手段22はアンテナ装置41を介して間引きデータDt等を送信するようになっている。
間引きデータDtは、例えば、二次元状に配列された各放射線検出素子7に対応して各画像データDを配列した場合に3×3画素や4×4画素ごとに1画素分の画像データDを抽出するようにして作成してもよく、或いは、走査線5の各ラインL1、L4、L7、…にそれぞれ接続された各放射線検出素子7からの画像データDのように、走査線5の所定の間隔ごとの各ラインLnに接続された各放射線検出素子7からの画像データDを抽出して作成するように構成することも可能である。
また、本実施形態では、制御手段22は、間引きデータDtを放射線画像処理装置70に送信した後、当該間引きデータDtの作成の基となった対応する画像データDを放射線画像処理装置70に自動的に送信するようになっている。
一方、本実施形態では、制御手段22は、1回の放射線画像撮影が終了するごとに、或いは、一連の放射線画像撮影が終了した時点で、これらの放射線画像撮影で得られた画像データDに重畳されている暗電荷等に起因するオフセット分を補正するためのオフセット補正値Oを得るために、いわゆるダーク読取処理を自動的に行うようになっている。なお、単独或いは一連の放射線画像撮影が開始される前に、ダーク読取処理を行うように構成することも可能である。
ダーク読取処理では、放射線画像撮影装置1の各TFT8をオフ状態とし、放射線画像撮影装置1に放射線が照射されない状態で所定時間放置した後、上記の読み出し処理と同様にして、各放射線検出素子7から、それらに蓄積された暗電荷等がいわゆるダーク読取値dとして読み出され、記憶手段23に保存される。
そして、制御手段22は、読み出した各放射線検出素子7ごとのダーク読取値dをオフセット補正値Oとしたり、或いは、ダーク読取処理を複数回行って、各放射線検出素子7ごとに得られた複数のダーク読取値dを平均化する等して、オフセット補正値Oを算出して放射線画像処理装置70に送信するようになっている。
次に、放射線画像撮影システム50における各装置等について説明する。
図1に示すように、ブッキー装置51は、本実施形態では、カセッテ保持部(カセッテホルダともいう。)51aに可搬型の放射線画像撮影装置1を装填して用いることができるようになっている。本実施形態では、図1に示すように、撮影室Raには、ブッキー装置51として立位撮影用と臥位撮影用のブッキー装置51A、51Bが設置されているが、いずれか一方のみが設けられていてもよい。
本実施形態では、放射線画像撮影装置1は、ブッキー装置51に装填して用いられる場合には、図8に示すように、放射線画像撮影装置1をブッキー装置51に装填する前に、ブッキー装置51から延びるケーブルの先端に設けられたコネクタ51bを放射線画像撮影装置1のコネクタ39に接続するようになっている。
なお、例えば、ブッキー装置51のカセッテ保持部51aの内部にコネクタ51bを設けておき、放射線画像撮影装置1をカセッテ保持部51aに装填すると放射線画像撮影装置1のコネクタ39とコネクタ51bが自動的に接続されるように構成することも可能である。
図1に示すように、撮影室Raには、被写体に放射線を照射する放射線源52が少なくとも1つ設けられており、本実施形態では、放射線源52のうち、1つの放射線源52Aは、例えば撮影室Raの天井からつり下げられて配置されており、その照射方向等を変えることにより、立位撮影用のブッキー装置51Aや臥位撮影用のブッキー装置51Bに装填された放射線画像撮影装置1に対して放射線を照射することができるようになっている。
また、本実施形態では、ポータブルの放射線発生装置52Bも設けられており、ポータブルの放射線発生装置52Bは、撮影室Ra内のいかなる場所にも持ち運びでき、任意の方向に放射線を照射できるようになっている。そして、ブッキー装置51に装填されていない単独の状態の放射線画像撮影装置1を被写体である患者の身体の部分にあてがったり、臥位撮影用のブッキー装置51Bや図示しないベッドと患者の身体との間に差し込んだりした状態で、ポータブルの放射線発生装置52Bから放射線を照射することができるようになっている。
撮影室Ra内には、放射線画像撮影装置1と放射線画像処理装置70との通信を中継する無線アンテナ53を備えた基地局54が設けられている。なお、本実施形態では、ブッキー装置51と基地局54とはケーブル等で接続されており、放射線画像撮影装置1をブッキー装置51に装填して用いる場合には、放射線画像撮影装置1から出力された間引きデータDt等は、ブッキー装置51や基地局54等を介して放射線画像処理装置70に送信される。
また、基地局54には、クレードル55が接続されている。クレードル55は撮影室Raと前室Rbのいずれに設置されてもよく、撮影室Raに設置される場合には、放射線発生装置52から照射される放射線が到達しない位置、すなわち、例えば撮影室Raのコーナーの位置等に設置される。
図1に示すように、前室(操作室等ともいう。)Rbには、放射線源52に対して放射線の照射開始等を指示するための曝射スイッチ56等を備えた放射線発生装置の操作卓57が設けられている。そして、放射線技師等の操作者は、曝射スイッチ56を操作することにより放射線源52から放射線画像撮影装置1に放射線を照射するようになっている。
また、前述したように、本実施形態では、放射線画像処理装置70は、前室Rbに設けられており、図示しないCPUやROM、RAM、入出力インターフェース等がバスに接続されたコンピュータ等で構成されている。ROMには所定のプログラムが格納されており、放射線画像処理装置70は、必要なプログラムを読み出してRAMの作業領域に展開してプログラムに従って各種処理を実行するようになっている。
放射線画像処理装置70には、LCD(Liquid Crystal Display)等からなる表示部70aが設けられており、その他、キーボードやマウス等の図示しない入力手段等が接続されている。また、放射線画像処理装置70には、ハードディスク等で構成された記憶手段70bが接続されており、その他、図示を省略するが、他のコンピュータや、放射線画像処理装置70から出力された画像データに基づいて放射線画像をフィルム等の画像記録媒体に記録して出力するイメージャ等の外部機器が、LAN等を介して接続されている。
本実施形態では、放射線画像処理装置70は、前述したように、放射線画像撮影装置1から間引きデータDtが送信されてくると、間引きデータDtに基づいて表示部70b上にプレビュー画像を表示するようになっている。このプレビュー画像の表示は、放射線技師等の操作者がプレビュー画像を見て放射線画像撮影が正常に行われたか否かを確認して再撮影の要否等を判断するために行われるものであり、操作者は、再撮影が必要であると判断すると、当該放射線画像撮影を再度行う。
また、操作者が再撮影は不要であると判断して、放射線画像処理装置70に対して当該放射線画像撮影について画像処理を行う旨の入力が行われると、放射線画像処理装置70は、当該放射線画像撮影で得られた画像データDに対して画像処理を行うようになっている。
本実施形態では、放射線画像処理装置70は、図9に示すフローチャートに従って画像データDに対する画像処理を行うようになっている。
そして、放射線画像処理装置70は、図10に示すように、機能的に、少なくとも分割手段71、平均値算出手段72、エッジ圧縮手段73、フィルタ処理手段74および補正データ作成手段75を備えて構成されている。
そして、図9のムラ除去処理(ステップS3)では、これらの各手段71~75で、後述するように前処理が行われた画像データD(すなわち後述する元画像データDori)から、前述した筋状のムラ(図29参照)の成分のみからなる補正データDsを作成して、各元画像データDoriからそれぞれ対応する補正データDsを減算処理することによって、筋状のムラが除去された各画像データD*を生成するようになっている。
以下、放射線画像処理装置70における画像データDに対する画像処理について詳しく説明する。また、本実施形態に係る放射線画像処理装置70の作用について説明する。
なお、本発明に特徴的な処理は、上記のようにムラ除去処理(ステップS3)であり、その他の処理は、放射線画像撮影で得られた画像データDに対する通常の画像処理で行われる公知の処理であるため、簡単に説明したり、或いは説明を省略する。
放射線画像処理装置70は、上記のように、放射線画像撮影装置1から各放射線検出素子7ごとの画像データD(すなわちいわゆるrawデータ)とオフセット補正値Oが送信されてくると、それらの各データを記憶手段70bに保存して、まず、ゲイン・オフセット補正処理(ステップS1)を行う。
ゲイン・オフセット補正処理(ステップS1)では、放射線画像処理装置70は、記憶手段70bに記憶されている当該放射線画像撮影装置1の各放射線検出素子7ごとのゲイン補正値Gを読み出して、下記(1)式に従って、画像データDに基づいて元画像データDoriを各放射線検出素子7ごとに算出する。
Dori=G×(D-O) …(1)
Dori=G×(D-O) …(1)
なお、元画像データDoriとは、下記のムラ除去処理(ステップS3)で、補正データDsを作成する基となる画像データという意味であり、作成された補正データが減算処理されて筋状のムラが除去された各画像データD*を生成する基となる画像データという意味である。
次に、放射線画像処理装置70は、欠陥補正処理(ステップS2)では、記憶手段70bに記憶されている欠陥画素マップを参照して、欠陥画素の補正を行う。ここで、欠陥画素とは、異常な画像データDが読み出される放射線検出素子7、或いはその放射線検出素子7から読み出された異常な画像データDのことをいう。
欠陥補正処理(ステップS2)では、放射線画像処理装置70は、欠陥画素マップを参照して、欠陥画素としてマップ上に示されている異常な放射線検出素子7の元画像データDoriを無効として廃棄し、当該異常な放射線検出素子7に検出部P上で隣接する正常な放射線検出素子7の元画像データDoriで線形補間する等して、異常な放射線検出素子7の元画像データDoriとする。
放射線画像処理装置70は、続いて、本発明に特徴的なムラ除去処理(ステップS3)を行う。なお、このムラ除去処理(ステップS3)は、上記のようにゲイン・オフセット補正処理(ステップS1)や欠陥補正処理(ステップS2)を行った後に行うことが好ましい。
また、後述する対数変換処理等(ステップS4)が行われた後に行うと、上記(1)式に従って算出された元画像データDoriのうち、値が相対的に小さい元画像データDoriは値の差異が増幅され、値が相対的に大きい元画像データDoriは値の差異が圧縮されるため、筋状のムラを補正しづらくなる。そのため、このムラ除去処理(ステップS3)は、対数変換処理等(ステップS4)を行う前に行うことが好ましい。
ムラ除去処理(ステップS3)は、分割手段71から補正データ作成手段75までの各手段で、前述したように元画像データDoriから、前述した筋状のムラ(図29参照)の成分のみからなる補正データDsを作成する。そして、各元画像データDoriからそれぞれ対応する補正データDsを減算処理することによって、筋状のムラが除去された各画像データD*を生成する処理である。
まず、補正データDsを作成するための前提として、図11に示すように、各元画像データDoriを、検出部P(図5や図7参照)上で二次元状に配列された複数の放射線検出素子7にそれぞれ対応するように、二次元状に配列した場合を考える。
なお、検出部P上で、m本目の信号線6とn本目の走査線5(すなわち走査線5のラインLn)との交点の位置にある放射線検出素子7(m,n)に対応する各元画像データDoriをDori(m,n)と表す。また、図11中では、元画像データDori(m,n)が単に(m,n)で表されている。
また、元画像データDori(m,n)を図11に示したように配列させた場合、検出部P上での走査線5の延在方向(以下、走査線方向という。)は図中の矢印Aで表され、信号線6の延在方向(以下、信号線方向という。)は図中の矢印Bで表される。
放射線画像処理装置70の分割手段71は、図11に示したように各元画像データDori(m,n)が二次元状に配列されて形成された領域Qを、図12に示すように、少なくとも走査線方向Aにおいて所定個数の元画像データDoriごとに区分して、領域Qを各区分領域Skに分割するようになっている。
区分領域Skの走査線方向Aにおける元画像データDoriの所定個数は、適宜の個数に設定され、例えば20個等に設定される。また、図12に示したように、区分領域Skの信号線方向Bは、領域Qの信号線方向Bの全ての元画像データDoriとすることも可能であり、区分領域Skを信号線方向Bにおいても所定個数の元画像データDoriごとに区分して、区分領域Skを信号線方向Bにも分割するように構成してもよい。
放射線画像処理装置70の平均値算出手段72は、図13に示すように、分割手段71が領域Qを分割して作成した各区分領域Sk内の各元画像データDori(m,n)のうち、同一の走査線5のラインLnの延在方向に並ぶ各元画像データDori(m,n)の平均値Daを算出する。そして、この処理を走査線5の各ラインL1~Lx(図7参照)ごとに行う。
このように、各元画像データDori(m,n)の走査線方向Aの平均値Daを算出すると、各平均値Daは、走査線5の各ラインL1~Lxごとの値となり、ライン番号nの変数になる。そのため、以下、走査線5のラインLnについての平均値DaをDa(n)と表す。また、それらの平均値Da(n)を信号線方向Bに並べると、図14に示すように、各平均値Da(n)の信号線方向Bのプロファイルが得られる。なお、各平均値Da(n)の算出は、各区分領域Skごとに行われる。
図29に示したように、本実施形態のようなFPD型の放射線画像撮影装置1では、走査線方向Aに延びる筋状のムラが特徴的に現れる。そのような状況において、前述したように、各元画像データDori(m,n)の平均値を算出する際に、仮に、図15に示すように、各元画像データDori(m,n)の信号線方向Bの平均値Db(m)を算出するように構成してしまうと、平均値Db(m)の算出の際に走査線方向Aに延びる筋状のムラの成分が互いに打ち消しあってしまうため、元画像データDoriから筋状のムラの成分を抽出することができない。
それに対し、図13に示したように、各元画像データDori(m,n)の走査線方向Aの平均値Daを算出すると、平均値Da(n)の算出の際に、走査線方向Aに延びる筋状のムラの成分が互いに打ち消しあうことはなく、算出された平均値Da(n)に、筋状のムラの成分がそのまま残存するため、元画像データDoriから筋状のムラの成分を抽出することが可能となる。
なお、ゲイン・オフセット処理(図9のステップS1)や欠陥補正処理(ステップS2)を施されただけの元画像データDoriには、筋状のムラの成分だけでなく、種々のノイズ成分が重畳されているが、上記のように、各元画像データDori(m,n)の走査線方向Aの平均値Daを算出することにより、筋状のムラの成分は残存するが、その他の種々のノイズ成分は互いに打ち消しあうため低減される。
そのため、各元画像データDori(m,n)の走査線方向Aの平均値Daを算出することにより、平均値Daは、ノイズ成分として、主に筋状のムラの成分を含むデータとなる。
また、このようにして、各平均値Da(n)の信号線方向Bのプロファイルを得ると、例えば患者の人体内の臓器や骨等が撮影された部分とそれ以外の部分とで、平均値Da(n)が異なる値になる。なお、図14では、臓器や骨等が撮影されている部分の平均値Da(n)がその周囲の部分の平均値Da(n)よりも大きな値とされる場合が示されている。
次に、放射線画像処理装置70のエッジ圧縮手段73は、上記のように、各区分領域Skごとに算出した平均値Da(n)の信号線方向Bのプロファイルにおいて、撮影された患者の人体内の臓器や骨等の対象物とその周囲との境界部分に生じる上記の平均値Da(n)の差分を圧縮したプロファイルを、各区分領域Skごとに作成するようになっている。
具体的には、エッジ圧縮手段73は、例えば、上記のように、各区分領域Skごとに算出した平均値Da(n)の信号線方向Bのプロファイルを、例えばnが小さい方から順に見ていった場合に、ある平均値Da(n)からその直前の平均値Da(n-1)を減算して算出した差分Da(n)-Da(n-1)の絶対値|Da(n)-Da(n-1)|が、予め設定された閾値ΔDath以上に大きい場合には、当該平均値Da(n)を含むそれ以降の各平均値Da(n)、Da(n+1)、Da(n+2)、…から前記差分を減算して、差分を圧縮したプロファイルを作成するように構成される。なお、以下、この差分を圧縮した後の平均値を、平均値Da(n)と区別するために、平均値DA(n)と表す。
すなわち、区分領域Skの各平均値Da(n)の信号線方向Bのプロファイルが、例えば図14に示したように算出された場合、そのプロファイルを例えばnが小さい方から順に見ていった場合、n=aやn=bの所で平均値Da(n)が大きく変化している。そして、例えばn=aの所では、平均値Da(a)からその直前の平均値Da(a-1)を減算して算出した差分Da(a)-Da(a-1)の絶対値|Da(a)-Da(a-1)|が、予め設定された閾値ΔDath以上に大きくなる。
そのため、エッジ圧縮手段73は、平均値Da(a)とそれ以降の各平均値Da(a+1)、Da(a+2)、…から、それぞれ差分Da(a)-Da(a-1)を減算して、圧縮後の平均値DA(a)を算出する。
この場合、例えばn=aでの平均値Da(a)から差分Da(a)-Da(a-1)を減算すると、圧縮後の平均値DA(a)は、
DA(a)=Da(a)-(Da(a)-Da(a-1))
∴DA(n)=Da(a-1) …(2)
となり、直前のn=a-1における平均値Da(a-1)に等しくなる。
DA(a)=Da(a)-(Da(a)-Da(a-1))
∴DA(n)=Da(a-1) …(2)
となり、直前のn=a-1における平均値Da(a-1)に等しくなる。
また、平均値Da(a+1)、Da(a+2)、…についても、
DA(a+1)=Da(a+1)-(Da(a)-Da(a-1)) …(3)
DA(a+2)=Da(a+2)-(Da(a)-Da(a-1)) …(4)
等の各減算処理が行われて、圧縮後の平均値DA(a+1)、DA(a+2)、…が算出される。また、平均値Da(b)を含むそれ以降の各平均値Da(b)、Da(b+1)、Da(b+2)、…についても同様の処理が行われて、圧縮後の平均値DA(b)、DA(b+1)、DA(b+2)、…が算出される。
DA(a+1)=Da(a+1)-(Da(a)-Da(a-1)) …(3)
DA(a+2)=Da(a+2)-(Da(a)-Da(a-1)) …(4)
等の各減算処理が行われて、圧縮後の平均値DA(a+1)、DA(a+2)、…が算出される。また、平均値Da(b)を含むそれ以降の各平均値Da(b)、Da(b+1)、Da(b+2)、…についても同様の処理が行われて、圧縮後の平均値DA(b)、DA(b+1)、DA(b+2)、…が算出される。
エッジ圧縮手段73でこのように処理が行われることで、例えば図14に示した区分領域Skの平均値Da(n)のプロファイルから、図16に示すような差分を圧縮後の平均値DA(n)のプロファイルが作成される。エッジ圧縮手段73は、この処理を、各区分領域Skごとに行うようになっている。
上記のように患者の人体内の臓器や骨等を撮影した画像では、臓器や骨等が撮影された部分と、それ以外の部分が撮影された部分との境界、すなわちいわゆるエッジの部分で、隣接する平均値Da(n)同士の値が大きく変化する。一方、筋状のムラの成分でも隣接する平均値Da(n)同士で異なる値になるが、これらの臓器や骨等のエッジ部分における隣接する平均値Da(n)同士の値の差分は、筋状のムラの成分の場合の差分より格段に大きくなる。
そのため、上記の閾値ΔDathとして、臓器や骨等のエッジ部分における隣接する平均値Da(n)同士の値の差分と、筋状のムラの成分における隣接する平均値Da(n)同士の値の差分とを的確に区別できるような適切な値が設定される。
そして、このように設定することで、筋状のムラの成分については、上記のエッジ圧縮処理が行われずに圧縮後の平均値DA(n)のプロファイル中に残り、臓器や骨等のエッジ部分については、上記のエッジ圧縮処理が行われて、圧縮後の平均値DA(n)のプロファイル中から除去することが可能となる。
このようにして、上記のようなエッジ圧縮処理を行うことにより、圧縮後の平均値DA(n)のプロファイルに基づいて作成される後述する補正データDs中に、臓器や骨等のエッジ部分が残らないようにすることが可能となり、補正データDsを、筋状のムラの成分のみのデータとすることが可能となる。
なお、放射線画像処理装置70のエッジ圧縮手段73におけるエッジ圧縮処理は、上記のように、ある平均値Da(n)からその直前の平均値Da(n-1)を減算して算出した差分Da(n)-Da(n-1)に基づいて行う代わりに、例えば、平均値Da(n)のプロファイルにおける低周波成分を圧縮する、いわゆるダイナミック圧縮処理によって行うように構成することも可能である。
また、例えば、注目するnを含む例えばn-5~n+5等の所定の範囲における平均値Da(n)の平均値を算出し、当該nにおける平均値Da(n)と当該所定範囲における平均値Da(n)の平均値との差分を算出し、その差分の絶対値が予め設定された閾値を越えた場合に、上記のようにエッジを圧縮するように構成することも可能である。
このように、エッジ圧縮手段73におけるエッジ圧縮処理は、上記のように、各区分領域Skごとに算出した平均値Da(n)の信号線方向Bのプロファイルにおいて、撮影された患者の人体内の臓器や骨等の対象物とその周囲との境界部分に生じる上記の平均値Da(n)の差分を圧縮したプロファイルを適切に作成することができるものであればよく、特定の手法に限定されない。
次に、放射線画像処理装置70のフィルタ処理手段74は、上記のように作成された圧縮後の平均値DA(n)のプロファイルに対して、適応フィルタを施す処理を、各区分領域Skごとに行うようになっている。
本実施形態では、上記の適応フィルタとしてウィナーフィルタが用いられるようになっている。ウィナーフィルタは、よく知られているように、ノイズが重畳されていない真のデータとノイズが重畳されているデータとの二乗誤差を最小にする適応フィルタであり、本実施形態の場合に適用する場合には、ウィナーフィルタは、例えば下記(5)式で表されるような形になる。
ここで、Df(n)はフィルタ処理後のデータの値を表し、DAaveは各区分領域Skごとの圧縮後の平均値DA(n)の平均値を表し、σv2は予め設定された分散の値を表し、σ2(n)はDA(n)の分散を表す。すなわち、DA(n)の分散σ2(n)は、nを含む微小範囲(例えばn-5~n+5の範囲)の(DA(n)-DAave)2の総和Σ(DA(n)-DAave)2である。
そして、DA(n)の分散σ2(n)が設定された分散σv2に比べて非常に大きい場合(すなわちσ2(n)>>σv2の場合)には、上記(5)式は、
Df(n)≒DA(n) …(6)
となり、フィルタ処理後のデータDf(n)は、元の圧縮後の平均値DA(n)に近い値になる。
Df(n)≒DA(n) …(6)
となり、フィルタ処理後のデータDf(n)は、元の圧縮後の平均値DA(n)に近い値になる。
また、DA(n)の分散σ2(n)が設定された分散σv2とほぼ等しい場合(すなわちσ2(n)≒σv2の場合)には、上記(5)式は、
Df(n)≒DAave …(7)
となり、フィルタ処理後のデータDf(n)は、圧縮後の平均値DA(n)の平均値DAaveに近い値になる。
Df(n)≒DAave …(7)
となり、フィルタ処理後のデータDf(n)は、圧縮後の平均値DA(n)の平均値DAaveに近い値になる。
上記(5)式に従って算出されるフィルタ処理後のデータDf(n)は、上記のように各区分領域Skにおける筋状のムラ等のノイズが重畳されていない真のデータであり、図17に示すようなプロファイルになる。
一方、このムラ除去処理(ステップS3)で得たいのは、逆に、前述した筋状のムラ(図29参照)の成分のみからなる補正データDsである。
そこで、放射線画像処理装置70のフィルタ処理手段74は、下記(8)式に従って、図16に示した圧縮後の平均値DA(n)から、上記のようにして算出した筋状のムラ等のノイズが重畳されていない真のデータであるフィルタ処理後のデータDf(n)(図17参照)を減算することで、当該区分領域Skにおける補正データDs作成の基となるデータds(n)を算出するようになっている。
ds(n)=DA(n)-Df(n) …(8)
以下、このds(n)を区分領域Skごとの補正データという。
ds(n)=DA(n)-Df(n) …(8)
以下、このds(n)を区分領域Skごとの補正データという。
そして、この区分領域Skごとの補正データds(n)を筋状のムラの成分のみのデータとするために、上記の(5)式において予め設定する分散σv2の値を的確に設定することが必要となる。そして、分散σv2の値を、筋状のムラに対応する分散の値と同程度の値に予め設定しておくと、この目的を達成することが可能となる。
すなわち、分散σv2の値を筋状のムラに対応する分散の値と同程度の値に設定した場合、上記(5)式によれば、ある圧縮後の平均値DA(n)の分散σ2(n)が分散σv2(すなわち筋状のムラに対応する分散)に比べて非常に大きい場合(すなわちσ2(n)>>σv2の場合)には、上記(6)式に示したようにDf(n)≒DA(n)となり、これを上記(8)式に代入すると、区分領域Skごとの補正データds(n)はほぼ0になる。
つまり、前述したようにnを含む例えばn-5~n+5の微小範囲について算出されるDA(n)の分散σ2(n)が、分散σv2すなわち筋状のムラに対応する分散に比べて非常に大きい場合とは、例えばその微小範囲において圧縮後の平均値DA(n)が圧縮後の平均値DA(n)の平均値DAaveよりも全体的に大きい、或いは全体的に小さい場合を表している。これは、その微小範囲では、筋状のムラの成分よりも、例えば患者の人体内の臓器や骨等に由来する長周期の成分が支配的であることを意味する。
そして、上記のように区分領域Skごとの補正データds(n)を筋状のムラの成分のみのデータとするためには、そのような筋状のムラの成分以外の長周期の成分が区分領域Skごとの補正データds(n)中に残らないことが望まれる。そして、上記のように、分散σv2の値を筋状のムラに対応する分散の値と同程度の値に設定すればこの微小範囲では区分領域Skごとの補正データds(n)はほぼ0になり、筋状のムラの成分以外の長周期の成分が区分領域Skごとの補正データds(n)中に残らない。このように、分散σv2の値を筋状のムラに対応する分散の値と同程度の値に設定することが望ましい。
また、分散σv2の値を筋状のムラに対応する分散の値と同程度の値に設定すると、上記(5)式によれば、ある圧縮後の平均値DA(n)の分散σ2(n)が分散σv2(すなわち筋状のムラに対応する分散)とほぼ等しい場合(すなわちσ2(n)≒σv2の場合)には、上記(7)式に示したようにDf(n)≒DAaveとなり、これを上記(8)式に代入すると、区分領域Skごとの補正データds(n)はほぼDA(n)-DAaveになる。
つまり、この場合には、区分領域Skごとの補正データds(n)が0ではなく、圧縮後の平均値DA(n)に含まれる筋状のムラの成分とほぼ等しい値となる。このように、分散σv2の値を筋状のムラに対応する分散の値と同程度の値に設定すれば、上記の微小範囲が筋状のムラを反映した部分であれば、区分領域Skごとの補正データds(n)が筋状のムラの成分とほぼ等しい値となる。
上記のように区分領域Skごとの補正データds(n)を筋状のムラの成分のみのデータとするためには、筋状のムラの成分が反映されたプロファイルの部分では筋状のムラの成分が的確に区分領域Skごとの補正データds(n)中に残ることが必要になるが、上記のように、分散σv2の値を筋状のムラに対応する分散の値と同程度の値に設定すれば、筋状のムラの成分が反映されたプロファイルの部分で、筋状のムラの成分を的確に区分領域Skごとの補正データds(n)中に残すことが可能となる。従って、この点においても、分散σv2の値を筋状のムラに対応する分散の値と同程度の値に設定することが望ましいのである。
なお、筋状のムラの成分がどの程度の大きさで出現するかは、個々の放射線画像撮影装置1で異なるため、分散σv2の値は、放射線画像撮影装置1ごとに、上記の目的を達成することが可能な適宜の値に設定される。
以上のように、ウィナーフィルタを用いることで、前述した従来の場合のように、圧縮後の平均値DA(n)のプロファイルに対してハイパスフィルタ処理やローパスフィルタ処理を施さなくても、図18に示すように、区分領域Skごとの補正データds(n)を的確に抽出することが可能となる。
なお、上記の(5)式では、圧縮後の平均値DA(n)をそのまま用いる場合を示したが、例えば、圧縮後の平均値DA(n)にSINCフィルタを施す等してローパスフィルタ処理を行い、そのように処理された圧縮後の平均値DA(n)に対して、上記の(5)式を適用するように構成することも可能である。
その場合、上記(5)式中のDA(n)にはローパスフィルタ処理が施された圧縮後の平均値DA(n)が代入され、DAaveはローパスフィルタ処理が施された圧縮後の平均値DA(n)の平均値となる。
次に、放射線画像処理装置70の補正データ作成手段75は、各区分領域Skごとに得られた区分領域Skごとの補正データds(n)のプロファイルに基づいて、補正データDs(m,n)を作成するようになっている。
具体的には、補正データ作成手段75は、図19に示すように、例えば、区分領域Skごとの補正データds(n)を、各区分領域Skにおける走査線方向Aの例えば中央の位置にそれぞれ信号線方向Bに並べ、隣接する区分領域Skの対応する補正データds(n)同士の間の位置の補正データDs(m,n)を例えば線形補間等により作成して、二次元状に配列された各補正データDs(m,n)を作成する。
このようにして、放射線画像処理装置70の分割手段71から補正データ作成手段75における各処理で、元画像データDori(m,n)から、筋状のムラ(図29参照)の成分のみからなる補正データDs(m,n)が抽出される。
そして、放射線画像処理装置70は、図10に示したように、元画像データDori(m,n)から、それぞれ対応する補正データDs(m,n)を減算処理することによって、元画像データDori(m,n)から筋状のムラを除去し、筋状のムラが除去された各画像データD*を生成する。
放射線画像処理装置70は、筋状のムラが除去された各画像データD*を生成すると、続いて、生成された各画像データD*に対して正規化処理や対数変換処理、階調処理等の公知の画像処理を行って、生成された各画像データD*を最終的な画像データ(すなわち診断用画像データ)に変換する(図9のステップS4)。
なお、放射線画像撮影装置1を用いて放射線画像撮影を行う際、放射線画像撮影装置1の放射線入射面R(図2等参照)にグリッドを取り付ける等して撮影が行われる場合がある。このような場合には、グリッドにより放射線画像上にムラ(以下、上記の筋状のムラと区別するためグリッドモアレという。)が現れるが、このグリッドモアレは、上記の筋状のムラに比べて長周期で画像上に現れるため、上記のムラ除去処理(ステップS3)ではグリッドモアレを除去し切れず、或いは除去できない。
そのため、このような場合には、放射線画像処理装置70は、生成された各画像データD*から生成した最終的な画像データに対して、公知のモアレ除去のための画像処理を施して、画像データからグリッドモアレを除去する処理を行う(ステップS5)。
そして、放射線画像処理装置70は、以上のようにして生成した最終的な画像データを記憶手段70bに保存する(ステップS6)。また、必要に応じて、最終的な画像データに基づく放射線画像をイメージャでフィルム等の画像記録媒体に記録して出力したり、病院や医院のメインコンピュータに送信して格納させたりする。
以上のように、本実施形態に係る放射線画像処理装置70によれば、分割手段71や平均値算出手段72、エッジ圧縮手段73、フィルタ処理手段74、補正データ作成手段75の各手段で、画像データDに対して前処理が行われて得られた元画像データDoriから筋状のムラの成分のみを的確に抽出し、筋状のムラの成分のみからなる補正データDsを作成する。
前述したように、FPD型の放射線画像撮影装置で放射線画像撮影を行う場合、ゲイン・オフセット補正処理や欠陥補正処理等の通常の前処理を行っただけの各元画像データDoriに基づいて放射線画像を生成すると、放射線画像上に、FPD型の放射線画像撮影装置特有の、走査線5の延在方向に延びる筋状のムラ(図29参照)が生じる。
しかし、平均値算出手段72で、各区分領域Sk内の、同一の走査線5の延在方向に並ぶ各元画像データDoriの平均値を算出することで、筋状のムラの成分が互いに打ち消しあうことなく、筋状のムラの成分を抽出することができる。また、フィルタ処理手段74で、例えば上記(5)式に示したようなウィナーフィルタを適用することで、筋状のムラの成分が強く現れているデータを的確に抽出することが可能となる。
そのため、分割手段71から補正データ作成手段75の各手段で前述した各処理を施すことにより、元画像データDoriから走査線5の延在方向に延びる筋状のムラの成分のみを的確に抽出して、筋状のムラの成分のみからなる補正データDsを作成することが可能となる。
そして、各元画像データDoriからそれぞれ対応する補正データDsを減算処理することによって、各元画像データDoriから、それに重畳されている筋状のムラを的確に除去することが可能となり、筋状のムラが除去された各画像データD*を的確に生成することが可能となる。
そのため、筋状のムラが除去された各画像データD*に基づいて放射線画像を生成すれば、放射線画像上から筋状のムラが的確に除去され、筋状のムラのない放射線画像を生成することが可能となる。
なお、前述したように、放射線画像処理装置70の平均値算出手段72で、各区分領域Skごとに、各元画像データDori(m,n)の走査線方向Aの平均値Da(n)を算出することにより、筋状のムラの成分以外のノイズ成分が打ち消しあって低減され、平均値Da(n)を、筋状のムラの成分をノイズ成分の主成分とするデータとすることができる。
そのため、これを拡張すれば、図12に示したように、分割手段71で、各元画像データDori(m,n)が二次元状に配列されて形成された領域Qを、少なくとも走査線方向Aにおいて所定個数の元画像データDoriごとに区分して、領域Qを各区分領域Skに分割せずに、領域Qの全域に対して各元画像データDori(m,n)の走査線方向Aの平均値Da(n)を算出するように構成すればよいとも考えられる。
しかし、領域Qの全域に対して平均値Da(n)を算出するように構成すると、図14にa、bで示したようなエッジ部分が平均値Da(n)のプロファイル中に多数現れるようになり、エッジ圧縮手段73で、上記のようにエッジ部分を適切に圧縮することが困難になる。
また、上記の所定個数を小さい個数とすると(すなわち簡単に言えば、各区分領域Skの走査線方向Aの幅を狭くすると)、上記のように平均値算出手段72で各元画像データDori(m,n)の走査線方向Aの平均値Da(n)を算出する際に、筋状のムラの成分以外のノイズ成分が打ち消しあって低減される程度が低くなり、平均値Da(n)中に、筋状のムラの成分以外のノイズ成分が多く含まれる状態になる。
そのため、そのような平均値Da(n)に対してエッジ圧縮手段73以降の各手段で各処理を行っても、筋状のムラの成分を的確に抽出することができず、筋状のムラ(図29参照)の成分のみからなる補正データDsを的確に作成することができない。そのため、分割手段71で、領域Qを、走査線方向Aにおいて所定個数の元画像データDoriごとに区分する際の所定個数は、適宜の個数に設定される。
その際、上記の所定個数を、領域Qにおける走査線方向Aの元画像データDoriの総数を自然数で除算した個数に設定すれば、平均値算出手段72で、各区分領域Skごとに上記の平均値Da(n)を算出する際の平均化の対象となる元画像データDoriの数が同じになり、演算処理の処理構成が単純になる。また、算出された各平均値Da(n)を均等に扱うことが可能となる。
一方、上記のように、領域Qを分割する際の上記の所定個数を大きく設定すると、平均値Da(n)のプロファイル中にエッジ部分が多数現れるようになり、エッジ圧縮手段73で、上記のようにエッジ部分を適切に圧縮することが困難になるが、エッジ部分を的確に圧縮して除去できないと、補正データDs中にエッジ部分が残る。そして、元画像データDoriから補正データDsを減算処理すると、元画像データDori中のエッジ部分から補正データDs中のエッジ部分が減算されてしまうため、筋状のムラが除去された各画像データD*は、エッジ部分がぼやけたデータになってしまう。
そして、このような各画像データD*に基づいて放射線画像を生成すると、放射線画像は、例えば人体内の臓器や骨等が撮影された部分とそれ以外の部分との境界がぼやけた画像となり、画質が劣化する。このように、エッジ圧縮手段73で、エッジ部分を適切に検出して的確に圧縮することが必要となる。
そこで、上記の分割手段71では、領域Qを細かく分割し、すなわち上記の所定個数を小さく設定して、エッジ部分が含まれない区分領域Skの数を多くし、また、区分領域Skにエッジ部分が含まれてもエッジ圧縮手段73でエッジ部分を適切に検出できるようにして、エッジ圧縮手段73やフィルタ処理手段74での各処理を行う。
そして、図20に示すように、フィルタ処理手段74と補正データ作成手段75との間に再平均化手段76を設け、上記のように細かく分割された各区分領域Skの各元画像データDoriの走査線方向Aの平均値Da(n)のプロファイルに対して、エッジ圧縮処理やフィルタ処理を行った後、再平均化手段76で、隣接する所定個数の各区分領域Skについて、区分領域Skごとの補正データds(n)の走査線方向Aの平均値ds(n)aveを算出するように構成することが可能である。
補正データ作成手段75は、この再平均化手段76が算出した各平均値ds(n)aveに基づいて、例えば、図19に示したように各平均値ds(n)aveを、上記のように隣接する所定個数の各区分領域Skの走査線方向Aの中央の位置にそれぞれ信号線方向Bに並べ、線形補間する等して、領域Qにおける他の位置の各データを復元して、二次元状に配列された補正データDsを作成する。
このように構成すれば、領域Qを細かく分割して生成された各区分領域Skのうち、エッジ部分が含まれない区分領域Skの数が多くなり、また、エッジ部分が含まれる区分領域Skでは、区分領域Sk中に含まれるエッジ部分の数が少なくなり、かつ明確に現れるようになるため、エッジ圧縮手段73でエッジ部分を的確に検出して圧縮することが可能となる。
また、エッジ部分を含まない区分領域Skや、エッジ部分が的確に圧縮された区分領域Skの各元画像データDoriの走査線方向Aの平均値DA(n)のプロファイルに対してフィルタ処理を行った後、再平均化手段76で、隣接する所定個数の各区分領域Skについて、区分領域Skごとの補正データds(n)の走査線方向Aの平均値ds(n)aveを算出することで、上記のように、筋状のムラ以外のノイズ成分が互いに打ち消しあって低減される。そのため、算出された平均値ds(n)aveが、筋状のムラの成分をノイズ成分の主成分とするデータになる。
ところで、上記のように、エッジ圧縮手段73で、各区分領域Skの各元画像データDoriの走査線方向Aの平均値Da(n)のプロファイルにおけるエッジ部分を圧縮する場合、仮にエッジ部分に筋状のムラの成分があったとしても、平均値Da(n)を圧縮することでその部分の筋状のムラの成分が消されてしまう。
すなわち、圧縮前の平均値Da(n)に筋状のムラの情報が残っていたとしても、圧縮後の平均値DA(n)中からはその情報が失われてしまう。そのため、筋状をムラの成分を抽出するという観点から言えば、エッジ部分の平均値DA(n)の値は必ずしも信頼がおけるデータとは言えない。
そこで、例えば、エッジ圧縮手段73が、各区分領域Skのプロファイル中からエッジ部分を検出した場合、その位置、すなわち、例えば図14に示した区分領域Skではn=aやn=bを、各区分領域Skごとに記憶手段70bに保存する。また、フィルタ処理手段74によりフィルタ処理が行われた後で、各区分領域Skごとに得られた区分領域Skごとの補正データds(n)のプロファイルに基づいて、全ての区分領域Skについて、各区分領域Skごとの補正データds(n)の走査線方向Aの平均値ds(n)allのプロファイルを作成するように構成する。
そして、各区分領域Skのエッジ部分の区分領域Skごとの補正データds(n)については、各区分領域Skごとに算出したデータds(n)を破棄して全区分領域Skについての平均値ds(n)allで置換することにより、エッジ圧縮手段73におけるエッジ圧縮処理で失われた各区分領域Skのエッジ部分の筋状のムラの成分を復元することが可能となる。
[第2の実施の形態]
前述したように、画像データDや元画像データDori上に重畳される筋状のムラは、前述したように、走査駆動手段15の電源回路15a(図7参照)やバイアス電源14等で発生する電圧のノイズが走査線5やバイアス線9を介して各放射線検出素子7やTFT8に伝わることに起因して生じるものである。
前述したように、画像データDや元画像データDori上に重畳される筋状のムラは、前述したように、走査駆動手段15の電源回路15a(図7参照)やバイアス電源14等で発生する電圧のノイズが走査線5やバイアス線9を介して各放射線検出素子7やTFT8に伝わることに起因して生じるものである。
そして、画像データDの読み出し処理の際に、走査線5の各ラインL1~Lxにオン電圧が印加されるタイミングが異なるため、同じ走査線5に接続されている各放射線検出素子7には同じノイズが重畳されるが、走査線5が異なるとノイズの大きさが異なる値になる。そのため、図29に拡大して示したように、放射線画像pすなわち各画像データD上に、走査線方向Aに延びる筋状のムラとなって現れるものである。
そして、上記の第1の実施形態では、図5等に示したように、少なくとも各走査線5が検出部Pの全域(すなわち図5の場合は左側から右側までの全ての領域)にわたって配列されているため、ゲイン・オフセット補正処理等が施された各元画像データDoriを、検出部P上で二次元状に配列された複数の放射線検出素子7にそれぞれ対応するように、二次元状に配列して、分割手段71以降の各手段で各処理を行った。
しかし、FPD型の放射線画像撮影装置の中には、例えば図21に示すように、検出部P上で、少なくとも各走査線5が、走査線方向Aにおいてそれぞれ走査線5aと走査線5bに分割されて構成されているような放射線画像撮影装置1*もある。
このような放射線画像撮影装置1*では、分割された各走査線5a、5bがそれぞれ存在する検出部Pの各領域Pa、Pbで、走査線5a、5b等を介して各放射線検出素子7等に伝わるノイズが各領域Pa、Pbで異なる場合があり、また、各領域Pa、Pbで必ずしも同じタイミングで読み出し処理等が行われない場合もある。
このような場合、仮に、放射線画像撮影装置1*で撮影された各画像データDに前処理が施されて作成された各元画像データDoriに基づいて放射線画像p1を生成すると、図22に拡大して示すように、検出部Pの各領域Pa、Pbに対応する各領域で、走査線方向Aに延びる筋状のムラがそれぞれ異なる状態で現れる。
なお、図22や後述する図23、図24において、放射線画像中に信号線方向Bに延びる筋状のムラは、各信号線6ごとに接続されている各読み出し回路17の読み出し特性に由来するものであり、本発明における走査線方向Aに延びる筋状のムラとは異なるものである。この信号線方向Bに延びる筋状のムラについては、別途、除去する処理が適切に行われる。
そして、上記の第1の実施形態と同様に、このような元画像データDoriを、検出部Pの全域で二次元状に配列された複数の放射線検出素子7にそれぞれ対応するように、二次元状に配列して、分割手段71以降の各手段で各処理を行うと、図23に拡大して示すように、放射線画像p2上の、特に各領域Pa、Pbの境界Bo付近で、各領域Pa、Pbで筋状のムラが異なる状態で現れることの影響が排除し切れず、筋状のムラが除去し切れずに残存するという現象が生じる場合がある。
そこで、このように、検出部P上で少なくとも各走査線5が走査線方向Aにおいて分割されて構成されている放射線画像撮影装置1*により撮影された各画像データDに対しては、分割された走査線5a、5bがそれぞれ存在する検出部Pの各領域Pa、Pbでそれぞれ得られた各画像データDについてそれぞれ別々に第1の実施形態で説明した画像処理を行うように構成することが望ましい。
このように構成すれば、検出部Pの各領域Pa、Pbでそれぞれ得られた各画像データDに基づく各元画像データDoriから、それぞれ別々に補正データDsが抽出されて作成され、それらの補正データDsがそれぞれ別々に各元画像データDoriから減算処理される。
そのため、検出部Pの領域Paで得られた各画像データDに基づく各元画像データDoriからはそれから抽出されて作成された補正データDsが減算処理され、検出部Pの領域Pbで得られた各画像データDに基づく各元画像データDoriからはそれから抽出されて作成された補正データDsが減算処理されるため、検出部Pの各領域Pa、Pbで、筋状のムラが除去された各画像データD*をそれぞれ的確に生成することが可能となる。
そのため、図24に拡大して示すように、上記のように的確に生成された各画像データD*に基づいて生成された放射線画像p3では、筋状のムラが的確に除去され、各領域Pa、Pbの境界Bo付近においても、筋状のムラが残存していない状態とすることが可能となる。
[第3の実施の形態]
第1および第2の実施形態では、放射線画像撮影装置1で読み出された各画像データDについて、ゲイン・オフセット補正処理等の前処理を行って各元画像データDoriを作成し、各元画像データDoriに対して分割手段71から補正データ作成処理75までの各処理を行って、補正データDsを作成し、各元画像データDoriから対応する補正データDsを減算処理して、筋状のムラが除去された各画像データD*を生成する場合について説明した。
第1および第2の実施形態では、放射線画像撮影装置1で読み出された各画像データDについて、ゲイン・オフセット補正処理等の前処理を行って各元画像データDoriを作成し、各元画像データDoriに対して分割手段71から補正データ作成処理75までの各処理を行って、補正データDsを作成し、各元画像データDoriから対応する補正データDsを減算処理して、筋状のムラが除去された各画像データD*を生成する場合について説明した。
一方、ゲイン・オフセット補正処理(図9のステップS1)で各画像データDから減算されるオフセット補正値Oを取得するためのダーク読取処理は、前述したように、放射線画像撮影装置1に放射線を照射しないことを除いて放射線画像撮影と同様の処理が行われて読み出されるものであるため、オフセット補正値Oにも、画像データDと同様の走査線方向Aに延びる筋状のムラが重畳されている。
また、ダーク読取処理は、放射線画像撮影と同時に行うことができず、放射線画像撮影の前や後に行われるため、ダーク読取処理で読み出される各オフセット補正値Oに重畳される筋状のムラと、放射線画像撮影における読み出し処理で読み出される各画像データDに重畳される筋状のムラは、全く異なる状態となる。
しかし、ゲイン・オフセット補正処理で、上記(1)式に従って各画像データDから対応する各オフセット補正値Oを減算すると、減算して得られる値D-Oには、各画像データDに重畳されている筋状のムラから各オフセット補正値Oに重畳されている筋状のムラが減算されて生成される筋状のムラが重畳される状態になる。
従って、この筋状のムラが重畳された値D-Oにゲイン補正値Gを乗算して各元画像データDoriを作成し、これらの各元画像データDoriに対して画像処理を行うように構成することも可能である。これは、第1の実施形態や第2の実施形態で説明した場合に相当する。
しかし、筋状のムラが重畳された各画像データDから、別の筋状のムラが重畳された各オフセット補正値Oを減算すると、減算して得られる値D-Oに重畳されている筋状のムラは、その大きさが、2の平方根倍すなわち約1.4倍に大きくなる。そして、画像データDに重畳される筋状のムラがもともと比較的大きいような場合には、減算して得られる値D-Oに重畳される筋状のムラの大きさが異常に大きくなり、元画像データDoriから筋状のムラを的確に除去できなくなる場合がある。
そこで、ゲイン・オフセット補正処理(図9のステップS1)で、上記(1)式に従って各画像データDからオフセット補正値Oを減算する前に、このオフセット補正値Oについて、第1および第2の実施形態で示した各元画像データDoriに対する各処理と同様に処理を行い、筋状のムラが除去された各オフセット補正値Oを生成しておく。
そして、ゲイン・オフセット補正処理(図9のステップS1)で、この筋状のムラが除去された各オフセット補正値Oを、筋状のムラが重畳されている画像データDから減算するように構成することが可能である。
このように構成すれば、元画像データDoriに重畳されている筋状のムラが異常に大きくならず、元の画像データDに重畳されていた筋状のムラの状態のままで各処理を行うことが可能となり、各元画像データDoriから筋状のムラを的確に除去することが可能となる。
また、オフセット補正値Oは、前述したように、放射線画像撮影装置1に放射線を照射しない状態で行われるダーク読取処理で取得される。そのため、各元画像データDoriの場合と異なり、オフセット補正値Oには、人体内の臓器や骨等の情報は含まれておらず、上記のようなエッジ部分が存在しない。
そのため、オフセット補正値Oに対する処理では、各オフセット補正値Oを、検出部P(図5や図7参照)上で二次元状に配列された複数の放射線検出素子7にそれぞれ対応するように、二次元状に配列した場合の領域Qを、分割手段71で分割せず、領域Qの全体に対して平均値算出手段72から補正データ作成手段75までの各手段での各処理を行うように構成することも可能である。
一方、オフセット補正値Oから筋状のムラをより厳密に除去するために、以下のように構成することも可能である。
すなわち、事前に、放射線画像撮影装置1でダーク読取処理を複数回行って、各放射線検出素子7ごとに複数のオフセット補正値Oを取得しておき、それらの平均値として各放射線検出素子7ごとのオフセット補正値Oを算出して、放射線画像処理装置70の記憶手段70bに保存させておく。
以下、この各放射線検出素子7ごとに予め取得されたオフセット補正値Oを、各回の放射線画像撮影ごとに得られるオフセット補正値Oと区別するため、基準となるオフセット補正値Ocという。
また、上記のように、基準となるオフセット補正値Ocを取得するためのダーク読取処理は複数回行われ、各回のダーク読取処理ごとにオフセット補正値Oに重畳される筋状のムラの状態が異なるため、すなわち、図29に示したように、筋状のムラが強く出る走査線5のラインLと弱く出る走査線5のラインLが各回のダーク読取処理ごとに異なるため、各回のダーク読取処理で取得された各オフセット補正値Oを平均化することで、各オフセット補正値Oに重畳されている筋状のムラが互いに打ち消しあう。そのため、基準となるオフセット補正値Ocは、筋状のムラが重畳されていない状態となる。
放射線画像処理装置70では、このような基準となるオフセット補正値Ocを有する状態で、図25に示すフローチャートに従って、各回の放射線画像撮影の前や後に行われるダーク読取処理で取得された実際のオフセット補正値Oから筋状のムラを除去して、筋状のムラが除去されたオフセット補正値Oを作成する。
放射線画像処理装置70は、まず、各回のダーク読取処理で取得された実際のオフセット補正値Oから基準となるオフセット補正値Ocを各放射線検出素子7ごとに減算して、各放射線検出素子7ごとの各差分ΔOを算出する(ステップS11)。上記のように、基準となるオフセット補正値Ocには筋状のムラが重畳されていないため、各差分ΔOには、実際のオフセット補正値Oに重畳されていた筋状のムラが重畳された状態となる。
そして、算出した各差分ΔOに対して、前述した分割手段71から補正データ作成手段75における各処理を実施する(ステップS12)。なお、この場合も、各差分ΔOには、人体内の臓器や骨等の情報は含まれておらず、上記のようなエッジ部分が存在しないため、必ずしも分割手段71における分割処理を行う必要はない。
そして、補正データ作成手段75で、各差分ΔOに関する補正データΔOsを作成すると、元の差分ΔOからそれぞれ対応する補正データΔOsを減算処理して、筋状のムラが除去された各差分ΔO*を生成する。
次に、放射線画像処理装置70は、生成した筋状のムラが除去された各差分ΔO*を、基準となる各オフセット補正値Ocに加算する(ステップS13)。前述したように、基準となるオフセット補正値Ocには筋状のムラが重畳されておらず、各差分ΔO*からも筋状のムラが除去されているため、これらの加算値も筋状のムラが除去されたものとなる。
そのため、放射線画像処理装置70は、筋状のムラが除去された各差分ΔO*と基準となる各オフセット補正値Ocの加算値を、筋状のムラが除去されたオフセット補正値O*として記憶手段70bに保存する。放射線画像処理装置70は、このようにして、実際のオフセット補正値Oから、筋状のムラが除去されたオフセット補正値O*を作成する(ステップS14)。
なお、このようにして作成された筋状のムラが除去されたオフセット補正値O*を、上記(1)式のOに代入して元画像データDoriが作成され、この元画像データDoriに対して、再度、分割手段71から補正データ作成手段75における各処理が施されて、作成された補正データDsを元画像データDoriから減算処理して、筋状のムラが除去された各画像データD*が生成されることは、改めて説明するまでもない。
このように構成すれば、オフセット補正値Oから筋状のムラがより的確に除去されるため、作成された筋状のムラが除去されたオフセット補正値O*を画像データDから減算して得られる元画像データDoriに対して適切に筋状のムラを除去する各処理を施すことが可能となり、筋状のムラが除去された各画像データD*を的確に生成することが可能となる。
放射線画像撮影を行う分野(特に医療分野)において利用可能性がある。
1 放射線画像撮影装置
5、L1~Lx 走査線
6 信号線
7 放射線検出素子
70 放射線画像処理装置
71 分割手段
72 平均値算出手段
73 エッジ圧縮手段
74 フィルタ処理手段
75 補正データ作成手段
76 再平均化手段
A 走査線方向(走査線の延在方向)
B 信号線方向(信号線の延在方向)
D 画像データ
D* ノイズが除去された画像データ
DAave 差分が圧縮されたプロファイル上の各平均値の各区分領域ごとの平均値
Da(n) 平均値
Da(n)-Da(n-1) 差分
Dori 元画像データ(画像データ)
Ds、Ds(m,n) 補正データ
ds(n) 区分領域ごとの補正データ(適応フィルタが施されたプロファイル上の各データ)
ds(n)ave 区分領域ごとの補正データの平均値
O オフセット補正値
O* ノイズが除去されたオフセット補正値
Oc 基準となるオフセット補正値
P 検出部
Pa、Pb 検出部の各領域
Q 領域
r 区画
Sk 各区分領域
ΔDath 閾値
ΔO 差分
ΔO* ノイズが除去された差分
ΔOs 補正データ
σ2(n) 差分が圧縮されたプロファイル上の平均値の分散
σv2 予め設定された分散の値
5、L1~Lx 走査線
6 信号線
7 放射線検出素子
70 放射線画像処理装置
71 分割手段
72 平均値算出手段
73 エッジ圧縮手段
74 フィルタ処理手段
75 補正データ作成手段
76 再平均化手段
A 走査線方向(走査線の延在方向)
B 信号線方向(信号線の延在方向)
D 画像データ
D* ノイズが除去された画像データ
DAave 差分が圧縮されたプロファイル上の各平均値の各区分領域ごとの平均値
Da(n) 平均値
Da(n)-Da(n-1) 差分
Dori 元画像データ(画像データ)
Ds、Ds(m,n) 補正データ
ds(n) 区分領域ごとの補正データ(適応フィルタが施されたプロファイル上の各データ)
ds(n)ave 区分領域ごとの補正データの平均値
O オフセット補正値
O* ノイズが除去されたオフセット補正値
Oc 基準となるオフセット補正値
P 検出部
Pa、Pb 検出部の各領域
Q 領域
r 区画
Sk 各区分領域
ΔDath 閾値
ΔO 差分
ΔO* ノイズが除去された差分
ΔOs 補正データ
σ2(n) 差分が圧縮されたプロファイル上の平均値の分散
σv2 予め設定された分散の値
Claims (9)
- 互いに交差するように配設された複数の走査線および複数の信号線により区画された各区画に各放射線検出素子が設けられ、前記各放射線検出素子が二次元状に配列された検出部を備える放射線画像撮影装置により撮影された前記各放射線検出素子ごとの各画像データに対して画像処理を行う放射線画像処理装置であって、
前記各画像データを前記二次元状に配列された複数の放射線検出素子にそれぞれ対応するように二次元状に配列した場合に、前記各画像データが二次元状に配列されて形成された領域を、少なくとも前記検出部上での前記走査線の延在方向において所定個数の前記画像データごとに区分して、前記領域を各区分領域に分割する分割手段と、
前記各区分領域内の前記各画像データのうち、同一の前記走査線の延在方向に並ぶ前記各画像データの平均値を算出する処理を前記各走査線ごとに行い、前記各区分領域ごとに行う平均値算出手段と、
前記各区分領域ごとに算出した前記各走査線ごとの前記各画像データの平均値の前記信号線の延在方向のプロファイルにおいて、撮影された対象物とその周囲との境界部分に生じる前記各画像データの平均値の差分を圧縮した前記プロファイルを前記各区分領域ごとに作成するエッジ圧縮手段と、
前記差分が圧縮された前記プロファイルに対して適応フィルタを施す処理を前記各区分領域ごとに行うフィルタ処理手段と、
前記適応フィルタが施された前記プロファイル上の各データに基づいて、前記領域における他の位置の各データを復元して、二次元状に配列された各補正データを作成する補正データ作成手段と、
を備え、
前記二次元状に配列された前記各画像データから、対応する位置の前記補正データを減算処理して、前記各放射線検出素子ごとにノイズが除去された前記各画像データを生成することを特徴とする放射線画像処理装置。 - 前記フィルタ処理手段は、前記差分が圧縮された前記プロファイルに対する前記適応フィルタとして、ウィナーフィルタを施すことを特徴とする請求の範囲第1項に記載の放射線画像処理装置。
- 前記フィルタ処理手段は、前記ウィナーフィルタを、下記(1)式で表される形で前記差分が圧縮された前記プロファイルに適用し、算出された下記の値Df(n)を、前記差分が圧縮された前記プロファイル上の前記平均値から減算した値を、前記適応フィルタが施された前記プロファイル上のデータとすることを特徴とする請求の範囲第2項に記載の放射線画像処理装置。
- 前記フィルタ処理手段は、前記差分が圧縮された前記プロファイル上の前記各平均値にローパスフィルタ処理を施して、前記ウィナーフィルタを施すことを特徴とする請求の範囲第2項または第3項に記載の放射線画像処理装置。
- 前記分割手段で前記領域を少なくとも前記検出部上での前記走査線の延在方向において所定個数の前記画像データごとに区分する際の前記所定個数は、前記領域における前記走査線の延在方向の前記画像データの総数を自然数で除算した個数に設定されることを特徴とする請求の範囲第1項から第4項のいずれか一項に記載の放射線画像処理装置。
- 互いに隣接する所定個数の前記各区分領域について、前記フィルタ処理手段により前記適応フィルタが施された前記プロファイル上の各データのうち、同一の前記走査線の延在方向に存在する前記各データの平均値を算出する処理を前記各走査線ごとに行う再平均化手段を備え、
前記補正データ作成手段は、前記再平均化手段が算出した前記各データの平均値に基づいて、前記領域における他の位置の各データを復元して、二次元状に配列された各補正データを作成することを特徴とする請求の範囲第1項から第5項のいずれか一項に記載の放射線画像処理装置。 - 前記検出部上で、少なくとも前記複数の走査線が走査線の延在方向において分割されて構成されている前記放射線画像撮影装置により撮影された前記各画像データに対しては、前記複数の走査線が分割された前記検出部の各領域でそれぞれ得られた前記各画像データについてそれぞれ別々に前記分割手段から前記補正データ作成手段までの各処理を行い、それぞれ別々に二次元状に配列された前記各画像データから、対応する位置の前記補正データをそれぞれ別々に減算処理して、ノイズが除去された前記各画像データを生成することを特徴とする請求の範囲第1項から第6項のいずれか一項に記載の放射線画像処理装置。
- 前記放射線画像撮影装置により取得された前記各放射線検出素子ごとの各オフセット補正値に対しても、前記分割手段から前記補正データ作成手段までの各処理を行い、二次元状に配列された前記各オフセット補正値から、対応する位置の前記補正データを減算処理して、ノイズが除去された前記各オフセット補正値を生成し、前記各画像データからノイズが除去された前記各オフセット補正値を減算することを特徴とする請求の範囲第1項から第7項のいずれか一項に記載の放射線画像処理装置。
- 前記放射線画像撮影装置の前記各放射線検出素子ごとに基準となる前記各オフセット補正値を予め有しておき、前記放射線画像撮影装置により取得された前記各オフセット補正値から前記基準となる各オフセット補正値を減算した各差分に対して、前記分割手段から前記補正データ作成手段までの各処理を行い、二次元状に配列された前記各差分から、対応する位置の前記補正データを減算処理して、ノイズが除去された前記各差分を生成し、ノイズが除去された前記各差分を前記基準となる各オフセット補正値に加算して、前記放射線画像撮影装置の前記各放射線検出素子ごとのノイズが除去された前記各オフセット補正値を作成し、前記各画像データからノイズが除去された前記各オフセット補正値を減算することを特徴とする請求の範囲第1項から第7項のいずれか一項に記載の放射線画像処理装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/698,424 US8885909B2 (en) | 2010-05-17 | 2011-02-25 | Radiographic-image processing apparatus |
JP2012515776A JP5747913B2 (ja) | 2010-05-17 | 2011-02-25 | 放射線画像処理装置 |
EP11783314.5A EP2572641B1 (en) | 2010-05-17 | 2011-02-25 | Radiographic-image processing device |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010112739 | 2010-05-17 | ||
JP2010-112739 | 2010-05-17 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011145377A1 true WO2011145377A1 (ja) | 2011-11-24 |
Family
ID=44991495
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/054280 WO2011145377A1 (ja) | 2010-05-17 | 2011-02-25 | 放射線画像処理装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8885909B2 (ja) |
EP (1) | EP2572641B1 (ja) |
JP (1) | JP5747913B2 (ja) |
WO (1) | WO2011145377A1 (ja) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103155434B (zh) * | 2010-06-18 | 2016-08-24 | 康普技术有限责任公司 | 在分布天线***的端点间传输数字数据的方法和端点元件 |
CN103337074B (zh) * | 2013-06-18 | 2016-01-13 | 大连理工大学 | 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 |
CN104739509B (zh) * | 2013-12-31 | 2017-11-07 | 深圳迈瑞生物医疗电子股份有限公司 | 一种生理信号压缩方法、***及医疗设备 |
JP6251147B2 (ja) * | 2014-09-29 | 2017-12-20 | 富士フイルム株式会社 | 電子カセッテおよび電子カセッテの作動方法 |
JP6072102B2 (ja) * | 2015-01-30 | 2017-02-01 | キヤノン株式会社 | 放射線撮影システム及び放射線撮影方法 |
US9953430B1 (en) * | 2015-10-29 | 2018-04-24 | Indoor Reality Inc. | Methods for detecting luminary fixtures |
JP2017192605A (ja) * | 2016-04-22 | 2017-10-26 | コニカミノルタ株式会社 | 放射線画像撮影システムおよび放射線画像撮影装置 |
JP6677100B2 (ja) * | 2016-06-24 | 2020-04-08 | コニカミノルタ株式会社 | 放射線画像撮影システム |
CN106254724A (zh) * | 2016-07-29 | 2016-12-21 | 努比亚技术有限公司 | 一种实现图像降噪的方法、装置和终端 |
JP6951158B2 (ja) * | 2017-09-05 | 2021-10-20 | キヤノン株式会社 | 放射線撮像装置、その駆動方法および放射線撮像システム |
JP7087435B2 (ja) * | 2018-02-19 | 2022-06-21 | コニカミノルタ株式会社 | 放射線画像撮影装置及び放射線画像撮影システム |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0772256A (ja) * | 1993-09-01 | 1995-03-17 | Fuji Photo Film Co Ltd | 画像信号補正装置および方法 |
JP2001099944A (ja) * | 1999-07-29 | 2001-04-13 | Toshiba Medical System Co Ltd | 放射線検出器のノイズリダクション方法、放射線検出器、放射線診断装置 |
JP2003087656A (ja) * | 2001-09-13 | 2003-03-20 | Toshiba Corp | X線画像診断装置 |
JP2005003578A (ja) * | 2003-06-13 | 2005-01-06 | Shimadzu Corp | フラットパネル型検出器におけるノイズ除去方法 |
JP2005129709A (ja) * | 2003-10-23 | 2005-05-19 | Shimadzu Corp | 2次元放射線検出器 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5065444A (en) | 1988-02-08 | 1991-11-12 | Northrop Corporation | Streak removal filtering method and apparatus |
FR2686471B1 (fr) * | 1992-01-17 | 1994-03-04 | Thomson Tubes Electroniques | Procede d'elaboration d'une image matricielle avec filtrage de valeurs moyennes de lignes de pixels. |
US5661309A (en) | 1992-12-23 | 1997-08-26 | Sterling Diagnostic Imaging, Inc. | Electronic cassette for recording X-ray images |
EP0712092A1 (en) * | 1994-11-10 | 1996-05-15 | Agfa-Gevaert N.V. | Image contrast enhancing method |
JP3486490B2 (ja) | 1995-09-04 | 2004-01-13 | キヤノン株式会社 | 放射線検出装置 |
DE69631126T2 (de) * | 1996-08-08 | 2004-09-16 | Agfa-Gevaert | Verfahren zur Verbesserung der Aufzeichnungsmaterialfehler von Strahlungsbildern |
JP3765155B2 (ja) | 1997-03-28 | 2006-04-12 | コニカミノルタホールディングス株式会社 | 放射線画像読取装置 |
JP3890163B2 (ja) | 1999-04-27 | 2007-03-07 | キヤノン株式会社 | 撮影システム |
US7221782B1 (en) * | 1999-06-24 | 2007-05-22 | General Electric Company | Method and apparatus for determining a dynamic range of a digital medical image |
JP2005303448A (ja) | 2004-04-07 | 2005-10-27 | Konica Minolta Medical & Graphic Inc | 放射線画像処理装置、放射線画像読取装置、放射線画像処理方法及び放射線画像読取方法 |
JP4012182B2 (ja) | 2004-08-19 | 2007-11-21 | キヤノン株式会社 | カセッテ型x線画像撮影装置 |
EP1875855B1 (en) * | 2005-04-27 | 2012-05-02 | Olympus Medical Systems Corp. | Image processing apparatus, image processing method, and image processing program |
US20070014468A1 (en) * | 2005-07-12 | 2007-01-18 | Gines David L | System and method for confidence measures for mult-resolution auto-focused tomosynthesis |
US7474731B2 (en) * | 2006-08-29 | 2009-01-06 | Siemens Medical Solutions Usa, Inc. | Systems and methods for adaptive image processing using acquisition data and calibration/model data |
-
2011
- 2011-02-25 JP JP2012515776A patent/JP5747913B2/ja not_active Expired - Fee Related
- 2011-02-25 WO PCT/JP2011/054280 patent/WO2011145377A1/ja active Application Filing
- 2011-02-25 EP EP11783314.5A patent/EP2572641B1/en not_active Not-in-force
- 2011-02-25 US US13/698,424 patent/US8885909B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0772256A (ja) * | 1993-09-01 | 1995-03-17 | Fuji Photo Film Co Ltd | 画像信号補正装置および方法 |
JP2001099944A (ja) * | 1999-07-29 | 2001-04-13 | Toshiba Medical System Co Ltd | 放射線検出器のノイズリダクション方法、放射線検出器、放射線診断装置 |
JP2003087656A (ja) * | 2001-09-13 | 2003-03-20 | Toshiba Corp | X線画像診断装置 |
JP2005003578A (ja) * | 2003-06-13 | 2005-01-06 | Shimadzu Corp | フラットパネル型検出器におけるノイズ除去方法 |
JP2005129709A (ja) * | 2003-10-23 | 2005-05-19 | Shimadzu Corp | 2次元放射線検出器 |
Non-Patent Citations (1)
Title |
---|
See also references of EP2572641A4 * |
Also Published As
Publication number | Publication date |
---|---|
EP2572641A4 (en) | 2013-12-11 |
US20130071000A1 (en) | 2013-03-21 |
JP5747913B2 (ja) | 2015-07-15 |
US8885909B2 (en) | 2014-11-11 |
EP2572641B1 (en) | 2018-01-24 |
JPWO2011145377A1 (ja) | 2013-07-22 |
EP2572641A1 (en) | 2013-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5747913B2 (ja) | 放射線画像処理装置 | |
JP7085043B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP6214226B2 (ja) | 画像処理装置、断層撮影装置、画像処理方法およびプログラム | |
JP5838702B2 (ja) | 放射線画像撮影システム | |
JP6609119B2 (ja) | 放射線撮影装置、放射線撮影方法、放射線撮影システム、及びプログラム | |
KR20150122392A (ko) | 엑스선 영상 장치 및 엑스선 영상 장치의 영상 처리 방법 | |
JP2011193306A (ja) | 放射線画像撮影装置および放射線画像撮影システム | |
US20090304151A1 (en) | Radiographic image data correction method and apparatus and radiography apparatus | |
CN111163698A (zh) | 图像处理设备、图像处理方法以及程序 | |
JP2008073342A (ja) | 放射線画像撮影システム及び放射線画像撮影方法 | |
JP4497615B2 (ja) | 画像処理装置、補正方法及び記録媒体 | |
JP7063199B2 (ja) | 放射線画像撮影システム | |
JP2020149237A (ja) | 画像処理装置、放射線撮影装置および画像処理方法 | |
JP2017077405A (ja) | 放射線撮影システム、画像処理装置及び画像処理方法 | |
JP2005028114A (ja) | 放射線撮影装置及び放射線撮影方法 | |
US7615756B2 (en) | Apparatus for and method of capturing radiation image | |
JP4980769B2 (ja) | 放射線撮像装置及び方法 | |
JP4754812B2 (ja) | X線撮影装置 | |
JP4258092B2 (ja) | 画像処理装置および画像処理方法 | |
US20230404510A1 (en) | Radiation image processing device, radiation image processing method, and radiation image processing program | |
JP2013138280A (ja) | 放射線画像撮影装置および放射線画像撮影システム | |
JP6739511B2 (ja) | 放射線画像撮影装置および放射線画像撮影システム | |
JP2001222704A (ja) | 画像処理方法及び画像処理装置 | |
JP2530223B2 (ja) | 放射線撮影の補正方法 | |
JP2019084307A (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: 11783314 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012515776 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011783314 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13698424 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |