WO2013005833A1 - Ｘ線撮影装置およびそのキャリブレーション方法 - Google Patents

Ｘ線撮影装置およびそのキャリブレーション方法 Download PDF

Info

Publication number
WO2013005833A1
WO2013005833A1 PCT/JP2012/067336 JP2012067336W WO2013005833A1 WO 2013005833 A1 WO2013005833 A1 WO 2013005833A1 JP 2012067336 W JP2012067336 W JP 2012067336W WO 2013005833 A1 WO2013005833 A1 WO 2013005833A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
ray
marker
ray source
calibration
Prior art date
Application number
PCT/JP2012/067336
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
Priority claimed from JP2011149912A external-priority patent/JP5702236B2/ja
Priority claimed from JP2011153815A external-priority patent/JP5702240B2/ja
Application filed by 富士フイルム株式会社 filed Critical 富士フイルム株式会社
Priority to EP12807484.6A priority Critical patent/EP2730226A4/en
Priority to CN201280033444.4A priority patent/CN103635139A/zh
Publication of WO2013005833A1 publication Critical patent/WO2013005833A1/ja
Priority to US14/147,113 priority patent/US9380985B2/en

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/025—Tomosynthesis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/04—Positioning of patients; Tiltable beds or the like
    • A61B6/0492—Positioning of patients; Tiltable beds or the like using markers or indicia for aiding patient positioning
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/42—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4233—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/582—Calibration
    • A61B6/585—Calibration of detector units
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/587—Alignment of source unit to detector unit
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/588—Setting distance between source unit and detector unit
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/42—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4283—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by a detector unit being housed in a cassette
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4452—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being able to move relative to each other
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4464—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit or the detector unit being mounted to ceiling
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/502—Clinical applications involving diagnosis of breast, i.e. mammography
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10016—Video; Image sequence
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30204—Marker
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30241—Trajectory

