JP6670560B2

JP6670560B2 - 画像処理方法、画像処理プログラム、画像処理装置及び画像処理システム

Info

Publication number: JP6670560B2
Application number: JP2015140289A
Authority: JP
Inventors: 金吾七戸; 正勝帆足; 根本　茂; 茂根本; ▲隆▼義山口
Original assignee: Nemoto Kyorindo Co Ltd
Current assignee: Nemoto Kyorindo Co Ltd
Priority date: 2015-07-14
Filing date: 2015-07-14
Publication date: 2020-03-25
Anticipated expiration: 2035-07-14
Also published as: JP2017018457A; WO2017010097A1

Description

本発明は、冠動脈ＣＴ画像と遅延造影ＣＴ画像とを用いた差分処理を行って差分画像を生成する画像処理方法、画像処理プログラム、画像処理装置及び画像処理システムに関する。

従来、ＣＴ画像から心筋梗塞部を抽出することが困難であることが知られていた。例えば、特許文献１の発明が解決しようとする課題には、責任血管を特定する技術では、心筋梗塞部が描出される領域を抽出できないと記載されている。この特許文献１では、ＣＴ画像から心筋梗塞部が描出される領域を抽出するために、ＭＲＩ（Magnetic Resonance Imaging）装置によって心臓が撮像されたＭＲ画像と、ＣＴ（Computed Tomography）装置によって心臓が撮像されたＣＴ画像とを用いている。

特開２０１０−１０４７１０号公報

心筋梗塞部では、心筋内膜側（内腔側）から線維化が進行する。造影剤は冠動脈を経て間質腔に浸透し、心筋梗塞部における造影剤の浸透速度は心筋正常部よりも遅い。そして、間質腔に浸透した造影剤の流出速度も、心筋正常部よりも心筋梗塞部において遅い。そのため、冠動脈に造影剤が到達してから所定の遅延時間が経過した遅延相で撮像した画像（遅延造影画像）においては、心筋正常部と心筋梗塞部とに画素値の差が生じる。ＭＲ画像においてはこの差が大きい（コントラストが高い）ため、ＣＴ画像と比較して心筋梗塞部の抽出がより容易である。

ＭＲＩ装置を用いる場合、ＡＴＰ薬剤（血管拡張剤）による負荷を利用して撮像が行われる。そして、ＭＲＩ装置は、ＣＴ装置の検査時間（１０〜１５分）よりも長い検査時間（４０〜６０分）を必要とする。そのため、被写体である患者の負担を少なくするために、ＣＴ装置で撮像したＣＴ画像から心筋梗塞部を識別することが望ましい。しかし、遅延相で撮像したＣＴ画像（遅延造影ＣＴ画像）では、心筋部と内腔との画素値の差が小さい（コントラストが低い）ため、医師が両者の境界を識別することが難しかった。さらに、遅延造影ＣＴ画像では、心筋梗塞部と心筋正常部との画素値の差も少ないため、医師が心筋梗塞部を識別することが難しかった。

そこで、ＣＴ画像に差分処理を行って、遅延造影ＣＴ画像における心筋梗塞部を強調することが考えられる。具体的には、造影された冠動脈を含む左心室を撮像して得られる冠動脈ＣＴ画像を、遅延造影ＣＴ画像から差し引く差分処理を行うことが考えられる。通常、差分処理は、撮像対象である左心室の平均画素値が高い（濃い）画像から、左心室の平均画素値が低い（淡い）画像を差し引くことによって行われる。そのため、一般に、淡い遅延造影ＣＴ画像から濃い冠動脈ＣＴ画像を差し引くという差分処理は行われない。

上記差分処理を行うことによって、左心室の内腔と心筋部との画素値の差を強調して、両者の境界を明確にすることができる。また、左心室の心筋正常部と心筋梗塞部の画素値の差を強調することもできる。しかし、この差分処理を行うことにより、新たな課題が生じてしまう。すなわち、通常は差分処理時に判別することを要しない遅延造影ＣＴ画像と冠動脈ＣＴ画像を判別する必要が生じてしまう。また、心臓の鼓動によって、非剛体である心臓、特に左心室が変形してしまう。そのため、左心室の心筋梗塞部を正確に識別するためには、遅延造影ＣＴ画像と冠動脈ＣＴ画像との位置合わせを行う必要がある。しかし、遅延造影ＣＴ画像では左心室の内腔と心筋部との境界が不明瞭であるため、冠動脈ＣＴ画像との位置合わせが困難であるという問題がある。

上記課題を解決するため、本発明の一例としての画像処理方法は、造影された冠動脈を含む左心室をＣＴ装置により撮像して得られた冠動脈ＣＴ画像と、前記冠動脈ＣＴ画像の撮像から所定の遅延時間が経過した後に、造影された心筋梗塞部を含む左心室をＣＴ装置により撮像して得られた遅延造影ＣＴ画像と、を用いた差分処理を行って差分画像を生成する画像処理方法であって、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを取得し、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを判別し、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との位置合わせを行い、前記遅延造影ＣＴ画像から前記冠動脈ＣＴ画像を差し引く前記差分処理を行うことを特徴とする。

また、本発明の他の例としての画像処理プログラムは、造影された冠動脈を含む左心室をＣＴ装置により撮像して得られた冠動脈ＣＴ画像と、前記冠動脈ＣＴ画像の撮像から所定の遅延時間が経過した後に、造影された心筋梗塞部を含む左心室をＣＴ装置により撮像して得られた遅延造影ＣＴ画像と、を用いた差分処理をコンピューターに実行させて差分画像を生成する画像処理プログラムであって、前記コンピューターに、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを判別する判別処理と、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との位置合わせを行う位置合わせ処理と、前記遅延造影ＣＴ画像から前記冠動脈ＣＴ画像を差し引く前記差分処理を実行させることを特徴とする。

また、本発明の他の例としての画像処理装置は、造影された冠動脈を含む左心室をＣＴ装置により撮像して得られた冠動脈ＣＴ画像と、前記冠動脈ＣＴ画像の撮像から所定の遅延時間が経過した後に、造影された心筋梗塞部を含む左心室をＣＴ装置により撮像して得られた遅延造影ＣＴ画像と、を用いた差分処理を行って差分画像を生成する画像処理装置であって、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを取得する取得部と、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを判別する判別部と、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との位置合わせを行う位置合わせ部と、前記遅延造影ＣＴ画像から前記冠動脈ＣＴ画像を差し引く前記差分処理を行う差分処理部を備えることを特徴とする。

また、本発明の他の例として画像処理システムは、上記画像処理装置と、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを撮像するＣＴ装置と、前記ＣＴ装置と接続され且つ造影剤を注入する注入装置とを備えることを特徴とする。

これにより、遅延造影ＣＴ画像と冠動脈ＣＴ画像とを判別することができる。また、遅延造影ＣＴ画像と冠動脈ＣＴ画像とを正確に位置合わせすると共に、心筋梗塞部を正確に識別することができる。

本発明のさらなる特徴は、添付図面を参照して例示的に示した以下の実施例の説明から明らかになる。

画像処理装置のブロック図である。冠動脈相のタイムデンシティカーブである。遅延相のタイムデンシティカーブである。非剛体変形位置合わせを説明する概略図である。差分処理を説明する概略図である。画像処理のフローチャートである。第２実施形態に係る画像処理システムの概略図である。第３実施形態に係る非剛体変形位置合わせを説明する概略図である。

以下、本発明を実施するための例示的な実施形態を、図面を参照して詳細に説明する。ただし、以下の実施形態で説明する寸法、材料、形状、構成要素の相対的な位置等は任意であり、本発明が適用される装置の構成又は様々な条件に応じて変更できる。また、特別な記載がない限り、本発明の範囲は、以下に具体的に記載された実施形態に限定されるものではない。

［第１実施形態］
図１に示す画像処理装置１は、ワークステーション又は検像システムであり、画像を処理する処理部１０としても機能するＣＰＵ（制御部）と、画像処理プログラムを記憶すると共に、画像データを記憶する記憶部１１とを有している。また、記憶部１１は、処理部１０が動作するためのシステムワークメモリであるＲＡＭ、画像処理プログラム、制御プログラム及びシステムソフトウェア等を格納するＲＯＭ、及びハードディスクドライブ等を有する。

ＣＰＵは、記憶部１１に記憶されたプログラムに基づいて画像処理装置１の全体を制御する。すなわち、ＣＰＵは、記憶された制御プログラムに従って、種々の演算、制御、判別等の処理動作を実行する。なお、ＣＰＵは、ＣＤ（Compact Disc）、又はインターネット上のサーバー等の外部記憶媒体に記憶されたプログラムに従って処理動作を実行することもできる。

画像処理装置１は、所定の指令あるいはデータを入力するキーボード、又は各種スイッチ等を含む入力部１２と、ＣＴ画像、装置の入力状態、設定状態、計測結果、及び各種情報等を表示する表示部１３とを備えている。なお、画像処理装置１は、入力部１２及び表示部１３の両者として機能するタッチパネルを備えていてもよい。

この画像処理装置１は、第１実施形態に係る画像処理プログラムに従って各種処理を実行するコンピューターである。そのため、処理部１０では、取得部１４、判別部１５、抽出部１６、算出部１７、位置合わせ部１８及び差分処理部１９等の各部が各種機能として論理的に実現されている。また、画像処理プログラムは、画像処理装置１を取得部１４、判別部１５、抽出部１６、算出部１７、位置合わせ部１８及び差分処理部１９として機能させる。この画像処理プログラムは、コンピューター読み取り可能な内部又は外部の記録媒体に記録することができる。

また、画像処理装置１は、外部記憶装置（サーバー）としてのＰＡＣＳ（Picture Archiving and Communication Systems）４と無線又は有線接続している。そして、画像処理装置１は、ＰＡＣＳ４へ情報を送信できると共に、ＰＡＣＳ４から情報を受信できる。なお、画像処理装置１は、ＲＩＳ（Radiology Information System）又はＨＩＳ（Hospital Information System）等の他の外部記憶装置と接続していてもよい。さらに、画像処理装置１はＣＴ装置３と無線又は有線接続しており、ＣＴ装置３は造影剤の注入装置２と無線又は有線接続している。そして、ＣＴ装置３及び注入装置２は、ＲＩＳ５と無線又は有線接続している。

注入装置２は、注入プロトコルを生成すると共に、該注入プロトコルに従って被写体である患者の体内に造影剤及び生理食塩水を注入する。ＣＴ装置３は、造影された冠動脈を含む左心室を撮像して冠動脈ＣＴ画像を得る。また、ＣＴ装置３は、冠動脈ＣＴ画像の撮像から所定の遅延時間が経過した後に、造影された心筋梗塞部を含む左心室を撮像して遅延造影ＣＴ画像を得る。さらに、ＣＴ装置３は、撮像したＣＴ画像を画像処理装置１又はＰＡＣＳ４に送信する。

ここでＣＴ画像の撮像タイミングを、図２及び図３を参照して説明する。図２は、造影剤によって冠動脈が造影されるとき（冠動脈相）を含む、概略タイムデンシティカーブ（ＴＤＣ）の一例を示す。図３は、冠動脈相から所定の遅延時間が経過したとき（遅延相）を含む、概略タイムデンシティカーブの一例を示す。図２及び図３では、縦軸が画素値を示し、横軸が造影剤の注入開始からの経過時間を示している。

図２では、冠動脈ＣＴ画像の撮像タイミングを縦棒で図示している。この撮像タイミングは、冠動脈が良好に造影されているタイミングであり、図２では注入開始から２８秒が経過した時点である。冠動脈は、注入開始から１５〜３０秒経過した時点で良好に造影される。図３では、遅延造影ＣＴ画像の撮像タイミングを縦棒で図示している。この撮像タイミングは、心筋正常部と心筋梗塞部の画素値の差が大きくなるタイミングであり、図３では冠動脈相から遅延時間３００秒（注入開始から３２８秒）が経過した時点である。この遅延時間は、５〜２０分が好ましい。

図１に戻り、画像処理装置１の取得部１４は、冠動脈ＣＴ画像と遅延造影ＣＴ画像とをＣＴ装置３又はＰＡＣＳ４から取得する。また、取得部１４は、画像処理装置１の記憶部１１に記憶された冠動脈ＣＴ画像と遅延造影ＣＴ画像とを取得することもできる。判別部１５は、取得部１４が取得した冠動脈ＣＴ画像と遅延造影ＣＴ画像とを判別し、必要に応じてＣＴ画像に判別情報（例えば、画像の種類を示すＩＤ）を付加する。

判別部１５は、ＣＴ画像の左心室を含む領域（例えば、中央領域）又はＣＴ画像全体における平均画素値を基準に、冠動脈ＣＴ画像と遅延造影ＣＴ画像とを判別することができる。具体的に判別部１５は、両ＣＴ画像を比較して、平均画素値が低い方が遅延造影ＣＴ画像であると判別する。そして、判別部１５は、遅延造影ＣＴ画像であることを示す判別情報を、遅延造影ＣＴ画像に付加する。また、判別部１５は、画素値とピクセル数（ボクセル数）とに基づくヒストグラムにおける、所定の画素値範囲のピクセル数（ボクセル数）を基準に判別してもよい。例えば、判別部１５は、両ＣＴ画像を比較して、２００ＨＵ〜５００ＨＵの範囲におけるピクセル数（ボクセル数）が少ない方が遅延造影ＣＴ画像であると判別する。

抽出部１６は、冠動脈ＣＴ画像における左心室の外輪郭（以下、単に「外輪郭」ともいう）と、遅延造影ＣＴ画像における左心室の外輪郭とを抽出する。外輪郭は心外膜の輪郭に対応し、後述する内輪郭は心内膜の輪郭に対応する。例えば、抽出部１６は、ＣＴ画像に対して任意の閾値による二値化を行い、画素値が大きい領域のみを抽出することによって心臓領域のみを抽出する。そして、抽出部１６は、抽出した心臓領域と他の領域との境界を外輪郭として抽出する。

位置合わせ部１８は、判別部１５が判別した冠動脈ＣＴ画像と遅延造影ＣＴ画像との位置合わせを行う。具体的に位置合わせ部１８は、冠動脈ＣＴ画像と遅延造影ＣＴ画像との剛体位置合わせを行う。例えば、位置合わせ部１８は、冠動脈ＣＴ画像と遅延造影ＣＴ画像とのそれぞれの左心室の重心を特定し、両重心が一致するように冠動脈ＣＴ画像の平行移動と回転移動とを行う。なお、位置合わせ部１８は、遅延造影ＣＴ画像の平行移動と回転移動とを行ってもよい。

また、位置合わせ部１８は、位置合わせを行う際に、冠動脈ＣＴ画像及び遅延造影ＣＴ画像の左心室の外輪郭同士の位置ズレが最小となるように、剛体位置合わせを行ってもよい。例えば、位置合わせ部１８は、両ＣＴ画像においてグリッド状に設定された複数の制御点のうち、外輪郭上の制御点を対応点として設定する。そして、位置合わせ部１８は、両外輪郭上の対応点同士の位置ズレが最小（又は所定の範囲内）となるように、冠動脈ＣＴ画像の平行移動と回転移動とを行う。外輪郭とその外側の領域との間では画素値の差が大きいため、淡い遅延造影ＣＴ画像及び濃い冠動脈ＣＴ画像のいずれにおいても外輪郭を正確に抽出することができる。そのため、正確に抽出可能な外輪郭を基準とすることで、正確に位置合わせを行うことができる。

さらに、位置合わせ部１８は、位置合わせを行う際に、冠動脈ＣＴ画像と遅延造影ＣＴ画像との非剛体変形位置合わせを行うこともできる。冠動脈ＣＴ画像の左心室の形状又は大きさは、心臓の鼓動又は体動等に起因して遅延造影ＣＴ画像と異なることがある。そのため、非剛体変形位置合わせを行うことにより、一方の画像の大きさ又は形状等を変形して、他方の画像の大きさ又は形状と一致させることによって、差分処理により適した位置合わせを行うことができる。

非剛体変形位置合わせを行う際、位置合わせ部１８は、遅延造影ＣＴ画像を変形せずに冠動脈ＣＴ画像を変形する。遅延造影ＣＴ画像においては、心筋正常部と心筋梗塞部との画素値の差が冠動脈ＣＴ画像よりも大きく、心筋梗塞部の領域がより正確に撮像されている。そのため、遅延造影ＣＴ画像を変形させないことによって、差分処理によって得られた差分画像においても、心筋梗塞部の領域を正確に描出することができる。非剛体変形位置合わせは、例えば、画素値の類似性（例えば、相互情報量又は相関係数）を基準とする尺度（類似度）を用いた方法によって行われる。そして、位置合わせ部１８は、類似度が最大又は所定の閾値以上となるように冠動脈ＣＴ画像を変形する。

非剛体変形位置合わせについて、左心室の概略断面を示す図４を参照して説明する。図４においては、冠動脈ＣＴ画像と遅延造影ＣＴ画像とは、それぞれの左心室の重心が一致するように位置合わせされた状態である。また、破線で示された遅延造影ＣＴ画像の左心室よりも、実線で示された冠動脈ＣＴ画像の左心室が小さい。さらに、冠動脈ＣＴ画像の内輪郭（内腔）の形状は、遅延造影ＣＴ画像と異なっている。

位置合わせ部１８は、両ＣＴ画像においてグリッド状に設定された複数の制御点のうち、冠動脈ＣＴ画像の外輪郭上の制御点を対応点Ｃとして設定する。また、位置合わせ部１８は、複数の制御点のうち、遅延造影ＣＴ画像の外輪郭上の制御点を対応点Ｄとして設定する。そして、位置合わせ部１８は、対応点Ｃと対応点Ｄとの位置が略一致するように、冠動脈ＣＴ画像の外輪郭を図４の点線で示す位置まで拡大する。その後、位置合わせ部１８は、両ＣＴ画像の類似度が最大又は所定の閾値以上となるように、冠動脈ＣＴ画像の外輪郭の形状を変形させる。さらに、位置合わせ部１８は、両ＣＴ画像の類似度が最大又は所定の閾値以上となるように、冠動脈ＣＴ画像の内輪郭の形状を変形させる。

一方、遅延造影ＣＴ画像の左心室よりも冠動脈ＣＴ画像の左心室が大きい場合、位置合わせ部１８は、対応点Ｃと対応点Ｄとの位置が略一致するように冠動脈ＣＴ画像の外輪郭を縮小する。その後、位置合わせ部１８は、両ＣＴ画像の類似度が最大又は所定の閾値以上となるように、冠動脈ＣＴ画像の外輪郭の形状を変形させる。さらに、位置合わせ部１８は、両ＣＴ画像の類似度が最大又は所定の閾値以上となるように、冠動脈ＣＴ画像の内輪郭の形状を変形させる。

ただし、遅延造影ＣＴ画像の内輪郭は、内腔と心筋部との画素値の差が少ないため抽出することが難しい。そのため、遅延造影ＣＴ画像の内輪郭に対応点を設定することは難しく、冠動脈ＣＴ画像の内輪郭を外輪郭と同様に変形することはできない。そこで、位置合わせ部１８は、心筋部の面積又は体積を利用して冠動脈ＣＴ画像の内輪郭の形状を以下のように変形させることができる。

すなわち、まず画像処理装置１の算出部１７が、冠動脈ＣＴ画像における左心室の心筋部（図４に斜線で示す領域）の面積を算出する。冠動脈ＣＴ画像においては、内腔と心筋部との画素値の差が大きいので、内輪郭と外輪郭とを抽出できる。そのため、算出部１７は、抽出された内輪郭と外輪郭との間の領域の面積を算出することによって心筋部の面積を算出できる。また、算出部１７は、取得部１４が取得したボリュームデータに含まれる多数のＣＴ画像に対して同じ処理を施すことによって、心筋部の体積を算出できる。なお、算出部１７は、心筋部に対する心筋梗塞部の割合、及び心筋梗塞部の心内膜（心筋壁部）からの広さ（深さ）を算出することもできる。

次に、位置合わせ部１８は、非剛体変形位置合わせを行う際に、算出した心筋部の面積又は体積を維持するように冠動脈ＣＴ画像（冠動脈ＣＴ画像の内輪郭）を変形する。この場合も、位置合わせ部１８は、両ＣＴ画像の類似度が最大又は所定の閾値以上となるように、冠動脈ＣＴ画像を変形させる。通常、冠動脈ＣＴ画像及び遅延造影ＣＴ画像は、心電同期を用いて同一心位相で撮像されている。そのため、両心筋部の面積及び体積は略同一であると仮定することができる。そこで、心筋部の面積又は体積を維持することにより、冠動脈ＣＴ画像の内輪郭をより実際の大きさ及び形状に近似させることができる。なお、面積又は体積は完全に一致させる必要はなく、所定の範囲内の誤差があっても面積又は体積が維持されているとみなすことができる。

左心室の概略断面を示す図５を参照して、遅延造影ＣＴ画像と冠動脈ＣＴ画像を用いた差分処理について説明する。図５において、左側には遅延造影ＣＴ画像が示され、中央には冠動脈ＣＴ画像が示され、右側には差分処理後の画像（差分画像）が示されている。なお、説明の便宜上、図５は位置合わせされた後のＣＴ画像を示している。

位置合わせ後に、画像処理装置１の差分処理部１９は、遅延造影ＣＴ画像から冠動脈ＣＴ画像を差し引く差分処理を行う。これにより、遅延造影ＣＴ画像から冠動脈ＣＴ画像の画素値を減算して、左心室内腔（内輪郭）を明瞭に描出できる。すなわち、心内膜を含む内腔の画素値が減算され、心筋部と内腔の画素値の差が強調される。また、心筋正常部と心筋梗塞部の画素値の差が強調され、心筋梗塞部を明瞭に描出できる。

遅延造影ＣＴ画像においては、例えば、心筋正常部の画素値が８０ＨＵであり、内腔の画素値が１００ＨＵであり、心筋梗塞部の画素値が１００ＨＵである。冠動脈ＣＴ画像においては、心筋正常部の画素値が５０ＨＵであり、内腔の画素値が１８０ＨＵであり、心筋梗塞部の画素値が２０ＨＵである。そのため、差分画像においては、心筋正常部の画素値が３０ＨＵ（＝８０ＨＵ−５０ＨＵ）であり、内腔の画素値が−８０ＨＵ（＝１００ＨＵ−１８０ＨＵ）であり、心筋梗塞部の画素値が８０ＨＵ（＝１００ＨＵ−２０ＨＵ）である。

差分処理後、画像処理装置１は生成した差分画像をＰＡＣＳ４に送信し、ＰＡＣＳ４は受信した差分画像を記憶する。そして、記憶された差分画像は、医師によって読影に用いられる。差分処理の結果、遅延造影ＣＴ画像においては差が無かったが、差分画像においては心筋梗塞部と内腔の画素値の差が強調されて１６０ＨＵになっている。また、遅延造影ＣＴ画像においては２０ＨＵであった心筋正常部と心筋梗塞部の画素値の差が強調されて、差分画像においては５０ＨＵになっている。この結果、差分画像においては、例えば比較的に画素値が低い内腔が黒く表示され、心筋部の心筋正常部が灰色に表示される。また、比較的に画素値が高い心筋梗塞部は白色に表示される。

遅延造影ＣＴ画像では、心筋梗塞部を描出して心筋正常部と区別させることが可能である。しかし、通常の遅延造影ＣＴ画像では、内腔と心筋梗塞部とが近い画素値を有する。そのため、心内膜から進行する心筋梗塞部を正確に識別することは難しい。一方、冠動脈ＣＴ画像では、内腔および冠動脈のみの画素値が高い。そのため、冠動脈ＣＴ画像では、内腔と心筋部とを容易に識別することが可能である。そこで、遅延造影ＣＴ画像から冠動脈ＣＴ画像を減算するという手法によれば、上記のように簡易且つ明瞭に心筋部を描出できる。これにより、医師は、正確且つ容易に心筋梗塞部を識別することができる。

次に、図６に示すフローチャートを参照して、画像処理方法について説明する。まず、画像処理装置１の取得部１４が、ＣＴ装置３又はＰＡＣＳ４（画像サーバー）から遅延造影ＣＴ画像と冠動脈ＣＴ画像を取得して、記憶部１１に一時的に記憶する（Ｓ１０１）。そして、判別部１５は、遅延造影ＣＴ画像と冠動脈ＣＴ画像とを判別する（Ｓ１０２）。このとき、オペレーターは表示部１３に表示された「（左心室）抽出ボタン」をクリックする。すると抽出部１６は、判別された冠動脈ＣＴ画像の外輪郭と内輪郭（及び外輪郭と内輪郭の間の心筋部）の形状（範囲）を抽出する（Ｓ１０３）。そして、表示部１３は抽出結果をＣＴ画像上に表示し、オペレーターは必要に応じて入力部１２を操作して外輪郭及び内輪郭の範囲を修正する。なお、抽出部１６は、オペレーターの指示によらずに自動で外輪郭と内輪郭の形状を抽出してもよい。

非剛体変形を行うために面積（体積）を算出する場合（Ｓ１０４でＹＥＳ）、算出部１７は、抽出された左心室の心筋部から、冠動脈ＣＴ画像の心筋部の面積（体積）を算出する（Ｓ１０５）。そして、抽出部１６は、判別部１５が判別した遅延造影ＣＴ画像の外輪郭の形状（範囲）を抽出する（Ｓ１０６）。この抽出結果を表示部１３がＣＴ画像上に表示するため、オペレーターは外輪郭及び内輪郭の範囲を修正できる。

また、面積を算出しない場合も（Ｓ１０４でＮＯ）、抽出部１６は遅延造影ＣＴ画像の外輪郭の範囲を抽出する（Ｓ１０６）。この抽出結果を表示部１３がＣＴ画像上に表示するため、オペレーターは外輪郭及び内輪郭の範囲を修正できる。そして、位置合わせ部１８は、冠動脈ＣＴ画像と遅延造影ＣＴ画像の位置合わせを行う（Ｓ１０７）。この位置合わせには、剛体位置合わせ及び非剛体変形位置合わせの少なくとも一方が含まれる。

剛体位置合わせの場合、位置合わせ部１８は、例えば、外輪郭で特定される左心室同士の重心が一致するように位置合わせを行う。また、位置合わせ部１８は、冠動脈ＣＴ画像と遅延造影ＣＴ画像のそれぞれの左心室の外輪郭を基準に、外輪郭同士のずれが最小又は所定範囲内となるように剛体位置合わせを行うこともできる。右心室は、冠動脈ＣＴ画像と遅延造影ＣＴ画像との間で変形の度合いが大きい。そのため、左心室の外輪郭を基準に剛体位置合わせをすることによって、より正確な位置合わせを行うことができる。

非剛体変形位置合わせの場合、位置合わせ部１８は、例えば、遅延造影ＣＴ画像の外輪郭と一致するように、冠動脈ＣＴ画像の外輪郭を変形して位置合わせを行う。このとき、位置合わせ部１８は、類似度が最大となるように冠動脈ＣＴ画像を変形する。心筋梗塞部がより明確に撮像された遅延造影ＣＴ画像を変形させないため、心筋梗塞部をより明確に描出することができる。

また、位置合わせ部１８は、ステップＳ１０５で算出した心筋部の面積（体積）が維持されるように、冠動脈ＣＴ画像を変形させることもできる。心電同期によって、いずれのＣＴ画像も拡張中期に撮像しているため、両ＣＴ画像における心筋部の断面積（体積）は大きく変わらないと仮定できる。また、左心室の外輪郭は、冠動脈ＣＴ画像と遅延造影ＣＴ画像とのいずれにおいても比較的明瞭に抽出できる。そのため、非剛体変形位置合わせ時に、冠動脈ＣＴ画像から算出した面積（体積）に基づいて冠動脈ＣＴ画像内を変形させることにより、より正確な位置合わせを行うことができる。

位置合わせ後、差分処理部１９は、遅延造影ＣＴ画像から冠動脈ＣＴ画像を差し引く差分処理を行い（Ｓ１０８）、作成された差分画像は記憶部１１に一時的に記憶される。これにより画像処理が終了し、表示部１３は差分画像を表示する。また、画像処理装置１は、差分画像をＰＡＣＳ４に送信して記憶させる。

なお、該フローにおいて、画像処理装置１の判別部１５は、ＣＴ装置３又はＰＡＣＳ４に記憶されている遅延造影ＣＴ画像と冠動脈ＣＴ画像とを判別してもよい。例えば、判別部１５は、ＣＴ画像の平均画素値を基準に冠動脈ＣＴ画像と遅延造影ＣＴ画像とを判別する。そして、取得部１４が遅延造影ＣＴ画像及び冠動脈ＣＴ画像を取得する際には、判別情報が両ＣＴ画像に付加されて記憶部１１に記憶される。差分処理時には、該判別情報に基づいて遅延造影ＣＴ画像と冠動脈ＣＴ画像とを判別することができる。また、遅延造影ＣＴ画像に対する抽出処理と、冠動脈ＣＴ画像に対する抽出処理及び算出処理とは順番を入れ替えることもできる。

この画像処理は、画像処理プログラムによって以下のように実行される。すなわち、冠動脈ＣＴ画像及び遅延造影ＣＴ画像の取得後に、オペレーターは、画像処理プログラムを起動して両ＣＴ画像を読み込む。画像処理プログラムは、画像処理装置１を判別部１５として機能させ、冠動脈ＣＴ画像と遅延造影ＣＴ画像とを判別する判別処理を実行させる（Ｓ１０２）。次に、画像処理プログラムは、画像処理装置１を抽出部１６として機能させ、冠動脈ＣＴ画像の心筋部を抽出させる抽出処理を実行させる（Ｓ１０３）。そして、画像処理プログラムは、画像処理装置１の表示部１３に抽出結果を表示させる。

面積を算出する場合（Ｓ１０４でＹＥＳ）、画像処理プログラムは、画像処理装置１を算出部１７として機能させ、冠動脈ＣＴ画像の心筋部の面積を算出させる算出処理を実行させる（Ｓ１０５）。その後、画像処理プログラムは、画像処理装置１を抽出部１６として機能させ、遅延造影ＣＴ画像の外輪郭を抽出する抽出処理を実行させる（Ｓ１０６）。なお、面積を算出しない場合も（Ｓ１０４でＮＯ）、同様に画像処理プログラムは、抽出処理を実行させる（Ｓ１０６）。そして、画像処理プログラムは、画像処理装置１を位置合わせ部１８として機能させ、冠動脈ＣＴ画像と遅延造影ＣＴ画像との位置合わせを行う位置合わせ処理を実行させる（Ｓ１０７）。

位置合わせ後、画像処理プログラムは、画像処理装置１を差分処理部１９として機能させ、遅延造影ＣＴ画像から冠動脈ＣＴ画像を差し引く差分処理を実行させる（Ｓ１０８）。そして、画像処理プログラムは、作成した差分画像を記憶部１１に一時的に記憶させる。これにより画像処理が終了し、画像処理プログラムは、表示部１３に差分画像を表示させる。その後、画像処理装置１は、差分画像をＰＡＣＳ４に送信して記憶させる。なお、画像処理プログラムは、心筋部に対する心筋梗塞部の割合及び心筋梗塞部の心内膜からの広さを算出するように、画像処理装置１に差分画像を用いた算出処理を実行させてもよい。

このような第１実施形態に係る発明によれば、遅延造影ＣＴ画像と冠動脈ＣＴ画像とを正確に判別して、差分処理を行うことができる。また、内腔と心筋部との画素値の差が小さい遅延造影ＣＴ画像であっても、冠動脈ＣＴ画像と正確に位置合わせをすることができる。また、第１実施形態に係る発明によれば、サブトラクションによって内腔と心筋部との画素値の差を強調して、境界を明確に描出することができる。また、心筋正常部と心筋梗塞部の画素値の差を強調して、心筋梗塞部を明確に描出することができる。これにより、ＭＲ画像よりも解像度が高いＣＴ画像において、医師がより正確に心筋梗塞部を識別できる。また、ＣＴ装置３による撮像はＭＲＩ装置よりも撮像時間が短く、心筋梗塞部を撮像する際に被写体である患者の負担を軽減できる。さらに、ＣＴ装置はＭＲＩ装置よりも１スキャンでの撮像幅が広いため、１回の撮像でより広範囲の画像を得ることができる。

なお、第１実施形態では冠動脈ＣＴ画像及び遅延造影ＣＴ画像が二次元画像である例を説明した。しかし、冠動脈ＣＴ画像及び遅延造影ＣＴ画像は、三次元画像（ボリュームデータ）であってもよい。また、上述した画像処理を複数の二次元画像に対して同様に行うことによって、三次元の差分画像を作成することもできる。例えば、３２０枚の二次元画像の全てにおいて心筋梗塞部を抽出する必要はない。冠動脈ＣＴ画像及び遅延造影ＣＴ画像のそれぞれ１枚の画像に対して、心筋梗塞部と心筋正常部との閾値を指定してやれば、その結果を全画像に反映させて三次元画像を作成することができる。

［第２実施形態］
図７を参照して、第２実施形態に係る画像処理システムについて説明する。なお、第２実施形態の説明においては、第１実施形態との相違点について説明し、第１実施形態で説明した構成要素については同じ参照番号を付し、その説明を省略する。特に説明した場合を除き、同じ参照符号を付した構成要素は略同一の動作及び機能を奏し、その作用効果も略同一である。

画像処理システム２００は、第１実施形態で説明した画像処理装置１と、冠動脈ＣＴ画像と遅延造影ＣＴ画像とを撮像するＣＴ装置３と、ＣＴ装置３と有線又は無線で接続され且つ造影剤を注入する注入装置２とを備えている。ＣＴ装置３は、撮像プロトコルに従って被写体の撮像を行う撮像部３１と、ＣＴ装置３全体を制御する制御装置３２と、ディスプレイ３３とを有している。撮像部３１は、寝台と、被写体にＸ線を照射するＸ線源と、被写体を透過したＸ線を検出するＸ線検出器等を有している。そして、撮像部３１は、被写体にＸ線を曝射し、被写体を透過したＸ線に基づいて被写体の体内を逆投影することで、被写体の透視画像を撮像する。なお、制御装置３２とディスプレイ３３とは、一体的に構成してもよい。

オペレーターは、医師が予め記憶させた検査オーダーをＲＩＳ５からＣＴ装置３へ読み出す。この検査オーダーは、撮像部位と、撮像方法（冠動脈相及び遅延相等の撮像タイミング）との少なくとも一つの情報を含んでいる。なお、検査オーダーは、被写体に関する情報として、被写体の体重、除脂肪体重、循環血液量、被写体番号（被写体ＩＤ）、被写体氏名、性別、生年月日、年齢、身長、血液量、血流速度、体表面積、被写体の疾病、副作用の履歴、クレアチニン値、心拍数、及び心拍出量等を含んでいてもよい。また、検査オーダーは、検査情報として、検査番号（検査ＩＤ）、検査部位、検査日時、薬液種類、及び薬液名称等を含んでいてもよい。

そして、オペレーターは、ＣＴ装置３に撮像部位及び撮像タイミング（冠動脈相及び遅延相）を入力する。そして、ＣＴ装置３は、冠動脈相及び遅延相の撮像プロトコルを作成する。この撮像プロトコルをオペレーターが確認して、修正が無ければ撮像プロトコルを確定させる。なお、オペレーターは、撮像部位又は撮像タイミングの一方のみを入力することもある。

撮像プロトコルには、例えば、撮像部位、実効管電圧、機種名、メーカー名、撮像時間、管電圧、撮像範囲、回転速度、ヘリカルピッチ、曝射時間、線量及び撮像方法等の情報が含まれている。そして、制御装置３２は、撮像プロトコルに従うように撮像部３１を制御して被写体を撮像する。また、制御装置３２は、ディスプレイ３３に接続されており、ディスプレイ３３には、装置の入力状態、設定状態、撮像結果、及び各種情報等が表示される。

ＣＴ装置３は、確定した撮像プロトコルを注入装置２に送信する。また、ＣＴ装置３は、基準注入プロトコルを注入装置２に送信することができる。例えば、ＣＴ装置３と注入装置２とは、ゲートウェイ装置（不図示）を介して接続されており、ＣＴ装置３は予め記憶していた基準注入プロトコルを注入装置２に送信する。

注入装置２は、シリンジに充填された造影剤を被写体の体内に注入するために、注入ヘッド２１を備えている。さらに、注入装置２は、注入ヘッド２１を保持するスタンド２２と、注入ヘッド２１に有線又は無線で接続されたコンソール２３とを有する。注入ヘッド２１は、造影剤が充填されたシリンジが搭載される第１保持部２１４と、造影剤を後押しするための生理食塩水が充填されたシリンジが搭載される第２保持部２１５とを有する。また、注入ヘッド２１は、第１保持部２１４に搭載されたシリンジ内の薬液を注入プロトコルに従って押し出す駆動機構（不図示）と、第２保持部２１５に搭載されたシリンジ内の薬液を注入プロトコルに従って押し出す駆動機構（不図示）とを有している。

コンソール２３は、注入ヘッド２１を制御する制御装置として機能すると共に、造影剤の注入プロトコルを生成する生成装置としても機能する。そして、コンソール２３は、ＣＴ装置３から受信した撮像プロトコルに基づいて注入プロトコルを生成する。なお、コンソール２３は、ＣＴ装置３から受信した基準注入プロトコルを修正することによって、注入プロトコルを生成することもできる。そして、オペレーターがコンソール２３に表示された注入プロトコルを確認すると、注入準備が完了する。

注入装置２は、コンソール２３に代えて、注入ヘッド２１に接続された制御装置と、該制御装置に接続され且つ薬液の注入状況等が表示されるタッチパネルディスプレイとを有していてもよい。また、注入ヘッド２１及び制御装置は、スタンド２２と一体的に構成することもできる。スタンド２２に代えて天吊部材を設け、該天吊部材を介して天井から注入ヘッド２１を天吊することもできる。

また、注入装置２は、電源又はバッテリー、コンソール２３に接続されるハンドスイッチ、及び注入ヘッド２１を遠隔操作する遠隔操作装置等を有していてもよい。この遠隔操作装置は、注入ヘッド２１を遠隔操作して注入を開始又は停止することができる。また、電源又はバッテリーは、注入ヘッド２１又はコンソール２３のいずれかに設けることができ、これらとは別に設けることもできる。

注入ヘッド２１は、制御装置が生成した注入プロトコルに従って造影剤を注入する。この注入プロトコルには、少なくとも注入速度及び注入時間が含まれている。また、注入プロトコルは、注入量、注入タイミング、造影剤濃度及び注入圧力等の注入条件に関する情報を含んでいてもよい。また、注入ヘッド２１は、注入条件、注入状況、装置の入力状態、設定状態、及び各種注入結果等が表示されるヘッドディスプレイ２１１と、駆動機構の動作を入力するための操作部２１２とを有している。ヘッドディスプレイ２１１は、タッチパネル等から構成することによって、操作部２１２として用いることもできる。なお、ヘッドディスプレイ２１１は省略することができる。

薬液が注入される際には、注入ヘッド２１に搭載されたシリンジの先端部にミキシングチューブ等が接続される。また、操作部２１２には、駆動機構の前進ボタン、駆動機構の後進ボタン、又は最終確認ボタン等が設けられている。そして、ミキシングチューブの接続等の注入準備が完了すると、オペレーターが最終確認ボタンを押す。これにより、注入ヘッド２１は、注入を開始できる状態で待機する。注入開始後、シリンジから押し出された薬液は、ミキシングチューブ等を介して、被写体の体内へ注入される。このミキシングチューブは、造影剤と希釈薬液との混合器としても機能する。このような混合器としては、株式会社根本杏林堂製の「ＳＰＩＲＡＬＦＬＯＷ（登録商標）」がある。

注入ヘッド２１には、例えば、ＲＦＩＤチップ、ＩＣタグ、バーコード等のデータキャリアを有するプレフィルドシリンジ等の、種々のシリンジを搭載することができる。そして、注入ヘッド２１は、シリンジに取り付けられたデータキャリアの読み取りを行う読取部を内蔵している。このデータキャリアには、薬液に関連する薬液情報が記憶されている。薬液情報としては、製品名称、製品ＩＤ、化学分類、含有成分、濃度、粘度、消費期限、シリンジ容量、シリンジ耐圧、シリンダ内径、ピストンストローク及びロット番号等がある。

