JP2008006211A

JP2008006211A - Heart beat phase determining apparatus, program and x-ray diagnostic apparatus

Info

Publication number: JP2008006211A
Application number: JP2006182309A
Authority: JP
Inventors: Takuya Sakaguchi; 卓弥坂口
Original assignee: Toshiba Corp; Toshiba Medical Systems Corp
Current assignee: Toshiba Corp; Canon Medical Systems Corp
Priority date: 2006-06-30
Filing date: 2006-06-30
Publication date: 2008-01-17
Anticipated expiration: 2026-06-30
Also published as: JP4777164B2

Abstract

<P>PROBLEM TO BE SOLVED: To highly accurately specify a heart beat phase and a window width relatively small in the movement of the heart. <P>SOLUTION: A storage part 10 stores the data of a plurality of projection images relating to an object gathered in a fixed cycle. A feature point tracking part 13 specifies a morphological feature point included in each projection image. A determination part 15 determines the cardiac phase and a time width on the basis of the time fluctuation of the moving distance of the feature point and the electrocardiogram of the object. <P>COPYRIGHT: (C)2008,JPO&INPIT

Description

本発明は、一定周期で収集された被検体に関する複数の投影画像のデータを処理する心拍位相決定装置、プログラム及びＸ線診断装置に関する。 The present invention relates to a heartbeat phase determining apparatus, a program, and an X-ray diagnostic apparatus that process data of a plurality of projection images related to a subject collected at a constant period.

周知の通り、近年では、Ｘ線コンピュータ断層撮影装置に限らず、Ｃアームシステムでも、多方向から撮影した複数の投影画像から三次元画像を再構成する技術が実用化されつつある。Ｃアームシステムでは、Ｘ線管をＸ線検出器とともに１８０°程度回転させ、その間に２０フレーム／秒程度で撮影を繰り返している。 As is well known, in recent years, not only an X-ray computed tomography apparatus but also a C-arm system, a technique for reconstructing a three-dimensional image from a plurality of projection images photographed from multiple directions is being put into practical use. In the C-arm system, the X-ray tube is rotated about 180 ° together with the X-ray detector, and imaging is repeated at about 20 frames / second during that time.

ここで、検査対象が心臓であれば、心臓の動きに対してＣアームの回転速度が遅いので、モーションアーチファクトの影響を軽減するために、心電同期再構成のように、観察者が指定した心拍の位相（以下単に位相又は心位相という）で、それを中心とした指定幅（ウインドウの幅又は時間幅という）の期間内に収集された投影画像を、上記１８０°程度の回転期間中に収集した１セットの投影画像から部分的に抽出して、再構成処理にかけることがなされている。 Here, if the subject to be examined is the heart, the rotation speed of the C-arm is slow relative to the motion of the heart, so that the observer specified as in ECG synchronization reconstruction in order to reduce the effects of motion artifacts. Projection images collected within a specified width (window width or time width) centered on the phase of the heartbeat (hereinafter simply referred to as phase or heart phase) are rotated during the rotation period of about 180 °. A part of the collected set of projection images is extracted and subjected to reconstruction processing.

心位相を、動きの比較的少ない拡張中末期や収縮末期に指定することである程度、モーションアーチファクトの影響を軽減できるが、その位相やウインドウ幅について様々に試行錯誤を繰り返して最適な心位相や幅を絞り込む必要があった。
特開２００４−１８１２４６号公報 Although the influence of motion artifacts can be reduced to some extent by specifying the cardiac phase to the end diastole or end systole with relatively little movement, the optimal cardiac phase and width can be reduced by repeating trial and error for the phase and window width. It was necessary to narrow down.
JP 2004-181246 A

本発明の目的は、心臓の動きの比較的少ない心位相とウインドウ幅を高精度に特定することにある。 An object of the present invention is to specify a cardiac phase and a window width with relatively little heart motion with high accuracy.

本発明は、第１局面において、一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて心拍位相と時間幅とを決定する決定部とを具備する心拍位相決定装置を提供する。
本発明は、第２局面において、一定周期で収集された被検体に関する複数の投影画像各々に含まれる形態上の特徴点を特定する手段と、前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて心拍位相と時間幅とを決定する手段とをコンピュータに実現させるためのプログラムを提供する。
本発明は、第３局面において、一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて、前記移動距離が閾値未満であって、前記移動距離に対応する時間幅が最大となる心拍位相を決定する決定部とを具備することを特徴とする心拍位相決定装置を提供する。
本発明は、第４局面において、一定周期で収集された被検体に関する複数の投影画像各々に含まれる形態上の特徴点を特定する手段と、前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて、前記移動距離が閾値未満であって、前記移動距離に対応する時間幅が最大となる心拍位相を決定する手段とをコンピュータに実現させるためのプログラムを提供する。
本発明は、第５局面において、Ｘ線を発生するＸ線管と、被検体を透過したＸ線を検出するＸ線検出器と、前記Ｘ線管と前記Ｘ線検出器とを移動する移動機構と、前記Ｘ線検出器を介して一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて心拍位相と時間幅とを決定する決定部と、前記決定された心拍位相と時間幅とに基づいて前記Ｘ線の発生／停止又は前記Ｘ線の強度変化を制御する制御部とを具備することを特徴とするＸ線診断装置を提供する。
本発明は、第６局面において、Ｘ線を発生するＸ線管と、被検体を透過したＸ線を検出するＸ線検出器と、前記Ｘ線管と前記Ｘ線検出器とを移動する移動機構と、前記Ｘ線検出器を介して一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて、前記移動距離が閾値未満であって、前記移動距離に対応する時間幅が最大となる心拍位相を決定する決定部と、前記決定された心拍位相と時間幅とに基づいて前記Ｘ線の発生／停止又は前記Ｘ線の強度変化を制御する制御部とを具備することを特徴とするＸ線診断装置を提供する。 According to the first aspect of the present invention, in the first aspect, a storage unit that stores data of a plurality of projection images related to a subject collected at a constant period, a specifying unit that specifies morphological feature points included in each of the projection images, There is provided a heartbeat phase determining apparatus comprising a determining unit that determines a heartbeat phase and a time width based on a temporal variation of a moving distance of the feature point and an electrocardiogram of the subject.
According to the second aspect of the present invention, in the second aspect, means for specifying a morphological feature point included in each of a plurality of projection images related to the subject collected at a constant period, time variation of the movement distance of the feature point, and the subject There is provided a program for causing a computer to realize means for determining a heartbeat phase and a time width based on the electrocardiogram.
According to a third aspect of the present invention, in the third aspect, a storage unit that stores data of a plurality of projection images related to a subject collected at a constant period, a specifying unit that specifies morphological feature points included in each of the projection images, A determination unit configured to determine a heartbeat phase in which the movement distance is less than a threshold and the time width corresponding to the movement distance is maximum based on a temporal variation of the movement distance of the feature point and an electrocardiogram of the subject; A heartbeat phase determining device is provided.
According to the fourth aspect of the present invention, in the fourth aspect, means for specifying a morphological feature point included in each of a plurality of projection images related to the subject collected at a constant period, temporal variation in the movement distance of the feature point, and the subject Based on the electrocardiogram, a program for causing a computer to realize means for determining a heartbeat phase in which the moving distance is less than a threshold and the time width corresponding to the moving distance is maximized is provided.
In the fifth aspect, the present invention provides an X-ray tube that generates X-rays, an X-ray detector that detects X-rays that have passed through a subject, and a movement that moves between the X-ray tube and the X-ray detector. A mechanism, a storage unit that stores data of a plurality of projection images related to the subject collected at a fixed period via the X-ray detector, and a specifying unit that specifies morphological feature points included in each of the projection images A determination unit that determines a heartbeat phase and a time width based on a time variation of a moving distance of the feature point and an electrocardiogram of the subject, and the X-ray based on the determined heartbeat phase and a time width An X-ray diagnostic apparatus comprising: a control unit that controls generation / stopping of X-rays or a change in intensity of the X-ray.
In the sixth aspect, the present invention relates to an X-ray tube that generates X-rays, an X-ray detector that detects X-rays that have passed through a subject, and a movement that moves between the X-ray tube and the X-ray detector. A mechanism, a storage unit that stores data of a plurality of projection images related to the subject collected at a fixed period via the X-ray detector, and a specifying unit that specifies morphological feature points included in each of the projection images And determining a heartbeat phase at which the movement distance is less than a threshold and the time width corresponding to the movement distance is the maximum, based on the time variation of the movement distance of the feature point and the electrocardiogram of the subject. And an X-ray diagnostic apparatus comprising: a control unit configured to control generation / stop of the X-ray or intensity change of the X-ray based on the determined heartbeat phase and time width To do.

本発明によれば、心臓の動きの比較的少ない心位相とウインドウ幅を高精度に特定することができる。 According to the present invention, a cardiac phase and a window width with relatively little heart movement can be specified with high accuracy.

以下図面を参照して本発明の実施形態に係るＸ線診断装置を説明する。なお、本実施形態のＸ線診断装置は、心臓の動きの比較的少ない心拍位相及びそれを中心とした時間ウインドウの幅を決定する心拍位相決定機能を特徴的に有している。この心拍位相決定機能は、Ｘ線診断装置から分離して、心拍位相決定装置として単体で提供され得る。また、当該心拍位相決定機能は、その処理をコンピュータに実現させるためのプログラムとして提供可能であり、またそれを記憶したコンピュータで読み取り可能な記憶媒体として提供され得る。 An X-ray diagnostic apparatus according to an embodiment of the present invention will be described below with reference to the drawings. Note that the X-ray diagnostic apparatus according to the present embodiment has a heartbeat phase determining function that determines a heartbeat phase with relatively little heart motion and a time window width centered on the heartbeat phase. This heartbeat phase determination function can be provided separately as a heartbeat phase determination device separately from the X-ray diagnostic apparatus. The heartbeat phase determination function can be provided as a program for causing a computer to execute the processing, and can be provided as a computer-readable storage medium that stores the program.

図１に示すように、Ｘ線診断装置は、Ｃアーム装置５を有する。Ｃアーム装置５は、Ｃアーム１６と、Ｃアーム１６を直交３軸に関して回転自在に支持する床置き又は天井吊り支持機構と、回転駆動源とを有する。Ｃアーム１６の一端にはＸ線管１が取り付けられる。Ｘ線制御部４は、システム制御部９の制御に従って、Ｘ線管１からＸ線を発生するために、Ｘ線管１の電極間に管電圧を印加し、またＸ線管１の陰極フィラメントに加熱電流を供給する。Ｃアーム１６の他端にはＸ線検出器２が取り付けられる。Ｘ線管１とＸ線検出器２とは天板３上の被検体Ｐを挟んで対向する。Ｘ線検出器２は、例えばイメージインテンシファイアとＴＶカメラとの組み合わせから構成される。または、Ｘ線検出器２は、マトリクス状に配列された半導体検出素子を有するフラットパネルディテクタ（ＦＰＤ：平面型Ｘ線検出器）で構成される。Ｃアーム／天板移動制御部６は、システム制御部９の制御に従って、Ｃアーム１６を回転するためにその駆動源に電力を供給する。またＣアーム／天板移動制御部６は、システム制御部９の制御に従って、天板を移動するためにその駆動源に電力を供給する。 As shown in FIG. 1, the X-ray diagnostic apparatus has a C-arm device 5. The C arm device 5 includes a C arm 16, a floor or ceiling suspension support mechanism that supports the C arm 16 so as to be rotatable about three orthogonal axes, and a rotational drive source. The X-ray tube 1 is attached to one end of the C arm 16. The X-ray control unit 4 applies a tube voltage between the electrodes of the X-ray tube 1 to generate X-rays from the X-ray tube 1 in accordance with the control of the system control unit 9, and the cathode filament of the X-ray tube 1 To supply a heating current. The X-ray detector 2 is attached to the other end of the C arm 16. The X-ray tube 1 and the X-ray detector 2 face each other with the subject P on the top 3 interposed therebetween. The X-ray detector 2 is composed of, for example, a combination of an image intensifier and a TV camera. Alternatively, the X-ray detector 2 includes a flat panel detector (FPD: planar X-ray detector) having semiconductor detection elements arranged in a matrix. The C arm / top movement control unit 6 supplies electric power to the drive source for rotating the C arm 16 according to the control of the system control unit 9. The C arm / top plate movement control unit 6 supplies power to the drive source for moving the top plate according to the control of the system control unit 9.