Definitions

  • the present invention relates to an X-ray imaging apparatus that performs imaging in a supine position, a standing position, a mammo, and the like, and its calibration method, and in particular, even if the moving axis of an X-ray source is shifted with respect to the X-ray detector.
  • the X-ray imaging apparatus capable of correcting the deviation, obtaining a high-accuracy composite image (tomosynthesis image) for a cross-section of the subject at a predetermined position, and further suppressing the deterioration of the deviation over time, and the same
  • the present invention relates to a calibration method.
  • the X-ray source is moved to irradiate the subject with X-rays from different angles, and the obtained images are added to obtain an image that emphasizes the desired tomographic plane.
  • Various X-ray imaging apparatuses that can perform tomosynthesis imaging have been proposed (see, for example, Patent Document 1).
  • tomosynthesis imaging for example, the X-ray source is moved in parallel with the X-ray detector, or is moved so as to draw an arc such as a circle or an ellipse, and a plurality of captured images obtained by imaging the subject at different irradiation angles are acquired. These tomographic images are reconstructed to create a tomographic image.
  • Japanese Patent Application Laid-Open No. 2004-228561 selects either a first mode in which only the radiation source is moved or a second mode in which both the radiation source and the detection unit are moved according to the imaging conditions during tomosynthesis imaging.
  • a radiation imaging apparatus capable of acquiring a plurality of captured images according to a selected mode is disclosed.
  • tomosynthesis imaging high accuracy is required for alignment of the X-ray detector (X-ray image sensor) and the movement axis of the X-ray source.
  • the X-ray source and the X-ray detector (X-ray image sensor) are not integrated, and are installed so that alignment can be achieved with high accuracy. It's not easy.
  • tomosynthesis imaging a method of generating calibration data using a predetermined marker assembly is known.
  • an X-ray image is corrected by detecting a positional deviation of an image using a marker (see, for example, Patent Document 2).
  • a supporter having a marker is attached to an imaging region of a subject, a mask image and a live image of the imaging region are captured, and an error between the marker position of the mask image and the marker position of the live image.
  • An X-ray imaging apparatus is disclosed in which subtraction is performed so that these values coincide with each other, and noise in the subtraction image caused by positional deviation of the catatization table is prevented.
  • Patent Document 3 discloses an X-ray CT apparatus that is suitable for operation in a clinical field and aims to acquire a high-accuracy CT image free from ring artifacts.
  • the air calibration data acquisition unit acquires new air calibration data
  • the air calibration data exchange unit Replace the air calibration data that has been used up to that point with new air calibration data.
  • Patent Document 3 discloses that the air calibration data is replaced, it relates to an X-ray CT apparatus, and there is currently no thing that can cope with deterioration of calibration with time in tomosynthesis imaging. is there.
  • the object of the present invention is to solve the problems based on the prior art, and even if the movement axis of the X-ray source is shifted with respect to the X-ray detector, the shift can be corrected, and the predetermined position of the subject can be corrected.
  • An X-ray imaging apparatus capable of obtaining a high-accuracy composite image (tomosynthesis image) for the cross section of the image, further correcting the deviation, and suppressing the deterioration of the deviation with time, and a calibration method thereof There is to do.
  • the first aspect of the present invention is directed to an X-ray source capable of moving along a predetermined movement path, a moving means for moving the X-ray source along a predetermined movement path, and facing the X-ray source.
  • the X-ray source is moved by a moving table, an X-ray detector provided on the imaging table, a flat X-ray detector provided on the imaging table, a marker provided on the imaging table, and a moving means.
  • a control unit that captures an image including a marker, and the position of the marker image for each captured image, and the inclination of the moving axis of the X-ray source with respect to the X-ray detector based on the relative relationship of the position of the marker image
  • An X-ray imaging apparatus characterized by having an image processing unit for obtaining the above.
  • the image processing unit preferably obtains an inclination with respect to the X-ray detector for each imaging position of the X-ray source.
  • the image processing unit preferably uses the inclination with respect to the X-ray detector for reconstruction of a tomographic image.
  • the inclination with respect to the X-ray detector is, for example, a first inclination in a plane parallel to the surface of the X-ray detector.
  • the image processing unit further calculates a variation in the size of the marker image for each captured image, obtains the distance between the X-ray source and the X-ray detector using the variation, and calculates the inclination with respect to the X-ray detector. It is preferable to obtain the second inclination in a plane orthogonal to the surface of the X-ray detector.
  • a calibration method for an X-ray imaging apparatus wherein the X-ray imaging apparatus includes an imaging table provided facing an X-ray source movable along a predetermined movement path, and imaging.
  • a flat plate-shaped X-ray detector provided on the stage and a marker; and a step of moving the X-ray source to capture images including the marker from at least two positions; and
  • a calibration method comprising a step of obtaining a position of a marker image with respect to an image, and obtaining an inclination of the moving axis of the X-ray source with respect to the X-ray detector based on a relative relationship between the positions of the marker images. To do. It is preferable to obtain an inclination with respect to the X-ray detector for each imaging position of the X-ray source.
  • the inclination with respect to the X-ray detector is, for example, a first inclination in a plane parallel to the surface of the X-ray detector.
  • the variation in the size of the marker image is calculated for each photographed image, the distance between the X-ray source and the X-ray detector is obtained using the variation, and the X-ray detector is used as the inclination with respect to the X-ray detector. It is preferable to have the process of calculating
  • an X-ray source movable along a predetermined movement path, a moving means for moving the X-ray source along a predetermined movement path, and an imaging table provided opposite to the X-ray source.
  • a calibration unit that performs calibration relating to a shift in the moving direction of an X-ray source disposed on the X-ray detector, a flat X-ray detector provided on the imaging table, a marker provided on the imaging table And a determination unit that determines whether or not the calibration unit needs to be updated, and causes the calibration unit to perform calibration based on the determination result, and an image including the marker is captured by the X-ray source, The position of the marker image of the image captured by the calibration unit is obtained, and the determination unit determines whether or not the calibration needs to be updated based on the position of the marker image.
  • an X-ray imaging apparatus There is provided an X-ray imaging apparatus.
  • the calibration unit moves the X-ray source so that images including markers are captured from at least two positions, and the positions of the marker images are obtained for the first and second images that have been captured, and the determination is made.
  • the unit obtains the position of the corresponding marker image corresponding to the marker image of the second image from the position of the marker image of the first image using the existing calibration data, and the position of the marker image of the second image; It is preferable to determine whether the calibration needs to be updated based on a deviation from the position of the corresponding marker image corresponding to the marker image of the second image.
  • the calibration unit moves the X-ray source to capture an image including a marker from at least two positions, obtains the position of the marker image for each captured image, and determines the relative relationship between the positions of the marker images. Based on this, a first inclination of the movement axis of the X-ray source in a plane parallel to the surface of the X-ray detector is obtained, and first calibration data is created based on the first inclination. It is preferable.
  • the calibration unit moves the X-ray source, captures an image including a marker from at least two positions, calculates a variation in the size of the marker image for each captured image, and uses the variation.
  • a distance between the X-ray source and the X-ray detector is obtained, a second inclination of the X-ray source in a plane orthogonal to the surface of the X-ray detector is obtained, and a second calibration is performed based on the second inclination. It is preferable to create the action data.
  • a calibration part calculates
  • the X-ray source is moved to a predetermined imaging position by the moving means and has a reconstruction unit for reconstructing a tomographic image using a plurality of images obtained by tomosynthesis imaging.
  • reconstructing a tomographic image using a plurality of captured images it is preferable to use at least the first calibration data among the first calibration data and the second calibration data.
  • a method for calibrating an X-ray imaging apparatus wherein the X-ray imaging apparatus includes an imaging table provided opposite to an X-ray source movable along a predetermined movement path, and imaging.
  • a flat plate X-ray detector provided on the table and a marker, a step of photographing an image including the marker with an X-ray source, a position of the marker image of the photographed image is obtained, and this marker image
  • the calibration step is a step of moving the X-ray source to capture images including markers from at least two positions, and obtaining the positions of the marker images for the captured first image and second image, respectively. Determining the position of the corresponding marker image corresponding to the marker image of the second image from the position of the marker image of the first image using the existing calibration data; and the position of the marker image of the second image; It is preferable to include a step of determining whether or not the calibration needs to be updated based on a deviation from the position of the corresponding marker image corresponding to the marker image of the second image.
  • the X-ray source is moved, and images including markers are captured from at least two positions, the positions of the marker images are obtained for the captured images, and the relative relationship between the positions of the marker images is determined.
  • the X-ray source is moved, and images including markers are captured from at least two positions, the variation in the size of the marker image is calculated for each captured image, and the variation is used.
  • a distance between the X-ray source and the X-ray detector is obtained, a second inclination of the X-ray source in a plane orthogonal to the surface of the X-ray detector is obtained, and a second calibration is performed based on the second inclination. It is preferable to create the action data. In addition, it is preferable to obtain at least first calibration data among the first calibration data and the second calibration data for each imaging position of the X-ray source.
  • the shift can be corrected. For this reason, when reconstructing tomosynthesis, reconstruction using each image can be performed with high accuracy. As a result, a highly accurate composite image (tomosynthesis image) can be obtained for the cross section of the subject at a predetermined position.
  • the present invention it is possible to easily determine whether or not the calibration data for calibrating the deviation between the X-ray detector and the moving axis of the X-ray source needs to be updated due to a change with time. Can do. Thereby, it is possible to easily know whether calibration is necessary, and calibration can be performed as necessary. Therefore, it is possible to always reconstruct tomosynthesis with high accuracy while suppressing the influence of changes over time, and to always obtain a high-accuracy synthesized image (tomosynthesis image) for a cross section of a predetermined position of the subject. it can.
  • 1 is a schematic diagram showing an X-ray imaging apparatus according to a first embodiment of the present invention. It is a schematic diagram which shows the arrangement
  • (A) And (b) is a schematic diagram which shows the fluctuation
  • (A) And (b) is a conceptual diagram which shows the mode at the time of the reconstruction of the X-ray tomographic image by tomosynthesis imaging
  • FIG. 1 is a schematic diagram showing an X-ray imaging apparatus according to the first embodiment of the present invention.
  • the X-ray imaging apparatus 10 of the first embodiment shown in FIG. 1 performs tomosynthesis imaging of a subject M such as a human body and obtains an X-ray tomographic image at a cross section of the subject M at an arbitrary height.
  • the X-ray imaging apparatus 10 can be used for, for example, recumbent position, standing position, and mammo tomosynthesis imaging.
  • the X-ray detector 18 can use built-in.
  • the X-ray imaging apparatus 10 can use a cassette type X-ray detector 18 as long as the cassette position does not change before and after calibration.
  • the X-ray imaging apparatus 10 includes an overhead traveling X-ray source 12, a moving unit 14, an X-ray source control unit 15, an imaging table 16, an X-ray detector 18, a marker 20, an image processing unit 22, a display unit 24, and an output. Part 26 and control part 30.
  • the control unit 30 is connected to the moving unit 14, the X-ray source control unit 15, the X-ray detector 18, the image processing unit 22, the display unit 24, the output unit 26, and the input unit 28.
  • the control unit 30 captures X-ray images by the X-ray source 12, the moving unit 14, and the X-ray source control unit 15, image processing in the image processing unit 22, and display unit 24.
  • the operation of the X-ray imaging apparatus 10 such as screen display and output processing in the output unit 26 is controlled.
  • the input unit 28 is a part for inputting various instructions including an instruction to start photographing and a switching instruction to be described later, and examples thereof include a mouse and a keyboard.
  • An instruction input is input to the control unit 30 via the input unit 28, and various operations of the X-ray imaging apparatus 10 are performed by the control unit 30.
  • the instruction input input via the input unit 28 is displayed on the display unit 24, for example.
  • the X-ray source 12 is movably attached to a rail 13 provided on the ceiling H so as to face the imaging table 16. As shown in FIG. 2, the rail 13 is arranged so that the moving direction (y direction) of the X-ray source 12 and the longitudinal direction (Y direction) of the X-ray detector 18 coincide, that is, the X-ray source 12.
  • the moving direction (y direction) of the X-ray detector 18 and the longitudinal direction (Y direction) of the X-ray detector 18 are arranged in parallel.
  • a moving portion 14 is provided on the rail 13.
  • the moving unit 14 includes, for example, a driving unit (not shown) such as a stepping motor connected to the X-ray source 12 via a belt (not shown) and a pulley (not shown).
  • the moving unit 14 allows the X-ray source 12 to move linearly along the rail 13 and can also stop at a predetermined position on the rail 13.
  • the configuration of the moving unit 14 is not particularly limited as long as the X-ray source 12 can be linearly moved along the rail 13 and the X-ray source 12 can be stopped at a predetermined imaging position. Absent.
  • the movement path of the X-ray source 12 is not limited to a straight line, but may be an arc shape. In this case, the rail 13 is provided in an arc shape.
  • the X-ray source 12 irradiates the subject M on the surface 16 a of the imaging table 16 with X-rays (radiation), and can change the X-ray irradiation angle in the direction of the subject M.
  • X-ray source 12 an X-ray source that can be used for a general X-ray imaging apparatus and that can change the X-ray irradiation angle can be used.
  • the X-ray source 12 is provided with a collimator (irradiation field stop) (not shown).
  • the X-ray source 12 is connected to the X-ray source control unit 15.
  • the X-ray source control unit 15 controls X-ray irradiation and X-ray irradiation angle from the X-ray source 12.
  • the X-ray source control unit 15 changes the X-ray irradiation angle in the direction of the subject M, and then X X-rays having a predetermined intensity are irradiated from the radiation source 12 for a predetermined time. That is, the X-ray source 12 irradiates the subject M with X-rays having a predetermined dose (dose).
  • the imaging stand 16 is a surface on which a subject M such as a person or an animal is placed on the surface 16a.
  • a subject M for example, a person is placed on the surface 16a of the photographing stand 16 in a supine position.
  • a flat type (flat plate) X-ray detector 18 (hereinafter referred to as FPD 18) is provided on the back surface 16 b of the imaging table 16 so as to face the rail 13.
  • FPD 18 flat type (flat plate) X-ray detector 18
  • the longitudinal direction of the FPD 18 Y direction, see FIG. 2
  • the longitudinal direction of the rail 13 that is, the moving direction of the X-ray source 12 (y direction) are parallel. Installed.
  • the FPD 18 is disposed on the back surface 16b of the imaging table 16 with the X-ray light receiving surface facing upward.
  • the FPD 18 detects X-rays transmitted through the subject M, performs photoelectric conversion, and obtains digital image data (projection data) corresponding to the photographed X-ray image of the subject M.
  • the movement direction (y direction) of the X-ray source 12 and the longitudinal direction (Y direction) of the FPD 18 are arranged in parallel, and the movement direction (y direction) of the X-ray source 12 is , And are arranged so as to coincide with the pixel arrangement direction of the FPD 18.
  • inclination refers to an inclination with respect to an arrangement of a plurality of pixels (a two-dimensional matrix arrangement direction).
  • the FPD 18 there are various methods such as a direct method for directly converting X-rays (radiation) into electric charges, an indirect method for once converting X-rays (radiation) into light, and further converting the converted light into electric signals. Is available. Further, the FPD 18 may be configured to be movable with respect to the moving direction of the X-ray source 12.
  • the imaging stand 16 is provided with, for example, one marker 20 within a range where the FPD 18 is provided.
  • the marker 20 is provided on the imaging stand 16 with a gap from the surface 18a of the FPD 18. Further, the marker 20 may be provided with the edge aligned with the outer edge 18 b of the FPD 18.
  • the marker 20 is used for calibrating a plurality of pieces of projection data (X-ray images) that have been photographed before the tomosynthesis photographing of the subject M, for example, when the apparatus is activated.
  • the marker 20 is, for example, a metal, preferably a metal having a large atomic number, and is made of a material that hardly transmits radiation. It is desirable that the marker 20 has a shape in which the center of the marker 20 is easy to understand, such as a shape with a cross hole.
  • the number of markers 20 is not particularly limited, but there may be a plurality of markers. In the case of a plurality, it is preferable that the markers are different from each other in order to facilitate identification of each marker. Furthermore, when there are a plurality of markers 20, for example, four are provided. In addition, when there are a plurality of markers 20, the arrangement pattern is not limited at all. For example, they are arranged in a rectangular shape.
  • the marker 20 may be detachable from the imaging table 16.
  • calibration can be performed only when calibration is performed, for example, when the marker 20 is appropriately placed on the imaging table 16 and imaged at the time of starting the apparatus.
  • calibration can be easily performed by making the marker 20 detachable.
  • the X-ray irradiation angle in the direction of the subject M at a predetermined imaging position while moving the X-ray source 12 in the y direction by the moving unit 14. And the subject M is irradiated with X-rays at different imaging angles (fixed time intervals). X-rays emitted from the X-ray source 12 pass through the subject M, enter the light receiving surface of the FPD 18, are detected and photoelectrically converted by the FPD 18, and projection data corresponding to the photographed X-ray image of the subject M ( Digital image data) is acquired.
  • a plurality of X-ray images (for example, 20 to 80 images) with different imaging angles of the subject M are captured by one imaging operation, and the plurality of X-ray images captured from the FPD 18 are captured.
  • Projection data (digital image data) corresponding to is sequentially output to the image processing unit 22 (storage unit 22a).
  • the image processing unit 22 receives projection data of a plurality of X-ray images acquired by the FPD 18 and performs image processing (correction processing (calibration), image synthesis) using the projection data of the plurality of X-ray images. This is a part for reconstructing an X-ray tomographic image in a cross section at an arbitrary height of the subject M.
  • the image processing unit 22 includes a storage unit 22a, a correction unit 22b (calibration unit), and a reconstruction unit 22c.
  • the storage unit 22a receives a plurality of projection data of X-ray images acquired by the FPD 18 and stores the projection data of these X-ray images.
  • the design size of the marker 20 in the projection data, the threshold value of the shift amount (shift amount ⁇ ), and the like are also stored in the storage unit 22a.
  • the correction unit 22b calculates a positional deviation amount (position information of the X-ray source 12) at each imaging position of the X-ray source 12, and when the positional deviation amount is larger than a predetermined value (for example, a threshold value), The position information (X-ray tube position information) of the X-ray source 12 of the projection data imaged at the imaging position with a large amount of deviation is corrected. That is, the deviation amount is calculated for all the projection data of the X-ray image obtained by tomosynthesis imaging, and when the deviation amount of the position is larger than a predetermined value (for example, a threshold value), the imaging position with the larger deviation amount is used. The position information (X-ray tube position information) of the X-ray source 12 of the captured projection data is corrected.
  • a predetermined value for example, a threshold value
  • the position information (X-ray tube position information) of the X-ray source 12 based on the shift amount may be corrected for all the projection data without setting a threshold for the shift amount.
  • the position information of the X-ray source 12 (X-ray tube position information) is based on the design value (the target value at which each of the plurality of images is taken) and the movement axis of the X-ray source 12 is tilted. By correcting the amount, the information is corrected to be close to the position of the true X-ray source 12 and used for the reconstruction calculation.
  • the correction unit 22b has an image recognition function for the marker 20 (marker image).
  • the inclination of the movement axis of the X-ray source 12 with respect to the FPD 18 (X-ray detector) is a positional shift amount of the X-ray source 12 at the imaging position.
  • the positional deviation amount at the imaging position of the X-ray source 12 includes the deviation amount of the distance between the X-ray source 12 and the surface 18a of the FPD 18 (distance L between the X-ray source and the FPD). Is also included.
  • the correction unit 22b when the amount of deviation is larger than a predetermined value (for example, a threshold value), based on the amount of deviation of the position, for example, in the case of deviation with respect to the longitudinal direction (Y direction) of the FPD 18, depending on the amount of deviation.
  • a predetermined value for example, a threshold value
  • the distance L between the X-ray source and the FPD is set with respect to the deviation amount of the distance between the X-ray source 12 and the surface 18a of the FPD 18 (the distance L between the X-ray source and the FPD).
  • the change rate of the marker 20 (marker image) that is, the enlargement rate or reduction rate is obtained with respect to the deviation amount of the distance L between the X-ray source and the FPD, and this is stored in, for example, the storage unit 22a.
  • the shift in the z direction shown in FIG. 1 can be corrected by enlarging or reducing the projection data in accordance with the enlargement rate or reduction rate. It goes without saying that the correction unit 22b may correct both the positional deviation amount and the distance deviation amount.
  • the reconstruction unit 22c corrects the projection data of a plurality of X-ray images including the projection data after correction processing by the correction unit 22b performed as necessary, and corrects the positional information of the X-ray source 12 as necessary.
  • An image composition process is performed using (X-ray tube position information) to reconstruct an X-ray tomographic image of a cross section of the subject M at an arbitrary height.
  • the image processing unit 22 (the storage unit 22a, the correction unit 22b (calibration unit), and the reconstruction unit 22c) is not limited to being configured by hardware (device).
  • a program for causing a computer to execute processing by the correction unit 22b (calibration unit) and the reconstruction unit 22c of the image processing unit 22 on the projection data stored in the storage unit 22a may be used.
  • the display unit 24 is a part that displays the content of the instruction input signal input from the input unit 28 to the control unit 30 and the X-ray tomographic image reconstructed by the image processing unit 22 in accordance with the instruction input signal.
  • a flat panel display such as a liquid crystal display can be exemplified.
  • the output unit 26 is a part that outputs the X-ray tomographic image reconstructed by the image processing unit 22 in accordance with an instruction input signal input from the input unit 28 to the control unit 30.
  • the X-ray tomographic image is printed out. Examples include various printers, storage devices that store digital image data of X-ray tomographic images in various recording media, and the like.
  • tomosynthesis imaging using the X-ray imaging apparatus 10 of the present embodiment high accuracy is required for alignment of the movement axes of the FPD 18 (X-ray detector 18) and the X-ray source 12 as described above.
  • the X-ray source 12 and the FPD 18 are installed with high accuracy because the X-ray source 12 and the FPD 18 are not integrated. Difficult to do. Therefore, when moving the X-ray source 12, as shown in FIG. 2, with respect to the longitudinal direction of the X-ray source 12 FPD18 (Y direction), and moves in the y theta direction displaced, resulting in that state It may be reconstructed using the captured image. In this case, a highly accurate composite image (tomosynthesis image) cannot be obtained. Therefore, in the present embodiment, calibration can be performed as follows to correct a shift in the movement direction of the X-ray source 12 with respect to the FPD 18.
  • the X-ray source 12 is moved in the moving direction (y direction) by the moving unit 14, and at least two images each including the marker 20 are included. Take the picture as you like. In this case, for example, as shown in FIG. 3, a first image 40 and a second image 42 are obtained. In this case, if the movement direction (y direction) of the X-ray source 12 and the length direction (Y direction) of the FPD 18 are parallel, the position of the marker 20 is orthogonal to the length direction (Y direction) of the FPD 18.
  • the position in the X direction is the same between the first image 40 and the second image 42.
  • the position in the X direction is, for example, ⁇ between the first image 40 and the second image 42. Just shift. That is, it is shifted by ⁇ .
  • the shift angle ⁇ of the X-ray source 12 can also be obtained from the shift amount ⁇ and the distance d between the imaging positions of the first image 40 and the second image 42.
  • the shift angle ⁇ (first inclination) of the X-ray source 12 can be obtained using the relative relationship between the positions of the captured images of the markers 20.
  • Reference numeral 44 a indicates an outer edge of the composite image 44.
  • the shift angle ⁇ (first inclination) refers to the above-described inclination and is relative to the arrangement of a plurality of pixels (two-dimensional matrix arrangement direction).
  • the relationship between the variation rate of the size of the marker 20 in the captured image and the distance L between the X-ray source and the FPD is obtained in advance, and stored in the storage unit 22a, for example.
  • the distance L between the X-ray source and the FPD can be obtained.
  • the shift angle (second inclination) of the X-ray source 12 in the z direction can be obtained using the relative relationship of the positions of the captured images of the markers 20.
  • the calculation of the shift amount ⁇ and the shift angle ⁇ of the X-ray source 12 and the shift amount and the shift angle in the z direction is a state where the subject M (patient) is not on the surface 16a of the imaging table 16. However, it may be performed in the imaging state of the subject M (patient). Further, in order to enable correction when the movement axis of the X-ray source 12 is not straight, that is, when the rail 13 is not parallel to the longitudinal direction (Y direction) of the FPD 18, The angle ⁇ and the shift amount and shift angle in the z direction are obtained for each position of the X-ray source 12, stored in the storage unit 22a, for example, in the form of a table, and used for tomosynthesis imaging. Also good.
  • tomosynthesis imaging of the X-ray imaging apparatus 10 will be described. First, after the subject M is positioned on the surface 16 a of the imaging table 16, when an instruction to start imaging is given from the input unit 28, tomosynthesis imaging is started under the control of the control unit 30.
  • the moving unit 14 moves the X-ray source 12 in the y direction, changes the irradiation angle of the X-ray source 12 in the direction of the subject M at a predetermined imaging position, and at different irradiation angles.
  • X-rays are irradiated onto the subject M, and a plurality of projection data (X-ray images) with different imaging angles are obtained by one imaging operation.
  • projection data corresponding to the captured X-ray image is output from the FPD 18 to the storage unit 22a of the image processing unit 22, and a plurality of images acquired by the FPD 18 are obtained.
  • the projection data of the X-ray image is stored in the storage unit 32.
  • the correction position 22b corrects the imaging position of the projection data (position information of the X-ray source 12 (X-ray tube position information)) for those having a positional deviation.
  • the corrected imaging position position information of the X-ray source 12 (X-ray tube position information)
  • the reconstruction unit 22c uses the projection data and the corrected imaging position (position information of the X-ray source 12 (X-ray tube position information)) to obtain an X-ray tomogram in a cross section at an arbitrary height of the subject M.
  • the image is reconstructed. In this way, since the displacement of the photographing position is corrected, a highly accurate composite image (tomosynthesis image) can be obtained for a cross section of the subject M at a predetermined position.
  • the obtained tomosynthesis image (X-ray tomographic image) is displayed on the display unit 24, for example. Further, the tomosynthesis image (X-ray tomographic image) is output to the output unit 26, and for example, the tomosynthesis image (X-ray tomographic image) is printed out and the digital image data of the X-ray tomographic image is recorded on the recording medium. Saved.
  • FIGS. 6A and 6B are conceptual diagrams showing a state when an X-ray tomographic image is reconstructed by tomosynthesis imaging.
  • the X-ray source 12 starts from the position S1 and moves to S3, and the subject M is irradiated with radiation at each of the source positions S1, S2, and S3. It is assumed that M X-ray images P1, P2, and P3 are obtained.
  • the photographing objects A and B exist at two positions where the height of the subject M is different.
  • S1, S2, and S3 the imaging position of the X-ray source 12 at the time of imaging
  • the X-rays emitted from the X-ray source 12 pass through the subject M and enter the FPD 18.
  • the two imaging objects A and B are projected with different positional relationships.
  • the imaging objects A and B are respectively projected to the positions of P1A and P1B which are shifted to the right side of the photographing objects A and B.
  • the projection is performed at the positions of P2A and P2B almost immediately below, and in the case of the X-ray image P3, the projection is performed at the positions of P3A and P3B shifted to the left side.
  • the projection positions P1A, P2A, and P3A of the imaging object A match. For example, as shown in FIG. 6B, the X-ray image P1 is shifted to the left, and the X-ray image P3 is shifted to the right for synthesis. Thereby, an X-ray tomographic image having a height at which the imaging object A exists is reconstructed. Similarly, an X-ray tomographic image in a cross section having an arbitrary height can be reconstructed.
  • a shift addition method is typically used as a method for reconstructing an X-ray tomographic image.
  • the shift addition method the positions of a plurality of photographed projection data (X-ray images) are sequentially shifted and added based on the position of the X-ray source 12 at the time of photographing each X-ray image.
  • FIG. 7 is a schematic diagram showing an X-ray imaging apparatus according to the second embodiment of the present invention
  • FIG. 8 is a flowchart showing a calibration method in the X-ray imaging apparatus according to the second embodiment of the present invention. It is.
  • the same components as those of the first embodiment shown in FIGS. 1 to 5 are denoted by the same reference numerals, and detailed description thereof is omitted.
  • the X-ray imaging apparatus 10 a of the second embodiment shown in FIG. Except for the difference in operation and calibration method, the configuration and operation are the same as those of the X-ray imaging apparatus 10 of the first embodiment, and detailed description thereof will be omitted.
  • the determination unit 23 is connected to the control unit 30.
  • the control unit 30 captures X-ray images by the X-ray source 12, the moving unit 14, and the X-ray source control unit 15, image processing in the image processing unit 22, and determination unit 23. It controls operations of the X-ray imaging apparatus 10a such as determination of necessity of calibration update, screen display on the display unit 24, and output processing on the output unit 26. Further, based on the determination result of the determination unit 23, at least the first calibration data and the second calibration data described later are transmitted to the image processing unit 22 (correction unit 22 b described later) via the control unit 30. 1 calibration data is created, and the created calibration data is stored in the image processing unit 22 (a storage unit 22a described later). In the present embodiment, the storage unit 22a, the correction unit 22b, and the control unit 30 constitute a calibration unit.
  • the determination unit 23 determines whether or not the calibration data needs to be updated in the calibration unit. In the determination unit 23, as will be described later, out of the position of the marker 20 (marker image) in the captured image obtained by the calibration unit (correction unit 22b) and the size of the marker 20 (marker image) in the captured image. The necessity of updating the calibration data is determined based on at least the position of the marker 20 (marker image) in the captured image. That is, the determination unit 23 determines whether or not to update the calibration data already stored in the storage unit 22a by causing the calibration unit to perform calibration.
  • the determination unit 23 also obtains the position of the image of the corresponding marker corresponding to the image of the marker 20 in the second image 42 using the first calibration data currently stored. Furthermore, the size of the image of the corresponding marker corresponding to the image of the marker 20 in the second image 42 is obtained using the second calibration data currently stored.
  • the marker 16 is provided with, for example, one marker 20 within the range where the FPD 18 is provided, as in the first embodiment.
  • the marker 20 is provided with an edge aligned with the outer edge 18 b of the FPD 18. Further, the marker 20 is provided on the imaging table 16 with a gap from the surface 18a of the FPD 18. Further, the marker 20 may be detachable from the imaging table 16. In this way, the calibration can be easily updated by making the marker 20 detachable.
  • the storage unit 22a receives a plurality of projection data of X-ray images acquired by the FPD 18 and stores the projection data of these X-ray images.
  • the storage unit 22a also stores a design size of the marker 20 in the projection data, a position of an image of the marker 20 to be described later, a threshold value regarding the degree of coincidence of the size of the marker 20 and the like.
  • the storage unit 22a stores first calibration data and second calibration data described later. The first calibration data and the second calibration data are appropriately updated according to the determination result of the determination unit 23 as described above.
  • the correction unit 22b calculates first calibration data, which will be described in detail later, at each imaging position of the X-ray source 12. Further, as will be described later, the second calibration data is calculated.
  • the correction unit 22b has a position detection function for detecting the image position of the marker 20 in the captured image, and a marker image is obtained from each of the first image and the second image captured from different X-ray tube positions. Is automatically used to generate first calibration data, and is used for determining whether or not the first calibration data needs to be updated in the determination unit 23.
  • the determination unit 23 obtains the corresponding marker image position where the marker image of the second image should be from the marker image position detected from the first image and the first calibration data, and this is detected from the second image.
  • the amount of deviation is obtained by comparing with the marker image position.
  • the necessity of calibration is determined by comparing the deviation amount with a preset threshold value.
  • the correction unit 22b calculates the first calibration data and updates the calibration data according to the determination result of whether or not to update the calibration data of the determination unit 23 described above.
  • First calibration data is created using at least two captured images.
  • the first calibration data is created using the first image 40 and the second image 42.
  • the position of the marker 20 is orthogonal to the length direction (Y direction) of the FPD 18.
  • the position in the X direction is the same between the first image 40 and the second image 42. That is, the shift amount is zero.
  • X-ray source 12 is displaced in the y theta direction as shown in FIG. 2 as in the first embodiment If the object moves, the position of the image of the marker 20 is shifted by, for example, ⁇ between the first image 40 and the second image 42 in the X direction. That is, the shift amount ⁇ .
  • the correction unit 22b creates first calibration data based on the deviation angle ⁇ (first inclination) and stores the first calibration data in the storage unit 22a.
  • the movement axis of the X-ray source 12 is inclined based on a design value (a target value at which each of a plurality of images is taken) by an amount corresponding to the deviation angle ⁇ .
  • the deviation in the y ⁇ direction shown in FIG. 2 can be corrected.
  • the first image 40 and the second image 42 are synthesized by combining the X-direction ends of the photographed image to obtain a synthesized image 44.
  • the shift amount ⁇ of the marker 20 is obtained using the outer edge 44a of the composite image 44 as a reference, and the X-ray source is determined by the shift amount ⁇ and the distance d between the imaging positions of the first image 40 and the second image 42. Twelve deviation angles ⁇ (first inclination) can also be obtained.
  • the shift angle ⁇ (first inclination) of the X-ray source 12 can be obtained using the relative relationship between the positions of the captured images of the markers 20.
  • the first calibration data described above can also be created based on the deviation angle ⁇ (first inclination).
  • the deviation angle ⁇ (first inclination) is the same as that described in the first embodiment, and a detailed description thereof will be omitted.
  • the distance L between the X-ray source and the FPD changes from the set distance.
  • the size of the marker 20 in the captured image changes.
  • the size of 20 and the enlargement rate or reduction rate of the marker 20 are obtained and stored in the storage unit 22a, for example.
  • a threshold is set for the size of the marker 20 in the captured image and the enlargement ratio or reduction ratio of the marker 20, and the threshold is stored in the storage unit 22a.
  • the correction unit 22b automatically detects the size of the marker image from each of the first image and the second image taken from different X-ray tube positions in addition to the position of the image of the marker 20,
  • the determination unit 23 can also determine whether the second calibration data needs to be updated. In this case, the determination unit 23 obtains the size of the image of the marker 20 in the first image 40 and the size of the image of the marker 20 in the second image 42.
  • the size of the corresponding marker image that should be the marker image of the second image 42 is obtained from the second calibration data, and is compared with the size of the marker 20 of the second image 42 to obtain a deviation amount.
  • the necessity of calibration is determined by comparing the deviation amount with a preset threshold value. As the amount of deviation, the difference in the diameter of the image of the marker 20, the enlargement rate or reduction rate of the image of the marker 20, and the like can be used.
  • Second calibration data is created using at least two captured images. For example, as a result of photographing at two photographing positions, when the size of the marker 20 shown in FIG. 5A is changed like the marker 21 shown in FIG. 5B as in the first embodiment, the correction unit In 22b, the distance L between the X-ray source and the FPD at each imaging position can be obtained. Further, in the correction unit 22b, a deviation angle (second inclination) in a plane orthogonal to the surface of the FPD 18 of the X-ray source 12 shown in FIG.
  • the shift angle (second inclination) of the X-ray source 12 in the z direction can be obtained using the relative relationship of the positions of the captured images of the markers 20.
  • the correction unit 22b creates second calibration data based on the deviation angle (second inclination) and stores the second calibration data in the storage unit 22a.
  • the second calibration data for example, enlargement rate data or reduction rate for enlarging or reducing the entire projection data so that the size of the marker 20 becomes a design size corresponding to the enlargement rate or reduction rate. Data, etc.
  • the deviation in the z direction shown in FIG. 7 can be corrected by the second calibration data.
  • the threshold value of the degree of coincidence of the image position of the marker 20 and the threshold value of the degree of coincidence of the image size of the marker 20 may be set via the input unit 28, for example.
  • the calibration data (first calibration data and second calibration data) created by the correction unit 22b may be set via the input unit 28, for example.
  • the first calibration data and the second calibration data are obtained for each position of the X-ray source 12, and stored in the storage unit 22a, for example, in the form of a table. It may be used for tomosynthesis photography.
  • the reconstruction unit 22c corrects the projection data of a plurality of X-ray images including the projection data after correction processing by the correction unit 22b performed as necessary, and corrects the positional information of the X-ray source 12 as necessary.
  • An image composition process is performed using (X-ray tube position information) to reconstruct an X-ray tomographic image of a cross section of the subject M at an arbitrary height.
  • the first calibration data and the second calibration data may be created.
  • a highly accurate tomosynthesis image can be obtained, and the first calibration data and the second calibration data are obtained.
  • a tomosynthesis image with higher accuracy can be obtained.
  • the reconstruction unit 22c reconstructs an X-ray tomographic image of a cross section at an arbitrary height of the subject M, the X-ray corrected by the correction unit 22b using at least the first calibration data.
  • the position information of the source 12 (X-ray tube position information) may be used.
  • the image processing unit 22 (the storage unit 22a, the correction unit 22b, and the reconstruction unit 22c) is not limited to being configured by hardware (device).
  • a program for causing a computer to execute processing by the correction unit 22b and the reconstruction unit 22c of the image processing unit 22 on the projection data stored in the storage unit 22a may be used.
  • the display unit 24 has the same configuration as that of the first embodiment, but further functions as a notification unit that notifies the determination result of whether the determination unit 23 needs to update calibration.
  • the control unit 30 displays the determination result of whether or not the calibration unit needs to be updated on the display unit 24 in the form of characters or the like.
  • the reconstructed image may be reconstructed using the captured image obtained in that state.
  • a highly accurate composite image tomosynthesis image
  • the calibration data may not be appropriate due to deterioration over time (change over time).
  • the necessity of calibration is determined as shown in FIG. 8, and the influence of deterioration with time (change with time) of the calibration data is suppressed.
  • the calibration of the X-ray imaging apparatus 10a is preferably performed in a state where the subject M (patient) is not present on the surface 16a of the imaging table 16, but may be performed in the imaging state of the subject M (patient).
  • the X-ray imaging apparatus 10a has already been calibrated, and the first calibration data and the second calibration data are stored in the storage unit 22a as initial calibration data (step S10).
  • the moving unit 14 moves the X-ray source 12 to at least two imaging positions, and the marker 20 is imaged from the X-ray source 12 at each imaging position.
  • X-rays are irradiated onto the FPD 18, and each image data of at least two captured images (for example, the first image 40 and the second image 42) including the marker 20 from the FPD 18 is stored in the image processing unit 22. Is output to the unit 22a.
  • the correction unit 22b for example, the position of the marker 20 image of the first image 40 and the position of the marker 20 image of the second image 42 are obtained.
  • the position of the corresponding marker image corresponding to the image of the marker 20 in the second image 42 is obtained using the first calibration data currently stored in the determination unit 23. Then, in the correction unit 22b, a deviation amount between the position of the marker 20 in the second image 42 and the position of the corresponding marker image obtained by using the first calibration data is obtained (step S12).
  • the determination unit 23 determines whether or not calibration (creation of first calibration data) is necessary based on the threshold value of the deviation amount at the image position of the marker 20 stored in the storage unit 22a (Ste S14). Note that the determination result regarding the necessity of calibration related to the first calibration data is displayed on the display unit 24 in the form of characters or the like, and is notified to the user or the like.
  • step S14 if the amount of deviation at the position of the image of the marker 20 is within the threshold, the first calibration data stored in the storage unit 22a is updated to the X-ray tomographic image without updating the initial calibration data. It is used for reconstruction.
  • step S18 first calibration data is created again (step S18).
  • step S12 for example, the size of the image of the marker 20 in the first image 40 and the size of the image of the marker 20 in the second image 42 are obtained. Further, the size of the corresponding marker image corresponding to the image of the marker 20 in the second image 42 is obtained using the second calibration data currently stored in the determination unit 23. Then, the correction unit 22b may obtain a deviation amount between the size of the marker 20 of the second image 42 and the size of the corresponding marker image obtained using the second calibration data (step S12). . In this case, the deviation amount of the image size of the marker 20 is then output to the determination unit 23.
  • the determination unit 23 determines whether or not calibration (creation of second calibration data) is necessary based on the threshold value of the deviation amount of the image size of the marker 20 stored in the storage unit 22a (step S14). ). Note that the determination result regarding the necessity of calibration regarding the second calibration data is displayed on the display unit 24 in the form of characters or the like, and is notified to the user or the like.
  • step S14 if the amount of deviation of the image size of the marker 20 is within the threshold value, the second calibration data stored in the storage unit 22a is updated to the X-ray tomography without updating the initial calibration data. Used for image reconstruction. On the other hand, if the amount of deviation of the image size of the marker 20 exceeds the threshold value in step S14, second calibration data is created again (step S18).
  • the first calibration data and the second calibration data are created by photographing an image including the marker 20 from at least two photographing positions and using each obtained photographed image data. This is a method, and a detailed description of the creation method is omitted.
  • step S20 the initial calibration data is updated.
  • the initial calibration data is updated based on the determination result. For this reason, even when the moving direction of the X-ray source 12 is shifted due to a change with time and the initial calibration data is deteriorated with time (change with time), it is possible to obtain appropriate calibration data.
  • the movement axis of the X-ray source 12 is not straight, that is, when the rail 13 is not parallel to the longitudinal direction (Y direction) of the FPD 18 (movement direction of the X-ray source 12 (y direction) ) And the pixel arrangement direction of the FPD 18 are not parallel to each other), and even if the movement axis of the X-ray source 12 is deviated, the deviation can be corrected and the influence of the change in the deviation of the movement axis of the X-ray source 12 over time.
  • the tomosynthesis can always be reconstructed with high accuracy by always using appropriate calibration data, and a high-accuracy composite image (tomosynthesis image) is always obtained for the cross-section at a predetermined position of the subject. Obtainable.
  • tomosynthesis imaging of the X-ray imaging apparatus 10a will be described. First, after the subject M is positioned on the surface 16 a of the imaging table 16, when an instruction to start imaging is given from the input unit 28, tomosynthesis imaging is started under the control of the control unit 30.
  • the moving unit 14 moves the X-ray source 12 in the y direction, changes the irradiation angle of the X-ray source 12 in the direction of the subject M at a predetermined imaging position, and at different irradiation angles.
  • X-rays are irradiated onto the subject M, and a plurality of projection data (X-ray images) with different imaging angles are obtained by one imaging operation.
  • projection data corresponding to the captured X-ray image is output from the FPD 18 to the storage unit 22a of the image processing unit 22, and a plurality of images acquired by the FPD 18 are obtained.
  • the projection data of the X-ray image is stored in the storage unit 32.
  • the first calibration data and the second calibration data described above are stored in advance in the storage unit 22a for each photographing position. Further, as in the method shown in FIG. 8, it is determined whether or not the calibration data needs to be updated before shooting, and the calibration data is updated according to the determination result. In addition, the determination result of necessity of this calibration is displayed on the display unit 24 in the form of characters or the like, for example.
  • the correction unit 22b uses the imaging position of the projection data (the position information of the X-ray source 12 (the X-ray tube) for those that are misaligned based on the first calibration data and the second calibration data. The position information)) is corrected.
  • the corrected imaging position position information of the X-ray source 12 (X-ray tube position information)
  • the reconstruction unit 22c uses the projection data and the corrected imaging position (position information of the X-ray source 12 (X-ray tube position information)) to obtain an X-ray tomogram in a cross section at an arbitrary height of the subject M.
  • the image is reconstructed. In this way, since the displacement of the photographing position is corrected, a highly accurate composite image (tomosynthesis image) can be obtained for a cross section of the subject M at a predetermined position.
  • the obtained tomosynthesis image (X-ray tomographic image) is displayed on the display unit 24, for example. Further, the tomosynthesis image (X-ray tomographic image) is output to the output unit 26, and for example, the tomosynthesis image (X-ray tomographic image) is printed out and the digital image data of the X-ray tomographic image is recorded on the recording medium. Saved.
  • the present invention is basically configured as described above. As described above, the X-ray imaging apparatus and the calibration method thereof according to the present invention have been described in detail. However, the present invention is not limited to the above-described embodiment, and various improvements or modifications can be made without departing from the gist of the present invention. Of course it is also good.

Abstract

　Ｘ線撮影装置は、所定の移動経路で移動可能なＸ線源と、Ｘ線源を所定の移動経路で移動させる移動手段と、Ｘ線源に対向して設けられた撮影台と、撮影台に設けられた平板状のＸ線検出器と、マーカと、Ｘ線源を移動させて少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させる制御部と、撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいてＸ線源の移動軸の、Ｘ線検出器に対する傾きを求める画像処理部とを有する。

Description

Ｘ線撮影装置およびそのキャリブレーション方法
　本発明は、臥位、立位、マンモ等で撮影を行うＸ線撮影装置およびそのキャリブレーション方法に関し、特に、Ｘ線検出器に対してＸ線源の移動軸がずれて設置されていても、そのずれを補正でき、被写体の所定の位置の断面について、高精度な合成画像（トモシンセシス画像）を得ることができ、更にはそのずれの経時劣化も抑制することができるＸ線撮影装置およびそのキャリブレーション方法に関する。
　現在、患部をより詳しく観察するために、Ｘ線源を移動させて異なる角度から被写体にＸ線を照射して撮影を行い、得た画像を加算して所望の断層面を強調した画像を得ることができるトモシンセシス撮影が可能なＸ線撮影装置が種々提案されている（例えば、特許文献１参照）。
　トモシンセシス撮影では、例えば、Ｘ線源をＸ線検出器と平行に移動させたり、円、楕円等の弧を描くように移動させて、異なる照射角で被写体を撮影した複数の撮影画像を取得して、これらの撮影画像を再構成して断層画像を作成する。
　特許文献１には、トモシンセシス撮影の際に、放射線源のみを移動させる第１のモードと、放射線源および検出手段の双方を移動させる第２のモードとのいずれかを撮影条件に応じて選択することができ、選択されたモードにしたがって複数の撮影画像を取得する放射線撮影装置が開示されている。
特開２０１１－６７５０３号公報 特開２０００－２７８６０６号公報 特開２０１０－２５２９５１号公報
　しかしながら、トモシンセシス撮影には、Ｘ線検出器（Ｘ線画像センサ）とＸ線源の移動軸のアライメントに高い精度が要求されている。特に、天井走行型のＸ線源を用いたＸ線撮影装置では、Ｘ線源とＸ線検出器（Ｘ線画像センサ）が一体化されておらず、高い精度でアライメントが取れるように設置するのは容易ではない。
　従来、トモシンセシス撮影においては、既定のマーカーアッセンブリを用いてキャリブレーションデータを生成する方法が知られている。しかし、この方法では、マーカの精度管理を要する点、ＳＩＤ（Ｘ線管焦点・受像面間距離）、Ｘ線源の角度、ショット数等の撮影条件毎にキャリブレーションデータを生成する必要がある点等の難点がある。
　さらには、既定のマーカーアッセンブリを用いてキャリブレーションデータを生成する方法では、事前に取得したキャリブレーションデータを用いるため、経時劣化（経時変化）に対応することができないという問題点もある。
　また、上述のもの以外にも、マーカを用いて画像の位置ずれを検出して、Ｘ線画像を補正することがなされている（例えば、特許文献２参照）。
　特許文献２には、被検体の撮影部位にマーカを有するサポータを装着し、この撮影部位のマスク像とライブ像とを撮影し、マスク像のマーカの位置とライブ像のマーカの位置との誤差を計算してこれが一致するようにしてサブトラクションを行い、カテタリゼーションテーブルの位置ずれにより生じるサブトラクション画像のノイズ発生を防止するＸ線画像撮影装置が開示されている。
　この場合においても、キャリブレーションデータを生成するために、マーカの精度管理を要する点、ＳＩＤ（Ｘ線管焦点・受像面間距離）、Ｘ線源の角度、ショット数等の撮影条件毎にキャリブレーションデータを生成する必要がある点等の難点がある。
　また、特許文献３には、臨床現場における運用に適し、リングアーチファクトのない精度の高いＣＴ画像を取得することを目的とするＸ線ＣＴ装置が開示されている。この特許文献３のＸ線ＣＴ装置においては、出力された画像にリングアーチファクトが検出された場合、エアキャリブレーションデータ取得部は新しいエアキャリブレーションデータを取得し、エアキャリブレーションデータ交換部は、その時点まで使用していたエアキャリブレーションデータを新しいエアキャリブレーションデータに差し替える。
　このように特許文献３は、エアキャリブレーションデータを差し替えることが開示されているものの、Ｘ線ＣＴ装置に関するものであり、トモシンセシス撮影において、キャリブレーションの経時劣化に対応できるものがないのが現状である。
　本発明の目的は、前記従来技術に基づく問題点を解消し、Ｘ線検出器に対してＸ線源の移動軸がずれて設置されていても、そのずれを補正でき、被写体の所定の位置の断面について、高精度な合成画像（トモシンセシス画像）を得ることができ、更にはそのずれを補正するとともに、そのずれの経時劣化も抑制することができるＸ線撮影装置およびそのキャリブレーション方法を提供することにある。
　上記目的を達成するために、本発明の第１の態様は、所定の移動経路で移動可能なＸ線源と、Ｘ線源を所定の移動経路で移動させる移動手段と、Ｘ線源に対向して設けられた撮影台と、撮影台に設けられた平板状のＸ線検出器と、撮影台に設けられるマーカと、移動手段によりＸ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させる制御部と、撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、Ｘ線源の移動軸の、Ｘ線検出器に対する傾きを求める画像処理部とを有することを特徴とするＸ線撮影装置を提供するものである。
　画像処理部は、Ｘ線源の撮影位置毎に、Ｘ線検出器に対する傾きを求めることが好ましい。また、画像処理部は、Ｘ線検出器に対する傾きを、断層画像の再構成に用いることが好ましい。
　Ｘ線検出器に対する傾きは、例えば、Ｘ線検出器の表面と平行な面内での第１の傾きである。
　画像処理部は、さらに撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いてＸ線源とＸ線検出器との距離を求めて、Ｘ線検出器に対する傾きとして、Ｘ線検出器の表面と直交する面内での第２の傾きを求めることが好ましい。
　本発明の第２の態様は、Ｘ線撮影装置のキャリブレーション方法であって、Ｘ線撮影装置は、所定の移動経路で移動可能なＸ線源に対向して設けられた撮影台と、撮影台に設けられた平板状のＸ線検出器と、マーカとを備えており、Ｘ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させる工程と、撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、Ｘ線源の移動軸の、Ｘ線検出器に対する傾きを求める工程を有することを特徴とするキャリブレーション方法を提供するものである。
　Ｘ線源の撮影位置毎に、Ｘ線検出器に対する傾きを求めることが好ましい。
　また、Ｘ線検出器に対する傾きは、例えば、Ｘ線検出器の表面と平行な面内での第１の傾きである。
　さらに撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いてＸ線源とＸ線検出器との距離を求めて、Ｘ線検出器に対する傾きとして、Ｘ線検出器の表面と直交する面内での第２の傾きを求める工程を有することが好ましい。
　本発明の第３の態様は、所定の移動経路で移動可能なＸ線源と、Ｘ線源を所定の移動経路で移動させる移動手段と、Ｘ線源に対向して設けられた撮影台と、撮影台に設けられた平板状のＸ線検出器と、撮影台に設けられるマーカと、Ｘ線検出器に対して配置されたＸ線源の移動方向のずれに関するキャリブレーションを行うキャリブレーション部と、キャリブレーション部によるキャリブレーションの更新の要否を判定し、判定結果に基づいてキャリブレーション部にキャリブレーションを行わせる判定部とを有し、Ｘ線源によりマーカを含む画像が撮影され、キャリブレーション部により撮影された画像のマーカ像の位置が求められ、判定部は、このマーカ像の位置に基づいてキャリブレーションの更新の要否を判定することを特徴とするＸ線撮影装置を提供するものである。
　キャリブレーション部により、Ｘ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させ、撮影された第１の画像、第２の画像についてマーカ像の位置が求められ、判定部は、既存のキャリブレーションデータを用いて第１の画像のマーカ像の位置から第２の画像のマーカ像に対応する対応マーカ像の位置を求め、第２の画像のマーカ像の位置と、第２の画像のマーカ像に対応する対応マーカ像の位置とのずれに基づいてキャリブレーションの更新の要否を判定することが好ましい。
　キャリブレーション部は、Ｘ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させ、撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、Ｘ線源の移動軸の、Ｘ線検出器の表面と平行な面内での第１の傾きを求め、第１の傾きに基づいて第１のキャリブレーションデータを作成するものであることが好ましい。
　キャリブレーション部は、Ｘ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させ、撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いてＸ線源とＸ線検出器との距離を求めて、Ｘ線検出器の表面と直交する面内でのＸ線源の第２の傾きを求め、第２の傾きに基づいて第２のキャリブレーションデータを作成するものであることが好ましい。
　また、キャリブレーション部は、Ｘ線源の撮影位置毎に、第１のキャリブレーションデータおよび第２のキャリブレーションデータのうち、少なくとも第１のキャリブレーションデータを求めることが好ましい。
　さらに、移動手段によりＸ線源が所定の撮影位置に移動されトモシンセシス撮影されて得られた複数の画像を用いて断層画像を再構成する再構成部を有し、再構成部は、トモシンセシス撮影で撮影された複数の画像を用いて断層画像を再構成する際に、第１のキャリブレーションデータおよび第２のキャリブレーションデータのうち、少なくとも第１のキャリブレーションデータを用いることが好ましい。
　さらに、判定部によるキャリブレーションの更新の要否の判定結果を通知する通知部を有することが好ましい。
　本発明の第４の態様は、Ｘ線撮影装置のキャリブレーション方法であって、Ｘ線撮影装置は、所定の移動経路で移動可能なＸ線源に対向して設けられた撮影台と、撮影台に設けられた平板状のＸ線検出器と、マーカとを備えており、Ｘ線源によりマーカを含む画像を撮影させる工程と、撮影された画像のマーカ像の位置を求め、このマーカ像の位置に基づいて、キャリブレーションの更新の要否を判定し、判定結果に基づいてキャリブレーションを行う工程とを有することを特徴とするキャリブレーション方法を提供するものである。
　キャリブレーション工程は、Ｘ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させ、撮影された第１の画像、第２の画像についてそれぞれマーカ像の位置を求める工程と、既存のキャリブレーションデータを用いて第１の画像のマーカ像の位置から第２の画像のマーカ像に対応する対応マーカ像の位置を求める工程と、第２の画像のマーカ像の位置と、第２の画像のマーカ像に対応する対応マーカ像の位置とのずれに基づいてキャリブレーションの更新の要否を判定する工程とを有することが好ましい。
　キャリブレーション工程は、Ｘ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させ、撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、Ｘ線源の移動軸の、Ｘ線検出器の表面と平行な面内での第１の傾きを求め、第１の傾きに基づいて第１のキャリブレーションデータを作成する工程を有することが好ましい。
　キャリブレーション工程は、Ｘ線源を移動させて、少なくとも２つの位置から、それぞれマーカを含む画像を撮影させ、撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いてＸ線源とＸ線検出器との距離を求めて、Ｘ線検出器の表面と直交する面内でのＸ線源の第２の傾きを求め、第２の傾きに基づいて第２のキャリブレーションデータを作成するものであることが好ましい。
　また、Ｘ線源の撮影位置毎に、第１のキャリブレーションデータおよび第２のキャリブレーションデータのうち、少なくとも第１のキャリブレーションデータを求めることが好ましい。
　本発明によれば、Ｘ線検出器に対してＸ線源の移動軸がずれて設置されていても、そのずれを補正することができる。このため、トモシンセシスの再構成の際に、高い精度で各画像を用いた再構成をすることができる。これにより、被写体の所定の位置の断面について、高精度な合成画像（トモシンセシス画像）を得ることができる。
　また、本発明によれば、Ｘ線検出器とＸ線源の移動軸とのずれを校正するためのキャリブレーションデータが、経時変化により、更新する必要があるか否かについて容易に判定することができる。これにより、キャリブレーションの要否を容易に知ることができ、必要に応じてキャリブレーションを行うことができる。このため、経時変化の影響を抑制して、トモシンセシスの再構成を、常に高い精度で行うことができ、被写体の所定の位置の断面について、高精度な合成画像（トモシンセシス画像）を常に得ることができる。
本発明の第１の実施形態に係るＸ線撮影装置を示す模式図である。 本発明の第１の実施形態に係るＸ線撮影装置のＸ線源とＸ線検出器との配置状態を示す模式図である。 本発明の第１の実施形態に係るＸ線撮影装置で撮影された２つの画像を示す模式図である。 本発明の第１の実施形態に係るＸ線撮影装置で撮影された２つの画像の合成画像を示す模式図である。 （ａ）および（ｂ）は、本発明の第１の実施形態に係るＸ線撮影装置で撮影されたマーカのマーカ像の変動を示す模式図である。 （ａ）および（ｂ）は、トモシンセシス撮影によるＸ線断層画像の再構成時の様子を示す概念図である。 本発明の第２の実施形態に係るＸ線撮影装置を示す模式図である。 本発明の第２の実施形態に係るＸ線撮影装置におけるキャリブレーション方法を示すフローチャートである。
　以下に、添付の図面に示す好適実施形態に基づいて、本発明のＸ線撮影装置およびそのキャリブレーション方法を詳細に説明する。
　図１は、本発明の第１の実施形態に係るＸ線撮影装置を示す模式図である。
　図１に示す第１の実施形態のＸ線撮影装置１０は、人体等の被写体Ｍをトモシンセシス撮影して、被写体Ｍの任意高さの断面におけるＸ線断層画像を得るものである。
　Ｘ線撮影装置１０は、例えば、臥位、立位、マンモのトモシンセシス撮影に用いることができる。
　Ｘ線撮影装置１０は、Ｘ線検出器１８は、ビルトインを用いることができる。また、Ｘ線撮影装置１０は、カセッテの位置がキャリブレーション前後で変わらなければ、カセッテタイプのＸ線検出器１８を用いることができる。
　Ｘ線撮影装置１０は、天井走行型のＸ線源１２、移動部１４、Ｘ線源制御部１５、撮影台１６、Ｘ線検出器１８、マーカ２０、画像処理部２２、表示部２４、出力部２６および制御部３０を有する。制御部３０は、移動部１４、Ｘ線源制御部１５、Ｘ線検出器１８、画像処理部２２、表示部２４、出力部２６および入力部２８に接続されている。
　制御部３０は、入力部２８から入力される指示入力信号に従って、Ｘ線源１２、移動部１４およびＸ線源制御部１５によるＸ線画像撮影、画像処理部２２における画像処理、表示部２４における画面表示、および出力部２６における出力処理等、Ｘ線撮影装置１０の動作を制御するものである。
　入力部２８は、撮影開始の指示および後述する切替指示等を含む、各種の指示を入力するための部位であり、例えば、マウスおよびキーボード等を例示することができる。入力部２８を介して指示入力が制御部３０に入力されて、Ｘ線撮影装置１０の各種の動作が制御部３０によりなされる。なお、入力部２８を介して入力される指示入力は、例えば、表示部２４に表示される。
　Ｘ線源１２は、撮影台１６に対向して天井Ｈに設けられたレール１３に、移動可能に取り付けられている。
　なお、レール１３は、図２に示すように、Ｘ線源１２の移動方向（ｙ方向）とＸ線検出器１８の長手方向（Ｙ方向）とが一致するように、すなわち、Ｘ線源１２の移動方向（ｙ方向）とＸ線検出器１８の長手方向（Ｙ方向）とが平行に配置されている。
　図１に示すように、レール１３には移動部１４が設けられている。この移動部１４は、例えば、Ｘ線源１２にベルト（図示せず）、プーリ（図示せず）を介して接続されたステッピングモータ等の駆動部（図示せず）を備える。この移動部１４により、Ｘ線源１２は、レール１３に沿って直線移動が可能であり、レール１３上の所定の位置に停止することもできる。
　なお、移動部１４は、Ｘ線源１２をレール１３に沿って直線移動させることができ、かつ所定の撮影位置にＸ線源１２を停止させることができれば、その構成は特に限定されるものではない。
　また、Ｘ線源１２の移動経路は、直線に限定されるものではなく、円弧状であってもよく、この場合、レール１３が円弧状に設けられる。
　Ｘ線源１２は、撮影台１６の表面１６ａ上の被写体ＭにＸ線（放射線）を照射するものであり、被写体Ｍの方向にＸ線の照射角度を変えることができるものである。Ｘ線源１２には、一般的なＸ線撮影装置に用いられる、Ｘ線の照射角度を変えることができるＸ線源を用いることができる。Ｘ線源１２には、例えば、図示はしないがコリメータ（照射野絞り）が設けられている。
　Ｘ線源１２は、Ｘ線源制御部１５に接続されている。このＸ線源制御部１５により、Ｘ線源１２によるＸ線の照射およびＸ線の照射角度が制御される。撮影条件に応じて、移動部１４によりＸ線源１２が所定の撮影位置に移動された後、Ｘ線源制御部１５により、被写体Ｍの方向にＸ線の照射角度を変えられた後、Ｘ線源１２から所定強度のＸ線を所定の時間だけ照射させる。すなわち、Ｘ線源１２から、所定の照射量（線量）のＸ線を被写体Ｍに向けて照射させる。
　撮影台１６は、表面１６ａに、人、動物などの被写体Ｍが載るものである。被写体Ｍとして、例えば、人が臥位でこの撮影台１６の表面１６ａに載る。
　また、撮影台１６の裏面１６ｂにフラットタイプ（平板状）のＸ線検出器１８（以下、ＦＰＤ１８という）が、レール１３に対向して設けられている。本実施形態においては、上述のように、ＦＰＤ１８の長手方向（Ｙ方向、図２参照）と、レール１３の長手方向、すなわち、Ｘ線源１２の移動方向（ｙ方向）とが平行になるように設置される。
　ＦＰＤ１８は、Ｘ線の受光面を上向きにして撮影台１６の裏面１６ｂに配設されている。ＦＰＤ１８は、被写体Ｍを透過したＸ線を検出して光電変換し、撮影された被写体ＭのＸ線画像に対応するデジタル画像データ（投影データ）を取得するものであり、検出器筐体の中に複数の画素が２次元マトリックス状に配置されてなる撮像面がある。ＦＰＤ１８から画像処理部２２にデジタル画像データが出力される。なお、上述のようにＸ線源１２の移動方向（ｙ方向）とＦＰＤ１８の長手方向（Ｙ方向）とが平行になるように配置されており、Ｘ線源１２の移動方向（ｙ方向）は、ＦＰＤ１８の画素配列方向と一致するように配置される。
　ここで、本発明において、“傾き”とは、複数の画素の並び（２次元マトリックスの配列方向）に対する傾きのことである。
　ＦＰＤ１８としては、Ｘ線（放射線）を電荷に直接変換する直接方式、Ｘ線（放射線）を一旦光に変換し、変換された光をさらに電気信号に変換する間接方式等、各種方式のものが利用可能である。また、Ｘ線源１２の移動方向に対してＦＰＤ１８を移動可能に構成してもよい。
　また、撮影台１６には、ＦＰＤ１８が設けられている範囲内にマーカ２０が、例えば、１個設けられている。このマーカ２０は、ＦＰＤ１８の表面１８ａに対して隙間をあけて撮影台１６に設けられている。また、マーカ２０は、ＦＰＤ１８の外縁１８ｂにエッジを一致させて設けてもよい。
　マーカ２０は、被写体Ｍのトモシンセシス撮影を行う前、例えば、装置起動時に撮影され、撮影された複数枚の投影データ（Ｘ線画像）のキャリブレーションのために用いられるものである。
　マーカ２０は、例えば、金属、好ましくは原子番号の大きい金属であって、放射線を透過しにくい材質のものである。マーカ２０の形状は、十字穴付き形状のように、マーカ２０の中心が分かりやすいものが望ましい。マーカ２０の数は、特に限定されるものではないが、複数あってもよい。複数の場合、各マーカの識別を容易にするために、マーカはそれぞれ異なることが好ましい。さらには、マーカ２０が複数の場合、例えば、４個設ける。なお、マーカ２０が複数の場合、その配置パターンは何ら限定されないが、例えば、矩形状に配置する。
　また、マーカ２０は、撮影台１６に対して着脱自在であってもよい。この場合、キャリブレーションのときにだけ、例えば、装置起動時等に、マーカ２０を撮影台１６上に適当に設置して撮影すれば、キャリブレーションすることができる。このように、マーカ２０を着脱自在にすることにより、容易にキャリブレーションすることができる。
　本実施形態のＸ線撮影装置１０において、トモシンセシス撮影を行う場合、移動部１４により、Ｘ線源１２をｙ方向に移動しつつ、所定の撮影位置で、被写体Ｍの方向にＸ線の照射角度を変えて、異なる撮影角度（一定の時間間隔）で被写体ＭにＸ線が照射される。Ｘ線源１２から照射されたＸ線は、被写体Ｍを透過してＦＰＤ１８の受光面に入射され、ＦＰＤ１８により検出されて光電変換され、撮影された被写体ＭのＸ線画像に対応する投影データ（デジタル画像データ）が取得される。
　トモシンセシス撮影の場合、１回の撮影操作により、被写体Ｍの、撮影角度の異なる複数枚（例えば、２０～８０枚）のＸ線画像が撮影され、ＦＰＤ１８から、撮影された複数枚のＸ線画像に対応する投影データ（デジタル画像データ）が順次画像処理部２２（記憶部２２ａ）に出力される。
　画像処理部２２は、ＦＰＤ１８で取得された複数枚のＸ線画像の投影データが入力されて、これら複数枚のＸ線画像の投影データを用いて画像処理（補正処理（キャリブレーション）、画像合成処理等を含む）を行って、被写体Ｍの任意高さの断面におけるＸ線断層画像を再構成する部位である。画像処理部２２は、記憶部２２ａ、補正部２２ｂ（キャリブレーション部）および再構成部２２ｃを有する。
　記憶部２２ａは、ＦＰＤ１８で取得されたＸ線画像の複数枚の投影データが入力されて、これらのＸ線画像の投影データを記憶するものである。
　また、記憶部２２ａには、後述するように、投影データにおけるマーカ２０の設計上の大きさ、ならびにずれ量（シフト量δ）の閾値等も記憶される。
　補正部２２ｂは、Ｘ線源１２の各撮影位置における位置のずれ量（Ｘ線源１２の位置情報）を算出し、その位置のずれ量が所定の値（例えば、閾値）よりも大きい場合、ずれ量が大きい撮影位置で撮影された投影データのＸ線源１２の位置情報（Ｘ線管球位置情報）を補正するものである。すなわち、トモシンセシス撮影で得られたＸ線画像の投影データの全てについて、ずれ量を算出し、その位置のずれ量が所定の値（例えば、閾値）よりも大きい場合、ずれ量が大きい撮影位置で撮影された投影データのＸ線源１２の位置情報（Ｘ線管球位置情報）を補正するものである。
　なお、ずれ量に対して、閾値を設定することなく、全ての投影データをずれ量に基づいてのＸ線源１２の位置情報（Ｘ線管球位置情報）を補正するようにしてもよい。
　なお、Ｘ線源１２の位置情報（Ｘ線管球位置情報）は、設計値（複数の画像をそれぞれどの位置で撮るかの目標値）をベースにしてＸ線源１２の移動軸が傾いている分を補正することによって、真のＸ線源１２の位置に近い情報に補正して、再構成演算に利用される。また、補正部２２ｂは、マーカ２０（マーカ像）の画像認識機能を有するものである。
　本実施形態において、Ｘ線源１２の移動軸のＦＰＤ１８（Ｘ線検出器）に対する傾きとは、Ｘ線源１２の撮影位置における位置のずれ量のことである。例えば、図２に示すＦＰＤ１８の長手方向（Ｙ方向）に対するＸ線源１２の移動方向（ｙθ方向）のずれ量である。更には、Ｘ線源１２の撮影位置における位置のずれ量には、図１に示すようにＸ線源１２とＦＰＤ１８の表面１８ａとの距離（Ｘ線源－ＦＰＤ間の距離Ｌ）のずれ量も含まれる。
　補正部２２ｂでは、所定の値（例えば、閾値）よりもずれ量が大きい場合、位置のずれ量に基づいて、例えば、ＦＰＤ１８の長手方向（Ｙ方向）に対するずれの場合には、ずれ量に応じた分だけ、Ｘ線源１２の撮影位置を補正する。
　また、補正部２２ｂでは、Ｘ線源１２とＦＰＤ１８の表面１８ａとの距離（Ｘ線源－ＦＰＤ間の距離Ｌ）のずれ量に対しては、Ｘ線源－ＦＰＤ間の距離Ｌが設定された距離から変化すると、撮影画像におけるマーカ２０（マーカ像）の大きさが変わる。このため、Ｘ線源－ＦＰＤ間の距離Ｌのずれ量に対して、マーカ２０（マーカ像）の変化率、すなわち、拡大率または縮小率を求めておき、これを、例えば、記憶部２２ａに記憶させておく。この拡大率または縮小率に応じて、投影データを拡大または縮小することにより、図１に示すｚ方向におけるずれを補正することができる。
　なお、補正部２２ｂにおいて、上述の位置のずれ量および距離のずれ量の両方に対して補正をしてもよいことはもちろんである。
　再構成部２２ｃは、必要に応じてなされた補正部２２ｂによる補正処理後の投影データを含む、複数枚のＸ線画像の投影データを、必要に応じて補正されたＸ線源１２の位置情報（Ｘ線管球位置情報）を用いて画像合成処理を施し、被写体Ｍの任意高さの断面におけるＸ線断層画像を再構成するものである。
　なお、画像処理部２２（記憶部２２ａ、補正部２２ｂ（キャリブレーション部）および再構成部２２ｃ）は、ハードウェア（装置）で構成することに限定されるものではない。例えば、記憶部２２ａに記憶された投影データに対して、画像処理部２２の補正部２２ｂ（キャリブレーション部）および再構成部２２ｃによる処理を、コンピュータに実行させるためのプログラムとすることもできる。
　表示部２４は、入力部２８から制御部３０に入力される指示入力信号の内容、およびこの指示入力信号に従って、画像処理部２２により再構成されたＸ線断層画像等を表示する部位であり、例えば、液晶ディスプレイ等のようなフラットパネル型ディスプレイを例示することができる。
　出力部２６は、入力部２８から制御部３０に入力される指示入力信号に従って、画像処理部２２により再構成されたＸ線断層画像を出力する部位であり、例えば、Ｘ線断層画像をプリント出力する各種のプリンタ、Ｘ線断層画像のデジタル画像データを各種の記録媒体に保存する記憶装置等を例示することができる。
　本実施形態のＸ線撮影装置１０を用いたトモシンセシス撮影においては、上述のように、ＦＰＤ１８（Ｘ線検出器１８）とＸ線源１２の移動軸のアライメントに高い精度が要求されている。本実施形態のように、天井走行型のＸ線源１２を用いたＸ線撮影装置１０では、Ｘ線源１２とＦＰＤ１８が一体化していないため、Ｘ線源１２とＦＰＤ１８とを高い精度で設置することが難しい。このため、Ｘ線源１２を移動させた場合、図２に示すように、Ｘ線源１２がＦＰＤ１８の長手方向（Ｙ方向）に対して、ずれてｙθ方向に移動し、その状態で得られた撮影画像を用いて再構成されることがある。この場合、高精度な合成画像（トモシンセシス画像）を得ることができない。
　そこで、本実施形態においては、以下のようにしてキャリブレーションを行い、Ｘ線源１２のＦＰＤ１８に対する移動方向のずれを補正することができる。
　以下、本実施形態のＸ線撮影装置１０のキャリブレーション方法について説明する。
　まず、被写体Ｍのトモシンセシス撮影を行う前、例えば、装置起動時に、Ｘ線源１２を、移動部１４により、移動方向（ｙ方向）に移動させつつ、少なくとも２つの画像を、それぞれマーカ２０を含むようにして撮影する。この場合、例えば、図３に示すように、第１の画像４０と第２の画像４２が得られる。
　この場合、Ｘ線源１２の移動方向（ｙ方向）と、ＦＰＤ１８の長さ方向（Ｙ方向）とが平行であれば、マーカ２０の位置は、ＦＰＤ１８の長さ方向（Ｙ方向）と直交するＸ方向における位置が第１の画像４０と第２の画像４２とで同じである。一方、Ｘ線源１２の移動方向（ｙ方向）が、ＦＰＤ１８の長さ方向に対してずれていれば、Ｘ方向における位置が第１の画像４０と第２の画像４２とで、例えば、δだけずれる。すなわち、δだけシフトする。
　このシフト量δと、第１の画像４０と第２の画像４２の撮影位置間の距離ｄとにより、Ｘ線源１２のずれ角θを求めることができる。すなわち、ずれ角θは、θ＝ｔａｎ－１（δ／ｄ）により算出することができる。
　なお、図４に示すように、第１の画像４０と第２の画像４２とを、撮影画像のＸ方向の端を合わせて合成して合成画像４４を得て、シフト量δを求め、このシフト量δと、第１の画像４０と第２の画像４２の撮影位置間の距離ｄとにより、Ｘ線源１２のずれ角θを求めることもできる。このように、マーカ２０の撮影画像の位置の相対関係を用いてＸ線源１２のずれ角θ（第１の傾き）を求めることができる。なお、符号４４ａは、合成画像４４の外縁を示す。
　ここで、ずれ角θ（第１の傾き）とは、上述の傾きのことから、複数の画素の並び（２次元マトリックスの配列方向）に対するものである。
　Ｘ線源１２のシフト量δ、ずれ角θを用いて、トモシンセシス画像を得る際に、各撮影画像の撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））を調整した後、再構成する。これにより、高精度なトモシンセシス画像を得ることができる。なお、Ｘ線源１２のシフト量δ、ずれ角θを用いて各撮影画像の位置を調整した後ではなく、Ｘ線源１２のシフト量δ、ずれ角θを用いて再構成してもよい。
　また、上述のように、撮影画像におけるマーカ２０の大きさの変動率とＸ線源－ＦＰＤ間の距離Ｌとの関係が予め求められており、例えば、記憶部２２ａに記憶されている。これにより、図５（ａ）に示すマーカ２０が図５（ｂ）に示すマーカ２１のように大きさが変動した場合でも、Ｘ線源－ＦＰＤ間の距離Ｌが得られる。この場合、Ｘ線源１２のＦＰＤ１８の表面に対して直交する面におけるずれ量およびずれ角、すなわち、図１に示すｚ方向におけるずれ量およびずれ角を求めることができる。このように、マーカ２０の撮影画像の位置の相対関係を用いてＸ線源１２のｚ方向におけるずれ角（第２の傾き）を求めることができる。
　Ｘ線源１２のＦＰＤ１８の表面に対するずれ量およびずれ角を用いて、トモシンセシス画像を得る際に、各画像の撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））での大きさを調整した後、再構成する。これにより、高精度なトモシンセシス画像を得ることができる。なお、ＦＰＤ１８の表面に対するずれ量およびずれ角を用いて各画像での撮影物の大きさを調整した後ではなく、ＦＰＤ１８の表面に対するずれ量およびずれ角を用いて再構成してもよい。
　本実施形態においては、以上のようにして、キャリブレーションを行うことができる。
　なお、上記Ｘ線源１２のシフト量δおよびずれ角θ、ならびにｚ方向におけるずれ量およびずれ角の算出、すなわち、キャリブレーションは、撮影台１６の表面１６ａ上に被写体Ｍ（患者）がいない状態で行うことが好ましいが、被写体Ｍ（患者）の撮影状態で行ってもよい。
　また、Ｘ線源１２の移動軸が真直でない場合、すなわち、レール１３がＦＰＤ１８の長手方向（Ｙ方向）と平行でない場合の補正を可能にするため、上記Ｘ線源１２のシフト量δ、ずれ角θならびにｚ方向におけるずれ量およびずれ角は、Ｘ線源１２の位置毎に求めておき、記憶部２２ａに、例えば、テーブルの形式で記憶しておいて、トモシンセシス撮影の際に利用してもよい。
　次に、Ｘ線撮影装置１０のトモシンセシス撮影について説明する。
　まず、被写体Ｍが撮影台１６の表面１６ａに位置決めされた後、入力部２８から撮影開始の指示が与えられると、制御部３０の制御によりトモシンセシス撮影が開始される。
　撮影が開始されると、移動部１４により、Ｘ線源１２をｙ方向に移動しつつ、所定の撮影位置で、被写体Ｍの方向にＸ線源１２の照射角度を変えて、異なる照射角度でＸ線が被写体Ｍに照射され、１回の撮影操作で撮影角度の異なる複数枚の投影データ（Ｘ線画像）が得られる。そして、被写体ＭのＸ線画像の撮影が行われる度に、ＦＰＤ１８から、撮影されたＸ線画像に対応する投影データが画像処理部２２の記憶部２２ａに出力され、ＦＰＤ１８で取得された複数枚のＸ線画像の投影データが記憶部３２に記憶される。
　次に、予め補正部２２ｂには、各撮影位置毎に、位置ずれがあれば、そのずれ量が算出されている。補正部２２ｂにより、位置ずれがあるものについては、投影データの撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））が補正される。補正された撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））が再構成部２２ｃに出力される。
　次に、再構成部２２ｃにおいて、投影データおよび補正された撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））を用いて、被写体Ｍの任意高さの断面におけるＸ線断層画像が再構成される。このようにして、撮影位置の位置ずれを補正しているため、被写体Ｍの所定の位置の断面について、高精度な合成画像（トモシンセシス画像）を得ることができる。
　例えば、得られたトモシンセシス画像（Ｘ線断層画像）は、例えば、表示部２４に表示される。また、トモシンセシス画像（Ｘ線断層画像）は出力部２６に出力され、出力部２６において、例えば、トモシンセシス画像（Ｘ線断層画像）がプリント出力され、Ｘ線断層画像のデジタル画像データが記録媒体に保存される。
　次に、トモシンセシス撮影のＸ線断層画像の再構成について説明する。
　図６（ａ）および（ｂ）は、トモシンセシス撮影によるＸ線断層画像の再構成時の様子を示す概念図である。
　トモシンセシス撮影時に、図６（ａ）に示すようにＸ線源１２が位置Ｓ１からスタートしてＳ３まで移動し、Ｓ１、Ｓ２、Ｓ３の各線源位置において被写体Ｍに放射線が照射され、それぞれ、被写体ＭのＸ線画像Ｐ１、Ｐ２、Ｐ３が得られるものとする。
　ここで、図６（ａ）に示すように、被写体Ｍの高さの異なる２つの位置に撮影対象物Ａ、Ｂが存在するとする。各撮影位置（撮影時のＸ線源１２の位置）Ｓ１、Ｓ２、Ｓ３において、Ｘ線源１２から照射されたＸ線は、被写体Ｍを透過してＦＰＤ１８に入射される。その結果、各撮影位置Ｓ１、Ｓ２、Ｓ３に対応するＸ線画像Ｐ１、Ｐ２、Ｐ３において、２つの撮影対象物Ａ、Ｂは、それぞれ異なる位置関係で投影される。
　例えば、Ｘ線画像Ｐ１の場合、Ｘ線源１２の位置Ｓ１が、Ｘ線源１２の移動方向に対して、撮影対象物Ａ、Ｂよりも左側に位置するため、撮影対象物Ａ、Ｂは、それぞれ、撮影対象物Ａ、Ｂよりも右側にずれたＰ１Ａ、Ｐ１Ｂの位置に投影される。同様に、Ｘ線画像Ｐ２の場合には、ほぼ直下のＰ２Ａ、Ｐ２Ｂの位置に、Ｘ線画像Ｐ３の場合には、左側にずれたＰ３Ａ、Ｐ３Ｂの位置に投影される。
　撮影対象物Ａが存在する高さの断面におけるＸ線断層画像を再構成する場合、Ｘ線源１２の位置に基づいて、撮影対象物Ａの投影位置Ｐ１Ａ、Ｐ２Ａ、Ｐ３Ａが一致するように、例えば、図６（ｂ）に示すように、Ｘ線画像Ｐ１を左へ、Ｘ線画像Ｐ３を右にシフトさせて合成する。これにより、撮影対象物Ａが存在する高さのＸ線断層画像が再構成される。同様にして、任意高さの断面におけるＸ線断層画像も再構成することができる。
　なお、Ｘ線断層画像の再構成方法としては、代表的にシフト加算法が利用される。シフト加算法は、それぞれのＸ線画像の撮影時におけるＸ線源１２の位置に基づいて、撮影した複数枚の投影データ（Ｘ線画像）の位置を順次シフトして加算するものである。
　次に、本発明の第２の実施形態について説明する。
　図７は、本発明の第２の実施形態に係るＸ線撮影装置を示す模式図であり、図８は、本発明の第２の実施形態に係るＸ線撮影装置におけるキャリブレーション方法を示すフローチャートである。
　なお、本実施形態においては、図１～図５に示す第１の実施形態の構成物と同一構成物には同一符号を付して、その詳細な説明は省略する。
　図７に示す第２の実施形態のＸ線撮影装置１０ａは、図１に示す第１の実施形態のＸ線撮影装置１０に比して、判定部２３を有する点、ならびに画像処理部２２の動作およびキャリブレーション方法が異なる点以外は、第１の実施形態のＸ線撮影装置１０と同様の構成および動作をするものであるため、その詳細な説明は省略する。
　本実施形態のＸ線撮影装置１０ａにおいて、判定部２３は制御部３０に接続されている。
　制御部３０は、入力部２８から入力される指示入力信号に従って、Ｘ線源１２、移動部１４およびＸ線源制御部１５によるＸ線画像撮影、画像処理部２２における画像処理、判定部２３におけるキャリブレーションの更新の要否などの判断、表示部２４における画面表示、および出力部２６における出力処理等のＸ線撮影装置１０ａの動作を制御するものである。
　また、判定部２３の判定結果に基づいて、制御部３０を介して画像処理部２２（後述する補正部２２ｂ）に後述する第１のキャリブレーションデータおよび第２のキャリブレーションデータのうち、少なくとも第１のキャリブレーションデータが作成させ、作成されたキャリブレーションデータが画像処理部２２（後述する記憶部２２ａ）に記憶される。
　本実施形態では、記憶部２２ａ、補正部２２ｂおよび制御部３０により、キャリブレーション部が構成される。
　判定部２３は、キャリブレーション部におけるキャリブレーションデータの更新の要否を判定するものである。判定部２３においては、後述するように、キャリブレーション部（補正部２２ｂ）で得られた撮影画像におけるマーカ２０（マーカ像）の位置、および撮影画像におけるマーカ２０（マーカ像）の大きさのうち、少なくとも撮影画像におけるマーカ２０（マーカ像）の位置に基づいて、キャリブレーションデータの更新の要否を判定する。すなわち、判定部２３は、キャリブレーション部にキャリブレーションを実施させて、既に記憶部２２ａに記憶されているキャリブレーションデータを更新するか否かを判定するものである。
　また、判定部２３は、現在保存している第１のキャリブレーションデータを用いて、第２の画像４２におけるマーカ２０の像に対応する対応マーカの像の位置を求めるものでもある。さらには、現在保存している第２のキャリブレーションデータを用いて、第２の画像４２におけるマーカ２０の像に対応する対応マーカの像の大きさを求めるものでもある。
　本実施形態においても、第１の実施形態と同じく撮影台１６には、ＦＰＤ１８が設けられている範囲内にマーカ２０が、例えば、１個設けられている。このマーカ２０は、ＦＰＤ１８の外縁１８ｂにエッジを一致させて設けられている。また、マーカ２０は、ＦＰＤ１８の表面１８ａに対して隙間をあけて撮影台１６に設けられている。
　また、マーカ２０は、撮影台１６に対して着脱自在であってもよい。このように、マーカ２０を着脱自在にすることにより、キャリブレーションの更新を容易にすることができる。
　次に、画像処理部２２の各部について詳細に説明する。
　記憶部２２ａは、ＦＰＤ１８で取得されたＸ線画像の複数枚の投影データが入力されて、これらのＸ線画像の投影データを記憶するものである。
　また、記憶部２２ａには、投影データにおけるマーカ２０の設計上の大きさ、ならびに後述するマーカ２０の像の位置、マーカ２０の像の大きさの一致度に関する閾値等も記憶される。
　さらには、記憶部２２ａには、後述する第１のキャリブレーションデータ、第２のキャリブレーションデータが記憶されるものでもある。この第１のキャリブレーションデータ、第２のキャリブレーションデータは、上述のように判定部２３の判定結果に応じて適宜更新される。
　補正部２２ｂは、Ｘ線源１２の各撮影位置における、後に詳細に説明する第１のキャリブレーションデータを算出するものである。また、さらには、後述するように、第２のキャリブレーションデータを算出するものである。
　補正部２２ｂは、撮影画像におけるマーカ２０の像位置を検出する位置検出機能を有しており、互いに異なるＸ線管球位置から撮影された第１の画像および第２の画像から、それぞれマーカ像の位置を自動検出し、第１のキャリブレーションデータの生成に用いるとともに、判定部２３における第１のキャリブレーションデータ更新要否の判定に供する。
　判定部２３は、第１の画像から検出されたマーカ像位置と第１のキャリブレーションデータから第２の画像のマーカ像のあるべき対応マーカ像位置を求め、これを第２の画像から検出されたマーカ像位置と比較し、ずれ量を求める。このずれ量を予め設定された閾値と比較することにより、キャリブレーションの要否が判定される。
　補正部２２ｂにおいては、上述の判定部２３のキャリブレーションデータを更新するか否かの判定結果に応じて、第１のキャリブレーションデータを算出し、キャリブレーションデータを更新するものである。
　判定部２３で、第１のキャリブレーションデータに関し、キャリブレーションが必要であると判定された場合、補正部２２ｂにおいては、移動部１４により移動方向（ｙ方向）にＸ線源１２を移動させつつ、撮影された少なくとも２つの画像を用いて第１のキャリブレーションデータを作成する。例えば、第１の実施形態と同様に図３に示すように、第１の画像４０と第２の画像４２を用いて第１のキャリブレーションデータを作成する。
　この場合、Ｘ線源１２の移動方向（ｙ方向）と、ＦＰＤ１８の長さ方向（Ｙ方向）とが平行であれば、マーカ２０の位置は、ＦＰＤ１８の長さ方向（Ｙ方向）と直交するＸ方向における位置が第１の画像４０と第２の画像４２とで同じである。すなわち、シフト量がゼロである。
　一方、Ｘ線源１２の移動方向（ｙ方向）がＦＰＤ１８の長さ方向に対して、例えば、第１の実施形態と同様に図２に示すようにｙθ方向にＸ線源１２がずれて移動するものであれば、Ｘ方向における位置が第１の画像４０と第２の画像４２とでマーカ２０の像の位置が、例えば、δだけずれる。すなわち、シフト量δである。
　このシフト量δと、第１の画像４０と第２の画像４２の撮影位置間の距離ｄとにより、Ｘ線源１２のずれ角θ（第１の傾き）を求めることができる。すなわち、ずれ角θは、θ＝ｔａｎ－１（δ／ｄ）により算出することができる。
　補正部２２ｂは、ずれ角θ（第１の傾き）に基づいて、第１のキャリブレーションデータを作成し、記憶部２２ａに記憶させる。
　第１のキャリブレーションデータは、例えば、ずれ角θに応じた分だけ、設計値（複数の画像をそれぞれどの位置で撮るかの目標値）をベースにしてＸ線源１２の移動軸が傾いている分を補正することによって、真のＸ線源１２の位置に近い情報に補正するためのＸ線源１２の位置情報（Ｘ線管球位置情報）である。この第１のキャリブレーションデータにより、図２に示すｙθ方向のずれを補正することができる。
　なお、第１の実施形態と同様に図４に示すように、第１の画像４０と第２の画像４２とを、撮影画像のＸ方向の端を合わせて合成して合成画像４４を得て、合成画像４４の外縁４４ａを基準として、マーカ２０のシフト量δを求め、このシフト量δと、第１の画像４０と第２の画像４２の撮影位置間の距離ｄとにより、Ｘ線源１２のずれ角θ（第１の傾き）を求めることもできる。このように、マーカ２０の撮影画像の位置の相対関係を用いてＸ線源１２のずれ角θ（第１の傾き）を求めることができる。このずれ角θ（第１の傾き）に基づいて、上述の第１のキャリブレーションデータを作成することもできる。
　なお、ずれ角θ（第１の傾き）については、第１の実施形態で説明した通りであるため、その詳細な説明は省略する。
　また、Ｘ線源１２とＦＰＤ１８の表面１８ａとの距離（Ｘ線源－ＦＰＤ間の距離Ｌ）のずれ量に対しては、Ｘ線源－ＦＰＤ間の距離Ｌが設定された距離から変化すると、撮影画像におけるマーカ２０の大きさが変わる。このため、本実施形態においては、Ｘ線源－ＦＰＤ間の距離Ｌのずれ量に対して、撮影画像におけるマーカ２０の大きさ、撮影画像におけるマーカ２０の変化率、すなわち、各距離Ｌにおけるマーカ２０の大きさ、マーカ２０の拡大率または縮小率を求めておき、これを、例えば、記憶部２２ａに記憶させておく。さらには、撮影画像におけるマーカ２０の大きさ、マーカ２０の拡大率または縮小率に対して、閾値を設定しておき、この閾値を記憶部２２ａに記憶させておく。
　補正部２２ｂは、マーカ２０の像の位置以外に、互いに異なるＸ線管球位置から撮影された第１の画像および第２の画像から、それぞれマーカ像の大きさを自動検出し、第２のキャリブレーションデータの生成に用いるとともに、判定部２３における第２のキャリブレーションデータ更新要否の判定に供することもできる。
　この場合、判定部２３は、第１の画像４０のマーカ２０の像の大きさおよび第２の画像４２のマーカ２０の像の大きさを求める。第２のキャリブレーションデータから第２の画像４２のマーカ像のあるべき対応マーカ像の大きさを求め、第２の画像４２のマーカ２０の大きさと比較し、ずれ量を求める。このずれ量を予め設定された閾値と比較することにより、キャリブレーションの要否が判定される。なお、ずれ量は、マーカ２０の像の直径の差、マーカ２０の像の拡大率または縮小率等を用いることができる。
　判定部２３で、第２のキャリブレーションデータに関し、キャリブレーションが必要であると判定された場合、補正部２２ｂは、移動部１４により移動方向（ｙ方向）にＸ線源１２を移動させつつ、撮影された少なくとも２つの画像を用いて第２のキャリブレーションデータを作成する。
　例えば、２つの撮影位置で撮影した結果、第１の実施形態と同様に図５（ａ）に示すマーカ２０が図５（ｂ）に示すマーカ２１のように大きさが変動した場合、補正部２２ｂにおいては、各撮影位置でのＸ線源－ＦＰＤ間の距離Ｌを得ることができる。更には、補正部２２ｂでは、各撮影位置の距離Ｌと撮影位置間の距離ｄとにより、図７に示すＸ線源１２のＦＰＤ１８の表面に対して直交する面におけるずれ角（第２の傾き）を求めることができる。このように、マーカ２０の撮影画像の位置の相対関係を用いてＸ線源１２のｚ方向におけるずれ角（第２の傾き）を求めることができる。
　補正部２２ｂは、このずれ角（第２の傾き）に基づいて、第２のキャリブレーションデータを作成し、記憶部２２ａに記憶させる。
　第２のキャリブレーションデータとしては、例えば、拡大率または縮小率に応じた分だけ、マーカ２０の大きさが設計上の大きさとなるように投影データ全体を拡大または縮小する拡大率データまたは縮小率データ等である。この第２のキャリブレーションデータにより、図７に示すｚ方向におけるずれを補正することができる。
　なお、判定部２３において、マーカ２０の像の位置の一致度の閾値、マーカ２０の像の大きさの一致度の閾値は、例えば、入力部２８を介して設定できるようにしてもよい。更には、補正部２２ｂにより作成されるキャリブレーションデータ（第１のキャリブレーションデータおよび第２のキャリブレーションデータ）についても、例えば、入力部２８を介して設定できるようにしてもよい。
　本実施形態においては、第１のキャリブレーションデータおよび第２のキャリブレーションデータについて、Ｘ線源１２の位置毎に求めておき、記憶部２２ａに、例えば、テーブルの形式で記憶しておいて、トモシンセシス撮影の際に利用してもよい。
　再構成部２２ｃは、必要に応じてなされた補正部２２ｂによる補正処理後の投影データを含む、複数枚のＸ線画像の投影データを、必要に応じて補正されたＸ線源１２の位置情報（Ｘ線管球位置情報）を用いて画像合成処理を施し、被写体Ｍの任意高さの断面におけるＸ線断層画像を再構成するものである。
　なお、補正部２２ｂにおいて、第１のキャリブレーションデータおよび第２のキャリブレーションデータを作成することが望ましいが、少なくとも第１のキャリブレーションデータを作成すればよい。第１のキャリブレーションデータおよび第２のキャリブレーションデータのうち、少なくとも第１のキャリブレーションデータを用いることにより、高精度なトモシンセシス画像を得ることができ、第１のキャリブレーションデータおよび第２のキャリブレーションデータの両方を用いることにより、更に高精度なトモシンセシス画像を得ることができる。このため、再構成部２２ｃで、被写体Ｍの任意高さの断面におけるＸ線断層画像を再構成する際には、補正部２２ｂにより、少なくとも第１のキャリブレーションデータを用いて補正されたＸ線源１２の位置情報（Ｘ線管球位置情報）を用いればよい。
　本実施形態においても、画像処理部２２（記憶部２２ａ、補正部２２ｂおよび再構成部２２ｃ）は、ハードウェア（装置）で構成することに限定されるものではない。例えば、記憶部２２ａに記憶された投影データに対して、画像処理部２２の補正部２２ｂおよび再構成部２２ｃによる処理を、コンピュータに実行させるためのプログラムとすることもできる。
　表示部２４は、第１の実施形態と同様の構成であるが、更に判定部２３によるキャリブレーションの更新の要否の判定結果を通知する通知部として機能するものである。判定部２３によるキャリブレーションの更新の要否の判定結果が、制御部３０により、例えば、文字等の形態で表示部２４に表示される。
　本実施形態のＸ線撮影装置１０ａを用いたトモシンセシス撮影においては、上述のように、ＦＰＤ１８（Ｘ線検出器１８）とＸ線源１２の移動軸のアライメントに高い精度が要求されている。本実施形態のように、天井走行型のＸ線源１２を用いたＸ線撮影装置１０ａでは、Ｘ線源１２とＦＰＤ１８が一体化していないため、Ｘ線源１２とＦＰＤ１８とを高い精度で設置することが難しい。このため、Ｘ線源１２を移動させた場合、図２に示すように、Ｘ線源１２がＦＰＤ１８の長手方向（Ｙ方向）に対してずれて、Ｘ線源１２がｙθ方向に移動し、その状態で得られた撮影画像を用いて再構成されることがある。この場合、高精度な合成画像（トモシンセシス画像）を得ることができない。また、キャリブレーションしても経時劣化（経時変化）により、そのキャリブレーションデータが適正ではないこともある。
　しかしながら、本実施形態においては、図８に示すようにしてキャリブレーションの要否を判定し、キャリブレーションデータの経時劣化（経時変化）の影響を抑制している。
　なお、Ｘ線撮影装置１０ａのキャリブレーションは、撮影台１６の表面１６ａ上に被写体Ｍ（患者）がいない状態で行うことが好ましいが、被写体Ｍ（患者）の撮影状態で行ってもよい。
　Ｘ線撮影装置１０ａについて、既にキャリブレーションがなされており、第１のキャリブレーションデータおよび第２のキャリブレーションデータが、初期キャリブレーションデータとして記憶部２２ａに記憶されている（ステップＳ１０）。
　次に、被写体Ｍのトモシンセシス撮影を行う前、例えば、装置起動時に、移動部１４によりＸ線源１２を少なくとも２つの撮影位置に移動させ、各撮影位置において、Ｘ線源１２からマーカ２０が撮影されるようにＸ線をＦＰＤ１８に照射して、ＦＰＤ１８からマーカ２０を含む少なくとも２つの撮影画像（例えば、第１の画像４０、第２の画像４２）の各画像データが画像処理部２２の記憶部２２ａに出力される。
　次に、補正部２２ｂにおいて、例えば、第１の画像４０のマーカ２０の像の位置および第２の画像４２のマーカ２０の像の位置を求める。さらに、判定部２３において現在保存している第１のキャリブレーションデータを用いて、第２の画像４２におけるマーカ２０の像に対応する対応マーカの像の位置が求められる。そして、補正部２２ｂにおいて、第２の画像４２のマーカ２０の位置と、第１のキャリブレーションデータを用いて得られた対応マーカの像の位置とのずれ量が求められる（ステップＳ１２）。
　次に、マーカ２０の像の位置におけるずれ量が、判定部２３に出力される。
　次に、判定部２３において、記憶部２２ａに記憶されたマーカ２０の像の位置におけるずれ量の閾値に基づいて、キャリブレーション（第１のキャリブレーションデータの作成）の要否が判定される（ステップＳ１４）。なお、第１のキャリブレーションデータに関するキャリブレーションの要否の判定結果は、例えば、表示部２４に文字等の形態で表示されて、使用者等に通知される。
　ステップＳ１４において、マーカ２０の像の位置におけるずれ量が閾値以内であれば、初期キャリブレーションデータを更新することなく、記憶部２２ａに記憶されている第１のキャリブレーションデータが、Ｘ線断層画像の再構成に利用される。
　一方、ステップＳ１４において、マーカ２０の像の位置におけるずれ量が閾値を越えた場合、再度、第１のキャリブレーションデータを作成する（ステップＳ１８）。
　なお、ステップＳ１２において、例えば、第１の画像４０のマーカ２０の像の大きさおよび第２の画像４２のマーカ２０の像の大きさを求める。さらに、判定部２３において現在保存している第２のキャリブレーションデータを用いて、第２の画像４２におけるマーカ２０の像に対応する対応マーカの像の大きさが求められる。そして、補正部２２ｂにおいて、第２の画像４２のマーカ２０の大きさと、第２のキャリブレーションデータを用いて得られた対応マーカの像の大きさとのずれ量を求めてもよい（ステップＳ１２）。
　この場合、次に、マーカ２０の像の大きさのずれ量が、判定部２３に出力される。判定部２３において、記憶部２２ａに記憶されたマーカ２０の像の大きさのずれ量の閾値に基づいて、キャリブレーション（第２のキャリブレーションデータの作成）の要否が判定される（ステップＳ１４）。なお、第２のキャリブレーションデータに関するキャリブレーションの要否の判定結果は、例えば、表示部２４に文字等の形態で表示されて、使用者等に通知される。
　ステップＳ１４において、マーカ２０の像の大きさのずれ量が閾値以内であれば、初期キャリブレーションデータを更新することなく、記憶部２２ａに記憶されている第２のキャリブレーションデータが、Ｘ線断層画像の再構成に利用される。
　一方、ステップＳ１４において、マーカ２０の像の大きさのずれ量が閾値を越えた場合、再度、第２のキャリブレーションデータを作成する（ステップＳ１８）。
　なお、第１のキャリブレーションデータおよび第２のキャリブレーションデータの作成方法は、既に説明したように、少なくとも２つの撮影位置からマーカ２０を含む画像を撮影し、得られた各撮影画像データを用いる方法であり、その作成方法についての詳細な説明は省略する。
　次に、作成された第１のキャリブレーションデータおよび第２のキャリブレーションデータと、初期キャリブレーションデータとして記憶部２２ａに記憶されていた第１のキャリブレーションデータおよび第２のキャリブレーションデータとを書き換え、記憶部２２ａに記憶させる。すなわち、初期キャリブレーションデータを更新する（ステップＳ２０）。
　このようにして、初期キャリブレーションデータの更新の要否が判定され、判定結果に基づいて初期キャリブレーションデータが更新される。このため、Ｘ線源１２の移動方向が経時変化によりずれて、初期キャリブレーションデータが経時劣化（経時変化）した場合でも、適正なキャリブレーションデータとすることができる。
　本実施形態においては、例えば、Ｘ線源１２の移動軸が真直ではない場合、すなわち、レール１３がＦＰＤ１８の長手方向（Ｙ方向）と平行ではない場合（Ｘ線源１２の移動方向（ｙ方向）とＦＰＤ１８の画素配列方向とが平行ではない場合）等、Ｘ線源１２の移動軸がずれていても、そのずれを補正できるとともに、Ｘ線源１２の移動軸のずれの経時変化の影響を抑制しつつ、常に適正なキャリブレーションデータを用いて、トモシンセシスの再構成を、常に高い精度で行うことができ、被写体の所定の位置の断面について、高精度な合成画像（トモシンセシス画像）を常に得ることができる。
　次に、Ｘ線撮影装置１０ａのトモシンセシス撮影について説明する。
　まず、被写体Ｍが撮影台１６の表面１６ａに位置決めされた後、入力部２８から撮影開始の指示が与えられると、制御部３０の制御によりトモシンセシス撮影が開始される。
　撮影が開始されると、移動部１４により、Ｘ線源１２をｙ方向に移動しつつ、所定の撮影位置で、被写体Ｍの方向にＸ線源１２の照射角度を変えて、異なる照射角度でＸ線が被写体Ｍに照射され、１回の撮影操作で撮影角度の異なる複数枚の投影データ（Ｘ線画像）が得られる。そして、被写体ＭのＸ線画像の撮影が行われる度に、ＦＰＤ１８から、撮影されたＸ線画像に対応する投影データが画像処理部２２の記憶部２２ａに出力され、ＦＰＤ１８で取得された複数枚のＸ線画像の投影データが記憶部３２に記憶される。
　次に、予め記憶部２２ａには、例えば、各撮影位置毎に、上述の第１のキャリブレーションデータおよび第２のキャリブレーションデータが記憶されている。また、図８に示す方法のように撮影前に、キャリブレーションデータの更新の要否が判定され、判定結果に応じて、キャリブレーションデータは更新されている。なお、このキャリブレーションの要否の判定結果は、例えば、表示部２４に文字等の形態で表示される。
　次に、補正部２２ｂにより、第１のキャリブレーションデータおよび第２のキャリブレーションデータに基づいて位置ずれがあるものについては、投影データの撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））が補正される。補正された撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））が再構成部２２ｃに出力される。
　次に、再構成部２２ｃにおいて、投影データおよび補正された撮影位置（Ｘ線源１２の位置情報（Ｘ線管球位置情報））を用いて、被写体Ｍの任意高さの断面におけるＸ線断層画像が再構成される。このようにして、撮影位置の位置ずれを補正しているため、被写体Ｍの所定の位置の断面について、高精度な合成画像（トモシンセシス画像）を得ることができる。
　例えば、得られたトモシンセシス画像（Ｘ線断層画像）は、例えば、表示部２４に表示される。また、トモシンセシス画像（Ｘ線断層画像）は出力部２６に出力され、出力部２６において、例えば、トモシンセシス画像（Ｘ線断層画像）がプリント出力され、Ｘ線断層画像のデジタル画像データが記録媒体に保存される。
　なお、本実施形態のトモシンセシス撮影のＸ線断層画像の再構成は、上述の第１の実施形態の再構成と同じ方法でなされるものであるため、その詳細な説明は省略する。
　本発明は、基本的に以上のように構成されるものである。以上、本発明のＸ線撮影装置およびそのキャリブレーション方法について詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良または変更をしてもよいのはもちろんである。
　１０、１０ａ　Ｘ線撮影装置
　１２　Ｘ線源
　１３　レール
　１４　移動部
　１５　Ｘ線源制御部
　１６　撮影台
　１８　Ｘ線検出器（ＦＰＤ）
　２０　マーカ
　２２　画像処理部
　２３　判定部
　２４　表示部
　２６　出力部
　２８　入力部
　３０　制御部
　Ｍ　被写体