造影剤の注入準備が完了し、オペレーターが注入装置２に注入開始を指示すると、注入装置２は造影剤のテストボーラス注入を行う。オペレーターは、予め関心領域（ＲＯＩ）を、例えば上行大動脈に設定する。そして、ＣＴ装置３は、心電同期を用いて、注入開始から例えば２８秒経過した時（冠動脈相）に冠動脈を撮像する。また、ＣＴ装置３は、撮像した冠動脈ＣＴ画像を画像処理装置１の取得部１４に送信する。さらに、ＣＴ装置３は、心電同期を用いて、冠動脈相から例えば遅延時間３００秒が経過した時に左心室を撮像する。そして、ＣＴ装置３は、撮像した遅延造影ＣＴ画像を取得部１４に送信する。なお、テストボーラス注入は省略することもできる。

ＣＴ装置３は、冠動脈ＣＴ画像及び遅延造影ＣＴ画像をＰＡＣＳ４に送信してもよい。さらに、ＣＴ装置３は、判別部１５が冠動脈ＣＴ画像と遅延造影ＣＴ画像とを判別できるように、ＣＴ画像に判別情報を付加してもよい。このような判別情報は、例えば、画像の種類を示すＩＤ、撮像の順番に応じた番号、及び撮像時間等である。

画像処理装置１においては、第１実施形態と同様に、判別部１５が冠動脈ＣＴ画像と遅延造影ＣＴ画像とを判別する。また、抽出部１６は、冠動脈ＣＴ画像における左心室の外輪郭と、遅延造影ＣＴ画像における左心室の外輪郭とを抽出する。そして、位置合わせ部１８が冠動脈ＣＴ画像と遅延造影ＣＴ画像との位置合わせを行い、差分処理部１９が遅延造影ＣＴ画像から冠動脈ＣＴ画像を差し引く。非剛体変形位置合わせを行う場合には、算出部１７が、冠動脈ＣＴ画像の左心室における心筋部の面積又は体積を算出する。

このような第２実施形態に係る発明によっても、遅延造影ＣＴ画像と冠動脈ＣＴ画像とを正確に判別して、差分処理を行うことができる。また、内腔と心筋部との画素値の差が小さい遅延造影ＣＴ画像であっても、冠動脈ＣＴ画像と正確に位置合わせをすることができる。また、サブトラクションによって内腔と心筋部との画素値の差を強調して、境界を明確に識別することができる。また、心筋正常部と心筋梗塞部の画素値の差を強調して、心筋梗塞部を明確に描出することができる。

［第３実施形態］
図４及び図８を参照して、第３実施形態に係る画像処理について説明する。この第３実施形態においては、非剛体変形位置合わせを行う際に、冠動脈ＣＴ画像における左心室の外輪郭の変位量に応じて、冠動脈ＣＴ画像における左心室の内輪郭を変形する。

なお、第３実施形態の説明においては、第１実施形態及び第２実施形態との相違点について説明し、第１実施形態及び第２実施形態で説明した構成要素については同じ参照番号を付し、その説明を省略する。特に説明した場合を除き、同じ参照符号を付した構成要素は略同一の動作及び機能を奏し、その作用効果も略同一である。

図４の冠動脈ＣＴ画像の左心室の外輪郭は、遅延造影ＣＴ画像よりも小さい。この場合、画像処理装置１の位置合わせ部１８は、冠動脈ＣＴ画像の外輪郭上の対応点Ｃが、遅延造影ＣＴ画像の外輪郭上の対応点Ｄと一致するように、冠動脈ＣＴ画像の外輪郭を拡大する。この拡大に伴う外輪郭の変位量及び変位方向は、図４において矢印で示すような変位ベクトルＡ，Ｂとして表すことができる。

変位量は、左心室の部位によって異なることがある。例えば、変位ベクトルＡ群は、変位ベクトルＢ群よりも変位量が大きい。そこで、位置合わせ部１８は、外輪郭の変位量に応じて、冠動脈ＣＴ画像における左心室の内輪郭を変形する。例えば、位置合わせ部１８は、対応点Ｃがどの位置に移動するかを示す変位ベクトルＡ，Ｂを算出する。そして、位置合わせ部１８は、第１実施形態で説明した類似度が最大又は所定の閾値以上となるように、算出した変位ベクトルＡ，Ｂに応じて非剛体変形位置合わせを行う。すなわち、変位量が大きい対応点Ｃの近傍に位置する部分の変位ベクトルＡ’（図８）が、変位量が小さい対応点の近傍に位置する部分の変位ベクトルＢ’よりも大きくなるように、内輪郭を変形させる。

その結果、冠動脈ＣＴ画像の内輪郭において、変位量が大きい対応点Ｃの近傍に位置する部分は、変位量が小さい対応点Ｃの近傍に位置する部分よりも変位量が大きくなる。一方、冠動脈ＣＴ画像の左心室の外輪郭が遅延造影ＣＴ画像よりも大きい場合、位置合わせ部１８は、冠動脈ＣＴ画像の外輪郭を縮小する。この場合も、位置合わせ部１８は、外輪郭の変位量に応じて、冠動脈ＣＴ画像の内輪郭を変形する。

このような第３実施形態に係る発明によれば、外輪郭の変位量に応じて内輪郭を変位させることができる。そのため、より正確に内輪郭を変形させて、非剛体変形位置合わせをすることができる。また、第３実施形態に係る発明によっても、遅延造影ＣＴ画像と冠動脈ＣＴ画像とを正確に判別して、差分処理を行うことができる。また、内腔と心筋部との画素値の差が小さい遅延造影ＣＴ画像であっても、冠動脈ＣＴ画像と正確に位置合わせをすることができる。また、サブトラクションによって内腔と心筋部との画素値の差を強調して、境界を明確に識別することができる。また、心筋正常部と心筋梗塞部の画素値の差を強調して、心筋梗塞部を明確に描出することができる。

［変形形態］
心筋正常部と心筋梗塞部との画素値の差が大きいＣＴ画像を用いて差分処理を行うと、差分画像における画素値の差が大きくなる。そのため、パーフュージョン画像を撮像する場合のように、１回の検査で心臓を複数回撮像する場合は、複数のＣＴ画像のうち画素値の差が大きいＣＴ画像を使用することが好ましい。ここで、画素値の差が最大となるのは、心筋正常部の画素値がピークとなる時点である。そのため、変形形態においては、画像処理装置１の取得部１４が、心筋正常部の画素値のピークに近い時点に撮像されたＣＴ画像を取得する。

例えば、心筋正常部の画素値は、冠動脈又は左心室内腔の画素値が１５０ＨＵに到達してから２秒経過後に上昇を開始する。そして、到達から７〜８秒経過時に、心筋正常部及び心筋梗塞部の画素値の差がピークに達する。その後、到達から１６〜１７秒経過時に、画素値の差が判別困難なレベルに達する。

