WO2019064632A1 - X線撮像装置およびx線撮像素子の画像処理方法 - Google Patents

X線撮像装置およびx線撮像素子の画像処理方法 Download PDF

Info

Publication number
WO2019064632A1
WO2019064632A1 PCT/JP2018/008420 JP2018008420W WO2019064632A1 WO 2019064632 A1 WO2019064632 A1 WO 2019064632A1 JP 2018008420 W JP2018008420 W JP 2018008420W WO 2019064632 A1 WO2019064632 A1 WO 2019064632A1
Authority
WO
WIPO (PCT)
Prior art keywords
light receiving
ray
pixel
receiving element
ray imaging
Prior art date
Application number
PCT/JP2018/008420
Other languages
English (en)
French (fr)
Inventor
桜井 健次
文洋 趙
Original Assignee
国立研究開発法人物質・材料研究機構
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 国立研究開発法人物質・材料研究機構 filed Critical 国立研究開発法人物質・材料研究機構
Priority to JP2019544217A priority Critical patent/JPWO2019064632A1/ja
Publication of WO2019064632A1 publication Critical patent/WO2019064632A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/22Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
    • G01N23/223Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material by irradiating the sample with X-rays or gamma-rays and by measuring X-ray fluorescence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/24Measuring radiation intensity with semiconductor detectors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36Measuring spectral distribution of X-rays or of nuclear radiation spectrometry

Definitions

  • the present invention relates to an X-ray imaging apparatus using a CCD imaging device, a CMOS imaging device or a semiconductor pixel radiation detector, and more specifically, a general purpose CCD imaging device for visible light or a CMOS imaging device, or a semiconductor pixel radiation detector
  • An X-ray imaging apparatus such as a fluorescent X-ray analyzer using the
  • the present invention also relates to an image processing method of an X-ray imaging device using a CCD imaging device, a CMOS imaging device, or a semiconductor pixel radiation detector.
  • Substances are composed of various elements, and their composition greatly affects their physical and chemical properties. For this reason, it is important to analyze the types and amounts of contained elements in order to understand substances and develop new materials. It is known that the type of an element is known from the energy of fluorescent X-rays emitted when the substance is irradiated with X-rays, and the amount is known from its intensity (Patent Documents 1 and 2). In order to perform this fluorescent X-ray analysis, a dedicated X-ray spectrometer and an X-ray detector are required together with the X-ray source.
  • Charge-coupled device (CCD) cameras usually have a very large active sensor area of more than 150 mm 2 .
  • Conventional CCD cameras are designed to record visible light images in a variety of applications, including scientific and industrial fields.
  • the quantum efficiency of the x-ray region is low due to the thin depletion of the CCD camera, but can also work in the x-ray wavelength region. Since the 1980's, advanced x-ray cameras such as pnCCD have been developed for astronomy. Since the thickness of the depletion layer is typically several hundred ⁇ m, the quantum efficiency of the X-ray region (1 to 10 keV) can be 90% or more.
  • the present invention solves the above-mentioned problems, and can separate X-ray energy information individually by processing the X-ray photon event occurring in a CCD image sensor, a CMOS image sensor or a semiconductor pixel radiation detector.
  • An object of the present invention is to provide an X-ray imaging apparatus capable of obtaining precise imaging pixel information.
  • Another object of the present invention is X-ray imaging suitable for utilizing a CCD imaging device or a CMOS imaging device designed for visible light wavelength, or a semiconductor pixel radiation detector as an imaging device of a large area X-ray spectrometer It is providing a device.
  • the X-ray imaging apparatus of the present invention is an X-ray imaging device (8) that receives an energy beam containing X-ray photons from an analysis object (5) as shown in, for example, FIG. 1 and FIG. 2A.
  • the X-ray imaging device has the two-dimensional arrangement of the unit light receiving devices (82) and the peripheral light receiving device (84) around the unit light receiving devices that causes charge sharing due to light reception of the X-ray photons.
  • An X-ray imaging device The light reception signal from the X-ray imaging device is read for each unit light receiving device, and an X-ray image of the analysis object captured according to the two-dimensional arrangement of the unit light receiving device in the X-ray imaging device An X-ray image generation circuit (10) to generate;
  • the X-ray image generation circuit 10 an effective event determination circuit (102) comparing the light reception signal read for each unit light receiving element and a first threshold value as a reference for recognizing an effective photon event; With respect to the peripheral light receiving element of the unit light receiving element recognized as an effective photon event by the effective event judging circuit, the light receiving signal read by the peripheral light receiving element is compared with a second threshold serving as a reference for recognizing noise.
  • a noise removing circuit (104) and a pixel for performing counting such that the pixel-by-pixel counter value of the X-ray image generating circuit (10) is incremented by 1 with respect to unit light receiving elements recognized as effective photon events without being removed by the noise removing circuit.
  • the peripheral light receiving element of the unit light receiving element recognized as a valid photon event and recognized by the valid event judging circuit (102) Isolated pixel extraction in which the light receiving signal of all the peripheral light receiving elements is compared with the first threshold and the unit light receiving element whose light receiving signal of all the peripheral light receiving elements is less than the first threshold is extracted as an isolated pixel
  • the noise removal circuit (104) has a circuit (103), and is a unit light receiving element recognized as a valid photon event and recognized by the valid event determination circuit (102), and the noise removal circuit (104) comprises the isolated pixel extraction circuit (103).
  • the peripheral light receiving element of the unit light receiving element extracted as an isolated pixel is characterized by comparing the light receiving signal read by the peripheral light receiving element with a second threshold value.
  • an X-ray imaging device (8) which receives an energy beam containing X-ray photons from an analysis object (5).
  • the X-ray imaging device has the two-dimensional arrangement of the unit light receiving devices (82) and the peripheral light receiving device (84) around the unit light receiving devices that causes charge sharing due to light reception of the X-ray photons.
  • An X-ray imaging device The light reception signal from the X-ray imaging device is read for each unit light receiving device, and an X-ray image of the analysis object captured according to the two-dimensional arrangement of the unit light receiving device in the X-ray imaging device And an X-ray image generation circuit (10) to generate The X-ray image generation circuit 10, an effective event determination circuit (102) comparing the light reception signal read for each unit light receiving element and a first threshold value as a reference for recognizing an effective photon event; With respect to the peripheral light receiving element of the unit light receiving element recognized as an effective photon event by the effective event judging circuit (102), the light receiving signal read by the peripheral light receiving element is compared with the first threshold value.
  • a charge division detection circuit (107) for detecting whether charge division has occurred between the unit light receiving element and the peripheral light receiving element; a unit light receiving element having charge division detected by the charge division detection circuit (107); Charge division processing circuit (108) for calculating median value, integral value and position information of peripheral light receiving elements for all light receiving elements exceeding the first threshold, and charge division processing circuit 108) counting the count value of the pixel-by-pixel counter corresponding to all the light receiving elements exceeding the first threshold value, using the median value, integral value and position information thereof calculated in 108). And an X-ray counter image generation circuit (108) for generating an X-ray image by the X-ray imaging device (8) based on the count value of the pixel-by-pixel counter circuit (106). It is characterized by
  • the X-ray imaging device (8) is a CCD imaging device, a CMOS imaging device, or a semiconductor pixel radiation detector.
  • the unit light receiving element (82) of the X-ray imaging element (8) has a two-dimensional area having a representative length of 1 to 10 ⁇ m as one pixel. It is good to do.
  • the first threshold is a value larger than background noise generated by each unit light receiving element (82) of the X-ray imaging element (8). Therefore, the value should be smaller than the energy level generated by the effective photon event.
  • the second threshold is a value larger than the background noise generated by each unit light receiving element (82) of the X-ray imaging element (8), and is a value smaller than the first threshold Good to have.
  • the peripheral light receiving element 84 is any of 3 ⁇ 3, 5 ⁇ 5 or 7 ⁇ 7 centered on the unit light receiving element (82) in which a photon event has occurred. It is good to have.
  • the energy level of the light receiving signal read for each unit light receiving element (82) of the X-ray image generation circuit (10) It is preferable to have an energy level determination circuit (105) for obtaining [10]
  • the pixel-by-pixel counter circuit of the X-ray image generation circuit (10) preferably corresponds to the energy level separated by the energy level determination circuit (105). It is good to be divided.
  • the pixel-by-pixel counter circuit of the X-ray image generation circuit (10) is an energy level of the light reception signal of the unit light receiving element (82) where a photon event occurred.
  • the energy levels of the light reception signals of the peripheral light reception elements (84) of the unit light reception element (82) may also be added.
  • the X-ray imaging apparatus of the present invention preferably further comprises an X-ray irradiation circuit that irradiates an energy beam containing X-rays to the object to be analyzed.
  • the X-ray imaging apparatus of the present invention further includes a pixel-by-pixel counter reset circuit that resets the pixel-by-pixel counter value of the X-ray image generation circuit (10).
  • the image processing method of the X-ray imaging device of the present invention is, for example, as shown in FIG. 5, FIG. 7A, and FIG. 7B, an X-ray imaging device that receives energy rays including X-ray photons from an analysis object.
  • An image processing method wherein the X-ray imaging device (8) has unit light receiving elements arranged in a two-dimensional arrangement, and the unit light receiving elements have a two-dimensional area which causes charge sharing by receiving the X-ray photons.
  • the light receiving signal read by the peripheral light receiving element is compared with the first threshold value, and the light receiving signals of all the peripheral light receiving elements are less than the first threshold value.
  • An isolated pixel extracting step (S200) for extracting a unit light receiving element as an isolated pixel, and the noise removing step (S300) is a unit light receiving element recognized as an effective photon event recognized in the effective event determining step (S100).
  • the light receiving signal of the unit light receiving element extracted as an isolated pixel in the isolated pixel extracting step (S200) is recognized as a light receiving signal read by the peripheral light receiving element and noise And comparing the second threshold value as a quasi.
  • a charge division processing step for calculating a median, an integral value, and position information thereof, and in the pixel-by-pixel counter step, the median, integral value calculated in the charge division processing step And the position information using, and performs counting to plus 1 the count value of the pixel-counter corresponding to all the light-receiving element exceeds the first threshold.
  • information on X-ray energy is obtained using a CCD imaging element, a CMOS imaging element or a semiconductor pixel radiation detector designed for visible light wavelength. Can be separated individually to obtain precise imaging pixel information.
  • a CCD imaging device, a CMOS imaging device, or a semiconductor pixel radiation detector can cope with a large-area X-ray spectrometer, the preconditions of spectral imaging can be established, and an X utilizing a CCD camera or a CMOS imaging device sufficiently. Line spectral imaging can be performed.
  • FIG. 2 is a block diagram of a configuration of a correction filtering circuit for isolated pixels, showing a first aspect of the present invention.
  • FIG. 7 is a block diagram of a correction filtering circuit for isolated pixels, which is a slight modification of the first aspect of the present invention.
  • FIG. 7 is a diagram specifically showing a method of obtaining X-ray energy (wavelength) from the brightness (charge amount) of the detected pixel, and particularly shows the X-ray imaging device and the distribution state of the charge amount.
  • FIG. 2 is a block diagram of a configuration of a correction filtering circuit for isolated pixels, showing a first aspect of the present invention.
  • FIG. 7 is a block diagram of a correction filtering circuit for isolated pixels, which is a slight modification of the first aspect of the present invention.
  • FIG. 7 is a diagram specifically showing a method of obtaining X-ray energy (wavelength) from the brightness (charge amount) of the detected pixel, and particularly shows the X-ray imaging device and the distribution state of the charge
  • FIG. 2 is a diagram specifically showing a method of obtaining X-ray energy (wavelength) from the brightness (charge amount) of a detected pixel, and particularly shows the charge amount and the count number. It is explanatory drawing about division (general charge division) of the electric charge which arose by the X-ray photon into several pixels. It is a flow chart explaining initial selection of a pixel which should be processed. It is a flowchart which shows the 1st type
  • FIG. 10 is a flow chart illustrating the strict threshold addition (second-order removal of the effect of charge division).
  • FIG. 17 illustrates the introduction of a strict threshold (secondary removal of the influence of charge division), showing the energy level of the threshold 120 for isolated pixel determination in the amount of charge.
  • FIG. 6 shows X-ray images by X-ray spectral energy using the method of strictly selecting and collecting the charges of isolated pixels (FIGS. 5 to 7B). It is an explanatory view showing a second type which is an object of the present invention and schematically showing a situation where charge division is the majority and no isolated pixel can be found. It is explanatory drawing of the calculation principle of the total charge amount including a surrounding pixel.
  • FIG. 6 is a block diagram of a correction filtering circuit for charge division, showing a second aspect of the present invention.
  • FIG. 10 shows X-ray images of different X-ray spectral energies using the method of correcting the influence of charge division (FIG. 9, 10) using the ratio of the median to the integral.
  • FIG. 10 is a block diagram of a configuration of a correction filtering circuit of a combined use type of isolated pixel and charge division showing a third aspect of the present invention. It is a principal part block diagram of the measuring apparatus for acquiring the X-ray image according to X-ray spectrum energy which shows the 2nd Example of this invention. It is explanatory drawing of the X-ray counter image by energy level produced
  • FIG. 1 is a schematic view of an X-ray analyzer showing an embodiment of the present invention.
  • the X-ray analyzer according to one embodiment of the present invention includes an X-ray tube 1, a crystal monochromator 2, a signal control mechanical shutter 3, a sample holder 4, a pinhole plate 6, and an X-ray camera 7.
  • the X-ray tube 1 is an electron tube for generating X-rays, and continuous X-rays are generated by colliding electrons generated by heating the filament (cathode) with a metal (target, anode) such as tungsten, molybdenum or copper. (Braking radiation) and characteristic X-rays are generated.
  • the monochromator 2 is also called a monochromator, and is a device for taking out only light of a specific wavelength from the light dispersed by a dispersion element such as a diffraction grating or a prism with a slit.
  • the monochromator includes, for example, an entrance slit, a collimator mirror, a dispersive element (diffraction grating or prism), a focusing mirror, and an exit slit.
  • graphite (002), 2d 6.72 ⁇ is used.
  • a curved graphite spectroscope (monochromator) is installed at a position in front of a detector to remove a diffraction line due to a CuK ⁇ ray.
  • a thin Ni foil filter may be used instead of the monochromator.
  • the signal control mechanical shutter 3 physically moves the shutter curtain, and in fact, it is necessary to block incident X-rays from being incident on the sample. Therefore, while using the same mechanism as the mechanical shutter for visible light, one using a material that can absorb X-rays and a thickness can be used, and UNIBLITZ (registered trademark) can be used, for example.
  • the necessity of the mechanical shutter is determined by the magnitude relationship between the actual imaging time and the time required to transfer the image. In the CCD sensor shown in Example 1 to be described later, the readout transfer time of one frame is affected, and the effective imaging time is different between the pixel at the beginning of reading and the pixel at the end of reading. In order to prevent this, it is necessary to stop the X-ray with a mechanical shutter. On the other hand, in the CMOS sensor of the second embodiment, since the readout is quick, the measurement is performed without using the mechanical shutter depending on the imaging time.
  • the sample holder 4 holds a sample 5 to be subjected to X-ray spectral analysis, and has, for example, a sample stage.
  • the sample 5 is a target of X-ray spectral analysis, and includes various test objects such as a metal material, an inorganic material, an organic material, and a biomaterial.
  • pinhole plates are used very often in XRF imaging.
  • the pinhole plate 6 is also referred to as a micro pinhole collimator and is used in the case of XRF imaging measurement.
  • the pinhole plate 6 is removed in non-imaging XRF spectrum measurement.
  • the X-ray camera 7 is provided with an imaging sensor 8.
  • the imaging sensor 8 is a two-dimensional semiconductor sensor, and for example, a CCD, a CMOS, or an X-ray pixel detector in which an element structure is specially formed is used.
  • the imaging sensor 8 has unit light receiving elements of a two-dimensional arrangement according to the number of pixels. For example, in the case where the number of pixels is 1,000,000, 1000 unit light receiving elements are two-dimensionally arranged in the vertical x horizontal directions. For example, when the imaging sensor 8 is 10 mm ⁇ 10 mm, each unit light receiving element has a size of 10 ⁇ m ⁇ 10 ⁇ m.
  • a CCD 47-10 which is a CCD sensor manufactured by e2V can be used.
  • the beam size at the sample position is about 1 mm (H), 5 mm (V).
  • the cooler 9 is used when the imaging sensor 8 is used in an X-ray application.
  • the cooling of the cooler 9 is generally electronic cooling by a Peltier element, but its exhaust heat needs to be removed by air cooling or water cooling.
  • air cooling a fan is used as the cooler 9.
  • water cooling the cooler 9 requires circulating cooling water.
  • FIG. 1 schematically shows the case of water cooling.
  • the cooling temperature is as wide as -50 ° C to 5 ° C depending on the sensor. Below 0 ° C., it is general to evacuate the area around the sensor to avoid freezing and the like.
  • a cooled CCD is used at -30.degree. C.
  • a CMOS camera is used at 5.degree.
  • the X-ray beam 101 a indicates that the X-ray generated by the X-ray tube 1 is irradiated to the sample 5 through the monochromator 2.
  • the X-ray beam 101a spreads horizontally on the sample 5, and the typical irradiation area is 10 mm ⁇ 5 mm.
  • Primary X-ray intensity is greater than 10 8 counts / sec.
  • the X-ray beam 101 b indicates one of the X-rays reflected on the sample 5 that passes through the pinhole plate 6 and the CCD sensor 7.
  • the operation of the device thus configured will now be described.
  • the X-ray beam 101 a strikes the sample at a low angle, the X-ray beam 101 a spreads horizontally on the sample 5.
  • a signal-controlled mechanical shutter 3 is placed in the beam path from the x-ray tube 1 to the sample 5.
  • the mechanical shutter 3 is open and allows the X-ray beam 101a to pass.
  • the mechanical shutter 3 is closed to shut off the X-ray beam 101a to avoid contamination on the camera image.
  • the X-ray beam 101 b reflected on the sample 5 is narrowed by the pinhole plate 6 and reaches the X-ray camera 7.
  • the centers of the sample holder 4, the pinhole plate 6, and the CCD sensor 7 are aligned.
  • the X-ray camera 7 can be moved back and forth in the optical axis direction, and may be provided with a position adjustment mechanism which can also be moved in a planar direction perpendicular to the optical axis.
  • the pixel detection signal received by the imaging sensor 8 is converted into an image signal by the X-ray image generation circuit 10.
  • FIG. 2A is a configuration block diagram of a correction filtering circuit for isolated pixels, showing the first aspect of the present invention.
  • the X-ray imaging device 8 receives energy beams including X-ray photons from the object to be analyzed 4 and has unit light receiving devices 82 in a two-dimensional arrangement.
  • a peripheral light receiving element 84 is generated around the unit light receiving element 82 to cause charge sharing due to the reception of the X-ray photon.
  • the peripheral light receiving element 84 may be any of 3 ⁇ 3, 5 ⁇ 5 or 7 ⁇ 7 centered on the unit light receiving element 82 in which the photon event has occurred.
  • the X-ray image generation circuit 10 reads the light reception signal from the X-ray imaging device 8 for each unit light receiving device, and picks up an image according to the two-dimensional arrangement of the unit light receiving device in the X-ray imaging device. It is generated as an X-ray image of the analyte.
  • the X-ray image generation circuit 10 includes an effective event determination circuit 102, an isolated pixel extraction circuit 103, a noise removal circuit 104, an energy level determination circuit 105, a pixel-by-pixel counter circuit 106, and an X-ray counter image generation circuit 109.
  • the components of the X-ray image generation circuit 10 may use hardware electronic circuits such as application specific integrated circuits for each of the components or may realize functions by computer software. FIG.
  • FIG. 2B shows a slight modification of the first aspect of the invention of FIG. 2A.
  • the pixel-by-pixel counter reset circuit 106 a is added below the pixel-by-pixel counter circuit 106.
  • the other configurations are the same, and thus redundant description will be omitted.
  • counting is performed to add 1 to the count value of the pixel-by-pixel counter corresponding to all light receiving elements exceeding the first threshold, but when the pixel to be measured is completed, the pixel-by-pixel counter
  • the pixel-by-pixel counter is reset by the reset circuit 106 a to end the pixel-by-pixel processing, and the X-ray counter image generation circuit 109 generates an image.
  • the valid event judging circuit 102 compares the light reception signal read for each unit light receiving element 82 with a first threshold value which is a reference for recognizing a valid photon event.
  • the first threshold is a value larger than the background noise generated in each unit light receiving element 82 of the X-ray imaging element 8 and smaller than the energy level generated in the effective photon event.
  • the first threshold is also called Lower Limit of Discrimination, and is used not only for CCD and CMOS sensors but also for signal processing circuits of any radiation detector. The purpose is to remove the kind of noise or dark current of the detector or the signal processing circuit that is also recognized when X-rays are not incident.
  • the signal is recognized only when the signal level is higher than the first threshold value, and the value of the luminance value (the brightness of the image, which corresponds to the charge amount generated by the X-ray) seen in the pixel is considered.
  • this luminance value corresponds to the information of X-ray energy, and in the aspect of elemental analysis, the information of elements.
  • the isolated pixel extraction circuit 103 detects the light receiving signal read by the peripheral light receiving element and the first threshold value for the peripheral light receiving element 84 of the unit light receiving element 82 recognized as an effective photon event by the effective event judging circuit 102. And the unit light receiving element 82 in which the light reception signals of all the peripheral light receiving elements 84 are equal to or less than the first threshold value is extracted as an isolated pixel.
  • the noise removing circuit 104 uses the light receiving signal read by the peripheral light receiving element 84 for the peripheral light receiving element 84 of the unit light receiving element 82 recognized as a valid photon event by the valid event judging circuit 102 and a reference for recognizing as noise. And a second threshold value.
  • the second threshold is a value larger than the background noise generated in each unit light receiving element 82 of the X-ray imaging device 8 and is a value smaller than the first threshold.
  • the energy level determination circuit 105 is also called a spectral circuit, and obtains the energy level of the light reception signal read from each unit light receiving element 82 of the X-ray image generation circuit 10.
  • the energy level is set, for example, to a value corresponding to an element to be analyzed, for example, a metal element such as calcium or chromium.
  • the energy level may have a region defined by the upper limit value and the lower limit value.
  • the pixel-by-pixel counter circuit 106 counts the pixel-by-pixel counter value of the X-ray image generation circuit 10 by 1 for unit light-receiving elements 82 that are not removed by the noise removal circuit 104 and are recognized as valid photon events.
  • the pixel-by-pixel counter circuit 106 is preferably divided according to the energy level separated by the energy level determination circuit 105, and the energy level of the light reception signal of the unit light receiving element 82 where a photon event has occurred.
  • the energy levels of the light reception signals of the peripheral light reception elements 84 of the unit light reception element 82 may also be added.
  • the X-ray counter image generation circuit 109 generates an X-ray image by the X-ray imaging device 8 based on the count value of the pixel-by-pixel counter circuit 106.
  • FIG. 16 is an explanatory diagram of an X-ray counter image by energy level generated by the X-ray counter image generation circuit, where (A1) is an element A image, (A2) is an element B image, and (A3) is an element C An image, (A4) represents an element D image, (B1) represents a region P spectrum, (B2) represents a region Q spectrum, and (B3) represents a region R spectrum.
  • the level of the brightness (charge amount I) of the (X, Y) pixel is in the range of X-ray photon energy corresponding to the element A ( AL ⁇ I ⁇ A H ) ,
  • the count of the (X, Y) position of the image of the A element is +1.
  • the count of the (X, Y) position of the image of the B element is incremented by 1 with respect to the X-ray photon energy range (B L ⁇ I ⁇ B H ) corresponding to the element B .
  • the count of the (X, Y) position of the image of the C element is incremented by 1 with respect to the X-ray photon energy range (C L ⁇ I ⁇ C H ) corresponding to the element C .
  • the count of the (X, Y) position of the D element image is incremented by +1. .
  • the photon event is displayed at the position of coordinate (X, Y) on FIG. 16 (A2). According to the region P spectrum shown in FIG.
  • X-ray photon energy E by the level of the brightness (charge amount I) of each pixel can be found for the (X, Y) region (or the entire region). That is, since the X-ray spectrum in the (X, Y) region is known, the count of the energy E of the spectrum of the region P spectral region is incremented by 1.
  • the region Q spectrum shown in FIG. 16 (B2) and the region R spectrum shown in FIG. 16 (B3) are also similar to the region P spectrum shown in FIG. 16 (B1).
  • the location of the brightness (charge amount) of the pixel corresponding to the photon event is +1 in the region Q spectrum in FIG. 16 (B2).
  • FIGS. 3A and 3B specifically show a method of obtaining X-ray energy (wavelength) from the brightness (charge amount) of the detected pixel
  • FIG. 3A shows the X-ray imaging device and the distribution state of the charge amount
  • FIG. 3B shows the charge amount and the count number.
  • the unit light receiving elements 110 are unit light receiving elements of the X-ray imaging element, and are two-dimensionally arranged.
  • the first threshold 112 is a value larger than background noise generated by the unit light receiving element 110 and smaller than an energy level generated by an effective photon event.
  • the noise equivalent signal 111 is such that the brightness (charge amount) of the pixel detected by the unit light receiving element 110 is lower than the first threshold 112.
  • the effective photon event signal 113 is such that the brightness (charge amount) of the pixel detected by the unit light receiving element 110 is higher than the first threshold 112.
  • the horizontal axis indicates the charge amount
  • the vertical axis indicates the count number.
  • the first threshold 115 is a threshold sufficiently larger than the amount of charge 117 at the level of dark current seen with or without X-ray irradiation.
  • the effective photon event signal 116 is a signal with a larger amount of charge than the first threshold 115, and can be summed to obtain information on the charge generated by the X-ray.
  • the amount of charge at the level of the dark current has a median value 118 of the dark current and has a noise width 119 with respect to the upper limit value of the actual dark current.
  • the imaging under the single photon counting condition is, in other words, imaging that is repeatedly and continuously performed imaging in a very short imaging time.
  • the imaging time is increased, the accumulated image becomes clear, but the meaning is that a large number of charges are accumulated in one pixel.
  • the imaging time is increased, the information of individual X-ray photons can not be known anymore, by which X-ray of X-ray photon energy produces the charge.
  • the X-ray energy can be determined by calculating the charge amount back.
  • FIG. 4 is an explanatory diagram of division of charge generated by X-ray photons into a plurality of pixels (general charge division, charge sharing).
  • an ordinary radiation detector semiconductor detector
  • the total charge amount 125 is shortened for a short time, as described in FIGS. 3A and 3B.
  • Energy analysis 122 can be measured.
  • the X-ray imaging device 124 is composed of a large number of unit image elements 123, the charge generated by one X-ray photon 121 does not necessarily go to one device.
  • FIG. 5 is a flowchart illustrating initial selection of pixels to be processed.
  • the initial selection (S100) of the pixel to be processed corresponds to the operation of the valid event determination circuit 102 of the X-ray image generation circuit 10.
  • the input data is a value recorded in all the pixels (unit light receiving elements) of the X-ray imaging device (S102).
  • the pixels (the unit light receiving elements 82) exceeding the first threshold (the lowest threshold) are sorted and counted (S104).
  • a histogram of charge amount (X-ray spectrum) is obtained (S106).
  • FIG. 6 shows a first type targeted by the present invention, and is a flow chart for explaining sorting of isolated pixels (first-order removal of the influence of charge division).
  • the selection of isolated pixels (S200) corresponds to the operation of the isolated pixel extraction circuit 103 of the X-ray image generation circuit 10.
  • the input image data is an X-ray image obtained by subtracting the background seen even when there is no X-ray (S202).
  • isolated pixels are selected (S204). That is, the peripheral pixels (peripheral light receiving element 84) of the pixel (unit light receiving element 82) for which a valid signal is recognized are checked to confirm that any peripheral pixels are equal to or less than the minimum threshold.
  • the output data is a histogram (X-ray spectrum) of the charge amount indicating the position of the isolated pixel and the charge amount (S206). That is, it is determined whether the pixel is an isolated pixel, and only the isolated pixel is selected. In this way, the influence of division of charge into multiple pixels in the X-ray spectrum is removed from the charge amount histogram.
  • X-ray spectrum X-ray spectrum
  • FIGS. 7A and 7B illustrate the introduction of a strict threshold (secondary removal of the influence of charge division), and FIG. 7A is a flowchart, and FIG. 7B is an energy of the threshold 120 for isolated pixel determination in charge amount. It shows the level.
  • the strict threshold addition (S300) corresponds to the operation of the noise removal circuit 104 of the X-ray image generation circuit 10.
  • the input image data is the position and charge amount of the isolated pixel (S302).
  • a threshold for isolated pixel determination (second threshold 120 shown in FIG. 7B) having a value smaller than the normally adopted minimum threshold is applied to pixels in the peripheral portion, and even one of the peripheral pixels is used. If there is something beyond this, the charge collection is regarded as incomplete and excluded (S304).
  • the output data is a histogram of charge (X-ray spectrum) indicating the position of the isolated pixel and the amount of charge, and the noise that is determined to be an incomplete charge collection noise due to the influence of charge division is removed.
  • the obtained data is obtained (S306).
  • the effective event determination circuit 102 and the isolated pixel extraction circuit 103 first search for an isolated pixel for determination.
  • the isolated pixel found has a luminance value larger than the first threshold, and eight pixels around the isolated pixel are smaller than the first threshold It is a matter of course from the definition of isolated pixels that the luminance value is obtained.
  • the noise removing circuit 104 imposes a more severe restriction on the peripheral eight pixels. It is a feature of the first type of the present invention.
  • the peripheral light receiving element 84 is adopted as an isolated pixel having eight effective peripheral pixels if it is smaller than the first threshold value but smaller than the more severe second threshold value, and if not applicable. , Charge division of the second type.
  • the (X, Y) pixel is a counter, and one is added to a specific point to brighten that point. Since the pixel-by-pixel counter in this case is in a position where it only needs to exceed the first threshold, it ignores the energy of the X-rays and gives an image incorporating all the X-ray spectra.
  • the energy level determination circuit 105 is used to determine whether the luminance value of the pixel is in the range corresponding to the specific X-ray energy. For example, in a CCD camera equipped with a 16-bit A / D converter, the brightness of a pixel of an image is expressed by an integer value between 0 and 65535.
  • the first threshold value Assuming that the brightness of the dark current level is 250, it is customary to set the first threshold value to, for example, 300 arbitrary values. Also, the second threshold is set, for example, between 250 and 300. Most X-ray signals have a level corresponding to an energy level much higher than the first threshold, and even a single photon has a value of, for example, about 800. At this time, the difference between X-ray energy and whether it was 600 or 800 is in the range of the characteristic value according to the element usually composed. Therefore, when creating an image of a specific element using the energy level determination circuit 105, for example, the setting is made to be 750 or more and 820 or less, and the luminance value is so at point (X, Y).
  • the (X, Y) point of the image of the element is added to 1. That is, when measuring, the image is only one image, but the image corresponding to the counter may be prepared and used independently for the number of elements to be focused and the type of X-ray energy focused on. it can.
  • FIG. 8 shows an X-ray image by X-ray spectral energy using the method of strictly selecting and collecting charges of isolated pixels (FIGS. 5 to 7 B), using the setup of FIG. There is.
  • FIG. 8 shows pinhole XRF imaging obtained at a spatial resolution of 20 ⁇ m
  • FIG. 8 (a) shows a photograph of the inspection object, and the observation area is indicated by a dashed square.
  • FIG. 8 (b) shows the obtained XRF spectrum.
  • the imaging results of Ca in the substrate and Cr of the black pattern are shown in FIGS. 8 (c) and 8 (d), respectively. Cr bars at 20 ⁇ m intervals are clearly distinguished.
  • the X-ray tube 1 is a sealed X-ray tube, and a copper target (TOSHIBA A26L-Cu, 1.5 kW) is used.
  • Characteristic X-rays generated from the copper target include K ⁇ rays of about 1.54 ⁇ and K ⁇ rays of 1.38 ⁇ .
  • the K ⁇ line is a mixture of the K ⁇ 1 line with a peak wavelength of 1.5406 ⁇ and the K ⁇ 2 line with 1.5444 ⁇ at an intensity ratio of 2: 1, but since the two peaks are in close proximity, a standard powder diffractometer The specification is to detect both without distinguishing them.
  • a tungsten foil having a thickness of 20 ⁇ m to 50 ⁇ m and a hole of 5 ⁇ m to 100 ⁇ m in diameter is used by a UV laser. Having a spatial resolution of more than 20 ⁇ m can be visualized directly, for example by testing a standard resolution target. Because pinholes have advantages such as low cost, easy adjustment, and energy dependence, obtaining high spatial resolution with pinholes is important to expand the application of full field of view XRF imaging.
  • the X-ray camera 7 is combined with the collimator of the micro pinhole plate 6 for full field of view XRF imaging to achieve a spatial resolution better than 20 ⁇ m. This result can be directly visualized without mathematical expression or model fitting.
  • the test object is a black chromium pattern coated on a transparent glass substrate [Fig. 8 (a)].
  • the thickness of the chromium layer is 1000 ⁇ .
  • the chrome bars are distributed equidistantly at 1 mm intervals, and the individual widths are equal to the respective intervals.
  • the numbers above indicate the number of bars in each group.
  • the distance from the target to the pinhole is 2 mm, so the magnification for pinhole imaging is 7.
  • the observation area is 1.9 mm ⁇ 1.9 mm, and is shown by the dashed square in FIG. 8 (a).
  • the separation of the central observed bar is 20 ⁇ m.
  • the spectrum can identify calcium and chromium peaks.
  • calcium is inferred from the glass substrate.
  • the glass substrate may contain other elements such as silicon and sodium. Since the fluorescent X-ray is strongly absorbed by the CCD window material, it can not be detected this time.
  • the calcium image [FIG. 8 (c)] and the chromium image [FIG. 8 (d)] are drawn by the signals in the regions (regions of interest) of each peak.
  • the XRF intensity distribution is homogeneous, indicating that calcium is uniformly distributed in the glass substrate.
  • the surface-coated chromium pattern does not leave a shadow on the calcium image, since the chromium layer is so thin as 1000 ⁇ thick.
  • the X-ray transmission of CaK ⁇ and K ⁇ is higher than 99.99%.
  • Chromium images can clearly distinguish chromium bars with a separation of 20 ⁇ m, indicating that the spatial resolution is better than 20 ⁇ m.
  • two black vertical lines are observed. They are due to linearly arranged dead pixels on the CCD sensor chip and do not affect the signals of other pixels.
  • the e2V CCD sensor used as the imaging sensor 8 in the first embodiment it takes 4 to 5 seconds to read and transfer one frame (a catalog data reading speed of 0.25 frame / second). Assuming that X-rays have been received for 4 to 5 seconds, there is a difference in effective imaging time between the reading start pixel and the reading end pixel, and the reading end pixel is the reading start pixel. As compared to the image, it took 4 to 5 seconds extra. In order to prevent this, it is necessary to stop the X-ray with a mechanical shutter. Some CCD sensors have a built-in mechanism for temporarily storing charge information to avoid such problems.
  • the temporary storage mechanism is such that an area of the same size as the area used for imaging is provided next to it, transferred at a speed on the order of microseconds (referred to as frame transfer), and then read out.
  • the area that is the transfer destination of the frame transfer is outside the opening of the X-ray window in order to shield the X-ray from hitting.
  • FIG. 9 shows a second type of object of the present invention, and is an explanatory view schematically showing a situation where charge division is the majority and isolated pixels can not be found.
  • the unit light receiving elements 155 are unit light receiving elements of the X-ray imaging element 156 and are two-dimensionally arranged.
  • the lowest threshold is also applied to one or more pixels around it (see FIG. A charge above the first threshold 115) shown in 3B is observed.
  • the charge division is remarkable, and the charge corresponding to the energy of the X-ray photon 151 is divided into a plurality of pixels.
  • the amount of charge generated by the X-ray photons 151 can be regarded as the sum of the peripheral portions of the largest unit light receiving element 152, for example, 25 pixels of 5 ⁇ 5.
  • the total charge amount for 25 pixels is equivalent to one pixel of the light receiving element 152.
  • FIG. 10 is an explanatory diagram of a method of determining what form of charge division should be numerically adopted in the situation shown in FIG.
  • the charges generated by the x-ray photons are distributed not only in one isolated pixel but in the periphery. That is, it extends around the largest unit light receiving element 161 (x, y: for example, about 5 pixels square).
  • the sum of the charge amounts in the wide area pixels gives information on the charge amount generated by the X-ray photon, that is, the energy of the X-ray.
  • the X-ray detection position should be regarded as 161 points.
  • FIG. 11 is a block diagram of a correction filtering circuit for charge division, showing a second aspect of the present invention. Note that, in FIG. 11, the same reference numerals are given to components having the same function as those of FIG. 2A, and the description will be omitted.
  • the X-ray image generation circuit 20 includes an effective event determination circuit 102, an energy level determination circuit 105, a pixel-by-pixel counter circuit 106, a charge division detection circuit 107, an X-ray counter image generation circuit 108, and an X-ray counter image generation circuit 109.
  • the charge division detection circuit 107 detects the light receiving signal read by the peripheral light receiving element 84 for the peripheral light receiving element 84 of the unit light receiving element 82 recognized as being recognized as a valid photon event by the effective event judging circuit 102; The threshold is compared to detect whether charge division occurs between the unit light receiving element 82 and the peripheral light receiving element 84.
  • the charge division processing circuit 108 calculates the median value of all light receiving elements exceeding the first threshold value for the unit light receiving element 82 and the peripheral light receiving element 84 in which charge division has occurred and detected by the charge division detection circuit 107. , Integral value and its position information are calculated.
  • the X-ray counter image generation circuit 109 generates an X-ray image by the X-ray imaging device 8 based on the count value of the pixel-by-pixel counter circuit 106.
  • the pixel-by-pixel counter circuit 106 uses the median, integral value, and position information thereof calculated by the charge division processing circuit 108. Then, counting is performed such that the count value of the pixel-by-pixel counter corresponding to all the light receiving elements exceeding the first threshold value is incremented by one.
  • FIG. 12 is an explanatory diagram of the sorting process of the pixel information in the case where the charge division is remarkable.
  • the sorting process of pixel information (S400) in the case where the charge division is remarkable corresponds to the operation of the charge division detection circuit 107 and the X-ray counter image generation circuit 108 of the X-ray image generation circuit 10.
  • the input image data exceeds the lowest threshold (first threshold 115 shown in FIG. 3B) for determining X-ray incidence / detection itself.
  • Position information of the center value, the integral value, and the pixel for which the center value is given (S402).
  • the ratio of the median value to the integral value is extracted and sorted out in a specified range (S404).
  • a histogram X-ray spectrum
  • an X-ray image showing the distribution and charge amount of pixels of specific X-ray energy are obtained (S406).
  • a histogram X-ray spectrum
  • the distribution can be given as an image for the X-ray energy or an X-ray energy range having a certain width.
  • the ratio of the median to the integral value is determined, and it is determined that the form of charge division that can be restored is effective only when it falls within a certain range.
  • 5 ⁇ 5 is empirically effective as the range of pixels for calculating the integral value, it is not limited to this, and depending on the size and structure of the pixel of the sensor, it may take a narrow range or a wider range. It may be valid.
  • a range in which the ratio of the median to the integral value is preferred is referred to as an appropriate range of the ratio
  • the lower limit of the appropriate range of the ratio is referred to as an appropriate lower limit of the ratio.
  • the upper limit of is called the appropriate upper limit of the ratio.
  • an appropriate range of pixels for which an integral value is calculated is referred to as an integration appropriate range.
  • an integration appropriate range is empirically effective when the ratio of the median to the integral is between 40% and 50%, it is not limited thereto, and it may be further different depending on the size and structure of the pixel of the sensor. Taking a range may be effective.
  • the X-ray irradiation position and the X-ray energy can be simultaneously recorded.
  • the data can be used to obtain an x-ray spectrum or an x-ray image group corresponding to any x-ray energy appearing in the spectrum.
  • the charge division detection circuit 107, the charge division processing circuit 108, and the pixel-by-pixel counter circuit 106 cooperate with the energy level determination circuit 105 in accordance with a CMOS camera in which it is difficult to find isolated pixels.
  • the operation of will be described.
  • the charge division detection circuit 107 and the charge division processing circuit 108 perform verification on all pixels at each image read out to obtain a median and an integral value, and determine that only when the ratio is within the predetermined range, it is valid. It is carried out. That is, it is interpreted that the total charge amount corresponding to the integral value is originally at the (X, Y) point giving the median value, and the reading is replaced.
  • the charge division processing circuit 108 cooperates with the energy level determination circuit 105 to determine whether or not the integral value is in a range corresponding to the specific X-ray energy.
  • the brightness of a pixel of an image is represented by an integer value between 0 and 65535. Assuming that the brightness of the dark current level is 250, it is customary to set the first threshold value to, for example, 300 arbitrary values.
  • Most X-ray signals have a level corresponding to an energy level much higher than the first threshold, and even a single photon has a value of, for example, about 800. At this time, the difference between X-ray energy and whether it was 600 or 800 is in the range of the characteristic value according to the element usually composed.
  • FIG. 13 shows an X-ray image by X-ray spectral energy using the method of correcting the influence of charge division (FIG. 9, 10) using the ratio of the median to the integral value. It uses a setup.
  • Fig. 13 shows a photograph and an X-ray fluorescence spectrum of the ceramic plate obtained by the same sCMOS camera
  • Fig. 13 (a) is a photograph of the front of the ceramic plate
  • Fig. 13 (b) is its X-ray fluorescence spectrum
  • FIG. 13 (c) is a photograph of the back side
  • FIG. 13 (d) is its spectrum. Photographs and X-ray fluorescence spectra are obtained with the same sCMOS camera. This is an embodiment using a method (FIGS. 9 and 10) of correcting the influence of charge division using the ratio of the median to the integral.
  • a visible light digital camera can measure the XRF spectrum of a sample, it can also analyze the composition of elements in the sample. For example, a ceramic plate with a white base (FIG. 13c, back) and a blue pattern (FIG. 13a, front) is tested.
  • the two photographs of FIGS. 13a and 13c are taken by the sCMOS camera before removing the optical lens system and the transparent glass cover. Since the camera does not have a color filter, the picture is monochrome.
  • the XRF spectra of the front Figure 13b
  • the back Figure 13d
  • the primary x-ray beam is a monochromatic copper K ⁇ ray.
  • the exposure time of one image is 100 ms.
  • the accumulation time is 30 minutes.
  • all peaks are identified, so the corresponding elements can be identified as well.
  • the elemental compositions of the front and back sides are very similar, it can be seen that cobalt appears only on the front and shows the relationship between cobalt and ceramic blue.
  • sCMOS cameras are used in the same way as other existing x-ray detectors.
  • concentrations of the elements in the sample can be used for quantitative XRF analysis by using experimentally obtained calibration curves or by using so-called reference free analysis based on X-ray fundamental parameters.
  • reference free analysis based on X-ray fundamental parameters can be used for quantitative XRF analysis by using experimentally obtained calibration curves or by using so-called reference free analysis based on X-ray fundamental parameters.
  • the use of sCMOS cameras is valuable in that it can be applied to the analysis of trace elements, especially when combined with total reflection geometry.
  • FIG. 14 is a configuration block diagram of a correction filtering circuit of a combined use type of an isolated pixel and charge division showing a third aspect of the present invention.
  • the X-ray image generation circuit 30 includes an effective event determination circuit 102, an isolated pixel extraction circuit 103, a noise removal circuit 104, an energy level determination circuit 105, a pixel-by-pixel counter circuit 106, a charge division detection circuit 107, and an X-ray counter image generation circuit. And an X-ray counter image generation circuit 109.
  • the correction filtering circuit for isolated pixel shown in FIG. 2A and the correction filtering circuit for charge division of FIG. Therefore, it is possible to properly cope with an event captured by an imaging sensor in which an isolated pixel and charge division coexist.
  • FIG. 15 is a block diagram of a main part of a measurement apparatus for obtaining an X-ray image by X-ray spectral energy according to a second embodiment of the present invention.
  • the X-ray analyzer according to the second embodiment of the present invention comprises a sample holder 204, a sample 205, an imaging sensor 208, and an X-ray beam guide 212.
  • the components corresponding to the X-ray tube 1, the crystal monochromator 2 and the signal control mechanical shutter 3 shown in FIG. 1 are not shown in FIG.
  • the X-ray beam guide 212 guides only the parallel component of the fluorescent X-ray as the X-ray beam 201 to the imaging sensor 208.
  • the present invention it is clear how to avoid charge division and correction of charge division, and it becomes possible to manufacture with a smaller element size than the present, and in CCD sensors and CMOS sensors designed for visible light, It enables x-ray imaging with energy discrimination. Furthermore, although the case where a CCD imaging device or a CMOS imaging device designed for visible light wavelength is used as an X-ray imaging device is shown as an embodiment of the present invention, the present invention is not limited to these. The present invention can also be applied to the case where a semiconductor pixel radiation detector specially manufactured for performing energy discrimination X-ray imaging is used as an X-ray imaging device.
  • FIG. 17 shows how to use the first threshold and the second threshold in filtering.
  • the first threshold is used to check the intensity of one light receiving element (pixel) to recognize a valid photon event.
  • the second threshold which is smaller than the first threshold, has surrounding pixels to ensure that there is no leakage current to avoid tailing on the low energy side of the peak. Used to check the Only the first threshold is used to remove the noise equivalent signal. Then, using the second threshold value, it is confirmed that the peripheral light receiving element has not caused the photon event (see S304 in FIG. 7A).
  • the functions of the components of the X-ray image generation circuit 10 or 20 can be realized by computer software. For example, in FIG. 18, after data is received by the computer via the interface and processed by software, necessary data is stored. About the principal part structure of the measuring apparatus for acquiring the X-ray image according to X-ray spectrum energy, since it is the same as that of the thing of FIG. 1, the overlapping description is abbreviate
  • the image data processed by the computer 54 is stored in a data storage 54 such as a hard disk, for example.
  • Lines 62, 64, 66, and 68 for transmitting control signals to the component parts 42, 43, and 47 of the apparatus and detection signals from the respective parts are temporarily connected to the interface 52 and can be captured by the computer 54.
  • the computer 54 sequentially fetches data from each part, performs processing, and transmits necessary instructions to each part.
  • the processed image data is stored in the data storage 54 either collectively or sequentially while performing the measurement.
  • the A / D converter can be placed on the camera side or on the camera controller.
  • the computer can receive the A / D converted image itself. More specifically, the computer can receive data of 1024 ⁇ 1024, and one pixel of 16 bits (65,536 tones).
  • the type may be diagnosed by paying attention to the shape distribution of the group.
  • Various diagnostic methods have been studied, but as a method with the highest degree of certainty and the quickest judgment, it is to use the quantitative relationship between the total charge amount and the central charge amount as an index. I understand.
  • the dispersed ones can be collected and summed.
  • the position and energy it is equivalent to the center of distribution of the dispersed charge or simply the pixel with the largest amount of charge being the detection position, and the total charge occurring there It can be considered easy.
  • this is not always the case.
  • electronic circuits for readout and the like are wired to every pixel. These are close to the places where X-rays enter and generate charges. As such, the charge distribution across pixels resulting from charge splitting may be in different patterns.
  • the ratio of the median to the sum is less than a predetermined value (specified value 1) (for example, 40%) is the occurrence of intense dissipation (something on the wall between pixels or wiring It is considered to be the cause of The amount of charge lost is so large that it is impossible to obtain information on the energy of X-rays that have already entered, and such inclusion may result in strange spectra being obtained. If the ratio between the median value and the sum value is larger than the prescribed value (for example, 50%), the total amount of charge lost is not large, but an abnormality occurs in the process of charge division It is considered to be a thing. Therefore, these are excluded.
  • a predetermined value for example, 40%
  • the specified value 1 and the specified value 2 described above are defined based on a large amount of experimental data, but the conditions on the detector side such as a CMOS element and the observation target side such as light to be detected and electromagnetic waves such as X-ray If the measurement conditions such as temperature and humidity are constant, the same value is obtained even if measurement is performed at different times and places. Therefore, the prescribed value 1 and the prescribed value 2 can be obtained for each of these conditions, and can be generalized as the standard value 1 and the standard value 2 under the conditions. From the above, it is preferable that the reference value is 1 or more such that severe dissipation does not occur. Further, it is preferable that the reference value is 2 or less so as not to cause an abnormality in the process of charge division.
  • an image obtained by correcting and converting information of X-ray energy (element) into a final state that can be easily extracted from the obtained original image is obtained.
  • the image looks like the original image you get an image that looks dark or not, so you hardly know what it is.
  • every pixel is completely isolated, and there may be differences in the brightness among the individual pixels that are bright in the dark.
  • the brightness of all the pixels of this image is checked, and those close to a specific brightness level are detected by comparison. If that is the case, at that (X, Y) point, the energy X-ray is determined to have one count, and the energy is added to the image counter.
  • the brightness of the image corresponds to energy (wavelength of light), not intensity (brightness of light).
  • the amount of charge generated corresponding to the energy of X-rays For example, for the X-ray of 6400 eV and the X-ray of 3200 eV, even if the same one photon enters, the charge generation amount differs by 2 times.
  • a threshold is generally set and used to distinguish a signal from electrical noise or the like.
  • a threshold is called a lower level discriminator (LLD).
  • LLD lower level discriminator
  • Such a threshold is also used in the embodiments of the present invention.
  • the determination that the pixel obtaining the signal is isolated is that the charge amount is above LLD for all pixels It becomes possible by checking whether it is.
  • information is acquired using only such an image that looks like only noise appears (referred to as “single photon counting mode”).
  • the charge amount corresponds to the brightness of the pixel.
  • the LLD used at this time may be a standard that is clearly brighter than noise. Since the values do not have to be strictly selected, they can be set for each of various measurements in consideration of the measurement conditions, the type of image to be obtained, and the like.
  • the newly introduced threshold is considered to be much closer to the noise level than LLD.
  • signal processing using such thresholds is not known. This is because it is determined that the signal is a clear signal distinguished from noise if it is LLD or more, and it is considered sufficient to look at only such a signal.
  • the second threshold which has been introduced, may be used to scrutinize the surrounding pixels, once for a pixel of interest once considered as isolated. Then, once the surrounding pixels that have been regarded as zero are examined again, if they exude a little next to the target pixel, even if it is a pixel that has once been regarded as isolated, it is isolated.
  • the idea is to exclude from the pixels that you This is because if the brightness of an isolated pixel is regarded as the total charge amount, and if a slight amount of bleeding is overlooked when determining the energy of the X-ray, then the energy of the X-ray may be underestimated. It is.
  • almost perfect X-rays are obtained in the case where it is recognized that one isolated pixel has a charge (in the case where the frequency of charge division is not so high, the case where the pixel size is relatively large, etc.) The effect that a spectrum is obtained can be recognized.
  • the X-ray imaging apparatus of the present invention it is possible to process X-ray photon events occurring in a CCD imaging device or a CMOS imaging device designed for visible light wavelength or a semiconductor pixel radiation detector.
  • Line energy information can be separated individually, and precise imaging pixel information can be obtained.
  • the signal processing of a digital camera detects (1) a light signal including X-rays and converts it into an electric signal (charge) and stores it (sensor), (2) it is stored for each pixel of the sensor
  • the process includes reading out the charge amount (readout circuit, A / D converter), (3) displaying the read-out as an image, or storing it as electronic data.
  • Each of these steps may be branched into further steps. There may also be detailed processes with certain features.
  • Reference Signs List 1 X-ray tube 2 crystal monochromator 3 signal control mechanical shutter 4 204 sample holder 5 205 sample 6 imaging optical system (pinhole plate) 7 X-ray camera (CCD camera, CMOS camera) 8, 208 Image sensor (two-dimensional semiconductor sensor) 82 unit light receiving element 84 peripheral light receiving element 9 cooling circuit 10, 20, 30 X-ray image generating circuit 101a, 101b X-ray beam 102 effective event judging circuit 103 isolated pixel extracting circuit 104 noise removing circuit 105 energy level judging circuit 106 pixel by pixel Counter circuit 107 Charge division detection circuit 108 Charge division processing circuit 109 X-ray counter image generation circuit 11 X-ray image display device 112, 115 First threshold 120 Second threshold

Abstract

可視光波長用に設計されたCCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器を大面積X線分光器の撮像素子として利用するのに適したX線撮像装置を提供する。分析対象物からのX線光子を含むエネルギー線を受光するX線撮像素子であって、当該X線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該X線光子の受光による電荷共有を生ずる二次元的領域を有する前記X線撮像素子と、このX線撮像素子からの受光信号を各単位受光素子毎に読み込むX線画像生成回路であって、前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第1のしきい値及びノイズと認識する基準となる第2のしきい値とを比較して、第1及び第2のしきい値を通過したものを有効なX線光子事象として、X線画像を生成する。

Description

X線撮像装置およびX線撮像素子の画像処理方法
 本発明は、CCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器を用いたX線撮像装置に関し、より詳しくは、汎用の可視光用のCCD撮像素子またはCMOS撮像素子、若しくは半導体ピクセル放射線検出器を用いた蛍光X線分析装置等のX線撮像装置に関する。
 また、本発明は、CCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器を用いたX線撮像素子の画像処理方法に関する。
 物質はさまざまな元素から構成されており、その組成によって物理的・化学的性質が大きく左右される。このため、物質の理解や新材料開発のために、含まれている元素の種類や量を分析することは重要である。物質にX線を照射した際に出てくる蛍光X線のエネルギーから元素の種類が、その強度から量がわかることが知られている(特許文献1、2)。この蛍光X線分析を行うためには、X線源とともに、専用のX線分光器やX線検出器が必要である。
 他方、X線蛍光(XRF)放射線が試料から等方的に広がる関係で、立体角が大きくなる課題がある。そこで、大きな立体角に対処することは、分光計にとって長年にわたり最も重要な要件の1つとなっている。蛍光強度が低い実験、例えばエネルギー分散型X線分光法(EDS)の微量分析では、できるだけ多くの蛍光X線を受けるには大面積X線エネルギー分散型検出器が常に好ましい。Si(Li)やSiドリフト検出器(SDD)などの通常の単一ピクセルX線分光器の面積は、100mmを超えることは珍しいことである。複雑なマルチピクセル分光器の中には、1000mmまでの非常に大きなアクティブセンサ領域を持つものがあるが、広範囲の大きさのX線分光器の代わりの選択肢を見つけることは、多様な作業環境で様々な機器をより多くの選択肢を提供するためには依然として必要である。
 電荷結合素子(CCD)カメラは、通常150mmを超える非常に大きなアクティブセンサ領域を備えている。従来のCCDカメラは、科学および産業分野を含む様々な用途で可視光画像を記録するように設計されている。X線領域の量子効率は、CCDカメラの薄い空乏層のために低いが、X線波長領域でも機能し得る。1980年代から、pnCCDのような高度なX線カメラが天文学用に開発された。空乏層の厚さは、典型的には数百μmであるので、X線領域(1~10keV)の量子効率は90%以上にできる。
 CCDカメラでは、X線光子によって生成された電荷を単純に蓄積する通常のX線撮影に加えて、このような蓄積を停止させることにより、異なるX線光子エネルギーに対応する多くの画像を得ることが可能である。単一のX線光子によって生成された電荷の数がX線のエネルギーに関する情報を与えるので、X線スペクトルおよび/またはカラフルなX線画像を得ることができる。
 しかしながら、CCDカメラで、X線光子によって生成された電荷の蓄積を停止すると、可視光で生じた信号と各種波長のX線光子による信号とが混在する為、ノイズと信号の峻別が困難になるという課題があった。
 従来の解決方法としては、CCD撮像素子またはCMOS撮像素子において、単一のX線光子によって生成された電荷が広がる領域よりも大きな画素とすることで、画像分解能を低下させて、ノイズの影響を積分することで、各種波長のX線光子による信号を抽出することが行われている(例えば、非特許文献1)。しかし、CCDカメラの画素を粗大にすると、本来単一のX線光子の有する空間領域の信号成分を捨象することに繋がり、ひいては緻密なX線分光撮像情報が得られなくなるという課題があった。
特許第3049313号公報 特許第3663439号公報
"R & D of a pixel sensor based on 0.15 μm fully depleted SOI technology", Toru Tsuboyama, et. al., Sep. 2006, Perugia, Italy, Vertex 2006, Nucl. Instr. and Meth. A. Volume 582, Issue 3, 1 December 2007, Pages 861-865. http://dx.doi.org/10.1016/j.nima.2007.07.111
 本発明は、上記の課題を解決するもので、CCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器で起こっているX線フォトン事象を処理することで、X線エネルギーの情報を個別に分離でき緻密な撮像画素情報の得られるX線撮像装置を提供することを目的としている。
 本発明の他の目的は、可視光波長用に設計されたCCD撮像素子またはCMOS撮像素子、若しくは半導体ピクセル放射線検出器を大面積X線分光器の撮像素子として利用するのに適したX線撮像装置を提供することである。
[1] 本発明のX線撮像装置は、例えば図1、図2Aに示すように、分析対象物(5)からのX線光子を含むエネルギー線を受光するX線撮像素子(8)であって、当該X線撮像素子は二次元配置の単位受光素子(82)を有すると共に、当該単位受光素子の周辺には当該X線光子の受光による電荷共有を生ずる周辺受光素子(84)を有する前記X線撮像素子と、
 このX線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該X線撮像素子内での二次元配置に応じて撮像された分析対象物のX線画像として生成するX線画像生成回路(10)と、を備え、
 前記X線画像生成回路10は、前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第1のしきい値を比較する有効事象判断回路(102)と、前記有効事象判断回路で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号およびノイズと認識する基準となる第2のしきい値を比較するノイズ除去回路(104)と、前記ノイズ除去回路で除去されないで有効な光子事象と認識された単位受光素子について、X線画像生成回路(10)の画素別カウンタ値をプラス1する計数を行う画素別カウンタ回路(106)と、この画素別カウンタ回路の計数値に基づいて、X線撮像素子8によるX線画像を生成するX線カウンタ画像生成回路(108)と、を備えることを特徴とする。
[2] 本発明のX線撮像装置においては、好ましくは、さらに、有効事象判断回路(102)で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、全ての周辺受光素子の受光信号が第1のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出回路(103)を有し、ノイズ除去回路(104)は、有効事象判断回路(102)で有効な光子事象と認識と認識された単位受光素子であって、孤立画素抽出回路(103)での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と第2のしきい値とを比較することを特徴とする。
[3] 本発明のX線撮像装置において、例えば図1、図11に示すように、分析対象物(5)からのX線光子を含むエネルギー線を受光するX線撮像素子(8)であって、当該X線撮像素子は二次元配置の単位受光素子(82)を有すると共に、当該単位受光素子の周辺には当該X線光子の受光による電荷共有を生ずる周辺受光素子(84)を有する前記X線撮像素子と、
 このX線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該X線撮像素子内での二次元配置に応じて撮像された分析対象物のX線画像として生成するX線画像生成回路(10)とを備え、
 前記X線画像生成回路10は、前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第1のしきい値を比較する有効事象判断回路(102)と、前記有効事象判断回路(102)で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出回路(107)と、電荷分割検出回路(107)で検出された電荷分割の生じた単位受光素子および周辺受光素子について、前記第1のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理回路(108)と、電荷分割処理回路(108)で演算された前記中央値、積分値およびその位置情報を用いて、前記第1のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス1する計数を行う画素別カウンタ回路(106)と、この画素別カウンタ回路の計数値に基づいて、X線撮像素子(8)によるX線画像を生成するX線カウンタ画像生成回路(108)と、を備えることを特徴とする。
[4] 本発明のX線撮像装置においては、好ましくは、X線撮像素子(8)は、CCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器であるとよい。
[5] 本発明のX線撮像装置においては、好ましくは、X線撮像素子(8)の単位受光素子(82)は、一画素として1乃至10μmの代表長さを有する二次元的領域を有するとよい。
[6] 本発明のX線撮像装置においては、好ましくは、第1のしきい値は、X線撮像素子(8)の各単位受光素子(82)で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値であるとよい。
[7] 第2のしきい値は、X線撮像素子(8)の各単位受光素子(82)で生ずるバックグランドノイズよりも大きな値であって、第1のしきい値よりも小さな値であるとよい。
[8] 本発明のX線撮像装置においては、好ましくは、周辺受光素子84は、光子事象のあった単位受光素子(82)を中心とする3×3、5×5又は7×7の何れかあるとよい。
[9] 本発明のX線撮像装置においては、好ましくは、さらに、X線画像生成回路(10)の各単位受光素子(82)毎に読み込まれた受光信号について、当該受光信号のエネルギー準位を求めるエネルギー準位判断回路(105)を有するとよい。
[10] 本発明のX線撮像装置においては、好ましくは、さらに、X線画像生成回路(10)の画素別カウンタ回路は、エネルギー準位判断回路(105)で分光されたエネルギー準位に応じて区分されているとよい。
[11] 本発明のX線撮像装置においては、好ましくは、X線画像生成回路(10)の画素別カウンタ回路は、光子事象のあった単位受光素子(82)の受光信号のエネルギー準位について、単位受光素子(82)の周辺受光素子(84)の受光信号のエネルギー準位も加算するとよい。
[12] 本発明のX線撮像装置においては、好ましくは、さらに、前記分析対象物にX線を含むエネルギー線を照射するX線照射回路を有するとよい。
[13] 本発明のX線撮像装置においては、好ましくは、さらに、X線画像生成回路(10)の画素別カウンタ値をリセットする画素別カウンタリセット回路を有するとよい。
[14] 本発明のX線撮像素子の画像処理方法は、例えば図5、図7A、図7Bに示すように、分析対象物からのX線光子を含むエネルギー線を受光するX線撮像素子の画像処理方法であって、X線撮像素子(8)は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該X線光子の受光による電荷共有を生ずる二次元的領域を有するものであって、前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第1のしきい値とを比較する有効事象判断工程(S100)と、有効事象判断工程(S100)で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第2のしきい値とを比較するノイズ除去工程(S300)と、前記ノイズ除去工程で除去されないで有効な光子事象と認識された単位受光素子について、X線画像生成回路(10)の画素別カウンタ値をプラス1する計数を行う画素別カウンタ工程とを備え、この画素別カウンタ工程での計数値に基づいて、X線撮像素子(8)によるX線画像を生成するものである。
[15] 本発明のX線撮像素子の画像処理方法においては、好ましくは、さらに、例えば図6に示すように、有効事象判断工程(S100)で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、全ての周辺受光素子の受光信号が第1のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出工程(S200)を有し、ノイズ除去工程(S300)は、有効事象判断工程(S100)で有効な光子事象と認識と認識された単位受光素子であって、孤立画素抽出工程(S200)での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第2のしきい値とを比較することを特徴とする。
[16] 本発明のX線撮像素子の画像処理方法においては、好ましくは、さらに、例えば図12に示すように、さらに、有効事象判断工程(S100)で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出工程と、前記電荷分割検出工程で検出された電荷分割の生じた単位受光素子と周辺受光素子について、前記第1のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理工程(S400)とを有し、前記画素別カウンタ工程では、前記電荷分割処理工程で演算された前記中央値、積分値およびその位置情報を用いて、前記第1のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス1する計数を行うことを特徴とする。
 本発明のX線撮像装置およびX線撮像素子の画像処理方法によれば、可視光波長用に設計されたCCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器を用いて、X線エネルギーの情報を個別に分離でき緻密な撮像画素情報が得られる。
 また、CCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器は、大面積のX線分光器に対応できるので、分光イメージングの前提条件を確立でき、CCDカメラやCMOS撮像素子を十分に活用したX線分光撮像が行える。
本発明の一実施例を示すX線スペクトルエネルギー別のX線画像を得るための測定装置の要部構成図である。 本発明の第1の態様を示す孤立画素用の補正フィルタリング回路の構成ブロック図である。 本発明の第1の態様を若干変更したもの示す孤立画素用の補正フィルタリング回路の構成ブロック図である。 検出された画素の明るさ(電荷量)からX線エネルギー(波長)を求める方法を具体的に示す図で、特に、X線撮像素子と電荷量の分布状態を示している。 検出された画素の明るさ(電荷量)からX線エネルギー(波長)を求める方法を具体的に示す図で、特に、電荷量とカウント数を示している。 X線光子により生じた電荷の複数画素への分割(一般的な電荷分割)についての説明図である。 処理対象とすべき画素の初期選別を説明する流れ図である。 本発明の対象となる第1の類型を示すもので、孤立している画素の選別(電荷分割の影響の第1次除去)を説明する流れ図である。 厳格な閾値の追加導入(電荷分割の影響の第2次除去)を説明する流れ図である。 厳格な閾値の追加導入(電荷分割の影響の第2次除去)を説明するもので、電荷量における孤立画素判定用の閾値120のエネルギー準位を示している。 孤立した画素の電荷を厳格に選定して収集する方法(図5~図7B)を用いたX線スペクトルエネルギー別のX線画像を示したものである。 本発明の対象となる第2の類型を示すもので、電荷分割が大多数を占め、孤立画素が見いだせない場合の状況を模式的に示した説明図である。 周辺画素を含めた総電荷量の算出原理の説明図である。 本発明の第2の態様を示す電荷分割用の補正フィルタリング回路の構成ブロック図である。 電荷分割が顕著な場合の画素情報の選別処理の説明図である。 中央値と積分値の比を用いて、電荷分割の影響を補正する方法(図9、10)を用いたX線スペクトルエネルギー別のX線画像を示したものである。 本発明の第3の態様を示す孤立画素と電荷分割の併用型の補正フィルタリング回路の構成ブロック図である。 本発明の第2の実施例を示すX線スペクトルエネルギー別のX線画像を得るための測定装置の要部構成図である。 X線カウンタ画像生成回路で生成されるエネルギー準位別のX線カウンタ画像の説明図である。 第1のしきい値及び第2のしきい値を使った例を示す図である。 図2の孤立画素用の補正フィルタリング回路の機能をコンピュータのソフトウェアで行う場合の、コンピュータの構成図である。 コンピュータのソフトウェアで行う場合の、流れ図である。 コンピュータのソフトウェアで行う場合の、流れ図である。
 以下、図面を用いて本発明を説明する。
 図1は、本発明の一実施例を示すX線分析装置の要部構成図である。本発明の一実施例のX線分析装置は、X線管1、結晶モノクロメータ2、信号制御メカニカルシャッター3、試料ホルダー4、ピンホール板6、X線カメラ7を有している。
 X線管1は、X線を発生させる電子管であり、フィラメント(カソード)を加熱することで発生した電子を、タングステン、モリブデン、銅などの金属(ターゲット、アノード)に衝突させることで連続X線(制動放射)と特性X線が発生するものである。
 モノクロメータ2は、単色計とも呼ばれるもので、回折格子やプリズムなどの分散素子で分散させた光の中から、ある特定の波長の光のみをスリットで取り出す装置のことである。モノクロメータは、例えば、入口スリット、コリメータ鏡、分散素子(回折格子またはプリズム)、集光鏡、出口スリットで構成されている。結晶モノクロメータ2は、例えばグラファイト(002)、2d=6.72Åが用いられる。銅ターゲットから発生する特性X線のうちKα線とKβ線とが混在していると、それぞれの特性X線に対応した異なる位置に回折ピークが現れてデータの解釈が困難になるので、通常の回折計では検出器の手前の位置に湾曲グラファイト分光器(モノクロメータ)を設置してCuKβ線による回折線を除去する。なお、モノクロメータに代えて薄いNi箔フィルターを用いても良い。
 信号制御メカニカルシャッター3は、物理的にシャッター幕を動かすものであり、実際に入射X線が試料に入射するのをブロックする必要がある。そこで、可視光用のメカニカルシャッターと同じ機構を用いながらも、X線吸収が可能な材質、厚さのものを使用するもので、例えばUNIBLITZ(登録商標)を用いることができる。
 メカニカルシャッターは、実際の撮像時間と画像の転送に要する時間の大小関係によって、その要否が決まる。後述の実施例1に示すCCDセンサでは、1フレームの読み出し転送時間が影響して、読み始めの画素と読み終わりの画素では、実効的な撮像時間に差異が生ずることになる。これを防ぐためには、メカニカルシャッターでX線を止める必要がある。
 これに対して、実施例2のCMOSセンサでは、読み出しが早いので、撮像時間によっては、メカニカルシャッターを使用せずに測定を行っている。
 試料ホルダー4は、X線分光分析の対象となる試料5を保持するもので、例えばサンプルステージを有している。
 試料5は、X線分光分析の対象となるもので、金属材料、無機材料、有機材料、生体材料等の各種の被検査対象物が含まれる。
 結像光学系6としては、ピンホール板、非球面ミラー、ゾーンプレートなどが用いられる。例えば、ピンホール板は、XRFイメージングにおいて非常に頻繁に使用される。ピンホール板6は、マイクロピンホールコリメータとも呼ばれるもので、XRFイメージング測定の場合に用いられる。なお、ピンホール板6は、非イメージングXRFスペクトル測定においては、除去される。
 X線カメラ7には、撮像センサ8が設けられている。撮像センサ8は、2次元の半導体センサで、例えばCCD、CMOS、もしくは専用に素子構造をつくりこんだX線ピクセル検出器などが用いられる。撮像センサ8は、画素数に応じた二次元配置の単位受光素子を有する。例えば、画素数が100万画素の場合には、縦x横として各1000個の単位受光素子が二次元的配置されている。例えば、撮像センサ8が10mmx10mmの場合には、各単位受光素子は10μmx10μmの大きさを有している。
 撮像センサ8としては、例えばe2V社製のCCDセンサであるCCD47-10を用いることができる。サンプル位置でのビームサイズは、約1mm(H)、5mm(V)である。
 冷却器9は、撮像センサ8をX線用途で使う場合に用いられる。冷却器9の冷却はペルチエ素子による電子冷却が一般的であるが、その排熱を空冷または水冷で除く必要がある。空冷の場合は、冷却器9としてファンが用いられる。水冷の場合は、冷却器9に循環冷却水が必要になる。図1では模式的に水冷の場合を示している。冷却温度はセンサによって-50℃~5℃と広範である。0℃以下では氷結などを避けるためセンサ周辺を真空排気することが一般的である。本実施例では、冷却CCDを―30℃、CMOSカメラを5℃で使用している。
 X線ビーム101aは、X線管1で発生したX線が、モノクロメータ2をへて、試料5に照射されるものを示している。X線ビーム101aは、試料5上で水平方向に広くなり、典型的な照射面積は10mmx5mmである。一次X線強度は10カウント/秒以上である。X線ビーム101bは、試料5上で反射したX線のうち、ピンホール板6とCCDセンサ7を通過するものを示している。
 このように構成された装置の動作を次に説明する。X線ビーム101aが低角度で試料に衝突すると、X線ビーム101aは、試料5上で水平方向に広くなる。信号制御された機械式シャッター3は、X線管1から試料5までのビーム経路に設置される。X線カメラ7が露光中の場合、機械式シャッター3は開いており、X線ビーム101aを通過させることができる。X線カメラ7が画像を読み取っているときには、機械式シャッター3を閉じてX線ビーム101aを遮断し、カメラ画像上の汚れを回避する。
 試料5上で反射したX線ビーム101bは、ピンホール板6で絞られて、X線カメラ7に至る。試料ホルダー4、ピンホール板6、CCDセンサ7の中心は、一直線に並んでいる。X線カメラ7は、光軸方向に前後に動かすことができるし、また光軸と垂直な平面方向にも動かせる位置調整機構を設けてもよい。撮像センサ8で受光された画素検出信号は、X線画像生成回路10で画像信号に変換される。
 次に、信号処理回路について説明する。
 図2Aは、本発明の第1の態様を示す孤立画素用の補正フィルタリング回路の構成ブロック図である。図において、X線撮像素子8は、分析対象物4からのX線光子を含むエネルギー線を受光するものであって、二次元配置の単位受光素子82を有する。単位受光素子82がX線光子を受光すると、当該単位受光素子82の周辺には当該X線光子の受光による電荷共有を生ずる周辺受光素子84が生ずる。周辺受光素子84は、光子事象のあった単位受光素子82を中心とする3×3、5×5又は7×7の何れかであるとよい。
 X線画像生成回路10は、X線撮像素子8からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該X線撮像素子内での二次元配置に応じて撮像された分析対象物のX線画像として生成する。
 X線画像生成回路10は、有効事象判断回路102、孤立画素抽出回路103、ノイズ除去回路104、エネルギー準位判断回路105、画素別カウンタ回路106およびX線カウンタ画像生成回路109を有する。
 X線画像生成回路10の構成要素は、各構成要素毎に、特定用途集積回路のようなハードウェア電子回路を用いても良く、またコンピュータソフトウェアで機能を実現させても良い。
 図2Bは、図2Aの本発明の第1の態様を若干変更したもの示す。ここでは、図2Aと異なり、画素別カウンタ回路106の下に、画素別カウンタリセット回路106aが加えられている。他の構成は、同じであるので、重複する説明は省略する。上述のように、第1のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス1する計数を行うが、測定すべき画素が終了した場合は、画素別カウンタリセット回路106aにより、画素別カウンタをリセットして、画素別の処理を終了させ、X線カウンタ画像生成回路109により、画像を生成する。
 有効事象判断回路102は、各単位受光素子82毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第1のしきい値とを比較する。第1のしきい値は、X線撮像素子8の各単位受光素子82で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値である。第1の閾値は、Lower Limit of Discriminationとも呼ばれ、CCDやCMOSセンサのみならず、いかなる放射線検出器の信号処理回路にも必ず使用されるものである。X線を入射しないときにも認められる検出器または信号処理回路のノイズ、暗電流の類を除去するのが目的である。従って、第1の閾値よりも信号レベルが高いときにのみ、信号と認め、その画素に見られる輝度値(画像の明るさ、X線によってつくられた電荷量に相当)の値を検討する。なお、単一光子計数モードでは、この輝度値がX線エネルギーの情報、元素分析の立場では元素の情報に対応する。
 孤立画素抽出回路103は、有効事象判断回路102で有効な光子事象と認識された単位受光素子82の周辺受光素子84について、当該周辺受光素子で読み込まれた受光信号と、第1のしきい値とを比較して、全ての周辺受光素子84の受光信号が第1のしきい値以下となっている単位受光素子82を孤立画素として抽出する。
 ノイズ除去回路104は、有効事象判断回路102で有効な光子事象と認識された単位受光素子82の周辺受光素子84について、当該周辺受光素子84で読み込まれた受光信号と、ノイズと認識する基準となる第2のしきい値とを比較する。第2のしきい値は、X線撮像素子8の各単位受光素子82で生ずるバックグランドノイズよりも大きな値であって、第1のしきい値よりも小さな値である。
 エネルギー準位判断回路105は、分光回路とも呼ばれるもので、X線画像生成回路10の各単位受光素子82毎に読み込まれた受光信号について、当該受光信号のエネルギー準位を求める。エネルギー準位は、例えば分析対象となる元素、例えばカルシウムやクロム等の金属元素に応じた値に定める。エネルギー準位は、上限値と下限値で定められる領域を有していてもよい。
 画素別カウンタ回路106は、ノイズ除去回路104で除去されないで有効な光子事象と認識された単位受光素子82について、X線画像生成回路10の画素別カウンタ値をプラス1する計数を行う。画素別カウンタ回路106は、エネルギー準位判断回路105で分光されたエネルギー準位に応じて区分されているとよく、また、光子事象のあった単位受光素子82の受光信号のエネルギー準位について、単位受光素子82の周辺受光素子84の受光信号のエネルギー準位も加算するとよい。
 X線カウンタ画像生成回路109は、この画素別カウンタ回路106の計数値に基づいて、X線撮像素子8によるX線画像を生成する。
 図16は、X線カウンタ画像生成回路で生成されるエネルギー準位別のX線カウンタ画像の説明図で、(A1)は元素A画像、(A2)は元素B画像、(A3)は元素C画像、(A4)は元素D画像、(B1)は領域Pスペクトル、(B2)は領域Qスペクトル、(B3)は領域Rスペクトルを表している。
 図16(A1)によれば、(X、Y)画素の明るさ(電荷量I)のレベルが、元素Aに相当するX線光子エネルギーの範囲(A<I<A)にあるものについて、A元素の画像の(X、Y)位置のカウントを+1する。
 図16(A2)によれば、元素Bに相当するX線光子エネルギーの範囲(B<I<B)にあるものについて、B元素の画像の(X、Y)位置のカウントを+1する。図16(A3)によれば、元素Cに相当するX線光子エネルギーの範囲(C<I<C)にあるものについて、C元素の画像の(X、Y)位置のカウントを+1する。図16(A4)によれば、元素Dに相当するX線光子エネルギーの範囲(D<I<D)にあるものについて、D元素の画像の(X、Y)位置のカウントを+1する。ここでは、光子事象が元素B画像に生じているので、図16(A2)上の座標(X、Y)位置に当該光子事象が表示される。
 図16(B1)に示す領域Pスペクトルによれば、(X、Y)領域(もしくは全領域)について、それぞれの画素の明るさ(電荷量I)のレベルによるX線光子エネルギーEがわかる。即ち、当該その(X、Y)領域内のX線スペクトルがわかるので、領域Pスペクトル領域のスペクトルのエネルギーEのカウントを+1する。
 図16(B2)に示す領域Qスペクトル、図16(B3)に示す領域Rスペクトルについても、図16(B1)に示す領域Pスペクトルと同様である。ここでは、光子事象が元素B画像に生じているので、図16(B2)上の領域Qスペクトルに当該光子事象に応じた画素の明るさ(電荷量)の箇所が+1される。
<検出された画素の明るさ(電荷量)からX線エネルギー(波長)を求める方法>
 図3A及び図3Bは、検出された画素の明るさ(電荷量)からX線エネルギー(波長)を求める方法を具体的に示す図で、図3AはX線撮像素子と電荷量の分布状態、図3Bは電荷量とカウント数を示している。図3Aにおいて、単位受光素子110は、X線撮像素子の単位受光素子で、二次元に配置されている。第1のしきい値112は、単位受光素子110で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値である。ノイズ相当信号111は、単位受光素子110で検出された画素の明るさ(電荷量)が、第1のしきい値112よりも低いものである。有効光子事象信号113は、単位受光素子110で検出された画素の明るさ(電荷量)が、第1のしきい値112よりも高いものである。
 図3Bでは、横軸に電荷量、縦軸にカウント数を示している。第1のしきい値115は、X線を照射してもしなくても見られる暗電流のレベルの電荷量117よりも十分大きな閾値である。有効光子事象信号116は、第1のしきい値115よりも大きな電荷量が得られたものであり、これを集計することで、X線により生じた電荷の情報を得ることができる。暗電流のレベルの電荷量は、暗電流の中央値118が存在し、事実上の暗電流の上限値に対してノイズ幅119を有している。
<単一光子計数条件での撮像>
 続いて、単一光子計数条件での撮像について説明する。
 単一光子計数条件での撮像は、換言すれば、非常に短い撮像時間での撮像を繰り返し、継続的に行う撮像ということである。撮像時間を長くすると、蓄積された画像は鮮明になるが、その意味するところは、1つの画素に多数の電荷が蓄積されているということである。ところが、撮像時間を長くしたのでは、どんなX線光子エネルギーのX線によって、その電荷が作られたのかという、個別X線光子の情報はもはやわからない。
 これに対して、1つの画素、もしくはその周辺に1つのX線光子しか来ていないのであれば、その電荷量から逆算して、X線エネルギーを求めることができる。例えば、仮に8000eVのX線と6000eVのX線が両方来ている場合、それぞれのX線光子が作り出す電荷量は4:3である。そのため、お互いを足し合わせたり、同じものを2回足し合わせたりしない限り、画素の明るさの情報から、どっちのX線エネルギーであるかを判別することができる。
<X線光子により生じた電荷の複数画素への分割>
 図4は、X線光子により生じた電荷の複数画素への分割(一般的な電荷分割、チャージシェアリング)についての説明図である。
 通常の放射線検出器(半導体検出器)では、もともと1素子であるので、図3A及び図3Bでの説明の通り、1個のX線光子121が来るたびに、その総電荷量125を短時間で収集すれば、エネルギー分析122が測定できる。しかし、X線撮像素子124は多数の単位画像素子123から成り立っているので、1個のX線光子121によって発生した電荷が必ずしも1つの素子にいくとは限らない。
 すると、X線撮像素子で、X線光子によって生成された電荷の蓄積を停止すると、可視光で生じた信号と各種波長のX線光子による信号とが混在する為、ノイズと信号の峻別が困難になるという課題があった。
<電荷分割が極端に多くはなく、孤立した画素による電荷収集が十分行える場合>
 図5は、処理対象とすべき画素の初期選別を説明する流れ図である。処理対象とすべき画素の初期選別(S100)は、X線画像生成回路10の有効事象判断回路102の動作に対応している。
 処理対象とすべき画素の初期選別(S100)では、入力されるデータは、X線撮像素子の全画素(単位受光素子)に記録されている値である(S102)。そして、第1のしきい値(最低閾値)を超えている画素(単位受光素子82)を選別して集計する(S104)。出力データとしては、電荷量のヒストグラム(X線スペクトル)となる(S106)。
 図6は、本発明の対象となる第1の類型を示すもので、孤立している画素の選別(電荷分割の影響の第1次除去)を説明する流れ図である。孤立している画素の選別(S200)は、X線画像生成回路10の孤立画素抽出回路103の動作に対応している。
 孤立している画素の選別(S200)では、入力される画像データは、X線がないときでもみられるバックグラウンドを差し引いたX線画像である(S202)。そして、孤立している画素の選別を行う(S204)。即ち、有効な信号が認められた画素(単位受光素子82)の周辺画素(周辺受光素子84)を点検し、どの周辺画素も最低閾値以下であることを確認する。出力データとしては、孤立している画素の位置と電荷量を示す、電荷量のヒストグラム(X線スペクトル)となる(S206)。
 即ち、その画素が孤立している画素であるかどうかを判定し、孤立しているものだけを選ぶということである。このようにして、電荷量のヒストグラムからX線スペクトルにおける電荷の複数画素への分割の影響を取り除く。
 図6に示す方法によって、電荷分割の影響はひとまず除外することができるが、図6のS206に示すように、得られたスペクトルには低エネルギー側に裾を引く現象132がよく見られる。これは、わずかながら周辺画素への散逸があること等、1画素への電荷収集が不完全であることに対応している。
 図7A及び図7Bは、厳格な閾値の追加導入(電荷分割の影響の第2次除去)を説明する図で、図7Aは流れ図、図7Bは電荷量における孤立画素判定用の閾値120のエネルギー準位を示している。なお、厳格な閾値の追加導入(S300)は、X線画像生成回路10のノイズ除去回路104の動作に対応している。
 厳格な閾値の追加導入(S300)では、入力される画像データは、孤立している画素の位置と電荷量である(S302)。周辺部分の画素に対しては、通常採用される最低閾値よりも小さな値の孤立画素判定用の閾値(図7Bに示す第2のしきい値120)を適用し、周辺画素のうち1つでも、これを超えたものがあれば、電荷収集が不完全とみなして、除外する(S304)。出力データとしては、孤立している画素の位置と電荷量を示す、電荷量のヒストグラム(X線スペクトル)であって、電荷分割の影響で電荷収集が不完全なノイズと判断されるものが除去されたデータが得られる(S306)。
 上記の図5から図7A及び図7Bの信号処理を採用することで、2次元検出器を用いて、X線の照射位置とX線エネルギーを同時に記録することができ、そのデータを用いて、X線スペクトルを得ることも、そのスペクトルに現れるどのX線エネルギーに対応するX線画像群も得ることもできるようになった。
 本発明の第1の類型では、有効事象判断回路102と孤立画素抽出回路103によって、まず孤立画素を探して、判定している。そこで、孤立画素の判定・探索の時点では、見つかった孤立画素は、第1のしきい値よりも大きな輝度値を持ち、孤立画素の周辺の8画素は、第1のしきい値よりも小さな輝度値になっているのは、孤立画素の定義から当然のことである。これに加え、周辺受光素子84が、光子事象のあった単位受光素子82を中心とする3×3とする場合に、ノイズ除去回路104によって、周辺の8画素に対して、いっそう厳しい制限を課しているのが、本発明の第1の類型の特徴部分である。周辺受光素子84は、第1のしきい値よりも小さいだけでなく、もっと厳しい第2のしきい値より小さい場合に、有効な周辺8画素を有する孤立画素として採用し、該当しない場合には、第2の類型の電荷分割で対処することになる。
 この場合、画素別カウンタ回路106の動作としては、第1の態様として、エネルギー準位判断回路105と協働しない類型がある。第1の態様では、(X,Y)の画素がカウンタになっていて、特定地点について1追加して、その地点が明るくなる。この場合の画素別カウンタは、第1のしきい値を超えてさえいればいいという立場であるので、X線のエネルギーを無視しており、すべてのX線スペクトルをすべて取り入れた画像を与える。
 次に、画素別カウンタ回路106の動作における第2の態様として、エネルギー準位判断回路105と画素別カウンタ回路106が協働して動作する類型について説明する。
 上述したように、本発明の第1の類型では、第1のしきい値を超えており、かつ、周辺の画素については、どの8画素も第2のしきい値よりも下回っており、厳格な孤立画素であることを確認したうえで、エネルギー準位判断回路105を用いて、その画素の輝度値が特定のX線エネルギーに対応する範囲にあるか判定する。例えば、16bitのA/D変換器を搭載しているCCDカメラなどでは、画像の画素の明るさは、0から65535の間の整数値で表現されいる。仮に暗電流のレベルの明るさが250だったとすると、第1のしきい値は、例えば300台の任意の値にセットすることが通例である。また、第2のしきい値は例えば250~300の間に設定される。ほとんどのX線の信号は、第1のしきい値よりも大幅に大きなエネルギー準位に相当するレベルであり、単一光子でも例えば800程度の値になる。このときに600だったのか、800だったのか、がX線エネルギーの違いであり、通常構成される元素に応じた特有の値の範囲にある。そこで、エネルギー準位判断回路105を用いて、特定の元素の画像を作成する場合には、例えば750以上820以下みたいな設定をし、もし地点(X,Y)で、その輝度値がそうだったときは、その元素の画像の(X,Y)地点をプラス1することになる。つまり、計測している際には、画像は1つの画像に過ぎないが、カウンタにあたる画像は、注目したい元素の数だけ、注目したX線エネルギーの種類だけ、独立に用意し、使用することができる。
(実施例1 孤立した画素の電荷を厳格に選定して収集する方法)
 図8は、孤立した画素の電荷を厳格に選定して収集する方法(図5~図7B)を用いたX線スペクトルエネルギー別のX線画像を示したもので、図1のセットアップを用いている。図8は20μmの空間分解能で得られたピンホールXRFイメージングで、図8(a)は検査対象の写真を示し、観測領域は破線の四角で表示される。図8(b)は、得られたXRFスペクトルを示す。基板内のCaと黒色パターンのCrの結像結果をそれぞれ図8(c)、図8(d)に示す。20μm間隔のCrバーははっきりと区別される。
 図1のセットアップにおいて、X線管1としては、密閉式X線管であって、銅ターゲット(TOSHIBA A26L-Cu、1.5kW)を用いている。銅ターゲットから発生する特性X線には、波長約1.54ÅのKα線と波長1.38ÅのKβ線が含まれる。Kα線はさらにピーク波長1.5406ÅのKα1線と1.5444ÅのKα2線が強度比2:1で混合したものであるが、二つのピークは近接した位置にあるので、標準的な粉末回折計ではこれらを区別せずに両方とも検出する仕様となっている。
 ピンホール板6として、UVレーザーによって厚さ20μm~50μmのタングステン箔に5μm~100μm径の穴をあけたものを用いている。20μm以上の空間分解能を有することは、例えば標準的な分解能目標をテストすることによって直接可視化できる。ピンホールは、低コスト、容易な調整、およびエネルギー依存性などの利点を有するので、ピンホールによる高空間分解能を得ることは、全視野XRF撮像の応用を拡大するために重要である。X線カメラ7は、全視野XRF撮像用のマイクロピンホール板6のコリメータと結合され、これにより20μmより良好な空間分解能が達成される。この結果は、数学的表現やモデルフィッティングなしで直接視覚化できる。
 試験対象は、透明なガラス基板上にコーティングされた黒色クロムパターンである[図8(a)]。クロム層の厚さは1000Åである。写真に示されているように、各グループにおいて、クロムバーは1mmの間隔で等距離に分布し、個々の幅はそれぞれの間隔に等しい。上の数字は各グループのバーの数を示している。この実施例では、ターゲットからピンホールまでの距離は2mmであり、従って、ピンホール撮像の倍率は7である。観測領域は1.9mm×1.9mmであり、図8(a)の破線の四角で示されている。中央の観察されたバーの分離は20μmである。
 ピンホールを通過したX線は弱いので、十分な統計を得るには長い蓄積時間が必要である。この実験では、蓄積時間が24時間に達すると、元素イメージングにおけるクロムパターンの輪郭が鮮明になり、区別できる。
 空間分解能の実用的な限界を見るために、高品質の要素画像を得るために蓄積時間が1週間に延長されている。CCDカメラ1枚の画像を読み取るのに4秒必要なので、合計測定時間は蓄積時間よりも20%長くなる。測定と並行して、画像は取得後直ちに実験室で用意されたPythonコードで処理される。XRFスペクトルおよび要素画像を含むすべてのX線データは、測定が終了した直後に測定および出力できる。XRFスペクトルを図8(b)に示す。スペクトルでは、カルシウムとクロムのピークを特定できる。ここで、カルシウムはガラス基板から推測される。ガラス基板は、シリコンおよびナトリウムなどの他の元素を含んでいてもよい。蛍光X線がCCD窓材料に強く吸収されるため、今回は検出できていない。
 カルシウムの画像[図8(c)]と、クロムの画像[図8(d)]はそれぞれのピークのROI(regions of interest)内の信号によって描かれる。カルシウム像では、XRF強度分布は均質であり、カルシウムがガラス基板中に均一に分布していることを示している。クロム層は1000Åの厚さと非常に薄いので、表面被覆されたクロムパターンはカルシウム像上に影を残さないことが分かる。CaKαとKβのX線透過率は99.99%より高い。クロム画像では、20μmの分離を有するクロムバーを明確に区別することができ、空間分解能が20μmより良好であることを示している。カルシウムとクロムの画像では、2本の黒い垂直線が観察される。それらは、CCDセンサチップ上の直線状に配置されたデッドピクセルに起因しているものであり、他のピクセルの信号に影響を与えない。
 なお、実施例1で撮像センサ8として使用するe2V社製CCDセンサでは、1フレームの読み出し転送にも4~5秒かかる(カタログデータで、0.25フレーム/秒の読み出し速度)。その4~5秒の間にもX線が来ているとすれば、読み始めの画素と読み終わりの画素では、実効的な撮像時間に差異があり、読み終わりの画素は読み始めの画素に比べ4~5秒余分に撮像したことになる。これを防ぐためには、メカニカルシャッターでX線を止める必要がある。CCDセンサのなかには、こうした問題を避けるために、別途電荷情報を一時蓄える機構を内蔵しているものがある。当該一時蓄える機構は、撮像に用いるエリアと同サイズのエリアを隣に設け、そこにマイクロ秒オーダーの速度で転送し(フレームトランスファーと呼ぶ)、その後に読み出すものである。フレームトランスファーの転送先にあたるエリアは、X線が当たらないように遮光するため、X線窓の開口の外にある。
<電荷分割が大多数を占め、孤立画素が見いだせない場合>
 上記の議論は、電荷分割が極端に多くはなく、孤立した画素による電荷収集が十分行える場合の説明である。しかし、センサの1画素あたりのサイズが小さくなり、またセンサの構造によっては、孤立した画素による電荷収集が困難になる。
 図9は、本発明の対象となる第2の類型を示すもので、電荷分割が大多数を占め、孤立画素が見いだせない場合の状況を模式的に示した説明図である。
 図9中、単位受光素子155は、X線撮像素子156の単位受光素子で、二次元に配置されている。第2の類型の場合には、1個のX線光子151が来ると、その検出電荷量が最大の単位受光素子152であっても、その周囲の1つ以上の画素にも最低閾値(図3Bに示す第1のしきい値115)を超える電荷が認められる。電荷分割が顕著であり、X線光子151のエネルギーに相当する電荷は複数画素に分割されている。X線光子151によって生じた電荷量は最大の単位受光素子152の周辺部分、例えば5×5の25画素の総和と見なすことができる。すなわち、その25画素分の総電荷量が受光素子152の1画素に来たのと等価である。現実には、電荷分割にも多くの形態があり、X線光子のエネルギーが散逸して電荷量の情報として復元できないものもある。
 図10は、図9に示す状況のなかで、どのような電荷分割の形態のみを数値的に採用すべきかを判定する方法の説明図である。第2の類型の場合には、X線光子によって生じた電荷は1つの孤立した画素にとどまらず、周辺に広がって分布している。すなわち、最大の単位受光素子161を中心として、その周囲(x、y:例えば、5画素四方程度)に広がっている。この広域の画素にある電荷量の総和が、X線光子によって生じた電荷量、すなわちX線のエネルギーの情報を与えるものである。また、X線の検出位置は161の地点と見なすべきである。電荷分布を数学関数でパラメータフィッティングを行い、その重心位置をもって位置とし、その関数の与える積分値をもって総電荷量とすることも可能であるが、その重心位置は、まず161の画素内にはいると考えられる。
 実際には、上記のような単純な総和をもって、電荷分割前の情報を復元できないケースは多く存在する。電荷分割と一口に言っても多くの形態があり、X線光子のエネルギーが散逸し、どの画素にも電荷量として記録されない場合もある。そこで、電荷分割の状態を識別するため、中央値と積分値という2つのパラメータを導入する。中央値とは、その周辺で最大の電荷量を持つ画素の電荷量である。積分値とは、中央値を示す画素の周辺の特定の範囲の画素に認められる電荷量の総和である。この両者の比に着目すると、復元可能な電荷分割の形態を判別することができる。
 図11は、本発明の第2の態様を示す電荷分割用の補正フィルタリング回路の構成ブロック図である。なお、図11において、前記図2Aと同一作用をするものには同一符号を付して、説明を省略する。
 X線画像生成回路20は、有効事象判断回路102、エネルギー準位判断回路105、画素別カウンタ回路106、電荷分割検出回路107、X線カウンタ画像生成回路108およびX線カウンタ画像生成回路109を有する。
 電荷分割検出回路107は、有効事象判断回路102で有効な光子事象と認識と認識された単位受光素子82の周辺受光素子84について、当該周辺受光素子84で読み込まれた受光信号と、第1のしきい値とを比較して、単位受光素子82と周辺受光素子84との間で電荷分割が生じているか検出する。
 電荷分割処理回路108は、電荷分割検出回路107で検出された電荷分割の生じた単位受光素子82と周辺受光素子84について、第1のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する。
 X線カウンタ画像生成回路109は、この画素別カウンタ回路106の計数値に基づいて、X線撮像素子8によるX線画像を生成する。単位受光素子82と周辺受光素子84との間で電荷分割が生じている場合には、画素別カウンタ回路106は、電荷分割処理回路108で演算された中央値、積分値およびその位置情報を用いて、第1のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス1する計数を行う。
 図12は、電荷分割が顕著な場合の画素情報の選別処理の説明図である。なお、電荷分割が顕著な場合の画素情報の選別処理(S400)は、X線画像生成回路10の電荷分割検出回路107、X線カウンタ画像生成回路108の動作に対応している。
 電荷分割が顕著な場合の画素情報の選別処理(S400)では、入力される画像データは、X線入射・検出そのものを判定する最低閾値(図3Bに示す第1のしきい値115)を超えている全画素の中央値、積分値、および、中央値を与える画素の位置情報である(S402)。そして、中央値と積分値の比が、指定範囲に入っているものを抽出・選別(S404)。出力データとしては、特定のX線エネルギーの画素の分布と電荷量を示す、電荷量のヒストグラム(X線スペクトル)および、X線画像が得られる(S406)。ヒストグラム(X線スペクトル)は、全画素を足し合わせたものに対しても、任意の画素範囲内についても与えることができる。また、そのようなスペクトルに現れる特定のピークに注目して、そのX線エネルギーまたは、ある程度の幅をもったX線エネルギー範囲について、その分布を画像として与えることができる。
 本発明の第2の類型においては、中央値と積分値の比を求め、それがある範囲に入っている場合のみ、復元可能な電荷分割の形態であるとして、有効と判断する。積分値を算出する画素の範囲として、経験的に5×5が有効であるが、これに限定されるものではなく、センサの画素のサイズや、構造によっては狭い範囲もしくはさらに広範囲をとることが有効である場合もある。ここで、中央値と積分値の比が入ると好ましいとされる範囲を、ここでは、比の適正範囲といい、比の適正範囲の下限値を比の適正下限値といい、比の適正範囲の上限値を比の適正上限値という。また、積分値を算出する画素の適正な範囲を積分適正範囲という。また、中央値と積分値の比が40%と50%の間のとき、経験的に有効であるが、これに限定されるものではなく、センサの画素のサイズや、構造によってはさらに別の範囲をとることが有効である場合もある。
 図12に示す中央値を積分値の比をとるという方法によれば、撮像センサ(二次元半導体センサ)を用いたカメラでも、X線の照射位置とX線エネルギーを同時に記録することができ、そのデータを用いて、X線スペクトルを得ることも、そのスペクトルに現れるどのX線エネルギーに対応するX線画像群も得ることもできる。
 この場合、孤立画素を見つけることが困難なCMOSカメラの場合に即して、電荷分割検出回路107と電荷分割処理回路108、および画素別カウンタ回路106がエネルギー準位判断回路105と協働する類型の動作について説明する。
 電荷分割検出回路107と電荷分割処理回路108によって、画像読み出しのつど、全画素についての検証を行い、中央値と積分値を求め、その比が既定の範囲にはいった場合のみを有効とする判定を行っている。つまり、積分値に対応する総電荷量が、中央値を与える(X,Y)地点にもともとあったのだと解釈して、読み替えを行う。読み替え後は、中央値のあった地点の輝度が積分値になり、周辺の5×5の画素の値はゼロにセットされる。同時に、電荷分割処理回路108は、エネルギー準位判断回路105と協働して、その積分値が、特定のX線エネルギーに対応する範囲にあるかどうかも判定する。例えば、16bitのA/D変換器を搭載しているCMOSカメラなどでは、画像の画素の明るさは、0から65535の間の整数値で表現される。仮に暗電流のレベルの明るさが250だったとすると、第1のしきい値は、例えば300台の任意の値にセットすることが通例である。ほとんどのX線の信号は、第1のしきい値よりも大幅に大きなエネルギー準位に相当するレベルであり、単一光子でも例えば800程度の値になる。このときに600だったのか、800だったのか、がX線エネルギーの違いであり、通常構成される元素に応じた特有の値の範囲にある。
 しかし、CMOSカメラでは、電荷分割があまりに頻繁に生じているので、孤立画素は見つけられない。5×5のグループで中央値と積分値をもとめ、その比が40~50%だった場合、有効と判定することにし、中央値のあった画素の地点に積分値を与え、周囲はゼロに置き換える。もし、特定の元素の画像を作ろうとするのであれば、例えば750以上820以下のような設定をし、もし地点(X,Y)で、その新しい輝度値である積分値がそうだったときは、その元素の画像の(X,Y)地点をプラス1することになる。つまり、計測している際には、画像は1つの画像に過ぎないが、カウンタにあたる画像は、注目したい元素の数だけ、注目したX線エネルギーの種類だけ、独立に用意し、使用することができる。
(実施例2 特徴的な化学元素への応用)
 図13は、中央値と積分値の比を用いて、電荷分割の影響を補正する方法(図9、10)を用いたX線スペクトルエネルギー別のX線画像を示したもので、図1のセットアップを用いている。図13は、同じsCMOSカメラで得られたセラミックプレートの写真とX線蛍光スペクトルで、図13(a)はセラミックプレートの前面の写真であり、図13(b)はそのX線蛍光スペクトルである。図13(c)は裏側の写真、図13(d)はそのスペクトルである。写真とX線蛍光スペクトルは同じsCMOSカメラで得られる。
これは、中央値と積分値の比を用いて、電荷分割の影響を補正する方法(図9、10)を用いている実施例になる。
 可視光デジタルカメラは、試料のXRFスペクトルを測定できるので、試料中の元素の組成も分析できる。例えば、白色ベース(図13c、裏面)および青色パターン(図13a、前面)を有するセラミックプレートを試験する。
 図13aおよび図13cの2つの写真は、光学レンズシステムおよび透明ガラスカバーを取り外す前にsCMOSカメラによって撮影される。カメラにカラーフィルタがないため、写真はモノクロである。レンズシステムを取り外し、ガラスカバーをX線窓で置き換えた後、前面(図13b)と背面(図13d)のXRFスペクトルをそれぞれ同じsCMOSカメラで測定する。実験では、一次X線ビームは単色の銅Kα線である。1枚の画像の露光時間は100msである。蓄積時間は30分である。
 スペクトルでは、すべてのピークが同定され、したがって対応する要素も同様に同定することができる。表側と裏側の元素組成は非常に似ているが、コバルトは正面のみに現れ、コバルトとセラミックの青色の関係を示していることが分かる。
 XRF分析では、sCMOSカメラは他の既存のX線検出器と同じように使用される。したがって、実験的に得られた較正曲線を用いるか、X線基本パラメータに基づいていわゆる基準フリー分析を使用することによって、サンプル中の元素の濃度を定量的XRF分析に使用可能である。また、sCMOSカメラの使用は、特に全反射ジオメトリと組み合わせた場合の微量元素の分析に適用できる点に価値がある。
 図14は、本発明の第3の態様を示す孤立画素と電荷分割の併用型の補正フィルタリング回路の構成ブロック図である。なお、図14において、前記図2A、図11と同一作用をするものには同一符号を付して、説明を省略する。
 X線画像生成回路30は、有効事象判断回路102、孤立画素抽出回路103、ノイズ除去回路104、エネルギー準位判断回路105、画素別カウンタ回路106、電荷分割検出回路107、X線カウンタ画像生成回路108およびX線カウンタ画像生成回路109を有する。
 図14に示す孤立画素と電荷分割の併用型の補正フィルタリング回路の構成ブロック図によれば、図2Aに示す孤立画素用の補正フィルタリング回路と、図11の電荷分割用の補正フィルタリング回路を統合してあるため、孤立画素と電荷分割が併存する様な撮像センサの撮影した事象にも的確に対応できる。
 図15は、本発明の第2の実施例を示すX線スペクトルエネルギー別のX線画像を得るための測定装置の要部構成図である。
 図において、本発明の第2の実施例のX線分析装置は、試料ホルダー204、試料205、撮像センサ208、X線ビームガイド212で構成されている。なお、図1に示したX線管1、結晶モノクロメータ2、信号制御メカニカルシャッター3に相当する構成要素に関しては、図15では図示を省略している。
 図15の実施例では、X線ビームガイド212によって、X線ビーム201として蛍光X線の平行成分のみを撮像センサ208に導いている。
 本発明によれば、電荷分割の避け方や電荷分割の補正が明らかとなり、現在よりももっと小さな素子サイズでの製造も可能になり、可視光用に設計されているCCDセンサ、CMOSセンサにおいて、エネルギー識別のできるX線イメージングを可能となる。
 さらに、本発明の実施例としては、可視光波長用に設計されたCCD撮像素子またはCMOS撮像素子をX線撮像素子として使用する場合を示したが、本発明はこれらに限定されるものではなく、エネルギー識別のできるX線イメージングを行なう目的で専用に製作されている半導体ピクセル放射線検出器をX線撮像素子として使用する場合にも適用できる。
 図17は、如何に第1のしきい値及び第2のしきい値をフィルタリングにおいて使用するかを示す。第1のしきい値は、有効な光子事象を認識するために1つの受光素子(pixel)の強度をチェックするために使用される。第1のしきい値より小さい第2のしきい値は、ピークの低エネルギー側に裾を引く現象を避けるために、リーク電流がないことを確かにするためにその周辺受光素子(surrounding pixels)をチェックするために使用される。第1のしきい値のみを使用して、ノイズ相当信号を取り除く。そして、第2のしきい値を使用して、周辺受光素子が光子事象を起こしていないことを確認する(図7AのS304参照)。
 X線画像生成回路10又は20の構成要素の機能を、コンピュータソフトウェアで実現させることができる。例えば、図18において、インターフェースを介してデータをコンピュータが受信し、ソフトウェアで処理した後に、必要なデータを保存する。X線スペクトルエネルギー別のX線画像を得るための測定装置の要部構成については、図1のものと同一であるので、重複する説明は省略する。X線管1を含むX線発生部42、信号制御メカニカルシャッター3、試料ホルダー4やピンホール板6を含む試料部44、及びX線カメラ7を含む撮像部46から、測定装置の要部は構成される。これら各部とは、インターフェース52を介してコンピュータ54に接続される。また、コンピュータ54により処理された画像データは、例えばハードディスクのようなデータストーレージ54に保存される。本装置の各構成部42、3、44、及び47への制御信号や各部からの検出信号等を伝達するライン62、64、66、及び68は、インターフェース52に一旦接続され、コンピュータ54により取り込める形の信号に変えた後、コンピュータ54にライン70で接続されるか、或いは、コンピュータからの制御命令を各部が制御信号として受信可能な信号に変えた後、各部に制御信号として送信される。コンピュータ54では、逐次的にデータを各部から取り込み、処理を行い、必要な命令を各部に送信する。そして、処理後の画像データは、一括して、或いは、測定を行いながら、逐次的に、データストーレージ54に保存される。より具体的には、例えば、A/D変換器をカメラ側若しくはカメラコントローラーに置くことができる。コンピュータは、A/D変換された画像そのものを受け取ることができる。より具体的な例示をすれば、1024×1024で、1画素が16bit (65,536諧調)のデータを転送してコンピュータが受け取ることができる。受け取ったデータの処理の内容には、個々の工程がある。例えば、まずは、全画素を1つずつ、どんな明るさになっているかを点検する。明るさがまったくゼロということはなく、何も映っていない場合でもダークカレント(電子ノイズ)に対応する明るさがある。そのダークレベルと比較してどの程度であるかを見渡せば、明らかに明るい画素は区別でき、その画素の周囲に明るいものがあるかないか(孤立して明るい画素であるかどうか)もわかる。具体的な例として、図19及び20において流れ図で示す。
(1)孤立した1個のピクセルに電荷があると認められるケース
(電荷分割の頻度がそれほど多くはないケース、比較的ピクセルサイズが大きい場合等)
まず、対象となるピクセルが孤立しているかを検査する(S500)。孤立していることを確認した(S502)後に、そのピクセルの周囲の8つのピクセルを点検する(S504)。このようなピクセルは、たぶんゼロだと普通は見なすところ、通常よりも厳しい基準で(新たな閾値を使って)点検する(S506)。新閾値は、通常よりも厳しいので、厳密な条件を満たすものだけが採用され(S506でYes)、そのデータが保存される(S508)。もし、厳密な条件を満たさなければ、そのデータは保存対象から外れる(S506でNo)。そして、全てのピクセルをチェックしたかどうかを確認し(S508)、そうであれば、このルーチンを終了させる。全てのピクセルをチェックしていない場合は(S510でNo)、再び、別のピクセルについて、孤立していることを確認し(S502)する。このようなケースは、CCDカメラのケースに相当することができる。上述のようにして、孤立して明るい画素を見出すことができる。孤立しておらず、隣同士も明るいといったものがあったとしても、それは採用しないようにする。このようにして、見つかった明るい孤立画素を更に厳しく選定してもよい。周囲の8画素が普通に暗いくらいでは容認せず、少しでも余分な電荷のしみだしがあった場合を除外するための厳選を行うことができる。
(2)画像を見渡してどのピクセルも孤立しておらず、その周囲のピクセルに電荷が分散しているケース(S600)(電荷分割の頻度が多いケース、ピクセルサイズが小さい場合、特にS¥CMOS素子の場合など)においては、分散している電荷の総和をとり(S602)、その値と、中央の電荷量の比をとり(S604)、規定の範囲におさまっているものだけを採用する(これは電荷の分布の仕方をみきわめることに対応し、散逸の激しいものをあらかじめ除外する意味がある)。規定の範囲におさまっているのは、その比(R)が、規定1及び規定2の間に入っているかどうかをチェックすることで行う(S606)。入っている場合(S606でYes)は、そのデータを保存する(S608)。入っていない場合(S606でNo)は、そのデータは保存されない。このチェック(S606)後には、全ピクセルをチェックしたかどうかを検査する(S610)。全てをチェックした場合(S610でYes)は、このルーチンを終了させる。全てをチェックしていない場合(S610でNo)は、別の対象ピクセルについて、分散している電荷の総和をとる(S602)。このようにして、対象となるピクセルが規定範囲に収まっているかを順次調べていく。このようなケースは、CMOSカメラのケースに相当する。CMOSカメラで、電荷分割問題を解決し、X線画像を得ることに成功している例は、本発明者らが知る限りにおいて、他になく、本発明者らが唯一成功している。このような場合においては、画素サイズが小さく、電荷分割は極めて高い頻度で生じる。その集合、散逸状態を画像の上でグループとして見出すことは可能である。そのグループの形状分布に着目して、タイプを診断してもよい。その診断の方法は種々検討されてきたが、1番確実性が高く、かつ最も迅速に判断できる方法として、総和の電荷量と中央の電荷量の量的な関係を指標にすることであることが分かった。
ここで、電荷が複数のピクセルに分散してしまう現象を補正する対策として、分散したものをかき集めて総和をとることができる。位置とエネルギー(電荷の総量)の両方を得るために、分散した電荷の分布の重心若しくは単純に最大の電荷量のあったピクセルを検出位置とし、そこに全電荷が生じたのと等価であるとみなせれば容易である。しかしながら、そうなるとは限らない。例えば、CMOS素子では、読み出しなどのための電子回路が、どのピクセルにも配線される。これらは、X線が入って電荷を発生している場所に近い。そのため、電荷分割の結果として生じる複数ピクセルにまたがる電荷分布はさまざまなパターンになる場合がある。中央値と総和値の比をとって所定の値(規定値1)(例えば、40%)より小さいのは、激しい散逸が生じたものである(ピクセルの間の壁でなにか生じたか、あるいは配線の影響を受けたかが原因と考えられる)。失われた電荷量が多すぎて、すでに入射したX線のエネルギーの情報を得ることはできないもので、このようなものを含めると奇妙なスペクトルが得られる結果を招くおそれがある。中央値と総和値の比をとって所定の値(規定値2)(例えば、50%)より大きいものは、失った電荷の総量はあまり多くはないものの、電荷分割の過程で異常を生じたものと考えられる。そのため、これらは除外される。上述する規定値1及び規定値2は、多くの実験データに基づいて規定されているが、CMOS素子等の検出器側の条件及び検出対象となる光やX線等の電磁波等の観測対象側の条件並びに温度や湿度等の測定条件を一定にすれば、時と場所を変えて測定を行っても同じ値となる。そのため、これらの条件毎に、規定値1及び規定値2を求め、その条件における基準値1及び基準値2として、一般化することができる。以上より、激しい散逸が生じることのないような基準値1以上であることが好ましい。また、電荷分割の過程で異常を生じたものとならないような基準値2以下であることが好ましい。
 上述するような処理を行うと、最終的に、取得した元の画像に対し、X線エネルギー(元素)の情報を抽出しやすい最終的な状態に補正、変換した画像が得られる。この時点では、画像を一見しても元の画像と同じく、何が写っているかほとんどわからないような暗い、見えるかどうかの画像が得られる。しかしながら、補正された後では、どの画素も完全に孤立しており、暗い中で明るくなっている個々の画素において、その明るさに差がある状態となり得る。次の段階において、この画像の全画素の明るさを点検し、特定の明るさのレベルに近いものを比較により検出する。もし、該当すれば、その(X、Y)地点で、そのエネルギーのX線は1カウントありと判定し、そのエネルギーの画像カウンタに追加する。これは、分光分析等の機器分析の信号処理でよく見られる波高分析と類似する。画像の明るさが、強度(光の明るさ)ではなく、エネルギー(光の波長)に対応する。1画素に1光子しか入らないと想定されるくらい暗い画像では、X線のエネルギーに対応して作られる電荷量の差異を見ることができ得る。(例えば、6400eVのX線と3200eVのX線では、同じ1光子が入っても、電荷発生量は2倍違う。従って、画像で見れば、その明るさが2倍違うことになる。)測定結果という画像を得るためには、上述するような工程において、短い時間で、休みなく読み続け、読みながら傍らで信号処理を行い、その結果を、並列している多数の画像カウンタに積算する。これにより、多数のエネルギー別の画像が作られる。
 ここで、カメラのような検出器において、一般的に電気ノイズのようなものと信号を区別するための閾値が設定され、使われている。多くの場合において、このような閾値は、lower level discriminator(LLD)と呼ばれる。本発明の実施例においても、このような閾値は用いられている。非常に短い時間で撮像した一見、全部がノイズ画像に見えるような画像において、信号を得ているピクセルが孤立しているという判定は、すべてのピクセルについて、その電荷量がLLDより上であるかどうかを検査することにより、可能となる。本発明の実施例において、そのような一見ノイズしか写っていないように見える画像(「single photon counting mode(一光子計数法)」と呼ぶ。)だけを使って情報を取得している。ここで、電荷量は、そのピクセルの明るさに相当するものである。このとき用いられるLLDは、ノイズより明らかに明るいといった基準のものであってよい。厳密に値を選択しなければならない訳ではないので、測定条件や得たい画像の種類等を考慮して、各種の測定ごとに設定することができる。
 本発明の実施例において、新規導入した閾値は、LLDよりもノイズレベルにずっと近いと考えられるものである。本発明者らが知る限りにおいて、このような閾値を使った信号処理は知られていない。LLD以上のものであれば、ノイズと区別された明瞭な信号と判定されているので、そのような信号だけを見れば十分であると考えられるからである。しかしながら、ここで、導入してきた第2の閾値と呼ぶものは、一旦は、孤立しているとみなした対象となるピクセルについて、その周囲のピクセルを念入りに検査するために用いられてよい。そして、一旦は、ゼロと見なしてきた周囲のピクセルについて、再度検査することにより、対象となるピクセルからわずかでも隣にしみだしていたら、一旦は孤立しているとみなしたピクセルであっても、孤立しているとするピクセルから除外するという考えである。これは、孤立したピクセルの明るさをもって総電荷量とみなし、それからX線のエネルギーを求めるにあたって、わずかでもしみだしを見逃してしまえば、X線のエネルギーを低めに見積もる結果を招くおそれがあるからである。このような基準を採用したため、孤立した1個のピクセルに電荷があると認められるケース(電荷分割の頻度がそれほど多くはないケース、比較的ピクセルサイズが大きい場合等)では、ほぼ完全なX線スペクトルが得られるという効果が認めらえる。
 本発明のX線撮像装置によれば、可視光波長用に設計されたCCD撮像素子またはCMOS撮像素子、若しくは半導体ピクセル放射線検出器で起こっているすべてのX線フォトン事象を処理することで、X線エネルギーの情報を個別に分離でき、緻密な撮像画素情報が得られる。一般に、デジタルカメラの信号処理は、(1)X線等を含む光の信号を検知して電気信号(電荷)に変えて貯蔵し(センサー)、(2)そのセンサーの画素毎に蓄えられた電荷量を読み出し(読み出し回路、A/D変換器)、(3)その読んだものを画像として表示する、若しくは、電子データとして保存するといった工程が含まれる。これらの各工程は、更に細部の工程に分岐しすることがある。また、一定の特色を持つ細部の工程もあり得る。上述するセンサーとしては、CCDやCMOS等が例示されるが、これらのセンサーの素子構造は物理的に大きく異なっているが、光の画像を得るという最終機能は類似する。電荷量を読み出す方式等については、大きな違いがあるものの、得られた画像の処理については差異がないともいえる。電気信号への変換や電荷量の読み出しの工程は、一般には、ハードウエアイメージで語られることが多く、読み込んだ画像等のデータの表示等においては、ソフトウエア(コンピュータ上で動作するプログラムによる処理)イメージで語られることが多い。しかしながら、本発明の実施例では、上述するように、ソフトウエアだけでなく、個々の特色ある信号処理(画像処理)の種々の機能を備える電子回路というように、ハードウエアイメージでも説明している。このため、行っている処理の機能を明確に説明することが容易であるばかりでなく、ハードウエアを使用する方が類似した機能をより高速に行うことが可能と考えられている。例えば、画像処理専用プロセッサを搭載している場合、コンピュータの画像処理は、高速にできるとされている。このように、機器構成として、コンピュータを含む場合も、含まない場合もあり、どちらの場合にも、個々の機能をハードウエアイメージで保有し、かつ、ソフトウエアでも代替機能を簡易的に実現できる仕組みを持つという特徴がある。そして、これらは、ハードウエアにおいても、また、ソフトウエアにおいても実行が可能であり、広く産業の発展に資する。
1 X線管
2 結晶モノクロメータ
3 信号制御メカニカルシャッター
4、204 試料ホルダー
5、205 試料
6 結像光学系(ピンホール板)
7 X線カメラ(CCDカメラ、CMOSカメラ)
8、208 撮像センサ(二次元半導体センサ)
82 単位受光素子
84 周辺受光素子
9 冷却回路
10、20、30 X線画像生成回路
101a、101b X線ビーム
102 有効事象判断回路
103 孤立画素抽出回路
104 ノイズ除去回路
105 エネルギー準位判断回路
106 画素別カウンタ回路
107 電荷分割検出回路
108 電荷分割処理回路
109 X線カウンタ画像生成回路
11 X線画像表示装置
112、115 第1のしきい値
120 第2のしきい値

 

Claims (16)

  1.  分析対象物からのX線光子を含むエネルギー線を受光するX線撮像素子であって、当該X線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該X線光子の受光による電荷共有を生ずる二次元的領域を有する前記X線撮像素子と、
     このX線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該X線撮像素子内での二次元配置に応じて撮像された分析対象物のX線画像として生成するX線画像生成回路と、を備え、
     前記X線画像生成回路は、
      前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第1のしきい値を比較する有効事象判断回路と、
      前記有効事象判断回路で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号およびノイズと認識する基準となる第2のしきい値を比較するノイズ除去回路と、
      前記ノイズ除去回路で除去されないで有効な光子事象と認識された単位受光素子について、前記X線画像生成回路の画素別カウンタ値をプラス1する計数を行う画素別カウンタ回路と、
     この画素別カウンタ回路の計数値に基づいて、前記X線撮像素子によるX線画像を生成するX線カウンタ画像生成回路と、
     を備えることを特徴とするX線撮像装置。
  2.  前記X線画像生成回路は、さらに、前記有効事象判断回路で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、全ての周辺受光素子の受光信号が第1のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出回路を有し、
     前記ノイズ除去回路は、前記有効事象判断回路で有効な光子事象と認識と認識された単位受光素子であって、前記孤立画素抽出回路での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第2のしきい値とを比較することを特徴とする請求項1に記載のX線撮像装置。
  3.  分析対象物からのX線光子を含むエネルギー線を受光するX線撮像素子であって、当該X線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該X線光子の受光による電荷共有を生ずる二次元的領域を有する前記X線撮像素子と、
     このX線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該X線撮像素子内での二次元配置に応じて撮像された分析対象物のX線画像として生成するX線画像生成回路と、を備え、
     前記X線画像生成回路は、
      前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第1のしきい値を比較する有効事象判断回路と、
      前記有効事象判断回路で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出回路と、
      前記電荷分割検出回路で検出された電荷分割の生じた単位受光素子および周辺受光素子について、前記第1のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理回路と、
     前記電荷分割処理回路で演算された前記中央値、積分値およびその位置情報を用いて、前記第1のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス1する計数を行う画素別カウンタ回路と、
     この画素別カウンタ回路の計数値に基づいて、前記X線撮像素子によるX線画像を生成するX線カウンタ画像生成回路と、
     を備えることを特徴とするX線撮像装置。
  4.  前記X線撮像素子は、CCD撮像素子、CMOS撮像素子または半導体ピクセル放射線検出器であることを特徴とする請求項1乃至3の何れか1項に記載のX線撮像装置。
  5.  前記X線撮像素子の単位受光素子は、一画素として1乃至10μmの代表長さを有する二次元的領域を有することを特徴とする請求項4に記載のX線撮像装置。
  6.  前記第1のしきい値は、前記X線撮像素子の各単位受光素子で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値であることを特徴とする請求項1乃至5の何れか1項に記載のX線撮像装置。
  7.  前記第2のしきい値は、前記X線撮像素子の各単位受光素子で生ずるバックグランドノイズよりも大きな値であって、前記第1のしきい値よりも小さな値であることを特徴とする請求項1又は2に記載のX線撮像装置。
  8.  前記周辺受光素子は、光子事象のあった前記単位受光素子を中心とする3x3、5x5又は7x7の何れかあることを特徴とする請求項1乃至7の何れか1項に記載のX線撮像装置。
  9.  さらに、前記X線画像生成回路の各単位受光素子毎に読み込まれた受光信号について、当該受光信号のエネルギー準位を求めるエネルギー準位判断回路を有することを特徴とする請求項1乃至8の何れか1項に記載のX線撮像装置。
  10.  さらに、前記X線画像生成回路の画素別カウンタ回路は、前記エネルギー準位判断回路で分光されたエネルギー準位に応じて区分されていることを特徴とする請求項9に記載のX線撮像装置。
  11.  前記X線画像生成回路の画素別カウンタ回路は、光子事象のあった前記単位受光素子の受光信号のエネルギー準位について、前記単位受光素子の周辺受光素子の受光信号のエネルギー準位も加算することを特徴とする請求項10に記載のX線撮像装置。
  12.  さらに、前記分析対象物にX線を含むエネルギー線を照射するX線照射回路を有することを特徴とする請求項1乃至11の何れか1項に記載のX線撮像装置。
  13.  さらに、前記X線画像生成回路の画素別カウンタ値をリセットする画素別カウンタリセット回路を有することを特徴とする請求項1乃至12の何れか1項に記載のX線撮像装置。
  14.  分析対象物からのX線光子を含むエネルギー線を受光するX線撮像素子の画像処理方法であって、前記X線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該X線光子の受光による電荷共有を生ずる二次元的領域を有するものであって、
     前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第1のしきい値とを比較する有効事象判断工程と、
     前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第2のしきい値とを比較するノイズ除去工程と、
     前記ノイズ除去工程で除去されないで有効な光子事象と認識された単位受光素子について、当該単位受光素子に対応する画素別カウンタのカウント値をプラス1する計数を行う画素別カウンタ工程と、
     この画素別カウンタ工程での計数値に基づいて、前記X線撮像素子によるX線画像を生成することを特徴とするX線撮像素子の画像処理方法。
  15.  請求項14に記載のX線撮像素子の画像処理方法において、さらに、
      前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、全ての周辺受光素子の受光信号が第1のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出工程を有し、
      前記ノイズ除去工程は、前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子であって、前記孤立画素抽出工程での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第2のしきい値とを比較することを特徴とするX線撮像素子の画像処理方法。
  16.  分析対象物からのX線光子を含むエネルギー線を受光するX線撮像素子の画像処理方法であって、前記X線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該X線光子の受光による電荷共有を生ずる二次元的領域を有するものであって、
     前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第1のしきい値とを比較する有効事象判断工程と、
     前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第1のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出工程と、
     前記電荷分割検出工程で検出された電荷分割の生じた単位受光素子と周辺受光素子について、前記第1のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理工程と、
     前記電荷分割処理工程で演算された前記中央値、積分値およびその位置情報を用いて、前記第1のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス1する計数を行う画素別カウンタ工程と、
     この画素別カウンタ工程での計数値に基づいて、前記X線撮像素子によるX線画像を生成することを特徴とするX線撮像素子の画像処理方法。

     
PCT/JP2018/008420 2017-09-28 2018-03-05 X線撮像装置およびx線撮像素子の画像処理方法 WO2019064632A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019544217A JPWO2019064632A1 (ja) 2017-09-28 2018-03-05 X線撮像装置およびx線撮像素子の画像処理方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017187361 2017-09-28
JP2017-187361 2017-09-28

Publications (1)

Publication Number Publication Date
WO2019064632A1 true WO2019064632A1 (ja) 2019-04-04

Family

ID=65901350

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/008420 WO2019064632A1 (ja) 2017-09-28 2018-03-05 X線撮像装置およびx線撮像素子の画像処理方法

Country Status (2)

Country Link
JP (1) JPWO2019064632A1 (ja)
WO (1) WO2019064632A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110596744A (zh) * 2019-09-12 2019-12-20 电子科技大学 一种椭圆弯晶谱仪数据采集与处理系统
CN115684223A (zh) * 2022-11-07 2023-02-03 南京航空航天大学 基于x射线康普顿散射及瑞利散射的成像装置及成像方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004325183A (ja) * 2003-04-23 2004-11-18 M & C:Kk 放射線検出方法、放射線検出器、及び、この検出器を搭載した放射線撮像システム
JP2017092820A (ja) * 2015-11-13 2017-05-25 キヤノン株式会社 放射線撮像装置およびフォトンカウンティングの方法
JP2017521682A (ja) * 2014-06-20 2017-08-03 ブルカー・エイエックスエス・インコーポレイテッドBruker AXS, Inc. 混合型光子計数/アナログ出力モードで動作可能なx線検出器

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011015235A1 (en) * 2009-08-04 2011-02-10 Organisation Europeenne Pour La Recherche Nucleaire Pixilated radiation sensing device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004325183A (ja) * 2003-04-23 2004-11-18 M & C:Kk 放射線検出方法、放射線検出器、及び、この検出器を搭載した放射線撮像システム
JP2017521682A (ja) * 2014-06-20 2017-08-03 ブルカー・エイエックスエス・インコーポレイテッドBruker AXS, Inc. 混合型光子計数/アナログ出力モードで動作可能なx線検出器
JP2017092820A (ja) * 2015-11-13 2017-05-25 キヤノン株式会社 放射線撮像装置およびフォトンカウンティングの方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110596744A (zh) * 2019-09-12 2019-12-20 电子科技大学 一种椭圆弯晶谱仪数据采集与处理系统
CN115684223A (zh) * 2022-11-07 2023-02-03 南京航空航天大学 基于x射线康普顿散射及瑞利散射的成像装置及成像方法
CN115684223B (zh) * 2022-11-07 2023-12-15 南京航空航天大学 基于x射线康普顿散射及瑞利散射的成像装置及成像方法

Also Published As

Publication number Publication date
JPWO2019064632A1 (ja) 2020-07-16

Similar Documents

Publication Publication Date Title
JP4960255B2 (ja) ガンマ線撮像装置
US8363781B2 (en) Nondestructive identification method and nondestructive identification device
EP3026429B1 (en) X-ray fluorescence analyzer and x-ray fluorescence analyzing method
Jacques et al. A laboratory system for element specific hyperspectral X-ray imaging
JP6769402B2 (ja) 電子線マイクロアナライザー及びデータ処理プログラム
Losko et al. New perspectives for neutron imaging through advanced event-mode data acquisition
US10209375B2 (en) X-ray data processing apparatus, X-ray data processing method, and X-ray data processing program
Zhao et al. CCD camera as feasible large-area-size x-ray detector for x-ray fluorescence spectroscopy and imaging
JP4853964B2 (ja) X線検査装置用のx線センサ
Zhao et al. Seeing elements by visible-light digital camera
WO2019064632A1 (ja) X線撮像装置およびx線撮像素子の画像処理方法
Romano et al. A new X-ray pinhole camera for energy dispersive X-ray fluorescence imaging with high-energy and high-spatial resolution
Naselli et al. Quantitative analysis of an ECR Ar plasma structure by X-ray spectroscopy at high spatial resolution
JP7150638B2 (ja) 半導体欠陥検査装置、及び、半導体欠陥検査方法
CN108140650A (zh) 具有高空间解析度的x射线检测器
Ponchut Characterization of X-ray area detectors for synchrotron beamlines
Vulovic et al. A toolkit for the characterization of CCD cameras for transmission electron microscopy
Goldstein et al. Energy dispersive X-ray spectrometry: physical principles and user-selected parameters
CN106290433A (zh) X射线数据处理装置及其方法以及程序
Cao et al. The measurement of the presampled MTF of a high spatial resolution neutron imaging system
JP2012058146A (ja) X線検出システム
KR20190071111A (ko) 엑스선 검사 장비 및 이를 이용하는 반도체 장치 제조 방법
Mor et al. Development of a CCD based thermal neutron imaging detector for the Israeli Research Reactor IRR-1 at Soreq NRC
JP2008089325A (ja) 電子プローブマイクロアナライザによる地質年代測定法
Uher et al. X-ray beam hardening based material recognition in micro-imaging

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: 18861708

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019544217

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18861708

Country of ref document: EP

Kind code of ref document: A1