JP4346297B2 - Ｘ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 - Google Patents
Ｘ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 Download PDFInfo
- Publication number
- JP4346297B2 JP4346297B2 JP2002307168A JP2002307168A JP4346297B2 JP 4346297 B2 JP4346297 B2 JP 4346297B2 JP 2002307168 A JP2002307168 A JP 2002307168A JP 2002307168 A JP2002307168 A JP 2002307168A JP 4346297 B2 JP4346297 B2 JP 4346297B2
- Authority
- JP
- Japan
- Prior art keywords
- polar
- heart
- model
- coronary artery
- unit
- 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.)
- Expired - Lifetime
Links
- 238000012545 processing Methods 0.000 title claims description 42
- 238000002591 computed tomography Methods 0.000 title claims description 14
- 238000003672 processing method Methods 0.000 title claims description 10
- 210000004351 coronary vessel Anatomy 0.000 claims description 72
- 230000004217 heart function Effects 0.000 claims description 45
- 230000002107 myocardial effect Effects 0.000 claims description 35
- 210000000056 organ Anatomy 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 claims description 14
- 230000000747 cardiac effect Effects 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 230000033001 locomotion Effects 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 7
- 230000002194 synthesizing effect Effects 0.000 claims description 5
- 230000007659 motor function Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 3
- 238000003786 synthesis reaction Methods 0.000 claims 3
- 238000000034 method Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 15
- 230000002861 ventricular Effects 0.000 description 14
- 210000004165 myocardium Anatomy 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 210000005240 left ventricle Anatomy 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 239000000470 constituent Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 210000001367 artery Anatomy 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 238000001990 intravenous administration Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/466—Displaying means of special interest adapted to display 3D data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/503—Clinical applications involving diagnosis of heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/504—Clinical applications involving diagnosis of blood vessels, e.g. by angiography
Description
【０００１】
【発明の属する技術分野】
本発明は、心臓の心機能指標を極座標で表現するポーラーマップの生成及び表示機能を備えたＸ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法に関する。
【０００２】
【従来の技術】
ヘリカルスキャンの高速化及び２次元検出器を使ったボリュームスキャンの開発は、Ｘ線コンピュータ断層撮影装置の概念を、単なる断面画像化装置からボリューム映像化化装置へ大きく移行させた。つまり、これらスキャン形式を備えたＸ線コンピュータ断層撮影装置により、被検体のボリュームを短い時間分解能で高速にスキャンすることが可能になった。このような機能を備えたＸ線コンピュータ断層撮影装置の実用化と共に、これまで使われていなかった多くの診断分野への進出が試みられ、また全く新しい診断方法を模索する動きが始まっている。
【０００３】
その一つの方向性として、心臓の心拍運動に関する評価診断への適用がある。ボリュームスキャンは、セグメントスキャン方式とあいまって、心臓全体を短いサイクルで連続的にスキャンすることを実質的に実現している。それにより例えば拡張末期（ＥＤ）や収縮末期（ＥＳ）に限定した心臓の立体像を生成することができ、さらに、拡張末期（ＥＤ）から収縮末期（ＥＳ）までの心壁厚の変化率等の心機能指標を心臓全体にわたって求めることも可能となる。
【０００４】
心臓は心尖部から心基部まで縦に長い。この縦に長い心臓全体の心機能指標を一度に観察できるように、その表示形式として一般的にポーラーマップ(Polar Map)と呼ばれる極座標画像で表示されることが多い。周知のとおり、極座標は、平面を（ｒ，θ）で表す表現方法であり、ここでは、θを心軸周りの角度、ｒを心尖部から心基部に向かって各スライスまでの距離にそれぞれ割り当てている。
【０００５】
このポーラーマップにより心機能低下を起こしている部分が容易に観察できるものであるが、実際に診断に用いる場合には、その部分と支配血管との関係が必要とされる。そのため、図１９に示すように、解剖学的に分割された既成の一般的な支配領域区分を重ねて表示することが多い。図１９において、“ＬＡＤ”表記の区分は左前下行枝の支配領域を表し、“ＬＣＸ”表記の区分は左回旋枝の支配領域を表し、“ＲＣＳ”表記の区分は右冠状動脈の支配領域を表している。
【０００６】
このように一般的な支配領域区分をポーラーマップに重ねて表示することにより、左前下行、枝左回旋枝、右冠状動脈のいずれの冠状動脈に障害が起きているか、つまりどの動脈の支配区域で障害が起きているかまでは把握することができるものの、その冠状動脈の果たして具体的にどの場所、どの位置で障害が起きているかは判断することはできない。
【０００７】
【発明が解決しようとする課題】
本発明の目的は、心臓の心機能指標を極座標で表現するポーラーマップの生成及び表示機能を備えたＸ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法において、障害が起きている冠状動脈の支配区域の判別だけでなく、その具体的な冠状動脈の障害位置まで特定又は推定できる情報の提供を実現することにある。
【０００８】
【課題を解決するための手段】
本発明の第１局面による画像処理装置は、被検体の特定臓器を含むボリュームデータから前記特定臓器の運動機能に関する指標を計算する指標計算部と、前記指標を極座標系に分布することによりポーラーマップを生成するポーラーマップ生成部と、前記ボリュームデータから前記特定臓器又はその周辺臓器の形態を、前記極座標系で表現したポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーモデルを前記ポーラーマップに合成する合成部と、前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備する。
本発明の第２局面によるＸ線コンピュータ断層撮影装置は、被検体の心臓を含む領域に関する投影データを収集する架台部と、前記投影データに基づいて画像データを再構成する再構成部と、前記画像データを、心軸に関して離散的な複数の短軸像のデータに変換する断面変換処理部と、前記短軸像各々から前記心臓の心筋領域を抽出する心筋抽出部と、前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算する心機能指標計算部と、前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成するポーラーマップ生成部と、前記短軸像各々から前記心臓の冠状動脈領域を抽出する冠状動脈抽出部と、前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーモデルを前記ポーラーマップに合成する合成部と、前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備する。
本発明の第３局面による画像処理装置は、被検体の心臓を含む領域に関する画像データを、心軸に関して離散的な複数の短軸像のデータに変換する断面変換処理部と、前記短軸像各々から前記心臓の心筋領域を抽出する心筋抽出部と、前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算する心機能指標計算部と、前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成するポーラーマップ生成部と、前記短軸像各々から前記心臓の冠状動脈領域を抽出する冠状動脈抽出部と、前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーモデルを前記ポーラーマップに合成する合成部と、前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備する。
本発明の第４局面による画像処理方法は、被検体の心臓を含む領域に関する画像データを、心軸に関して離散的な複数の短軸像のデータに変換し、前記短軸像各々から前記心臓の心筋領域を抽出し、前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算し、前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成し、前記短軸像各々から前記心臓の冠状動脈領域を抽出し、前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成し、前記ポーラーモデルを前記ポーラーマップに合成し、前記ポーラーモデルが合成されたポーラーマップを表示する。
本発明の第５局面によるＸ線コンピュータ断層撮影装置は、被検体の心臓を含む領域に関する投影データを収集する架台部と、前記投影データに基づいて画像データを再構成する再構成部と、前記画像データから前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成するポーラーマップ生成部と、前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーマップを前記ポーラーモデルとともに表示する表示部とを具備する。
本発明の第６局面による画像処理装置は、被検体の心臓を含む領域に関する画像データから、前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成するポーラーマップ生成部と、前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーマップを前記ポーラーモデルとともに表示する表示部とを具備する。 本発明の第７局面による画像処理方法は、被検体の心臓を含む領域に関する画像データから、前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成し、前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成し、前記ポーラーマップを前記ポーラーモデルとともに表示する。
【０００９】
【発明の実施の形態】
以下、図面を参照して本発明によるＸ線コンピュータ断層撮影装置（Ｘ線ＣＴ装置）、画像処理装置及び画像処理方法の実施形態を説明する。なお、以下では、本実施形態は、Ｘ線ＣＴ装置として説明する。しかし、以下のＸ線ＣＴ装置は、本実施形態の画像処理方法を実現する機能を有する画像処理装置を装備している。
【００１０】
ここで、Ｘ線ＣＴ装置には、Ｘ線管と放射線検出器とが１体として被検体の周囲を回転する回転／回転（ROTATE/ROTATE）タイプと、リング状に多数の検出素子がアレイされ、Ｘ線管のみが被検体の周囲を回転する固定／回転（STATIONARY/ROTATE）タイプ等様々なタイプがあり、いずれのタイプでも本発明を適用可能である。ここでは、現在、主流を占めている回転／回転タイプとして説明する。また、１スライスの断層像データを再構成するには、被検体の周囲１周、約３６０°分の投影データが、またハーフスキャン法でも１８０°＋ビュー角分の投影データが必要とされる。いずれの再構成方式にも本発明を適用可能である。ここでは、ハーフスキャン法を例に説明する。また、入射Ｘ線を電荷に変換するメカニズムは、シンチレータ等の蛍光体でＸ線を光に変換し更にその光をフォトダイオード等の光電変換素子で電荷に変換する間接変換形と、Ｘ線による半導体内の電子正孔対の生成及びその電極への移動すなわち光導電現象を利用した直接変換形とが主流である。Ｘ線検出素子としては、それらのいずれの方式を採用してもよいが、ここでは、前者の間接変換形として説明する。また、近年では、Ｘ線管とＸ線検出器との複数のペアを回転リングに搭載したいわゆる多管球型のＸ線ＣＴ装置の製品化が進み、その周辺技術の開発が進んでいる。本発明では、従来からの一管球型のＸ線ＣＴ装置であっても、多管球型のＸ線ＣＴ装置であってもいずれにも適用可能である。ここでは、一管球型として説明する。
【００１１】
また、本実施形態は、被検体の特定臓器を含むボリュームデータから特定臓器の運動機能に関する指標を計算し、その指標を極座標系に分布することによりポーラーマップを生成し、それと共に、同じボリュームデータから特定臓器又はその周辺臓器の形態を極座標系で表現したポーラーモデルを生成し、これらポーラーモデルとポーラーマップとを合成して表示するものであり、この実施形態で最も効果的なのは、心臓、特に左心室の心筋を特定臓器に適用し、またポーラーモデルとして形態を表現する対象に冠状動脈を適用した場合である。以下の説明では、この最も効果的な左心室の心筋に関するポーラーマップを生成し、冠状動脈のポーラーモデルを生成する場合を例に説明する。しかし、特定臓器は心臓、左心室の心筋には限定されることはなく、同様に、その周辺臓器は、冠状動脈に限定されるものではない。
【００１２】
図１に、本実施形態に係るコンピュータ断層撮影装置の構成をブロック図により示している。架台部１００は、架台駆動部１０７により回転駆動される回転リング１０２を有する。回転リング１０２には、コーンビーム形Ｘ線管１０１とマルチスライス形又は２次元アレイ形のＸ線検出器１０３とが、撮影時に被検体が挿入される空洞Ｓを隔てて、対向して搭載されている。Ｘ線管４は、高電圧発生装置１０９から連続的又は周期的に発生される高電圧をスリップリング１０８を経由して印加される。それにより、Ｘ線が四角錐形に放射される。
【００１３】
Ｘ線検出器１０３は、同時に複数スライス分の投影データを検出できるように、Ｘ線管１０１の焦点（コーンビームの頂点）を中心として円弧状に形成された多チャンネル型の検出素子列が、回転リング１０２の回転軸と略平行な方向（スライス方向）に沿って複数器並列されてなり、または平面的又は球部分形にコーンビームの頂点を中心として複数のＸ線検出素子がマトリクス状に配列されている。
【００１４】
Ｘ線検出器１０３の出力には、一般的にＤＡＳ(data acquisition system) と呼ばれているデータ収集回路１０４が接続されている。このデータ収集回路１０４には、Ｘ線検出器１０３の各チャンネルの電流信号を電圧に変換するＩ−Ｖ変換器と、この電圧信号をＸ線の曝射周期に同期して周期的に積分する積分器と、この積分器の出力信号を増幅するアンプと、このプリアンプの出力信号をディジタル信号に変換するアナログ・ディジタル・コンバータとが、チャンネルごとに設けられている。前処理装置１０６は、データ収集回路１０４の出力（ディジタル信号（純生データと呼ばれる））を光又は磁気を媒介して非接触のデータ伝送を実現する非接触データ伝送装置１０５を経由して受け取り、純生データに対して、チャンネル間の感度不均一を補正したり、またＸ線強吸収体、主に金属部による極端な信号強度の低下又は信号脱落を補正する等の前処理を実行する。前処理装置１０６から出力されるデータは、生データ又は投影データと呼ばれる。投影データは補助記憶装置１１２に記憶され、ホストコントローラ１１０の制御に従って再構成装置１１４に読み出される。再構成装置１１４は、補助記憶装置１１２からの投影データのセットに基づいてフェルドカンプ法又は他の再構成法に従って断層又はボリュームに関する画像データを再構成する。この画像データは補助記憶装置１１２に記憶され、また表示装置１１６に可視画像として表示される。
【００１５】
ホストコントローラ１１０には、Ｘ線ＣＴ装置としての基本的な構成要素としての上記補助記憶装置１１２、再構成装置１１４、表示装置１１６、キーボードやポインティングデバイス等を有する入力装置１１５の他に、本実施形態で特徴的な画像処理装置を構成する及び画像処理方法を実現するために、断面変換処理部１１７、心筋／冠状動脈抽出処理部１１８、心機能指標計算部１１９、心機能指標の極座標分布、つまりポーラーマップを生成するポーラーマップ生成部１２０、冠状動脈の構造（形態）を極座標で表現した画像（ポーラーモデル）を生成するポーラーモデル生成部１２１、ポーラーマップにポーラーモデルを合成して表示画像データを生成する表示画像生成部１２２、さらにボリュームレンダリングを始め各種の３次元処理を実現する３Ｄ処理部１２３を有している。以下、これら構成要素の詳細について処理手順に従って説明する。
【００１６】
図２は、ポーラーマップ生成手順を示している。まず、造影剤の静脈投与後に、被検体の心臓を含む胸部領域が、少なくとも拡張末期と収縮末期とを含む心拍位相にスキャン（ボリュームスキャン）が繰り返される（Ｓ１）。それにより各心拍位相の投影データが収集される。この投影データに基づいて拡張末期に対応するボリュームデータ、及び収縮末期に対応するボリュームデータが再構成される（Ｓ２）。
【００１７】
次に、断面変換処理部１１７は、図４に示すように、入力装置１１５を介して手動で又は拡張若しくは収縮末期のボリュームデータから自動認識した左心室の心軸に沿って、その心軸に直交する複数のスライスを、入力装置１１５を介して設定されたスライス枚数（後述の同心円数、ここでは４０と仮定する）に応じた一定の間隔で設定し、それら各スライスについて断層像データ（短軸像と呼ばれる）を、拡張末期のボリュームデータから生成する（Ｓ３）。同様に、拡張末期のそれと同じ位置のスライス各々に対応する断層像データを、収縮末期のボリュームデータから生成する。
【００１８】
心筋／冠状動脈抽出処理部１１８は、ＣＴ値に対するしきい値処理（階調処理）等により、拡張末期の各スライスの断層像データから、図５に示すように、左心室の心筋領域を抽出する（Ｓ４）。同様に、収縮末期の各スライスの断層像データから、左心室の心筋領域を抽出する。
【００１９】
なお、造影剤投与後の特定の時刻には、造影効果は、左心室の心筋で強く発揮される。従って、特定の時刻に収集したボリュームデータを使って、しきい値処理により、左心室の心筋領域を高精度で抽出することができる。上記説明では、ボリュームデータから各スライスについて断層像データを生成し（Ｓ３）、各断層像データから左心室の心筋領域を抽出（Ｓ４）したが、まず、ボリュームデータから左心室の心筋領域を抽出し、その抽出した左心室の心筋領域から各スライス像を生成するようにしてもよい。この場合、左心室の心軸の特定が簡易化される可能性がある。
【００２０】
心機能指標計算部１１９は、拡張末期と収縮末期とで、同じスライスどうしの心筋領域の変動に基づいて、心臓の収縮機能を定量化する心機能指標を、図６に示すように、心軸を中心として放射状に、入力装置１１５を介して設定された各同心円の分割区域（セグメント）の数に一致する複数、ここでは７２の方向に関して個々に計算する（Ｓ５）。心機能指標には、例えば、
・拡張末期の左心室の心軸から、左心室の心筋外壁までの距離から、収縮末期の左心室心軸から左心室心筋外壁までの距離を引き算して得られる容積半径変化(Wall Motion)、
・収縮末期の左心室心筋の壁厚から拡張末期の左心室心筋の壁厚を引き算した厚さの差を、拡張末期の左心室心筋の壁厚で割り算（正規化）して得られる左心室心筋壁厚変化率(W all Thickening)、
・拡張末期の左心室心筋内径の二乗から収縮末期の左心室心筋内径の二乗を引き算して得られる容積変化を、収縮末期の左心室心筋内径の二乗で割り算して得られる容積変化率(Rerional EF)
等があり、任意の種類の指標が入力装置１１５を介して選択される。
【００２１】
ポーラーマップ生成部１２０は、心機能指標計算部１１９で計算された心機能指標を、ポーラーマップテンプレート上において、対応するスライス番号に応じた同心円であって、その中の方向（セグメント番号）に応じたセグメントに、プロットすることによりポーラーマップデータを生成する（Ｓ６）。以上の処理により、心機能指標を極座標系で表現したポーラーマップデータが生成される。
【００２２】
次に、ポーラーマップデータＴ同じ被検体に関する実際の又は実際に近似した冠状動脈の構造を極座標系で表現したポーラーモデルデータの生成について説明する。図３にその手順を示している。まず、ポーラーモデルで冠状動脈構造を表現する心拍位相が、入力装置１１５を介して設定される（Ｓ７）。本実施形態では、冠状動脈の位置計算方法として、自動と手動（半自動）とを提供する。まず、自動方法について説明する。
【００２３】
まず、Ｓ７で設定された心拍位相に対応する上記Ｓ３で既生成された４０スライスの断層像データセットが心筋／冠状動脈抽出処理部１１８に供給される。心筋／冠状動脈抽出処理部１１８は、図８に示すように、Ｓ４で抽出済みの心筋領域の外側例えば１ｃｍの範囲を探索範囲として、各スライスの断層像データから冠状動脈領域を抽出する（Ｓ９）。
【００２４】
次に、ポーラーモデル生成部１２１は、図９に示すように、抽出された冠状動脈領域の位置（ここでは心軸回りの角度）を、スライスごとに計算する（Ｓ１０）。次に、３Ｄ処理部１２３は、計算した離散点を、図１５に示すように、心尖部側から心基部に向かって及びその逆向きの両方向で、最も近い離散点どうしを直線で連結することにより、冠状動脈の３Ｄ構造モデルを生成する（Ｓ１８）。ポーラーモデル生成部１２１は、生成された冠状動脈の３Ｄ構造モデルに対して、図１６に示すように、短軸像と同じ間隔で同じ枚数（４０枚）のスライスを設定し、各スライスごとに冠状動脈領域の位置（心軸回りの角度）を計算する（Ｓ１９）。
【００２５】
さらに図１０に示すように、ポーラーマップテンプレートと同じ様式のポーラーモデルテンプレートのスライス番号に対応する同心円の外縁上であって、角度に対応する位置に、点（動脈点という）をプロットする（Ｓ１１）。図１１に全スライスの動脈点をプロットした状態のポーラーモデルテンプレートを示している。
【００２６】
ポーラーモデル生成部１２１は、図１２に示すように、円中心（心尖部）側から外側（心基部）に向かって各同心円を挟んで最も近い動脈点どうしを直線で連結し、さらに、ポーラーモデル生成部１２１は、今度は逆向きに、外側（心基部）から円中心（心尖部）側に向かって各同心円を挟んで最も近い動脈点どうしを直線で連結することにより、冠状動脈のポーラーモデルデータを生成する（Ｓ１２）。このように両方向から動脈点を連結することにより、冠状動脈の分岐に対応することができる。この連結処理は、冠状動脈の解剖学的な知識を利用して自動又は手動で行うことが好ましい。
【００２７】
表示画像生成部１２２は、生成された冠状動脈のポーラーモデルを、Ｓ６で生成した心機能指標のポーラーマップに円角度を合わせて合成する（Ｓ１３）。冠状動脈のポーラーモデルが合成された心機能指標のポーラーマップを、図１３に示すように表示する（Ｓ１４）。
【００２８】
次に図３に戻り手動方法について説明する。まず、Ｓ１５において、３Ｄ処理部１２３は、拡張末期又は収縮末期のボリュームデータから、またはＳ３で生成した４０スライスの断層像データ（短軸像データ）から、ボリュームレンダリング処理、ＭＩＰ処理（最大値投影処理）により、図１４に示すように、心臓の３次元画像を生成する。心臓の３次元画像が、表示装置１１６の画面に、図１４に示すように、表示される（Ｓ１６）。操作者は、入力装置１１５のポインティングデバイスを介して、図１４に点線で示すように、３次元画像上に表示された冠状動脈の像上の複数点を離散的に指定する。その離散点の位置が３Ｄ処理部１２３で計算される（Ｓ１７）。次に、３Ｄ処理部１２３は、計算した離散点を、図１５に示すように、心尖部側から心基部に向かって及びその逆向きの両方向で、最も近い離散点どうしを直線で連結することにより、冠状動脈の３Ｄ構造モデルを生成する（Ｓ１８）。
【００２９】
ポーラーモデル生成部１２１は、生成された冠状動脈の３Ｄ構造モデルに対して、図１６に示すように、短軸像と同じ間隔で同じ枚数（４０枚）のスライスを設定し、各スライスごとに冠状動脈領域の位置（心軸回りの角度）を計算する（Ｓ１９）。以後、Ｓ１１乃至Ｓ１４は上述した通りである。
【００３０】
以上のように、心機能指標を極座標系で表現したポーラーマップに、その実際の被検体に関する冠状動脈の構造を表すポーラーモデルを重ねて表示することにより、観察者は、心臓の心機能指標を極座標で表現するポーラーマップで、障害が起きている冠状動脈の支配区域の判別だけでなく、その具体的な冠状動脈の障害位置まで特定又は推定することができる。
【００３１】
なお、ポーラーマップ生成部１２０は、４０の同心円×７２（セグメント（基本セグメントという））で生成されたポーラーマップを、図１７に示すような解剖学的に意義のあるセグメンテーションに再配置する機能を備えている。図１７に示す各セグメントには、それぞれのセグメントに含まれる複数の基本セグメントが含まれ、各セグメントの値を、そのセグメントに含まれる複数の基本セグメントの心機能指標値の平均値、最大値、最小値、最大頻度値をそのセグメントの特徴値として計算し、割り当てる。それにより様々な臨床診断上意義のあるポーラーマップを再生成することができる。もちろん、同じ極座標表現であるので、再生成されたポーラーマップにも冠状動脈のポーラーモデルを合成することができる。
【００３２】
また、３Ｄ処理部１２３は、図１８に示すように、入力装置１１５を介してポーラーマップ上で指定された冠状動脈モデル上の任意の位置に対応するＣｕｒｖｅｄＭＰＲ像を生成する機能を備えている。このＣｕｒｖｅｄＭＰＲ像は、心軸の直交する線が、指定点を含む冠状動脈に沿って移動することにより定義される局面に関する局面断層画像であり、任意の心拍位相に関するボリュームデータから再構成される。このＣｕｒｖｅｄＭＰＲ像により、心筋を含む心臓の動脈に沿った長軸像様の画像を観察することができる。例えば障害の可能性のある個所を指定することによりその周辺の心筋等の組織構造を縦に観察することができ、障害個所を特定することができる。
【００３３】
（変形例）
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
【００３４】
【発明の効果】
本発明によれば、観察者は、心臓の心機能指標を極座標で表現するポーラーマップで、障害が起きている冠状動脈の支配区域の判別だけでなく、その具体的な冠状動脈の障害位置まで特定又は推定することができる。
【図面の簡単な説明】
【図１】本発明の好ましい実施形態によるコンピュータ断層撮影装置の構成図。
【図２】本実施形態において、心機能指標のポーラーマップ作成手順を示すフローチャート。
【図３】本実施形態において、冠状動脈構造のポーラーモデル作成手順を示すフローチャート。
【図４】図２の断面変換処理Ｓ３の補足図。
【図５】図２の断面変換処理Ｓ３により生成された短軸像の例を示す図。
【図６】図２の心機能指標計算処理Ｓ５の補足図。
【図７】図２の心機能指標プロット処理Ｓ６の補足図。
【図８】図３の冠状動脈領域抽出処理Ｓ９の補足図。
【図９】図３の冠状動脈領域位置計算処理Ｓ１０の補足図。
【図１０】図３の動脈点プロット処理Ｓ１１の第１の補足図。
【図１１】図３の動脈点プロット処理Ｓ１１の第２の補足図。
【図１２】図３の動脈点連結処理Ｓ１２の補足図。
【図１３】本実施形態において、ポーラーモデルが合成されたポーラーマップの実際の表示例を示す中間調画像。
【図１４】図３の冠状動脈指定処理Ｓ１７の補足図。
【図１５】図３の冠状動脈構造モデル生成処理Ｓ１８の補足図。
【図１６】図３の冠状動脈位置計算処理Ｓ１９の補足図。
【図１７】本実施形態において、リセグメンテーションのテンプレートの例を示す図。
【図１８】本実施形態において、指定血管に対応するＣｕｒｖｅｄＭＰＲ像の実際の表示例を示す中間調画像。
【図１９】従来において、ポーラーマップに重ねられる支配領域区分図を示す図。
【符号の説明】
１０１…Ｘ線管、
１０２…回転リング、
１０３…２次元検出器、
１０４…ＤＡＳ、
１０５…非接触データ電装装置、
１０６…前処理装置、
１０７…架台駆動部、
１０８…スリップリング、
１０９…高電圧発生装置、
１１０…ホストコントローラ、
１１２…補助記憶装置、
１１４…再構成装置、
１１５…入力装置、
１１６…表示装置、
１１７…断面変換処理部、
１１８…心筋／冠状動脈抽出処理部、
１１９…心機能指標計算部、
１２０…ポーラーマップ生成部、
１２１…ポーラーモデル生成部、
１２２…表示画像生成部、
１２３…３Ｄ処理部。
【発明の属する技術分野】
本発明は、心臓の心機能指標を極座標で表現するポーラーマップの生成及び表示機能を備えたＸ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法に関する。
【０００２】
【従来の技術】
ヘリカルスキャンの高速化及び２次元検出器を使ったボリュームスキャンの開発は、Ｘ線コンピュータ断層撮影装置の概念を、単なる断面画像化装置からボリューム映像化化装置へ大きく移行させた。つまり、これらスキャン形式を備えたＸ線コンピュータ断層撮影装置により、被検体のボリュームを短い時間分解能で高速にスキャンすることが可能になった。このような機能を備えたＸ線コンピュータ断層撮影装置の実用化と共に、これまで使われていなかった多くの診断分野への進出が試みられ、また全く新しい診断方法を模索する動きが始まっている。
【０００３】
その一つの方向性として、心臓の心拍運動に関する評価診断への適用がある。ボリュームスキャンは、セグメントスキャン方式とあいまって、心臓全体を短いサイクルで連続的にスキャンすることを実質的に実現している。それにより例えば拡張末期（ＥＤ）や収縮末期（ＥＳ）に限定した心臓の立体像を生成することができ、さらに、拡張末期（ＥＤ）から収縮末期（ＥＳ）までの心壁厚の変化率等の心機能指標を心臓全体にわたって求めることも可能となる。
【０００４】
心臓は心尖部から心基部まで縦に長い。この縦に長い心臓全体の心機能指標を一度に観察できるように、その表示形式として一般的にポーラーマップ(Polar Map)と呼ばれる極座標画像で表示されることが多い。周知のとおり、極座標は、平面を（ｒ，θ）で表す表現方法であり、ここでは、θを心軸周りの角度、ｒを心尖部から心基部に向かって各スライスまでの距離にそれぞれ割り当てている。
【０００５】
このポーラーマップにより心機能低下を起こしている部分が容易に観察できるものであるが、実際に診断に用いる場合には、その部分と支配血管との関係が必要とされる。そのため、図１９に示すように、解剖学的に分割された既成の一般的な支配領域区分を重ねて表示することが多い。図１９において、“ＬＡＤ”表記の区分は左前下行枝の支配領域を表し、“ＬＣＸ”表記の区分は左回旋枝の支配領域を表し、“ＲＣＳ”表記の区分は右冠状動脈の支配領域を表している。
【０００６】
このように一般的な支配領域区分をポーラーマップに重ねて表示することにより、左前下行、枝左回旋枝、右冠状動脈のいずれの冠状動脈に障害が起きているか、つまりどの動脈の支配区域で障害が起きているかまでは把握することができるものの、その冠状動脈の果たして具体的にどの場所、どの位置で障害が起きているかは判断することはできない。
【０００７】
【発明が解決しようとする課題】
本発明の目的は、心臓の心機能指標を極座標で表現するポーラーマップの生成及び表示機能を備えたＸ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法において、障害が起きている冠状動脈の支配区域の判別だけでなく、その具体的な冠状動脈の障害位置まで特定又は推定できる情報の提供を実現することにある。
【０００８】
【課題を解決するための手段】
本発明の第１局面による画像処理装置は、被検体の特定臓器を含むボリュームデータから前記特定臓器の運動機能に関する指標を計算する指標計算部と、前記指標を極座標系に分布することによりポーラーマップを生成するポーラーマップ生成部と、前記ボリュームデータから前記特定臓器又はその周辺臓器の形態を、前記極座標系で表現したポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーモデルを前記ポーラーマップに合成する合成部と、前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備する。
本発明の第２局面によるＸ線コンピュータ断層撮影装置は、被検体の心臓を含む領域に関する投影データを収集する架台部と、前記投影データに基づいて画像データを再構成する再構成部と、前記画像データを、心軸に関して離散的な複数の短軸像のデータに変換する断面変換処理部と、前記短軸像各々から前記心臓の心筋領域を抽出する心筋抽出部と、前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算する心機能指標計算部と、前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成するポーラーマップ生成部と、前記短軸像各々から前記心臓の冠状動脈領域を抽出する冠状動脈抽出部と、前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーモデルを前記ポーラーマップに合成する合成部と、前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備する。
本発明の第３局面による画像処理装置は、被検体の心臓を含む領域に関する画像データを、心軸に関して離散的な複数の短軸像のデータに変換する断面変換処理部と、前記短軸像各々から前記心臓の心筋領域を抽出する心筋抽出部と、前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算する心機能指標計算部と、前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成するポーラーマップ生成部と、前記短軸像各々から前記心臓の冠状動脈領域を抽出する冠状動脈抽出部と、前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーモデルを前記ポーラーマップに合成する合成部と、前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備する。
本発明の第４局面による画像処理方法は、被検体の心臓を含む領域に関する画像データを、心軸に関して離散的な複数の短軸像のデータに変換し、前記短軸像各々から前記心臓の心筋領域を抽出し、前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算し、前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成し、前記短軸像各々から前記心臓の冠状動脈領域を抽出し、前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成し、前記ポーラーモデルを前記ポーラーマップに合成し、前記ポーラーモデルが合成されたポーラーマップを表示する。
本発明の第５局面によるＸ線コンピュータ断層撮影装置は、被検体の心臓を含む領域に関する投影データを収集する架台部と、前記投影データに基づいて画像データを再構成する再構成部と、前記画像データから前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成するポーラーマップ生成部と、前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーマップを前記ポーラーモデルとともに表示する表示部とを具備する。
本発明の第６局面による画像処理装置は、被検体の心臓を含む領域に関する画像データから、前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成するポーラーマップ生成部と、前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成するポーラーモデル生成部と、前記ポーラーマップを前記ポーラーモデルとともに表示する表示部とを具備する。 本発明の第７局面による画像処理方法は、被検体の心臓を含む領域に関する画像データから、前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成し、前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成し、前記ポーラーマップを前記ポーラーモデルとともに表示する。
【０００９】
【発明の実施の形態】
以下、図面を参照して本発明によるＸ線コンピュータ断層撮影装置（Ｘ線ＣＴ装置）、画像処理装置及び画像処理方法の実施形態を説明する。なお、以下では、本実施形態は、Ｘ線ＣＴ装置として説明する。しかし、以下のＸ線ＣＴ装置は、本実施形態の画像処理方法を実現する機能を有する画像処理装置を装備している。
【００１０】
ここで、Ｘ線ＣＴ装置には、Ｘ線管と放射線検出器とが１体として被検体の周囲を回転する回転／回転（ROTATE/ROTATE）タイプと、リング状に多数の検出素子がアレイされ、Ｘ線管のみが被検体の周囲を回転する固定／回転（STATIONARY/ROTATE）タイプ等様々なタイプがあり、いずれのタイプでも本発明を適用可能である。ここでは、現在、主流を占めている回転／回転タイプとして説明する。また、１スライスの断層像データを再構成するには、被検体の周囲１周、約３６０°分の投影データが、またハーフスキャン法でも１８０°＋ビュー角分の投影データが必要とされる。いずれの再構成方式にも本発明を適用可能である。ここでは、ハーフスキャン法を例に説明する。また、入射Ｘ線を電荷に変換するメカニズムは、シンチレータ等の蛍光体でＸ線を光に変換し更にその光をフォトダイオード等の光電変換素子で電荷に変換する間接変換形と、Ｘ線による半導体内の電子正孔対の生成及びその電極への移動すなわち光導電現象を利用した直接変換形とが主流である。Ｘ線検出素子としては、それらのいずれの方式を採用してもよいが、ここでは、前者の間接変換形として説明する。また、近年では、Ｘ線管とＸ線検出器との複数のペアを回転リングに搭載したいわゆる多管球型のＸ線ＣＴ装置の製品化が進み、その周辺技術の開発が進んでいる。本発明では、従来からの一管球型のＸ線ＣＴ装置であっても、多管球型のＸ線ＣＴ装置であってもいずれにも適用可能である。ここでは、一管球型として説明する。
【００１１】
また、本実施形態は、被検体の特定臓器を含むボリュームデータから特定臓器の運動機能に関する指標を計算し、その指標を極座標系に分布することによりポーラーマップを生成し、それと共に、同じボリュームデータから特定臓器又はその周辺臓器の形態を極座標系で表現したポーラーモデルを生成し、これらポーラーモデルとポーラーマップとを合成して表示するものであり、この実施形態で最も効果的なのは、心臓、特に左心室の心筋を特定臓器に適用し、またポーラーモデルとして形態を表現する対象に冠状動脈を適用した場合である。以下の説明では、この最も効果的な左心室の心筋に関するポーラーマップを生成し、冠状動脈のポーラーモデルを生成する場合を例に説明する。しかし、特定臓器は心臓、左心室の心筋には限定されることはなく、同様に、その周辺臓器は、冠状動脈に限定されるものではない。
【００１２】
図１に、本実施形態に係るコンピュータ断層撮影装置の構成をブロック図により示している。架台部１００は、架台駆動部１０７により回転駆動される回転リング１０２を有する。回転リング１０２には、コーンビーム形Ｘ線管１０１とマルチスライス形又は２次元アレイ形のＸ線検出器１０３とが、撮影時に被検体が挿入される空洞Ｓを隔てて、対向して搭載されている。Ｘ線管４は、高電圧発生装置１０９から連続的又は周期的に発生される高電圧をスリップリング１０８を経由して印加される。それにより、Ｘ線が四角錐形に放射される。
【００１３】
Ｘ線検出器１０３は、同時に複数スライス分の投影データを検出できるように、Ｘ線管１０１の焦点（コーンビームの頂点）を中心として円弧状に形成された多チャンネル型の検出素子列が、回転リング１０２の回転軸と略平行な方向（スライス方向）に沿って複数器並列されてなり、または平面的又は球部分形にコーンビームの頂点を中心として複数のＸ線検出素子がマトリクス状に配列されている。
【００１４】
Ｘ線検出器１０３の出力には、一般的にＤＡＳ(data acquisition system) と呼ばれているデータ収集回路１０４が接続されている。このデータ収集回路１０４には、Ｘ線検出器１０３の各チャンネルの電流信号を電圧に変換するＩ−Ｖ変換器と、この電圧信号をＸ線の曝射周期に同期して周期的に積分する積分器と、この積分器の出力信号を増幅するアンプと、このプリアンプの出力信号をディジタル信号に変換するアナログ・ディジタル・コンバータとが、チャンネルごとに設けられている。前処理装置１０６は、データ収集回路１０４の出力（ディジタル信号（純生データと呼ばれる））を光又は磁気を媒介して非接触のデータ伝送を実現する非接触データ伝送装置１０５を経由して受け取り、純生データに対して、チャンネル間の感度不均一を補正したり、またＸ線強吸収体、主に金属部による極端な信号強度の低下又は信号脱落を補正する等の前処理を実行する。前処理装置１０６から出力されるデータは、生データ又は投影データと呼ばれる。投影データは補助記憶装置１１２に記憶され、ホストコントローラ１１０の制御に従って再構成装置１１４に読み出される。再構成装置１１４は、補助記憶装置１１２からの投影データのセットに基づいてフェルドカンプ法又は他の再構成法に従って断層又はボリュームに関する画像データを再構成する。この画像データは補助記憶装置１１２に記憶され、また表示装置１１６に可視画像として表示される。
【００１５】
ホストコントローラ１１０には、Ｘ線ＣＴ装置としての基本的な構成要素としての上記補助記憶装置１１２、再構成装置１１４、表示装置１１６、キーボードやポインティングデバイス等を有する入力装置１１５の他に、本実施形態で特徴的な画像処理装置を構成する及び画像処理方法を実現するために、断面変換処理部１１７、心筋／冠状動脈抽出処理部１１８、心機能指標計算部１１９、心機能指標の極座標分布、つまりポーラーマップを生成するポーラーマップ生成部１２０、冠状動脈の構造（形態）を極座標で表現した画像（ポーラーモデル）を生成するポーラーモデル生成部１２１、ポーラーマップにポーラーモデルを合成して表示画像データを生成する表示画像生成部１２２、さらにボリュームレンダリングを始め各種の３次元処理を実現する３Ｄ処理部１２３を有している。以下、これら構成要素の詳細について処理手順に従って説明する。
【００１６】
図２は、ポーラーマップ生成手順を示している。まず、造影剤の静脈投与後に、被検体の心臓を含む胸部領域が、少なくとも拡張末期と収縮末期とを含む心拍位相にスキャン（ボリュームスキャン）が繰り返される（Ｓ１）。それにより各心拍位相の投影データが収集される。この投影データに基づいて拡張末期に対応するボリュームデータ、及び収縮末期に対応するボリュームデータが再構成される（Ｓ２）。
【００１７】
次に、断面変換処理部１１７は、図４に示すように、入力装置１１５を介して手動で又は拡張若しくは収縮末期のボリュームデータから自動認識した左心室の心軸に沿って、その心軸に直交する複数のスライスを、入力装置１１５を介して設定されたスライス枚数（後述の同心円数、ここでは４０と仮定する）に応じた一定の間隔で設定し、それら各スライスについて断層像データ（短軸像と呼ばれる）を、拡張末期のボリュームデータから生成する（Ｓ３）。同様に、拡張末期のそれと同じ位置のスライス各々に対応する断層像データを、収縮末期のボリュームデータから生成する。
【００１８】
心筋／冠状動脈抽出処理部１１８は、ＣＴ値に対するしきい値処理（階調処理）等により、拡張末期の各スライスの断層像データから、図５に示すように、左心室の心筋領域を抽出する（Ｓ４）。同様に、収縮末期の各スライスの断層像データから、左心室の心筋領域を抽出する。
【００１９】
なお、造影剤投与後の特定の時刻には、造影効果は、左心室の心筋で強く発揮される。従って、特定の時刻に収集したボリュームデータを使って、しきい値処理により、左心室の心筋領域を高精度で抽出することができる。上記説明では、ボリュームデータから各スライスについて断層像データを生成し（Ｓ３）、各断層像データから左心室の心筋領域を抽出（Ｓ４）したが、まず、ボリュームデータから左心室の心筋領域を抽出し、その抽出した左心室の心筋領域から各スライス像を生成するようにしてもよい。この場合、左心室の心軸の特定が簡易化される可能性がある。
【００２０】
心機能指標計算部１１９は、拡張末期と収縮末期とで、同じスライスどうしの心筋領域の変動に基づいて、心臓の収縮機能を定量化する心機能指標を、図６に示すように、心軸を中心として放射状に、入力装置１１５を介して設定された各同心円の分割区域（セグメント）の数に一致する複数、ここでは７２の方向に関して個々に計算する（Ｓ５）。心機能指標には、例えば、
・拡張末期の左心室の心軸から、左心室の心筋外壁までの距離から、収縮末期の左心室心軸から左心室心筋外壁までの距離を引き算して得られる容積半径変化(Wall Motion)、
・収縮末期の左心室心筋の壁厚から拡張末期の左心室心筋の壁厚を引き算した厚さの差を、拡張末期の左心室心筋の壁厚で割り算（正規化）して得られる左心室心筋壁厚変化率(W all Thickening)、
・拡張末期の左心室心筋内径の二乗から収縮末期の左心室心筋内径の二乗を引き算して得られる容積変化を、収縮末期の左心室心筋内径の二乗で割り算して得られる容積変化率(Rerional EF)
等があり、任意の種類の指標が入力装置１１５を介して選択される。
【００２１】
ポーラーマップ生成部１２０は、心機能指標計算部１１９で計算された心機能指標を、ポーラーマップテンプレート上において、対応するスライス番号に応じた同心円であって、その中の方向（セグメント番号）に応じたセグメントに、プロットすることによりポーラーマップデータを生成する（Ｓ６）。以上の処理により、心機能指標を極座標系で表現したポーラーマップデータが生成される。
【００２２】
次に、ポーラーマップデータＴ同じ被検体に関する実際の又は実際に近似した冠状動脈の構造を極座標系で表現したポーラーモデルデータの生成について説明する。図３にその手順を示している。まず、ポーラーモデルで冠状動脈構造を表現する心拍位相が、入力装置１１５を介して設定される（Ｓ７）。本実施形態では、冠状動脈の位置計算方法として、自動と手動（半自動）とを提供する。まず、自動方法について説明する。
【００２３】
まず、Ｓ７で設定された心拍位相に対応する上記Ｓ３で既生成された４０スライスの断層像データセットが心筋／冠状動脈抽出処理部１１８に供給される。心筋／冠状動脈抽出処理部１１８は、図８に示すように、Ｓ４で抽出済みの心筋領域の外側例えば１ｃｍの範囲を探索範囲として、各スライスの断層像データから冠状動脈領域を抽出する（Ｓ９）。
【００２４】
次に、ポーラーモデル生成部１２１は、図９に示すように、抽出された冠状動脈領域の位置（ここでは心軸回りの角度）を、スライスごとに計算する（Ｓ１０）。次に、３Ｄ処理部１２３は、計算した離散点を、図１５に示すように、心尖部側から心基部に向かって及びその逆向きの両方向で、最も近い離散点どうしを直線で連結することにより、冠状動脈の３Ｄ構造モデルを生成する（Ｓ１８）。ポーラーモデル生成部１２１は、生成された冠状動脈の３Ｄ構造モデルに対して、図１６に示すように、短軸像と同じ間隔で同じ枚数（４０枚）のスライスを設定し、各スライスごとに冠状動脈領域の位置（心軸回りの角度）を計算する（Ｓ１９）。
【００２５】
さらに図１０に示すように、ポーラーマップテンプレートと同じ様式のポーラーモデルテンプレートのスライス番号に対応する同心円の外縁上であって、角度に対応する位置に、点（動脈点という）をプロットする（Ｓ１１）。図１１に全スライスの動脈点をプロットした状態のポーラーモデルテンプレートを示している。
【００２６】
ポーラーモデル生成部１２１は、図１２に示すように、円中心（心尖部）側から外側（心基部）に向かって各同心円を挟んで最も近い動脈点どうしを直線で連結し、さらに、ポーラーモデル生成部１２１は、今度は逆向きに、外側（心基部）から円中心（心尖部）側に向かって各同心円を挟んで最も近い動脈点どうしを直線で連結することにより、冠状動脈のポーラーモデルデータを生成する（Ｓ１２）。このように両方向から動脈点を連結することにより、冠状動脈の分岐に対応することができる。この連結処理は、冠状動脈の解剖学的な知識を利用して自動又は手動で行うことが好ましい。
【００２７】
表示画像生成部１２２は、生成された冠状動脈のポーラーモデルを、Ｓ６で生成した心機能指標のポーラーマップに円角度を合わせて合成する（Ｓ１３）。冠状動脈のポーラーモデルが合成された心機能指標のポーラーマップを、図１３に示すように表示する（Ｓ１４）。
【００２８】
次に図３に戻り手動方法について説明する。まず、Ｓ１５において、３Ｄ処理部１２３は、拡張末期又は収縮末期のボリュームデータから、またはＳ３で生成した４０スライスの断層像データ（短軸像データ）から、ボリュームレンダリング処理、ＭＩＰ処理（最大値投影処理）により、図１４に示すように、心臓の３次元画像を生成する。心臓の３次元画像が、表示装置１１６の画面に、図１４に示すように、表示される（Ｓ１６）。操作者は、入力装置１１５のポインティングデバイスを介して、図１４に点線で示すように、３次元画像上に表示された冠状動脈の像上の複数点を離散的に指定する。その離散点の位置が３Ｄ処理部１２３で計算される（Ｓ１７）。次に、３Ｄ処理部１２３は、計算した離散点を、図１５に示すように、心尖部側から心基部に向かって及びその逆向きの両方向で、最も近い離散点どうしを直線で連結することにより、冠状動脈の３Ｄ構造モデルを生成する（Ｓ１８）。
【００２９】
ポーラーモデル生成部１２１は、生成された冠状動脈の３Ｄ構造モデルに対して、図１６に示すように、短軸像と同じ間隔で同じ枚数（４０枚）のスライスを設定し、各スライスごとに冠状動脈領域の位置（心軸回りの角度）を計算する（Ｓ１９）。以後、Ｓ１１乃至Ｓ１４は上述した通りである。
【００３０】
以上のように、心機能指標を極座標系で表現したポーラーマップに、その実際の被検体に関する冠状動脈の構造を表すポーラーモデルを重ねて表示することにより、観察者は、心臓の心機能指標を極座標で表現するポーラーマップで、障害が起きている冠状動脈の支配区域の判別だけでなく、その具体的な冠状動脈の障害位置まで特定又は推定することができる。
【００３１】
なお、ポーラーマップ生成部１２０は、４０の同心円×７２（セグメント（基本セグメントという））で生成されたポーラーマップを、図１７に示すような解剖学的に意義のあるセグメンテーションに再配置する機能を備えている。図１７に示す各セグメントには、それぞれのセグメントに含まれる複数の基本セグメントが含まれ、各セグメントの値を、そのセグメントに含まれる複数の基本セグメントの心機能指標値の平均値、最大値、最小値、最大頻度値をそのセグメントの特徴値として計算し、割り当てる。それにより様々な臨床診断上意義のあるポーラーマップを再生成することができる。もちろん、同じ極座標表現であるので、再生成されたポーラーマップにも冠状動脈のポーラーモデルを合成することができる。
【００３２】
また、３Ｄ処理部１２３は、図１８に示すように、入力装置１１５を介してポーラーマップ上で指定された冠状動脈モデル上の任意の位置に対応するＣｕｒｖｅｄＭＰＲ像を生成する機能を備えている。このＣｕｒｖｅｄＭＰＲ像は、心軸の直交する線が、指定点を含む冠状動脈に沿って移動することにより定義される局面に関する局面断層画像であり、任意の心拍位相に関するボリュームデータから再構成される。このＣｕｒｖｅｄＭＰＲ像により、心筋を含む心臓の動脈に沿った長軸像様の画像を観察することができる。例えば障害の可能性のある個所を指定することによりその周辺の心筋等の組織構造を縦に観察することができ、障害個所を特定することができる。
【００３３】
（変形例）
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
【００３４】
【発明の効果】
本発明によれば、観察者は、心臓の心機能指標を極座標で表現するポーラーマップで、障害が起きている冠状動脈の支配区域の判別だけでなく、その具体的な冠状動脈の障害位置まで特定又は推定することができる。
【図面の簡単な説明】
【図１】本発明の好ましい実施形態によるコンピュータ断層撮影装置の構成図。
【図２】本実施形態において、心機能指標のポーラーマップ作成手順を示すフローチャート。
【図３】本実施形態において、冠状動脈構造のポーラーモデル作成手順を示すフローチャート。
【図４】図２の断面変換処理Ｓ３の補足図。
【図５】図２の断面変換処理Ｓ３により生成された短軸像の例を示す図。
【図６】図２の心機能指標計算処理Ｓ５の補足図。
【図７】図２の心機能指標プロット処理Ｓ６の補足図。
【図８】図３の冠状動脈領域抽出処理Ｓ９の補足図。
【図９】図３の冠状動脈領域位置計算処理Ｓ１０の補足図。
【図１０】図３の動脈点プロット処理Ｓ１１の第１の補足図。
【図１１】図３の動脈点プロット処理Ｓ１１の第２の補足図。
【図１２】図３の動脈点連結処理Ｓ１２の補足図。
【図１３】本実施形態において、ポーラーモデルが合成されたポーラーマップの実際の表示例を示す中間調画像。
【図１４】図３の冠状動脈指定処理Ｓ１７の補足図。
【図１５】図３の冠状動脈構造モデル生成処理Ｓ１８の補足図。
【図１６】図３の冠状動脈位置計算処理Ｓ１９の補足図。
【図１７】本実施形態において、リセグメンテーションのテンプレートの例を示す図。
【図１８】本実施形態において、指定血管に対応するＣｕｒｖｅｄＭＰＲ像の実際の表示例を示す中間調画像。
【図１９】従来において、ポーラーマップに重ねられる支配領域区分図を示す図。
【符号の説明】
１０１…Ｘ線管、
１０２…回転リング、
１０３…２次元検出器、
１０４…ＤＡＳ、
１０５…非接触データ電装装置、
１０６…前処理装置、
１０７…架台駆動部、
１０８…スリップリング、
１０９…高電圧発生装置、
１１０…ホストコントローラ、
１１２…補助記憶装置、
１１４…再構成装置、
１１５…入力装置、
１１６…表示装置、
１１７…断面変換処理部、
１１８…心筋／冠状動脈抽出処理部、
１１９…心機能指標計算部、
１２０…ポーラーマップ生成部、
１２１…ポーラーモデル生成部、
１２２…表示画像生成部、
１２３…３Ｄ処理部。
Claims (8)
- 被検体の特定臓器を含むボリュームデータから前記特定臓器の運動機能に関する指標を計算する指標計算部と、
前記指標を極座標系に分布することによりポーラーマップを生成するポーラーマップ生成部と、
前記ボリュームデータから前記特定臓器又はその周辺臓器の形態を、前記極座標系で表現したポーラーモデルを生成するポーラーモデル生成部と、
前記ポーラーモデルを前記ポーラーマップに合成する合成部と、
前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備することを特徴とする画像処理装置。 - 被検体の心臓を含む領域に関する投影データを収集する架台部と、
前記投影データに基づいて画像データを再構成する再構成部と、
前記画像データを、心軸に関して離散的な複数の短軸像のデータに変換する断面変換処理部と、
前記短軸像各々から前記心臓の心筋領域を抽出する心筋抽出部と、
前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算する心機能指標計算部と、
前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成するポーラーマップ生成部と、
前記短軸像各々から前記心臓の冠状動脈領域を抽出する冠状動脈抽出部と、
前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成するポーラーモデル生成部と、
前記ポーラーモデルを前記ポーラーマップに合成する合成部と、
前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備することを特徴とするＸ線コンピュータ断層撮影装置。 - 前記冠状動脈抽出部は、前記抽出した心筋領域の外縁から所定範囲を探索することを特徴とする請求項２記載のＸ線コンピュータ断層撮影装置。
- 被検体の心臓を含む領域に関する画像データを、心軸に関して離散的な複数の短軸像のデータに変換する断面変換処理部と、
前記短軸像各々から前記心臓の心筋領域を抽出する心筋抽出部と、
前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算する心機能指標計算部と、
前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成するポーラーマップ生成部と、
前記短軸像各々から前記心臓の冠状動脈領域を抽出する冠状動脈抽出部と、
前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成するポーラーモデル生成部と、
前記ポーラーモデルを前記ポーラーマップに合成する合成部と、
前記ポーラーモデルが合成されたポーラーマップを表示する表示部とを具備することを特徴とする画像処理装置。 - 被検体の心臓を含む領域に関する画像データを、心軸に関して離散的な複数の短軸像のデータに変換し、
前記短軸像各々から前記心臓の心筋領域を抽出し、
前記抽出した心筋領域から前記心臓の心拍運動に関わる心機能指標を心軸周りの多方向にわたって計算し、
前記計算された心機能指標を、心軸周りの角度と心尖部又は心基部からの距離とで定義される極座標系にプロットすることにより、ポーラーマップを生成し、
前記短軸像各々から前記心臓の冠状動脈領域を抽出し、
前記抽出した冠状動脈領域の位置を前記極座標系にプロットし、プロットした点どうしを連結することにより、冠状動脈のポーラーモデルを生成し、
前記ポーラーモデルを前記ポーラーマップに合成し、
前記ポーラーモデルが合成されたポーラーマップを表示することを特徴とする画像処理方法。 - 被検体の心臓を含む領域に関する投影データを収集する架台部と、
前記投影データに基づいて画像データを再構成する再構成部と、
前記画像データから前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成するポーラーマップ生成部と、
前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成するポーラーモデル生成部と、
前記ポーラーマップを前記ポーラーモデルとともに表示する表示部とを具備することを特徴とするＸ線コンピュータ断層撮影装置。 - 被検体の心臓を含む領域に関する画像データから、前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成するポーラーマップ生成部と、
前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成するポーラーモデル生成部と、
前記ポーラーマップを前記ポーラーモデルとともに表示する表示部とを具備することを特徴とする画像処理装置。 - 被検体の心臓を含む領域に関する画像データから、前記心臓の局所心機能指標を極座標で表現するポーラーマップを生成し、
前記画像データから前記心臓の冠状動脈の構造を極座標で表現するポーラーモデルを生成し、
前記ポーラーマップを前記ポーラーモデルとともに表示することを特徴とする画像処理方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002307168A JP4346297B2 (ja) | 2002-10-22 | 2002-10-22 | Ｘ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 |
| US10/685,489 US7260252B2 (en) | 2002-10-22 | 2003-10-16 | X-ray computed tomographic apparatus, image processing apparatus, and image processing method |
| CNB200310102708XA CN1292711C (zh) | 2002-10-22 | 2003-10-22 | X射线计算机层析装置、图像处理装置以及图像处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002307168A JP4346297B2 (ja) | 2002-10-22 | 2002-10-22 | Ｘ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004141245A JP2004141245A (ja) | 2004-05-20 |
| JP4346297B2 true JP4346297B2 (ja) | 2009-10-21 |
Family
ID=32170931
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002307168A Expired - Lifetime JP4346297B2 (ja) | 2002-10-22 | 2002-10-22 | Ｘ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7260252B2 (ja) |
| JP (1) | JP4346297B2 (ja) |
| CN (1) | CN1292711C (ja) |
Families Citing this family (42)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7574247B2 (en) * | 2003-11-17 | 2009-08-11 | Siemens Medical Solutions Usa, Inc. | Automatic coronary isolation using a n-MIP ray casting technique |
| US8010175B2 (en) * | 2004-05-05 | 2011-08-30 | Siemens Medical Solutions Usa, Inc. | Patient-specific coronary territory mapping |
| CN1989527B (zh) * | 2004-07-27 | 2011-10-05 | 皇家飞利浦电子股份有限公司 | 成像几何参数的自动确定 |
| JP5060719B2 (ja) * | 2004-11-01 | 2012-10-31 | 富士フイルムＲｉファーマ株式会社 | 生体臓器の画像処理装置、方法及びコンピュータプログラム |
| JP5060720B2 (ja) * | 2004-11-01 | 2012-10-31 | 富士フイルムＲｉファーマ株式会社 | 診断支援装置、方法、及びコンピュータプログラム |
| DE102005002949A1 (de) * | 2005-01-21 | 2006-08-03 | Siemens Ag | Verfahren zur Visualisierung von Schädigungen im Myokard |
| US20060247544A1 (en) * | 2005-02-03 | 2006-11-02 | Maleeha Qazi | Characterization of cardiac motion with spatial relationship |
| US7715627B2 (en) * | 2005-03-25 | 2010-05-11 | Siemens Medical Solutions Usa, Inc. | Automatic determination of the standard cardiac views from volumetric data acquisitions |
| US7813535B2 (en) * | 2005-04-19 | 2010-10-12 | Siemens Medical Solutions Usa, Inc. | System and method for fused PET-CT visualization for heart unfolding |
| DE102005018067A1 (de) * | 2005-04-19 | 2006-11-02 | Siemens Ag | System zur Erzeugung, Auswertung und Verteilung computertomographischer 4D-Darstellungen des Herzens eines Patienten |
| US8285012B2 (en) * | 2005-04-28 | 2012-10-09 | Hitachi Medical Corporation | Image display apparatus and program |
| CN100562290C (zh) * | 2005-05-31 | 2009-11-25 | 柯尼卡美能达医疗印刷器材株式会社 | 图像处理方法以及图像处理装置 |
| CN100444798C (zh) * | 2005-07-22 | 2008-12-24 | 清华大学 | 一种用于ct重建的旋转对称体素离散化方法 |
| JP4777007B2 (ja) * | 2005-08-03 | 2011-09-21 | 東芝メディカルシステムズ株式会社 | Ｘ線コンピュータ断層撮影装置 |
| DE102005036998B4 (de) * | 2005-08-05 | 2014-11-20 | Siemens Aktiengesellschaft | Vorrichtung zur automatischen Detektion von Auffälligkeiten in medizinischen Bilddaten |
| DE102006026695A1 (de) * | 2006-06-08 | 2007-12-13 | Tomtec Imaging Systems Gmbh | Verfahren, Vorrichtung und Computerprogrammprodukt zum Auswerten von dynamischen Bildern einer Kavität |
| US20080008366A1 (en) * | 2006-06-20 | 2008-01-10 | Vladimir Desh | Simultaneous visualization, analysis and navigation of multi-modality medical imaging data |
| WO2008012755A2 (en) * | 2006-07-26 | 2008-01-31 | Koninklijke Philips Electronics N.V. | Visualization of stress level cardiac functional analysis results |
| US20080049999A1 (en) * | 2006-08-28 | 2008-02-28 | Siemens Medical Solutions Usa, Inc. | Computer Aided Detection of Bone Metastasis |
| US7978191B2 (en) * | 2007-09-24 | 2011-07-12 | Dolphin Imaging Systems, Llc | System and method for locating anatomies of interest in a 3D volume |
| US8571288B2 (en) * | 2007-12-07 | 2013-10-29 | Kabushiki Kaisha Toshiba | Image display apparatus and magnetic resonance imaging apparatus |
| US9064300B2 (en) * | 2008-02-15 | 2015-06-23 | Siemens Aktiengesellshaft | Method and system for automatic determination of coronory supply regions |
| JP5558672B2 (ja) * | 2008-03-19 | 2014-07-23 | 株式会社東芝 | 画像処理装置及びｘ線コンピュータ断層撮影装置 |
| WO2009147605A1 (en) * | 2008-06-04 | 2009-12-10 | Koninklijke Philips Electronics N.V. | Reconstruction of dynamical cardiac spect for measuring tracer uptake and redistribution |
| US8229192B2 (en) * | 2008-08-12 | 2012-07-24 | General Electric Company | Methods and apparatus to process left-ventricle cardiac images |
| EP2373218B1 (en) | 2008-12-03 | 2016-07-13 | Koninklijke Philips N.V. | Reparametrized bull's eye plots |
| JP5361410B2 (ja) * | 2009-01-22 | 2013-12-04 | 株式会社東芝 | 画像処理装置 |
| US8659603B2 (en) | 2009-03-31 | 2014-02-25 | General Electric Company | System and method for center point trajectory mapping |
| JP4926199B2 (ja) * | 2009-04-16 | 2012-05-09 | 富士フイルム株式会社 | 診断支援装置、診断支援プログラムおよび診断支援方法 |
| RU2584127C2 (ru) | 2010-03-23 | 2016-05-20 | Конинклейке Филипс Электроникс Н.В. | Данные ультразвукового объемного изображения, переформатированные в виде последовательности плоских изображений |
| JP5844093B2 (ja) * | 2010-09-15 | 2016-01-13 | 株式会社東芝 | 医用画像処理装置及び医用画像処理方法 |
| CN103764038A (zh) * | 2012-02-21 | 2014-04-30 | 株式会社东芝 | X射线ct装置、图像显示装置、图像显示方法 |
| WO2014084382A1 (ja) | 2012-11-30 | 2014-06-05 | 株式会社 東芝 | 医用画像処理装置 |
| EP2936364B1 (en) * | 2012-12-18 | 2019-09-25 | Koninklijke Philips N.V. | Method and apparatus for simulating blood flow under patient-specific boundary conditions derived from an estimated cardiac ejection output |
| US10371771B2 (en) * | 2014-07-30 | 2019-08-06 | Koninklijke Philips N.V. | Gradient amplifier system for driving a gradient coil and configuration method |
| US9949643B2 (en) * | 2014-10-18 | 2018-04-24 | International Business Machines Corporation | Automatic visualization of regional functional parameters of left ventricle from cardiac imaging |
| JP6734028B2 (ja) * | 2014-11-14 | 2020-08-05 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置、画像処理装置及び画像生成方法 |
| EP3245932A1 (en) * | 2015-07-06 | 2017-11-22 | Olympus Corporation | Medical device, medical image generating method, and medical image generating program |
| US10176604B2 (en) * | 2016-04-20 | 2019-01-08 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image reconstruction |
| JP6657132B2 (ja) * | 2017-02-27 | 2020-03-04 | 富士フイルム株式会社 | 画像分類装置、方法およびプログラム |
| CN108549530B (zh) * | 2018-04-17 | 2021-03-12 | 兰州理工大学 | 基于自然模式的Logo生成方法、智能Logo生成器 |
| EP3671649A1 (en) * | 2018-12-19 | 2020-06-24 | Siemens Healthcare GmbH | Method and computer system for generating a combined tissue-vessel representation |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60218166A (ja) | 1984-04-14 | 1985-10-31 | Fujitsu Ltd | 画像処理装置 |
| JPS6344103A (ja) | 1986-08-12 | 1988-02-25 | Nissan Motor Co Ltd | 位置姿勢認識装置 |
| JP2562047B2 (ja) | 1988-12-21 | 1996-12-11 | 株式会社不二越 | 対象物体の位置姿勢認識方法 |
| JPH0820207B2 (ja) | 1991-08-19 | 1996-03-04 | 川崎重工業株式会社 | 光学式３次元位置計測方法 |
| JPH06147830A (ja) | 1992-11-11 | 1994-05-27 | Mitsubishi Electric Corp | ３次元位置測定装置及び３次元位置測定結果補正方法 |
| JP3059602B2 (ja) | 1993-01-29 | 2000-07-04 | オリンパス光学工業株式会社 | 微小線幅測定装置及びその方法 |
| JP3604467B2 (ja) * | 1995-09-27 | 2004-12-22 | 株式会社東芝 | 心筋のねじれ補正方法 |
| US5970182A (en) * | 1995-11-15 | 1999-10-19 | Focus Imaging, S. A. | Registration process for myocardial images |
| US7194117B2 (en) * | 1999-06-29 | 2007-03-20 | The Research Foundation Of State University Of New York | System and method for performing a three-dimensional virtual examination of objects, such as internal organs |
| US5971767A (en) * | 1996-09-16 | 1999-10-26 | The Research Foundation Of State University Of New York | System and method for performing a three-dimensional virtual examination |
| JP2002109550A (ja) * | 2000-09-29 | 2002-04-12 | Fuji Photo Film Co Ltd | 心胸郭の輪郭検出方法および心胸郭比算出方法 |
| AU2002360478A1 (en) * | 2001-12-05 | 2003-06-17 | The Trustees Of The University Of Pennsylvania | Virtual bone biopsy |
| US6628743B1 (en) * | 2002-11-26 | 2003-09-30 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for acquiring and analyzing cardiac data from a patient |
-
2002
- 2002-10-22 JP JP2002307168A patent/JP4346297B2/ja not_active Expired - Lifetime
-
2003
- 2003-10-16 US US10/685,489 patent/US7260252B2/en active Active
- 2003-10-22 CN CNB200310102708XA patent/CN1292711C/zh not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| US20040087853A1 (en) | 2004-05-06 |
| JP2004141245A (ja) | 2004-05-20 |
| US7260252B2 (en) | 2007-08-21 |
| CN1541617A (zh) | 2004-11-03 |
| CN1292711C (zh) | 2007-01-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4346297B2 (ja) | Ｘ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 | |
| US6252924B1 (en) | Method and apparatus for motion-free cardiac CT imaging | |
| JP4777007B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| US6269140B1 (en) | X-ray computerized tomography apparatus control method therefor and image generating method using the apparatus | |
| US7180976B2 (en) | Rotational angiography based hybrid 3-D reconstruction of coronary arterial structure | |
| US7315605B2 (en) | Method and device for reconstructing a 3D image data set of a moving object | |
| JP5091644B2 (ja) | Ｘ線コンピュータ断層撮像装置及び医用画像処理装置 | |
| CN105934199B (zh) | 用于处理体腔的医学图像的设备和方法 | |
| JP2000083944A (ja) | 心臓用ｃｔシステム | |
| US8498462B2 (en) | Image processing apparatus and program | |
| US20100074490A1 (en) | Image processing apparatus and x-ray computer tomography apparatus | |
| US6718004B2 (en) | Methods and apparatus for coronary-specific imaging reconstruction | |
| JP4828839B2 (ja) | Ｘ線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 | |
| JP2004121840A (ja) | 周期的に運動する器官のｃｔ画像形成方法およびこの方法を実施するためのｃｔ装置 | |
| JP2007144139A (ja) | Ｘ線コンピュータ断層撮影装置及び画像処理装置 | |
| CN100424725C (zh) | 产生时间内插断层摄影图像的方法和装置 | |
| KR20170105876A (ko) | 단층 촬영 장치 및 그에 따른 단층 영상 재구성 방법 | |
| JP2005161032A (ja) | 画像処理装置 | |
| JP4309621B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| JP5624350B2 (ja) | 医用画像処理装置 | |
| JP5426783B2 (ja) | Ｘ線コンピュータ断層撮影装置、画像処理装置及びプログラム | |
| JP5455290B2 (ja) | 医用画像処理装置及び医用画像診断装置 | |
| JP2007068726A (ja) | 心機能診断装置 | |
| JP5984237B2 (ja) | 医用画像処理装置及びｘ線コンピュータ断層撮影装置 | |
| JP2003010173A (ja) | Ｘ線コンピュータ断層撮影装置及び画像処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051014 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20081024 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081111 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090623 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090714 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 4346297 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120724 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130724 Year of fee payment: 4 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |