Nothing Special   »   [go: up one dir, main page]

WO2020129213A1 - Ｘ線撮影装置 - Google Patents

Ｘ線撮影装置 Download PDF

Info

Publication number
WO2020129213A1
WO2020129213A1 PCT/JP2018/046996 JP2018046996W WO2020129213A1 WO 2020129213 A1 WO2020129213 A1 WO 2020129213A1 JP 2018046996 W JP2018046996 W JP 2018046996W WO 2020129213 A1 WO2020129213 A1 WO 2020129213A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
pixel
ray
resolution
ray fluoroscopic
Prior art date
Application number
PCT/JP2018/046996
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 JP2020560725A priority Critical patent/JP7205552B2/ja
Priority to US17/414,069 priority patent/US11821851B2/en
Priority to CN201880099878.1A priority patent/CN113226180B/zh
Priority to PCT/JP2018/046996 priority patent/WO2020129213A1/ja
Priority to TW108124892A priority patent/TWI708937B/zh
Publication of WO2020129213A1 publication Critical patent/WO2020129213A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating 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/02—Investigating 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 transmitting the radiation through the material
    • G01N23/04—Investigating 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 transmitting the radiation through the material and forming images of the material
    • G01N23/046—Investigating 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 transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
    • G06T3/4069—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution by subpixel displacements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/401—Imaging image processing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/42—Imaging image digitised, -enhanced in an image processor
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/427—Imaging stepped imaging (selected area of sample is changed)
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/66—Specific applications or type of materials multiple steps inspection, e.g. coarse/fine