そのため、判別部１５は、以下の条件で撮像されたＣＴ画像を判別して差分処理に使用することが好ましい。すなわち、上行大動脈にＲＯＩを設定した状態で、ＲＯＩの画素値が１５０ＨＵに到達した後に撮像されたＣＴ画像、ＲＯＩの画素値が１５０ＨＵに到達してから所定時間（例えば７秒〜１７秒）内に撮像されたＣＴ画像、ＲＯＩの画素値が上昇して下降を開始した（ピーク）後に撮像されたＣＴ画像、及びＲＯＩの画素値が上昇して下降を開始してから２秒〜９秒の範囲内に撮像されたＣＴ画像が使用されることが好ましい。なお、ＲＯＩは、冠動脈又は左心室内腔に設定することもできる。

これにより、画像処理装置１は、正確に心筋梗塞部を描出することができる。なお、遅延相についても、複数のＣＴ画像のうち心筋正常部と心筋梗塞部との画素値の差が最大となるＣＴ画像を用いることが好ましい。また、判別部１５は、ＣＴ装置３がＣＴ画像に付加した撮像時間情報に基づいて、複数のＣＴ画像から画素値の差が大きいＣＴ画像を判別することができる。さらに、判別部１５は、各ＣＴ画像の心筋正常部と心筋梗塞部との画素値の差を算出して、画素値の差が最大となるＣＴ画像を判別してもよい。

以上、各実施形態を参照して本発明について説明したが、本発明は上記実施形態に限定されるものではない。本発明に反しない範囲で変更された発明、及び本発明と均等な発明も本発明に含まれる。また、上述の各実施形態及び各変形形態は、本発明に反しない範囲で適宜組み合わせることができる。

なお、剛体位置合わせ、非剛体変形位置合わせ又は外輪郭の抽出後に、オペレーターは、画像の位置若しくは外輪郭の形状等を手動で修正することができる。また、オペレーターは、剛体位置合わせ、非剛体位変形置合わせ、外輪郭の特定又は差分処理後に、公知の方法を用いて、ノイズ除去等の画像補正を行うことができる。

さらに、剛体位置合わせ、非剛体変形位置合わせ及び外輪郭の抽出は、上記の方法には限定されない。例えば、非剛体変形位置合わせでは、冠動脈ＣＴ画像と遅延造影ＣＴ画像とのそれぞれの外輪郭に沿って対応点を配置し、外輪郭同士の類似度を向上させるように、冠動脈ＣＴ画像を変形させてもよい。この場合、対応点以外の冠動脈ＣＴ画像内の部分の移動量は、対応点の移動量から補完することができる。また、非剛体変形位置合わせを行う際には、遅延造影ＣＴ画像を変形せずに冠動脈ＣＴ画像を変形することが望ましいが、遅延造影ＣＴ画像の心筋梗塞部が変形されなければよい。そのため、遅延造影ＣＴ画像の少なくとも左心室の外輪郭の内側領域が変形されなければよく、外輪郭の外側領域は変形されてもよい。つまり、非剛体変形位置合わせを行う際には、遅延造影ＣＴ画像における少なくとも左心室の外輪郭の内側領域を変形せずに、冠動脈ＣＴ画像を変形することもできる。

１：画像処理装置、２：注入装置、３：ＣＴ装置、１４：取得部、１５：判別部、１６：抽出部、１７：算出部、１８：位置合わせ部、１９：差分処理部、２００：画像処理システム

Claims

造影された冠動脈を含む左心室をＣＴ装置により撮像して得られた冠動脈ＣＴ画像と、前記冠動脈ＣＴ画像の撮像から所定の遅延時間が経過した後に、造影された心筋梗塞部を含む左心室をＣＴ装置により撮像して得られた遅延造影ＣＴ画像と、を用いた差分処理を行って差分画像を生成する画像処理方法であって、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを取得し、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを判別し、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との位置合わせを行い、
前記遅延造影ＣＴ画像から前記冠動脈ＣＴ画像を差し引く前記差分処理を行う、画像処理方法。
前記判別を行う際に、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との間で左心室を含む領域の平均画素値を比較し、又は前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との間で画像全体の平均画素値を比較し、前記平均画素値が低い画像が前記遅延造影ＣＴ画像であると判別する、請求項１に記載の画像処理方法。
前記位置合わせを行う際に、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との一方の大きさ又は形状を変形して位置合わせを行う、請求項１又は２に記載の画像処理方法。
前記位置合わせを行う際に、前記遅延造影ＣＴ画像を変形せずに、前記冠動脈ＣＴ画像を変形する、請求項３に記載の画像処理方法。
前記冠動脈ＣＴ画像における左心室の心筋部の面積又は体積を算出し、
前記位置合わせを行う際に、前記心筋部の面積又は体積を維持するように前記冠動脈ＣＴ画像を変形する、請求項４に記載の画像処理方法。
前記位置合わせを行う際に、前記冠動脈ＣＴ画像における左心室の外輪郭の変位量に応じて、前記冠動脈ＣＴ画像における左心室の内輪郭を変形する、請求項４又は５に記載の画像処理方法。
前記冠動脈ＣＴ画像における左心室の外輪郭と、前記遅延造影ＣＴ画像における左心室の外輪郭とを抽出し、
前記位置合わせを行う際に、前記冠動脈ＣＴ画像及び前記遅延造影ＣＴ画像の前記外輪郭同士の位置ズレが最小となるように、剛体位置合わせを行う、請求項１または２に記載の画像処理方法。
造影された冠動脈を含む左心室をＣＴ装置により撮像して得られた冠動脈ＣＴ画像と、前記冠動脈ＣＴ画像の撮像から所定の遅延時間が経過した後に、造影された心筋梗塞部を含む左心室をＣＴ装置により撮像して得られた遅延造影ＣＴ画像と、を用いた差分処理をコンピューターに実行させて差分画像を生成する画像処理プログラムであって、前記コンピューターに、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを判別する判別処理と、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との位置合わせを行う位置合わせ処理と、
前記遅延造影ＣＴ画像から前記冠動脈ＣＴ画像を差し引く前記差分処理を実行させる、画像処理プログラム。
造影された冠動脈を含む左心室をＣＴ装置により撮像して得られた冠動脈ＣＴ画像と、前記冠動脈ＣＴ画像の撮像から所定の遅延時間が経過した後に、造影された心筋梗塞部を含む左心室をＣＴ装置により撮像して得られた遅延造影ＣＴ画像と、を用いた差分処理を行って差分画像を生成する画像処理装置であって、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを取得する取得部と、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを判別する判別部と、
前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像との位置合わせを行う位置合わせ部と、
前記遅延造影ＣＴ画像から前記冠動脈ＣＴ画像を差し引く前記差分処理を行う差分処理部を備える、画像処理装置。
請求項９に記載の画像処理装置と、前記冠動脈ＣＴ画像と前記遅延造影ＣＴ画像とを撮像するＣＴ装置と、前記ＣＴ装置と接続され且つ造影剤を注入する注入装置とを備える、画像処理システム。