JP2012055475A - Radiographic image capturing apparatus, method and program - Google Patents

Abstract

PROBLEM TO BE SOLVED: To prevent the occurrence of a body motion again during a retake operation which is conducted after a body motion has occurred during long-length imaging.SOLUTION: A local motion vector calculation unit 34 calculates local motion vectors in an overlapping area between two adjacent radiographic images. A body motion index value calculation unit 36 calculates a body motion index value using the local motion vectors. Further, a body motion determination unit 38 determines whether or not there is a body motion using the body motion index value. When it is determined that there is a body motion, a retake auxiliary information generation unit 46 generates retake auxiliary information, and an image display unit 60 displays the retake auxiliary information.

Description

  The present invention relates to a radiographic image capturing apparatus, method, and program for capturing a plurality of radiographic images and connecting the plurality of radiographic images to obtain one long captured image.

  2. Description of the Related Art Conventionally, in the medical field and the like, various types of radiation detectors (so-called “Flat Panel Detector”, hereinafter referred to as FPD) that record a radiation image related to a subject by irradiation with radiation that has passed through the subject have been proposed and put into practical use. As such an FPD, for example, there is an FPD using a semiconductor such as amorphous selenium that generates electric charge upon irradiation of radiation, and as such an FPD, a so-called optical reading type or TFT reading type is proposed. ing.

  On the other hand, in simple X-ray photography, there is a case where long photographing is performed on a long region such as the entire spine (whole spine) or the entire foot (whole leg), and long photographing using FPD is also performed. Has been. Here, since the FPD may have a narrower shooting range than the subject to be shot, in order to perform long shooting, the same subject is moved each time the FPD is moved along a predetermined movement axis and the position is changed. To receive radiation that has passed through. Then, a reading operation is performed from the FPD for each radiation irradiation (recording of a radiation image), and image data indicating a radiation image is acquired for each reading operation. Then, image data indicating a long portion of the subject can be obtained by combining the image data so as to be connected later.

  By the way, in the long shooting using the FPD, since the shot interval is 3 to 5 seconds, there is a possibility that the body movement of the subject may occur between shots. Since the line images cannot be properly connected and accurate measurement is hindered, it is necessary to re-photograph. In addition, even if body movement occurs in the middle, there are cases in which it is not noticed for the first time after seeing the composite image, and it is inefficient to re-shoot after shooting to the end. Furthermore, since the shot after the body movement is meaningless, the exposure of the subject also increases.

  For this reason, a method has been proposed in which the body movement of the subject is detected to stop photographing or warn of the body movement. For example, Patent Document 1 proposes a method of detecting a body motion of a subject at the time of photographing using a sensor and giving a warning when a body motion is detected. Patent Document 2 proposes a method for calculating the amount of body motion of a subject in consideration of an assembly error of an FPD and an installation error of a radiographic image capturing apparatus. Specifically, by performing template matching in the overlapping area of multiple radiographic images, a plurality of local shift amounts are obtained, the average value of the shift amounts is calculated as a body movement amount, and this is compared with a threshold value. Thus, a method for detecting the presence or absence of body movement has been proposed. Further, Patent Documents 1 and 2 also propose a method for correcting body motion by nonlinear warping.

JP 2009-240656 A Japanese Patent Application No. 2010-074023

  The method described in Patent Document 1 gives a warning when body movement occurs, and the operator re-photographs when there is a warning. However, since it is not possible to know how the body movement has occurred only by the warning, there is a possibility that the body movement may occur again when re-imaging is performed. If body movement occurs again at the time of re-shooting as described above, it is necessary to take another shot, resulting in poor shooting efficiency and an increased exposure of the subject.

  The present invention has been made in view of the above circumstances, and it is an object of the present invention to prevent reoccurrence of body movement when performing re-photographing after occurrence of body movement when performing long imaging.

A radiographic imaging device according to the present invention moves a radiation detector and irradiates the radiation detector with radiation that has passed through a subject each time the position is changed by the movement, so that a plurality of radiations at least partially overlapped with each other. A radiographic imaging device for acquiring an image,
An imaging means for moving the radiation detector along a predetermined movement axis and irradiating the radiation detector with radiation transmitted through the subject each time the position is changed by the movement;
A radiological image acquisition means for acquiring a plurality of radiographic images of the subject by reading a signal from the radiation detector each time the movement and the irradiation of the radiation are performed;
Body motion detection means for detecting body motion of the subject at the time of capturing the plurality of radiation images;
And an information generating unit configured to generate re-imaging assistance information for assisting re-imaging of the plurality of radiographic images when the body motion is detected.

  The radiographic image capturing apparatus according to the present invention may further include display means for displaying the re-imaging auxiliary information.

  The radiographic image capturing apparatus according to the present invention may further include re-imaging control means for controlling the re-imaging based on the re-imaging auxiliary information.

  In the radiographic image capturing apparatus according to the present invention, the re-imaging control unit may be a unit that divides the plurality of radiographic images into frames while avoiding a position where the body movement of the subject occurs.

  In the radiographic image capturing apparatus according to the present invention, the re-imaging control means may be means for increasing the fixing strength of the subject fixing tool.

Further, in the radiographic imaging device according to the present invention, the body movement detection unit includes a local movement vector calculation unit that calculates a local movement vector that represents a local shift of the subject in an overlapping region in the plurality of radiographic images. ,
Body motion index value calculating means for calculating the body motion index value based on the local movement vector;
Body movement determining means for determining the presence or absence of the body movement based on the body movement index value may be provided.

  In this case, the body motion index value calculating means may be a means for calculating an index value of the parallel movement amount of the subject as the body motion index value.

  Further, in this case, the body motion index value calculation means further uses at least one of the index value of the three-dimensional movement amount of the subject and the index value of the two-dimensional movement amount of the subject as the body motion index value. It is good also as a means to calculate.

  In the radiographic image capturing apparatus according to the present invention, the body motion detection means may be a sensor that detects the movement of the subject.

In the radiographic imaging method according to the present invention, the radiation detector is moved, and each time the position is changed by the movement, the radiation that has passed through the subject is irradiated to the radiation detector, and a plurality of radiations at least partially overlapped with each other. A radiographic imaging method for acquiring an image,
Moving the radiation detector along a predetermined movement axis, and irradiating the radiation detector with radiation transmitted through the subject each time the position is changed by the movement;
Acquiring a plurality of radiation images of the subject by reading a signal from the radiation detector each time the movement and irradiation of the radiation are performed;
Detecting body movement of the subject at the time of capturing the plurality of radiation images;
Generating re-imaging assistance information for assisting re-imaging of the plurality of radiographic images when the body movement is detected.

  In addition, you may provide as a program for making a computer perform the radiographic imaging method by this invention.

  According to the present invention, since the re-shooting auxiliary information is generated when the body movement is detected, the operator performs the shooting according to the re-shooting instruction information so that the body movement does not occur. be able to. Therefore, re-photographing can be performed efficiently, and the exposure amount to the subject can be reduced.

  In particular, by displaying the re-shooting auxiliary information, the operator can confirm the re-shooting auxiliary information at a glance, so that the re-shooting can be performed more efficiently.

  In addition, by performing re-shooting control based on the re-shooting auxiliary information, re-shooting can be performed without causing any body movement without bothering the operator.

  In addition, for each local region in the overlapping region in a plurality of radiographic images, a local movement vector representing a subject displacement in the overlapping region is calculated, a body motion index value is calculated based on the local movement vector, and the body motion index value is calculated. By determining the presence or absence of body movement based on this, body movement can be detected by a relatively simple calculation.

  In addition, by calculating the index value of the parallel movement amount of the subject as the body motion index value, it is possible to detect the body motion due to the parallel movement of the subject.

  Further, by calculating at least one of the index value of the three-dimensional movement amount of the subject and the index value of the two-dimensional movement amount of the subject as the body motion index value, a three-dimensional body of twisting and falling of the subject is calculated. It is possible to detect two-dimensional body movements such as movement and rotational movement and enlargement / reduction of a subject.

Schematic which shows the structure of the radiographic imaging apparatus by the 1st Embodiment of this invention. Diagram for explaining calculation of local movement vector Diagram for explaining template matching using multi-resolution conversion Diagram for explaining another method of template matching Diagram for explaining another method of template matching Figure showing an example of a histogram The figure for demonstrating calculation of the body movement index value of a three-dimensional movement Diagram for explaining body movement correction The flowchart which shows the process performed in 1st Embodiment. The figure which shows the example of a display in the display screen of the image display part in the case of no body movement The figure which shows the other example of the display in the display screen of the image display part in the case of no body movement The figure which shows the example of a display in the display screen of the image display part in the case of body movement The figure which shows the other example of the display in the display screen of the image display part in the case of body movement The figure which shows the example of a display in the display screen of the image display part in the case of body movement The figure which shows the other example of the display in the display screen of the image display part in the case of body movement Schematic which shows the structure of the radiographic imaging apparatus by the 2nd Embodiment of this invention. The flowchart which shows the process performed in 2nd Embodiment.

  Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a schematic diagram showing the configuration of a radiographic imaging apparatus according to the first embodiment of the present invention. As shown in FIG. 1, the radiographic image capturing apparatus 150 according to the first embodiment sequentially uses a single radiation source 100 and a single FPD 110 to sequentially form a plurality of adjacent areas N1, N2,. A long radiographing is performed in which a long radiographic image showing a large part of the subject N is obtained by synthesizing a plurality of radiographic images acquired by the radiographing, thereby obtaining a long radiographic image showing the majority of the subject N.

  Of course, the radiographic image capturing apparatus 150 according to the present embodiment can be used not only for long-length imaging but also for normal imaging in which only a specific part of the subject N such as the chest and lower limbs is imaged. In this case, energy subtraction imaging and temporal subtraction imaging can be performed, but only the long imaging will be described in detail below.

  Specifically, the radiographic imaging device 150 receives the radiation 104 that emits the radiation 104 from the exit window 111 to the irradiation range defined by the collimator 112 and the radiation 104 that has passed through the subject N, and detects the radiation 104. An FPD 110 having an imaging surface (radiation detection surface) 102, a detector moving unit 20 that moves the FPD 110 along the subject N, and a line for arranging the radiation source 100 so that the position and orientation of the exit window 111 are in a desired state. A source arrangement unit 25. In FIG. 1, the radiation center axis of the radiation 104 whose irradiation range is determined by the collimator 112 is indicated by Cr.

  The FPD 110 detects the radiation 104 transmitted through the subject N, converts it into an electrical signal, and outputs image data representing the radiation image of the subject N. Note that the FPD 110 can use either a direct FPD that directly converts radiation into electric charges, or an indirect FPD that converts radiation once into light and further converts the converted light into an electrical signal. The direct type FPD includes a photoconductive film such as amorphous selenium, a capacitor, a TFT (Thin Film Transistor) as a switch element, and the like. For example, when radiation such as X-rays is incident, electron-hole pairs (e-h pairs) are emitted from the photoconductive film. The electron-hole pair is accumulated in the capacitor, and the electric charge accumulated in the capacitor is read out as an electric signal through the TFT.

  On the other hand, an indirect FPD is composed of a scintillator layer made of a phosphor, a photodiode, a capacitor, a TFT, and the like. For example, when radiation such as “CsI: Tl” is incident, the scintillator layer emits light (fluoresces). Light emitted by the scintillator layer is photoelectrically converted by a photodiode and accumulated in a capacitor, and the electric charge accumulated in the capacitor is read out as an electrical signal through the TFT.

  The detector moving unit 20 stands up from the floor surface 5F in the vertical direction (in the direction of the arrow Y in the figure) and moves to move the FPD 110 in the vertical direction, which is the long direction, and the two columns 21 holding the FPD 110 therebetween. And a mechanism 22. As the moving mechanism 22, a mechanism in which the FPD 110 is supported by a conventionally known linear slide mechanism or the like and moved using a driving source such as a motor can be employed.

  When performing imaging for acquiring a radiographic image to be used for synthesis, the subject N is arranged along the moving direction of the FPD 110. That is, shooting is performed with the subject N standing on the floor.

  The radiation source arrangement unit 25 holds and moves the radiation source 100 so as to face the imaging surface 102 of the FPD 110 with the subject N in between, that is, to face the direction of the FPD 110. The radiation source arrangement portion 25 is engaged with the support column 26 extending vertically from the ceiling 5E, a ceiling base table 27 that moves the support column 26 along the ceiling 5E in the arrow Z direction in the figure, and the support column 26. It has a turntable 28 that can move in the middle arrow Y direction and can rotate about an axis perpendicular to the paper surface. The radiation source 100 is mounted on the turntable 28. Thereby, the radiation source 100 is movable in the vertical direction (arrow Y direction in the figure) and the horizontal direction (arrow Z direction in the figure), and has an axis parallel to the X axis in the figure passing through the approximate center of the radiation source 100. It can be rotated around. In addition, the radiation source arrangement | positioning part 25 can also be formed using drive sources, such as a linear slide mechanism and rotation mechanism known conventionally, and a motor.

  The radiographic image capturing apparatus 150 includes a long imaging control unit 50 that controls operations of the detector moving unit 20 and the radiation source arranging unit 25. The long imaging control unit 50 controls the operation of the detector moving unit 20 to sequentially move the FPD 110 to each position Q1, Q2,... For performing radiography along the direction in which the subject N extends. At the same time, the long imaging control unit 50 controls the operation of the radiation source arrangement unit 25 so that the irradiation direction of the radiation 104 emitted from the radiation source 100 faces the direction of the imaging surface 102 of the FPD 110 arranged at each position. The radiation source 100 is arranged as described above. Then, when the radiation source 100 is driven, adjacent regions N1, N2,... In the subject N are sequentially photographed, and a plurality of partial radiations for representing the entire subject N at each photographing. Image data representing an image is acquired.

  The radiographic image capturing apparatus 150 further includes a body motion detecting unit 30 that detects body motion of the subject N at the time of capturing, a body motion correcting unit 40 that corrects the shift amount of the subject N based on the body motion, and the above-described radiography. Each of the obtained image data is synthesized to include an image synthesis unit 42 that generates a long radiation image representing the entire subject N, a warning unit 44, and a re-imaging auxiliary information generation unit 46. The long radiation image generated by the image composition unit 42 is displayed on the image display unit 60 including, for example, a CRT display device, a liquid crystal display device, or the like.

  The body motion detection unit 30 includes an image acquisition unit 31, an imaging information acquisition unit 32, a local movement vector calculation unit 34, a body motion index value calculation unit 36, and a body motion determination unit 38.

  The image acquisition unit 31 includes various interfaces for acquiring a radiation image from the FPD 110. Note that when the body motion detection unit 30 is provided separately from the radiographic image capturing apparatus 150 and these are connected via a network, the image acquisition unit 31 serves as a network interface.

  The imaging information acquisition unit 32 acquires imaging information from the console 70 that controls the overall operation of the radiographic imaging device 150. Here, as the shooting information, information related to the shooting status and information related to the shooting target are acquired. Examples of the information regarding the imaging state include information such as an imaging time required for imaging all radiographic images during long imaging, an imaging interval between two adjacent radiographic images, and whether or not the subject N is fixed during imaging. . The information about the subject to be photographed includes information such as a photographing part of the subject N (chest, lower limbs, whole spine, all lower limbs, etc.) and a symptom of the subject N (whether severe or mild, preoperative or postoperative). It is done. The shooting information acquisition unit 32 acquires at least one of these pieces of information as shooting information.

  The imaging information may be information input by the operator to the console 70. The console 70 measures the imaging time, detects whether the subject N is fixed, and the acquired radiographic image. You may make it acquire what the console 70 calculated automatically by performing site | part recognition. Also, since the shooting information differs depending on the shooting menu, the shooting information automatically determined by the console 70 according to the shooting menu selected by the operator may be acquired.

  The local movement vector calculation unit 34 calculates a local movement vector in an overlapping region where two adjacent radiographic images overlap each other. FIG. 2 is a diagram for explaining the calculation of the local movement vector. As shown in FIG. 2, the local movement vector calculation unit 34 performs the local movement amount of the image that appears due to the body movement of the subject N, that is, local movement, in the overlapping regions K1 and K2 of the two adjacent radiation images S1 and S2. The same vector as the one image (template T with reference to the control point existing in the overlapping region K1 of the radiation image S1) is present in any part of the other image (overlapping region K2 of the radiation image S2). It is calculated by template matching that searches whether to do.

  Specifically, a correlation value between the template T and the search target image I of the same size sequentially searched from the overlapping region K2 of the radiation image S2 is calculated within a predetermined search range R. Thereby, a correlation distribution having the same size as the search range R is obtained. Then, in the correlation distribution, the pixel shift amount (movement amount) from the reference position of the template T when the correlation value becomes maximum (the position where the correlation becomes maximum when there is no body movement) is defined as the local movement vector V0. calculate.

  In the present embodiment, a plurality of control points are set. Thereby, a plurality of local movement vectors are calculated in the overlapping region K1. The control points may be all pixel positions of the image in the overlapping region K1, or may be pixel positions that are thinned out at a predetermined pixel interval. Further, it may be a representative pixel position such as an intersection of edges in the overlapping region K1 or a pixel position having a large variance. Here, the representative pixel position may be detected automatically, or may be manually set by the operator by displaying the overlapping region K1 of the radiation image S1.

  Further, not only the local movement vector V0 is calculated based on the position where the correlation value becomes maximum, but may be calculated based on the barycentric position of the correlation value in the correlation distribution. The barycentric position of the correlation value is obtained by calculating a weighted average position based on the correlation value of the pixel position in the search range R with reference to the origin in the search range R (for example, the reference position of the template described above). it can. At this time, the center-of-gravity position may be calculated using only pixel positions where the correlation value is equal to or greater than a predetermined threshold value. This barycentric position is referred to as a highly correlated barycentric position.

  Further, when calculating the local movement vector V0, as shown in FIG. 3, by performing multi-resolution conversion on the images of the overlapping regions K1 and K2, a plurality of overlapping region images having different resolutions are generated, and overlapping of each resolution is performed. The local movement vector V0 may be calculated sequentially from the low resolution to the high resolution between the region images. Note that FIG. 3 shows a state in which multi-resolution conversion is performed up to 1/4 resolution by performing resolution conversion three times. Also, in FIG. 3, the overlapping area images before multi-resolution conversion are denoted as K1-1 and K2-1, the overlapping area images of the next resolution are K1-2 and K2-2, and the overlapping area images of the next resolution are further displayed. K1-3 and K2-3.

  Here, the overlapping area images of each resolution differ by 1/2. Therefore, when the resolution conversion is performed three times and the multi-resolution conversion is performed up to a resolution of 1/4, if 16 local movement vectors V0 are calculated, first, the overlapping region images K1-3 and K2 with the lowest resolution are calculated. -3 is used to calculate one local movement vector V0. Note that, in the portion explaining the state where the local movement vector V0 of FIG. 3 is calculated, the size of the overlapping region image is the same size as the overlapping regions K1 and K2.

  Next, four local movement vectors V0 are calculated using the overlapping region images K1-2 and K2-2 of the next higher resolution. At this time, the local movement vector V0 can be efficiently calculated by using the local movement vector V0 calculated in the overlapping area images K1-3 and K2-3. For example, as illustrated in FIG. 3, when the local movement vector V0 that is directed diagonally right upward is calculated in the overlapping region images K1-3 and K2-3, in the overlapping region images K1-2 and K2-2 of the next resolution. If the local movement vector V0 calculated at the previous resolution is set as an initial value, template matching is performed only in a region near the upper right corner of the search range R (shown by diagonal lines) as shown in FIG. Since the calculation amount for calculating the correlation value can be reduced, the local movement vector V0 can be calculated efficiently.

  In the above, when performing template matching, the search range R is set in the overlap region K2, but as shown in FIG. 5, the search range R may be set beyond the overlap region K2. . Thereby, since the search range can be enlarged, the local movement vector V0 can be calculated with higher accuracy.

  In the above description, the template T is set in the overlapping region K1 of the radiation image S1 and the local movement vector V0 is calculated. In addition to this, the template T is set in the overlapping region K2 of the radiation image S2, and the overlapping is performed. You may make it calculate the local movement vector V0 on the basis of the area | region K2.

  In the above, the local movement vector V0 is calculated using the correlation value. However, the index is not limited to the correlation value as long as it is an index representing the similarity between the template T and the search target image I in the search range R. For example, any other index such as residual and mean square error can be used.

  Here, the local movement vector calculation unit 34 can calculate the local movement vector V0 using various methods as described above. In the present embodiment, the shooting information acquired by the shooting information acquisition unit 32 is used. Based on the above, the local movement vector V0 is calculated. Specifically, the method of calculating the local movement vector V0 is switched according to the shooting information. For example, when the information on the imaging status includes information such as the number of images to be taken and the imaging time is long, the imaging interval is long, and the subject is not fixed, or the imaging target information is the chest part (if the imaging is not long) If the patient, who is the subject, who is the subject or the entire spine (in the case of long imaging) is serious and cannot hold the body stationary during imaging, the body motion is assumed to be large. For this reason, when it is assumed that the body movement is large from the imaging information, the local movement vector calculation unit 34 increases the search range R or performs template matching by performing multi-resolution conversion on the overlapping region. Perform multi-resolution conversion to low resolution. Thereby, the local movement vector V0 can be calculated with higher accuracy. Note that whether or not the shooting time is long and whether or not the shooting interval is long may be determined using a threshold value.

  Conversely, if the shooting information includes information that is assumed to have a small body movement, the search range R should be reduced, the number of multi-resolution conversions should be reduced, or multi-resolution conversion should not be performed. That's fine. Here, the local movement vector calculation unit 34 may not only determine whether the body movement is large, but may also determine the magnitude of the body movement step by step. In this case, the size of the search range R may be changed stepwise according to the size of body movement, or the number of multi-resolution conversions may be changed stepwise.

  In addition, when the imaging information includes information of imaging regions (for example, the chest, the entire spine, etc.) that are assumed to have many two-dimensional structures, the local movement vector is based on the position where the correlation value is maximized during template matching. In the case where the calculation unit 34 is to be calculated and an imaging region (for example, the lower limbs, all lower limbs, etc.) assumed to have a large one-dimensional structure, the correlation value is likely to be shaken along the one-dimensional direction. The local movement vector V0 may be calculated based on the correlation centroid position.

  Here, the local movement vector calculation unit 34 calculates the local movement vector V0 based on the photographing information as described above, but calculates the local movement vector V0 using any one of the above-described techniques. It may be determined in advance. In this case, the photographing information acquisition unit 32 is not necessary.

  The body motion index value calculation unit 36 calculates a body motion index value using the plurality of local movement vectors V0 calculated by the local movement vector calculation unit 34. First, the body motion index value calculation unit 36 calculates a body motion index value for translation. The parallel movement is a linear movement parallel to the imaging surface 102 of the FPD 110 of the subject N. For this reason, the body motion index value calculation unit 36 creates a three-dimensional histogram having the vertical shift, the horizontal shift, and the frequency as three axes for the plurality of local movement vectors V0 calculated by the local movement vector V0. Here, the vertical shift and the horizontal shift correspond to the vertical and horizontal directions of the overlapping area of the radiographic image. By using the coordinate system shown in FIG. 1, the local movement vector V0 has the horizontal direction as the X coordinate and the vertical direction as the vertical direction. It is represented by a two-dimensional vector on the XY coordinate as the Y coordinate. Therefore, the vertical shift and the horizontal shift are respectively the magnitudes of the local movement vector V0 in the Y direction and the X direction.

  FIG. 6 is a diagram showing an example of a histogram. As shown in FIG. 6, in the histogram, the frequency of the local movement vector V0 according to the values of vertical shift and horizontal shift is represented. When creating a histogram, only the local movement vector V0 whose correlation value serving as a reference for calculating the local movement vector V0 is a predetermined threshold value or more may be used. Also, the variance value of the pixel values in the template T used when calculating the local movement vector V0 is calculated, and the local value calculated using the template T having a variance value smaller than a predetermined threshold and many existing edges. Only the movement vector V0 may be used.

  Then, the body motion index value calculation unit 36 calculates a body motion index value based on the histogram. As a method for calculating the body motion index value, a local movement vector V0 having the maximum histogram frequency is determined, and the position (frequency maximum position) of the local movement vector V0 on the histogram and the reference position (that is, the origin of the histogram) are determined. Can be used as a body motion index value. Alternatively, a method may be used in which the centroid position (frequency centroid position) of the local movement vector V0 having a frequency equal to or higher than a predetermined threshold is calculated, and the distance between the frequency centroid position and the reference position is calculated as a body movement index value. Good. In addition, the average value of the vertical shift and the horizontal shift of the local movement vector V0 are calculated to calculate the average position of the local movement vector V0, and the distance between the average position and the reference position is calculated as the body movement index value. May be. Further, the median of the vertical displacement and the lateral displacement of the local movement vector V0 are calculated to calculate the central position of the local movement vector V0, and the distance between the central position and the reference position is calculated as the body movement index value. May be.

  Further, the body motion index value calculation unit 36 calculates a body motion index value of the three-dimensional movement of the subject N. Note that the three-dimensional movement is a twisting movement of the subject N and a tilting movement in the front-rear direction with respect to the imaging surface 102 of the FPD 110. The body motion index value of the three-dimensional motion is calculated using the histogram of the local movement vector V0 calculated as described above. FIG. 7 is a diagram for explaining calculation of a body motion index value of a three-dimensional motion. Here, when the body movement is only parallel movement, the distribution of the histogram is concentrated in one place, but when the body movement is a three-dimensional movement, the distribution of the histogram is widened. For this reason, the body motion index value calculation unit 36 calculates the standard deviation or variance of the distribution of the histogram as the body motion index value of the three-dimensional movement. Note that the larger the standard deviation or variance, the greater the three-dimensional movement. As the center for calculating the standard deviation or variance, any of the above-described maximum frequency position, frequency center of gravity position, average position, and center position can be used.

  Further, the body motion index value calculation unit 36 calculates a body motion index value of the two-dimensional motion of the subject N. Note that the two-dimensional movement is a movement of rotation of the subject N on a plane parallel to the imaging surface 102 of the FPD 110, translation in the front-rear direction of the imaging plane 102, and translation in a direction parallel to the imaging plane 102. . In addition, since the parallel movement of the imaging surface in the front-rear direction appears as a body movement of enlargement / reduction in the radiographic image, the movement in the front-rear direction will be described as a body movement of enlargement / reduction in the following description. Further, as described above, the movement of the parallel movement can be calculated using the local movement vector V0. Here, the two-dimensional movement is described as a movement including all of rotation, enlargement / reduction, and parallel movement. .