Claims (21)

  1. 　所定の移動経路で移動可能なＸ線源と、
    　前記Ｘ線源を所定の移動経路で移動させる移動手段と、
    　Ｘ線源に対向して設けられた撮影台と、
    　前記撮影台に設けられた平板状のＸ線検出器と、
    　前記撮影台に設けられるマーカと、
    　前記移動手段により前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させる制御部と、
    　前記撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、前記Ｘ線源の移動軸の、前記Ｘ線検出器に対する傾きを求める画像処理部とを有することを特徴とするＸ線撮影装置。
  2. 　前記Ｘ線検出器に対する傾きは、前記Ｘ線検出器の表面と平行な面内での第１の傾きである請求項１に記載のＸ線撮影装置。
  3. 　前記画像処理部は、さらに前記撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いて前記Ｘ線源と前記Ｘ線検出器との距離を求めて、前記Ｘ線検出器に対する傾きとして、前記Ｘ線検出器の表面と直交する面内での第２の傾きを求める請求項１または２に記載のＸ線撮影装置。
  4. 　前記画像処理部は、前記Ｘ線源の撮影位置毎に、前記Ｘ線検出器に対する前記傾きを求める請求項１～３のいずれか１項に記載のＸ線撮影装置。
  5. 　前記画像処理部は、前記移動手段により前記Ｘ線源が所定の撮影位置に移動されトモシンセシス撮影されて得られた複数の画像を用いて断層画像を再構成する機能を有し、
    　前記画像処理部は、前記Ｘ線検出器に対する前記傾きを、前記断層画像の再構成に用いる請求項１～４のいずれか１項に記載のＸ線撮影装置。
  6. 　Ｘ線撮影装置のキャリブレーション方法であって、
    　前記Ｘ線撮影装置は、所定の移動経路で移動可能なＸ線源に対向して設けられた撮影台と、前記撮影台に設けられた平板状のＸ線検出器と、マーカとを備えており、
    　前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させる工程と、
    　前記撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、前記Ｘ線源の移動軸の、前記Ｘ線検出器に対する傾きを求める工程を有することを特徴とするキャリブレーション方法。
  7. 　前記Ｘ線検出器に対する前記傾きは、前記Ｘ線検出器の表面と平行な面内での第１の傾きである請求項６に記載のキャリブレーション方法。
  8. 　さらに前記撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いて前記Ｘ線源と前記Ｘ線検出器との距離を求めて、前記Ｘ線検出器に対する傾きとして、前記Ｘ線検出器の表面と直交する面内での第２の傾きを求める工程を有する請求項６または７に記載のキャリブレーション方法。
  9. 　前記Ｘ線源の撮影位置毎に、前記Ｘ線検出器に対する前記傾きを求める請求項６～８のいずれか１項に記載のキャリブレーション方法。
  10. 　所定の移動経路で移動可能なＸ線源と、
    　前記Ｘ線源を所定の移動経路で移動させる移動手段と、
    　Ｘ線源に対向して設けられた撮影台と、
    　前記撮影台に設けられた平板状のＸ線検出器と、
    　前記撮影台に設けられるマーカと、
    　前記Ｘ線検出器に対して配置された前記Ｘ線源の移動方向のずれに関するキャリブレーションを行うキャリブレーション部と、
    　前記キャリブレーション部によるキャリブレーションの更新の要否を判定し、判定結果に基づいて前記キャリブレーション部にキャリブレーションを行わせる判定部とを有し、
    　前記Ｘ線源により前記マーカを含む画像が撮影され、前記キャリブレーション部により前記撮影された画像のマーカ像の位置が求められ、前記判定部は、このマーカ像の位置に基づいて前記キャリブレーションの更新の要否を判定することを特徴とするＸ線撮影装置。
  11. 　前記キャリブレーション部により、前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させ、前記撮影された第１の画像、第２の画像についてマーカ像の位置が求められ、
    　前記判定部は、既存のキャリブレーションデータを用いて前記第１の画像のマーカ像の位置から前記第２の画像のマーカ像に対応する対応マーカ像の位置を求め、前記第２の画像のマーカ像の位置と、前記第２の画像のマーカ像に対応する対応マーカ像の位置とのずれに基づいて前記キャリブレーションの更新の要否を判定する請求項１０に記載のＸ線撮影装置。
  12. 　前記キャリブレーション部は、前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させ、前記撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、前記Ｘ線源の移動軸の、前記Ｘ線検出器の表面と平行な面内での第１の傾きを求め、前記第１の傾きに基づいて第１のキャリブレーションデータを作成するものである請求項１０に記載のＸ線撮影装置。
  13. 　前記キャリブレーション部は、前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させ、前記撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いて前記Ｘ線源と前記Ｘ線検出器との距離を求めて、前記Ｘ線検出器の表面と直交する面内での前記Ｘ線源の第２の傾きを求め、前記第２の傾きに基づいて第２のキャリブレーションデータを作成するものである請求項１１または１２に記載のＸ線撮影装置。
  14. 　前記キャリブレーション部は、前記Ｘ線源の撮影位置毎に、前記第１のキャリブレーションデータおよび前記第２のキャリブレーションデータのうち、少なくとも前記第１のキャリブレーションデータを求める請求項１３に記載のＸ線撮影装置。
  15. 　さらに、前記移動手段により前記Ｘ線源が所定の撮影位置に移動されトモシンセシス撮影されて得られた複数の画像を用いて断層画像を再構成する再構成部を有し、
    　前記再構成部は、トモシンセシス撮影で撮影された複数の画像を用いて断層画像を再構成する際に、前記第１のキャリブレーションデータおよび前記第２のキャリブレーションデータのうち、少なくとも前記第１のキャリブレーションデータを用いる請求項１３または１４に記載のＸ線撮影装置。
  16. 　さらに、前記判定部による前記キャリブレーションの更新の要否の判定結果を通知する通知部を有する請求項１０～１５のいずれか１項に記載のＸ線撮影装置。
  17. 　Ｘ線撮影装置のキャリブレーション方法であって、
    　前記Ｘ線撮影装置は、所定の移動経路で移動可能なＸ線源に対向して設けられた撮影台と、前記撮影台に設けられた平板状のＸ線検出器と、マーカとを備えており、
    　前記Ｘ線源により前記マーカを含む画像を撮影させる工程と、
    　前記撮影された画像のマーカ像の位置を求め、このマーカ像の位置に基づいて、キャリブレーションの更新の要否を判定し、判定結果に基づいてキャリブレーションを行う工程とを有することを特徴とするキャリブレーション方法。
  18. 　前記キャリブレーション工程は、前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させ、前記撮影された第１の画像、第２の画像についてそれぞれマーカ像の位置を求める工程と、
    　既存のキャリブレーションデータを用いて前記第１の画像のマーカ像の位置から前記第２の画像のマーカ像に対応する対応マーカ像の位置を求める工程と、
    　前記第２の画像のマーカ像の位置と、前記第２の画像のマーカ像に対応する対応マーカ像の位置とのずれに基づいて前記キャリブレーションの更新の要否を判定する工程とを有する請求項１７に記載のキャリブレーション方法。
  19. 　前記キャリブレーション工程は、前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させ、前記撮影された各画像についてマーカ像の位置を求め、このマーカ像の位置の相対関係に基づいて、前記Ｘ線源の移動軸の、前記Ｘ線検出器の表面と平行な面内での第１の傾きを求め、前記第１の傾きに基づいて第１のキャリブレーションデータを作成する工程を有する請求項１７に記載のキャリブレーション方法。
  20. 　前記キャリブレーション工程は、前記Ｘ線源を移動させて、少なくとも２つの位置から、それぞれ前記マーカを含む画像を撮影させ、前記撮影された各画像についてマーカ像の大きさの変動を算出し、この変動を用いて前記Ｘ線源と前記Ｘ線検出器との距離を求めて、前記Ｘ線検出器の表面と直交する面内での前記Ｘ線源の第２の傾きを求め、前記第２の傾きに基づいて第２のキャリブレーションデータを作成するものである請求項１７または１９に記載のキャリブレーション方法。
  21. 　前記Ｘ線源の撮影位置毎に、前記第１のキャリブレーションデータおよび前記第２のキャリブレーションデータのうち、少なくとも前記第１のキャリブレーションデータを求める請求項２０に記載のキャリブレーション方法。
PCT/JP2012/067336 2011-07-06 2012-07-06 Ｘ線撮影装置およびそのキャリブレーション方法 WO2013005833A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP12807484.6A EP2730226A4 (en) 2011-07-06 2012-07-06 X-RAY IMAGING DEVICE AND CORRESPONDING CALIBRATION METHOD
CN201280033444.4A CN103635139A (zh) 2011-07-06 2012-07-06 X射线成像设备和用于其的校准方法
US14/147,113 US9380985B2 (en) 2011-07-06 2014-01-03 X-ray tomosynthesis imaging device and calibration method of an X-ray tomosynthesis imaging device

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2011149912A JP5702236B2 (ja) 2011-07-06 2011-07-06 Ｘ線撮影装置およびそのキャリブレーション方法
JP2011-149912 2011-07-06
JP2011-153815 2011-07-12
JP2011153815A JP5702240B2 (ja) 2011-07-12 2011-07-12 Ｘ線撮影装置およびそのキャリブレーション方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/147,113 Continuation US9380985B2 (en) 2011-07-06 2014-01-03 X-ray tomosynthesis imaging device and calibration method of an X-ray tomosynthesis imaging device

Publications (1)

Publication Number Publication Date
WO2013005833A1 true WO2013005833A1 (ja) 2013-01-10