本Ｘ線診断装置では、被検体Ｐを計測して心電図を発生するために心電計７が装備される。心電計７で発生された心電図は、後述するように最適心位相及び最適ウインドウ幅の決定処理、心電同期撮影等に用いられる。画像データ記憶部１０は、Ｘ線検出器２で発生された複数の投影画像に関するデータを、心電図からシステム制御部９により取得される心位相のデータを関連付けて記憶する。なお、心位相とは、心電図のＲ波から次のＲ波までの期間内の各時点を規定する尺度として定義され、典型的には、Ｒ波から次のＲ波までの期間を１００％に規格化し、当該期間内の各時点をパーセントの単位で表記する。Ｃアーム１６を回転しながら撮影を繰り返すいわゆる回転撮影においては、各投影画像には心位相とともに撮影角度（Ｘ線管１の回転角度）が関連付けられる。画像再構成処理部１１は、画像データ記憶部１０に記憶された撮影角度が相違する複数の投影画像から断層画像を再構成する。モニタ１２は、主に、投影画像、断層画像、後述する時間ウインドウ幅及び心位相に対する特徴点の移動距離の関係性を表す立体的図表や等高線図表を表示するために設けられている。操作部８は、ユーザからの各種指令をシステム制御部９に伝達するために設けられ、キーボード、マウス等の各種入力デバイスを有する。 In the present X-ray diagnostic apparatus, an electrocardiograph 7 is equipped to measure the subject P and generate an electrocardiogram. The electrocardiogram generated by the electrocardiograph 7 is used for determining the optimum cardiac phase and optimum window width, electrocardiographic synchronization imaging, and the like, as will be described later. The image data storage unit 10 stores data regarding a plurality of projection images generated by the X-ray detector 2 in association with cardiac phase data acquired by the system control unit 9 from the electrocardiogram. The cardiac phase is defined as a scale that defines each time point in the period from the R wave of the electrocardiogram to the next R wave. Typically, the period from the R wave to the next R wave is 100%. Normalize and indicate each point in the period in percent. In so-called rotational imaging in which imaging is repeated while rotating the C-arm 16, an imaging angle (a rotation angle of the X-ray tube 1) is associated with each projection image together with a cardiac phase. The image reconstruction processing unit 11 reconstructs a tomographic image from a plurality of projection images having different shooting angles stored in the image data storage unit 10. The monitor 12 is provided mainly for displaying a projection image, a tomographic image, a three-dimensional chart or a contour chart showing the relationship between the movement distance of the feature point with respect to a time window width and a cardiac phase, which will be described later. The operation unit 8 is provided to transmit various commands from the user to the system control unit 9, and includes various input devices such as a keyboard and a mouse.

本実施形態では、心臓の動きの比較的少ない心位相とウインドウ幅を高精度に特定するために、特徴点追跡部１３、移動距離計算部１４、最適心位相／最適ウインドウ幅決定部１５が装備されている。特徴点追跡部１３は、操作部８を介して任意の投影画像上に指定された解剖学上形態的な特徴を有する特徴点を他の各投影画像上で同定する。移動距離計算部１４は、撮影時刻が前後する一対の投影画像上で心臓の動きに主に起因して変位する特徴点間距離、つまり一対の投影画像の撮影間隔で特徴点が移動する移動距離を計算する。最適心位相／最適ウインドウ幅決定部１５は、計算された移動距離の時間変化と、ユーザ指定された移動距離の上限値（閾値）とに基づいて、移動距離が閾値未満であって、ウインドウ幅が最大となる心位相を最適心位相として、またその最大のウインドウ幅を最適ウインドウ幅として決定する。 In the present embodiment, a feature point tracking unit 13, a movement distance calculation unit 14, and an optimal cardiac phase / optimum window width determination unit 15 are provided in order to specify a cardiac phase and a window width with relatively little heart motion. Has been. The feature point tracking unit 13 identifies feature points having anatomical morphological features specified on any projection image via the operation unit 8 on each other projection image. The movement distance calculation unit 14 is a distance between feature points that are displaced mainly due to the movement of the heart on a pair of projection images whose shooting times are before and after, that is, a movement distance that the feature points move at the shooting interval of the pair of projection images. Calculate The optimum cardiac phase / optimum window width determination unit 15 determines that the movement distance is less than the threshold based on the calculated time variation of the movement distance and the upper limit value (threshold value) of the movement distance designated by the user, and the window width Is determined as the optimal cardiac phase, and the maximum window width is determined as the optimal window width.

以下、最適心位相、最適ウインドウ幅の決定について詳細に説明する。まず、使用する各種用語について、図２、図３を参照して、定義をする。回転撮影とは、Ｃアーム１６の回転によりＸ線管１をＸ線検出器２とともに被検体の周囲を回転させ、その間、複数の投影画像を撮影する。本実施形態では、回転期間の全期間にわたってＸ線を継続的に被検体に照射して、撮影を繰り返すのではなく、特定の期間に限定してＸ線を被検体に照射して、その期間だけ撮影をする。特定の期間を心電図に対応付けた撮影法が心電同期撮影と呼ばれる。撮影が行われる特定の期間を、ここでは「時間ウインドウ」と定義し、その時間長を「時間ウインドウの幅」と定義する。また、時間ウインドウの中心に対応する心位相を関心位相と称する。本実施形態では、関心位相と、時間ウインドウの幅とを最適化して、心臓の動きにより断層画像上に顕在化するモーションアーチファクトの発生をできるだけ抑えることを目的としている。 Hereinafter, determination of the optimum cardiac phase and the optimum window width will be described in detail. First, various terms used will be defined with reference to FIGS. In the rotation imaging, the X-ray tube 1 and the X-ray detector 2 are rotated around the subject by the rotation of the C-arm 16, and a plurality of projection images are captured during the rotation. In the present embodiment, X-rays are continuously irradiated to the subject over the entire rotation period, and imaging is not repeated, but the subject is irradiated with X-rays limited to a specific period. Just shoot. An imaging method in which a specific period is associated with an electrocardiogram is called electrocardiogram synchronous imaging. A specific period during which photographing is performed is defined as a “time window”, and the time length is defined as a “time window width”. The cardiac phase corresponding to the center of the time window is called the phase of interest. The object of the present embodiment is to optimize the phase of interest and the width of the time window so as to suppress the occurrence of motion artifacts that appear on a tomographic image due to the motion of the heart as much as possible.

図４には、本実施形態によるシステム制御部９の制御のもとで実行される動作全体の手順を示している。図５にはその補足図を示している。まず、第１段階において、ユーザ指令にしたがってＣアーム１６が被検体に対して適当な角度で固定され、その方向から一連の複数の投影画像が撮像される（Ｓ１）。複数の投影画像は、動画を構成し、例えば２０フレーム／秒、つまり１秒あたり２０枚を撮影する速度で、複数の心拍周期、好ましくは少なくとも３心拍周期にわたって継続的に撮影される。撮影された投影画像の枚数は、Ｎと仮定する。複数の投影画像の撮影と並行して心電図が計測され、システム制御部９で特定された心位相（％）のデータとともに複数の投影画像のデータが、画像データ記憶部１０に記憶される。 FIG. 4 shows a procedure of the entire operation executed under the control of the system control unit 9 according to the present embodiment. FIG. 5 shows a supplementary diagram thereof. First, in the first stage, the C-arm 16 is fixed at an appropriate angle with respect to the subject in accordance with a user command, and a series of projection images are taken from that direction (S1). The plurality of projection images constitute a moving image, and are continuously captured over a plurality of heartbeat cycles, preferably at least three heartbeat cycles, at a speed of shooting 20 frames / second, that is, 20 frames per second. Assume that the number of photographed projection images is N. The electrocardiogram is measured in parallel with the photographing of the plurality of projection images, and the data of the plurality of projection images is stored in the image data storage unit 10 together with the data of the cardiac phase (%) specified by the system control unit 9.

次に、第２段階において、まず最初に撮影した第１フレームの投影画像が、画像データ記憶部１０から読み出され、モニタ１２に表示される。観察者は操作部８を介して表示画像上に特徴点を指定する（Ｓ２）。図６Ａは狭窄部分、図６Ｂは分岐血管、図６Ｃはカテーテル等のワイヤ先端、図６Ｄはステントに装備されたＸ線不透過性を有するマーカーをそれぞれ特徴点として指定した様子を例示している。特徴点には同定しやすいと考えられる形態的特徴が明瞭な部位が指定される。 Next, in the second stage, the first frame of the projected image of the first frame is read from the image data storage unit 10 and displayed on the monitor 12. The observer designates feature points on the display image via the operation unit 8 (S2). 6A illustrates a stenosis portion, FIG. 6B illustrates a branch vessel, FIG. 6C illustrates a distal end of a wire such as a catheter, and FIG. 6D illustrates a state in which a radiopaque marker mounted on a stent is designated as a feature point. . As the feature point, a region with a clear morphological feature that is considered to be easily identified is designated.

特徴点追跡部１３により、２フレーム目以降の各投影画像から特徴点が抽出され、各位置が特定される（Ｓ３）。特徴点の位置としては、心臓の収縮方向に近い画像のＸ軸又はＹ軸いずれか一方の座標が特定される。しかし、特徴点の位置としては、Ｘ軸又はＹ軸いずれか一方の一次元位置に限定されず、ＸＹ座標で表される二次元位置として特定してもよい。図７に特徴点の位置が時間とともに変位する様子を示している。拍動にしたがって特徴点の位置が周期的に変位する様子が理解される。 The feature point tracking unit 13 extracts feature points from the projection images of the second and subsequent frames, and specifies each position (S3). As the position of the feature point, the coordinates of either the X axis or the Y axis of the image close to the contraction direction of the heart are specified. However, the position of the feature point is not limited to the one-dimensional position of either the X axis or the Y axis, and may be specified as a two-dimensional position represented by XY coordinates. FIG. 7 shows how the positions of feature points are displaced with time. It is understood that the position of the feature point is periodically displaced according to the beat.

次に、移動距離計算部１４では、図７、図１１Ａに示した複数心拍にわたる特徴点の位置の時間変化を、典型的にはＲ波を基準として心拍周期毎に区分し、各区分を心位相を揃えて規格化して並列する。なお、各区分を心位相を揃えて並列する手法に限定されず、図１１Ｂに示すように各区分の後端を合わせてそこからの戻し時間（絶対時間）で揃えて並列しても良いし、図１１Ｃに示すように各区分の前端を合わせてそこからの遅れ時間（絶対時間）を揃えて並列しても良い。 Next, the movement distance calculation unit 14 divides temporal changes in the positions of feature points over a plurality of heartbeats shown in FIG. 7 and FIG. 11A, typically for each heartbeat cycle based on the R wave. Standardize and align in parallel. Note that the method is not limited to the method in which the respective sections are aligned in parallel with each other, and the rear ends of the respective sections may be aligned and aligned in return time (absolute time) as shown in FIG. 11B. As shown in FIG. 11C, the front ends of the respective sections may be combined and the delay time (absolute time) from there may be aligned and paralleled.

この特徴点の位置の時間変化に対して、移動距離計算部１４では、時間ウインドウの候補を適用する。この時間ウインドウの候補の時間幅（％で表す）の中で特徴点の最大移動距離（以下単に移動距離という）を計算する。つまり、図８に示すように、時間ウインドウ候補の時間幅の中で、最も離れた２つの特徴点の位置を選択し、選択した２つの特徴点の位置の間の長さを移動距離として計算する。この移動距離の計算は、時間ウインドウの候補の中心位相を１％間隔で移動しながら繰り返される。それにより図９に示すように、移動距離と心位相との関係が得られる。さらに、この関係の中で、最小移動距離が特定される。この最小移動距離は、当該処理に使用した時間ウインドウ候補の時間幅に関する固有値とされる。つまり、当該時間幅では、心臓の動きが最も少ない心位相がその移動距離（最小移動距離）とともに特定される。上記処理は、時間ウインドウ候補の時間幅を、予め決められた範囲を予め決められた間隔で変えながら繰り返される。それにより、図１０、図１２に示すように、時間ウインドウ候補の時間幅と最小移動距離との関係が得られる。また、図１３に示すように、心位相と時間ウインドウ候補の時間幅とに対して移動距離を分布させ、同じ移動距離を等高線のように結んで表現した図表が得られる。 The movement distance calculation unit 14 applies a time window candidate to the time change of the position of the feature point. Within this time window candidate time width (expressed in%), the maximum moving distance of feature points (hereinafter simply referred to as moving distance) is calculated. That is, as shown in FIG. 8, the position of the two most distant feature points in the time window candidate time width is selected, and the length between the positions of the two selected feature points is calculated as the movement distance. To do. This calculation of the moving distance is repeated while moving the center phase of the time window candidate at 1% intervals. Thereby, as shown in FIG. 9, the relationship between the movement distance and the cardiac phase is obtained. Furthermore, the minimum movement distance is specified in this relationship. This minimum moving distance is an eigenvalue related to the time width of the time window candidate used for the processing. In other words, in the time width, the cardiac phase with the least movement of the heart is specified together with the movement distance (minimum movement distance). The above processing is repeated while changing the time width of the time window candidate by changing a predetermined range at a predetermined interval. Thereby, as shown in FIGS. 10 and 12, the relationship between the time width of the time window candidate and the minimum moving distance is obtained. Further, as shown in FIG. 13, a chart is obtained in which the moving distance is distributed with respect to the cardiac phase and the time width of the time window candidate, and the same moving distance is connected like a contour line.

最適心位相／最適ウインドウ幅決定部１５は、例えば図１３に示すように、移動距離が閾値未満であって、時間ウインドウの時間幅が最大となる心位相（最適心位相）と時間ウインドウの時間幅（時間ウインドウの最適な時間幅）とを決定する（Ｓ５）。時間ウインドウの時間幅を拡大することで、画像再構成に用いることのできる投影画像の枚数を増加させることができ、それにより再構成画像（断層画像）の画質向上を図ることができる。 For example, as shown in FIG. 13, the optimum cardiac phase / optimum window width determination unit 15 has a cardiac phase (optimum cardiac phase) and a time window time at which the moving distance is less than the threshold and the time width of the time window is the maximum. The width (the optimal time width of the time window) is determined (S5). By expanding the time width of the time window, it is possible to increase the number of projection images that can be used for image reconstruction, thereby improving the image quality of the reconstructed image (tomographic image).

なお、最適心位相／最適ウインドウ幅決定部１５は、心位相の変化及びウインドウ時間幅の変化に対する特徴点の移動距離（ウインドウ時間幅内の最大移動距離）の変化を表す図表を作成する機能を有する。当該図表として典型的には、図１５に示すように、心位相、ウインドウ候補の時間幅及び移動距離を直交３軸に対応させた立体的な図表であり、また図１６に示すように、心位相とウインドウ候補の時間幅との直交２軸上で移動距離が同じ点を結んだ等高線を含む図表である。 The optimum cardiac phase / optimum window width determination unit 15 has a function of creating a chart that represents changes in the movement distance of the feature point (maximum movement distance within the window time width) with respect to changes in the cardiac phase and the window time width. Have. Typically, as shown in FIG. 15, the chart is a three-dimensional chart in which the cardiac phase, the window candidate time width and the movement distance correspond to the three orthogonal axes, and as shown in FIG. It is a graph containing the contour line which connected the point with the same movement distance on two orthogonal axes of a phase and the time width of a window candidate.

上述の通り決定した最適な心位相を中心とした最適な時間幅を有する時間ウインドウにしたがってシステム制御部９は撮像条件を設定し（Ｓ６）、心電計７からの心電図を使って心電同期により回転撮像を実行する（Ｓ７）。つまり、Ｘ線管１及びＸ線検出器２をＣアーム１６とともに一定速度で回転し（図１７参照）、その回転期間において、上記決定した最適な心位相を中心とした最適な時間幅を有する時間ウインドウの期間に限定してＸ線を発生し、被検体を透過したＸ線に関する投影画像を収集する（図１８、図１９、図２０、図２１参照）。 In accordance with a time window having an optimal time width centered on the optimal cardiac phase determined as described above, the system control unit 9 sets imaging conditions (S6) and uses the electrocardiogram from the electrocardiograph 7 to synchronize the electrocardiogram. Thus, rotational imaging is executed (S7). That is, the X-ray tube 1 and the X-ray detector 2 are rotated at a constant speed together with the C-arm 16 (see FIG. 17), and have an optimal time width centered on the determined optimal cardiac phase during the rotation period. X-rays are generated only during the period of the time window, and projection images relating to X-rays transmitted through the subject are collected (see FIGS. 18, 19, 20, and 21).

上記決定された最適な心位相を中心とした最適な時間幅を有する時間ウインドウの期間は、心臓の動きが通常は上記閾値未満を示すので、収集された投影画像から断層画像を再構成したとき（Ｓ８）、その断層画像は心臓の動きに起因するアーチファクトが軽減されている。 When the tomographic image is reconstructed from the acquired projection image, the duration of the time window having the optimal time width centered on the determined optimal cardiac phase is usually less than the threshold value. (S8) In the tomographic image, artifacts due to the motion of the heart are reduced.

なお、上述では、ウインドウ候補の時間幅内での移動距離を、特徴点の１次元又は２次元の位置から計算したが、図２２に示すように、バイプレーン撮影により２方向の投影画像を収集し、２方向の投影画像からエピポーラ理論により特徴点の３次元的な位置を求め、３次元空間上の特徴点の移動距離を計算するようにしてもよい。 In the above description, the movement distance within the time width of the window candidate is calculated from the one-dimensional or two-dimensional position of the feature point. However, as shown in FIG. 22, two-direction projection images are collected by biplane shooting. Then, the three-dimensional position of the feature point may be obtained from the projection image in two directions by the epipolar theory, and the movement distance of the feature point in the three-dimensional space may be calculated.

また、上述の説明では、心臓の動きを特徴点の移動距離によりその指標としたが、移動距離に代えて、図２３に示すように、左室造影撮影により得られた画像内の黒化部分（造影領域）の面積（ピクセル数）を用いてもよい。 In the above description, the motion of the heart is used as an index based on the moving distance of the feature point. Instead of the moving distance, as shown in FIG. 23, the blackened portion in the image obtained by left ventricular contrast imaging is used. The area (number of pixels) of (contrast region) may be used.

また、上述では、回転撮像の例を説明したが、図２４に示すように、一方向に固定した動画撮影においても適用できる。 In the above description, the example of rotational imaging has been described. However, as illustrated in FIG. 24, the present invention can also be applied to moving image shooting fixed in one direction.

また、上述では、最適な心位相を中心とした最適な時間幅を有する時間ウインドウの期間に限定してＸ線を発生したが、その期間を他の期間と比較して、パルスレートを変え（遅く）、Ｘ線強度を変え（強くする又は弱くする）、Ｘ線照射時間を変え（長くする）、Ｘ線量を変え（増やす）、またはリカーシブを変える（過去画像を多くする）ようにしてもよい。 In the above description, X-rays are generated only in the period of the time window having the optimum time width centered on the optimum cardiac phase. However, the period is compared with other periods, and the pulse rate is changed ( (Slow), change X-ray intensity (intensify or weaken), change X-ray irradiation time (increase), change X-ray dose (increase), or change recursive (increase past images) Good.

また、上述では、心拍位相決定機能を、回転撮影型のＸ線診断装置へ適用する例を説明したが、２方向からの立体画像構成装置にもそれ特有の高い効果とともに適用することができる。２方向からの立体画像構成は、図２５に示すように、２方向それぞれで被検体を静止画もしくは動画で撮影する。撮影した２方向の画像間の同じ部位の位置ずれから立体的な位置を特定する。この装置で心臓血管のように動く被検体を撮像する場合、２枚の画像は心臓が同じ位相にある２画像を使う必要がある。この際、もし心臓の動きが激しい心位相付近を用いると、撮像タイミングが少しずれただけで心臓が大きく動いてしまうため、２枚の画像が同じ心位相にならず、よって立体画像の精度が悪くなる。そこで本実施形態による心拍位相決定法を適応し、心臓の動きが比較的小さい心位相を選択し、その付近の画像を用いるようにする。これにより撮像タイミングが少しずれても、心臓はほとんど動いていないため、立体画像の精度を良く保てる。 In the above description, the example in which the heartbeat phase determination function is applied to the rotational imaging type X-ray diagnostic apparatus has been described. In the stereoscopic image configuration from two directions, as shown in FIG. 25, the subject is photographed as a still image or a moving image in each of the two directions. A three-dimensional position is specified from the positional shift of the same part between the two images taken. When a subject moving like a cardiovascular vessel is imaged with this apparatus, it is necessary to use two images in which the heart is in the same phase. At this time, if the vicinity of the heart phase where the heart motion is intense is used, the heart moves greatly only with a slight shift in the imaging timing, so the two images do not have the same heart phase, so the accuracy of the stereoscopic image is improved. Deteriorate. Therefore, the heartbeat phase determination method according to the present embodiment is applied, a heart phase with relatively small heart motion is selected, and an image in the vicinity thereof is used. As a result, even if the imaging timing is slightly deviated, the heart hardly moves, so that the accuracy of the stereoscopic image can be maintained well.

なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。 Note that the present invention is not limited to the above-described embodiment as it is, and can be embodied by modifying the constituent elements without departing from the scope of the invention in the implementation stage. In addition, various inventions can be formed by appropriately combining a plurality of components disclosed in the embodiment. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, constituent elements over different embodiments may be appropriately combined.

本発明の実施形態に係るＸ線診断装置の構成を示す図。The figure which shows the structure of the X-ray diagnostic apparatus which concerns on embodiment of this invention. 本実施形態による心電同期スキャンの概要を示す図。The figure which shows the outline | summary of the electrocardiogram synchronous scan by this embodiment. 図２の補足図。The supplementary figure of FIG. 本実施形態の動作手順を示す流れ図。The flowchart which shows the operation | movement procedure of this embodiment. 本実施形態の全体動作の概要を示す図。The figure which shows the outline | summary of the whole operation | movement of this embodiment. 図５の特徴点として狭窄部分を例示する図。The figure which illustrates a constriction part as a feature point of FIG. 図５の特徴点として分岐血管を例示する図。The figure which illustrates a branch blood vessel as a feature point of FIG. 図５の特徴点としてワイヤ先端を例示する図。The figure which illustrates a wire tip as a feature point of FIG. 図５の特徴点としてステントマーカを例示する図。The figure which illustrates a stent marker as a feature point of FIG. 図５のＳ３の特徴点の位置の時間変化を示す図。The figure which shows the time change of the position of the feature point of S3 of FIG. 図５のＳ４における特徴点の移動距離の計算方法を補足する図。The figure which supplements the calculation method of the movement distance of the feature point in S4 of FIG. 図５のＳ４における特徴点の移動距離に関する時間変化を示す図。The figure which shows the time change regarding the movement distance of the feature point in S4 of FIG. 図５のＳ４における特徴点の最小移動距離に関する時間幅変化を示す図。The figure which shows the time width change regarding the minimum moving distance of the feature point in S4 of FIG. 図５のＳ３の特徴点位置の時間変化を平均化する第１の例を示す図。The figure which shows the 1st example which averages the time change of the feature point position of S3 of FIG. 図５のＳ３の特徴点位置の時間変化を平均化する第２の例を示す図。The figure which shows the 2nd example which averages the time change of the feature point position of S3 of FIG. 図５のＳ３の特徴点位置の時間変化を平均化する第３の例を示す図。The figure which shows the 3rd example which averages the time change of the feature point position of S3 of FIG. 図５のＳ４における特徴点の最小移動距離に関する時間幅変化の作成処理を示す図。The figure which shows the creation process of the time width change regarding the minimum moving distance of the feature point in S4 of FIG. 図５のＳ４における特徴点の移動距離に関する時間変化を複数の時間ウインドウ幅について同時に示す図。The figure which shows the time change regarding the moving distance of the feature point in S4 of FIG. 5 simultaneously about several time window width | variety. 図５のＳ４における最適な心位相と最適な時間ウインドウ幅の選択法を示す図。The figure which shows the selection method of the optimal cardiac phase in S4 of FIG. 5, and the optimal time window width. 図５のＳ４において、時間ウインドウ幅及び心位相に対する特徴点の移動距離の関係性を表す立体的図表を示す図。The figure which shows the three-dimensional chart showing the relationship of the movement distance of the feature point with respect to time window width and a cardiac phase in S4 of FIG. 図５のＳ４において、時間ウインドウ幅及び心位相に対する特徴点の移動距離の関係性を表す等高線図表を示す図。The figure which shows the contour-line chart showing the relationship of the movement distance of the feature point with respect to time window width and a cardiac phase in S4 of FIG. 本実施形態によるＸ線管の焦点軌道を示す図。The figure which shows the focal track | orbit of the X-ray tube by this embodiment. 本実施形態において、図１７の焦点軌道に対して撮影位置を示す図。In this embodiment, the figure which shows an imaging position with respect to the focal track | orbit of FIG. 本実施形態において、心電波形に対する撮影位置を示す図。The figure which shows the imaging | photography position with respect to an electrocardiogram waveform in this embodiment. 本実施形態において、焦点軌道及び心電波形に対する撮影位置を示す図。The figure which shows the imaging | photography position with respect to a focal track | orbit and an electrocardiogram waveform in this embodiment. 図２０の焦点位置に対する検出器の位置を示す図。The figure which shows the position of the detector with respect to the focus position of FIG. 本実施形態において、特徴点の移動距離を３次元に拡張する変形例を示す図。The figure which shows the modification which expands the moving distance of a feature point to three dimensions in this embodiment. 本実施形態において、特徴点の移動距離に代えて左室造影部分の面積を適用する変形例を示す図。The figure which shows the modification which replaces with the movement distance of a feature point in this embodiment, and applies the area of the left ventricular contrast part. 本実施形態の変形例において、動作手順を示す流れ図。The flowchart which shows an operation | movement procedure in the modification of this embodiment. 本実施形態において、立体画像構成装置への適用例を示す図。The figure which shows the example applied to a three-dimensional image structure apparatus in this embodiment.

符号の説明Explanation of symbols

１…Ｘ線管、２…Ｘ線検出器、３…天板、４…Ｘ線制御部、５…Ｃアーム機構、６…移動制御部、７…心電計、８…操作部、９…システム制御部、１０…画像データ記憶部、１１…画像再構成処理部、１２…モニタ、１３…特徴点追跡部、１４…移動距離計算部、１５…最適心位相／最適ウインドウ幅決定部。 DESCRIPTION OF SYMBOLS 1 ... X-ray tube, 2 ... X-ray detector, 3 ... Top plate, 4 ... X-ray control part, 5 ... C arm mechanism, 6 ... Movement control part, 7 ... Electrocardiograph, 8 ... Operation part, 9 ... System control unit, 10 ... image data storage unit, 11 ... image reconstruction processing unit, 12 ... monitor, 13 ... feature point tracking unit, 14 ... moving distance calculation unit, 15 ... optimum cardiac phase / optimum window width determination unit.

Claims

一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、
前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、
前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて心拍位相と時間幅とを決定する決定部とを具備することを特徴とする心拍位相決定装置。 A storage unit for storing data of a plurality of projection images related to the subject collected at a fixed period;
A specifying unit for specifying a feature point on the form included in each of the projection images;
A heartbeat phase determination apparatus comprising: a determination unit that determines a heartbeat phase and a time width based on a temporal variation of a moving distance of the feature point and an electrocardiogram of the subject.

一定周期で収集された被検体に関する複数の投影画像各々に含まれる形態上の特徴点を特定する手段と、
前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて心拍位相と時間幅とを決定する手段とをコンピュータに実現させるためのプログラム。 Means for identifying morphological feature points included in each of a plurality of projection images relating to a subject collected at a fixed period;
A program for causing a computer to realize means for determining a heartbeat phase and a time width based on a time variation of a moving distance of the feature point and an electrocardiogram of the subject.

一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、
前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、
前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて、前記移動距離が閾値未満であって、前記移動距離に対応する時間幅が最大となる心拍位相を決定する決定部とを具備することを特徴とする心拍位相決定装置。 A storage unit for storing data of a plurality of projection images related to the subject collected at a fixed period;
A specifying unit for specifying a feature point on the form included in each of the projection images;
A determination unit configured to determine a heartbeat phase in which the movement distance is less than a threshold and the time width corresponding to the movement distance is maximum based on a temporal variation of the movement distance of the feature point and an electrocardiogram of the subject; A heartbeat phase determining device comprising:

一定周期で収集された被検体に関する複数の投影画像各々に含まれる形態上の特徴点を特定する手段と、
前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて、前記移動距離が閾値未満であって、前記移動距離に対応する時間幅が最大となる心拍位相を決定する手段とをコンピュータに実現させるためのプログラム。 Means for identifying morphological feature points included in each of a plurality of projection images relating to a subject collected at a fixed period;
Means for determining a heartbeat phase at which the moving distance is less than a threshold and the time width corresponding to the moving distance is maximum based on a temporal variation of the moving distance of the feature point and an electrocardiogram of the subject; A program to be realized on a computer.

Ｘ線を発生するＸ線管と、
被検体を透過したＸ線を検出するＸ線検出器と、
前記Ｘ線管と前記Ｘ線検出器とを移動する移動機構と、
前記Ｘ線検出器を介して一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、
前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、
前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて心拍位相と時間幅とを決定する決定部と、
前記決定された心拍位相と時間幅とに基づいて前記Ｘ線の発生／停止又は前記Ｘ線の強度変化を制御する制御部とを具備することを特徴とするＸ線診断装置。 An X-ray tube that generates X-rays;
An X-ray detector for detecting X-rays transmitted through the subject;
A moving mechanism for moving the X-ray tube and the X-ray detector;
A storage unit for storing data of a plurality of projection images related to the subject collected at a constant cycle via the X-ray detector;
A specifying unit for specifying a feature point on the form included in each of the projection images;
A determination unit that determines a heartbeat phase and a time width based on a time variation of a moving distance of the feature point and an electrocardiogram of the subject;
An X-ray diagnostic apparatus comprising: a control unit that controls generation / stop of the X-rays or change in intensity of the X-rays based on the determined heartbeat phase and time width.

Ｘ線を発生するＸ線管と、
被検体を透過したＸ線を検出するＸ線検出器と、
前記Ｘ線管と前記Ｘ線検出器とを移動する移動機構と、
前記Ｘ線検出器を介して一定周期で収集された被検体に関する複数の投影画像のデータを記憶する記憶部と、
前記投影画像各々に含まれる形態上の特徴点を特定する特定部と、
前記特徴点の移動距離の時間変動と前記被検体の心電図とに基づいて、前記移動距離が閾値未満であって、前記移動距離に対応する時間幅が最大となる心拍位相を決定する決定部と、
前記決定された心拍位相と時間幅とに基づいて前記Ｘ線の発生／停止又は前記Ｘ線の強度変化を制御する制御部とを具備することを特徴とするＸ線診断装置。 An X-ray tube that generates X-rays;
An X-ray detector for detecting X-rays transmitted through the subject;
A moving mechanism for moving the X-ray tube and the X-ray detector;
A storage unit for storing data of a plurality of projection images related to the subject collected at a constant cycle via the X-ray detector;
A specifying unit for specifying a feature point on the form included in each of the projection images;
A determination unit configured to determine a heartbeat phase in which the movement distance is less than a threshold and the time width corresponding to the movement distance is maximum based on a temporal variation of the movement distance of the feature point and an electrocardiogram of the subject; ,
An X-ray diagnostic apparatus comprising: a control unit that controls generation / stop of the X-rays or change in intensity of the X-rays based on the determined heartbeat phase and time width.

前記心拍位相の変化及び前記時間幅の変化に対する前記特徴点の移動距離の変化を表す図表を作成する作成部をさらに備えることを特徴とする請求項６記載のＸ線診断装置。 The X-ray diagnostic apparatus according to claim 6, further comprising a creation unit that creates a chart representing a change in a moving distance of the feature point with respect to a change in the heartbeat phase and a change in the time width.

前記作成部は、前記心拍位相、前記時間幅及び前記移動距離を直交３軸に対応させた立体的な図表を作成することを特徴とする請求項７記載のＸ線診断装置。 The X-ray diagnostic apparatus according to claim 7, wherein the creating unit creates a three-dimensional chart in which the heartbeat phase, the time width, and the movement distance are associated with three orthogonal axes.

前記作成部は、前記移動距離が同じ点を結んだ等高線を含む図表を作成することを特徴とする請求項７記載のＸ線診断装置。 The X-ray diagnostic apparatus according to claim 7, wherein the creating unit creates a chart including contour lines connecting the same movement distances.