The body motion index value calculation unit 36 calculates the two-dimensional movement of the subject N based on the formula described in, for example, “Graphic Processing Engineering, Fujio Yamaguchi, Nikkan Kogyo Shimbun, 68-82, 1981”. A body movement index value is calculated. Specifically, the local movement vector V0 calculated as described above is applied to the expression of the two-dimensional affine transformation, and each element of translation, rotation, and enlargement / reduction is calculated from the plurality of local movement vectors V0 using the least square method. Is calculated as a body movement index value. Here, in the two-dimensional affine transformation, when the horizontal movement amount is tx and the vertical movement amount is ty, the parallel movement is expressed by the following equation (1). In Equation (1) and the following equations, (x, y) is a two-dimensional coordinate position before movement, and (x ^* , y ^* ) is a two-dimensional coordinate position after movement.

In the two-dimensional affine transformation, if the rotation angle is θ, the rotational movement is expressed by the following equation (2).

In the two-dimensional affine transformation, if the enlargement / reduction parameters for enlargement / reduction are a (X direction) and d (Y direction), the enlargement / reduction is expressed by the following equation (3).

Summarizing formulas (1 ′) to (3 ′),

H is a matrix of nine elements. Here, since the local movement vector V0 can be expressed as (x ^* −x, y ^* −y), a plurality of local movement vectors V0 are applied to the above equation (4), and nine are obtained by the least square method. By obtaining these elements, the parallel movement amount, the rotation angle, and the enlargement / reduction parameter can be calculated. The parallel movement amount, the rotation angle, and the enlargement / reduction parameter obtained in this way become the body motion index value of the two-dimensional movement. In the present embodiment, in the subsequent processing, the parallel movement body motion index value is calculated based on the histogram as described above, and the two-dimensional motion body motion index value is Only the rotation angle and the enlargement / reduction parameters are used.

  Here, the body motion index value calculation unit 36 can calculate the body motion index value using various methods for the parallel movement and the three-dimensional motion as described above, but in the present embodiment, The body motion index value is calculated based on the imaging information acquired by the imaging information acquisition unit 32. Specifically, the method for calculating the body motion index value is switched according to the photographing information. For example, regarding the body movement index value of translation, when the imaging region is the lower limb (not long imaging) or all the lower limbs (long imaging), the body movement of the subject N such as gas accumulated in the heart or intestine Since the moving structure is not included separately, if the imaging information includes information that the imaging region is the lower limb, the body motion index value is calculated using the average position and the central position, and as much as possible An index value using the local movement vector V0 is calculated. As for the body motion index value of the three-dimensional motion, the body motion index value is calculated by using the average position and the center position as the center when calculating the standard deviation or the variance.

  On the other hand, when the imaging region is the chest (not long imaging) or the entire spine (long imaging), a structure that moves separately from the body movement of the subject N, such as gas accumulated in the heart or intestine, is included. Therefore, when the imaging information includes information indicating that the imaging region is the chest, the body motion index value for the parallel movement is calculated using the frequency maximum position or the frequency center of gravity position. Thus, the index value can be calculated so as not to be affected by local movement. As for the body motion index value of the three-dimensional motion, the body motion index value is calculated by using the frequency maximum position or the frequency centroid position as the center when calculating the standard deviation or variance.

  Note that the body motion index value calculation unit 36 calculates all body motion index values of parallel movement, three-dimensional motion, and two-dimensional motion, but the parallel motion, three-dimensional motion, and two-dimensional motion are calculated. At least one body motion index value may be calculated. However, it is preferable to include a body motion index value of translation.

  Here, the body motion index value calculation unit 36 calculates the body motion index value based on the imaging information as described above. However, the body motion index value is calculated using any one of the above methods. You may make it decide beforehand whether to do. In this case, the photographing information acquisition unit 32 is not necessary.

  The body movement determination unit 38 compares the body movement index value calculated by the body movement index value calculation unit 36 with a threshold value, determines that there is body movement when the body movement index value is equal to or greater than the threshold value, and determines Output the result. In addition, when a plurality of body motion index values of parallel movement, three-dimensional motion, and two-dimensional motion are calculated, the body motion index values for each of them are compared with a threshold value corresponding to each. Then, the presence or absence of body movement is determined. In this case, if any one of the body motion index values is equal to or greater than the threshold value, it is determined that there is a body motion. Further, when any two body motion index values are equal to or greater than the threshold value, or when all the body motion index values are equal to or greater than the threshold value, it may be determined that there is a body motion. Further, the presence or absence of body movement may be individually determined for each of the parallel movement, the three-dimensional movement, and the two-dimensional movement. Further, for two-dimensional movement, threshold values may be prepared for each of rotation and enlargement / reduction, and the presence or absence of body movement may be determined for each of rotation and enlargement / reduction.

  Here, in the present embodiment, the body movement determination unit 38 determines the presence or absence of body movement by comparing the body movement index value with a threshold value based on the shooting information calculated by the shooting information acquisition unit 32. . Specifically, the magnitude of the threshold is switched according to the shooting information. For example, when the imaging region is the lower limb (not long imaging) or the entire lower limb (long imaging), a structure that moves separately from the body movement of the subject N, such as gas accumulated in the heart or intestine, is not included. On the other hand, when the imaging region is the chest (not long imaging) or the entire spine (long imaging), a structure that moves separately from the body movement of the subject N, such as gas accumulated in the heart or intestine, is included. . For this reason, when the imaging information includes information indicating that the imaging region is the chest, the threshold for determination is set larger than when the imaging region is the lower limb. Thereby, it is possible to determine the presence or absence of body movement without being affected by local movement.

  In addition, when the imaging part is an imaging part with a small amount of detected body movement, such as when the imaging part is the entire spine (long imaging), the threshold value is set small, and the imaging part is set to all the lower limbs (long). In the case of an imaging region where the detected body movement amount is large, such as in the case of (scale shooting), the threshold value is set large.

  Note that the presence or absence of body movement may be determined using a predetermined threshold without being based on the photographing information. In this case, the photographing information acquisition unit 32 is not necessary.

  Next, correction of body movement by the body movement correction unit 40 will be described. The body motion correction unit 40 detects image distortion caused by body motion of the subject N based on the body motion index value detected by the body motion detection unit 30 when the determination result by the body motion determination unit 38 is no body motion. Correction for eliminating the above is performed on the radiation image. When the determination result is that there is a body motion, the correction is not performed because the image is too distorted when the body motion is corrected. Here, the case where the image distortion of two adjacent radiographic images is corrected will be described. As shown in FIG. 8A, when the body motion index value of the parallel movement is detected, FIG. ), Body motion correction is performed by relatively shifting the two radiation images S1 and S2 based on the body motion index value, that is, the local movement vector V0. In addition, when a plurality of body motion index values are detected, the position of the template for which the amount of deviation is detected matches the position of the search target image, that is, the position of each of the local movement vectors V0 is positioned. Body motion correction is performed by relatively nonlinear warping of a plurality of radiation images.

  In the case where all of the body movement index values of the parallel movement, the three-dimensional movement, and the two-dimensional movement are calculated, all these body movement index values are not 0 (that is, the parallel movement, the three-dimensional movement, and the When at least two of the two-dimensional movements exist), body motion correction may be performed by nonlinear warping. Further, when only one of the body motion index values of the parallel movement, the three-dimensional motion, and the two-dimensional motion has a value, the body motion correction according to the motion may be performed. . For example, when only the translational movement index value has a value, it is only necessary to correct the body movement by shifting the radiation images from each other, and when only the three-dimensional movement index value has a value. The body motion may be corrected by nonlinear warping. Further, when only the body motion index value of the two-dimensional motion has a value, the body motion index value (that is, the rotation and enlargement / reduction elements other than the translation in Expressions (1) to (3)) is used. The body motion may be corrected by the two-dimensional affine transformation. In addition, when only one of the rotation and the enlargement / reduction among the body motion index values of the two-dimensional motion has a value, the body motion correction is performed by two-dimensional affine transformation using only one of the body motion index values. Can be done. Further, when a two-dimensional movement is generated, nonlinear warping may be performed. Furthermore, body motion correction may be performed according to the motion that is the largest body motion index value among the three-dimensional motion and the two-dimensional motion body motion index values.

  The image synthesizing unit 42 generates a synthesized image C1 by synthesizing the radiographic images corrected for body movement so as to be joined together.

  The warning unit 44 issues a warning when the body movement is large as will be described later.

  The re-shooting auxiliary information generating unit 46 generates re-shooting auxiliary information for assisting re-shooting when the body movement determining unit 38 determines that there is a body movement. Here, in the case of performing long shooting, re-shooting is required when the body movement between each shot is large. For this reason, the re-imaging auxiliary information generation unit 46 performs re-imaging with information for specifying two adjacent radiographic images in which body movement has occurred and body movement reduced based on the determination result by the body movement determination unit 38. Therefore, information to be displayed on the image display unit 60 is generated as re-shooting auxiliary information. The information displayed on the image display unit 60 includes information such as text, images (including still images and moving images), audio, and the like that prompts the user to fix the location where the subject N has moved. . In addition, when it is determined that there is a body motion, it is information that recommends one-shot imaging using a single large FPD or CR long cassette instead of long imaging that captures a plurality of radiation images. There may be. Note that text and images that are created in advance and stored in a storage unit (not shown) may be used.

  The overall operation of the radiation image capturing apparatus 150 is controlled by the console 70. Therefore, information relating to the subject N, imaging conditions for obtaining a long radiation image, and the like are input to the console 70, and the information sets a radiation irradiation range defined by the long imaging control unit 50 and the collimator 112. The image is input to a photographing adjustment unit (not shown) or the like. This imaging adjustment unit adjusts the position of the radiation source 100, the state of the collimator 112, the position of the FPD 110, etc. at the time of each radiation imaging so that a composite radiation image of a predetermined size is obtained every four times of radiation imaging, for example. Perform frame allocation. Thereafter, in response to a command from the console 70, an operation for capturing four radiation images is executed.

  Here, in order to determine the size of the radiographic image captured a plurality of times, in addition to defining the radiation irradiation range by the collimator 112 or the like as described above, a part of the radiographic image obtained by each imaging is cut out and each image is extracted. You may make it adjust the length and width of a part.

  Next, processing performed in the first embodiment will be described. FIG. 9 is a flowchart showing processing performed in the first embodiment. First, long imaging is performed while moving the FPD 110, and a radiation image at each position of movement is acquired (step ST1). Then, the local movement vector calculation unit 34 calculates a local movement vector V0 in an overlapping region of two adjacent radiographic images (step ST2), and the body movement index value calculation unit 36 calculates a body movement index value (step ST3). ). Subsequently, the body movement determination unit 38 determines the presence or absence of body movement (step ST4).

  When there is no body movement, the body movement determination unit 38 outputs the determination result to the body movement correction unit 40 and the image display unit 60 together with the body movement index value (step ST5). The body motion correction unit 40 corrects the body motion of the plurality of radiographic images based on the body motion index value (step ST6), and the image synthesis unit 42 synthesizes the radiographic images after the body motion correction to generate a composite image C1. (Step ST7). Then, the image display unit 60 displays the composite image C1 together with the body movement index value (step ST8), and the process ends.

  FIG. 10 is a diagram illustrating an example of display on the display screen of the image display unit 60. As shown in FIG. 10, the display screen 61 displays an image display area 62 for displaying the composite image C1 and a body movement index value display area 63 for displaying the body movement index value. When a plurality of body motion index values are calculated, the body motion index value is displayed in the body motion index value display area 63 in association with each of parallel movement, three-dimensional motion, and two-dimensional motion. That's fine. Further, only the body movement index value of the parallel movement may be displayed, or only the largest body movement index value among the body movement index values of the parallel movement, the three-dimensional movement, and the two-dimensional movement may be displayed. Good. In this case, since the body motion index values of the parallel movement, the three-dimensional motion, and the two-dimensional motion have different units, the unit may be normalized to determine the size of the body motion index value. As for the body motion index value of the two-dimensional movement, only the body motion index value of rotation or enlargement / reduction may be displayed, or all the body motion index values may be displayed, and the largest body Only the dynamic index value may be displayed. As shown in FIG. 11, when the cursor is brought close to the overlapping area of the composite image C1 displayed in the image display area 62, the body movement index value may be displayed in a pop-up manner.

  On the other hand, when the determination in step ST4 is that there is body movement, the body movement determination unit 38 outputs the determination result to the warning unit 44 and the re-shooting auxiliary information generation unit 46 (step ST9). The reshooting auxiliary information generation unit 46 generates reshooting auxiliary information (step ST10) and outputs it to the image display unit 60. The image display unit 60 displays re-shooting auxiliary information (step ST11). On the other hand, the warning unit 44 issues a warning by a voice message (including a warning voice) or a buzzer sound (warning sound) (step ST12), and ends the process. In this case, the voice message may be a message that prompts the operator to re-shoot, for example, “Because body movement is large, please take a picture again”.

  FIG. 12 is a diagram illustrating an example of display of reshooting auxiliary information. As shown in FIG. 12, an image display area 64 schematically showing a plurality of radiographic images (four in this case) acquired by long imaging is displayed on the display screen 61 of the image display unit 60. As shown in FIG. 12, in the four radiographic images S <b> 1 to S <b> 4 displayed on the display screen 61, an arrow is given as re-imaging auxiliary information at a position where a body movement has occurred. In FIG. 12, since an arrow is provided between the radiographic image S1 and the radiographic image S2, it can be seen that body movement has occurred between the radiographic image S1 and the radiographic image S2. Further, symbols or the like may be used instead of the arrows. In FIG. 12, a text 65 “A body movement has occurred in the first and second images” is displayed as the re-shooting assistance information. By viewing the re-shooting auxiliary information, the operator urges attention to the subject N so that body movement does not occur between the first and second images during re-shooting, and the first and second images. A portion corresponding to the first sheet can be fixed by a fixing tool such as a belt, or a cushion can be inserted and fixed.

  Note that, as shown in FIG. 13, text may be displayed in a pop-up when the cursor is brought close to a location where an arrow is given. Further, as shown in FIG. 14, instead of the four radiation images S1 to S4, a human body image 66 is schematically displayed in the image display area 64, and an arrow is given to a place where a body motion has occurred in the human body image 66. You may do it.

  Further, instead of the text or together with the text, a guide image for prompting fixing of the part where the body movement has occurred may be displayed. FIG. 15 is a diagram showing a state in which a guide image is displayed. FIG. 15 shows a state in which a guide image is displayed together with text. As shown in FIG. 15, the image display area 64 is provided with an arrow between the radiation image S3 and the radiation image S4. The space between the radiation image S3 and the radiation image S4 corresponds to the knee position of the subject N. For this reason, the text 65 has the content “Please fix the knee with a belt”, and the guide image 67 is in a state where the knee is fixed with the belt, and the belt 68 included in the guide image 67 reciprocates in the direction of arrow A. The moving state is indicated by an animation.

  If the re-shooting assistance information is information that recommends one-shot shooting using a large FPD or CR long cassette, etc., the text should be "Please take one shot" Should be.

  As described above, in this embodiment, since the re-shooting auxiliary information is generated when it is determined that there is a body movement, the operator does not move by performing the shooting according to the re-shooting instruction information. So that re-shooting can be performed. Therefore, re-photographing can be performed efficiently, and the exposure amount to the subject can be reduced. In particular, by displaying the re-shooting auxiliary information, the operator can confirm the re-shooting auxiliary information at a glance, so that the re-shooting can be performed more efficiently.

  Next, a second embodiment of the present invention will be described. FIG. 16 is a schematic diagram showing a configuration of a radiographic image capturing apparatus to which the body movement detecting apparatus according to the second embodiment of the present invention is applied. In the second embodiment, the same components as those in the first embodiment are denoted by the same reference numerals, and detailed description thereof is omitted here. The radiographic image capturing apparatus 150A according to the second embodiment is different from the first embodiment in that a re-imaging control unit 48 is further provided.

  Based on the re-shooting auxiliary information generated by the re-shooting auxiliary information generation unit 46, the re-shooting control unit 48 instructs the above-described shooting adjustment unit to perform frame division again so as to avoid a place where body movement has occurred. . In the radiographic imaging apparatus 150A, a fixing tool that fixes the subject N with a belt or the like is used, and when the fixing tool can set the tightening strength, the fixing tool is strongly tightened. It may be an instruction to perform. Hereinafter, this instruction is referred to as re-photographing control.

  Next, processing performed in the second embodiment of the present invention will be described. FIG. 17 is a flowchart showing processing performed in the second embodiment. First, long imaging is performed while the FPD 110 is moved, and a radiation image at each position of movement is acquired (step ST21). Then, the local movement vector calculation unit 34 calculates a local movement vector V0 in an overlapping region of two adjacent radiographic images (step ST22), and the body movement index value calculation unit 36 calculates a body movement index value (step ST23). ). Subsequently, the body movement determination unit 38 determines the presence or absence of body movement (step ST24).

  When there is no body movement, the body movement determination unit 38 outputs the determination result together with the body movement index value to the body movement correction unit 40 and the image display unit 60 (step ST25). The body motion correction unit 40 corrects the body motion of the plurality of radiographic images based on the body motion index value (step ST26), and the image synthesis unit 42 synthesizes the radiographic images after the body motion correction to generate a composite image C1. (Step ST27). Then, the image display unit 60 displays the composite image C1 together with the body motion index value (step ST28), and the process ends.

  On the other hand, when the determination in step ST24 is that there is a body movement, the body movement determination unit 38 outputs the determination result to the warning unit 44 and the re-shooting auxiliary information generation unit 46 (step ST29). The re-shooting auxiliary information generation unit 46 generates re-shooting auxiliary information (step ST30) and outputs it to the image display unit 60 and the re-shooting control unit 48. The image display unit 60 displays the re-shooting information (step ST31). Further, the re-imaging control unit 48 performs the above-described re-imaging control (step ST32). On the other hand, the warning unit 44 issues a warning by a voice message (including a warning voice) or a buzzer sound (warning sound) (step ST33), and ends the process.

  In the second embodiment, the reshooting assistance information is displayed on the image display unit 60. However, only the reshooting control may be performed and the reshooting assistance information may not be displayed.

  Further, in the first and second embodiments, two adjacent FPDs 110 are moved for each position of the second and subsequent movements, that is, for each photographing after the second photographing. Calculation of body motion index value and determination of presence / absence of body motion based on the overlapping area of the radiographic image, a warning is generated when it is determined that there is body motion, and re-shooting auxiliary information is generated Good.

  In this case, when the operator confirms the body movement of the subject N by a warning or the like during the long shooting, the operator can stop shooting halfway by operating the emergency stop switch provided on the console 70. Thereby, it is possible to prevent the subject N from being unnecessarily exposed by continuing shooting even though there is a body movement amount that cannot be combined. In addition, the re-shooting auxiliary information can take measures so that body movement does not occur during re-shooting.

  In addition, a warning is generated when it is determined that the subject is moving during long shooting, but the subsequent shooting may be automatically stopped when the warning is generated. This eliminates the need for the emergency stop switch to be operated immediately after the operator confirms the body movement of the subject by a warning or the like, or the next shooting is performed before the emergency stop switch is operated. As a result, it is possible to prevent the subject from being unnecessarily exposed without burdening the operator.

  In the first and second embodiments, the body movement index value is detected by the local movement vector calculation unit 34 and the body movement index value calculation unit 36, but a sensor attached to the detector movement unit 20. The body movement may be detected by. As the sensor, for example, a sensor for detecting shaking and twisting provided on a gripping bar (not shown), a weight scale provided at a position where the left and right feet are respectively mounted, a sensor for detecting weight shift, or a subject A position sensor or the like attached to N itself can be considered.

  As described above, the apparatus 150 according to the embodiment of the present invention has been described. However, the computer includes the above-described imaging information acquisition unit 32, local movement vector calculation unit 34, body motion index value calculation unit 36, body motion determination unit 38, and re-imaging. A program that functions as means corresponding to the auxiliary information generation unit 46 and the re-imaging control unit 48 and performs processing as shown in FIGS. 9 and 17 is also one embodiment of the present invention. A computer-readable recording medium in which such a program is recorded is also one embodiment of the present invention.

DESCRIPTION OF SYMBOLS 30 Body motion detection part 31 Image acquisition part 32 Shooting information acquisition part 34 Local movement vector calculation part 36 Body motion index value calculation part 38 Body motion discrimination | determination part 40 Body motion correction part 42 Image composition part 44 Warning part 46 Re-shooting auxiliary information generation | occurrence | production Unit 48 re-imaging control unit 60 image display unit 100 radiation source 102 imaging surface 104 radiation 110 FPD
150,150A Radiation imaging device

Claims (11)

A radiographic imaging device that moves a radiation detector and irradiates the radiation detector with radiation that has passed through a subject each time the position is changed by the movement, and acquires a plurality of radiation images at least partially overlapping each other. There,
An imaging means for moving the radiation detector along a predetermined movement axis and irradiating the radiation detector with radiation transmitted through the subject each time the position is changed by the movement;
A radiological image acquisition means for acquiring a plurality of radiographic images of the subject by reading a signal from the radiation detector each time the movement and the irradiation of the radiation are performed;
Body motion detection means for detecting body motion of the subject at the time of capturing the plurality of radiation images;
An apparatus for generating a radiographic image, comprising: information generating means for generating re-imaging auxiliary information for assisting in re-imaging of the plurality of radiographic images when the body movement is detected.   The radiographic image capturing apparatus according to claim 1, further comprising display means for displaying the re-imaging auxiliary information.   The radiographic image capturing apparatus according to claim 1, further comprising a reimaging control unit configured to control the reimaging based on the reimaging auxiliary information.   The radiographic image capturing apparatus according to claim 3, wherein the re-imaging control unit is a unit that divides the plurality of radiographic images into frames while avoiding a position where the body movement of the subject occurs.   5. The radiographic image capturing apparatus according to claim 3, wherein the re-imaging control unit is a unit that increases a fixing strength of the subject by a fixture. The body motion detecting means; a local movement vector calculating means for calculating a local movement vector representing a local shift of the subject in an overlapping region in the plurality of radiographic images;
Body motion index value calculating means for calculating the body motion index value based on the local movement vector;
The radiographic imaging apparatus according to claim 1, further comprising body movement determination means for determining presence or absence of the body movement based on the body movement index value.   The radiographic image capturing apparatus according to claim 6, wherein the body motion index value calculating unit is a unit that calculates an index value of the parallel movement amount of the subject as the body motion index value.   The body movement index value calculating means is a means for calculating at least one of an index value of the three-dimensional movement amount of the subject and an index value of the two-dimensional movement amount of the subject as the body movement index value. The radiographic imaging apparatus according to claim 7.   The radiographic image capturing apparatus according to claim 1, wherein the body movement detection unit is a sensor that detects a movement of the subject. A radiation image capturing method for moving a radiation detector and irradiating the radiation detector with radiation that has passed through a subject each time the position is changed by the movement, thereby acquiring a plurality of radiation images at least partially overlapping each other. There,
Moving the radiation detector along a predetermined movement axis, and irradiating the radiation detector with radiation transmitted through the subject each time the position is changed by the movement;
Acquiring a plurality of radiation images of the subject by reading a signal from the radiation detector each time the movement and irradiation of the radiation are performed;
Detecting body movement of the subject at the time of capturing the plurality of radiation images;
And a step of generating re-imaging support information for assisting re-imaging of the plurality of radiographic images when the body motion is detected. A radiation image capturing method for moving a radiation detector and irradiating the radiation detector with radiation that has passed through a subject each time the position is changed by the movement, and acquiring a plurality of radiation images at least partially overlapping each other. A radiographic imaging program for causing a computer to execute,
A procedure of moving the radiation detector along a predetermined movement axis and irradiating the radiation detector with radiation transmitted through the subject each time the position is changed by the movement;
A procedure for acquiring a plurality of radiation images of the subject by reading a signal from the radiation detector each time the movement and the irradiation of the radiation are performed,
Detecting a body movement of the subject at the time of capturing the plurality of radiographic images;
A radiographic imaging program that causes a computer to execute a procedure of generating re-imaging assistance information for assisting re-imaging of the plurality of radiographic images when the body movement is detected.

Images