(57)【要約】【課題】  主に気管支や血管などの管状構造における分岐位置を特定するにあたって、その処理工数を低減し処理時間を短縮すること。【解決手段】  医用画像処理装置は、被検体の３次元領域を対象とするボリュームデータを記憶する記憶部を有する。断面画像発生部は、ボリュームデータファイルから所定の基準軸に略直交する複数の断面にそれぞれ対応する複数の断面画像のデータを発生する。領域抽出部は、複数の断面画像から閾値処理を用いて目的部位に関する複数の領域を抽出する。位置特定部は、抽出された領域の数が変化する基準軸上の位置を特定する。【選択図】図１
2012096024.tif 000002
【特許請求の範囲】
【請求項１】
被検体の３次元領域を対象とするボリュームデータを記憶する記憶部と、
前記ボリュームデータファイルから所定の基準軸に略直交する複数の断面にそれぞれ対応する複数の断面画像のデータを発生する断面画像発生部と、
前記複数の断面画像から閾値処理を用いて目的部位に関する複数の領域を抽出する領域抽出部と、
前記抽出された領域の数が変化する前記基準軸上の位置を特定する位置特定部とを具備することを特徴とする医用画像処理装置。
【請求項２】
前記基準軸は前記被検体の体軸であることを特徴とする請求項１記載の医用画像処理装置。
【請求項３】
前記基準軸を前記被検体の体軸に初期的に設定し、前記抽出した領域の位置に基づいて前記基準軸を修正するための基準軸設定部をさらに備えることを特徴とする請求項１記載の医用画像処理装置。
【請求項４】
前記基準軸を操作者指示に従って修正するための基準軸設定部をさらに備えることを特徴とする請求項１記載の医用画像処理装置。
【請求項５】
前記ボリュームデータファイルから前記閾値処理により抽出した前記目的部位に関する３次元領域に基づいて前記基準軸を設定するための基準軸設定部をさらに備えることを特徴とする請求項１記載の医用画像処理装置。
【請求項６】
前記目的部位は管状部位であり、
前記基準軸設定部は、前記３次元領域から特定した前記管状部位の軸線に基づいて前記基準軸を設定することを特徴とする請求項５記載の医用画像処理装置。
【請求項７】
前記領域の重心に基づいて前記基準軸を修正する基準軸設定部をさらに備えることを特徴とする請求項１記載の医用画像処理装置。
【請求項８】
前記領域抽出部は、前記閾値処理により抽出した複数の領域候補から画素数に基づいて前記目的部位の領域を抽出することを特徴とする請求項１記載の医用画像処理装置。
【請求項９】
前記領域抽出部は、前記閾値処理により抽出した複数の領域候補から領域周囲長に基づいて前記目的部位の領域を抽出することを特徴とする請求項１記載の医用画像処理装置。
【請求項１０】
前記領域抽出部は、前記閾値処理により抽出した複数の領域候補から直径に基づいて前記目的部位の領域を抽出することを特徴とする請求項１記載の医用画像処理装置。
【請求項１１】
前記領域抽出部は、前記閾値処理により抽出した複数の領域候補から最大径又は最小径に基づいて前記目的部位の領域を抽出することを特徴とする請求項１記載の医用画像処理装置。
【請求項１２】
前記目的部位は気管支又は血管であり、前記位置特定部は前記気管支又は血管の分岐位置を前記特定された断面位置により特定することを特徴とする請求項１記載の医用画像処理装置。
【請求項１３】
前記位置特定部は、前記領域の数が１から２に、又は２から１に変化する位置を前記気管支の第１分岐の位置として特定することを特徴とする請求項１２記載の医用画像処理装置。
【請求項１４】
前記位置特定部は、前記領域の数が１だけ増減する位置を特定することを特徴とする請求項１記載の医用画像処理装置。
【請求項１５】
前記発生された複数の断面画像から選択した、前記特定された位置を含む一の断面画像を表示する表示部をさらに備えることを特徴とする請求項１記載の医用画像処理装置。
【請求項１６】
前記表示部は、前記表示された断面画像に前記特定された位置を示すマークを重畳することを特徴とする請求項１５記載の医用画像処理装置。
【請求項１７】
前記発生された複数の断面画像から選択した、前記特定された位置を中心とした所定距離内の複数の断面画像を表示する表示部をさらに備えることを特徴とする請求項１記載の医用画像処理装置。
【請求項１８】
前記ボリュームデータファイルから前記目的部位に関する３次元画像を発生する３次元画像発生部と、
前記３次元画像を前記特定された位置を示すマークを重ねて表示する表示部をさらに備えることを特徴とする請求項１記載の医用画像処理装置。
【請求項１９】
前記断面画像発生部は、操作者指示に従って前記断面画像を発生する範囲を限局することを特徴とする請求項１記載の医用画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、医用画像処理装置に関する。
【背景技術】
【０００２】
例えばＸ線コンピュータ断層画像撮影装置ではＸ線検出器の列数が格段に増加し、肺野全体の投影データを一度に収集することが可能になりつつあり、それによる呼吸や拍動による動きアーチファクトの影響の低いボリュームデータを用いて診断支援情報の自動提供への応用が様々に試みられている。例えば肺野全体に関するボリュームデータを得て、気管支の分岐部を抽出し、その位置の断面画像を表示する機能がある。
【０００３】
その処理では、３次元上で気管支領域を抽出し、その管状構造に対して芯線を同定し、芯線構造から分岐部を特定している。
【０００４】
気管支領域抽出処理、芯線化処理、分岐部特定処理は全て３次元上での処理であり、非常に多い処理工数、非常に長い処理時間を要するものであった。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開2010-136765号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
目的は、主に気管支や血管などの管状構造における分岐位置を特定するにあたって、その処理工数を低減し処理時間を短縮することにある。
【課題を解決するための手段】
【０００７】
本実施形態による医用画像処理装置は、被検体の３次元領域を対象とするボリュームデータを記憶する記憶部を有する。断面画像発生部は、ボリュームデータファイルから所定の基準軸に略直交する複数の断面にそれぞれ対応する複数の断面画像のデータを発生する。領域抽出部は、複数の断面画像から閾値処理を用いて目的部位に関する複数の領域を抽出する。位置特定部は、抽出された領域の数が変化する基準軸上の位置を特定する。
【図面の簡単な説明】
【０００８】
【図１】本実施形態に係る医用画像処理装置を含むＸ線コンピュータ断層画像撮影装置の全体構成を示す図である。
【図２】本実施形態において気管支分岐位置特定処理の手順を示す流れ図である。
【図３】図１のＳ１１における探索範囲の例を示す図である。
【図４】図１のＳ１２で発生される断面画像の断面配列を示す図である。
【図５】図１のＳ１３で発生される二値画像例を示す図である。
【図６】図１のＳ１７で判定される気管支領域数の変化を示す図である。
【図７】本実施形態において気管支分岐位置特定処理の他の手順を示す流れ図である。
【発明を実施するための形態】
【０００９】
以下、図面を参照しながら本実施形態に係る医用画像処理装置を説明する。  
なお、医用画像処理装置は、Ｘ線コンピュータ断層撮影装置、磁気共鳴映像装置（ＭＲＩ）、超音波診断装置、Ｘ線診断装置などの被検体の３次元領域に関するボリュームデータを発生可能な医用画像発生装置に適合する。本実施形態に係わる医用画像処理装置は、これら医用画像発生装置に組み込まれ、又は単独で機能する。単独で機能するとき、本実施形態に係わる医用画像処理装置は、ＬＡＮなどの電気的通信回線に接続され、電気的通信回線を介して医用画像発生装置から又は病院内若しくは外部の医用画像補間通信システム（ＰＡＣＳ）から処理対象のボリュームデータを受信する。ここでは、本実施形態に係わる医用画像処理装置は、Ｘ線コンピュータ断層撮影装置に組み込まれたものとして説明する。
【００１０】
図１に、本実施形態に係る医用画像処理装置を装備したＸ線コンピュータ断層撮影装置の構成をブロック図により示している。架台部１００は、回転自在に支持される回転フレーム１０２を有する。回転フレーム１０２の回転中心軸をＺ軸、水平方向をＸ軸、垂直方向をＹ軸として説明する。撮影に際して回転フレーム１０２内側の撮影領域Ｓ内に挿入される被検体の体軸はＺ軸に略一致する。
【００１１】
架台駆動部１０７はホストコントローラ１１０の制御の下で回転フレーム１０２を回転駆動するための駆動信号を発生する。回転フレーム１０２にはコーンビーム形Ｘ線管１０１と２次元検出器（エリア検出器とも呼ばれる）１０３とがＺ軸を中心とした撮影領域Ｓを挟んで対向して搭載される。高電圧発生器１０９はホストコントローラ１１０の制御の下でＸ線管１０１に管電流を供給し、また高電圧を両極間に印加する。それによりＸ線管１０１の焦点ＦからＸ線絞り１１１で四角錐形に製形されたＸ線が被検体に照射される。２次元検出器１０３は、個々にチャンネルを構成する複数のＸ線検出素子を有する。複数のＸ線検出素子はＸ線焦点Ｆを中心としてＸＺ各軸を中心として略円弧状に配列される。
【００１２】
２次元検出器１０３には一般的にＤＡＳ(data acquisition system) と呼ばれているデータ収集装置１０４が接続される。データ収集装置１０４には、２次元検出器１０３の各チャンネルの電流信号を電圧に変換するＩ−Ｖ変換器と、この電圧信号をＸ線の曝射周期に同期して周期的に積分する積分器と、この積分器の出力信号を増幅するアンプと、このプリアンプの出力信号をディジタル信号に変換するアナログ・ディジタル・コンバータとが、チャンネルごとに設けられている。データ収集装置１０４の出力には光学的又は磁気的要素を媒介させる非接触データ伝送装置１０５を介して前処理装置１０６が接続される。前処理装置１０６は、データ収集装置１０４で検出されたデータに対して、チャンネル間の感度不均一を補正し、またＸ線強吸収体、主に金属部による極端な信号強度の低下又は信号脱落を補正する等の前処理を実行する。前処理装置１０６で前処理を受けたデータ（投影データという）はデータ記憶部１１６を介してコーンビーム再構成処理部１１２に供給される。
【００１３】
コーンビーム再構成処理部１１６は、ホストコントローラ１１０の制御のもとで、３６０°又は（１８０素＋ファン角）の角度範囲分の投影データに基づいて例えばコーンビーム再構成法によりＣＴ値分布を３次元座標系（ｘｙｚ）で表現してなるボリュームデータを再構成する。ボリュームデータの３次元座標系（ｘｙｚ）は、実空間座標系（ＸＹＺ）に対応する。再構成されたボリュームデータはデータ記憶部１１６に記憶される。
【００１４】
ボリュームデータから気管支や血管の芯線を特定してから分岐位置を同定する従来の３次元処理よりも処理工数、処理量を減らし、且つ一定の精度を維持するために、本実施形態では以下の通りマルチスライスで断面画像を生成し、その２次元画像上での処理で分岐位置を同定する。
【００１５】
断面画像発生部１１５は、表示部（ディスプレイ）１１３への表示に適合するいわゆるＭＰＲ処理（断面変換処理）により、ボリュームデータから、基準線設定部１２５により設定された基準線に略直交する複数の断面にそれぞれ対応する複数の断面画像を発生する。複数の断面は、基準軸に沿って一定の間隔（スライスピッチ）で略平行に配列され、いわゆるマルチスライスを構成する。複数の断面画像のデータはデータ記憶部１１６に記憶される。断面変換条件としてスライスピッチ、断面画像の空間分解能を操作者は操作デバイス１１５を介して任意に設定することができる。
【００１６】
基準線設定部１２５は、基準線を初期的に被検体の体軸に設定する。通常、被検体の体軸がＺ軸にほぼ一致するように被検体の天板上での位置が調整される。基準線設定部１２５は、操作デバイス１１５を介して入力された操作者指示に従って基準線の位置及び向きを修正する。また基準線設定部１２５は、操作デバイス１１５を介して入力された操作者による第１の基準線自動修正トリガに従って後述する気管支領域判定部１１９により複数の断面画像から抽出された気管支領域の位置に基づいて基準軸を修正する。典型的には、基準線は、複数の断面画像の中の気管支上部を横断する少なくとも２枚の断面画像から抽出された気管支領域の重心位置を連結する線に修正される。また基準線設定部１２５は、操作デバイス１１５を介して入力された操作者による第２の基準線自動修正トリガに従って後述する二値化処理部１１７により抽出された気管支に関する３次元領域から気管支上部の芯線に基づいて基準軸を修正する。
【００１７】
二値化処理部１１７は、複数の断面画像各々に対して閾値処理により複数の二値画像を発生する。この閾値は処理対象の部位に固有のＣＴ値に応じて決められる。例えば処理対象部位が気管支であれば空気ガスを抽出対象として−８００から−１０００の範囲の何れかの値に設定され、当該閾値未満のＣＴ値を有する一塊の画素群が気管支領域候補として抽出される。処理部位が造影血管であれば造影剤を抽出対象として＋８００から＋１０００の範囲の何れかの値に設定され、当該閾値を超過するＣＴ値を有する一塊の画素群が血管領域候補として抽出される。なおここでは処理対象部位を気管支として説明する。
【００１８】
このような領域候補を抽出するための二値化条件として、閾値、閾値に対する未満／超過の区別を操作者は操作デバイス１１５を介して任意に設定することができる。なお二値化処理は、下限閾値と上限閾値とから決まる範囲のＣＴ値を有する画素を抽出するようにしてもよく、その場合、二値化条件としては抽出方法の選択、下限閾値、上限閾値が設定される。
【００１９】
気管支領域判定部１１９は、各二値画像内の気管支領域候補各々が気管支領域であるか否かを所定の判定規則に従って判定する。操作者は、操作デバイス１１５を介して、次の１）−４）の判定方法のいずれか、又は２以上の組み合わせを選択する。１）気管支領域候補の面積に基づいて当該判定を行う。面積は画素数に画素あたりの単位面積を乗算することにより得られる。従って画素数は面積に実質的に等価である。具体的には、各気管支領域候補を構成する画素数を計数し、計数した画素数が閾値、例えば１００ピクセルを超過しているとき、その気管支領域候補を気管支領域として判定する。２）気管支領域候補の周囲長に基づいて当該判定を行う。周囲長は周囲画素数に画素あたりの単位長を乗算することにより得られる。従って周囲画素数は周囲長に実質的に等価である。具体的には、各気管支領域候補の外縁を構成する画素数を計数し、計数した画素数が閾値、例えば５０ピクセルを超過しているとき、その気管支領域候補を気管支領域として判定する。３）気管支領域候補の直径に基づいて当該判定を行う。「直径」は、気管支領域候補の重心を通り、かつ所定方向と平行な直線が、気管支領域候補を横断する長さとして定義される。直径は当該横断線上の画素数に画素あたりの単位長を乗算することにより得られる。従って当該画素数は直径に実質的に等価である。具体的には、当該画素数を計数し、計数した画素数が閾値、例えば２０ピクセルを超過しているとき、その気管支領域候補を気管支領域として判定する。４）気管支領域候補の最大径又は最小径に基づいて当該判定を行う。「最大径（又は最小径）」は、気管支領域候補の重心を通る複数の直線各々が気管支領域候補を横断する長さ（画素数）を計数し、その最大画素数（又は最小画素数）により得られる。当該最大画素数（又は最小画素数）が閾値、例えば２０ピクセルを超過しているとき、その気管支領域候補を気管支領域として判定する。なお、ここでは判定方法１）を例に説明する。
【００２０】
気管支領域計数部１２１は、判定した気管支領域を各二値画像（各スライス）ごとに、１、２・・・等の連番でラベリングする。ラベリングにより各二値画像ごとに気管支領域の数が決定される。なお、対象部位領域の判定処理は、その対象部位が気管支であれば、上記の通り閾値超過する領域サイズを有する領域候補を領域として判定することが好ましいが、他の部位であればそれぞれの部位に応じて閾値未満の領域サイズを有する領域候補を当該部位領域として判定するようにしてもよいし、所定範囲の領域サイズを有する領域候補を当該部位領域として判定するようにしてもよい。このような領域判定のための判定条件として、閾値、閾値に対する未満／超過の区別を操作者は操作デバイス１１５を介して任意に設定することができ、また判定方法の選択、下限閾値、上限閾値が設定される。
【００２１】
分岐位置判定部１２３は、気管支領域計数部１２１により各スライスごとに判定した気管支領域の数が近接スライスに対して変化するスライス位置により気管支の第１分岐位置を同定する。つまり、分岐位置判定部１２３は、気管支領域の数が１だけ増減する基準軸上の位置を分岐位置として特定する。具体的には、分岐位置判定部１２３は、被検体の頭部側から下肢側に向かって気管支領域の数を探索し、その領域数が１から２に変化するスライス位置を特定し、そのスライスのＺ位置、直前スライスのＺ位置、又は特定したスライスとその直前のスライスとの中間位置を気管支の第１分岐位置（Ｚ位置）として同定する。逆向きに探索しても良く、被検体の下肢側から頭部側に向かって気管支領域の数を探索し、その領域数が２から１に変化するスライス位置を特定するようにしてもよい。判定条件は、判定対象部位及び判定分岐状態に応じて任意に設定可能であり、変化前後の数、探索方向は操作デバイス１１５を介して任意に設定することができる。
【００２２】
断面画像発生部１１５は、ホストコントローラ１１０の制御のもとで、同定された分岐位置を含む断面画像をボリュームデータから自動的に発生する。同定された分岐位置を含む断面画像の生成条件、例えば空間分解能、補間処理の有無、ＸＹＺ各軸に対するオブリーク角等は操作者が操作デバイス１１５を介して予め任意に設定可能である。
【００２３】
表示部１１３は当該分岐位置を含む断面に関する断面画像を、分岐位置を示すマークを重畳して表示する。また表示部１１３は当該分岐位置を中心とした所定距離内の複数断面、典型的には３断面の断面画像を同一画面に並べて表示する。当該分岐位置を含む断面画像には、分岐位置を示すマークが重畳される。表示部１１３は、３次元処理部１２７によりボリュームデータファイルから気管支に関する３次元画像をレンダリング処理等により発生された３次元画像を表示する。この３次元画像には、分岐位置を示すマークが重ねられる。
【００２４】
図２には本実施形態において気管支分岐位置特定処理の手順を示している。操作者により操作デバイス１１５を介して処理対象のボリュームデータが指定され、各種処理条件が設定される（Ｓ１１）。処理条件としては、上述したとおり、断面変換条件、二値化条件、領域判定条件、分岐位置判定が含まれる。ここでは、処理対象部位が気管支であり、断面変換条件としてスライスピッチをΔＳＰ、断面画像の空間分解能をボリュームデータオリジナルの空間分解能に等価であると設定し、二値化条件として閾値ｔｈ１未満との条件に設定し、領域判定条件としてその閾値ｔｈ２を１００ピクセル未満との条件に設定し、分岐位置判定条件として頭部側から下肢側に向かって領域数が１から２に変化するとの条件に設定されるものとして説明する。
【００２５】
なお、処理開始に先立って、図３に示すように、図示しないボリュームレンダリング処理部による３次元画像を表示部１１３に表示し、操作者が操作デバイス１１５を介して気管支の主要部分を含む探索範囲を限局するようにしてもよい。この探索範囲内に限局して、断面画像発生から分岐位置探索までの処理が実行される。
【００２６】
図４に示すように、断面画像発生部１１５により、ボリュームデータからスライスピッチΔＳＰの間隔で基準軸に略直交するＸＹ面に平行な断面（スライス）に関する断面画像が発生される（Ｓ１２）。説明の便宜上、複数のスライスを頭部側から順番にＳ１，Ｓ２，Ｓ３，・・・と表記する。まず、探索順において最初のスライスＳ１に関する断面画像ＳＩ１が生成される。
【００２７】
断面画像ＳＩ１は、二値化処理部１１７により、図５に示すように閾値ｔｈ１により二値画像ＢＩ１に変換される（Ｓ１３）。二値画像ＢＩ１には閾値ｔｈ１より低いＣＴ値を有する画素が閾値ｔｈ１以上のＣＴ値を有する画素と区別されている。図５に示すように閾値ｔｈ１より低いＣＴ値を有する画素が連結する複数の領域が複数の気管支領域候補として気管支領域判定部１１９により抽出され、複数の気管支領域候補各々を構成する画素数が計数される（Ｓ１４）。複数の気管支領域候補全てについて、それぞれ計数された画素数が閾値ｔｈ２、例えば１００画素（ピクセル）と気管支領域判定部１１９で比較される（Ｓ１５）。気管支領域候補全てが閾値ｔｈ２以下の画素数からなるとき、工程Ｓ１２にリターンし、次のスライスＳ２についてＳ１２−Ｓ１５の処理が同様に実施される。画素数が閾値ｔｈ２より多い領域サイズを有する気管支領域候補は、気管支領域として気管支領域判定部１１９により判定される。
【００２８】
判定された気管支領域の数が気管支領域計数部１２１により計数される（Ｓ１６）。分岐位置判定部１２３では、気管支領域の数が分岐位置判定条件に従って変化するスライスを特定し、そのスライスに従って分岐位置（Ｚ位置）を同定する。図５、図６に示すように頭部側のスライスＳ１では、気管支領域の数は“１”である。分岐より下肢側では“２”に変化する。分岐位置は、頭部側から探索して１から２に最初に変化するスライスに従って分岐位置（Ｚ位置）が同定される。
【００２９】
図５、図６に示すように気管支領域の数が“１”から“２”に変化するｎ番目のスライスＳnの直前スライスＳn-1まで、領域数が１に維持されるので、Ｓ１２−Ｓ１７の処理が繰り返される。スライスＳnでは、領域数が“２”で最初に計数される。領域数が“２”で最初に計数されたスライスＳnを特定し歌時点で、Ｓ１２−Ｓ１７の処理は終了する。
【００３０】
分岐位置判定部１２３は、領域数が１から２に変化したスライスＳnのＺ位置を分岐位置として同定する（Ｓ１８）。スライスＳnの直前スライスＳn-1のＺ位置、又は特定したスライスＳnとその直前のスライスＳn-1との中間位置を気管支の分岐位置（Ｚ位置）として同定してもよい。それらのいずれの位置を分岐位置とするかは操作者の任意であり、スライスピッチΔＳＰが１ボクセルピッチから数ボクセルピッチ程度の比較的短いときはスライスＳnのＺ位置を分岐位置として同定することが好ましく、逆にスライスピッチΔＳＰが数十ボクセルピッチ程度の比較的長いときは特定したスライスＳnとその直前のスライスＳn-1との中間位置を分岐位置として同定することが好ましい。
【００３１】
工程Ｓ１９において、断面画像発生部１１５により、同定された分岐位置を含む断面画像がボリュームデータから発生され、表示部１１３に表示される。上述したとおり、この最終的に発生される分岐位置を含む断面画像は、操作者により予め設定された空間分解能でＸＹＺ各軸に対するオブリーク角に対応する断面に関して発生される。
【００３２】
本実施形態では芯線を抽出せずに気管支の第1分岐を抽出する。芯線を抽出する処理をなくすことができるため、処理時間を短縮することができる。
【００３３】
図２の処理手順は図７の処理手順に変形可能である。図２では領域候補からその領域サイズにより領域を判定してから、その領域数の変化により分岐位置を判定したが、図７の例では、領域候補の数の変化を判定してから（Ｓ１７）、それら全ての領域候補の領域サイズが領域判定条件を満たしているか否かを判定し（Ｓ２０）、満たしているときにそのスライスを分岐位置として同定する。この手順であっても、図２の処理手順より処理量は多くなるものの、図２の処理手順と同様に分岐位置を精度良く判定することができ、しかもその処理量は従来の芯線抽出処理を含む処理の処理量より少なくなるという効果を奏することができる。
【００３４】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００３５】
１００…架台部、１０１…Ｘ線管、１０２…回転フレーム、１０３…２次元検出器、１０４…データ収集装置、１０６…前処理装置、１０７…架台駆動部、１０９…高電圧発生器、１１０…ホストコントローラ、１１２…投影データ記憶部、１１３…ボリュームデータファイル記憶部、１１５…断面画像発生部、１１６…コーンビーム再構成処理部、１１７…表示部（ディスプレイ）、１１８…再生コントローラ、１２０…呼吸センサ、１２１…呼吸波形記憶部、１２３…最大値／最小値特定部、１２５…断面位置決定部。
【図１】
2012096024.tif 000003
【図２】
2012096024.tif 000004
【図３】
2012096024.tif 000005
【図４】
2012096024.tif 000006
【図５】
2012096024.tif 000007
【図６】
2012096024.tif 000008
【図７】
2012096024.tif 000009