Family

ID=47437176

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/067336 WO2013005833A1 (ja) 2011-07-06 2012-07-06 Ｘ線撮影装置およびそのキャリブレーション方法

Country Status (4)

Country Link
US (1) US9380985B2 (ja)
EP (1) EP2730226A4 (ja)
CN (1) CN103635139A (ja)
WO (1) WO2013005833A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104414662A (zh) * 2013-09-04 2015-03-18 江苏瑞尔医疗科技有限公司 成像设备的位置校准和误差补偿装置及其补偿方法
WO2015079120A1 (en) * 2013-11-29 2015-06-04 Planmed Oy Mammography imaging arrangement for tomosynthesis
CN107427274A (zh) * 2014-02-21 2017-12-01 三星电子株式会社 断层扫描设备及其用于重构断层扫描图像的方法
CN110098102A (zh) * 2018-01-31 2019-08-06 帕罗德克斯集团公司 减少成像板扫描器中的部件的校准
US10748293B2 (en) 2014-02-21 2020-08-18 Samsung Electronics Co., Ltd. Tomography apparatus and method for reconstructing tomography image thereof

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6105903B2 (ja) * 2012-11-09 2017-03-29 キヤノン株式会社 画像処理装置、画像処理方法、放射線撮影システム及びプログラム
JP6379785B2 (ja) * 2014-07-18 2018-08-29 コニカミノルタ株式会社 断層画像生成システム
KR102393294B1 (ko) * 2014-09-26 2022-05-03 삼성전자주식회사 의료 영상 장치 및 의료 영상 장치의 제어 방법
DE102014223103A1 (de) 2014-11-12 2016-05-12 Siemens Aktiengesellschaft Tomographisches Gerät sowie Verfahren und Computerprogramm zur modellbasierten Positionierung
KR102501712B1 (ko) * 2014-12-12 2023-02-21 삼성전자주식회사 엑스선 장치 및 그 동작 방법
US9872663B2 (en) * 2015-02-04 2018-01-23 Dentsply Sirona Inc. Methods, systems, apparatuses, and computer programs for removing marker artifact contribution from a tomosynthesis dataset
JP6631624B2 (ja) * 2015-04-24 2020-01-15 株式会社ニコン Ｘ線検査装置、ｘ線検査方法および構造物の製造方法
US11071504B2 (en) 2016-06-07 2021-07-27 Koninklijke Philips N.V. Calibration-free tomosynthesis
DE102017205113A1 (de) * 2017-03-27 2018-09-27 Siemens Aktiengesellschaft Ermitteln der Pose einer Röntgeneinheit relativ zu einem Objekt anhand eines digitalen Modells des Objekts
WO2019040056A1 (en) * 2017-08-23 2019-02-28 Carestream Dental Technology Topco Limited TOMOSYNTHESIS SYSTEM ON THE SIDE OF A DENTAL ARMCHAIR
US10531850B2 (en) * 2017-09-07 2020-01-14 General Electric Company Mobile X-ray imaging with detector docking within a spatially registered compartment
US10893842B2 (en) 2018-02-08 2021-01-19 Covidien Lp System and method for pose estimation of an imaging device and for determining the location of a medical device with respect to a target
JP7118750B2 (ja) * 2018-06-04 2022-08-16 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
US10743822B2 (en) * 2018-06-29 2020-08-18 Carestream Health, Inc. Fiducial marker for geometric calibration of bed-side mobile tomosynthesis system
JP6988732B2 (ja) * 2018-08-02 2022-01-05 株式会社島津製作所 放射線画像処理装置および放射線画像処理方法
KR20200050515A (ko) * 2018-11-01 2020-05-12 삼성디스플레이 주식회사 표시 장치 및 그 제어 방법
WO2020129384A1 (ja) * 2018-12-21 2020-06-25 株式会社島津製作所 放射線画像処理装置、放射線画像処理方法及び放射線画像処理プログラム
CN109924994B (zh) * 2019-04-02 2023-03-14 晓智未来(成都)科技有限公司 一种x光拍摄过程中的检测位置自动校准方法及系统
CN111514476B (zh) * 2020-04-30 2022-03-15 江苏瑞尔医疗科技有限公司 一种用于x射线图像引导系统中的校准方法
JP2022131969A (ja) * 2021-02-26 2022-09-07 富士フイルム株式会社 放射線画像撮影システム、画像処理装置、画像処理方法、及び画像処理プログラム
CN113892959B (zh) * 2021-09-26 2024-03-26 有方（合肥）医疗科技有限公司 X射线成像系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06142093A (ja) * 1992-11-09 1994-05-24 Toshiba Corp Ｘ線断層撮影装置
JPH08317287A (ja) * 1995-05-18 1996-11-29 Shimadzu Corp Ｘ線撮像装置
JP2000278606A (ja) 1999-03-24 2000-10-06 Hitachi Medical Corp Ｘ線画像撮影装置
JP2010252951A (ja) 2009-04-23 2010-11-11 Hitachi Medical Corp Ｘ線ｃｔ装置およびこれを用いたデータ取得方法
JP2011067503A (ja) 2009-09-28 2011-04-07 Fujifilm Corp 放射線撮影装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4629313A (en) * 1982-10-22 1986-12-16 Nippon Kogaku K.K. Exposure apparatus
JPH03232215A (ja) * 1990-02-08 1991-10-16 Toshiba Corp 位置合せ方法
US6379043B1 (en) * 1998-12-08 2002-04-30 U.S. Philips Corporation X-ray examination apparatus and method for generating distortion-free X-ray images
FR2812295B1 (fr) 2000-07-27 2003-01-31 Snf Sa Polymeres cationiques de haut poids moleculaire, procede pour leur preparation, et leurs applications
US6671349B1 (en) * 2000-11-13 2003-12-30 Olganix Corporation Tomosynthesis system and registration method
DE60319462T2 (de) * 2002-06-11 2009-03-12 Asml Netherlands B.V. Lithographischer Apparat und Verfahren zur Herstellung eines Artikels
US7656998B2 (en) * 2005-11-14 2010-02-02 Accuray Incorporated Unified quality assurance for a radiation treatment delivery system
US7950849B2 (en) * 2005-11-29 2011-05-31 General Electric Company Method and device for geometry analysis and calibration of volumetric imaging systems
US8000522B2 (en) * 2007-02-02 2011-08-16 General Electric Company Method and system for three-dimensional imaging in a non-calibrated geometry
US8038347B2 (en) 2007-12-21 2011-10-18 General Electric Company Portable tomographic diagnostic system with open gantry
US8007173B2 (en) * 2009-10-14 2011-08-30 Siemens Medical Solutions Usa, Inc. Calibration of imaging geometry parameters

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06142093A (ja) * 1992-11-09 1994-05-24 Toshiba Corp Ｘ線断層撮影装置
JPH08317287A (ja) * 1995-05-18 1996-11-29 Shimadzu Corp Ｘ線撮像装置
JP2000278606A (ja) 1999-03-24 2000-10-06 Hitachi Medical Corp Ｘ線画像撮影装置
JP2010252951A (ja) 2009-04-23 2010-11-11 Hitachi Medical Corp Ｘ線ｃｔ装置およびこれを用いたデータ取得方法
JP2011067503A (ja) 2009-09-28 2011-04-07 Fujifilm Corp 放射線撮影装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2730226A4

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104414662A (zh) * 2013-09-04 2015-03-18 江苏瑞尔医疗科技有限公司 成像设备的位置校准和误差补偿装置及其补偿方法
WO2015079120A1 (en) * 2013-11-29 2015-06-04 Planmed Oy Mammography imaging arrangement for tomosynthesis
CN105960204A (zh) * 2013-11-29 2016-09-21 普兰梅德有限公司 用于断层合成的乳房x射线摄影成像布置
US10881374B2 (en) 2013-11-29 2021-01-05 Planmeca Oy Mammography imaging arrangement for tomosynthesis
CN107427274A (zh) * 2014-02-21 2017-12-01 三星电子株式会社 断层扫描设备及其用于重构断层扫描图像的方法
US10748293B2 (en) 2014-02-21 2020-08-18 Samsung Electronics Co., Ltd. Tomography apparatus and method for reconstructing tomography image thereof
CN110098102A (zh) * 2018-01-31 2019-08-06 帕罗德克斯集团公司 减少成像板扫描器中的部件的校准
CN110098102B (zh) * 2018-01-31 2022-02-18 帕罗德克斯集团公司 减少成像板扫描器中的部件的校准

Also Published As

Publication number Publication date
EP2730226A4 (en) 2015-04-08
US20140119500A1 (en) 2014-05-01
US9380985B2 (en) 2016-07-05
EP2730226A1 (en) 2014-05-14
CN103635139A (zh) 2014-03-12

Similar Documents

Publication Publication Date Title
WO2013005833A1 (ja) Ｘ線撮影装置およびそのキャリブレーション方法
JP5702236B2 (ja) Ｘ線撮影装置およびそのキャリブレーション方法
JP2012045278A (ja) Ｘ線画像撮影装置およびｘ線画像撮影方法
US9636079B2 (en) Motion layer decomposition calibration of x-ray CT imagers
US10761418B2 (en) Imaging method and imaging system
US20110064193A1 (en) Long length multiple detector imaging apparatus and method
JP5019879B2 (ja) Ｘ線ｃｔ装置、画像処理プログラム、及び画像処理方法
JP2006075236A (ja) Ｘ線撮影装置
WO2015022888A1 (ja) 放射線断層像撮影装置
US20110075792A1 (en) Radiography apparatus
JP2005013738A (ja) トモシンセシス用途における対象物を走査するためのシステム及び方法
JP2009022464A (ja) Ｘ線アンギオ撮影装置
JP5702240B2 (ja) Ｘ線撮影装置およびそのキャリブレーション方法
JP2008023241A (ja) デジタルｘ線撮影装置及び方法
JP5016231B2 (ja) 撮像の幾何学的パラメータを決定する方法及び装置
JP6875954B2 (ja) 医用画像診断装置、及び画像処理方法
JP2010187812A (ja) 医用寝台装置
JP2010239997A (ja) 放射線画像撮影装置および方法
JP2011072404A (ja) 放射線撮影システム
JP5564385B2 (ja) 放射線画像撮影装置、放射線画像撮影方法およびプログラム
JP2004081275A (ja) Ｘ線診断装置およびその制御方法
JP5559648B2 (ja) 放射線撮影装置、方法およびプログラム
KR102246090B1 (ko) Fov 확장기능을 갖는 방사선 촬영 장치
JP5676883B2 (ja) Ｘ線ｃｔ装置
JP2010081996A (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: 12807484

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2012807484

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE