JP2012075862A - Body motion detector, method, and program - Google Patents

Body motion detector, method, and program Download PDF

Info

Publication number
JP2012075862A
JP2012075862A JP2011178725A JP2011178725A JP2012075862A JP 2012075862 A JP2012075862 A JP 2012075862A JP 2011178725 A JP2011178725 A JP 2011178725A JP 2011178725 A JP2011178725 A JP 2011178725A JP 2012075862 A JP2012075862 A JP 2012075862A
Authority
JP
Japan
Prior art keywords
body motion
movement
index value
subject
imaging
Prior art date
Legal status (The legal status 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 status listed.)
Abandoned
Application number
JP2011178725A
Other languages
Japanese (ja)
Other versions
JP2012075862A5 (en
Inventor
Yoshitaka Yamaguchi
義隆 山口
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Corp
Original Assignee
Fujifilm Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujifilm Corp filed Critical Fujifilm Corp
Priority to JP2011178725A priority Critical patent/JP2012075862A/en
Priority to EP11180443A priority patent/EP2428164A3/en
Priority to EP13173587.0A priority patent/EP2653105B1/en
Priority to EP13173594.6A priority patent/EP2659837A1/en
Priority to US13/227,184 priority patent/US20120059239A1/en
Priority to EP13198106.0A priority patent/EP2727533A1/en
Priority to CN201410575508.4A priority patent/CN104306009A/en
Priority to CN201110264400XA priority patent/CN102429672A/en
Publication of JP2012075862A publication Critical patent/JP2012075862A/en
Priority to US14/014,965 priority patent/US20140064450A1/en
Priority to US14/140,897 priority patent/US20140140481A1/en
Publication of JP2012075862A5 publication Critical patent/JP2012075862A5/ja
Abandoned legal-status Critical Current

Links

Images

Abstract

PROBLEM TO BE SOLVED: To achieve accurate detection of a body motion of a subject in the case where the subject may possibly move between imaging operations, such as during long-length imaging.SOLUTION: A local motion vector calculation unit 34 calculates local motion vectors in an overlapping area between each two adjacent radiographic images. A body motion index value calculation unit 36 calculates a body motion index value of the parallel translation, the three-dimensional movement, and the two-dimensional movement of a subject with using the local motion vectors. Further, a post treatment choice section 39, based on these body motion index values, performs the choice of the presence or absence of the enforcement of body motion correction, the choice of the body motion correction method, and the choice of the body motion index value displayed on an image display section 60.

Description

本発明は、同一被写体に関して複数回の撮影を行うことにより得られた複数の放射線画像を用いて、撮影時における被写体の体動を検出する体動検出装置、方法およびプログラムに関するものである。   The present invention relates to a body motion detection apparatus, method, and program for detecting body motion of a subject at the time of photographing using a plurality of radiographic images obtained by photographing the same subject a plurality of times.

従来、医療分野等において、被写体を透過した放射線の照射により被写体に関する放射線画像を記録する放射線検出器（いわゆる「Flat Panel Detector」 以下、ＦＰＤとする）が各種提案、実用化されている。このようなＦＰＤとしては、例えば、放射線の照射により電荷を発生するアモルファスセレン等の半導体を利用したＦＰＤがあり、そのようなＦＰＤとして、いわゆる光読取方式のものやＴＦＴ読取方式のものが提案されている。   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.

なお、体動は、被写体を構成する物質によって透過した放射線の減衰量が違うことを利用して、エネルギーが異なる２種類のエネルギーの放射線を被写体に照射して得られた２枚の放射線画像を用いたエネルギーサブトラクション処理を行うためのエネルギーサブトラクション撮影、および撮影時期が異なる２つの時系列的な放射線画像間の差異を表す差分画像を取得する、経時サブトラクション撮影を行う場合にも生じる。また、患部をより詳しく観察するために、Ｘ線管を移動させて異なる角度から被写体にＸ線を照射して撮影を行い、これにより取得した画像を加算して所望の断層面を強調した画像を得るトモシンセシス撮影、あるいは連続して被写体に放射線を照射して複数の画像を取得する連写撮影時にも生じる。   Note that body movement uses two radiation images obtained by irradiating a subject with radiations of two different energies using the fact that the attenuation amount of the transmitted radiation is different depending on the substance constituting the subject. This also occurs when performing energy subtraction imaging for performing the energy subtraction processing used and temporal subtraction imaging for acquiring a difference image representing a difference between two time-series radiographic images having different imaging timings. Also, in order to observe the affected area in more detail, the X-ray tube is moved to irradiate the subject with X-rays from different angles, and an image obtained by adding the acquired images to emphasize the desired tomographic plane It also occurs during tomosynthesis imaging that obtains image quality, or continuous shooting that continuously irradiates a subject with radiation and acquires a plurality of images.

このため、被写体の体動を検出して撮影を停止したり、体動を警告したりする手法が提案されている。例えば特許文献１には、センサを用いて撮影時の被写体の体動を検出し、体動が検出された場合に警告を行う手法が提案されている。また、特許文献２には、ＦＰＤの組み立て誤差および放射線画像撮影装置の設置誤差を考慮して、被写体の体動量を算出する手法が提案されている。具体的には、複数の放射線画像の重複領域においてテンプレートマッチングを行うことにより、局所的な複数のズレ量を求め、ズレ量の平均値を体動量として算出し、これをしきい値と比較して体動の有無を検出する手法が提案されている。なお、特許文献１，２に記載の手法は、いずれもＦＰＤの撮像面に対する被写体の平行移動を体動として検出している。さらに、特許文献１，２には、非線形ワーピングにより体動補正を行う手法も提案されている。   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. Note that the methods described in Patent Documents 1 and 2 both detect the parallel movement of the subject relative to the imaging surface of the FPD as a body movement. 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

しかしながら、体動には、特許文献１，２に記載されたような平行移動のみならず、被写体のねじれおよびＦＰＤの撮像面に対する前後方向の倒れという３次元的な動きも含む。さらには、撮像面に平行な面における被写体の回転および前後方向への平行移動による拡大縮小という２次元的な動きも含む。特許文献１，２に記載された手法は、体動を撮像面に平行な方向のズレ量としてのみ算出しているため、被写体の３次元的な動きおよび２次元的な動きは、平行移動に近似されてしまう。このように被写体の３次元的な動きおよび２次元的な動きが平行移動に近似されてしまうと、平行移動以外の被写体の動きが大きくても体動が小さいものとして検出される、あるいは体動が検出されなくなってしまうという問題がある。   However, the body motion includes not only parallel movement as described in Patent Documents 1 and 2, but also three-dimensional movement such as twisting of the subject and tilting in the front-rear direction with respect to the imaging surface of the FPD. Furthermore, it also includes a two-dimensional movement such as rotation of the subject on a plane parallel to the imaging plane and enlargement / reduction by translation in the front-rear direction. Since the methods described in Patent Documents 1 and 2 calculate body movement only as a shift amount in a direction parallel to the imaging surface, the three-dimensional movement and two-dimensional movement of the subject are translated. It will be approximated. Thus, if the three-dimensional movement and two-dimensional movement of the subject are approximated to parallel movement, even if the movement of the subject other than the parallel movement is large, it is detected that the body movement is small, or the body movement Is no longer detected.

一方、体動が検出されない（すなわち体動が小さい）と、非線形ワーピングにて体動補正を行うこととなるが、上述したように平行移動以外の被写体の動きが大きくても体動が小さくなる場合があるため、３次元的な動きおよび２次元的な動きが大きくても体動補正が行われてしまうこととなる。このように被写体の動きが大きい場合に体動補正が行われると、過度に画像が補正されて歪みが大きくなるため、正確な診断に支障をきたすおそれがある。   On the other hand, if body movement is not detected (that is, body movement is small), body movement correction is performed by non-linear warping. However, as described above, body movement decreases even if the movement of the subject other than parallel movement is large. In some cases, body motion correction is performed even if the three-dimensional motion and the two-dimensional motion are large. Thus, if body motion correction is performed when the movement of the subject is large, the image is excessively corrected and distortion increases, which may hinder accurate diagnosis.

本発明は上記事情に鑑みなされたものであり、長尺撮影のように撮影毎に被写体が移動する可能性がある場合に、被写体の体動を精度良く検出できるようにすることを目的とする。   The present invention has been made in view of the above circumstances, and an object of the present invention is to enable accurate detection of body movement of a subject when there is a possibility that the subject moves for each photographing as in long photographing. .

本発明による体動検出装置は、同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得する画像取得手段と、
前記複数の放射線画像の撮影時における前記被写体の動きの種類に応じた複数種類の体動指標値を取得する体動指標値取得手段とを備えたことを特徴とするものである。 The body motion detection device according to the present invention includes an image acquisition unit that acquires a plurality of radiographic images obtained by performing imaging a plurality of times on the same subject and at least some of the regions overlap.
It is characterized by comprising body motion index value acquisition means for acquiring a plurality of types of body motion index values corresponding to the types of movement of the subject at the time of capturing the plurality of radiographic images.

「少なくとも一部が重複した」とは、一部が重複する場合のみならず、全部が重複する場合も含む。   “At least partially overlapped” includes not only partially overlapping but also completely overlapping.

なお、本発明による体動検出装置においては、前記体動指標値取得手段を、前記複数の放射線画像における重複領域内の局所的な前記被写体のズレを表す局所移動ベクトルを算出する局所移動ベクトル算出手段と、
前記局所移動ベクトルに基づいて前記複数種類の体動指標値を算出する体動指標値算出手段とを備えてなるものとしてもよい。 In the body motion detection device according to the present invention, the body motion index value acquisition unit calculates a local motion vector that calculates a local motion vector that represents a local shift of the subject in an overlapping region in the plurality of radiation images. Means,
Body movement index value calculating means for calculating the plurality of types of body movement index values based on the local movement vector may be provided.

この場合、前記体動指標値算出手段を、前記被写体の平行移動量の指標値、前記被写体の３次元的な移動量の指標値および該被写体の２次元的な移動量の指標値の少なくとも２つを前記体動指標値として算出する手段としてもよい。   In this case, the body motion index value calculating means is configured to perform at least two of an index value of the parallel movement amount of the subject, 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. It is good also as a means to calculate one as said body movement index value.

また、本発明による体動検出装置においては、前記複数種類の体動指標値または該複数種類の体動指標値のうちの主たる体動指標値に基づいて、前記複数の放射線画像に関連する後処理を選択する後処理選択手段をさらに備えるものとしてもよい。   In the body movement detection device according to the present invention, the plurality of types of body movement index values or the main body movement index values of the plurality of types of body movement index values may be Post-processing selection means for selecting a process may be further provided.

また、本発明による体動検出装置においては、前記後処理選択手段を、前記複数種類の体動指標値または前記主たる体動指標値に基づいて体動がある前記動きの種類を判定し、該判定結果に応じて前記後処理を選択する手段としてもよい。   Further, in the body motion detection device according to the present invention, the post-processing selection means determines the type of motion having body motion based on the plurality of types of body motion index values or the main body motion index values, The post-processing may be selected according to the determination result.

また、本発明による体動検出装置においては、前記後処理選択手段を、前記複数の放射線画像に体動補正処理を施すか否かを選択する手段としてもよい。   In the body movement detection device according to the present invention, the post-processing selection unit may be a unit that selects whether or not to perform body movement correction processing on the plurality of radiation images.

また、本発明による体動検出装置においては、前記後処理選択手段を、前記複数の放射線画像に施す体動補正処理の種類を選択する手段としてもよい。   In the body motion detection device according to the present invention, the post-processing selection unit may be a unit that selects a type of body motion correction processing to be performed on the plurality of radiation images.

また、本発明による体動検出装置においては、前記後処理選択手段を、表示する体動指標値の種類を選択する手段としてもよい。   In the body motion detection device according to the present invention, the post-processing selection means may be means for selecting the type of body motion index value to be displayed.

本発明による放射線画像撮影装置は、放射線検出器を移動させ、該移動により位置を変える毎に被写体を透過した放射線を該放射線検出器に照射して、少なくとも一部の領域が重複した複数の放射線画像を取得する放射線画像撮影装置であって、
前記放射線検出器を所定の移動軸に沿って移動させ、該移動により位置を変える毎に前記被写体を透過した放射線を前記放射線検出器に照射する撮影手段と、
前記移動および前記放射線の照射を行う毎に前記放射線検出器から信号を読み出すことにより、前記被写体の複数の放射線画像を取得する放射線画像取得手段と、
本発明による体動検出装置とを備えたことを特徴とするものである。 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;
The body motion detecting device according to the present invention is provided.

本発明による体動検出方法は、同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得し、
前記複数の放射線画像の撮影時における前記被写体の動きの種類に応じた複数種類の体動指標値を取得することを特徴とするものである。 The body motion detection method according to the present invention obtains a plurality of radiation images obtained by performing imaging a plurality of times on the same subject, at least a part of which overlaps,
A plurality of types of body motion index values corresponding to the types of movement of the subject at the time of capturing the plurality of radiation images are acquired.

なお、本発明による体動検出方法をコンピュータに実行させるためのプログラムとして提供してもよい。   In addition, you may provide as a program for making a computer perform the body movement detection method by this invention.

本発明の体動検出装置および方法によれば、複数の放射線画像の撮影時における被写体の体動に応じた複数種類の体動指標値が取得される。このため、被写体の体動に応じた適切な体動指標値を取得することができ、その結果、精度よく体動を検出することができる。また、検出した体動指標値に基づいて、後段で行われる複数の放射線画像に関連する後処理を適切に選択することができる。   According to the body motion detection apparatus and method of the present invention, a plurality of types of body motion index values corresponding to the body motion of the subject at the time of capturing a plurality of radiation images are acquired. For this reason, an appropriate body motion index value according to the body motion of the subject can be acquired, and as a result, the body motion can be detected with high accuracy. Further, based on the detected body motion index value, it is possible to appropriately select post-processing related to a plurality of radiation images performed in the subsequent stage.

また、複数の放射線画像における重複領域内の局所領域毎に、重複領域における被写体のズレを表す局所移動ベクトルを算出し、局所移動ベクトルに基づいて体動指標値を算出することにより、比較的簡易な演算により体動指標値を取得することができる。   In addition, for each local region in the overlapping region in a plurality of radiographic images, a local movement vector that represents the displacement of the subject in the overlapping region is calculated, and a body motion index value is calculated based on the local movement vector, so that it is relatively simple. The body motion index value can be acquired by simple calculation.

また、被写体の平行移動量の指標値、被写体の３次元的な移動量の指標値および被写体の２次元的な移動量の指標値の少なくとも２つを体動指標値として算出することにより、被写体の平行移動に起因する体動、被写体のねじれおよび倒れという３次元的な体動、並びに被写体の回転移動および拡大縮小という２次元的な体動の少なくとも２つを検出することができる。   Further, by calculating at least two of the index value of the parallel movement amount of the subject, 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 movement index value, Can detect at least two of the three-dimensional body movements caused by the parallel movement of the subject, the three-dimensional body movements of the subject twisting and tilting, and the two-dimensional body movements of the subject rotational movement and enlargement / reduction.

また、本発明による体動検出装置および方法を、長尺撮影に適用することにより、長尺撮影により取得される複数の放射線画像間に生じる被写体の体動を精度良く検出することができる。   In addition, by applying the body motion detection device and method according to the present invention to long imaging, it is possible to accurately detect body motion of a subject that occurs between a plurality of radiographic images acquired by long imaging.

本発明の第１の実施形態による体動検出装置を適用した放射線画像撮影装置の構成を示す概略図Schematic which shows the structure of the radiographic imaging apparatus to which the body movement detection apparatus by the 1st Embodiment of this invention is applied. 局所移動ベクトルの算出を説明するための図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 the display in the display screen of an image display part 画像表示部の表示画面における表示の他の例を示す図The figure which shows the other example of the display in the display screen of an image display part 本発明の第２の実施形態による体動検出装置を適用した放射線画像撮影装置の構成を示す概略図Schematic which shows the structure of the radiographic imaging apparatus to which the body movement detection apparatus by the 2nd Embodiment of this invention is applied. 第２の実施形態において行われる処理を示すフローチャートThe flowchart which shows the process performed in 2nd Embodiment. 複数の合成画像を表示した状態を示す図The figure which shows the state which displayed the plural synthetic images 複数の合成画像を表示した状態を示す図The figure which shows the state which displayed the plural synthetic images 第３の実施形態において行われる処理を示すフローチャートThe flowchart which shows the process performed in 3rd Embodiment

以下、図面を参照して本発明の実施形態について説明する。図１は本発明の第１の実施形態による体動検出装置を適用した放射線画像撮影装置の構成を示す概略図である。図１に示すように第１の実施形態による放射線画像撮影装置１５０は、１つの放射線源１００と１つのＦＰＤ１１０とを用いて、被写体Ｎ中の隣接する複数の領域Ｎ１、Ｎ２・・・を順次撮影し、これにより取得した複数の放射線画像を合成して被写体Ｎの大部分を示す長尺の放射線画像を得る、長尺撮影を行うことができるように構成されたものである。   Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a schematic diagram showing a configuration of a radiographic imaging apparatus to which a body motion detection apparatus according to a first embodiment of the present invention is applied. 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.

ＦＰＤ１１０は、被写体Ｎを透過した放射線１０４を検出して電気信号に変換し、その被写体Ｎの放射線画像を表す画像データを出力する。なお、ＦＰＤ１１０は、放射線を電荷に直接変換する直接方式のＦＰＤ、または放射線を一旦光に変換し、変換された光をさらに電気信号に変換する間接方式のＦＰＤのいずれも利用可能である。直接方式のＦＰＤは、アモルファスセレン等の光導電膜、キャパシタおよびスイッチ素子としてのＴＦＴ（Thin Film Transistor）等によって構成される。例えば、Ｘ線等の放射線が入射されると、光導電膜から電子−正孔対（ｅ−ｈペア）が発せられる。その電子−正孔対はキャパシタに蓄積され、キャパシタに蓄積された電荷が、ＴＦＴを介して電気信号として読み出される。   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 combined, and an image combining unit 42 that generates a long radiation image representing the entire subject N and a warning unit 44 are provided. 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 post-processing selection unit 39. The local movement vector calculation unit 34 and the body motion index value calculation unit 36 constitute a body motion index value acquisition unit.

画像取得部３１は、ＦＰＤ１１０から放射線画像を取得するための各種インターフェースからなる。なお、体動検出部３０を放射線画像撮影装置１５０と別個に設け、これらをネットワークにより接続した場合、画像取得部３１はネットワークインターフェースとなる。   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). Can do. 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 twice. 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. Accordingly, in the case where the resolution conversion is performed twice and the multi-resolution conversion is performed up to ¼ resolution, assuming that 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 is calculated in FIG. 3, the size of the overlapping region image is the same as that of the overlapping regions K1-1 and K2-1.

次いで、次に高い解像度の重複領域画像Ｋ１−２，Ｋ２−２を用いて、４つの局所移動ベクトルＶ０を算出する。この際、重複領域画像Ｋ１−３，Ｋ２−３において算出した局所移動ベクトルＶ０を用いることにより、効率よく局所移動ベクトルＶ０を算出できる。例えば、図３に示すように、重複領域画像Ｋ１−３，Ｋ２−３において右斜め上を向く局所移動ベクトルＶ０が算出された場合、次の解像度の重複領域画像Ｋ１−２，Ｋ２−２においては、前の解像度にて算出した局所移動ベクトルＶ０を初期値として、図４に示すように、探索範囲Ｒの右斜め上付近の領域においてのみ（斜線で示す）テンプレートマッチングを行うようにすれば、相関値算出のための演算量を低減できるため、局所移動ベクトルＶ０の算出を効率よく行うことができる。   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. When V0 is calculated and the imaging region is assumed to have many one-dimensional structures (for example, lower limbs, all lower limbs, etc.), the correlation value is likely to be shaken along the one-dimensional direction. The local movement vector V0 may be calculated based on the 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 calculation unit 34. 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. .

体動指標値算出部３６は、例えば「図形処理工学、山口富士夫、日刊工業新聞社刊、68-82頁、1981年」に記載された式に基づいて、被写体Ｎの２次元的な動きの体動指標値を算出する。具体的には、上述したように算出した局所移動ベクトルＶ０を２次元アフィン変換の式に適用し、複数の局所移動ベクトルＶ０から最小二乗法を用いて、平行移動、回転および拡大縮小の各要素を体動指標値として算出する。ここで、２次元アフィン変換において、横方向の移動量をｔｘ、縦方向の移動量をｔｙとすると、平行移動は下記の式（１）により表される。なお、式（１）および以降の式において、（ｘ，ｙ）は移動前の２次元座標位置、（ｘ*，ｙ*）は移動後の２次元座標位置である。

Figure 2012075862
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.
Figure 2012075862

また２次元アフィン変換において、回転角をθとすると、回転移動は下記の式（２）により表される。

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

また、２次元アフィン変換において、拡大縮小のための拡大縮小パラメータをａ（Ｘ方向），ｄ（Ｙ方向）とすると、拡大縮小は下記の式（３）により表される。

Figure 2012075862
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).
Figure 2012075862

式（１′）〜（３′）をまとめると、

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

となり、Ｈは９個の要素からなる行列となる。ここで、局所移動ベクトルＶ０は、（ｘ*−ｘ，ｙ*−ｙ）と表すことができるため、複数の局所移動ベクトルＶ０を上記の式（４）に適用し、最小二乗法により９個の要素を求めることにより、平行移動量、回転角度および拡大縮小パラメータを算出することができる。このようにして求めた平行移動量、回転角度および拡大縮小パラメータが２次元的な動きの体動指標値となる。なお、本実施形態においては、以降の処理において、平行移動の体動指標値は上述したようにヒストグラムに基づいて算出されたものを用いるものとし、２次元的な動きの体動指標値としては、回転角度および拡大縮小パラメータのみを用いるものとする。 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 two body motion index values 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 post-processing selection unit 39 selects a process to be performed in the subsequent stage based on a plurality of types of body motion index values corresponding to the movement of the subject calculated as described above. Therefore, the post-processing selection unit 39 first compares the body motion index value calculated by the body motion index value calculation unit 36 with a threshold value, and there is a body motion when the body motion index value is equal to or greater than the threshold value. Is determined. In the present embodiment, since a plurality of body motion index values for parallel movement, three-dimensional motion, and two-dimensional motion are calculated, the body motion index value for each of the body motion index values is set to a threshold value corresponding to each. Compared with the value, the presence or absence of body movement is individually determined. Further, for two-dimensional movement, threshold values are prepared for each of rotation and enlargement / reduction, and the presence / absence of body movement is determined for each of rotation and enlargement / reduction.

ここで、本実施形態においては、後処理選択部３９は、撮影情報取得部３２が算出した撮影情報に基づいて、体動指標値をしきい値と比較することにより体動の有無を判別する。具体的には、撮影情報に応じてしきい値の大きさを切り換える。例えば、撮影部位が下肢（長尺撮影でない）や全下肢（長尺撮影）の場合、心臓または腸にたまったガスのように被写体Ｎの体動とは別に動く構造物が含まれない。これに対して、撮影部位が胸部（長尺撮影でない）や全脊椎（長尺撮影）の場合、心臓または腸にたまったガスのように被写体Ｎの体動とは別に動く構造物が含まれる。このため、撮影情報が、撮影部位が胸部である旨の情報を含む場合には、撮影部位が下肢の場合よりも判別のためのしきい値を大きくする。これにより、局所的な動きに影響されることなく、体動の有無を判別できる。   Here, in the present embodiment, the post-processing selection unit 39 determines the presence / 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.

さらに、後処理選択部３９は、体動補正部４０において体動補正を行うか否かを選択する。具体的には、３次元的な動きの体動があると判別した場合には、平行移動および２次元的な動きの体動の有無に拘わらず体動補正を行わない。また、２次元的な動きの体動があると判別した場合にも、平行移動および３次元的な動きの体動の有無に拘わらず体動補正を行わない。なお、平行移動の体動については、３次元的な動きの体動および２次元的な動きの体動がないと判別された場合にのみ、体動の有無に拘わらず体動補正を行う。すなわち、３次元的な動きおよび２次元的な動きのいずれかの体動があると判別された場合以外は体動補正を行う。   Further, the post-processing selection unit 39 selects whether or not the body motion correction unit 40 performs body motion correction. Specifically, when it is determined that there is a three-dimensional body movement, no body movement correction is performed regardless of the presence of parallel movement and two-dimensional body movement. Even when it is determined that there is a two-dimensional movement, the body movement correction is not performed regardless of whether there is a parallel movement or a three-dimensional movement. For translational body motion, body motion correction is performed regardless of the presence or absence of body motion only when it is determined that there is no three-dimensional motion body motion or two-dimensional motion body motion. That is, body motion correction is performed except when it is determined that there is any one of three-dimensional motion and two-dimensional motion.

また、後処理選択部３９は、体動の有無の判別結果および体動指標値に基づいて、体動補正を行うと選択した場合の体動補正の方法を選択する。ここで、本実施形態においては、体動指標値をしきい値と比較し、体動指標値がしきい値未満となった場合に体動無しと判別しているため、体動無しと判別された場合でも、体動指標値は０でないことから体動は生じているものである。このため、３次元的な動きの体動が生じている場合には、非線形ワーピングを選択する。また、２次元的な動きにおける回転の体動が生じている場合には、上記式（２）を用いて回転により体動補正を行う方法を選択する。また、２次元的な動きにおける拡大縮小の体動が生じている場合には、上記式（３）を用いて拡大縮小により体動補正を行う方法を選択する。なお、２次元的な動きの体動が生じている場合には、非線形ワーピングを行う方法を選択してもよい。また、３次元的な動きおよび２次元的な動きの体動の双方が生じている場合には、非線形ワーピングを行う方法を選択する。なお、これ以外の場合には、局所移動ベクトルＶ０の方向に基づいて、平行移動の体動補正を行う方法を選択する。   Further, the post-processing selection unit 39 selects a body motion correction method when it is selected that body motion correction is to be performed based on the determination result of the presence or absence of body motion and the body motion index value. Here, in the present embodiment, the body motion index value is compared with a threshold value, and when the body motion index value falls below the threshold value, it is determined that there is no body motion. Even in such a case, the body motion index value is not 0, so that the body motion has occurred. For this reason, when a three-dimensional movement is generated, nonlinear warping is selected. In addition, when a rotational body motion occurs in a two-dimensional motion, a method of correcting body motion by rotation is selected using the above equation (2). In addition, when enlargement / reduction body movement in a two-dimensional movement occurs, a method of performing body movement correction by enlargement / reduction is selected using the above equation (3). Note that when a two-dimensional movement is occurring, a method of performing nonlinear warping may be selected. When both three-dimensional motion and two-dimensional motion are occurring, a method for performing nonlinear warping is selected. In other cases, a method of correcting the body movement of the parallel movement is selected based on the direction of the local movement vector V0.

また、後処理選択部３９は、体動の有無の判別結果および体動指標値に基づいて、後述する画像表示部６０に表示する体動指標値を選択する。具体的には、３次元的な動きの体動があると判別した場合あるいは３次元的な動きの体動がないと判別したが、３次元的な動きの体動が生じている場合には、３次元的な動きの体動指標値（すなわち標準偏差または分散）を表示することを選択する。また、２次元的な動きにおける回転の体動があると判別した場合あるいは回転の体動がないと判別したが回転の体動が生じている場合には、回転の体動指標値（すなわち回転角θ）を表示することを選択する。また、２次元的な動きにおける拡大縮小の体動があると判別した場合あるいは拡大縮小の体動がないと判別したが拡大縮小の体動が生じている場合には、拡大縮小の体動指標値（すなわちパラメータａ，ｄ）を表示することを選択する。なお、平行移動の体動については、平行移動の体動指標値を表示することを選択する。   Further, the post-processing selection unit 39 selects a body motion index value to be displayed on the image display unit 60 described later based on the determination result of the presence or absence of body motion and the body motion index value. Specifically, if it is determined that there is a three-dimensional movement or if there is no three-dimensional movement, but there is a three-dimensional movement, It is selected to display a body motion index value (that is, standard deviation or variance) of a three-dimensional motion. When it is determined that there is a rotational body motion in a two-dimensional motion, or when it is determined that there is no rotational body motion, but there is a rotational body motion, a rotational body motion index value (that is, a rotational body motion) Select to display the angle θ). Also, when it is determined that there is an enlargement / reduction body movement in a two-dimensional movement or when there is no enlargement / reduction body movement, but an enlargement / reduction body movement occurs, an enlargement / reduction body movement index Choose to display values (ie parameters a, d). For translational body motion, it is selected to display a translational body motion index value.

次に、体動補正部４０による体動の補正について説明する。体動補正部４０は、後処理選択部３９が体動補正を行うと選択した場合において、さらに後処理選択部３９が選択した体動補正の方法を用いて、被写体Ｎの体動に起因する画像歪みを解消するための補正を放射線画像に対して行う。ここで、平行移動の体動補正が選択された場合について説明すると、図８（ａ）に示すように、平行移動の体動指標値が検出されていることとなるため、図８（ｂ）に示すように、体動指標値すなわち局所移動ベクトルＶ０に基づいて、２つの放射線画像Ｓ１，Ｓ２を相対的にずらすことにより体動補正を行う。また、非線形ワーピングが選択された場合には、複数の放射線画像を相対的に非線形ワーピングすることにより体動補正を行う。また、回転または拡大縮小の体動補正が選択されている場合には、上記式（２）または式（３）を用いて体動補正を行う。   Next, correction of body movement by the body movement correction unit 40 will be described. When the post-processing selecting unit 39 selects to perform body motion correction, the body motion correcting unit 40 is caused by the body motion of the subject N using the body motion correcting method selected by the post-processing selecting unit 39. Correction for eliminating image distortion is performed on the radiation image. Here, the case where the translational movement correction is selected will be described. Since the translational movement index value is detected as shown in FIG. 8A, FIG. 8B. As shown in FIG. 4, 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. When nonlinear warping is selected, body motion correction is performed by relatively nonlinearly warping a plurality of radiation images. Further, when the body motion correction of rotation or enlargement / reduction is selected, the body motion correction is performed using the above formula (2) or formula (3).

画像合成部４２は、体動補正済みの放射線画像をつなぎ合わせるように合成して合成画像Ｃ１を生成する。   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 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 post-processing selection unit 39 determines the presence or absence of body movement based on the body movement index value (step ST4). Then, the post-processing selection unit 39 selects the body motion index value to be displayed on the image display unit 60 based on the determination result of the presence / absence of body motion and the body motion index value (step ST5), and determines the presence / absence of body motion. It is selected whether or not to perform body motion correction according to the result (step ST6). When step ST6 is affirmed, the post-processing selection unit 39 selects a body motion correction method (step ST7), and stores the body motion index value used for the body motion correction and information on the selected body motion correction method. When displayed on the motion correction unit 40, the selected body motion index value is output to the image display unit 60 (information output, step ST8).

次いで、後処理選択部３９が選択した体動補正の方法を用いて、体動補正部４０が複数の放射線画像の体動を補正し（ステップＳＴ９）、画像合成部４２が体動補正済みの放射線画像を合成して合成画像Ｃ１を生成する（ステップＳＴ１０）。そして、画像表示部６０が合成画像Ｃ１を、表示すると選択された体動指標値とともに表示し（ステップＳＴ１１）、処理を終了する。   Next, using the body motion correction method selected by the post-processing selection unit 39, the body motion correction unit 40 corrects the body motion of the plurality of radiographic images (step ST9), and the image composition unit 42 has corrected the body motion. The radiation images are combined to generate a combined image C1 (step ST10). Then, the image display unit 60 displays the composite image C1 together with the selected body motion index value when displayed (step ST11), 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. In the body motion index value display area 63, only the body motion index value selected when displayed is displayed. FIG. 10 shows a state in which body movement index values for parallel movement, three-dimensional movement, and two-dimensional movement rotation and enlargement / reduction are 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, if step ST6 is negative, the post-processing selection unit 39 outputs a selection result indicating that the body motion correction is not performed to the warning unit 44, and outputs the selected body motion index value to the image display unit 60 when the result is displayed ( Information output, step ST12). The warning unit 44 issues a warning by a voice message (including a warning voice) or a buzzer sound (warning sound) (step ST13). Further, when the image display unit 60 displays, the selected body motion index value is displayed (step ST14), and the process ends. 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”. Further, a warning mark or a warning message may be displayed on the display screen of the image display unit 60. Moreover, you may make it generate a warning by both a sound and a display.

このように、本実施形態においては、複数の放射線画像の撮影時における被写体の体動に応じた複数種類の体動指標値を取得するようにしたため、平行移動、３次元的な動きおよび２次元的な動きという、被写体の体動に応じた適切な体動指標値を取得することができ、その結果、精度よく体動を検出することができる。また、検出した体動指標値に基づいて、後段で行われる体動補正実行の有無の選択、体動補正方法の選択および表示する体動指標値を選択する処理を適切に行うことができる。   As described above, in the present embodiment, since a plurality of types of body motion index values corresponding to the body motion of the subject at the time of capturing a plurality of radiographic images are acquired, parallel movement, three-dimensional motion, and two-dimensional motion are obtained. It is possible to acquire an appropriate body motion index value corresponding to the body motion of the subject, and to detect the body motion with high accuracy. In addition, based on the detected body motion index value, it is possible to appropriately perform processing for selecting whether or not to perform body motion correction performed in a later stage, selecting a body motion correction method, and selecting a body motion index value to be displayed.

次いで、本発明の第２の実施形態について説明する。図１２は本発明の第２の実施形態による体動検出装置を適用した放射線画像撮影装置の構成を示す概略図である。なお、第２の実施形態において第１の実施形態と同一の構成については同一の参照番号を付与し、ここでは詳細な説明は省略する。第２の実施形態による放射線画像撮影装置１５０Ａは、後処理選択部３９Ａが、平行移動、３次元の動きおよび２次元の動きの体動指標値のうちの、最も大きい体動指標値を主たる体動指標値に決定し、主たる体動指標値に基づいて体動補正の実施の選択、体動補正方法の選択および表示する体動指標値の選択を行うようにした点が第１の実施形態と異なる。   Next, a second embodiment of the present invention will be described. FIG. 12 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. In the radiographic image capturing apparatus 150A according to the second embodiment, the post-processing selection unit 39A has the largest body motion index value among the body motion index values of the parallel movement, the three-dimensional motion, and the two-dimensional motion. The first embodiment is that the movement index value is determined, and the selection of the body movement correction, the selection of the body movement correction method, and the selection of the body movement index value to be displayed are performed based on the main body movement index value. And different.

すなわち、第２の実施形態においては、後処理選択部３９Ａは、体動指標値算出部３６が算出した平行移動、３次元の動きおよび２次元の動きの体動指標値の大きさを判定する。この場合、平行移動、３次元の動きおよび２次元の動きの体動指標値は単位が異なるため、単位を正規化して体動指標値の大きさを判定すればよい。そして最も大きい体動指標値を主たる体動指標値に決定し、主たる体動指標値に基づいて、体動補正の実施の選択を行う。すなわち、主たる体動指標値が３次元的な動きまたは２次元的な動きの体動指標値である場合、その体動指標値をしきい値と比較して、体動補正を行うか否かを選択する。なお、主たる体動指標値が平行移動である場合は、３次元的な動きの体動および２次元的な動きの体動がないと判別された場合に、体動補正を行うことを選択する。   That is, in the second embodiment, the post-processing selection unit 39A determines the magnitude of the body motion index value of the parallel movement, the three-dimensional motion, and the two-dimensional motion calculated by the body motion index value calculation unit 36. . 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. Then, the largest body motion index value is determined as the main body motion index value, and selection of the body motion correction is performed based on the main body motion index value. That is, when the main body motion index value is a three-dimensional motion or a two-dimensional motion body motion index value, whether or not to perform body motion correction by comparing the body motion index value with a threshold value. Select. When the main body motion index value is parallel movement, it is selected to perform body motion correction when it is determined that there is no body motion of 3D motion and 2D motion. .

また、体動補正を実施すると選択した場合には、体動補正の方法を選択する。すなわち、主たる体動指標値が３次元的な動きの体動指標値である場合には、非線形ワーピングを選択する。また、主たる体動指標値が２次元的な動きにおける回転の体動指標値である場合には、上記式（２）を用いて回転により体動補正を行う方法を選択する。また、主たる体動指標値が２次元的な動きにおける拡大縮小の体動指標値である場合には、上記式（３）を用いて拡大縮小により体動補正を行う方法を選択する。なお、２次元的な動きの体動が生じている場合には、非線形ワーピングを行う方法を選択してもよい。さらに、主たる体動指標値が平行移動の体動指標値である場合には、局所移動ベクトルＶ０の方向に基づいて、平行移動の体動補正を行う方法を選択する。   If it is selected that body motion correction is to be performed, a method for body motion correction is selected. That is, when the main body motion index value is a three-dimensional motion body motion index value, nonlinear warping is selected. In addition, when the main body motion index value is a rotational body motion index value in a two-dimensional motion, a method of correcting body motion by rotation is selected using the above equation (2). Further, when the main body motion index value is a body motion index value of enlargement / reduction in a two-dimensional motion, a method of performing body motion correction by enlargement / reduction is selected using the above formula (3). Note that when a two-dimensional movement is occurring, a method of performing nonlinear warping may be selected. Further, when the main body motion index value is a translational body motion index value, a method of performing translational body motion correction is selected based on the direction of the local motion vector V0.

次いで、本発明の第２の実施形態において行われる処理について説明する。図１３は第２の実施形態において行われる処理を示すフローチャートである。まず、ＦＰＤ１１０を移動させつつ長尺撮影を行い、移動の各位置における放射線画像を取得する（ステップＳＴ２１）。そして、局所移動ベクトル算出部３４が、隣接する２つの放射線画像の重複領域における局所移動ベクトルＶ０を算出し（ステップＳＴ２２）、体動指標値算出部３６が体動指標値を算出する（ステップＳＴ２３）。   Next, processing performed in the second embodiment of the present invention will be described. FIG. 13 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 post-processing selection unit 39 determines the main body movement index value by comparing the magnitudes of the body movement index values (step ST24). Furthermore, the post-processing selection unit 39 determines the presence or absence of body movement based on the main body movement index value (step ST25). Then, the post-processing selection unit 39 selects the main body motion index value as the body motion index value to be displayed on the image display unit 60 (step ST26), and according to the determination result of the presence or absence of the body motion index value. It is selected whether or not to perform body movement correction (step ST27). If step ST27 is positive, the post-processing selection unit 39 selects a body motion correction method based on the main body motion index value (step ST28), and determines the main body motion index value and the selected body motion correction method. When the information is displayed on the body motion correction unit 40, the selected main body motion index value is output to the image display unit 60 (information output, step ST29).

次いで、後処理選択部３９が選択した体動補正の方法を用いて、体動補正部４０が複数の放射線画像の体動を補正し（ステップＳＴ３０）、画像合成部４２が体動補正済みの放射線画像を合成して合成画像Ｃ１を生成する（ステップＳＴ３１）。そして、画像表示部６０が合成画像Ｃ１を、表示すると選択された体動指標値とともに表示し（ステップＳＴ３２）、処理を終了する。   Next, using the body movement correction method selected by the post-processing selection unit 39, the body movement correction unit 40 corrects the body movements of the plurality of radiographic images (step ST30), and the image composition unit 42 has corrected the body movements. The radiation images are combined to generate a combined image C1 (step ST31). Then, the image display unit 60 displays the composite image C1 together with the body motion index value selected when displayed (step ST32), and the process is terminated.

一方、ステップＳＴ２７が否定されると、後処理選択部３９は、体動補正を行わない旨の選択結果を警告部４４に、表示すると選択した主たる体動指標値を画像表示部６０に出力する（情報出力、ステップＳＴ３３）。警告部４４は、音声メッセージ（警告音声を含む）やブザー音（警告音）により警告を行う（ステップＳＴ３４）。さらに画像表示部６０が表示すると選択された体動指標値を表示し（ステップＳＴ３５）、処理を終了する。   On the other hand, if step ST27 is negative, the post-processing selection unit 39 outputs the selection result indicating that the body motion correction is not performed to the warning unit 44, and outputs the selected main body motion index value to the image display unit 60. (Information output, step ST33). The warning unit 44 issues a warning by a voice message (including a warning voice) or a buzzer sound (warning sound) (step ST34). Further, when the image display unit 60 displays, the selected body motion index value is displayed (step ST35), and the process ends.

なお、上記第１の実施形態においては、後処理選択部３９が複数種類の体動補正方法を選択する場合がある。例えば、平行移動、３次元的な動きおよび２次元的な動きのすべてが生じている場合には、平行移動、非線形ワーピング、回転および拡大縮小のすべての体動補正方法が選択される。なお、２次元的な動きについては、非線形ワーピングが選択される場合もある。このため、体動補正部４０においては、後処理選択部３９が選択した複数の体動補正方法について、取り得るすべての組み合わせまたは任意の組み合わせの体動補正方法を用いて体動補正を行うようにしてもよい。例えば、平行移動、非線形ワーピング、回転および拡大縮小のすべての体動補正方法が選択されている場合は、体動補正の方法は４種類あるため、４種類すべての体動補正方法（１通り）、４種類のうちの３種類の体動補正方法（４通り）、４種類のうちの２種類の体動補正方法（６通り）、４種類のうちの１種類の体動補正方法（４通り）の合計１５通りの体動補正方法により体動補正を行うことが可能である。   In the first embodiment, the post-processing selection unit 39 may select a plurality of types of body motion correction methods. For example, when all of translation, three-dimensional movement, and two-dimensional movement have occurred, all body motion correction methods of translation, nonlinear warping, rotation, and scaling are selected. Note that non-linear warping may be selected for two-dimensional movement. Therefore, the body motion correction unit 40 performs body motion correction using all possible combinations or any combination of body motion correction methods for the plurality of body motion correction methods selected by the post-processing selection unit 39. It may be. For example, when all body motion correction methods of translation, nonlinear warping, rotation and enlargement / reduction are selected, there are four types of body motion correction methods, so all four types of body motion correction methods (one way) 4 types of body motion correction methods (4 types), 4 types of body motion correction methods (6 types), 4 types of body motion correction methods (4 types) ) Body motion correction can be performed by a total of 15 body motion correction methods.

この場合、体動補正部４０は、１５通りのすべての体動補正方法またはこれらのうちの任意の数の体動補正方法を用いて体動補正を行う。画像合成部４２は行われたすべての体動補正方法により補正された放射線画像を用いて複数の合成画像を生成する。そして、画像表示部６０は、複数の合成画像を表示する。図１４は複数の合成画像を表示した状態を示す図である。なお、図１４においては、６通りの体動補正が行われた場合における６つの合成画像を表示した状態を示している。図１４に示すように、表示画面６１の画像表示領域６２には、複数の合成画像の縮小画像が一覧として表示される。この場合、複数の合成画像のうち、選択された合成画像を拡大して表示することが好ましい。これにより、操作者は好ましい補正方法により補正された合成画像を選択することができる。なお、この場合、図１５に示すように、最も好ましい体動補正方法により補正された合成画像以外の画像をグレーアウトさせるようにしてもよい。なお、図１５においてはグレーアウトした状態を斜線を付与することにより示している。   In this case, the body motion correction unit 40 performs body motion correction using all the 15 body motion correction methods or any number of these body motion correction methods. The image synthesizing unit 42 generates a plurality of synthesized images using the radiographic images corrected by all the performed body motion correction methods. Then, the image display unit 60 displays a plurality of composite images. FIG. 14 is a diagram showing a state in which a plurality of composite images are displayed. FIG. 14 shows a state in which six composite images are displayed when six types of body motion correction are performed. As shown in FIG. 14, reduced images of a plurality of synthesized images are displayed as a list in the image display area 62 of the display screen 61. In this case, it is preferable to enlarge and display the selected composite image among the plurality of composite images. Thereby, the operator can select the composite image corrected by a preferable correction method. In this case, as shown in FIG. 15, an image other than the synthesized image corrected by the most preferable body motion correction method may be grayed out. In FIG. 15, the grayed out state is indicated by hatching.

ここで、上記第１および第２の実施形態においては、平行移動、３次元的な動きおよび２次元的な動きの体動指標値を算出しているが、とくに長尺撮影時における体動は、ＦＰＤ１１０に対する被写体の縦横方向の移動、すなわち平行移動が主なものであることが分かっている。また、平行移動の体動は、複数の放射線画像を相対的にずらすことにより補正することができるため、放射線画像の画質を劣化させることがない。これにより、計測および診断に影響がない範囲の平行移動の体動であることが分かれば、平行移動の体動を補正しても、放射線画像を合成することにより得られる合成画像の読影に影響を与えない。したがって、長尺撮影中に体動が平行移動であり、かつ計測および診断に影響がない範囲の体動であることが分かれば、放射線画像を補正しても合成画像の読影に影響を与えないことから、再度の撮影は不要となり、その結果、被写体である患者が無用に被曝されることを防止できる。   Here, in the first and second embodiments, the body movement index value of the parallel movement, the three-dimensional movement, and the two-dimensional movement is calculated. It has been found that the movement of the subject relative to the FPD 110 in the vertical and horizontal directions, that is, the parallel movement is the main one. Further, the translational body movement can be corrected by relatively shifting a plurality of radiographic images, so that the image quality of the radiographic image is not deteriorated. As a result, if it is known that the body movement is a translational movement in a range that does not affect the measurement and diagnosis, even if the translational movement is corrected, it affects the interpretation of the synthesized image obtained by combining the radiographic images. Not give. Therefore, if it is known that the body movement is a parallel movement during long imaging and that it does not affect the measurement and diagnosis, correction of the radiation image will not affect the interpretation of the composite image. Therefore, it is not necessary to perform another imaging, and as a result, it is possible to prevent the patient as the subject from being unnecessarily exposed.

このため、体動指標値を算出するに際し、平行移動の体動指標値のみを算出するようにしてもよい。以下、これを第３の実施形態として説明する。なお、第３の実施形態においては、上記第１および第２の実施形態における体動指標値算出部３６において行われる処理のみが異なるため、ここでは構成についての詳細な説明は省略する。   For this reason, when calculating the body movement index value, only the parallel movement body movement index value may be calculated. Hereinafter, this will be described as a third embodiment. Note that in the third embodiment, only the process performed in the body motion index value calculation unit 36 in the first and second embodiments is different, and thus detailed description of the configuration is omitted here.

第３の実施形態においては、体動指標値算出部３６は、平行移動の体動指標値のみを算出する。具体的には、局所移動ベクトル算出部３４が算出した複数の局所移動ベクトルＶ０を用いて図６に示すヒストグラムを作成し、ヒストグラムに基づいて平行移動の体動指標値を算出する。ここで、撮影部位が胸部（長尺撮影でない）や全脊椎（長尺撮影）の場合、心臓または腸にたまったガスのように被写体Ｎの体動とは別に動く構造物が含まれる。すなわち、被写体Ｎを剛体と見なした場合における被写体Ｎの姿勢に関する動き以外の、局所的な動きが含まれる。体動指標値算出部３４は、作成したヒストグラムにおける頻度最大位置または頻度重心位置を用いて体動指標値を算出する。なお、頻度が大きいほど重みづけを大きくしたヒストグラムの重みづけ平均値を体動指標値として算出するようにしてもよい。   In the third embodiment, the body motion index value calculation unit 36 calculates only the translational body motion index value. Specifically, a histogram shown in FIG. 6 is created using a plurality of local movement vectors V0 calculated by the local movement vector calculation unit 34, and a body movement index value for parallel movement is calculated based on the histogram. Here, 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. That is, local motion other than the motion related to the posture of the subject N when the subject N is regarded as a rigid body is included. The body motion index value calculation unit 34 calculates a body motion index value using the frequency maximum position or the frequency centroid position in the created histogram. Note that a weighted average value of a histogram in which the weight is increased as the frequency is increased may be calculated as the body movement index value.

これにより、被写体Ｎを剛体と見なした場合における被写体Ｎの姿勢に関する動きである平行移動の体動指標値のみを、局所的な動きに影響されないように算出することができる。   As a result, it is possible to calculate only the translational movement index value that is a movement related to the posture of the subject N when the subject N is regarded as a rigid body so as not to be affected by the local motion.

なお、平行移動の体動指標値を算出する手法は、ヒストグラムを用いる手法に限定されるものではない。例えば、複数の局所移動ベクトルＶ０を方向および大きさ毎にクラスタリングすることにより局所移動ベクトルＶ０を分離し、分離した局所移動ベクトルＶ０のうち、大きさが所定のしきい値未満の局所移動ベクトルＶ０を除去する、あるいは局所移動ベクトルＶ０の大きさが大きいほど重みづけを大きくする等の手法を用いることにより、局所的な動きを抑制して、平行移動の体動指標値を算出することが可能である。   Note that the method for calculating the body motion index value of the parallel movement is not limited to the method using the histogram. For example, the local movement vector V0 is separated by clustering a plurality of local movement vectors V0 for each direction and size, and among the separated local movement vectors V0, the local movement vector V0 having a magnitude less than a predetermined threshold value. It is possible to calculate a body motion index value for parallel movement while suppressing local movement by using a technique such as removing weights or increasing weighting as the size of the local movement vector V0 increases. It is.

次いで、第３の実施形態において行われる処理について説明する。図１６は第３の実施形態において行われる処理を示すフローチャートである。まず、ＦＰＤ１１０を移動させつつ長尺撮影を行い、移動の各位置における放射線画像を取得する（ステップＳＴ４１）。そして、局所移動ベクトル算出部３４が、隣接する２つの放射線画像の重複領域における局所移動ベクトルＶ０を算出し（ステップＳＴ４２）、体動指標値算出部３６が平行移動の体動指標値を算出する（ステップＳＴ４３）。   Next, processing performed in the third embodiment will be described. FIG. 16 is a flowchart showing processing performed in the third embodiment. First, long imaging is performed while moving the FPD 110, and a radiation image at each position of movement is acquired (step ST41). Then, the local movement vector calculation unit 34 calculates the local movement vector V0 in the overlapping region of two adjacent radiographic images (step ST42), and the body movement index value calculation unit 36 calculates the body movement index value of the parallel movement. (Step ST43).

続いて、後処理選択部３９が、体動指標値に基づいて体動の有無を判別する（ステップＳＴ４４）。そして、後処理選択部３９が、体動の有無の判別結果に応じて体動補正を行うか否かを選択する（ステップＳＴ４５）。ステップＳＴ４５が肯定されると、後処理選択部３９は、平行移動の体動指標値の情報を体動補正部４０および画像表示部６０にそれぞれ出力する（情報出力、ステップＳＴ４６）。   Subsequently, the post-processing selection unit 39 determines the presence or absence of body movement based on the body movement index value (step ST44). Then, the post-processing selection unit 39 selects whether or not to perform body motion correction according to the determination result of the presence or absence of body motion (step ST45). When step ST45 is affirmed, the post-processing selection unit 39 outputs information on the parallel movement body motion index value to the body motion correction unit 40 and the image display unit 60 (information output, step ST46).

次いで、体動補正部４０が複数の放射線画像の平行移動の体動を補正し（ステップＳＴ４７）、画像合成部４２が体動補正済みの放射線画像を合成して合成画像Ｃ１を生成する（ステップＳＴ４８）。そして、画像表示部６０が合成画像Ｃ１を平行移動の体動指標値とともに表示し（ステップＳＴ４９）、処理を終了する。   Next, the body motion correction unit 40 corrects the translational body motion of the plurality of radiation images (step ST47), and the image composition unit 42 synthesizes the body motion corrected radiation images to generate a composite image C1 (step ST47). ST48). Then, the image display unit 60 displays the composite image C1 together with the translational body motion index value (step ST49), and the process ends.

一方、ステップＳＴ４５が否定されると、後処理選択部３９は、体動補正を行わない旨の選択結果を警告部４４に、平行移動の体動指標値を画像表示部６０に出力する（情報出力、ステップＳＴ５０）。警告部４４は、音声メッセージ（警告音声を含む）やブザー音（警告音）により警告を行う（ステップＳＴ５１）。さらに画像表示部６０が平行移動の体動指標値を表示し（ステップＳＴ５２）、処理を終了する。   On the other hand, if step ST45 is negative, the post-processing selection unit 39 outputs the selection result indicating that the body motion correction is not performed to the warning unit 44 and the parallel movement body motion index value to the image display unit 60 (information). Output, step ST50). The warning unit 44 issues a warning by a voice message (including a warning sound) or a buzzer sound (warning sound) (step ST51). Further, the image display unit 60 displays the body movement index value of the parallel movement (step ST52), and the process is terminated.

これにより、第３の実施形態においては、被写体の平行移動の動きを精度よく体動を検出することができる。また、検出した平行移動の体動指標値に基づいて、後段で行われる体動補正実行の有無の選択を選択する処理を適切に行うことができる。   Thereby, in the third embodiment, it is possible to accurately detect the body movement of the parallel movement of the subject. In addition, based on the detected body movement index value of the parallel movement, it is possible to appropriately perform the process of selecting whether or not to perform body movement correction performed in the subsequent stage.

なお、上記第１から第３の実施形態においては、長尺撮影により取得した放射線画像についての体動指標値を取得しているが、体動検出部３０のみを単独で使用することにより、エネルギーサブトラクション撮影、経時サブトラクション撮影、およびトモシンセシス撮影および連写撮影等、同一被写体についてショット間で体動が生じる可能性がある撮影により取得した放射線画像についても、上記第１および第２の実施形態と同様に、撮影情報に基づく局所移動ベクトルＶ０の算出、体動指標値の算出および体動の有無の判別を行うことができる。   In the first to third embodiments, the body motion index value for the radiographic image obtained by the long imaging is acquired. However, by using only the body motion detection unit 30 alone, energy can be obtained. Radiation images acquired by imaging that may cause body movement between shots of the same subject, such as subtraction imaging, temporal subtraction imaging, tomosynthesis imaging, and continuous imaging, are the same as in the first and second embodiments. In addition, it is possible to calculate the local movement vector V0 based on the imaging information, calculate the body motion index value, and determine the presence or absence of body motion.

また、上記第１から第３の実施形態においては、ＦＰＤ１１０の移動中、すなわち撮影中において、２回目以降の移動の各位置毎、すなわち２回目の撮影以降の各撮影毎に、隣接する２つの放射線画像の重複領域に基づいて体動指標値の算出および体動の有無の判別を行い、体動ありと判別された場合に警告を発生するようにしてもよい。   Further, in the first to third embodiments, during the movement of the FPD 110, that is, during shooting, two adjacent positions are moved for each position after the second shooting, that is, for each shooting after the second shooting. A body motion index value may be calculated and whether or not there is a body motion based on the overlapping region of the radiographic images, and a warning may be generated when it is determined that there is a body motion.

この場合、操作者は、長尺撮影中に警告等によって被写体Ｎの体動を確認すると、コンソール７０に設けられた緊急停止スイッチを操作することによって、撮影を途中で停止させることができる。これにより、合成ができないほどの体動量があるにも拘わらず撮影を続けることにより、被写体Ｎが無用に被曝されることを防止できる。   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, 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 post-processing selection unit 39 selects whether to perform body motion correction, selects a body motion correction method, and sets the body motion index value displayed on the image display unit 60. Although the selection is performed, any one of these or any combination may be selected.

以上、本発明の実施形態に係る装置１５０について説明したが、コンピュータを、上記の撮影情報取得部３２、局所移動ベクトル算出部３４、体動指標値算出部３６、および後処理選択部３９に対応する手段として機能させ、図９，１３，１６に示すような処理を行わせるプログラムも、本発明の実施形態の１つである。また、そのようなプログラムを記録したコンピュータ読取り可能な記録媒体も、本発明の実施形態の１つである。   Although the apparatus 150 according to the embodiment of the present invention has been described above, the computer corresponds to the imaging information acquisition unit 32, the local movement vector calculation unit 34, the body motion index value calculation unit 36, and the post-processing selection unit 39. A program that functions as a means for performing the processing as shown in FIGS. 9, 13, and 16 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.

以下、本発明の実施態様項について説明する。   Hereinafter, embodiments of the present invention will be described.

［実施態様項１］
同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得する画像取得手段と、
前記複数の放射線画像における重複領域内の局所的な前記被写体のズレを表す局所移動ベクトルを算出する局所移動ベクトル算出手段と、
前記局所移動ベクトルに基づいて前記被写体の姿勢に起因する体動指標値を算出する体動指標値算出手段とを備えたことを特徴とする体動検出装置。 [Embodiment Item 1]
An image acquisition means for acquiring a plurality of radiation images obtained by performing a plurality of times of imaging on the same subject, at least a part of the region overlapping;
A local movement vector calculating means for calculating a local movement vector representing a deviation of the subject in the overlapping region in the plurality of radiation images;
A body motion detection device comprising: a body motion index value calculating means for calculating a body motion index value resulting from the posture of the subject based on the local movement vector.

「被写体の姿勢に起因する体動指標値」とは、被写体を剛体と見なした場合における、被写体の姿勢に起因する体動の指標値を意味し、被写体内部に含まれる心臓、肺、腸等の内臓の動きに起因する指標値は含まない。   The “body motion index value resulting from the subject's posture” means the body motion index value resulting from the subject's posture when the subject is regarded as a rigid body, and includes the heart, lungs, intestines contained within the subject. The index value resulting from the movement of internal organs such as is not included.

［実施態様項２］
前記体動指標値算出手段は、前記被写体の平行移動量の指標値を前記体動指標値として算出する手段であることを特徴とする実施態様項１記載の体動検出装置。 [Embodiment 2]
2. The body motion detection device according to claim 1, wherein the body motion index value calculation means is a means for calculating an index value of the parallel movement amount of the subject as the body motion index value.

［実施態様項３］
前記体動指標値算出手段は、さらに前記被写体の３次元的な移動量の指標値および該被写体の２次元的な移動量の指標値の少なくとも一方を前記体動指標値として算出する手段であることを特徴とする実施態様項２記載の体動検出装置。 [Embodiment 3]
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 body motion detection device according to claim 2, wherein the body motion detection device is characterized in that

［実施態様項４］
前記体動指標値に基づいて、前記複数の放射線画像に関連する後処理を選択する後処理選択手段をさらに備えたことを特徴とする実施態様項３項記載の体動検出装置。 [Embodiment 4]
The body motion detection device according to claim 3, further comprising post-processing selection means for selecting post-processing related to the plurality of radiation images based on the body motion index value.

［実施態様項５］
前記後処理選択手段は、前記複数の放射線画像に体動補正処理を施すか否かを選択する手段であることを特徴とする実施態様項４記載の体動検出装置。 [Embodiment 5]
5. The body motion detection device according to claim 4, wherein the post-processing selection means is means for selecting whether or not to perform body motion correction processing on the plurality of radiation images.

［実施態様項６］
前記後処理選択手段は、前記複数の放射線画像に施す体動補正処理の種類を選択する手段であることを特徴とする実施態様項４または５記載の体動検出装置。 [Embodiment 6]
6. The body motion detection device according to claim 4, wherein the post-processing selection unit is a unit that selects a type of body motion correction processing to be performed on the plurality of radiation images.

［実施態様項７］
前記後処理選択手段は、表示する体動指標値の種類を選択する手段であることを特徴とする実施態様項４から６のいずれか１項記載の体動検出装置。 [Embodiment 7]
The body motion detection device according to any one of embodiments 4 to 6, wherein the post-processing selection unit is a unit that selects a type of body motion index value to be displayed.

［実施態様項８］
前記撮影時の撮影状況および撮影対象を表す撮影情報を取得する撮影情報取得手段をさらに備えたことを特徴とする実施態様項１から７のいずれか１項記載の体動検出装置。 [Embodiment 8]
8. The body movement detection device according to any one of claims 1 to 7, further comprising imaging information acquisition means for acquiring imaging information representing an imaging situation and an imaging object at the time of imaging.

［実施態様項９］
前記局所移動ベクトル算出手段は、前記撮影情報に基づいて前記局所移動ベクトルを算出する手段であることを特徴とする実施態様項８記載の体動検出装置。 [Embodiment 9]
9. The body movement detecting device according to claim 8, wherein the local movement vector calculating means is means for calculating the local movement vector based on the photographing information.

［実施態様項１０］
前記体動指標値算出手段は、前記撮影情報に基づいて、前記局所移動ベクトルを統合して前記平行移動の体動指標値を算出する手段であることを特徴とする実施態様項８または９記載の体動検出装置。 [Embodiment 10]
10. The embodiment according to claim 8 or 9, wherein the body motion index value calculating means is means for calculating the body motion index value of the parallel movement by integrating the local movement vectors based on the photographing information. Body motion detection device.

［実施態様項１１］
前記撮影情報および前記体動指標値に基づいて、前記被写体の体動の有無を判別する体動判別手段をさらに備えたことを特徴とする実施態様項８から１０のいずれか１項記載の体動検出装置。 [Embodiment 11]
11. The body according to any one of embodiments 8 to 10, further comprising body movement determining means for determining presence or absence of body movement of the subject based on the photographing information and the body movement index value. Motion detection device.

［実施態様項１２］
前記体動指標値に基づいて、前記被写体の体動の有無を判別する体動判別手段をさらに備えたことを特徴とする実施態様項１から１０のいずれか１項記載の体動検出装置。 [Embodiment 12]
11. The body motion detection device according to any one of embodiments 1 to 10, further comprising body motion determination means for determining presence or absence of body motion of the subject based on the body motion index value.

［実施態様項１３］
放射線検出器を移動させ、該移動により位置を変える毎に被写体を透過した放射線を該放射線検出器に照射して、少なくとも一部の領域が重複した複数の放射線画像を取得する放射線画像撮影装置であって、
前記放射線検出器を所定の移動軸に沿って移動させ、該移動により位置を変える毎に前記被写体を透過した放射線を前記放射線検出器に照射する撮影手段と、
前記移動および前記放射線の照射を行う毎に前記放射線検出器から信号を読み出すことにより、前記被写体の複数の放射線画像を取得する放射線画像取得手段と、
実施態様項１から１２のいずれか１項記載の体動検出装置とを備えたことを特徴とする放射線画像撮影装置。 [Embodiment 13]
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;
A radiographic imaging apparatus comprising the body movement detection device according to any one of the embodiments 1 to 12.

［実施態様項１４］
前記体動が検出された場合に、前記複数の放射線画像の再撮影を補助するための再撮影補助情報を生成する情報生成手段をさらに備えたことを特徴とする実施態様項１３記載の放射線画像撮影装置。 [Embodiment 14]
14. The radiographic image according to claim 13, further comprising information generating means for generating re-imaging support information for assisting re-imaging of the plurality of radiographic images when the body movement is detected. Shooting device.

［実施態様項１５］
前記再撮影補助情報を表示する表示手段をさらに備えたことを特徴とする実施態様項１４記載の放射線画像撮影装置。 [Embodiment 15]
The radiographic image capturing apparatus according to claim 14, further comprising display means for displaying the re-imaging auxiliary information.

［実施態様項１６］
前記再撮影補助情報に基づいて、前記再撮影の制御を行う再撮影制御手段をさらに備えたことを特徴とする実施態様項１４または１５記載の放射線画像撮影装置。 [Embodiment 16]
16. The radiographic image capturing apparatus according to claim 14, further comprising a re-imaging control unit that controls the re-imaging based on the re-imaging auxiliary information.

［実施態様項１７］
同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得し、
前記複数の放射線画像における重複領域内の局所的な前記被写体のズレを表す局所移動ベクトルを算出し、
前記局所移動ベクトルに基づいて前記被写体の姿勢に起因する体動指標値を算出することを特徴とする体動検出方法。 [Embodiment 17]
Obtain a plurality of radiation images obtained by performing multiple times of imaging on the same subject, at least some of the regions overlap,
Calculating a local movement vector representing a local shift of the subject in an overlapping region in the plurality of radiation images;
A body motion detection method, wherein a body motion index value resulting from the posture of the subject is calculated based on the local movement vector.

［実施態様項１８］
同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得する手順と、
前記複数の放射線画像における重複領域内の局所的な前記被写体のズレを表す局所移動ベクトルを算出する手順と、
前記局所移動ベクトルに基づいて前記被写体の姿勢に起因する体動指標値を算出する手順とを有することを特徴とする体動検出方法をコンピュータに実行させるためのプログラム。 [Embodiment 18]
A procedure for obtaining a plurality of radiographic images obtained by performing imaging a plurality of times on the same subject and overlapping at least some of the regions;
A procedure for calculating a local movement vector representing a local shift of the subject in an overlapping region in the plurality of radiation images;
A program for causing a computer to execute a body motion detection method, comprising: calculating a body motion index value resulting from the posture of the subject based on the local movement vector.

３０ 体動検出部
３１ 画像取得部
３２ 撮影情報取得部
３４ 局所移動ベクトル算出部
３６ 体動指標値算出部
３９，３９Ａ 後処理選択部
４０ 体動補正部
４２ 画像合成部
４４ 警告部
６０ 画像表示部
１００ 放射線源
１０２ 撮像面
１０４ 放射線
１１０ ＦＰＤ
１５０，１５０Ａ 放射線画像撮影装置 DESCRIPTION OF SYMBOLS 30 Body motion detection part 31 Image acquisition part 32 Imaging | photography information acquisition part 34 Local movement vector calculation part 36 Body motion index value calculation part 39,39A Post-processing selection part 40 Body motion correction part 42 Image composition part 44 Warning part 60 Image display part 100 Radiation source 102 Imaging surface 104 Radiation 110 FPD
150,150A Radiation imaging device

Claims (11)

同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得する画像取得手段と、
前記複数の放射線画像の撮影時における前記被写体の動きの種類に応じた複数種類の体動指標値を取得する体動指標値取得手段とを備えたことを特徴とする体動検出装置。 An image acquisition means for acquiring a plurality of radiation images obtained by performing a plurality of times of imaging on the same subject, at least a part of the region overlapping;
A body motion detection apparatus comprising: a body motion index value acquisition unit configured to acquire a plurality of types of body motion index values corresponding to types of motion of the subject at the time of capturing the plurality of radiation images. 前記体動指標値取得手段は、前記複数の放射線画像における重複領域内の局所的な前記被写体のズレを表す局所移動ベクトルを算出する局所移動ベクトル算出手段と、
前記局所移動ベクトルに基づいて前記複数種類の体動指標値を算出する体動指標値算出手段とを備えてなることを特徴とする請求項１記載の体動検出装置。 The body movement index value acquisition means includes a local movement vector calculation means for calculating a local movement vector representing a local deviation of the subject in an overlapping region in the plurality of radiation images;
2. The body motion detection device according to claim 1, further comprising body motion index value calculating means for calculating the plurality of types of body motion index values based on the local movement vector. 前記体動指標値算出手段は、前記被写体の平行移動量の指標値、前記被写体の３次元的な移動量の指標値および該被写体の２次元的な移動量の指標値の少なくとも２つを前記体動指標値として算出する手段であることを特徴とする請求項２項記載の体動検出装置。   The body motion index value calculating means is configured to obtain at least two of an index value of the parallel movement amount of the subject, 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. The body motion detection device according to claim 2, wherein the body motion detection device is a means for calculating the body motion index value. 前記複数種類の体動指標値または該複数種類の体動指標値のうちの主たる体動指標値に基づいて、前記複数の放射線画像に関連する後処理を選択する後処理選択手段をさらに備えたことを特徴とする請求項１から３のいずれか１項記載の体動検出装置。   And further comprising post-processing selection means for selecting post-processing related to the plurality of radiation images based on the plurality of types of body motion index values or main body motion index values of the plurality of types of body motion index values. The body motion detection device according to claim 1, wherein the body motion detection device is a body motion detection device. 前記後処理選択手段は、前記複数種類の体動指標値または前記主たる体動指標値に基づいて体動がある前記動きの種類を判定し、該判定結果に応じて前記後処理を選択する手段であることを特徴とする請求項４記載の体動検出装置。   The post-processing selection unit is configured to determine the type of motion having a body motion based on the plurality of types of body motion index values or the main body motion index value, and select the post-processing according to the determination result. The body motion detection device according to claim 4, wherein 前記後処理選択手段は、前記複数の放射線画像に体動補正処理を施すか否かを選択する手段であることを特徴とする請求項４または５記載の体動検出装置。   6. The body motion detection device according to claim 4, wherein the post-processing selection means is means for selecting whether or not to perform body motion correction processing on the plurality of radiation images. 前記後処理選択手段は、前記複数の放射線画像に施す体動補正処理の種類を選択する手段であることを特徴とする請求項４から６のいずれか１項記載の体動検出装置。   The body motion detection device according to claim 4, wherein the post-processing selection unit is a unit that selects a type of body motion correction processing to be performed on the plurality of radiation images. 前記後処理選択手段は、表示する体動指標値の種類を選択する手段であることを特徴とする請求項４から７のいずれか１項記載の体動検出装置。   The body motion detection device according to claim 4, wherein the post-processing selection unit is a unit that selects a type of body motion index value to be displayed. 放射線検出器を移動させ、該移動により位置を変える毎に被写体を透過した放射線を該放射線検出器に照射して、少なくとも一部の領域が重複した複数の放射線画像を取得する放射線画像撮影装置であって、
前記放射線検出器を所定の移動軸に沿って移動させ、該移動により位置を変える毎に前記被写体を透過した放射線を前記放射線検出器に照射する撮影手段と、
前記移動および前記放射線の照射を行う毎に前記放射線検出器から信号を読み出すことにより、前記被写体の複数の放射線画像を取得する放射線画像取得手段と、
請求項１から８のいずれか１項記載の体動検出装置とを備えたことを特徴とする放射線画像撮影装置。 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;
A radiographic imaging apparatus comprising the body movement detection device according to claim 1. 同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得し、
前記複数の放射線画像の撮影時における前記被写体の動きの種類に応じた複数種類の体動指標値を取得することを特徴とする体動検出方法。 Obtain a plurality of radiation images obtained by performing multiple times of imaging on the same subject, at least some of the regions overlap,
A body motion detection method for acquiring a plurality of types of body motion index values corresponding to types of motion of the subject at the time of capturing the plurality of radiation images. 同一被写体に関して複数回の撮影を行うことにより得られた、少なくとも一部の領域が重複した複数の放射線画像を取得する手順と、
前記複数の放射線画像の撮影時における前記被写体の動きの種類に応じた複数種類の体動指標値を取得する手順とをコンピュータに実行させることを特徴とする体動検出プログラム。 A procedure for obtaining a plurality of radiographic images obtained by performing imaging a plurality of times on the same subject and overlapping at least some of the regions;
A body motion detection program for causing a computer to execute a procedure of acquiring a plurality of types of body motion index values according to types of motion of the subject at the time of capturing the plurality of radiographic images.
JP2011178725A 2010-09-08 2011-08-18 Body motion detector, method, and program Abandoned JP2012075862A (en)

Priority Applications (10)

Application Number Priority Date Filing Date Title
JP2011178725A JP2012075862A (en) 2010-09-08 2011-08-18 Body motion detector, method, and program
EP13198106.0A EP2727533A1 (en) 2010-09-08 2011-09-07 Body motion detection device, as well as radiographic imaging apparatus
EP13173587.0A EP2653105B1 (en) 2010-09-08 2011-09-07 Body motion detection device and method, as well as radiographic imaging apparatus and method
EP13173594.6A EP2659837A1 (en) 2010-09-08 2011-09-07 Body motion detection device, as well as radiographic imaging apparatus
US13/227,184 US20120059239A1 (en) 2010-09-08 2011-09-07 Body motion detection device and method, as well as radiographic imaging apparatus and method
EP11180443A EP2428164A3 (en) 2010-09-08 2011-09-07 Body motion detection device, as well as radiographic imaging apparatus
CN201410575508.4A CN104306009A (en) 2010-09-08 2011-09-08 Radiographic image capturing apparatus
CN201110264400XA CN102429672A (en) 2010-09-08 2011-09-08 Body motion detection device and method, as well as radiographic imaging apparatus and method
US14/014,965 US20140064450A1 (en) 2010-09-08 2013-08-30 Body motion detection device and method, as well as radiographic imaging apparatus and method
US14/140,897 US20140140481A1 (en) 2010-09-08 2013-12-26 Body motion detection device and method, as well as radiographic imaging apparatus and method

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2010201075 2010-09-08
JP2010201075 2010-09-08
JP2011178725A JP2012075862A (en) 2010-09-08 2011-08-18 Body motion detector, method, and program

Publications (2)

Publication Number Publication Date
JP2012075862A true JP2012075862A (en) 2012-04-19
JP2012075862A5 JP2012075862A5 (en) 2014-08-14

Family

ID=46236786

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011178725A Abandoned JP2012075862A (en) 2010-09-08 2011-08-18 Body motion detector, method, and program

Country Status (1)

Country Link
JP (1) JP2012075862A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014136415A1 (en) * 2013-03-06 2014-09-12 富士フイルム株式会社 Body-movement detection device and method
WO2014136641A1 (en) * 2013-03-06 2014-09-12 富士フイルム株式会社 Body motion display device and body motion display method
JP2018153297A (en) * 2017-03-16 2018-10-04 コニカミノルタ株式会社 X-ray dynamic image processing device
JP2020000313A (en) * 2018-06-25 2020-01-09 富士フイルム株式会社 Imaging control device, method, and program
JP2020068870A (en) * 2018-10-29 2020-05-07 キヤノン株式会社 Image processing device, image processing method, and program
JP2022125356A (en) * 2018-09-27 2022-08-26 富士フイルム株式会社 Positional deviation amount derivation device, method and program

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002272726A (en) * 2001-03-21 2002-09-24 Konica Corp Method for identifying image and radiation image photographing apparatus
JP2006000376A (en) * 2004-06-17 2006-01-05 Konica Minolta Medical & Graphic Inc Medical image photography system, and information processor
JP2007007411A (en) * 2005-06-30 2007-01-18 Siemens Ag Method or x-ray device for creating a series of recordings of medical x-ray images
JP2007082907A (en) * 2005-09-26 2007-04-05 Canon Inc Apparatus and method for photographing image
JP2007143952A (en) * 2005-11-29 2007-06-14 Shimadzu Corp X-ray imaging apparatus
JP2009240656A (en) * 2008-03-31 2009-10-22 Fujifilm Corp X-ray radiographic apparatus
JP2010115349A (en) * 2008-11-13 2010-05-27 Fujifilm Corp Radiation tomographic apparatus

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002272726A (en) * 2001-03-21 2002-09-24 Konica Corp Method for identifying image and radiation image photographing apparatus
JP2006000376A (en) * 2004-06-17 2006-01-05 Konica Minolta Medical & Graphic Inc Medical image photography system, and information processor
JP2007007411A (en) * 2005-06-30 2007-01-18 Siemens Ag Method or x-ray device for creating a series of recordings of medical x-ray images
JP2007082907A (en) * 2005-09-26 2007-04-05 Canon Inc Apparatus and method for photographing image
JP2007143952A (en) * 2005-11-29 2007-06-14 Shimadzu Corp X-ray imaging apparatus
JP2009240656A (en) * 2008-03-31 2009-10-22 Fujifilm Corp X-ray radiographic apparatus
JP2010115349A (en) * 2008-11-13 2010-05-27 Fujifilm Corp Radiation tomographic apparatus

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014136415A1 (en) * 2013-03-06 2014-09-12 富士フイルム株式会社 Body-movement detection device and method
WO2014136641A1 (en) * 2013-03-06 2014-09-12 富士フイルム株式会社 Body motion display device and body motion display method
JP2014171487A (en) * 2013-03-06 2014-09-22 Fujifilm Corp Body motion display device and method
US9595116B2 (en) 2013-03-06 2017-03-14 Fujifilm Corporation Body motion detection device and method
US10194880B2 (en) 2013-03-06 2019-02-05 Fujifilm Corporation Body motion display device and body motion display method
JP2018153297A (en) * 2017-03-16 2018-10-04 コニカミノルタ株式会社 X-ray dynamic image processing device
JP2020000313A (en) * 2018-06-25 2020-01-09 富士フイルム株式会社 Imaging control device, method, and program
JP7122886B2 (en) 2018-06-25 2022-08-22 富士フイルム株式会社 Imaging control device, method and program
JP2022125356A (en) * 2018-09-27 2022-08-26 富士フイルム株式会社 Positional deviation amount derivation device, method and program
JP7275363B2 (en) 2018-09-27 2023-05-17 富士フイルム株式会社 Positional deviation amount derivation device, method and program
JP2020068870A (en) * 2018-10-29 2020-05-07 キヤノン株式会社 Image processing device, image processing method, and program
JP7229721B2 (en) 2018-10-29 2023-02-28 キヤノン株式会社 Image processing device, image processing method and program

Similar Documents

Publication Publication Date Title
JP5572040B2 (en) Radiography equipment
US20140140481A1 (en) Body motion detection device and method, as well as radiographic imaging apparatus and method
US7555100B2 (en) Long length imaging using digital radiography
JP5501443B2 (en) Radiation image capturing apparatus, radiation image capturing method, body movement amount measuring method, and program
JP5890598B2 (en) Method for enhancing an image of a fluoroscopic image
JP5948275B2 (en) Radiographic apparatus, radiographic method, and radiographic control program
US20070071171A1 (en) Image photographing apparatus and method
CN103705265A (en) Radiographic image generation device and method
US8184765B2 (en) Radiography apparatus
JP2012075862A (en) Body motion detector, method, and program
KR101665513B1 (en) Computer tomography apparatus and method for reconstructing a computer tomography image thereof
KR20160103518A (en) Medical image processing apparatus and method for the same
JP3897925B2 (en) Cone beam CT system
CN111803110B (en) X-ray fluoroscopic apparatus
KR101768520B1 (en) A method of integrated operation of chest X-ray digital radiography and chest digital tomosynthesis
JPH11221206A (en) Device and method for x-ray ct imaging utilizing c-arm
JP5562184B2 (en) Radiation imaging apparatus, method and program
JP2010273931A (en) X-ray imaging device
JP2012055474A (en) Body motion detection device, method and program
JP5283882B2 (en) X-ray diagnostic apparatus, image processing apparatus, and filter coefficient calculation program used for image reconstruction processing
JP6750742B2 (en) Medical X-ray image processing apparatus and X-ray imaging apparatus
EP4162882A1 (en) Computer-implemented method for positioning a patient for x-ray image recording
Chawla et al. Design and development of a new multi-projection x-ray system for chest imaging
JP2005253572A (en) Image processor, x-ray diagnosis apparatus, medical image information system and calibration table attaching method
JP2017006595A (en) Image processing apparatus, tomographic image generation system, and program

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140701

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140701

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150227

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150324

A762 Written abandonment of application

Free format text: JAPANESE INTERMEDIATE CODE: A762

Effective date: 20150501