特許 7236747 コンピュータプログラム、画像処理装置、および画像処理方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-03-02

(45)【発行日】2023-03-10

(54)【発明の名称】コンピュータプログラム、画像処理装置、および画像処理方法

(51)【国際特許分類】

   A61B 6/03 20060101AFI20230303BHJP

【ＦＩ】

A61B6/03 370A

A61B6/03 375

A61B6/03 360Z

【請求項の数】 8

(21)【出願番号】P 2020522222

(86)(22)【出願日】2019-05-28

(86)【国際出願番号】 JP2019021147

(87)【国際公開番号】W WO2019230738

(87)【国際公開日】2019-12-05

【審査請求日】2022-04-22

(31)【優先権主張番号】P 2018102516

(32)【優先日】2018-05-29

(33)【優先権主張国・地域又は機関】JP

(73)【特許権者】

【識別番号】504147254

【氏名又は名称】国立大学法人愛媛大学

(74)【代理人】

【識別番号】100181250

【弁理士】

【氏名又は名称】田中 信介

(72)【発明者】

【氏名】城戸 輝仁

【審査官】櫃本 研太郎

(56)【参考文献】

【文献】国際公開第２０１３／０５７９８２（ＷＯ，Ａ１）

【文献】国際公開第２０１６／００９９５７（ＷＯ，Ａ１）

【文献】米国特許出願公開第２０１４／０１６３４０３（ＵＳ，Ａ１）

【文献】国際公開第２０１４／１６８２０６（ＷＯ，Ａ１）

【文献】国際公開第２０１０／１０１１８４（ＷＯ，Ａ１）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ６／００－６／１４

(57)【特許請求の範囲】

【請求項1】

造影剤が投与された後の被験者の臓器を撮像したものであって、時系列に沿った所定数のフレームのＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像からなる画像データを記憶する記憶手段を備えた画像処理装置のためのコンピュータプログラムにおいて、
前記記憶手段は、
第１所定数のフレームのＣＴ画像からなる第１画像データと、
前記第１所定数に比べて少ない第２所定数のフレームのＣＴ画像からなる第２画像データと、
を記憶しており、
前記画像処理装置に、
前記第１画像データに基づいて、ＣＴ値の経時変化を特定する第１ステップと、
前記第１画像データに基づき前記第１ステップで特定された前記ＣＴ値の経時変化において、前記ＣＴ値が最大値となるタイミングである最大タイミングを特定する第２ステップと、
前記第２ステップで特定された前記最大タイミングに基づく前記第２所定数のフレームのＣＴ画像からなる前記第２画像データに基づいて、ＣＴ値の経時変化を特定する第３ステップと、
前記第２画像データに基づき前記第３ステップで特定された前記ＣＴ値の経時変化に基づいて、所定時間に対する前記ＣＴ値の傾きであって所定の傾きを特定する第４ステップと、
を実行させることを特徴とするコンピュータプログラム。

【請求項2】

前記所定の傾きは、前記所定時間に対する前記ＣＴ値の傾きであって所定値より大きい傾きである
ことを特徴とする請求項１に記載のコンピュータプログラム。

【請求項3】

前記画像処理装置に、
前記第２画像データに基づき前記第４ステップで特定された前記所定の傾きに基づいて、前記ＣＴ値の経時変化を所定の関数で近似する第５ステップをさらに実行させる
ことを特徴とする請求項１に記載のコンピュータプログラム。

【請求項4】

前記第２画像データは、前記第２ステップで特定された前記最大タイミングに基づいた前記臓器の１回の拍動に相当する前記第２所定数のフレームのＣＴ画像からなる
ことを特徴とする請求項１に記載のコンピュータプログラム。

【請求項5】

造影剤が投与された後の被験者の臓器を撮像したものであって、時系列に沿った第１所定数のフレームのＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像からなる第１画像データと前記第１所定数に比べて少ない第２所定数のフレームのＣＴ画像からなる第２画像データとを記憶する記憶手段と、
前記第１画像データに基づいて、ＣＴ値の経時変化を特定する第１手段と、
前記第１画像データに基づき前記第１手段で特定された前記ＣＴ値の経時変化において、前記ＣＴ値が最大値となるタイミングである最大タイミングを特定する第２手段と、
前記第２手段で特定された前記最大タイミングに基づく前記第２所定数のフレームのＣＴ画像からなる前記第２画像データに基づいて、ＣＴ値の経時変化を特定する第３手段と、
前記第２画像データに基づき前記第３手段で特定された前記ＣＴ値の経時変化に基づいて、所定時間に対する前記ＣＴ値の傾きであって所定の傾きを特定する第４手段と、
を備えることを特徴とする画像処理装置。

【請求項6】

前記第２画像データに基づき前記第４手段で特定された前記所定の傾きに基づいて、前記ＣＴ値の経時変化を所定の関数で近似する第５手段をさらに備える
ことを特徴とする請求項５に記載の画像処理装置。

【請求項7】

造影剤が投与された後の被験者の臓器を撮像したものであって、時系列に沿った所定数のフレームのＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像からなる画像データを記憶する記憶手段を備えた画像処理装置が実行する画像処理方法であって、
前記記憶手段は、
第１所定数のフレームのＣＴ画像からなる第１画像データと、
前記第１所定数に比べて少ない第２所定数のフレームのＣＴ画像からなる第２画像データと、
を記憶しており、
前記画像処理装置は、
前記第１画像データに基づいて、ＣＴ値の経時変化を特定する第１ステップと、
前記第１画像データに基づき前記第１ステップで特定された前記ＣＴ値の経時変化において、前記ＣＴ値が最大値となるタイミングである最大タイミングを特定する第２ステップと、
前記第２ステップで特定された前記最大タイミングに基づく前記第２所定数のフレームのＣＴ画像からなる前記第２画像データに基づいて、ＣＴ値の経時変化を特定する第３ステップと、
前記第２画像データに基づき前記第３ステップで特定された前記ＣＴ値の経時変化に基づいて、所定時間に対する前記ＣＴ値の傾きであって所定の傾きを特定する第４ステップと、
を実行することを特徴とする画像処理方法。

【請求項8】

前記画像処理装置は、
前記第２画像データに基づき前記第４ステップで特定された前記所定の傾きに基づいて、前記ＣＴ値の経時変化を所定の関数で近似する第５ステップをさらに実行する
ことを特徴とする請求項７に記載の画像処理方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、画像処理技術に関し、特に時系列のダイナミック画像を解析する技術に関する。

【背景技術】

【0002】

被験者に造影剤を投与した後、その被験者の心臓を撮像した画像を使って、心臓の血流を解析する手法が知られている。このような心筋血流の解析手法の一つとして、冠動脈ＣＴ検査時に追加して行われるＣＴＰ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ Ｐｅｒｆｕｓｉｏｎ）がある。このＣＴＰは、例えば、造影剤の心筋へのファーストパスを造影効果として観察することで、心筋血流の状態を評価する方法である。また、このＣＴＰにおける評価法としては、徐々に染まっていく心筋を複数心拍で撮像するＤｙｎａｍｉｃ撮影によって得られたＴＤＣ（Ｔｉｍｅ ｄｅｎｓｉｔｙ ｃｕｒｖｅ）を解析することで、心筋血流を定量評価するものがある。

【0003】

このようにＴＤＣの解析によって心筋血流等を定量評価する手法としては、例えば、Ｄｙｎａｍｉｃ撮影によって得られたＣＴ画像を解析して、入力関数および出力関数を導出して、導出した入力関数および出力関数に基づいて所定の心筋領域での造影剤が到達した時刻であるアライバルタイム（以下、「ＡＴ」とする。）や臓器内画素位置の画素のベースとなるＣＴ値であるベース値を算出したうえで、心筋血流を定量評価する方法がある（例えば、特許文献１を参照。）。

【先行技術文献】

【特許文献】

【0004】

【文献】国際公開２０１６／００９９５７号

【発明の概要】

【発明が解決しようとする課題】

【0005】

しかしながら、上述した特許文献１に記載された方法では、少なくとも３０～４０心拍数におけるＤｙｎａｍｉｃ撮影によって得られた複数フレームのＣＴ画像を解析する必要があり、解析処理が煩雑であるため、より簡便な手法で時系列画像を解析できるような解析方法が望まれている。また、少なくとも３０～４０心拍数に応じた時間のＣＴＰ検査が必要であるため、被験者に対しての被ばく量が多くなってしまう、という問題があった。そのため、撮像タイミングを的確に把握し、撮像時間を短くして、被験者に対する被ばく量を少なくする、ということも望まれている。

【0006】

そこで、本発明の目的は、従来に比べてより簡便な手法で時系列画像を解析できる解析手法を提供することである。また、本発明の別の目的は、従来に比べて被験者に対する被ばく量を少なくするとともに、時系列画像の解析において客観性及び定量性を担保した解析手法を提供することである

【課題を解決するための手段】

【0007】

本発明は、上述の課題の少なくとも一部を解決するためになされたものであり、以下の適用例として実現することが可能である。なお、本欄における括弧内の参照符号や補足説明等は、本発明の理解を助けるために、後述する実施例との対応関係を示したものであって、本発明を何ら限定するものではない。

【0008】

本発明における適用例１のコンピュータプログラムは、造影剤が投与された後の被験者の臓器を撮像したものであって、時系列に沿った所定数のフレームのＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像からなる画像データを記憶する記憶手段（１１０）を備えた画像処理装置（１００）のためのコンピュータプログラムであって、前記画像処理装置に、前記所定数のフレームのＣＴ画像からなる画像データに基づいて、ＣＴ値の経時変化を特定する第１ステップ（Ｓ３５０、Ｓ５４０）と、前記第１ステップで特定された前記ＣＴ値の経時変化に基づいて、所定時間に対する前記ＣＴ値の傾きであって所定の傾きを特定する第２ステップ（Ｓ１１００）と、前記第２ステップで特定された前記所定の傾きに基づいて、前記ＣＴ値の経時変化を所定の関数で近似する第３ステップ（Ｓ１３００）と、を実行させることを要旨とする。

【0009】

本発明における適用例２のコンピュータプログラムは、適用例１のコンピュータプログラムであって、前記所定の傾きは、前記所定時間に対する前記ＣＴ値の傾きであって所定値より大きい傾きであることを要旨とする。

【0010】

本発明における適用例３のコンピュータプログラムは、適用例１または適用例２のコンピュータプログラムであって、 前記記憶手段は、第１所定数のフレームのＣＴ画像からなる第１画像データと、前記第１所定数に比べて少ない第２所定数のフレームのＣＴ画像からなる第２画像データと、を記憶しており、前記画像処理装置に、前記第１画像データに基づき前記第１ステップで特定された前記ＣＴ値の経時変化において、前記ＣＴ値が最大値となるタイミングである最大タイミングを特定する第４ステップ（Ｓ４８０）、をさらに実行させ、前記第４ステップで特定された前記最大タイミングに基づく前記第２所定数のフレームのＣＴ画像からなる前記第２画像データに基づいて、前記第１ステップ、前記第２ステップ、および前記第３ステップを実行させる、ことを要旨とする。

【0011】

本発明における適用例４のコンピュータプログラムは、適用例３のコンピュータプログラムであって、前記画像処理装置に、前記第４ステップで特定された前記最大タイミングに基づいた前記臓器の１回の拍動に相当する前記第２所定数のフレームのＣＴ画像からなる前記第２画像データに基づいて、前記第１ステップ、前記第２ステップ、および前記第３ステップを実行させる、ことを要旨とする。

【0012】

本発明における画像処理装置は、造影剤が投与された後の被験者の臓器を撮像したものであって、時系列に沿った所定数のフレームのＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像からなる画像データを記憶する記憶手段（１１０）と、前記所定数のフレームのＣＴ画像からなる画像データに基づいて、ＣＴ値の経時変化を特定する第１手段（１２３、１３３）と、前記第１手段によって特定された前記ＣＴ値の経時変化に基づいて、所定時間に対する前記ＣＴ値の傾きであって所定の傾きを特定する第２手段（１２５、１３７）と、前記第２手段によって特定された前記所定の傾きに基づいて、前記ＣＴ値の経時変化を所定の関数で近似する第３手段（１２５、１３７）と、を備えることを要旨とする。

【0013】

本発明における画像処理方法は、造影剤が投与された後の被験者の臓器を撮像したものであって、時系列に沿った所定数のフレームのＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像からなる画像データを記憶する記憶手段（１１０）を備えた画像処理装置（１００）が実行する画像処理方法であって、前記画像処理装置が、前記所定数のフレームのＣＴ画像からなる画像データに基づいて、ＣＴ値の経時変化を特定する第１ステップ（Ｓ３５０、Ｓ５４０）と、前記第１ステップで特定された前記ＣＴ値の経時変化に基づいて、所定時間に対する前記ＣＴ値の傾きであって所定の傾きを特定する第２ステップ（Ｓ１１００）と、前記第２ステップで特定された前記所定の傾きに基づいて、前記ＣＴ値の経時変化を所定の関数で近似する第３ステップ（Ｓ１３００）と、を実行することを要旨とする。

【図面の簡単な説明】

【0014】

【図1】本発明の一つの実施例である実施例１における画像処理装置１００の全体構成図である。

【図2】ＣＴ画像の一例を例示した説明図である。

【図3】入力関数導出処理を示すフローチャートである。

【図4】入力関数特定処理を示すフローチャートである。

【図5】出力関数導出処理を示すフローチャートである。

【図6】出力関数特定処理を示すフローチャートである。

【図7】（Ａ）は入力関数データ記憶部１２９の入力関数テーブル７００のデータ構成の一例を示した説明図であり、（Ｂ）は出力関数データ記憶部１３９の出力関数テーブル７５０のデータ構成の一例を示した説明図である。

【図8】（Ａ）は実施例１における大動脈の第１ＴＤＣの一例を示した説明図であり、（Ｂ）は実施例１における心筋内冠動脈の第２ＴＤＣの一例を示した説明図である。

【図9】（Ａ）は実施例２における大動脈の第１ＴＤＣの一例を示した説明図であり、（Ｂ）は実施例２における心筋内冠動脈の第２ＴＤＣの一例を示した説明図である。

【発明を実施するための形態】

【0015】

以下、本発明が適用された実施例について図面を用いて説明する。なお、本発明の実施の形態は、下記の実施例に何ら限定されることはなく、本発明の技術的範囲に属する限り種々の形態を採りうる。

【実施例1】

【0016】

＜画像処理装置１００の構成＞
図１を参照して、まずは、本実施例の画像処理装置１００の構成について説明する。図１は、画像処理装置１００の全体構成図である。画像処理装置１００は、例えば、汎用的なコンピュータシステムにより構成され、以下に説明する画像処理装置１００内の個々の構成要素又は機能を、コンピュータが読み取り可能な記録媒体等に保存されたコンピュータプログラムを実行することにより実現される。本実施例の画像処理装置１００は、画像データ記憶部１１０と、入力関数導出部１２０と、入力関数データ記憶部１２９と、出力関数導出部１３０と、出力関数データ記憶部１３９と、血流解析処理部１５０とを備える。また、画像処理装置１００には、入力装置１６０と、表示装置１７０と、が接続されている。

【0017】

画像処理装置１００は、被験者に造影剤を投与したあとに撮像された被験者の臓器のＣＴ画像を用いて、その臓器の血流量の定量解析を行う装置である。本実施形態では、画像処理装置１００は、被験者の心臓を心電図に同期して撮像した同位相の複数枚のＣＴのフレーム画像を対象とし、これら心臓のＣＴ画像を画素単位で解析することにより、心臓および心臓近傍の血流量の定量値を算出するために用いる入力関数及び出力関数等を特定する。本実施形態において、画像処理装置１００は、静脈に投与された造影剤によるＣＴ画像の画素値の変化に基づいて心臓等の血流量を推定し、例えば、対象となる部位に造影剤が到達することによってＣＴ画像の画素値が急激に変化した場合、画素値が変化し始める前の画素値をベース値としたり、画素値が変化し始める時刻をＡＴとしたりする。なお、本造影剤としては、種々のものを利用することができるが、本実施形態においては、例えば、イオプロミドを利用するものとする。

【0018】

画像データ記憶部１１０は、コンピュータシステムにおける一般的なメモリであって、本実施形態では、心臓の撮像によって得られた心臓及び心臓近傍のＣＴ画像の画像データを記憶する。ＣＴ画像は、造影剤が投与された後の被験者の臓器を撮像した時系列の複数フレームのＣＴ画像であって、複数のスライス画像（短軸断層画像）で構成される３次元画像（３次元ボクセルデータ）であり、詳細については後述する。

【0019】

入力関数導出部１２０は、画像データ記憶部１１０に保存されている画像データに基づいて、被験者の臓器に流入する造影剤がＣＴ画像の画素値（以下、「ＣＴ値」と称する。）に及ぼす変化に関係する入力関数等を導出する。なお、入力関数導出部１２０は、以下に説明するＲＯＩ設定部１２１と、ＲＯＩ経時変化特定部１２３と、入力関数化処理部１２５と、を有している。入力関数データ記憶部１２９は、入力関数導出部１２０で導出された入力関数等に係るデータを後述する入力関数テーブル７００として記憶する。

【0020】

ＲＯＩ設定部１２１は、医師等の解析者による入力装置１６０を介した操作に従って、ＣＴ画像において入力関数を特定するための領域であるＲＯＩ（Ｒｅｇｉｏｎ Ｏｆ Ｉｎｔｅｒｅｓｔ；関心領域）を設定する。具体的に、ＲＯＩ設定部１２１は、画像データ記憶部１１０から画像データを読み出して、対象とする領域が写っているスライスのフレーム画像を表示装置１７０に表示させ、解析者が入力装置３を用いて行った入力に従って、表示装置１７０に表示されている画像上にＲＯＩを設定する。なお、このＲＯＩの位置は、同一スライスの全フレームに対して共通であるものとする。

【0021】

ＲＯＩ経時変化特定部１２３は、時系列の複数フレームのＣＴ画像の画像データに基づいて、ＲＯＩ内のＣＴ値の経時変化を特定する。具体的に、ＲＯＩ経時変化特定部１２３は、ＣＴ画像における所定の領域のＣＴ値の経時変化として、ＲＯＩ内のＣＴ値から定まるＲＯＩ値のＴＤＣを生成する。なお、ＲＯＩ値とは、ＲＯＩ内のＣＴ値を統計的アルゴリズムで処理した値（統計値）であり、例えば、平均値、最頻値、中央値、最大値、最小値などのいずれかを採り得るものとする。

【0022】

入力関数化処理部１２５は、ＲＯＩ経時変化特定部１２３により生成されたＴＤＣに基づいて、ＲＯＩ内のＣＴ値が最大値Ｍａｘとなる時刻であるＭＴ、ＲＯＩが設定された領域に造影剤が到達した時刻であるＡＴ、後述する第２ＴＤＣを得るための撮像する期間の撮像区間であるＭＢ、入力関数等を特定する。なお、入力関数化処理部１２５は、ＲＯＩ経時変化特定部１２３によって生成されたＴＤＣにおいて、ＴＤＣがなめらかな曲線になるように平滑化処理を実行したうえで、入力関数等を特定する。また、入力関数化処理部１２５は、ＣＴ値が急激に上昇した後の所定のフレームである上限フレームＦａや、後述するベース値であるＢＬも特定する。

【0023】

出力関数導出部１３０は、画像データ記憶部１１０に保存されている画像データに基づいて、被験者の臓器の血管に造影剤が流入したことによりＣＴ画像のＣＴ値に及ぼす変化に関係する出力関数等を導出する。なお、出力関数導出部１３０は、ＲＯＩ単位ではなく画素単位に出力関数を導出するものであり、対象画素抽出部１３１と、ＣＴ経時変化特定部１３３と、出力関数化処理部１３７と、を有している。出力関数データ記憶部１３９は、出力関数導出部１３０で導出された出力関数に係るデータを後述する出力関数テーブル７５０として記憶する。

【0024】

対象画素抽出部１３１は、出力関数の導出対象となる画素を抽出し、フレーム画像上でその位置を特定する。具体的に、対象画素抽出部１３１は、抽出処理として、時系列の複数フレームのＣＴ画像のすべてにおいて、スライス単位にＣＴ値が所定の範囲内（例えば、５０～１６０）である画素を選択する。なお、この抽出処理において、目的とする対象画素の領域内で小さな対象外となる領域が含まれる場合には、その対象外の領域の大きさに応じて（例えば、画素数が所定数以下）、この領域を処理上の欠損とみなし、対象画素に転換する処理を実行することができる。また、対象画素の領域外で孤立した対象となる画素がある場合は、孤立した対象領域の大きさに応じて（例えば、画素数が所定数以下）、この領域を対象外とする処理を実行することもできる。そして、対象画素抽出部１３１は、このようにスライス毎に抽出された画素の位置を、出力関数を算出する対象である対象画素位置（臓器内画素位置）として特定する。

【0025】

ＣＴ経時変化特定部１３３は、時系列の複数フレームのＣＴ画像の画像データに基づいて、対象画素位置（臓器内画素位置）の画素のＣＴ値の経時変化を特定する。具体的に、ＣＴ経時変化特定部１３３は、スライス毎に、対象画素位置のＣＴ値の経時変化としてＴＤＣを生成する。

【0026】

出力関数化処理部１３７は、ＣＴ経時変化特定部１３３により生成されたＴＤＣに基づいて、出力関数等を特定する。具体的に、出力関数化処理部１３７は、ＣＴ経時変化特定部１３３で特定された経時変化であるＴＤＣに基づいて各ＣＴ値間の傾きの値を算出し、算出された傾きの値のうちで所定値以上の傾きに基づき出力関数を特定する。

【0027】

血流解析処理部１５０は、入力関数及び出力関数に基づいて、血流量の定量解析を行う。具体的に、血流解析処理部１５０は、入力関数データ記憶部１２９の入力関数テーブル７００及び出力関数データ記憶部１３９の出力関数テーブル７５０を参照して、所定の解析処理を行う。この血流量の定量解析の手法としては、例えば、Ｐａｔｌａｋ Ｐｌｏｔ法や、Ｄｅｃｏｎｖｏｌｕｔｉｏｎ法などを用いて、ＭＢＦ（Ｍｙｏｃａｒｄｉａｌ Ｂｌｏｏｄ Ｆｌｏｗ：心筋血流量）、ＭＢＶ（Ｍｙｏｃａｒｄｉａｌ Ｂｌｏｏｄ Ｖｏｌｕｍｅ：心筋血液量）、ＭＴＴ（Ｍｅａｎ Ｔｒａｎｓｉｔ Ｔｉｍｅ：平均通過時間）等を特定する。具体的には、後述する入力関数テーブル７００の関数Ｆの傾きａ又は関数Ｆのデータと、出力関数テーブル７５０の傾きα又は関数Ｆのデータと、に基づいて所定の解析処理が実行される。また、血流解析処理部１５０による定量解析の結果は、表示装置１７０に表示される。このとき、例えば、定量解析の結果をＣＴ画像と並べてまたは重ねて表示して、二つの画像を比較することにより、心筋梗塞または狭心症などの際の冠動脈狭窄の評価と心筋虚血評価を同時に行うようにすることができる。また、例えば、特定されたＭＢＦと正常なＭＢＦとを比較し、特定されたＭＢＦが正常なＭＢＦよりも低い場合には、虚血状態であるものと判断でき、延いては、心筋の冠動脈における心筋梗塞または狭心症などの際の冠動脈狭窄の評価と心筋虚血評価をすることができる。

【0028】

＜ＣＴ画像データの構成＞
図２を参照して、画像データ記憶部１１０に記憶されるＣＴ画像について説明する。図２は、ＣＴ画像の一例を例示する説明図である。図２に示したＣＴ画像は、心臓を心電図に同期して撮像したものであって、心拍が同位相の所定枚数（例えば３０枚）のフレーム画像２１０である。１フレームは、複数枚のスライス画像によって構成されているとともに、１フレームのデータは、３次元の画像データで構成されている。なお、図２では、１つのスライスのフレーム画像を例示している。ＣＴ画像は、被験者に造影剤を投与した直後から撮像が開始されたものでもよく、撮像順に第１フレームから第３０フレームまでが含まれてもよい。

【0029】

＜入力関数導出処理の内容＞
図３および図４を参照して、次に、入力関数導出部１２０によって実行される入力関数導出処理について説明する。図３は、入力関数導出処理を示すフローチャートである。図４は、入力関数特定処理を示すフローチャートである。

【0030】

図３に示すように、入力関数導出処理において、まず、ＲＯＩ設定部１２０は、解析者による選択処理に沿って、画像データ記憶部１１０から特定の領域が写っているスライス画像のフレーム画像の画像データを読み出し（Ｓ３１０）、選択されたスライス画像の中で、特定の領域が鮮明に写っているフレーム画像が表示装置５に表示されているとき、解析者の操作に従ってＲＯＩを設定する（Ｓ３３０）。次に、ＲＯＩ経時変化特定部１２３は、ＲＯＩが設定された全フレーム画像についてそれぞれＲＯＩ値を算出し、ＲＯＩ値に基づくＴＤＣを描画する（Ｓ３５０）。そして、入力関数化処理部１２５は、後述する入力関数特定処理を実行する（Ｓ３９０）。

【0031】

図４に示すように、入力関数特定処理において、まず、入力関数化処理部１２５は、ＲＯＩ経時変化特定部１２３により生成されたＴＤＣに基づいて、入力関数の導出対象とするフレームの範囲を規定する上限フレームＦａを特定し（Ｓ４０５）、フレームナンバーを示す変数ｎを１に設定する（Ｓ４１０）とともに、ｎに１を加算する（Ｓ４２０）。次に、入力関数化処理部１２５は、第ｎフレームより以前のフレームのＲＯＩ値を一つの直線Ｌに近似し（Ｓ４３０）、第ｎフレーム以降で上限フレームＦａ以前のフレームのＲＯＩ値を２次関数Ｆに近似する（Ｓ４４０）。次に、入力関数化処理部１２５は、直線Ｌと第ｎフレームより以前のフレーム のＲＯＩ値の最小二乗誤差を算出するとともに、２次関数Ｆと第ｎフレーム以降で上限フレームＦａ以前のフレームのＲＯＩ値との最小二乗誤差を算出して、これらの誤差の総和を算出する（Ｓ４５０）。そして、入力関数化処理部１２５は、ステップＳ４１０～Ｓ４５０までの処理を、ｎが上限フレームＦａ－１よりも小さい間は繰り返し（Ｓ４６０でＹＥＳ）、ｎ＝２から上限フレームＦａ－１までの最小二乗誤差の総和をすべて算出する。次に、入力関数化処理部１２５は、ｎ＝２から上限フレームＦａ－１までの最小二乗誤差の総和をすべて算出すると（Ｓ４６０でＮＯ）、その中から最小二乗誤差の総和が最も小さいｎを特定する（Ｓ４７０）。そして、入力関数化処理部１２５は、ＴＤＣにおけるＲＯＩ値が最大値ＭａｘとなるフレームをＭＴとし、ステップＳ４７０で特定されたｎのときの直線Ｌ及び２次関数Ｆの交点となるフレームをＡＴとし、ＡＴとしたフレームにおけるＲＯＩ値（直線Ｌの高さ（Ｙ切片））をＢＬとし、ＭＴの前後の所定の心拍数に相当する期間を撮像区間ＭＢとするとともに、入力関数および入力関数の傾きａを特定する。

【0032】

＜出力関数導出処理の内容＞
図５および図６を参照して、出力関数導出部１３０において実行される出力関数導出処理について説明する。図５は、出力関数導出処理を示すフローチャートである。図６は、出力関数特定処理を示すフローチャートである。

【0033】

図５に示すように、出力関数導出処理において、まず、対象画素抽出部１３１は、一つのスライス画像を選択し、画像データ記憶部１１０からそのスライスのフレーム画像の画像データを読み出し（Ｓ５１０）、選択されたスライスの全画素の中から、全フレーム画像におけるＣＴ値が前述の条件を満たす画素が特定領域の画素とみなして、これらの画素を対象画素として抽出し（Ｓ５２０）、抽出された対象画素から一つの画素を選択する（Ｓ５３０）。次に、ＣＴ経時変化特定部１３３は、対象画素抽出部１３１によって選択された画素についてのＴＤＣを描画する（Ｓ５４０）。次に、出力関数化処理部１３７は、後述する出力関数特定処理をＣＴ値に基づき実行して出力関数等を特定する（Ｓ５５０）。そして、出力関数導出部１３０は、ステップＳ５３０～Ｓ５５０までの処理をすべての対象画素について行う（Ｓ５８０）とともに、ステップＳ５１０～Ｓ５８０までの処理をすべてのスライスについて行う（Ｓ５９０）。

【0034】

図６に示すように、出力関数特定処理において、まず、出力関数化処理部１３７は、ＣＴ経時変化特定部１３３によって生成されたＴＤＣにおける各プロット間の傾きαを算出し（Ｓ６１０）、算出された傾きにおいて所定値以上（例えば、０．２以上）となる傾きのみを抽出する（Ｓ６２０）。そして、出力関数化処理部１３７は、抽出された傾きαに基づき所定の関数で近似し（Ｓ６３０）、出力関数を特定する（Ｓ６５０）。

【0035】

＜入力関数テーブル７００および出力関数テーブル７５０の構成＞
図７（Ａ），（Ｂ）を参照して、入力関数データ記憶部１２９に記憶される入力関数テーブル７００および出力関数データ記憶部１３９に記憶される出力関数テーブル７５０のデータ構成について説明する。図７（Ａ），（Ｂ）は、入力関数テーブル７００及び出力関数テーブル７５０のデータ構造の一例を示した説明図であり、（Ａ）は入力関数データ記憶部１２９の入力関数テーブル７００のデータ構成の一例を示した説明図であり、（Ｂ）は出力関数データ記憶部１３９の出力関数テーブル７５０のデータ構成の一例を示した説明図である。

【0036】

図７（Ａ）に示すように、入力関数テーブル７００は、ＡＴ７１１と、ＭＴ７１２、ＢＬ７１３と、関数Ｆの傾きａ７１４と、関数Ｌ及びＦ７１５と、ＲＯＩ値７１７と、ＭＢ７１８と、をデータ項目として含んでいる。なお、関数Ｌ及びＦ７１５は、入力関数化処理部１２５で求められたｎ＝１から上限フレームＦａまでを近似した直線Ｌの式及び２次関数Ｆの式であり、ＲＯＩ値７１７は、ＴＤＣの元となるフレーム毎のＲＯＩ値である。また、血流解析処理部１５０は、入力関数テーブル７００における入力関数を血流解析に使用する場合に、ＡＴ７１１及びＢＬ７１５が原点となるようにＲＯＩ値７１７もしくは関数Ｌ及びＦ７１５を補正したものを入力関数として使用されるようにしてもよい。

【0037】

図７（Ｂ）に示すように、出力関数テーブル７５０は、スライス番号７６１と、画素位置７６２と、傾きα７６３と、関数Ｆ７６５と、をデータ項目として含んでいる。なお、画素位置７６２は、傾きα７６３を算出するための基となる２つのＣＴ値の位置を示すデータであり、関数Ｆ７６５は、傾きα７６３から算出された近似式である関数Ｆの式である。

【0038】

＜心臓のＣＴ画像に対する画像処理の内容＞
次に、本実施例における心臓のＣＴ画像の画像処理について具体的に説明する。本実施例における心臓のＣＴ画像に基づくＴＤＣに対する画像処理においては、心臓に対して、造影剤を生理食塩水で所定の濃度に希釈したもの（例えば、造影剤を生理食塩水で４倍に希釈したもの）を利用するテストインジェクション法での撮影を実施した後に、希釈しない造影剤を利用した撮影を実施する、という流れである。具体的に、まず、テストインジェクション法で撮影が実施されることによって、被験者の大動脈を含む心臓のＣＴ画像である第１ＣＴ画像が撮像される。そして、この第１ＣＴ画像に基づいて、画像処理装置１００による上述のような入力関数導出処理が実行されることによって、大動脈の第１ＴＤＣが生成されるとともに、ＭＢや入力関数等が特定される。次に、第１ＴＤＣに基づいて特定されたＭＢにおいて、希釈しない造影剤を利用した撮影が実行されることによって、被験者の大動脈および心筋内冠動脈を含む心臓のＣＴ画像である第２ＣＴ画像が撮像される。そして、この第２ＣＴ画像に基づいて、画像処理装置１００による上述のような出力関数導出処理が実行されることによって、心筋内冠動脈の第２ＴＤＣが生成されるとともに、出力関数等が特定される。そして、このように特定された入力関数や出力関数等に基づいて、画像処理装置１００によって心筋内冠動脈の血流量の定量解析が実行されることとなる。

【0039】

図８（Ａ），（Ｂ）を参照して、ここで、第１ＴＤＣおよび第２ＴＤＣに対する画像処理の概要について説明する。図８（Ａ），（Ｂ）は、本実施例におけるＴＤＣの一例を例示した説明図であり、（Ａ）は本実施例における大動脈の第１ＴＤＣの一例を示した説明図であり、（Ｂ）は本実施例における心筋内冠動脈の第２ＴＤＣの一例を示した説明図である。なお、図８（Ａ）は、心臓の大動脈近傍の第１ＣＴ画像において大動脈にＲＯＩを設定した場合のものであり、図８（Ｂ）は、対象画素を所定のＣＴ値（本実施例では、５０～１５０ＨＵとする。）に設定した場合のものである。

【0040】

図８（Ａ）に示すように、まずは、テストインジェクション法での撮影で得られた心臓の第１ＣＴ画像に基づく画像処理装置１００の画像処理によって、大動脈の第１ＴＤＣが生成されるとともに、この第１ＴＤＣに基づいてＭＢ等が特定される。具体的には、入力関数導出部１２０のＲＯＩ経時変化特定部１２３は、まず、ＲＯＩ設定部１２１で設定されたＲＯＩを、同一スライスの他のフレーム画像２１０にも適用し、フレーム画像毎に、ＲＯＩ内の画素のＣＴ値に基づいてＲＯＩを代表するＲＯＩ値を決定する。次に、ＲＯＩ経時変化特定部１２３は、同一スライスの全フレームについて定められたＲＯＩ値をプロットすることで、第１ＴＤＣを生成する。次に、ＲＯＩ経時変化特定部１２３は、第１ＴＤＣにおけるＣＴ値の最大値Ｍａｘ、またはカーブが大きく立ち上がった （ＣＴ値が急激に上昇した）後の最初のピーク値を検出する。そして、ＲＯＩ経時変化特定部１２３は、ＴＤＣにおけるＲＯＩ値の最小値Ｍｉｎから、この最大値Ｍａｘまたはピーク値までの一定割合（例えば、７０％、８０％、９０％などを設定することができ、本実施例では、７０％とする。）の値をとる。

【0041】

そして、入力関数化処理部１２５は、最大値Ｍａｘまたはピーク値が検出されたフレームよりも以前のフレームであって、上述のようにＲＯＩ経時変化特定部１２３によって特定された一定割合の値をとるフレームを上限フレームＦａとする。次に、入力関数化処理部１２５は、第ｎ（ｎは、２から上限フレームＦａ－１まで）フレームＦｎ以前のＲＯＩ値に対して最小二乗法等を適用して直線近似を行い、直線Ｌの式を導出するとともに、第ｎフレームＦｎから上限フレームＦａまでのＲＯＩ値に対して最小二乗法等を適用して２次関数に近似し、２次関数Ｆの式を導出する。次に、入力関数化処理部１２５は、直線Ｌと第ｎフレームＦｎ以前のＲＯＩ値との誤差の二乗和を算出するとともに、２次関数Ｆと第ｎフレームＦｎ以降のＲＯＩ値との誤差の二乗和（残差平方和）を算出する。そして、入力関数化処理部２５は、すべてのｎについて上記の処理を行い、これら誤差の二乗和の合計が最も小さくなるｎを特定し、この誤差が最小となるｎの場合の直線Ｌ及び２次関数Ｆをそれぞれ近似された入力関数として特定するとともに、入力関数の傾きａを特定する。ここで、ＲＯＩ値が最大値ＭａｘとなるフレームがＭＴとなり、直線Ｌ及び２次関数Ｆの交点（境界）となるフレームがＡＴとなり、直線Ｌ及び２次関数Ｆの交点（境界）となるフレームのＲＯＩ値がＢＬとなり、ＭＴの前後で３心拍に相当する区間がＭＢとなる。このように入力関数化処理部１２５は、第１ＴＤＣに基づいて、大動脈の領域のＭＴ、ＡＴ、ＢＬ、ＭＢ、傾きａ、入力関数を特定する。なお、上述のＢＬの特定において、ＲＯＩ値そのものをＢＬとする以外に、例えば、直線Ｌの高さ（Ｙ切片）をＢＬとしても良い。

【0042】

ここで、第２ＣＴ画像を得るための撮影を実施するタイミングとなる撮像区間であるＭＢについて説明する。一般的には、大動脈に造影剤が流れ込んだ後に心筋内冠動脈に造影剤が流れ込むこととなるため、ＭＴの前後であって所定の心拍数に相当する区間における撮影でＣＴ画像を取得すれば、造影剤が流れ込んだ心筋内冠動脈に対するＣＴ画像を取得することができる。特に、第２ＣＴ画像において、心筋内冠動脈に造影剤が確実に流れ込んだときのフレームにおける心筋内冠動脈のＣＴ値を得るために、本実施例における撮像区間ＭＢは、ＭＴの前の１心拍とＭＴの後の２心拍とに相当する区間としている。

【0043】

図８（Ｂ）に示すように、次に、ＭＢにおける撮影で得られた第２ＣＴ画像に基づく画像処理装置１００の画像処理によって、心筋内冠動脈の第２ＴＤＣが生成されるとともに、この第２ＴＤＣに基づいて出力関数等が特定される。ここで、心筋内冠動脈に血液が主に流れるのは、心拍動に伴う心臓の拡張期であるため、心臓の拡張期と非拡張期とではＣＴ値に変化が生じてしまい、第２ＴＤＣは段差が形成されたようなプロットとなってしまう。そこで、第２ＴＤＣに基づき出力関数を特定する場合には、所定値以上となる傾きの値に基づいて近似することによって精度の高い出力関数を特定することができる。具体的には、撮像区間ＭＢでの撮影に基づく第２ＴＤＣにおける各プロット間の傾きの値を算出して、全ての傾きの値のうち所定値以上となる傾きα１～α３の値を抽出し、この傾きα１～α３に基づく出力関数特定処理によって出力関数が特定されることとなる。

【0044】

ここで、上述のように心臓の拡張期と非拡張期との関係に沿った第２ＴＤＣにおける１心拍中の傾きの値の変化の態様は、前半は小さい傾きの値であり、後半は大きい傾きの値となるような態様である。そして、出力関数特定処理において利用する傾きの値は、１心拍中における後半の大きい傾きの値である傾きα１～α３であるものとしている。つまり、言い換えれば、出力関数特定処理は、第２ＴＤＣにおいて１心拍中の傾きのうち後半部分の傾きの値を利用した処理である。なお、出力関数を特定するときにおいて、傾きα１～α３の平均値を利用して出力関数を特定してもよいし、傾きα１～α３の各々の近似式を特定してから種々の特定方法によって出力関数を特定してもよい。

【実施例2】

【0045】

図９（Ａ），（Ｂ）を参照して、次に、実施例２における心臓のＣＴ画像の画像処理について具体的に説明する。図９（Ａ），（Ｂ）は、実施例２におけるＴＤＣの一例を例示した説明図であり、（Ａ）は実施例２における大動脈の第１ＴＤＣの一例を示した説明図であり、（Ｂ）は実施例２における心筋内冠動脈の第２ＴＤＣの一例を示した説明図である。なお、実施例２における画像処理装置１００の構成は実施例１と同様のものであり、実施例２における入力関数導出処理および出力関数導出処理のフローは実施例１と略同様のものであるため、＜画像処理装置１００の構成＞、＜入力関数導出処理の内容＞、および＜出力関数導出処理の内容＞等の説明は省略する。

【0046】

図９に示すように、実施例２における心臓のＣＴ画像の画像処理においても、上述の実施例１と同様に、心臓に対して、テストインジェクション法での撮影を実施した後に、希釈しない造影剤を利用した撮影を実施する、という流れである。但し、上述の実施例１とは異なり、実施例２における第１ＴＤＣに基づき特定されるＭＢは１心拍に相当区間であり、第２ＴＤＣも１心拍に相当するものが生成される。そして、撮像区間ＭＢでの撮影に基づく第２ＴＤＣにおける傾きα４を抽出し、この傾きα４に基づく出力関数特定処理によって出力関数が特定されることとなる。なお、実施例２における撮像区間であるＭＢは、心筋内冠動脈に造影剤が確実に流れ込んだときのフレームにおける心筋内冠動脈のＣＴ値を得るために、ＭＴの後の１心拍に相当する区間としている。

【0047】

＜実施例１および実施例２の特徴＞
以上のように、上述の実施例１および実施例２の画像処理装置１００によれば、Ｓ５４０において、撮像区間であるＭＢでの撮影による時系列の複数フレームの第２ＣＴ画像からなる画像データに基づきＣＴ値の経時変化である第２ＴＤＣを特定し、Ｓ６２０において、第２ＴＤＣにおける各プロット間のＣＴ値の傾きであって所定の傾きα１～α３、α４を抽出し、Ｓ６３０において、これら所定の傾きα１～α３、α４に基づき所定の関数で近似することができる。特に、撮像区間であるＭＢにおける算出された複数の傾きに対して、所定値より大きい傾きα１～α３、α４を抽出した処理を実行することができる。したがって、このような実施例１および実施例２の画像処理装置１００であれば、例えば、心筋内冠動脈のＣＴ画像に対する画像処理に利用すると、本来、心筋血流が流れていない心臓の非拡張期のデータを排除して、拡張期のデータに基づいた精度の高い画像処理、延いては、精度の高い血流量の定量解析を実行することができる。また、従来の画像処理に比べて簡便な方法で解析することができるため、画像処理による画像処理装置１００における負荷を軽減することもできる。

【0048】

また、上述の実施例１および実施例２の画像処理装置１００によれば、画像データ記憶部１１０は、第１所定数のフレームの第１ＣＴ画像からなる第１画像データと、第１所定数に比べて少ない第２所定数のフレームの第２ＣＴ画像からなる第２画像データと、を記憶しており、希釈した造影剤を利用したテストインジェクション法の撮影による第１ＣＴ画像から特定されたＭＴに基づいて、所定の心拍数に相当する撮像区間であるＭＢを特定し、ＭＢでの希釈していない造影剤を利用した撮影が実施されて、得られた第２ＣＴ画像に基づいて、画像処理装置１００による画像処理を実行することができる。したがって、従来のように２０～３０心拍に応じた期間のＣＴＰ検査をする場合に比べて、少ない心拍数に応じた期間である撮像区間であるＭＢでの撮影によるＣＴ画像で画像処理を実行することができるため、従来に比べて被験者における被ばく量を軽減するとともに簡便な検査に基づく画像処理を実行することができる。

【0049】

さらに、上述の実施例２の画像処理装置１００によれば、撮像区間であるＭＢに基づいた１回の拍動に相当するフレーム数の第２ＣＴ画像および第２ＴＤＣに基づいて、画像処理装置１００による画像処理を実行することができる。したがって、従来のように２０～３０心拍に応じた期間のＣＴＰ検査をする場合に比べて、少ない心拍数に応じた期間である撮像区間ＭＢのＣＴＰ検査によるＣＴ画像で画像処理を実行することができるため、従来に比べて被験者における被ばく量を軽減するとともに簡便な検査に基づく画像処理を実行することができる。

【0050】

＜その他の実施例＞
上述の実施例において、第１ＣＴ画像に基づくＭＴの前後で３心拍又は１心拍に相当する区間を撮像区間であるＭＢとしたが、本発明ではこれに限定されず、ＭＢを以下のようにしてもよい。例えば、ＭＴの前後で２心拍や４～１０心拍に相当する期間を撮像区間であるＭＢとしてもよい。このようなＭＢや撮像タイミングに基づく撮影であっても、従来に比べて被験者における被ばく量を軽減するとともに簡便な検査に基づく画像処理を実行することができる。

【0051】

上述の実施例において、図６に示した出力関数特定処理では、算出された傾きにおいて０．２以上となる傾きのみを抽出するものとしたが、本発明ではこれに限定されず、０．２以外の数値を閾値としてもよい。例えば、複数回の撮影による結果から算出された統計値に基づく数値を閾値としてもよい。このような閾値であれば、より精度の高い画像処理、延いては、精度の高い血流量の定量解析を実行することができる。

【0052】

上述の実施例において、図６に示した出力関数特定処理では、算出された傾きにおいて０．２以上となる傾きのみを抽出し、つまり、心臓の拡張期における傾き（１心拍中の傾きのうち後半の傾き）の値のみを抽出するものとしたが、本発明ではこれに限定されず、心臓の非拡張期におけるＣＴ値の変化に基づき出力関数を特定するようにしてもよい。具体的には、例えば、第２ＴＤＣにおける０．２未満となる傾きを抽出し、抽出された傾きに基づいた各ＣＴ値の変化（心臓の非拡張期におけるＣＴ値の変化）の差分に基づく差分法等によって出力関数を特定する、といった処理であってもよい。

【0053】

上述の実施例において、図６に示した出力関数特定処理は、出力関数導出処理において利用するようにしたが、本発明はこれに限定されず、入力関数導出処理において図６に示した出力関数特定処理のような処理を利用してもよい。例えば、第１ＴＤＣに基づく入力関数を特定するときに出力関数特定処理のように傾きから入力関数を特定するようにしてもよい。このような画像処理であれば、従来の画像処理に比べてより簡便な方法で解析することができるため、画像処理による画像処理装置１００における負荷をより軽減することができる。

【0054】