Definitions

  • the present invention relates to an X-ray imaging apparatus, and in particular, generates a super-resolution image having a higher resolution than the first image and the second image generated based on the X-ray detected at the first position and the second position.
  • the present invention relates to an X-ray imaging apparatus.
  • an X-ray generator that emits an X-ray beam, a radiation detector that has a plane that is arranged perpendicular to the emission direction of the X-ray beam, and a radiation detector that holds the radiation detector
  • An X-ray imaging apparatus including a translation stage that translates a detector in the pixel direction is disclosed.
  • the radiation detector is configured to move to different detection positions by being translated by the translation stage in a distance smaller than the pixel size. Then, it is configured to generate an image (super-resolution image) having a finer resolution than the acquired radiograph based on the plurality of radiographs (images) acquired at different detection positions.
  • a super-resolution image is generated based on a plurality of acquired images like the conventional X-ray imaging apparatus described in Patent Document 1 described above.
  • a sequential calculation is performed using a large number of parameters.
  • the calculation time tends to be relatively long.
  • the quality of the generated super-resolution image is easily influenced by the parameters, and a certain super-resolution effect (the effect of increasing the resolution without degrading the image quality) cannot be obtained. There is.
  • the present invention has been made to solve the above problems, and an object of the present invention is to suppress the increase in the calculation time when generating a super-resolution image and to maintain a constant super-resolution.
  • An object of the present invention is to provide an X-ray imaging apparatus capable of generating a super-resolution image while ensuring a resolution effect.
  • an X-ray imaging apparatus provides an X-ray source and an X-ray emitted from the X-ray source from a first position and from one pixel from the first position. Is detected at the second position, which is translated in the first direction by a small amount of movement, and a first image is generated based on the X-ray detected at the first position, and is detected at the second position. And an image processing unit that generates a second image based on the X-rays, and a moving mechanism that moves the detection position of the X-rays by the detector between the first position and the second position.
  • the first image in the other of the first image and the second image is divided into the first direction based on the pixel values of two pixels that overlap with the first pixel. It is configured to generate a super-resolution image having a higher resolution in the first direction than the first image and the second image.
  • the image processing unit determines the pixel value of the first pixel of one of the first image and the second image as the first image and the second image. Based on the pixel values of the two pixels that overlap the first pixel in the other of the first image and the second image when the images are displayed in a shifted manner in the first direction by the amount corresponding to the movement amount, It is configured to generate a super-resolution image having a higher resolution in the first direction than the first image and the second image, based on the divided image that has been subjected to the first-direction division processing for dividing in one direction.
  • the first pixel in the other of the first image and the second image is overlaid.
  • the two pixels to be wrapped are pixels that are partially displaced to one side and the other side of the first pixel in the first direction, respectively. That is, two pixels that overlap the first pixel in the other of the first image and the second image are a guideline for dividing the pixel value of the first pixel into one side and the other side in the first direction. Can be thought of as the value of.
  • the pixel value of the first pixel is divided based on the pixel values of the two pixels, and the distribution of the actual pixel value is taken into consideration in consideration of the distribution of the actual pixel value. It can be divided in one direction. That is, since the calculation is performed only by dividing the pixel value based on the pixel values of the two pixels, it is not necessary to use the sequential calculation and the large number of parameters that are performed in the conventional super-resolution processing.
  • the pixels are divided in consideration of the actual distribution of pixel values, a super-resolution image having a higher resolution in at least the first direction than that of the first image and the second image and a constant image quality is generated. be able to. As a result, it is possible to prevent the calculation time for generating the super-resolution image from increasing and to maintain a certain super-resolution effect (the effect of increasing the resolution without degrading the image quality) while maintaining the super-resolution image. Can be generated.
  • the image processing unit sets the pixel value of the first pixel of one of the first image and the second image to the movement amount of the first image and the second image. Based on the pixel value of two pixels that overlap with the first pixel in the other of the first image and the second image and the overlapping area ratio when they are overlapped and displayed in the first direction by a corresponding amount.
  • the super-resolution image is generated based on the divided image subjected to the division process.
  • the first pixel is added to the pixel values of the two pixels that overlap the first pixel, and the overlapping area ratio of the two pixels is taken into consideration. It is possible to divide in the first direction by more accurately considering the distribution of values. As a result, compared with the case where the overlapping area ratio of two pixels is not considered, a super-resolution image having higher resolution and better image quality in at least the first direction than the first image and the second image is generated. be able to.
  • the image processing unit sets the pixel value of the first pixel in the first image to a first value corresponding to the movement amount of the first image and the second image.
  • the first divided image obtained by performing the division processing based on the pixel values of the two pixels that overlap the first pixel in the second image when the images are displayed in a superimposed manner by shifting in the direction, and the second divided image in the second image.
  • the super-resolution image is generated based on the second divided image obtained by performing the division processing based on the pixel value of one pixel.
  • both the first divided image and the second divided image are considered in comparison with the case where the super-resolution image is generated by considering only one of the first divided image and the second divided image.
  • the resolution in at least the first direction is higher than that of the first image and the second image as compared with the case of generating the super-resolution image in consideration of only one of the first divided image and the second divided image. It is possible to generate a super-resolution image with higher image quality.
  • the image processing unit is configured to generate a super-resolution image by generating an average image obtained by averaging the first divided image and the second divided image in pixels corresponding to each other.
  • the image processing unit is configured to generate a super-resolution image by generating an average image obtained by averaging the first divided image and the second divided image in pixels corresponding to each other.
  • the detector is, in addition to the first position and the second position, moved from the first position and the second position by a movement amount smaller than one pixel in the first direction.
  • the image processing unit is configured to detect at the first position, the second position, the third position, and the fourth position.
  • the first image, the second image, the third image, and the fourth image are generated based on the detected X-rays, respectively, and the division processing in the first direction is performed based on the first image and the second image.
  • a post-process for displaying a super-resolution image on the display device such as an image interpolation process for adjusting the pixel size in the direction and the second direction, can be simplified.
  • the image processing unit has the first image and the second image corresponding to the movement amount at the other end of the first image or the second image.
  • the pixel of one pixel that overlaps Based on the value, it is configured to virtually generate another overlapping pixel. According to this structure, by virtually generating another overlapping pixel, even if there is only one overlapping pixel with respect to the first pixel located at the end, the first pixel is The calculation for dividing in the first direction can be reliably performed.
  • the detector preferably translates the X-ray emitted from the X-ray source in the first direction in parallel with the first position by a half length of one pixel from the first position. And the second position.
  • the overlapping area ratios of the two pixels become equal (become 1:1), so that the pixel value of the first pixel is the pixel value of the two pixels that overlap the first pixel. It is possible to divide in the first direction by the ratio of the values themselves.
  • the calculation in the division processing can be simplified as compared with the case where the overlapping area ratio of two pixels is divided into a ratio other than two. Therefore, calculation when generating a super-resolution image It is possible to further prevent the time from increasing.
  • the detector detects the X-rays emitted from the X-ray source from a plurality of directions while rotating around the third direction as a rotation axis. It is configured to perform tomography by alternately repeating parallel movements in three directions. With such a configuration, in the X-ray imaging apparatus that performs the above-described (so-called non-helical scan type) tomography, it is possible to suppress an increase in calculation time when generating a super-resolution image and to maintain a constant value. A super-resolution image can be generated while ensuring the super-resolution effect.
  • FIG. 1 is a diagram showing the overall configuration of an X-ray imaging apparatus according to an embodiment of the present invention, viewed from the Y direction. It is the figure seen from the X direction showing the whole X-ray imaging system composition by one embodiment of the present invention.
  • FIG. 6 is a diagram for explaining an X-ray detection position in the X-ray imaging apparatus according to the embodiment of the present invention.
  • FIG. 3 is a diagram for explaining an X-ray fluoroscopic image captured by the X-ray imaging apparatus according to the embodiment of the present invention.
  • FIG. 7 is a diagram for explaining super-resolution processing in the X-ray imaging apparatus according to the embodiment of the present invention.
  • FIG. 1 is a diagram showing the overall configuration of an X-ray imaging apparatus according to an embodiment of the present invention, viewed from the Y direction. It is the figure seen from the X direction showing the whole X-ray imaging system composition by one embodiment of the present invention.
  • FIG. 6 is a diagram for explaining an X-ray detection position
  • FIG. 6 is a diagram showing a part of a first X-ray fluoroscopic image and a part of a second X-ray fluoroscopic image on which super-resolution processing in the X direction is performed. It is a figure for demonstrating the 1st division image and the 2nd division image which were generated based on the 1st X-ray fluoroscopic image and the 2nd X-ray fluoroscopic image. It is a figure for demonstrating the average image which averaged the 1st division image and the 2nd division image. It is a figure for demonstrating the super-resolution process in the X direction with respect to a 3rd X-ray fluoroscopic image and a 4th X-ray fluoroscopic image.
  • the X-ray imaging apparatus 100 includes an X-ray source 1, a detector 2, a processing unit 3, a rotary stage 4, and a detector moving mechanism 5.
  • the detector moving mechanism 5 is an example of the "moving mechanism" in the claims.
  • the X-ray source 1, the rotary stage 4, and the detector 2 are arranged side by side in this order in the direction (Z direction) connecting the X-ray source 1 and the detector 2. ..
  • the direction from the X-ray source 1 to the detector 2 is the Z2 direction, and the opposite direction is the Z1 direction.
  • the direction in which the rotary stage 4 translates is referred to as the Y direction.
  • a direction orthogonal to the Z direction and the Y direction is defined as the X direction.
  • the X direction is an example of the “first direction” and the “row direction” in the claims.
  • the Y direction is an example of the “second direction” and the “column direction” in the claims.
  • the X-ray source 1 is an X-ray generator capable of generating X-rays by applying a high voltage.
  • the X-ray source 1 is configured to emit the generated X-rays in the Z2 direction.
  • the detector 2 detects X-rays emitted from the X-ray source 1 and converts the detected X-rays into electric signals.
  • the detector 2 has a detection surface 2b including a plurality of detection elements 2a arranged side by side in a matrix in a predetermined cycle L1 in the X and Y directions.
  • the detector 2 is arranged so as to correspond to each of the plurality of detection elements 2a, and includes a plurality of conversion elements (not shown) that convert the detected X-rays into electric signals.
  • the detector 2 is, for example, an FPD (Flat Panel Detector).
  • the detection signal (electrical signal) converted by the detector 2 is sent to an image processing unit 3b (described later) included in the processing unit 3.
  • the processing unit 3 includes a control unit 3a and an image processing unit 3b.
  • the control unit 3a is configured to control the operations of the rotary stage 4 and the detector moving mechanism 5.
  • the control unit 3a includes, for example, a CPU (Central Processing Unit), a ROM (Read Only Memory), a RAM (Random Access Memory), and the like.
  • the image processing unit 3b is configured to generate an X-ray fluoroscopic image Ia (see FIG. 4) based on the detection signal sent from the detector 2. Further, the image processing unit 3b is configured to generate a reconstructed image 70x (see FIG. 8) described later based on the plurality of X-ray fluoroscopic images Ia (see FIG. 4).
  • the image processing unit 3b includes, for example, a processor such as a GPU (Graphics Processing Unit) or an FPGA (Field-Programmable Gate Array) configured for image processing.
  • the data of each pixel Ea(E) in the X-ray fluoroscopic image Ia (see FIG. 4) corresponds to the X-ray dose detected by each detection element 2a of the detector 2.
  • the rotary stage 4 has a mounting surface 4a (see FIG. 2) on which the subject S is mounted.
  • the rotary stage 4 is configured to be capable of rotating 360 degrees around the rotation axis 90 along the Y direction in the XZ plane. Further, as shown in FIG. 2, the rotary stage 4 is configured to be capable of translational movement in the Y direction.
  • the detector 2 detects the X-rays emitted from the X-ray source 1 from a plurality of directions (in the XZ plane) while rotating around the Y direction as a rotation axis, and The tomography can be performed by alternately repeating the parallel movement in the Y direction. That is, the X-ray imaging apparatus 100 is an imaging apparatus capable of performing so-called non-helical scan tomography.
  • the detector moving mechanism 5 is configured to move the detector 2 in the X direction and the Y direction by a moving amount smaller than that of the detecting element 2a.
  • the detector moving mechanism 5 moves the detector 2 among the first position P1, the second position P2, the third position P3, and the fourth position P4.
  • the second position P2 is a position translated in the X direction by a distance L2 from the first position P1.
  • the third position P3 is a position translated in the Y direction by a distance L3 from the first position P1.
  • the fourth position is a position translated from the first position P1 by a distance L3 to the same side as the third position in the Y direction.
  • the distance L2 and the distance L3 are smaller than the size of one detection element 2a (the period in which the detection elements 2a are arranged) L1 in the X direction and the Y direction, respectively.
  • the distance L2 and the distance L3 are half the length (length) of the size L1 of one detection element 2a in the X and Y directions, respectively. That is, in the present embodiment, the detector 2 applies the X-ray emitted from the X-ray source 1 to the first position P1 and one pixel Ea from the first position P1 (the length L4 (see FIG. 5)). Of the second position P2 that is translated in the X direction by a length L5 (see FIG.
  • the image processing unit 3b based on the X-rays detected at the first position P1, the second position P2, the third position P3, and the fourth position P4, respectively. It is possible to generate the first X-ray fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40 as the fluoroscopic image Ia.
  • the first X-ray fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40 are the “first image”, the “second image”, and the “second image” in the claims, respectively. It is an example of a "3rd image” and a "4th image.”
  • the first X-ray fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40 are each a plurality of X-ray fluoroscopic images Ia acquired by tomography using the rotary stage 4. including. That is, the first X-ray fluoroscopic image 10 includes first X-ray fluoroscopic images 11, 12, 13, 14,... Which are different from each other in at least one of the imaging position in the XZ plane and the imaging position in the Y direction by the rotary stage 4.
  • the second X-ray fluoroscopic image 20 includes second X-ray fluoroscopic images 21, 22, 23, 24,...
  • the third X-ray fluoroscopic image 30 includes third X-ray fluoroscopic images 31, 32, 33, 34,... Which differ from each other in at least one of the imaging direction in the XZ plane and the imaging position in the Y direction due to the operation of the rotary stage 4.
  • the fourth X-ray fluoroscopic image 40 includes fourth X-ray fluoroscopic images 41, 42, 43, 44,... Which differ from each other in at least one of the imaging position in the XZ plane and the imaging position in the Y direction due to the operation of the rotary stage 4.
  • the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 are data in which the imaging positions differ in the X direction by a size L5 that is half the size L4 of one pixel Ea.
  • the third X-ray fluoroscopic image 30 and the fourth X-ray fluoroscopic image 40 are data in which imaging positions differ in the X direction by a size L5 that is half the size L4 of one pixel Ea.
  • the first X-ray fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40 have a size L5 that is half the size L4 of one pixel Ea in the Y direction.
  • FIG. 5 This is the data in which the shooting position differs by the minute.
  • the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 are superimposed and displayed while being shifted in the X direction by an amount corresponding to the movement amount (of the X-ray detection position), and the third X-ray image.
  • the perspective image 30 and the fourth X-ray fluoroscopic image 40 are shown in a state of being overlapped with each other by being shifted in the X direction by an amount corresponding to the movement amount (of the X-ray detection position).
  • the image processing unit 3b performs the first X-ray fluoroscopic image based on the first X-ray fluoroscopic image 10 (X-ray fluoroscopic image Ia) and the second X-ray fluoroscopic image 20 (X-ray fluoroscopic image Ia). It is configured to generate a super-resolution image Ib having a higher resolution in the X direction than the image 10 and the second X-ray fluoroscopic image 20.
  • the process of generating an image having a higher resolution than the original image based on a plurality of images (original image) may be referred to as “super-resolution process”.
  • the image processing unit 3b performs a division image in the X direction in which the pixel value of the pixel E1 (see FIG. 6A) in the first X-ray fluoroscopic image 10 is divided in the X direction.
  • 51 see FIG. 6B
  • a divided image 52 obtained by performing a division process in the X direction for dividing the pixel value of the pixel E2 (see FIG. 6A) in the second X-ray fluoroscopic image 20 into the X direction.
  • the super-resolution image Ib (see FIG. 5) having a higher resolution in the X direction than the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 is generated.
  • the image processing unit 3b generates an average image 60 (see FIG. 6C) obtained by averaging the divided image 51 (see FIG. 6B) and the divided image 52 (see FIG. 6B) in the pixels E corresponding to each other. It is configured to generate the super-resolution image Ib (see FIG. 5).
  • the pixel E1 and the pixel E2 are examples of the “first pixel” and the “second pixel” in the claims, respectively.
  • the divided image 51 and the divided image 52 are examples of the “first divided image” and the “second divided image” in the claims, respectively.
  • the pixel value of the pixel E1 (see FIG. 6A) in the first X-ray fluoroscopic image 10 is changed to the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (as shown in FIG. 5).
  • the pixel E1 (see FIG. 6A) in the second X-ray fluoroscopic image 20 overlaps the pixel E1 (see FIG. 6A) in the second X-ray fluoroscopic image 20.
  • the divided image 52 overlaps the pixel value of the pixel E2 (see FIG. 6A) in the second X-ray fluoroscopic image 20 with the pixel E2 (see FIG. 6A) in the first X-ray fluoroscopic image 2.
  • FIG. 6B is a divided image 50 (see FIG. 6B) that has been divided in the X direction based on the pixel value of one pixel E1 and the overlapping area ratio. Note that, in FIG. 6, for convenience of illustration, only one row in the X-ray fluoroscopic image Ia is extracted and described. 6A and 6B, similarly to FIG.
  • the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 are shifted in the X direction by an amount corresponding to the movement amount (of the X-ray detection position).
  • the display is shifted in the vertical direction of the paper for convenience of illustration.
  • the image processing unit 3b divides the pixel E1b in the first X-ray fluoroscopic image 10 into a pixel E51c and a pixel E51d that are pixels of the divided image 51.
  • the image processing unit 3b sets the pixel values of the pixel E51c and the pixel E51d at the ratio of the pixel value of the pixel E2a and the pixel value of the pixel E2b that overlap the pixel E1b in the second X-ray fluoroscopic image 20, respectively.
  • the pixel value of E1b is a divided value.
  • the image processing unit 3b performs the same process on the other pixel E1 in the first X-ray fluoroscopic image 10. That is, the image processing unit 3b divides the pixel E1c in the first X-ray fluoroscopic image 10 into the pixel E51e and the pixel E51f which become the pixel E of the divided image 51. The image processing unit 3b determines the pixel values of the pixel E51e and the pixel E51f by the ratio of the pixel value of the pixel E2b and the pixel value of the pixel E2c that overlap the pixel E1c in the second X-ray fluoroscopic image 20, respectively. The pixel value of E1c is a divided value.
  • the image processing unit 3b divides the pixel E1d in the first X-ray fluoroscopic image 10 into a pixel E51g and a pixel E51h which are pixels of the divided image 51.
  • the image processing unit 3b sets the pixel values of the pixel E51g and the pixel E51h at the ratio of the pixel value of the pixel E2c and the pixel value of the pixel E2d that overlap the pixel E1d in the second X-ray fluoroscopic image 20, respectively.
  • the pixel value of E1d is a divided value.
  • the image processing unit 3b corresponds to the movement amount (of the X-ray detection position) between the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (as shown in FIG. 5).
  • the image processing unit 3b corresponds to the movement amount (of the X-ray detection position) between the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (as shown in FIG. 5).
  • the image processing unit 3b divides the pixel E1a in the first X-ray fluoroscopic image 10 into a pixel E51a and a pixel E51b that become the pixel E of the divided image 51.
  • the image processing unit 3b generates a virtual pixel E2e adjacent to the pixel E2a (at the end of the second X-ray fluoroscopic image 20) that overlaps the pixel E1a in the second X-ray fluoroscopic image 20.
  • the image processing unit 3b then divides the pixel values of the pixel E51a and the pixel E51b by dividing the pixel value of the pixel E1a by the ratio between the pixel value of the pixel E2a and the pixel value of the (virtual) pixel E2d. And In the X-ray imaging apparatus 100, the image processing unit 3b is configured to generate the pixel value of the virtual pixel E so as to have the same pixel value as the pixel value of the adjacent pixel E.
  • the image processing unit 3b is configured to generate the pixel E at the end of the divided image 51 so that the divided image 51 and the divided image 52 have the same number of pixels. Specifically, the image processing unit 3b generates the pixel E51i based on the (virtual) pixel E1e adjacent to the (end) pixel E1d in the first X-ray fluoroscopic image 10. With the above configuration, the divided image 51 is generated by dividing the pixel E1 of the first X-ray fluoroscopic image 10 in the X direction.
  • the image processing unit 3b also generates a divided image 52 obtained by dividing the pixel E2 of the second X-ray fluoroscopic image 20 in the X direction by the same method (as the generation of the divided image 51). That is, the image processing unit 3b divides the pixel E2a in the second X-ray fluoroscopic image 20 into the pixel E52b and the pixel E52c that become the pixel E of the divided image 52. Further, the image processing unit 3b divides the pixel E2b in the second X-ray fluoroscopic image 20 into a pixel E52d and a pixel E52e which are the pixels E of the divided image 52.
  • the image processing unit 3b divides the pixel E2c in the second X-ray fluoroscopic image 20 into a pixel E52f and a pixel E52g that are pixels of the divided image 52.
  • the pixel E2d in the second X-ray fluoroscopic image 20 is divided into a pixel E52h and a pixel E52i that become the pixel E of the divided image 52.
  • the image processing unit 3b is based on the (virtual) pixel E2e adjacent to the (end) pixel E2a in the second X-ray fluoroscopic image 20 so that the divided images 51 and 52 have the same number of pixels. Then, the pixel E52a is generated. With the above configuration, the divided image 52 is generated by dividing the pixel E2 of the second X-ray fluoroscopic image 20 in the X direction.
  • the image processing unit 3b generates an average image 61 (60) obtained by averaging the pixel values of the corresponding pixels E of the divided image 51 and the divided image 52. That is, the image processing unit 3b averages the pixel value (1.5) of the pixel E51a of the divided image 51 and the pixel value (2.5) of the pixel E52a of the divided image 52 to obtain the pixel E61a of the average image 61. Pixel value (2).
  • the image processing unit 3b performs the same process for other pixels.
  • the image processing unit 3b determines the pixel values of the pixels E51b, E51c, E51d, E51e, E51g, E51h, and E51i of the divided image 51 and the pixels E52b, E52c, E52d, E52e, E52g, E52h, and E52i of the divided image 52.
  • the respective pixel values are averaged to obtain the pixel values of the pixels E61b, E61c, E61d, E61e, E61g, E61h, and E61i of the average image 61.
  • the average image 61 (60) has a higher resolution (doubled) in the X direction than the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20.
  • the length L5 of one pixel Eb of the average image 61 (60) is half the length L4 of one pixel Ea of the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20.
  • the average image 61 (60) is a super-resolution image Ib having a higher resolution in the X direction than the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (super-resolution processing has been performed in the X direction). is there.
  • the image processing unit 3b causes the first X-ray fluoroscopic image 10 (X-ray fluoroscopic image Ia), the second X-ray fluoroscopic image 20 (X-ray fluoroscopic image Ia), and the third X-ray fluoroscopic image 30( Based on the X-ray fluoroscopic image Ia) and the fourth X-ray fluoroscopic image 40 (X-ray fluoroscopic image Ia), the first X-ray fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray. It is configured to generate a super-resolution image Ic having higher resolution in the X and Y directions than the perspective image 40.
  • the image processing unit 3b includes the divided image 50 (divided image 51 and divided image 52) that has been divided in the X direction based on the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20, and the third X image.
  • the divided image 50 divided image 53 and divided image 54
  • the same process as the division process in the X direction is performed. It is configured to generate the super-resolution image Ic by performing it in the direction.
  • the divided image 53 and the divided image 54 are examples of the “first divided image” and the “second divided image” in the claims, respectively.
  • the image processing unit 3b is similar to the divided image 51 and the divided image obtained by performing the division processing in the X direction based on the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20. Then, the divided image 53 and the divided image 54 are generated by performing the division processing in the X direction based on the third X-ray transparent image 30 and the fourth X-ray transparent image 40. That is, the image processing unit 3b performs division processing in the X direction that divides the pixel value of the pixel E3 in the third X-ray fluoroscopic image 31 (30) in the X direction to generate the divided image 53, and also the fourth X-ray fluoroscopic image.
  • a division image 54 is generated by performing division processing in the X direction that divides the pixel value of the pixel E4 in 41 (40) in the X direction.
  • the pixel E3 and the pixel E4 are examples of the “first pixel” and the “second pixel” in the claims, respectively.
  • the image processing unit 3b generates an average image 62 (60) obtained by averaging the pixel values of the corresponding pixels E of the divided image 53 and the divided image 54. That is, the image processing unit 3b averages the pixel value of the pixel E53 of the divided image 53 and the pixel value of the pixel E54 of the divided image 54 (corresponding to the pixel E53) to obtain the pixel E62 of the average image 62.
  • the resolution is higher in the X direction than the third X-ray fluoroscopic image 30 and the fourth X-ray fluoroscopic image 40 (doubled). Average image 62 (60) is generated.
  • the average image 62 (60) is a super-resolution image Ib having a higher resolution in the X direction than the third X-ray fluoroscopic image 30 and the fourth X-ray fluoroscopic image 40 (super-resolution processing has been performed in the X direction). is there.
  • the image processing unit 3b divides in the Y direction based on the average image 61 and the average image 62 (which are the super-resolution images Ib subjected to the super-resolution processing in the X direction). Processing is performed to generate the divided image 55 and the divided image 56. That is, the image processing unit 3b performs the Y-direction division process of dividing the pixel value of the pixel E61 in the average image 61 in the Y direction to generate the divided image 55, and sets the pixel value of the pixel E62 in the average image 62 to Y.
  • a division image 56 is generated by performing division processing in the Y direction for division in the direction.
  • the image processing unit 3b generates an average image 71 (70) obtained by averaging the pixel values of the corresponding pixels E of the divided image 55 and the divided image 56. That is, the image processing unit 3b averages the pixel value of the pixel E55 of the divided image 55 and the pixel value of the pixel E56 of the divided image 56 (corresponding to the pixel E55) to obtain the pixel E71 of the average image 71.
  • an average image 71 (70) having a higher resolution (doubled) in the Y direction than the average image 61 (60) and the average image 62 (60) is generated. To be done.
  • the length L5 of one pixel Ec of the average image 71 (70) is half the length L4 of one pixel Eb of the average image 61 (60) and the average image 62 (60).
  • the average image 71 (70) is closer to the X direction and Y than the X-ray fluoroscopic image Ia (the first X-ray fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40).
  • the image processing unit 3b reconstructs the super-resolution image Ic (which has been subjected to super-resolution processing in the X direction and the Y direction) to generate a reconstructed image 70x.
  • the image processing unit 3b generates the super-resolution image Ic for the plurality of X-ray fluoroscopic images Ia acquired by tomography using the rotary stage 4.
  • the image processing unit 3b reconstructs the (two-dimensional) average images 71, 72, 73, 74,... Which are the super-resolution images Ic to generate a three-dimensional reconstructed image 70x.
  • the average image 72 is the super-resolution image Ic generated by the first X-ray fluoroscopic image 12, the second X-ray fluoroscopic image 22, the third X-ray fluoroscopic image 32, and the fourth X-ray fluoroscopic image 42, which is the average image 72 (70). ).
  • the average image 73 is the average image 73 (70) that is the super-resolution image Ic generated by the first X-ray fluoroscopic image 13, the second X-ray fluoroscopic image 23, the third X-ray fluoroscopic image 33, and the fourth X-ray fluoroscopic image 43. is there.
  • the average image 74 is the average image 74 (70) that is the super-resolution image Ic generated by the first X-ray fluoroscopic image 14, the second X-ray fluoroscopic image 24, the third X-ray fluoroscopic image 34, and the fourth X-ray fluoroscopic image 44. is there.
  • step 110 the image processing unit 3b performs the super-resolution processing Ib in the X direction on the basis of the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20. To generate.
  • step 111 the image processing unit 3b converts the pixel E1 of the first X-ray fluoroscopic image 10 into the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (X Based on the pixel values of two pixels E2 that overlap the pixel E1 in the second X-ray fluoroscopic image 20 when they are overlapped and displayed in the X direction by an amount corresponding to the movement amount (of the detection position of the line). , In the X direction to generate a divided image 51.
  • step 112 the image processing unit 3b sets the pixel E2 of the second X-ray fluoroscopic image 20 to the movement amount (of the X-ray detection position) between the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20.
  • the divided image is divided in the X direction based on the pixel values of two pixels E1 overlapping the pixel E2 in the first X-ray fluoroscopic image 10. 52 is generated.
  • the order of step 111 and step 112 may be reversed.
  • step 113 the image processing unit 3b averages the divided image 51 and the divided image 52 between the pixels E corresponding to each other to generate an average image 61 (super-resolution processed by super-resolution processing in the X direction). Image Ib) is generated.
  • step 120 the image processing unit 3b performs super-resolution processing in the X direction based on the third X-ray fluoroscopic image 30 and the fourth X-ray fluoroscopic image 40. Generate Ib.
  • the order of step 110 and step 120 may be reversed.
  • step 121 the image processing unit 3b sets the pixel E3 of the third X-ray fluoroscopic image 30 to the third X-ray fluoroscopic image 30 and the fourth X-ray fluoroscopic image 40 (X Based on the pixel values of two pixels E4 that overlap with the pixel E3 in the fourth X-ray fluoroscopic image 40 when they are displayed in a superimposed manner by being shifted in the X direction by the amount corresponding to the movement amount (of the detection position of the line). , In the X direction to generate a divided image 53.
  • step 122 the image processing unit 3b sets the pixel E4 of the fourth X-ray fluoroscopic image 40 to the movement amount (of the X-ray detection position) between the third X-ray fluoroscopic image 30 and the fourth X-ray fluoroscopic image 40.
  • the divided image is divided in the X direction based on the pixel values of two pixels E3 that overlap the pixel E4 in the third X-ray fluoroscopic image 30. 54 is generated.
  • the order of steps 121 and 122 may be reversed.
  • step 123 the image processing unit 3b averages the divided image 53 and the divided image 54 between the pixels E corresponding to each other to generate an average image 62 (super-resolution processed by super-resolution processing in the X direction). Image Ib) is generated.
  • step 130 the image processing unit 3b, based on the average image 61 and the average image 62 (which is the super-resolution image Ib subjected to the super-resolution processing in the X direction). , And a super-resolution image Ic that has undergone super-resolution processing in the Y direction is generated.
  • step 131 the image processing unit 3b shifts the average image 61 and the average image 62 in the Y direction by an amount corresponding to the movement amount (of the X-ray detection position).
  • the pixel E61 of the average image 61 is divided in the Y direction based on the pixel values of two pixels E62 overlapping the pixel E61 in the average image 62 to generate a divided image 55. To do.
  • step 132 the image processing unit 3b shifts the average image 61 and the average image 62 in the Y direction by the amount corresponding to the movement amount (of the X-ray detection position), and displays them in an averaged manner.
  • the pixel E62 of the image 62 is divided in the Y direction based on the pixel values of the two pixels E61 overlapping the pixel E62 in the average image 61 to generate the divided image 56.
  • the order of step 131 and step 132 may be reversed.
  • step 133 the image processing unit 3b averages the divided images 55 and 56 between the pixels E corresponding to each other, and the average image 70 (super-resolution processing is performed in the X direction and the Y direction).
  • a super-resolution image Ic) is generated.
  • the image processing unit 3b uses the pixel value of the pixel E1 of the first X-ray fluoroscopic image 10, the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (X-ray detection). In the case of overlapping and displaying by shifting in the X direction by an amount corresponding to the movement amount (position), in the X direction based on the pixel values of two pixels E2 overlapping with the pixel E1 in the second X-ray fluoroscopic image 20.
  • a super-resolution image Ib having a higher resolution in the X direction than the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 is generated based on the divided image 50 (divided image 51) that has been divided in the X direction.
  • the pixel value of the pixel E1 can be divided in the X direction in consideration of the distribution of the actual pixel value by a relatively simple operation of dividing the pixel value of the pixel E1 based on the pixel values of the two pixels E2. ..
  • the calculation is performed only by dividing the pixel value based on the pixel values of the two pixels E2, it is not necessary to use the sequential calculation and the large number of parameters that are performed in the conventional super-resolution processing.
  • the pixel E1 is divided in consideration of the distribution of the actual pixel value E2, at least the resolution in the X direction is higher than that of the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20, and a certain image quality is secured.
  • the super-resolution image Ib can be generated. As a result, it is possible to prevent the calculation time when generating the super-resolution image Ib from increasing and to maintain a certain super-resolution effect (an effect of increasing the resolution without degrading the image quality).
  • the image Ib can be generated.
  • the image processing unit 3b uses the pixel value of the pixel E1 of the first X-ray fluoroscopic image 10 to determine the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (X-ray fluoroscopic image).
  • the super-resolution image Ib is generated based on the divided image 50 (divided image 51) that has been divided based on the area ratio.
  • the pixel E1 is added to the pixel values of the two pixels E2 that overlap with the pixel E1, and the overlapping area ratio of the two pixels E2 is taken into consideration. It can be divided in the X direction by considering it more accurately. As a result, as compared with the case where the overlapping area ratio of the two pixels E2 is not considered, the resolution in at least the X direction is higher than that of the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20, and the image quality is superior.
  • the resolution image Ib can be generated.
  • the image processing unit 3b uses the pixel value of the pixel E1 in the first X-ray fluoroscopic image 10 to determine the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 (X-ray fluoroscopic image).
  • the pixels are shifted in the X direction by an amount corresponding to the amount of movement of (detection position of) and are displayed in an overlapping manner, division is performed based on the pixel values of two pixels E2 that overlap the pixel E1 in the second X-ray fluoroscopic image 20.
  • the pixel value of the pixel E2 in the processed divided image 51 and the second X-ray fluoroscopic image 20 corresponds to the movement amount (of the X-ray detection position) between the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20.
  • a divided image 52 obtained by performing division processing based on the pixel values of two pixels E1 that overlap the pixel E2 in the first X-ray fluoroscopic image 10 when they are overlapped and displayed in the X direction. Based on this, the super-resolution image Ib is generated.
  • the super-resolution image Ib is considered by considering both the divided image 51 and the divided image 52.
  • the generation can improve the accuracy of dividing the pixel value in the X direction.
  • the super-resolution image Ib is generated by considering only one of the divided image 51 and the divided image 52, at least in the X direction than the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20. It is possible to generate a super-resolution image Ib having high resolution and high image quality.
  • the image processing unit 3b generates the average image 61 obtained by averaging the divided images 51 and 52 in the pixels E corresponding to each other, thereby generating the super-resolution image Ib. To generate.
  • the error generated between the divided images 51 and 52 can be reduced as compared with the case where the average image 61 is not generated, and thus the accuracy of dividing pixel values in the X direction can be easily improved. be able to.
  • the detector 2 moves from the first position P1 and the second position P2 by less than one pixel each.
  • the detection is performed at the third position P3 and the fourth position P4 that are translated in the Y direction orthogonal to the X direction.
  • the image processing unit 3b is configured to perform the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image, respectively, based on the X-rays detected at the first position P1, the second position P2, the third position P3, and the fourth position P4. It is configured to generate the image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40.
  • the image processing unit 3b performs division processing in the X direction on the basis of the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20, divided images 50 (51, 52), the third X-ray fluoroscopic image 30, and the third X-ray fluoroscopic image 30.
  • the same processing as the division processing in the X direction is performed in the Y direction, and thus the first X
  • the super-resolution image Ic having a higher resolution in the X direction and the Y direction than the fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40 is configured to be generated. This makes it possible to generate a super-resolution image Ic that has a high resolution not only in the X direction but also in the Y direction and has a certain image quality.
  • the resolution is high only in the X direction
  • the pixels E in a general display device have the same size in the row direction (X direction) and the column direction (Y direction) of the detection elements 2a arranged in a matrix
  • the resolutions in the X direction and the Y direction are equal.
  • it is possible to simplify the post-process for displaying the super-resolution image Ic such as the interpolation processing of the image for adjusting the pixel size in the X direction and the Y direction on the display device (not shown). ..
  • the image processing unit 3b causes the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 to correspond to the movement amount (of the X-ray detection position) in the X direction.
  • the two X-ray fluoroscopic images 20 are overlapped with each other and overlapped with each other, if there is only one pixel E2 that overlaps the pixel E1 of the first X-ray fluoroscopic image 10, the pixels overlap each other. Based on the pixel value of one pixel E2, another overlapping pixel E2 is virtually generated.
  • the detector 2 receives the X-rays emitted from the X-ray source 1 at the first position P1 (third position P3) and the first position P1 (third position P3).
  • the second position P2 (fourth position P4), which is translated in the X direction by a length L5 that is half the length of one pixel Ea (the length L4 thereof).
  • the overlapping area ratios of the two pixels E2 (pixel E4) become equal (one-to-one), so that the pixel value of the pixel E1 (pixel E3) exceeds that of the pixel E1 (pixel E3).
  • the detector 2 detects the X-rays emitted from the X-ray source 1 from a plurality of directions while rotating the X-rays about the X-axis as a rotation axis, and the detection operation in the X-direction.
  • the tomography is performed by alternately repeating the parallel movement and the parallel movement.
  • the image processing unit 3b is configured to generate the virtual pixel E so as to have a pixel value equal to the pixel value of the adjacent pixel E, but the present invention is not limited to this. It is not limited to this.
  • the image processing unit may be configured to generate a pixel value of a virtual pixel so as to have a pixel value different from that of an adjacent pixel.
  • the X-ray detection position is moved in parallel by the length L5 which is half the length L4 of one pixel Ea, and the image processing unit 3b is operated by the first X-ray fluoroscopic image 10 (third X-ray image).
  • the X-ray detection position may be moved in parallel by a length longer than half the length of one pixel, or a position shorter than half the length of one pixel. You may comprise so that it may be moved in parallel.
  • the image processing unit sets the pixel value of “the first pixel (of the claims)” of “the first image of the claims” to “the second image of the claims”. So as to generate a divided image based on the pixel value of two "(second claim)" pixels that overlap with the "first pixel” and the overlapping area ratio. Just configure it.
  • the X-ray fluoroscopic images Ia based on the X-rays detected at the position where the image processing unit 3b is moved in parallel in either the X direction or the Y direction are divided, and the divided image 50 is obtained.
  • the X-ray fluoroscopic images based on the X-rays detected at the position where the image processing unit is translated in both the X direction and the Y direction May be divided to generate divided images.
  • the image processing unit 3b is configured to generate the super-resolution image Ib by generating the average image 60 obtained by averaging the two divided images 50 in the pixels E corresponding to each other.
  • the present invention is not limited to this.
  • the image processing unit may be configured to process one divided image itself as a super-resolution image.
  • the X-ray imaging apparatus 100 performs scan imaging at the first position P1, the second position P2, the third position P3, and the fourth position P4 prior to the generation of the super-resolution image Ic.
  • the first X-ray fluoroscopic image 10, the second X-ray fluoroscopic image 20, the third X-ray fluoroscopic image 30, and the fourth X-ray fluoroscopic image 40 are configured to be generated as described above, but the present invention is not limited to this.
  • the X-ray fluoroscopic image necessary for generating one super-resolution image is generated and before the X-ray fluoroscopic image necessary for generating another super-resolution image is generated, It may be configured to start generation of the super-resolution image.
  • the third X-ray fluoroscopic image 30 and the fourth X-ray fluoroscopic image 40 are not generated at the time when the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 are generated.
  • the generation of the super-resolution image Ib that has been subjected to the super-resolution processing in the X direction based on the first X-ray fluoroscopic image 10 and the second X-ray fluoroscopic image 20 may be started.
  • the image processing unit 3b is configured to generate the super-resolution image Ib based on the two X-ray fluoroscopic images Ia.
  • the present invention is not limited to this. ..
  • the image processing unit may be configured to generate a super-resolution image based on three or more X-ray fluoroscopic images.
  • the image processing unit 3b is configured to reconstruct the super-resolution image Ic that has undergone the super-resolution processing in the X direction and the Y direction to generate the reconstructed image 70x.
  • the present invention is not limited to this.
  • the image processing unit may be configured to reconstruct a super-resolution image that has undergone super-resolution processing in only the X direction or the Y direction, and generate a reconstructed image. .. In this case, X-ray fluoroscopic images with different X-ray detection positions in the direction in which the super-resolution processing is not necessary are unnecessary.
  • the X-ray imaging apparatus 100 that performs non-helical scan tomography is configured to generate the super-resolution images Ib and Ic, but the present invention is not limited to this. Absent.
  • a super-resolution image may be generated in an X-ray imaging apparatus that performs helical scan tomography (which operates in a spiral shape).
  • an X-ray imaging apparatus that does not perform tomography may be configured to generate a super-resolution image.
  • Detector moving mechanism (moving mechanism) 10 (11, 12, 13, 14) 1st X-ray fluoroscopic image (1st image) 20 (21, 22, 23, 24) Second X-ray fluoroscopic image (second image) 30 (31, 32, 33, 34) Third X-ray fluoroscopic image (third image) 40 (41, 42, 43, 44) 4th X-ray fluoroscopic image (4th image) 50 divided images 51, 53 divided images (first divided image) 52, 54 split image (second split image) 60 (61, 62) Average image 90 Rotation axis 100 X-ray imaging device E (Ea, Eb, Ec) Pixel E1, E3 Pixel (first pixel) E2 and E4 pixels (second pixel) Ib, Ic Super-resolution image P1 1st position P2 2nd position P3 3rd position P4 4th position

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Pulmonology (AREA)
  • Immunology (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Optics & Photonics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Image Processing (AREA)

Abstract

このＸ線撮影装置（１００）では、画像処理部（３ｂ）は、第１Ｘ線透視画像（１０）の第１画素（Ｅ１）の画素値を、第１Ｘ線透視画像と第２Ｘ線透視画像とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像（２０）における第１画素（Ｅ１）に対してオーバラップする２つの画素（Ｅ）の画素値に基づいてＸ方向に分割することにより、第１Ｘ線透視画像および第２Ｘ線透視画像よりもＸ方向における解像度が高い超解像画像（Ｉｂ）を生成するように構成されている。

Description

Ｘ線撮影装置
　本発明は、Ｘ線撮影装置に関し、特に、第１位置および第２位置において検出されたＸ線に基づいて生成された第１画像および第２画像よりも解像度が高い超解像画像を生成するＸ線撮影装置に関する。
　従来、第１位置および第２位置において検出されたＸ線に基づいて生成された第１画像および第２画像よりも解像度が高い超解像画像を生成するＸ線撮影装置が知られている（たとえば、特許文献１参照）。
　上記特許文献１には、Ｘ線ビームを放出するＸ線発生器と、Ｘ線ビームの放出方向に対して垂直に配置された平面を有する放射線検出器と、放射線検出器を保持するとともに、放射線検出器をピクセル方向に平行移動させる平行移動ステージと、を備えたＸ線撮影装置が開示されている。上記特許文献１のＸ線撮影装置では、放射線検出器は、平行移動ステージによりピクセルサイズよりも微細な距離だけ平行移動されることによって、互いに異なる検出位置に移動するように構成されている。そして、互いに異なる検出位置において取得された複数のラジオグラフ（画像）に基づいて、取得されたラジオグラフよりも微細な解像度を有する画像（超解像画像）を生成するように構成されている。
特表２０１６－５１７９６１号公報
　ここで、上記特許文献１には記載されていないが、上記特許文献１に記載のような従来のＸ線撮影装置のように、取得された複数の画像に基づいて超解像画像を生成する場合には、多数のパラメータを用いた逐次的な演算が行われることが知られている。このように、逐次的な演算が行われる場合、計算時間が比較的長くなり易い。また、多数のパラメータが用いられる場合、生成される超解像画像の画質がパラメータに左右され易く、一定の超解像効果（画質を低下させることなく解像度を高くする効果）を得られない場合がある。このため、上記特許文献１に記載のような従来のＸ線撮影装置では、取得された複数の画像に基づいて超解像画像を生成するための計算時間が長くなるとともに、生成された超解像画像において一定の超解像効果を得られない場合があるという問題点がある。
　この発明は、上記のような課題を解決するためになされたものであり、この発明の１つの目的は、超解像画像を生成する際の計算時間が長くなるのを抑制するとともに一定の超解像効果を確保しながら超解像画像を生成することが可能なＸ線撮影装置を提供することである。
　上記目的を達成するために、この発明の一の局面におけるＸ線撮影装置は、Ｘ線源と、Ｘ線源から照射されたＸ線を、第１位置と、第１位置から１画素分よりも小さい移動量だけ第１方向へ平行移動させた第２位置とにおいて検出する検出器と、第１位置において検出されたＸ線に基づいて第１画像を生成するとともに、第２位置において検出されたＸ線に基づいて第２画像を生成する画像処理部と、検出器によるＸ線の検出位置を、第１位置と第２位置との間で移動させる移動機構と、を備え、画像処理部は、第１画像と第２画像とを上記移動量に対応する分だけ第１方向にずらして重ねて表示した場合に、第１画像および第２画像の一方の第１画素の画素値を、第１画像および第２画像の他方における第１画素に対してオーバラップする２つの画素の画素値に基づいて第１方向に分割する第１方向の分割処理を行った分割画像に基づいて、第１画像および第２画像よりも第１方向における解像度が高い超解像画像を生成するように構成されている。
　この発明の一の局面によるＸ線撮影装置では、上記のように、画像処理部は、第１画像および第２画像の一方の第１画素の画素値を、第１画像と第２画像とを上記移動量に対応する分だけ第１方向にずらして重ねて表示した場合に、第１画像および第２画像の他方における第１画素に対してオーバラップする２つの画素の画素値に基づいて第１方向に分割する第１方向の分割処理を行った分割画像に基づいて、第１画像および第２画像よりも第１方向における解像度が高い超解像画像を生成するように構成されている。ここで、第１画像と第２画像とを上記移動量に対応する分だけ第１方向にずらして重ねて表示した場合に、第１画像および第２画像の他方における第１画素に対してオーバラップする２つの画素は、第１方向において、それぞれ、第１画素の一方側および他方側に部分的にずれた画素である。すなわち、第１画像および第２画像の他方における第１画素に対してオーバラップする２つの画素は、第１画素の画素値を、第１方向の一方側と他方側とに分割する場合の目安の値と考えることができる。したがって、上記のように構成すれば、第１画素の画素値を、２つの画素の画素値に基づいて分割するだけの比較的単純な演算により、実際の画素値の分布を考慮して、第１方向に分割することができる。すなわち、２つの画素の画素値に基づいて画素値を分割するだけの演算なので、従来の超解像処理で行っていた逐次的な演算や多数のパラメータを用いる必要がない。また、実際の画素値の分布を考慮して画素を分割するので、第１画像および第２画像よりも少なくとも第１方向における解像度が高いとともに一定の画質が確保された超解像画像を生成することができる。その結果、超解像画像を生成する際の計算時間が長くなるのを抑制するとともに一定の超解像効果（画質を低下させることなく解像度を高くする効果）を確保しながら、超解像画像を生成することができる。
　上記一の局面によるＸ線撮影装置において、好ましくは、画像処理部は、第１画像および第２画像の一方の第１画素の画素値を、第１画像と第２画像とを上記移動量に対応する分だけ第１方向にずらして重ねて表示した場合に、第１画像および第２画像の他方における第１画素に対してオーバラップする２つの画素の画素値およびオーバラップする面積比率に基づいて分割処理を行った分割画像に基づいて、超解像画像を生成するように構成されている。このように構成すれば、第１画素を、第１画素に対してオーバラップする２つの画素の画素値に加えて、２つの画素のオーバラップする面積比率が考慮されているので、実際の画素値の分布をより精度よく考慮して、第１方向に分割することができる。その結果、２つの画素のオーバラップする面積比率を考慮しない場合と比較して、第１画像および第２画像よりも少なくとも第１方向における解像度が高いとともにより画質の良い超解像画像を生成することができる。
　上記一の局面によるＸ線撮影装置において、好ましくは、画像処理部は、第１画像における第１画素の画素値を、第１画像と第２画像とを上記移動量に対応する分だけ第１方向にずらして重ねて表示した場合に、第２画像における第１画素に対してオーバラップする２つの画素の画素値に基づいて分割処理を行った第１分割画像と、第２画像における第２画素の画素値を、第１画像と第２画像とを上記移動量に対応する分だけ第１方向にずらして重ねて表示した場合に、第１画像における第２画素に対してオーバラップする２つの画素の画素値に基づいて分割処理を行った第２分割画像とに基づいて、超解像画像を生成するように構成されている。このように構成すれば、第１分割画像または第２分割画像の一方のみを考慮して超解像画像を生成する場合と比較して、第１分割画像および第２分割画像の両方を考慮して超解像画像を生成することにより、第１方向における画素値の分割の精度を向上することができる。その結果、第１分割画像または第２分割画像の一方のみを考慮して超解像画像を生成する場合と比較して、第１画像および第２画像よりも少なくとも第１方向における解像度が高いとともにより画質の良い超解像画像を生成することができる。
　この場合、好ましくは、画像処理部は、第１分割画像と第２分割画像とを、互いに対応する画素において平均した平均画像を生成することにより、超解像画像を生成するように構成されている。このように構成すれば、平均画像を生成しない場合と比較して、第１分割画像と第２分割画像との間に生じる誤差を低減することができるので、第１方向における画素値の分割の精度を容易に向上することができる。
　上記一の局面によるＸ線撮影装置において、好ましくは、検出器は、第１位置および第２位置に加えて、第１位置および第２位置からそれぞれ１画素分よりも小さい移動量だけ第１方向に直交する第２方向へ平行移動させた第３位置および第４位置とにおいて検出するように構成されており、画像処理部は、第１位置、第２位置、第３位置および第４位置において検出されたＸ線に基づいて、それぞれ、第１画像、第２画像、第３画像および第４画像を生成するとともに、第１画像および第２画像に基づいて第１方向の分割処理を行った分割画像と、第３画像および第４画像に基づいて第１方向の分割処理を行った分割画像とに基づいて、第１方向における分割処理と同様の処理を第２方向に対して行うことにより、第１画像、第２画像、第３画像および第４画像よりも第１方向および第２方向における解像度が高い超解像画像を生成するように構成されている。このように構成すれば、第１方向に加えて第２方向においても解像度が高いとともに一定の画質が確保された超解像画像を生成することができるので、第１方向のみにおいて解像度が高い場合と比較して、第１方向と第２方向とにおいて、解像度の程度に偏りのない超解像画像を生成することができる。また、一般的な表示装置における画素は、マトリクス状に配列された検出素子の行方向および列方向におけるサイズが等しいので、第１方向と第２方向との解像度が異なる場合のように、第１方向と第２方向とにおける画素サイズを調整する画像の補間処理等の超解像画像を表示装置に表示させるための後工程を簡略化することができる。
　上記一の局面によるＸ線撮影装置において、好ましくは、画像処理部は、第１画像または第２画像の他方の端部において、第１画像と第２画像とを上記移動量に対応する分だけ第１方向にずらして重ねて表示した場合に、第１画像および第２画像の一方の第１画素に対してオーバラップする画素が１つしかない場合には、オーバラップする１つの画素の画素値に基づいて、もう１つのオーバラップする画素を仮想的に生成するように構成されている。このように構成すれば、もう１つのオーバラップする画素を仮想的に生成することにより、端部に位置する第１画素に対してオーバラップする画素が１つしかない場合でも、第１画素を第１方向に分割する演算を確実に行うことができる。
　上記画像処理部が第１画素の画素値を第１画素に対してオーバラップする２つの画素の画素値およびオーバラップする面積比率に基づいて分割処理を行った分割画像に基づいて超解像画像を生成する構成において、好ましくは、検出器は、Ｘ線源から照射されたＸ線を、第１位置と、第１位置から１画素分の半分の長さ分だけ第１方向へ平行移動させた第２位置とにおいて検出するように構成されている。このように構成すれば、２つの画素のオーバラップする面積比率が等しくなる（１対１になる）ので、第１画素の画素値を、第１画素に対してオーバラップする２つの画素の画素値自体の比率により、第１方向に分割することができる。その結果、２つの画素のオーバラップする面積比率を２分割以外の比率に分割する場合と比較して、分割処理における演算を簡略化することができるので、超解像画像を生成する際の計算時間が長くなるのをより抑制することができる。
　上記一の局面によるＸ線撮影装置において、好ましくは、検出器は、Ｘ線源から照射されたＸ線を、第３方向を回転軸として回転しながら複数の方向から検出する検出動作と、第３方向への平行移動とを、交互に繰り返し行うことにより、断層撮影を行うように構成されている。このように構成すれば、上記のような（いわゆるノンヘリカルスキャン型の）断層撮影を行うＸ線撮影装置において、超解像画像を生成する際の計算時間が長くなるのを抑制するとともに一定の超解像効果を確保しながら超解像画像を生成することができる。
　本発明によれば、上記のように、超解像画像を生成する際の計算時間が長くなるのを抑制するとともに一定の超解像効果を確保しながら超解像画像を生成することができる。
本発明の一実施形態によるＸ線撮影装置の全体構成を示したＹ方向から見た図である。 本発明の一実施形態によるＸ線撮影装置の全体構成を示したＸ方向から見た図である。 本発明の一実施形態によるＸ線撮影装置におけるＸ線の検出位置を説明するための図である。 本発明の一実施形態によるＸ線撮影装置により撮影されたＸ線透視画像を説明するための図である。 本発明の一実施形態によるＸ線撮影装置における超解像処理を説明するための図である。 Ｘ方向の超解像処理が行われる第１Ｘ線透視画像の一部および第２Ｘ線透視画像の一部を示した図である。 第１Ｘ線透視画像および第２Ｘ線透視画像に基づいて生成された第１分割画像および第２分割画像を説明するための図である。 第１分割画像と第２分割画像とを平均した平均画像を説明するための図である。 第３Ｘ線透視画像および第４Ｘ線透視画像に対するＸ方向の超解像処理を説明するための図である。 Ｘ方向の超解像処理が行われた２つの画像に基づくＹ方向の超解像処理を説明するための図である。 Ｘ方向およびＹ方向の超解像処理が行われた複数の超解像画像を再構成した再構成画像の生成を説明するための図である。 超解像画像の生成フローである。 第１Ｘ線透視画像および第２Ｘ線透視画像に基づいて、Ｘ方向の超解像処理を行うフローである。 第３Ｘ線透視画像および第４Ｘ線透視画像に基づいて、Ｘ方向の超解像処理を行うフローである。 Ｘ方向の超解像処理が行われた２つの画像に基づいて、さらにＹ方向の超解像処理を行うフローである。
　以下、本発明を具体化した実施形態を図面に基づいて説明する。
　まず、図１～図９を参照して、本発明の一実施形態によるＸ線撮影装置１００の構成について説明する。
　図１に示すように、Ｘ線撮影装置１００は、Ｘ線源１と、検出器２と、処理ユニット３と、回転ステージ４と、検出器移動機構５と、を備えている。なお、検出器移動機構５は、特許請求の範囲の「移動機構」の一例である。
　Ｘ線撮影装置１００では、Ｘ線源１と、回転ステージ４と、検出器２とが、Ｘ線源１と検出器２とを結ぶ方向（Ｚ方向）に、この順に並んで配置されている。本明細書では、Ｘ線源１から検出器２に向かう方向をＺ２方向、その逆方向をＺ１方向とする。また、Ｚ方向と直交する面内方向のうちの回転ステージ４が並進移動する方向をＹ方向とする。また、Ｚ方向およびＹ方向と直交する方向をＸ方向とする。なお、Ｘ方向は、特許請求の範囲の「第１方向」および「行方向」の一例である。また、Ｙ方向は、特許請求の範囲の「第２方向」および「列方向」の一例である。
　Ｘ線源１は、高電圧が印加されることにより、Ｘ線を発生させることが可能なＸ線発生装置である。Ｘ線源１は、発生させたＸ線をＺ２方向に照射するように構成されている。
　検出器２は、Ｘ線源１から照射されたＸ線を検出するとともに、検出されたＸ線を電気信号に変換する。検出器２は、マトリクス状にＸ方向およびＹ方向に所定の周期Ｌ１で並んで配列された複数の検出素子２ａを含む検出面２ｂを有する。検出器２は、複数の検出素子２ａのそれぞれに対応するように配置されるとともに、検出されたＸ線を電気信号に変化する複数の変換素子（図示せず）を含む。検出器２は、たとえば、ＦＰＤ（Ｆｌａｔ　Ｐａｎｅｌ　Ｄｅｔｅｃｔｏｒ）である。検出器２で変換された検出信号（電気信号）は、処理ユニット３が備える画像処理部３ｂ（後述する）に送られる。
　処理ユニット３は、制御部３ａと、画像処理部３ｂと、を備えている。
　制御部３ａは、回転ステージ４および検出器移動機構５の動作を制御するように構成されている。制御部３ａは、たとえば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）およびＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）などを含む。
　画像処理部３ｂは、検出器２から送られた検出信号に基づいて、Ｘ線透視画像Ｉａ（図４参照）を生成するように構成されている。また、画像処理部３ｂは、複数のＸ線透視画像Ｉａ（図４参照）に基づいて、後述する再構成画像７０ｘ（図８参照）を生成するように構成されている。画像処理部３ｂは、たとえば、ＧＰＵ（Ｇｒａｐｈｉｃｓ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）や画像処理用に構成されたＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ　Ｇａｔｅ　Ａｒｒａｙ）などのプロセッサを含む。なお、Ｘ線透視画像Ｉａ（図４参照）における各画素Ｅａ（Ｅ）のデータは、それぞれ、検出器２の各検出素子２ａで検出されたＸ線量に対応する。
　回転ステージ４は、被写体Ｓを載置する載置面４ａ（図２参照）を有する。回転ステージ４は、ＸＺ平面内において、Ｙ方向に沿った回転軸９０周りに３６０度回転移動することが可能に構成されている。また、図２に示すように、回転ステージ４は、Ｙ方向に並進移動することが可能に構成されている。これにより、本実施形態では、検出器２は、Ｘ線源１から照射されたＸ線を、Ｙ方向を回転軸として回転しながら（ＸＺ平面内の）複数の方向から検出する検出動作と、Ｙ方向への平行移動とを、交互に繰り返し行うことにより、断層撮影を行うことが可能に構成されている。すなわち、Ｘ線撮影装置１００は、いわゆるノンヘリカルスキャン型の断層撮影を行うことが可能な撮影装置である。
　図１および図２に示すように、検出器移動機構５は、検出器２を、検出素子２ａよりも小さい移動量だけＸ方向およびＹ方向にそれぞれ移動させることが可能に構成されている。これにより、図３に示すように、検出器移動機構５は、検出器２を、第１位置Ｐ１と、第２位置Ｐ２と、第３位置Ｐ３と、第４位置Ｐ４との間で移動させることが可能である。第２位置Ｐ２は、第１位置Ｐ１から距離Ｌ２だけＸ方向へ平行移動させた位置である。第３位置Ｐ３は、第１位置Ｐ１から距離Ｌ３だけＹ方向へ平行移動させた位置である。第４位置は、第１位置Ｐ１から距離Ｌ３だけＹ方向における第３位置と同じ側へ平行移動させた位置である。距離Ｌ２および距離Ｌ３は、それぞれ、Ｘ方向およびＹ方向において、検出素子２ａの１つ分の大きさ（検出素子２ａが並ぶ周期）Ｌ１よりも小さい距離である。なお、Ｘ線撮影装置１００では、距離Ｌ２および距離Ｌ３は、それぞれ、Ｘ方向およびＹ方向において、検出素子２ａの１つ分の大きさＬ１の半分の距離（長さ）である。すなわち、本実施形態では、検出器２は、Ｘ線源１から照射されたＸ線を、第１位置Ｐ１と、第１位置Ｐ１から１画素Ｅａ分（の長さＬ４（図５参照））の半分の長さＬ５（図５参照）分だけＸ方向へ平行移動させた第２位置Ｐ２と、第１位置Ｐ１および第２位置Ｐ２からそれぞれ１画素Ｅａ分（の長さＬ４（図５参照））の半分の長さＬ５（図５参照）分だけＹ方向へ平行移動させた第３位置Ｐ３および第４位置Ｐ４とにおいて検出することが可能に構成されている。
　以上の構成により、図４に示すように、画像処理部３ｂは、第１位置Ｐ１、第２位置Ｐ２、第３位置Ｐ３および第４位置Ｐ４において検出されたＸ線に基づいて、それぞれ、Ｘ線透視画像Ｉａとして、第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０を生成することが可能である。なお、第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０は、それぞれ、特許請求の範囲の「第１画像」、「第２画像」、「第３画像」および「第４画像」の一例である。
　第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０は、それぞれ、回転ステージ４を用いた断層撮影により取得された複数のＸ線透視画像Ｉａを含む。すなわち、第１Ｘ線透視画像１０は、回転ステージ４によりＸＺ平面内の撮影方向およびＹ方向の撮影位置の少なくとも一方が互いに異なる第１Ｘ線透視画像１１、１２、１３、１４、…を含む。第２Ｘ線透視画像２０は、回転ステージ４の動作によりＸＺ平面内の撮影方向およびＹ方向の撮影位置の少なくとも一方が互いに異なる第２Ｘ線透視画像２１、２２、２３、２４、…を含む。第３Ｘ線透視画像３０は、回転ステージ４の動作によりＸＺ平面内の撮影方向およびＹ方向の撮影位置の少なくとも一方が互いに異なる第３Ｘ線透視画像３１、３２、３３、３４、…を含む。第４Ｘ線透視画像４０は、回転ステージ４の動作によりＸＺ平面内の撮影方向およびＹ方向の撮影位置の少なくとも一方が互いに異なる第４Ｘ線透視画像４１、４２、４３、４４、…を含む。
　図５に示すように、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とは、Ｘ方向において、１画素Ｅａの大きさＬ４の半分の大きさＬ５分だけ撮影位置が異なるデータである。また、第３Ｘ線透視画像３０と第４Ｘ線透視画像４０とは、Ｘ方向において、１画素Ｅａの大きさＬ４の半分の大きさＬ５分だけ撮影位置が異なるデータである。また、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０と、第３Ｘ線透視画像３０と第４Ｘ線透視画像４０とは、Ｙ方向において、１画素Ｅａの大きさＬ４の半分の大きさＬ５分だけ撮影位置が異なるデータである。なお、図５では、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示するとともに、第３Ｘ線透視画像３０と第４Ｘ線透視画像４０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した状態を示している。
　ここで、本実施形態では、画像処理部３ｂは、第１Ｘ線透視画像１０（Ｘ線透視画像Ｉａ）と第２Ｘ線透視画像２０（Ｘ線透視画像Ｉａ）とに基づいて、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりもＸ方向における解像度が高い超解像画像Ｉｂを生成するように構成されている。なお、以下の説明では、複数の画像（元画像）に基づいて、元画像よりも解像度が高い画像を生成する処理を「超解像処理」と呼ぶ場合がある。
　詳細には、図６に示すように、画像処理部３ｂは、第１Ｘ線透視画像１０における画素Ｅ１（図６Ａ参照）の画素値をＸ方向に分割するＸ方向の分割処理を行った分割画像５１（図６Ｂ参照）と、第２Ｘ線透視画像２０における画素Ｅ２（図６Ａ参照）の画素値をＸ方向に分割するＸ方向の分割処理を行った分割画像５２（図６Ｂ参照）とに基づいて、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりもＸ方向における解像度が高い超解像画像Ｉｂ（図５参照）を生成するように構成されている。また、画像処理部３ｂは、分割画像５１（図６Ｂ参照）と分割画像５２（図６Ｂ参照）とを、互いに対応する画素Ｅにおいて平均した平均画像６０（図６Ｃ参照）を生成することにより、超解像画像Ｉｂ（図５参照）を生成するように構成されている。なお、画素Ｅ１および画素Ｅ２は、それぞれ、特許請求の範囲の「第１画素」および「第２画素」の一例である。また、分割画像５１および分割画像５２は、それぞれ、特許請求の範囲の「第１分割画像」および「第２分割画像」の一例である。
　分割画像５１（図６Ｂ参照）は、第１Ｘ線透視画像１０における画素Ｅ１（図６Ａ参照）の画素値を、（図５に示すように）第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像２０における画素Ｅ１（図６Ａ参照）に対してオーバラップする２つの画素Ｅ２の画素値およびオーバラップする面積比率に基づいて、Ｘ方向の分割処理を行った分割画像５０（図６Ｂ参照）である。分割画像５２（図６Ｂ参照）は、第２Ｘ線透視画像２０における画素Ｅ２（図６Ａ参照）の画素値を、第１Ｘ線透視画像１０における画素Ｅ２（図６Ａ参照）に対してオーバラップする２つの画素Ｅ１の画素値およびオーバラップする面積比率に基づいて、Ｘ方向の分割処理を行った分割画像５０（図６Ｂ参照）である。なお、図６では、図示の都合上、Ｘ線透視画像Ｉａの内の１つの行だけを取り出して説明している。また、図６Ａおよび図６Ｂでは、図５と同様に、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして表示するとともに、（重ねて表示した）図５とは異なり、図示の都合上、紙面の上下方向にずらして表示している。
　具体的には、図６Ｂに示すように、画像処理部３ｂは、第１Ｘ線透視画像１０における画素Ｅ１ｂを、分割画像５１の画素となる画素Ｅ５１ｃと画素Ｅ５１ｄとに分割する。画像処理部３ｂは、画素Ｅ５１ｃおよび画素Ｅ５１ｄの画素値を、それぞれ、第２Ｘ線透視画像２０において画素Ｅ１ｂに対してオーバラップする画素Ｅ２ａの画素値と画素Ｅ２ｂの画素値との比率で、画素Ｅ１ｂの画素値を分割した値とする。たとえば、画像処理部３ｂは、画素Ｅ１ｂの画素値（７）を、画素Ｅ２ａの画素値（５）と画素Ｅ２ｂの画素値（４）との比率で分割した値（７×（５／（５＋４））＝３．９）および値（７×（４／（５＋４））＝３．１）を、それぞれ、画素Ｅ５１ｃの画素値および画素Ｅ５１ｄの画素値とする。すなわち、Ｘ線撮影装置１００では、画像処理部３ｂは、第１Ｘ線透視画像１０の画素Ｅ１の画素値を、第２Ｘ線透視画像２０における画素Ｅ１に対してそれぞれ半分ずつがオーバラップする２つの画素Ｅ２の画素値に基づいて半分に分割するように構成されている。
　画像処理部３ｂは、同様の処理を、第１Ｘ線透視画像１０における他の画素Ｅ１についても行う。すなわち、画像処理部３ｂは、第１Ｘ線透視画像１０における画素Ｅ１ｃを、分割画像５１の画素Ｅとなる画素Ｅ５１ｅと画素Ｅ５１ｆとに分割する。画像処理部３ｂは、画素Ｅ５１ｅおよび画素Ｅ５１ｆの画素値を、それぞれ、第２Ｘ線透視画像２０において画素Ｅ１ｃに対してオーバラップする画素Ｅ２ｂの画素値と画素Ｅ２ｃの画素値との比率で、画素Ｅ１ｃの画素値を分割した値とする。また、画像処理部３ｂは、第１Ｘ線透視画像１０における画素Ｅ１ｄを、分割画像５１の画素となる画素Ｅ５１ｇと画素Ｅ５１ｈとに分割する。画像処理部３ｂは、画素Ｅ５１ｇおよび画素Ｅ５１ｈの画素値を、それぞれ、第２Ｘ線透視画像２０において画素Ｅ１ｄに対してオーバラップする画素Ｅ２ｃの画素値と画素Ｅ２ｄの画素値との比率で、画素Ｅ１ｄの画素値を分割した値とする。
　ここで、本実施形態では、画像処理部３ｂは、（図５に示すように）第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像２０の端部において、第１Ｘ線透視画像１０の画素Ｅ１に対してオーバラップする画素Ｅが１つしかない場合には、オーバラップする１つの画素Ｅ２の画素値に基づいて、もう１つのオーバラップする画素Ｅ２を仮想的に生成するように構成されている。具体的には、画像処理部３ｂは、第１Ｘ線透視画像１０における画素Ｅ１ａを、分割画像５１の画素Ｅとなる画素Ｅ５１ａと画素Ｅ５１ｂとに分割する。画像処理部３ｂは、第２Ｘ線透視画像２０において画素Ｅ１ａに対してオーバラップする（第２Ｘ線透視画像２０における端部の）画素Ｅ２ａに隣接した仮想的な画素Ｅ２ｅを生成する。そして、画像処理部３ｂは、画素Ｅ５１ａおよび画素Ｅ５１ｂの画素値を、それぞれ、画素Ｅ２ａの画素値と（仮想的な）画素Ｅ２ｄの画素値との比率で、画素Ｅ１ａの画素値を分割した値とする。なお、Ｘ線撮影装置１００では、画像処理部３ｂは、仮想的な画素Ｅの画素値を、隣接する画素Ｅの画素値と等しい画素値となるように生成するように構成されている。
　画像処理部３ｂは、分割画像５１と分割画像５２との画素数を等しくように、分割画像５１の端部の画素Ｅを生成するように構成されている。具体的には、画像処理部３ｂは、第１Ｘ線透視画像１０における（端部の）画素Ｅ１ｄに隣接した（仮想的な）画素Ｅ１ｅに基づいて、画素Ｅ５１ｉを生成する。以上の構成により、第１Ｘ線透視画像１０の画素Ｅ１をＸ方向に分割した分割画像５１が生成される。
　そして、画像処理部３ｂは、第２Ｘ線透視画像２０の画素Ｅ２をＸ方向に分割した分割画像５２も（分割画像５１の生成と）同様の方法により生成する。すなわち、画像処理部３ｂは、第２Ｘ線透視画像２０における画素Ｅ２ａを、分割画像５２の画素Ｅとなる画素Ｅ５２ｂと画素Ｅ５２ｃとに分割する。また、画像処理部３ｂは、第２Ｘ線透視画像２０における画素Ｅ２ｂを、分割画像５２の画素Ｅとなる画素Ｅ５２ｄと画素Ｅ５２ｅとに分割する。また、画像処理部３ｂは、第２Ｘ線透視画像２０における画素Ｅ２ｃを、分割画像５２の画素となる画素Ｅ５２ｆと画素Ｅ５２ｇとに分割する。また、第２Ｘ線透視画像２０における画素Ｅ２ｄを、分割画像５２の画素Ｅとなる画素Ｅ５２ｈと画素Ｅ５２ｉとに分割する。また、画像処理部３ｂは、分割画像５１と分割画像５２との画素数を等しくように、第２Ｘ線透視画像２０における（端部の）画素Ｅ２ａに隣接した（仮想的な）画素Ｅ２ｅに基づいて、画素Ｅ５２ａを生成する。以上の構成により、第２Ｘ線透視画像２０の画素Ｅ２をＸ方向に分割した分割画像５２が生成される。
　そして、図６Ｂおよび図６Ｃに示すように、画像処理部３ｂは、分割画像５１と分割画像５２との互いに対応する画素Ｅの画素値同士を平均した平均画像６１（６０）を生成する。すなわち、画像処理部３ｂは、分割画像５１の画素Ｅ５１ａの画素値（１．５）と、分割画像５２の画素Ｅ５２ａの画素値（２．５）とを平均して、平均画像６１の画素Ｅ６１ａの画素値（２）とする。画像処理部３ｂは、他の画素についても同様の処理を行う。すなわち、画像処理部３ｂは、分割画像５１の画素Ｅ５１ｂ、Ｅ５１ｃ、Ｅ５１ｄ、Ｅ５１ｅ、Ｅ５１ｇ、Ｅ５１ｈおよびＥ５１ｉの画素値と、分割画像５２の画素Ｅ５２ｂ、Ｅ５２ｃ、Ｅ５２ｄ、Ｅ５２ｅ、Ｅ５２ｇ、Ｅ５２ｈおよびＥ５２ｉの画素値とを、それぞれ平均して、平均画像６１の画素Ｅ６１ｂ、Ｅ６１ｃ、Ｅ６１ｄ、Ｅ６１ｅ、Ｅ６１ｇ、Ｅ６１ｈおよびＥ６１ｉの画素値とする。
　図５に示すように、平均画像６１（６０）は、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりもＸ方向に解像度が高くなっている（２倍になっている）。具体的には、Ｘ方向において、平均画像６１（６０）の１画素Ｅｂ分の長さＬ５は、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０の１画素Ｅａ分の長さＬ４の半分となっている。すなわち、平均画像６１（６０）は、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりもＸ方向における解像度が高い（Ｘ方向に超解像処理が行われた）超解像画像Ｉｂである。
　また、本実施形態では、画像処理部３ｂは、第１Ｘ線透視画像１０（Ｘ線透視画像Ｉａ）と、第２Ｘ線透視画像２０（Ｘ線透視画像Ｉａ）と、第３Ｘ線透視画像３０（Ｘ線透視画像Ｉａ）と、第４Ｘ線透視画像４０（Ｘ線透視画像Ｉａ）とに基づいて、第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０よりもＸ方向およびＹ方向における解像度が高い超解像画像Ｉｃを生成するように構成されている。詳細には、画像処理部３ｂは、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０に基づいてＸ方向の分割処理を行った分割画像５０（分割画像５１および分割画像５２）と、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０に基づいてＸ方向の分割処理を行った分割画像５０（分割画像５３および分割画像５４）とに基づいて、Ｘ方向における分割処理と同様の処理をＹ方向に対して行うことにより、超解像画像Ｉｃを生成するように構成されている。なお、分割画像５３および分割画像５４は、それぞれ、特許請求の範囲の「第１分割画像」および「第２分割画像」の一例である。
　具体的には、図７に示すように、画像処理部３ｂは、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０に基づいてＸ方向の分割処理を行った分割画像５１および分割画像と同様に、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０に基づいてＸ方向の分割処理を行って分割画像５３および分割画像５４を生成する。すなわち、画像処理部３ｂは、第３Ｘ線透視画像３１（３０）における画素Ｅ３の画素値をＸ方向に分割するＸ方向の分割処理を行って分割画像５３を生成するとともに、第４Ｘ線透視画像４１（４０）における画素Ｅ４の画素値をＸ方向に分割するＸ方向の分割処理を行って分割画像５４を生成する。なお、画素Ｅ３および画素Ｅ４は、それぞれ、特許請求の範囲の「第１画素」および「第２画素」の一例である。
　そして、画像処理部３ｂは、分割画像５３と分割画像５４との互いに対応する画素Ｅの画素値同士を平均した平均画像６２（６０）を生成する。すなわち、画像処理部３ｂは、分割画像５３の画素Ｅ５３の画素値と、（画素Ｅ５３に対応する）分割画像５４の画素Ｅ５４の画素値とを平均して、平均画像６２の画素Ｅ６２とする。これにより、図５に示すように、平均画像６１（６０）と同様に、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０よりもＸ方向に解像度が高くなっている（２倍になっている）平均画像６２（６０）が生成される。すなわち、平均画像６２（６０）は、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０よりもＸ方向における解像度が高い（Ｘ方向に超解像処理が行われた）超解像画像Ｉｂである。
　そして、図８に示すように、画像処理部３ｂは、（Ｘ方向に超解像処理が行われた超解像画像Ｉｂである）平均画像６１および平均画像６２に基づいて、Ｙ方向の分割処理を行って分割画像５５および分割画像５６を生成する。すなわち、画像処理部３ｂは、平均画像６１における画素Ｅ６１の画素値をＹ方向に分割するＹ方向の分割処理を行って分割画像５５を生成するとともに、平均画像６２における画素Ｅ６２の画素値をＹ方向に分割するＹ方向の分割処理を行って分割画像５６を生成する。
　そして、画像処理部３ｂは、分割画像５５と分割画像５６との互いに対応する画素Ｅの画素値同士を平均した平均画像７１（７０）を生成する。すなわち、画像処理部３ｂは、分割画像５５の画素Ｅ５５の画素値と、（画素Ｅ５５に対応する）分割画像５６の画素Ｅ５６の画素値とを平均して、平均画像７１の画素Ｅ７１とする。これにより、図５に示すように、平均画像６１（６０）および平均画像６２（６０）よりもＹ方向に解像度が高くなっている（２倍になっている）平均画像７１（７０）が生成される。具体的には、Ｙ方向において、平均画像７１（７０）の１画素Ｅｃ分の長さＬ５は、平均画像６１（６０）および平均画像６２（６０）の１画素Ｅｂ分の長さＬ４の半分となっている。すなわち、平均画像７１（７０）は、Ｘ線透視画像Ｉａ（第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０）よりもＸ方およびＹ方向における解像度が高い（Ｘ方向およびＹ方向に超解像処理が行われた）超解像画像Ｉｃである。
　なお、図９に示すように、画像処理部３ｂは、（Ｘ方向およびＹ方向に超解像処理が行われた）超解像画像Ｉｃを再構成して、再構成画像７０ｘを生成するように構成されている。すなわち、画像処理部３ｂは、回転ステージ４を用いた断層撮影により取得された複数のＸ線透視画像Ｉａに対して、超解像画像Ｉｃを生成する。具体的には、画像処理部３ｂは、超解像画像Ｉｃである（２次元の）平均画像７１、７２、７３、７４、…を再構成して３次元の再構成画像７０ｘを生成する。なお、平均画像７２は、第１Ｘ線透視画像１２、第２Ｘ線透視画像２２、第３Ｘ線透視画像３２および第４Ｘ線透視画像４２により生成された超解像画像Ｉｃである平均画像７２（７０）である。平均画像７３は、第１Ｘ線透視画像１３、第２Ｘ線透視画像２３、第３Ｘ線透視画像３３および第４Ｘ線透視画像４３により生成された超解像画像Ｉｃである平均画像７３（７０）である。平均画像７４は、第１Ｘ線透視画像１４、第２Ｘ線透視画像２４、第３Ｘ線透視画像３４および第４Ｘ線透視画像４４により生成された超解像画像Ｉｃである平均画像７４（７０）である。
　（超解像画像の生成フロー）
　次に、図１０～１３を参照して、Ｘ線透視画像ＩａよりもＸ方向およびＹ方向における解像度が高い超解像画像Ｉｃの生成フローを説明する。なお、超解像画像Ｉｃの生成に先立って、第１位置Ｐ１、第２位置Ｐ２、第３位置Ｐ３および第４位置Ｐ４において、スキャン撮影が行われることにより、第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０が生成されているものとする。
　まず、図１０に示すように、ステップ１１０において、画像処理部３ｂは、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０に基づいてＸ方向に超解像処理を行った超解像画像Ｉｂを生成する。
　具体的には、図１１に示すように、ステップ１１１において、画像処理部３ｂは、第１Ｘ線透視画像１０の画素Ｅ１を、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像２０において画素Ｅ１に対してオーバラップする２つの画素Ｅ２の画素値に基づいて、Ｘ方向に分割して分割画像５１を生成する。
　次に、ステップ１１２において、画像処理部３ｂは、第２Ｘ線透視画像２０の画素Ｅ２を、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第１Ｘ線透視画像１０において画素Ｅ２に対してオーバラップする２つの画素Ｅ１の画素値に基づいて、Ｘ方向に分割して分割画像５２を生成する。なお、ステップ１１１およびステップ１１２の順序は逆でもよい。
　次に、ステップ１１３において、画像処理部３ｂは、分割画像５１および分割画像５２を、互いに対応する画素Ｅ同士で平均して平均画像６１（Ｘ方向において超解像処理が行われた超解像画像Ｉｂ）を生成する。
　次に、図１０に示すように、ステップ１２０において、画像処理部３ｂは、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０に基づいてＸ方向に超解像処理を行った超解像画像Ｉｂを生成する。なお、ステップ１１０およびステップ１２０の順序は逆でもよい。
　具体的には、図１２に示すように、ステップ１２１において、画像処理部３ｂは、第３Ｘ線透視画像３０の画素Ｅ３を、第３Ｘ線透視画像３０と第４Ｘ線透視画像４０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第４Ｘ線透視画像４０において画素Ｅ３に対してオーバラップする２つの画素Ｅ４の画素値に基づいて、Ｘ方向に分割して分割画像５３を生成する。
　次に、ステップ１２２において、画像処理部３ｂは、第４Ｘ線透視画像４０の画素Ｅ４を、第３Ｘ線透視画像３０と第４Ｘ線透視画像４０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第３Ｘ線透視画像３０において画素Ｅ４に対してオーバラップする２つの画素Ｅ３の画素値に基づいて、Ｘ方向に分割して分割画像５４を生成する。なお、ステップ１２１およびステップ１２２の順序は逆でもよい。
　次に、ステップ１２３において、画像処理部３ｂは、分割画像５３および分割画像５４を、互いに対応する画素Ｅ同士で平均して平均画像６２（Ｘ方向において超解像処理が行われた超解像画像Ｉｂ）を生成する。
　次に、図１０に示すように、ステップ１３０において、画像処理部３ｂは、（Ｘ方向において超解像処理が行われた超解像画像Ｉｂである）平均画像６１および平均画像６２に基づいて、さらにＹ方向に超解像処理を行った超解像画像Ｉｃを生成する。
　具体的には、図１３に示すように、ステップ１３１において、画像処理部３ｂは、平均画像６１と平均画像６２とを（Ｘ線の検出位置の）移動量に対応する分だけＹ方向にずらして重ねて表示した場合に、平均画像６１の画素Ｅ６１を、平均画像６２において画素Ｅ６１に対してオーバラップする２つの画素Ｅ６２の画素値に基づいて、Ｙ方向に分割して分割画像５５を生成する。
　次に、ステップ１３２において、画像処理部３ｂは、平均画像６１と平均画像６２とを（Ｘ線の検出位置の）移動量に対応する分だけＹ方向にずらして重ねて表示した場合に、平均画像６２の画素Ｅ６２を、平均画像６１において画素Ｅ６２に対してオーバラップする２つの画素Ｅ６１の画素値に基づいて、Ｙ方向に分割して分割画像５６を生成する。なお、ステップ１３１およびステップ１３２の順序は逆でもよい。
　次に、ステップ１３３において、画像処理部３ｂは、分割画像５５および分割画像５６を、互いに対応する画素Ｅ同士で平均して平均画像７０（Ｘ方向およびＹ方向において超解像処理が行われた超解像画像Ｉｃ）を生成する。
　（実施形態の効果）
　本実施形態では、以下のような効果を得ることができる。
　本実施形態では、上記のように、画像処理部３ｂを、第１Ｘ線透視画像１０の画素Ｅ１の画素値を、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像２０における画素Ｅ１に対してオーバラップする２つの画素Ｅ２の画素値に基づいてＸ方向に分割するＸ方向の分割処理を行った分割画像５０（分割画像５１）に基づいて、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりもＸ方向における解像度が高い超解像画像Ｉｂを生成するように構成する。これにより、画素Ｅ１の画素値を、２つの画素Ｅ２の画素値に基づいて分割するだけの比較的単純な演算により、実際の画素値の分布を考慮して、Ｘ方向に分割することができる。すなわち、２つの画素Ｅ２の画素値に基づいて画素値を分割するだけの演算なので、従来の超解像処理で行っていた逐次的な演算や多数のパラメータを用いる必要がない。また、実際の画素値Ｅ２の分布を考慮して画素Ｅ１を分割するので、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりも少なくともＸ方向における解像度が高いとともに一定の画質が確保された超解像画像Ｉｂを生成することができる。その結果、超解像画像Ｉｂを生成する際の計算時間が長くなるのを抑制するとともに一定の超解像効果（画質を低下させることなく解像度を高くする効果）を確保しながら、超解像画像Ｉｂを生成することができる。
　また、本実施形態では、上記のように、画像処理部３ｂを、第１Ｘ線透視画像１０の画素Ｅ１の画素値を、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像２０における画素Ｅ１に対してオーバラップする２つの画素Ｅ２の画素値およびオーバラップする面積比率に基づいて分割処理を行った分割画像５０（分割画像５１）に基づいて、超解像画像Ｉｂを生成するように構成する。これにより、画素Ｅ１を、画素Ｅ１に対してオーバラップする２つの画素Ｅ２の画素値に加えて、２つの画素Ｅ２のオーバラップする面積比率が考慮されているので、実際の画素値の分布をより精度よく考慮して、Ｘ方向に分割することができる。その結果、２つの画素Ｅ２のオーバラップする面積比率を考慮しない場合と比較して、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりも少なくともＸ方向における解像度が高いとともにより画質の良い超解像画像Ｉｂを生成することができる。
　また、本実施形態では、上記のように、画像処理部３ｂを、第１Ｘ線透視画像１０における画素Ｅ１の画素値を、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像２０における画素Ｅ１に対してオーバラップする２つの画素Ｅ２の画素値に基づいて分割処理を行った分割画像５１と、第２Ｘ線透視画像２０における画素Ｅ２の画素値を、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第１Ｘ線透視画像１０における画素Ｅ２に対してオーバラップする２つの画素Ｅ１の画素値に基づいて分割処理を行った分割画像５２とに基づいて、超解像画像Ｉｂを生成するように構成する。これにより、分割画像５１または分割画像５２の一方のみを考慮して超解像画像Ｉｂを生成する場合と比較して、分割画像５１および分割画像５２の両方を考慮して超解像画像Ｉｂを生成することにより、Ｘ方向における画素値の分割の精度を向上することができる。その結果、分割画像５１または分割画像５２の一方のみを考慮して超解像画像Ｉｂを生成する場合と比較して、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０よりも少なくともＸ方向における解像度が高いとともにより画質の良い超解像画像Ｉｂを生成することができる。
　また、本実施形態では、上記のように、画像処理部３ｂを、分割画像５１と分割画像５２とを、互いに対応する画素Ｅにおいて平均した平均画像６１を生成することにより、超解像画像Ｉｂを生成するように構成する。これにより、平均画像６１を生成しない場合と比較して、分割画像５１と分割画像５２との間に生じる誤差を低減することができるので、Ｘ方向における画素値の分割の精度を容易に向上することができる。
　また、本実施形態では、上記のように、検出器２は、第１位置Ｐ１および第２位置Ｐ２に加えて、第１位置Ｐ１および第２位置Ｐ２からそれぞれ１画素分よりも小さい移動量だけＸ方向に直交するＹ方向へ平行移動させた第３位置Ｐ３および第４位置Ｐ４とにおいて検出するように構成する。そして、画像処理部３ｂを、第１位置Ｐ１、第２位置Ｐ２、第３位置Ｐ３および第４位置Ｐ４において検出されたＸ線に基づいて、それぞれ、第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０を生成するように構成する。また、画像処理部３ｂを、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０に基づいてＸ方向の分割処理を行った分割画像５０（５１、５２）と、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０に基づいてＸ方向の分割処理を行った分割画像５０（５３、５４）とに基づいて、Ｘ方向における分割処理と同様の処理をＹ方向に対して行うことにより、第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０よりもＸ方向およびＹ方向における解像度が高い超解像画像Ｉｃを生成するように構成する。これにより、Ｘ方向に加えてＹ方向においても解像度が高いとともに一定の画質が確保された超解像画像Ｉｃを生成することができるので、Ｘ方向のみにおいて解像度が高い場合と比較して、Ｘ方向とＹ方向とにおいて、解像度の程度に偏りのない超解像画像Ｉｃを生成することができる。また、一般的な表示装置における画素Ｅは、マトリクス状に配列された検出素子２ａの行方向（Ｘ方向）および列方向（Ｙ方向）におけるサイズが等しいので、Ｘ方向とＹ方向との解像度が異なる場合のように、Ｘ方向とＹ方向とにおける画素サイズを調整する画像の補間処理等の超解像画像Ｉｃを表示装置（図示しない）に表示させるための後工程を簡略化することができる。
　また、本実施形態では、上記のように、画像処理部３ｂを、第１Ｘ線透視画像１０と第２Ｘ線透視画像２０とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第２Ｘ線透視画像２０の端部において、第１Ｘ線透視画像１０の画素Ｅ１に対してオーバラップする画素Ｅ２が１つしかない場合には、オーバラップする１つの画素Ｅ２の画素値に基づいて、もう１つのオーバラップする画素Ｅ２を仮想的に生成するように構成する。これにより、もう１つのオーバラップする画素Ｅを仮想的に生成することにより、第１Ｘ線透視画像１０の端部に位置する画素Ｅ１に対してオーバラップする画素Ｅ２が１つしかない場合でも、画素Ｅ１を分割する演算を確実に行うことができる。また、画像処理部３ｂを、第１Ｘ線透視画像１０の端部において、第２Ｘ線透視画像２０の画素Ｅ２に対してオーバラップする画素Ｅ１が１つしかない場合には、オーバラップする１つの画素Ｅ１の画素値に基づいて、もう１つのオーバラップする画素Ｅ１を仮想的に生成するように構成する。これにより、もう１つのオーバラップする画素Ｅ１を仮想的に生成することにより、第２Ｘ線透視画像２０の端部に位置する画素Ｅ２に対してオーバラップする画素Ｅ１が１つしかない場合でも、画素Ｅ２をＸ方向に分割する演算を確実に行うことができる。
　また、本実施形態では、上記のように、検出器２を、Ｘ線源１から照射されたＸ線を、第１位置Ｐ１（第３位置Ｐ３）と、第１位置Ｐ１（第３位置Ｐ３）から１画素Ｅａ分（の長さＬ４）の半分の長さＬ５分だけＸ方向へ平行移動させた第２位置Ｐ２（第４位置Ｐ４）とにおいて検出するように構成する。これにより、２つの画素Ｅ２（画素Ｅ４）のオーバラップする面積比率が等しくなる（１対１になる）ので、画素Ｅ１（画素Ｅ３）の画素値を、画素Ｅ１（画素Ｅ３）に対してオーバラップする２つの画素Ｅ２（画素Ｅ４）の画素値自体の比率により、Ｘ方向に分割することができる。その結果、２つの画素Ｅ２（画素Ｅ４）のオーバラップする面積比率を２分割以外の比率に分割する場合と比較して、分割処理における演算を簡略化することができるので、超解像画像Ｉｂを生成する際の計算時間が長くなるのをより抑制することができる。
　また、本実施形態では、上記のように、検出器２を、Ｘ線源１から照射されたＸ線をＸ方向を回転軸として回転しながら複数の方向から検出する検出動作と、Ｘ方向への平行移動とを、交互に繰り返し行うことにより、断層撮影を行うように構成する。これにより、ノンヘリカルスキャン型の断層撮影を行うＸ線撮影装置１００において、超解像画像Ｉｂを生成する際の計算時間が長くなるのを抑制するとともに一定の超解像効果を確保しながら超解像画像Ｉｂ、Ｉｃを生成することができる。
　［変形例］
　なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく特許請求の範囲によって示され、さらに特許請求の範囲と均等の意味および範囲内でのすべての変更（変形例）が含まれる。
　たとえば、上記実施形態では、画像処理部３ｂを、仮想的な画素Ｅを、隣接する画素Ｅの画素値と等しい画素値となるように生成するように構成した例を示したが、本発明はこれに限られない。本発明では、画像処理部を、仮想的な画素の画素値を、隣接する画素と異なる画素値となるように生成するように構成してもよい。
　また、上記実施形態では、Ｘ線の検出位置を、１画素Ｅａ分の長さＬ４の半分の長さＬ５分だけ平行移動させるとともに、画像処理部３ｂを、第１Ｘ線透視画像１０（第３Ｘ線透視画像３０）と第２Ｘ線透視画像２０（第４Ｘ線透視画像４０）とを（Ｘ線の検出位置の）移動量に対応する分だけＸ方向にずらして重ねて表示した場合に、第１Ｘ線透視画像１０（第３Ｘ線透視画像３０）の画素Ｅ１（画素Ｅ３）の画素値を、第２Ｘ線透視画像２０（第４Ｘ線透視画像４０）における画素Ｅ１（画素Ｅ３）に対してそれぞれ半分ずつがオーバラップする２つの画素Ｅ２（画素Ｅ４）の画素値に基づいて半分に分割するように構成した例を示したが、本発明はこれに限られない。本発明では、Ｘ線の検出位置を、１画素分の長さの半分よりも長い長さ分だけ平行移動させるように構成してもよいし、１画素分の長さの半分よりも短い長さ分だけ平行移動させるように構成してもよい。この場合、画像処理部を、「（特許請求の範囲の）第１画像」の「（特許請求の範囲の）第１画素」の画素値を、「（特許請求の範囲の）第２画像」における「第１画素」に対してオーバラップする２つの「（特許請求の範囲の）第２画素」の画素値およびオーバラップする面積比率に基づいて分割処理を行って分割画像を生成するように構成すればよい。
　また、上記実施形態では、画像処理部３ｂを、Ｘ方向またはＹ方向のいずれかに平行移動させた位置において検出されたＸ線に基づくＸ線透視画像Ｉａ同士を分割して、分割画像５０を生成するように構成した例を示したが、本発明はこれに限られない。本発明では、画像処理部を、Ｘ方向およびＹ方向のいずれにも平行移動（Ｘ方向およびＹ方向に対して斜めに平行移動）させた位置において検出されたＸ線に基づくＸ線透視画像同士を分割して、分割画像を生成するように構成してもよい。
　また、上記実施形態では、画像処理部３ｂを、２つの分割画像５０を互いに対応する画素Ｅにおいて平均した平均画像６０を生成することにより、超解像画像Ｉｂを生成するように構成した例を示したが、本発明はこれに限られない。本発明では、画像処理部を、１つの分割画像自体を超解像画像として処理するように構成してもよい。
　また、上記実施形態では、Ｘ線撮影装置１００を、超解像画像Ｉｃの生成に先立って、第１位置Ｐ１、第２位置Ｐ２、第３位置Ｐ３および第４位置Ｐ４において、スキャン撮影が行われることにより、第１Ｘ線透視画像１０、第２Ｘ線透視画像２０、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０が生成されているものとして構成した例を示したが、本発明はこれに限られない。本発明では、一の超解像画像の生成に必要なＸ線透視画像が生成された時点で、他の超解像画像の生成に必要なＸ線透視画像が生成される前に、一の超解像画像の生成を開始するように構成してもよい。たとえば、上記実施形態の例に当てはめると、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０が生成された時点で、第３Ｘ線透視画像３０および第４Ｘ線透視画像４０が生成されていない場合でも、第１Ｘ線透視画像１０および第２Ｘ線透視画像２０に基づいてＸ方向に超解像処理を行った超解像画像Ｉｂの生成を開始してもよい。
　また、上記実施形態では、画像処理部３ｂを、２つのＸ線透視画像Ｉａに基づいて、超解像画像Ｉｂを生成するように構成した例を示したが、本発明はこれに限られない。本発明では、画像処理部を、３つ以上のＸ線透視画像に基づいて、超解像画像を生成するように構成してもよい。
　また、上記実施形態では、画像処理部３ｂを、Ｘ方向およびＹ方向において超解像処理が行われた超解像画像Ｉｃを再構成して、再構成画像７０ｘを生成するように構成した例を示したが、本発明はこれに限られない。本発明では、画像処理部を、Ｘ方向またはＹ方向のいずれかのみにおいて超解像処理が行われた超解像画像を再構成して、再構成画像を生成するように構成してもよい。この場合、超解像処理を行わない方向における、Ｘ線の検出位置が異なるＸ線透視画像が不要となる。
　また、上記実施形態では、ノンヘリカルスキャン型の断層撮影を行うＸ線撮影装置１００において、超解像画像Ｉｂ、Ｉｃを生成するように構成した例を示したが、本発明はこれに限られない。本発明では、（螺旋状に動作する）ヘリカルスキャン型の断層撮影を行うＸ線撮影装置において、超解像画像を生成するように構成してもよい。また、断層撮影を行わないＸ線撮影装置において、超解像画像を生成するように構成してもよい。
　１　Ｘ線源
　２　検出器
　３ｂ　画像処理部
　５　検出器移動機構（移動機構）
　１０（１１、１２、１３、１４）　第１Ｘ線透視画像（第１画像）
　２０（２１、２２、２３、２４）　第２Ｘ線透視画像（第２画像）
　３０（３１、３２、３３、３４）　第３Ｘ線透視画像（第３画像）
　４０（４１、４２、４３、４４）　第４Ｘ線透視画像（第４画像）
　５０　分割画像
　５１、５３　分割画像（第１分割画像）
　５２、５４　分割画像（第２分割画像）
　６０（６１、６２）　平均画像
　９０　回転軸
　１００　Ｘ線撮影装置
　Ｅ（Ｅａ、Ｅｂ、Ｅｃ）　画素
　Ｅ１、Ｅ３　画素（第１画素）
　Ｅ２、Ｅ４　画素（第２画素）
　Ｉｂ、Ｉｃ　超解像画像
　Ｐ１　第１位置
　Ｐ２　第２位置
　Ｐ３　第３位置
　Ｐ４　第４位置

Claims (8)

  1. 　Ｘ線源と、
    　前記Ｘ線源から照射されたＸ線を、第１位置と、前記第１位置から１画素分よりも小さい移動量だけ第１方向へ平行移動させた第２位置とにおいて検出する検出器と、
    　前記第１位置において検出されたＸ線に基づいて第１画像を生成するとともに、前記第２位置において検出されたＸ線に基づいて第２画像を生成する画像処理部と、
    　前記検出器によるＸ線の検出位置を、前記第１位置と前記第２位置との間で移動させる移動機構と、
    を備え、
    　前記画像処理部は、前記第１画像および前記第２画像の一方の第１画素の画素値を、前記第１画像と前記第２画像とを前記移動量に対応する分だけ前記第１方向にずらして重ねて表示した場合に、前記第１画像および前記第２画像の他方における前記第１画素に対してオーバラップする２つの画素の画素値に基づいて、前記第１方向に分割する前記第１方向の分割処理を行った分割画像に基づいて、前記第１画像および前記第２画像よりも前記第１方向における解像度が高い超解像画像を生成するように構成されている、Ｘ線撮影装置。
  2. 　前記画像処理部は、前記第１画像および前記第２画像の一方の前記第１画素の画素値を、前記第１画像と前記第２画像とを前記移動量に対応する分だけ前記第１方向にずらして重ねて表示した場合に、前記第１画像および前記第２画像の他方における前記第１画素に対してオーバラップする２つの画素の画素値およびオーバラップする面積比率に基づいて前記分割処理を行った前記分割画像に基づいて、前記超解像画像を生成するように構成されている、請求項１に記載のＸ線撮影装置。
  3. 　前記画像処理部は、前記第１画像における前記第１画素の画素値を、前記第１画像と前記第２画像とを前記移動量に対応する分だけ前記第１方向にずらして重ねて表示した場合に、前記第２画像における前記第１画素に対してオーバラップする２つの画素の画素値に基づいて前記分割処理を行った第１分割画像と、前記第２画像における第２画素の画素値を、前記第１画像と前記第２画像とを前記移動量に対応する分だけ前記第１方向にずらして重ねて表示した場合に、前記第１画像における前記第２画素に対してオーバラップする２つの画素の画素値に基づいて前記分割処理を行った第２分割画像とに基づいて、前記超解像画像を生成するように構成されている、請求項１または２に記載のＸ線撮影装置。
  4. 　前記画像処理部は、前記第１分割画像と前記第２分割画像とを、互いに対応する画素において平均した平均画像を生成することにより、前記超解像画像を生成するように構成されている、請求項３に記載のＸ線撮影装置。
  5. 　前記検出器は、前記第１位置および前記第２位置に加えて、前記第１位置および前記第２位置からそれぞれ１画素分よりも小さい前記移動量だけ前記第１方向に直交する第２方向へ平行移動させた第３位置および第４位置とにおいて検出するように構成されており、
    　前記画像処理部は、前記第１位置、前記第２位置、前記第３位置および前記第４位置において検出されたＸ線に基づいて、それぞれ、前記第１画像、前記第２画像、第３画像および第４画像を生成するとともに、前記第１画像および前記第２画像に基づいて前記第１方向の前記分割処理を行った前記分割画像と、前記第３画像および前記第４画像に基づいて前記第１方向の前記分割処理を行った前記分割画像とに基づいて、前記第１方向における前記分割処理と同様の処理を前記第２方向に対して行うことにより、前記第１画像、前記第２画像、前記第３画像および前記第４画像よりも前記第１方向および前記第２方向における解像度が高い前記超解像画像を生成するように構成されている、請求項１～４のいずれか１項に記載のＸ線撮影装置。
  6. 　前記画像処理部は、前記第１画像または前記第２画像の他方の端部において、前記第１画像と前記第２画像とを前記移動量に対応する分だけ前記第１方向にずらして重ねて表示した場合に、前記第１画像および前記第２画像の一方の前記第１画素に対してオーバラップする画素が１つしかない場合には、前記オーバラップする１つの画素の画素値に基づいて、もう１つの前記オーバラップする画素を仮想的に生成するように構成されている、請求項１～５のいずれか１項に記載のＸ線撮影装置。
  7. 　前記検出器は、前記Ｘ線源から照射されたＸ線を、前記第１位置と、前記第１位置から１画素分の半分の長さ分だけ前記第１方向へ平行移動させた前記第２位置とにおいて検出するように構成されている、請求項２に記載のＸ線撮影装置。
  8. 　前記検出器は、前記Ｘ線源から照射されたＸ線を、第３方向を回転軸として回転しながら複数の方向から検出する検出動作と、前記第３方向への平行移動とを、交互に繰り返し行うことにより、断層撮影を行うように構成されている、請求項１～７のいずれか１項に記載のＸ線撮影装置。
PCT/JP2018/046996 2018-12-20 2018-12-20 Ｘ線撮影装置 WO2020129213A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2020560725A JP7205552B2 (ja) 2018-12-20 2018-12-20 Ｘ線撮影装置
US17/414,069 US11821851B2 (en) 2018-12-20 2018-12-20 X-ray imaging apparatus
CN201880099878.1A CN113226180B (zh) 2018-12-20 2018-12-20 X射线摄影装置
PCT/JP2018/046996 WO2020129213A1 (ja) 2018-12-20 2018-12-20 Ｘ線撮影装置
TW108124892A TWI708937B (zh) 2018-12-20 2019-07-15 X射線攝影裝置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2018/046996 WO2020129213A1 (ja) 2018-12-20 2018-12-20 Ｘ線撮影装置

Publications (1)

Publication Number Publication Date
WO2020129213A1 true WO2020129213A1 (ja) 2020-06-25

Family

ID=71102656

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/046996 WO2020129213A1 (ja) 2018-12-20 2018-12-20 Ｘ線撮影装置

Country Status (5)

Country Link
US (1) US11821851B2 (ja)
JP (1) JP7205552B2 (ja)
CN (1) CN113226180B (ja)
TW (1) TWI708937B (ja)
WO (1) WO2020129213A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2023006563A (ja) * 2021-06-30 2023-01-18 キヤノンメディカルシステムズ株式会社 Ｘ線診断装置、ｘ線診断方法、およびプログラム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10155158A (ja) * 1996-11-20 1998-06-09 Sony Corp 撮像装置及びカラー画像信号の処理方法
WO2017007024A1 (ja) * 2015-07-08 2017-01-12 国立大学法人香川大学 分光測定装置
JP2017045431A (ja) * 2015-08-28 2017-03-02 パナソニックＩｐマネジメント株式会社 画像出力装置、画像送信装置、画像受信装置、画像出力方法およびプログラム
JP2018108364A (ja) * 2016-12-28 2018-07-12 キヤノンメディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7492967B2 (en) * 2003-09-24 2009-02-17 Kabushiki Kaisha Toshiba Super-resolution processor and medical diagnostic imaging apparatus
JP4519434B2 (ja) 2003-09-24 2010-08-04 株式会社東芝 超解像処理装置及び医用画像診断装置
JP4720299B2 (ja) * 2005-06-07 2011-07-13 株式会社島津製作所 断層撮影装置
JP4806433B2 (ja) * 2008-08-11 2011-11-02 株式会社吉田製作所 Ｘ線撮影装置およびその信号処理方法
US8559590B2 (en) * 2010-01-28 2013-10-15 Varian Medical Systems, Inc. Imaging breast cancerous lesions with microcalcifications
JP2011217264A (ja) * 2010-04-01 2011-10-27 Toshiba Corp Ｘ線平面検出器、画像データ転送システム及び画像データ転送方法
FI20110106L (fi) * 2011-03-21 2012-04-13 Planmeca Oy Hammaslääketieteellinen kuvauslaitteisto
JP2013158444A (ja) * 2012-02-03 2013-08-19 Toshiba Corp Ｘ線診断装置
US9459217B2 (en) * 2013-04-12 2016-10-04 Illinois Tool Works, Inc. High-resolution computed tomography
JP2016099308A (ja) * 2014-11-26 2016-05-30 株式会社日立ハイテクサイエンス 蛍光ｘ線分析装置及び蛍光ｘ線分析方法
US10531855B2 (en) * 2016-12-28 2020-01-14 Canon Medical Systems Corporation X-ray computed tomography apparatus

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10155158A (ja) * 1996-11-20 1998-06-09 Sony Corp 撮像装置及びカラー画像信号の処理方法
WO2017007024A1 (ja) * 2015-07-08 2017-01-12 国立大学法人香川大学 分光測定装置
JP2017045431A (ja) * 2015-08-28 2017-03-02 パナソニックＩｐマネジメント株式会社 画像出力装置、画像送信装置、画像受信装置、画像出力方法およびプログラム
JP2018108364A (ja) * 2016-12-28 2018-07-12 キヤノンメディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YONEYAMA, AKIO ET AL.: "Non-official translation: Examination of high-resolution CT by S-17 X-ray sub-pixel shifting (SPS) method", SUNBEAM ANNUAL REPORT WITH RESEARCH RESULTS, vol. 4, no. 3, 2014, pages 122 - 123 *

Also Published As

Publication number Publication date
JP7205552B2 (ja) 2023-01-17
TWI708937B (zh) 2020-11-01
US11821851B2 (en) 2023-11-21
TW202024620A (zh) 2020-07-01
JPWO2020129213A1 (ja) 2021-10-14
CN113226180A (zh) 2021-08-06
US20220057340A1 (en) 2022-02-24
CN113226180B (zh) 2023-08-04

Similar Documents

Publication Publication Date Title
JP5537074B2 (ja) コンピュータ断層撮影装置及び方法
US7396162B1 (en) Scatter correction for CT method and apparatus
US9532755B2 (en) X-ray computed tomographic imaging apparatus and method for same
CN105828719A (zh) X射线成像装置
JP2008012319A (ja) トモシンセシス・イメージング・システムでのアーティファクトを低減する方法及びシステム
JP2010259653A (ja) Ｘ線コンピュータ断層撮影装置
US20180055459A1 (en) Image processing apparatus, method, and program
US20040105528A1 (en) Method and system for tomosynthesis image enhancement using transverse filtering
KR101284986B1 (ko) 고해상도 토모신세시스 단면 영상의 재구성 방법 및 그 장치
US20090207964A1 (en) Method and system for reconstructing cone-beam projection data with reduced artifacts
WO2020129213A1 (ja) Ｘ線撮影装置
JP7120442B2 (ja) Ｘ線撮影装置
JP2010104480A (ja) Ｘ線ｃｔ装置およびプログラム
JP5872212B2 (ja) マルチスライスｃｔ装置及びデータ前処理方法
US20180306735A1 (en) X-ray phase imaging apparatus
JP2018143574A (ja) Ｘ線ｃｔ装置及び画像処理方法
JP2016059610A (ja) 医用画像処理装置および医用画像処理装置における画像データ表示方法
US10258296B2 (en) X-ray CT apparatus including processing circuitry to improve a spatial resolution in a row direction and a channel direction
US20200268326A1 (en) X-ray examination arrangement and method for operating an x-ray examination arrangement
EP2433265B1 (en) Interpolation free fan-to-parallel beam rebinning
JP7487683B2 (ja) 放射線画像生成方法および放射線画像撮影装置
JP7233865B2 (ja) Ｘ線ｃｔ装置及び補正方法
CN111093508B (zh) 医用x射线图像处理装置和x射线图像摄影装置
WO2014038588A1 (ja) Ｘ線ｃｔ装置及び医用画像処理方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18943971

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2020560725

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

Country of ref document: EP

Kind code of ref document: A1