上述の実施例において、第１ＣＴ画像を撮像するために、造影剤を生理食塩水で希釈した造影剤を利用したテストインジェクション法でＤｙｎａｍｉｃ撮影を実施していたが、本発明ではこれに限定されない。例えば、対象部位にＲＯＩを設定してリアルタイムでＲＯＩのＣＴ値（ＲＯＩ値）をモニタリングするボーラストラッキング法で撮影を実施してもよい。また、希釈していない造影剤を利用したテストインジェクション法でＤｙｎａｍｉｃ撮影を実施してもよく、さらに、テストインジェクション法とボーラストラッキング法とを組み合わせた、所謂、テストボーラストラッキング法でＤｙｎａｍｉｃ撮影を実施してもよい。つまり、上述の実施例において、第２ＣＴ画像を撮像するための撮影を実施するタイミングを最適化するためには、種々の手法を利用して第１ＣＴ画像を撮像することができる。

【0055】

上述の実施例においては、希釈した造影剤を利用したテストインジェクション法での撮影で第１ＣＴ画像を得て、その後に希釈していない造影剤を利用した撮影で第２ＣＴ画像を得る、といった手順であるものとしたが、本発明はこれに限定されない。例えば、第１ＣＴ画像を撮像するためのテストインジェクション法での撮影の実行中に、所定の心拍においてのみ複数心位相を撮影することによって第２ＣＴ画像を撮像する、といった手法であってもよい。

【0056】

上述の実施例において、対象画素抽出部１３１は、上述のような抽出処理を実行していたが、本発明ではこれに限定されず、画素単位に以下の処理を実行してもよい。例えば、対象画素抽出部１３１が、時系列の複数フレームのＣＴ画像において、各画素のＣＴ値の変化量、例えば、最大値と最小値との差分の値に基づいて画素を選択してもよい。具体的に、対象画素抽出部１３１は、一つのスライスにおいて、全フレーム画像の中でＣＴ値の最大値と最小値との差を求め、その差が所定の値（例えば５０～１５０）である画素を抽出してもよい。さらに、対象画素抽出部１３１は、このような抽出処理を上述の抽出処理に加えて実行してもよいし、いずれか一方のみを実行してもよい。

【0057】

上述の実施例において、入力関数特定処理では上限フレームＦａを上述のような方法で特定していたが、本発明はこれに限定されない。例えば、第１ＴＤＣの変化率を求めて、変化率が最大となるフレームと変化率が０となるフレームとの間のいずれかのフレームを上限フレームＦａとしてもよい。

【0058】

上述の実施例において、入力関数特定処理では直線と２次曲線との２つの関数で第１ＴＤＣを関数化しているが、本発明はこれに限定されない。例えば、第ｎフレームＦｎ以前を直線Ｌで近似し、第ｎフレームＦｎ以降のＲＯＩ値又はＣＴ値の分布を直線または３次以上の多次関数で近似してもよい。あるいは、上限フレームＦａまでのすべてのＲＯＩ値又はＣＴ値の分布を多次元多項式で表される関数に近似してもよい。また、第１ＴＤＣを３つ以上の関数で近似してもよく、例えば、第２フレームから上限フレームＦａ－１までを３つ以上の区間に分けて、それぞれの区間を所定の関数で近似してもよい。第１ＴＤＣにおいては、カーブが大きく立ち上がる（ＲＯＩ値又はＣＴ値が急激に上昇する）直前に、一旦ＲＯＩ値又はＣＴ値が低下する場合がある。このような現象が現れた場合には、ＲＯＩ値又はＣＴ値がほぼ一定で直線近似できる区間（第１区間）と、ＲＯＩ値又はＣＴ値が減少する区間（第２区間）と、その後のＲＯＩ値又はＣＴ値が急激に上昇する区間（第３区間）とに分けてもよい。そして、例えば、第１区間を直線、第２区間を２次以上の関数、第３区間を別の２次以上の関数で近似するようにしてもよい。このような場合には、第１区間と第２区間の境界となるフレーム、または第２区間と第３区間の境界となるフレームの何れかをＡＴとし、このＡＴとなるフレームのＲＯＩ値又はＣＴ値をベース値としてもよい。

【0059】

さらに、第１ＴＤＣを一つの関数で近似するようにしてもよく、例えば、正規累積分布関数または累積分布関数を第１ＴＤＣにフィッティングさせて、第１ＴＤＣを単一の関数で近似してもよい。このように第１ＴＤＣを正規累積分布関数で関数化する場合、第１ＴＤＣの立ち上がりのカーブに最もよくフィットするように正規分布の標準偏差（ＳＤ）及び平均値を選択するようにしてもよい。このとき、例えば、第１ＴＤＣに近似された正規累積分布関数の－３ＳＤに最も近いフレームをＡＴとし、このＡＴとなるフレームのＲＯＩ値又はＣＴ値をベース値としてもよい。

【0060】

上述の実施例において、関数特定処理では平滑化した第１ＴＤＣに対して関数近似していたが、本発明はこれに限定されず、上述のようにｍ次関数にフィッティングする以外の方法で第１ＴＤＣが平滑化されてもよい。例えば、第１ＴＤＣを生成するときに、周辺の画素との平均値を求めてこれを利用して平滑化してもよいし、第１ＴＤＣを移動平均によって平滑化してもよい。また、平滑化していない第１ＴＤＣに対して関数近似してもよい。

【0061】

上述の実施例において、上述のように血流解析処理部１５０による定量解析の結果が表示装置１７０に表示されていたが、本発明はこれに限定されない。例えば、心臓領域を任意の複数のセグメントに分けて表示する画像を表示装置１７０に表示させるとき、各セグメントはそれぞれの血流量に応じた表示態様で表示されるようにしてもよい。あるいは、画像データ記憶部１１０に保存されている画像データに基づいた臓器の３Ｄ画像を表示させるときは、その３Ｄ画像における各画素の表示態様が血流量に応じたものになるようにしてもよい。なお、ＣＴ画像は、個々の画素の座標情報を基に３Ｄ画像表示を行ってもよい。

【0062】

以上、実施例、変形例に基づき本発明について説明してきたが、上記した発明の実施の形態は、本発明の理解を容易にするためのものであり、本発明を限定するものではない。本発明は、その趣旨並びに特許請求の範囲を逸脱することなく、変更、改良され得ると共に、本発明にはその等価物が含まれる。

【符号の説明】

【0063】

１００…画像処理装置、１１０…画像データ記憶部、１２０…入力関数導出部、１２１…ＲＯＩ設定部、１２３…経時変化特定部、１２５…関数化処理部、１２９…入力関数データ記憶部、１３０…出力関数導出部、１３１…対象画素抽出部、１３３…ＣＴ経時変化特定部、１３７…出力関数化処理部、１３９…出力関数データ記憶部、１５０…血流解析処理部、１６０…入力装置、１７０…表示装置。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】