JP2010035806A - 生体情報イメージング装置、画像構成方法 - Google Patents

Abstract

【課題】取得する情報の不確実性が増大しても、アーティファクトの発生を可及的に抑制することのできる生体情報イメージング装置および画像構成方法を提供する。
【解決手段】生体情報イメージング装置が、生体１に光２を照射する光源４と、光吸収によって発生する超音波７を受信し信号を出力する超音波検出器８と、超音波検出器８から出力される信号からタイムドメイン法により第１の画像を再構成するタイムドメイン処理部６１と、超音波検出器から出力される信号からフーリエドメイン法により第２の画像を再構成するフーリエドメイン処理部６２と、第１の画像と第２の画像とから診断画像を生成する画像演算部６３と、を備える。
【選択図】図１

Description

本発明は、光吸収によって発生する超音波を受信して画像化する生体情報イメージング装置、およびその画像構成方法に関する。

従来、乳癌などの診断では、エックス線（マンモグラフィー）、超音波、ＭＲＩ（核磁気共鳴法）を用いたイメージング装置が多く使われている。近年、パルス光を生体内に伝播させ、その伝播光の吸収によって発生した光音響波を生体表面で検知し、生体内の初期発生圧力分布あるいは光吸収係数値分布をイメージングする光音響波イメージング装置が注目を浴びている（特許文献１）。

一般的に、光音響波イメージング（光音響波トモグラフィー）の利点としては、以下のようなものが挙げられる。

第一に、エックス線、超音波、ＭＲＩのような体内の形態のイメージングではなく、酸素代謝や脳の活性のような機能イメージングが可能であることが挙げられる。例えば、酸素化ヘモグロビンと脱酸素化ヘモグロビンでは光の吸収スペクトルが異なっているため、吸収スペクトルを測定することで、血液内の酸素飽和度を求めることができ、酸素代謝測定などを行うことができる。なお、酸素代謝がイメージングできれば、そこから癌の腫瘍や脳の活性部位を調べることができるなどの利点がある。

第二に、光はガンマ線などの放射線やエックス線などと比べ被曝がないため、非侵襲で繰り返し連続して診断が可能であることが挙げられる。

第三に、ＭＲＩやＰＥＴ（ポジトロン放射断層撮影）などと比較し、医療機器の小型化・低価格化が可能であることが挙げられる。

第四に、光拡散を検知し、生体内のイメージングを行う光拡散イメージング装置と比較すると、光音響波の生体内での散乱が小さいため、比較的、解像度を上げることが可能であるという利点がある。

生体内の光吸収体の画像を取得するためには、受信した光音響波信号から、画像を再構成する必要があり、この画像再構成の手法としては、タイムドメイン法とフーリエドメイン法といった各種手法が提案されている（非特許文献１）。
Jonathan I. Sperl et al., "Photoacoustic Image Reconstruction - A Quantitative Analysis", SPIE-OSA Vol. 6631, 663103, 2007 米国特許第５７１３３５６号明細書

光音響波イメージングにおいては、生体内での光拡散のために、生体内深部に到達するフォトンの数が大きく減少する。よって、光吸収体が生体内深部に存在する場合は、生体表面に配置した音波検出器で光音響波信号を観測するのが困難になる。そのため、たとえ人体に対して３６０度の全方位に音波検出器が設置されていても、信号を取得可能な部位は限定されることになり、取得される情報の不確実性が増すことになる。

また、測定対象部位によっては３６０度の全方位に音波検出器を設置できないこともあ
る。このような制約がある場合には、例えば人体の片側に音波検出器を配置し、信号を測定することになる。しかし、光吸収体と音波検出器との間の距離が大きくなるほど、音波検出器と光吸収体とがなす立体角は小さくなることになり、この場合も取得する情報の不確実性が増すことになる。

このように、取得する情報の不確実性が増加すると、画像再構成処理において実際には存在していない虚像（アーティファクト）が発生してしまい、診断精度が低下するという課題がある。

そこで、本発明は、取得する情報の不確実性が増大しても、アーティファクトの発生を可及的に抑制することのできる生体情報イメージング装置および画像構成方法を提供することを目的とする。

上記目的を達成するために本発明は、以下の構成を採用する。

本発明の第１態様は、生体に光を照射する光源と、光吸収によって発生する超音波を受信し信号を出力する超音波検出器と、前記超音波検出器から出力される信号からタイムドメイン法により第１の画像を再構成するタイムドメイン処理部と、前記超音波検出器から出力される信号からフーリエドメイン法により第２の画像を再構成するフーリエドメイン処理部と、前記第１の画像と前記第２の画像とから診断画像を生成する画像演算部と、を備える生体情報イメージング装置である。

本発明の第２態様は、生体に光を照射し、光吸収によって発生する超音波を受信し、画像化する生体情報イメージング装置における、画像構成方法であって、前記超音波を受信して得られる信号からタイムドメイン法により第１の画像を再構成する工程と、前記超音波を受信して得られる信号からフーリエドメイン法により第２の画像を再構成する工程と、前記第１の画像と前記第２の画像とから診断画像を生成する工程と、を備える画像構成方法である。

本発明によれば、人体に対して３６０度の全方位から音波を検出できない場合など、取得する情報の不確実性が増した場合であっても、アーティファクトの発生を可及的に抑制できる。

次に、本発明の実施の形態について図面を参照しながら説明する。まず最初に、本発明の実施形態で使用する画像再構成手法について説明する。以下、説明を容易にするために２次元の場合で概要を説明する（３次元の場合も同様に考えることができる）。

画像再構成手法は、大まかにはタイムドメイン法とフーリエドメイン法とに分類することができる。

タイムドメイン法は、音波が球面波であるという性質を利用して、実空間において音波源の候補を描画していく方法である。代表的な手法として、ＤｎＳ（Delay and Sum）法
やＣＢＰ（Circular Back Projection）法、ＨＴＡ（Hough Transform Algorithm）法と
いった手法がある。

ＤｎＳ法による画像再構成においては、球面波をある一辺で観測した場合に、その時系列データの軌跡が、２次曲線をなすということを利用する。例えば、図４Ａの様に、横軸
に信号を観測する辺の座標を、縦軸に時間をとった時空間内においては、点音源の信号は符号１０１の様な２次曲線を描くことになる。

従って、その時空間内において、その２次曲線と相似な曲線上のデータを２次曲線の頂点に足し合わせていけば、音波源の候補が描画されることになる。

最後に、音速を用いて、時空間を実空間に変換すれば（ｚ＝ｃｔ、ｃ：音速）、図４Ｂの様に音波源１０２が描画される。

実際にＤｎＳ法を用いた画像再構成の例を示したのが、図４Ｃ、図４Ｄである。図４Ｃには、二つの円状音波源があるが、それを再構成すると図４Ｄの様な再構成画像が得られる。

このように、ＤｎＳ法は、２次曲線上の信号を足し合わせるという特性上、測定辺の近傍で測定辺に垂直なアーティファクト１０３が生じる傾向がある。

また、タイムドメイン法のひとつであるＣＢＰ法による画像再構成においては、図５Ａの様に、音波源の候補を実空間内において、信号の測定点を中心とする半径ｃｔ（ｃ：音速、ｔ：光を照射してからの経過時間）の円として足し合わせていく。この結果、音波源としての可能性の高い部分が描画されることになる。

実際にＣＢＰ法を用いた画像再構成の例を示したのが図５Ｂである。音波源としては、前記ＤｎＳ法の場合と同じ図４Ｃのデータを用いている。このＣＢＰ法を用いた場合には、円弧を描画するという手法の特性上、図５Ｂに示す様に、ＤｎＳ法と同様に信号を測定する辺に対し垂直方向のアーティファクト１０３が発生する傾向がある。

さらに、タイムドメイン法のひとつであるＨＴＡ法による画像再構成は、音波源を最初に描画する空間が時空間であるという点以外は、本質的にＣＢＰ法と同様である。そのため、ＣＢＰ法と同様に、信号を測定する辺に対し垂直方向のアーティファクトが発生する傾向がある。

一方、フーリエドメイン法（ＦＴＡ法）による画像再構成においては、フーリエ変換と音波の分散関係を用いて、測定された時空間の情報（図６Ａ）を、実空間の情報（図６Ｂ）に変換できることを利用する。すなわち、測定された音波の時系列に対し、まず順方向フーリエ変換を行い、時間スケールでのスペクトルを獲得する。その後、音波の分散関係（ω＝ｃｋ、ω：時間周波数、ｋ：空間周波数）から時間・空間スケールを変換した後に、逆方向フーリエ変換を行うことによって音波源の空間分布を画像として得る。実際にＦＴＡ法を用いた画像再構成の例を示したのが図６Ｃである。音波源としては、前記ＤｎＳ法の場合と同じ図４Ｃのデータを用いている。このＦＴＡ法を用いた場合に得られる画像には、図６Ｃに示す様な平行方向のアーティファクト１０４が発生する傾向がある。

前述した様に、タイムドメイン法は信号を測定する面に対し、垂直方向のアーティファクトを生じる傾向があるのに対し、フーリエドメイン法は平行方向のアーティファクトも生じる傾向がある。本発明者らはこのような傾向の違いに着目し、２つの画像再構成手法で得られた結果を組み合わせることでアーティファクトを減少させるという着想を得た。すなわち、フーリエドメイン法のみならずタイムドメイン法も併せて用いることにより、フーリエドメイン法のみを用いた画像よりもアーティファクトを減少させた画像が得られるということに至った。具体的には、本発明の生体情報イメージング装置では、タイムドメイン法により再構成された第１の画像とフーリエドメイン法により再構成された第２の画像とを演算処理（合成）することにより、アーティファクトの少ない診断画像を取得す
る。

一般に、タイムドメイン法による画像再構成処理は、フーリエドメイン法による画像再構成処理に比べて処理時間がかかる。そこで、タイムドメイン法による処理に用いる信号のデータ点数を（フーリエドメイン法による処理に用いるデータ点数よりも）少なくするとよい。これによりタイムドメイン法の処理速度を向上することができる。なお、タイムドメイン法により得られる画像は、主としてフーリエドメイン法で得られた画像のアーティファクトを減少させる目的で利用される。

［実施形態］
図１を参照して、本発明の実施形態に係る生体情報イメージング装置の概要について説明する。図１は、本実施形態に係る生体情報イメージング装置の構成例を示す図である。

本実施形態の生体情報イメージング装置は、生体内の光学特性値分布、特に光吸収係数値分布に関する情報、及び、それらの情報から得られる生体組織を構成する物質の濃度分布に関する情報の画像化を可能とするものである。この装置は、例えば、悪性腫瘍や血管疾患などの診断、化学治療の経過観察などに利用される。特に、この装置は、乳癌などの診断画像を得るための装置に好ましく適用することができる。

本実施形態の生体情報イメージング装置は、概略、光源４、超音波検出器８、信号処理装置６、及び、画像表示装置５を備える。

光源４は、被検体である生体１に照射するためのパルス光２を発生する手段である。なお、ここでいう生体とは乳房・指・足や首など人体の任意の場所を意味している。通常は、光源４から発生したパルス光は光ファイバや液体ライトガイドなどの光伝播装置３を通して、生体表面に照射される。ただし、パルス光のエネルギーが光ファイバなど光伝播路の限界を超える場合などは、ミラーやレンズを用いて、パルス光を生体に照射することも可能である。なお、生体表面での単位面積あたりの照射エネルギーを下げるために、パルス光は集光させるよりも広げる方が好ましい。典型的には数ｃｍ^２程度であるが、その面積は、パルス光のエネルギーと人体への光照射限界により最適な値に設計される。

超音波検出器８は、生体内の光吸収体９が光のエネルギーの一部を吸収して発生した超音波７を検出し電気信号に変換する手段である。光吸収体９は、例えば、腫瘍、血管、またはこれらに類する組織である。複数の超音波受信エレメントを有する超音波検出器８を用いることが好ましい。その場合は、２次元アレイ状に超音波受信エレメントを配列するとよい。信号処理装置６は、超音波検出器８で得られた電気信号を解析し、画像再構成処理を行う手段である。信号処理装置６は、タイムドメイン処理部６１とフーリエドメイン処理部６２と画像演算部６３とを備えている。画像表示装置５は、信号処理装置６の処理信号に基づいた画像を表示するための手段である。

（光源）
光源４は、生体を構成する成分のうち特定の成分に吸収される特定の波長の光を照射する。また、光源４は、数ナノから数百ナノ秒オーダーのパルス光を発生可能なパルス光源であることが望ましい。光源としてはレーザーが好ましいが、レーザーのかわりに発光ダイオードなどを用いることも可能である。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。光学特性値分布、特に、光吸収係数値分布の波長による違いを測定するために、単一波長の光のみ発振するものではなく、異なる波長の光を発振できる光源を用いても良い。そのときの光源としては、発振する波長を変更可能な色素、ＯＰＯ（Optical Parametric Oscillators）、チタンサファイヤ、アレキサンドライトなどの光学結晶を用いたレーザーなどが使
用できる。使用する波長に関しては、生体内において吸収が少ない７００ｎｍから１１００ｎｍの領域が好ましい。比較的生体表面付近の生体組織の光学特性値分布を求める場合は、上記の波長領域よりも範囲の広い、例えば４００ｎｍから１６００ｎｍの波長領域を使用することも可能である。

（超音波検出器）
超音波検出器８は、生体内を伝播した光のエネルギーの一部を吸収した光吸収体９が発生する超音波７を検出し、電気信号に変換する必要がある。そのため、超音波検出器８の受信可能な周波数帯は、生体内にある光吸収体９の大きさに応じて最適化されることが望ましい。超音波検出器８としては、圧電現象を用いたトランスデューサー、光の共振を用いたトランスデューサー、容量の変化を用いたトランスデューサーなど音響波信号を検知できるものであれば、どのような音波検出器を用いてもよい。様々な大きさの光吸収体から発生した超音波を受信する場合には、検出周波数帯の広い、容量の変化を用いたトランスデューサーが好ましい。

本実施形態では、一つの超音波検出器８を生体表面近傍に配置した例を示しているが、このような配置に限らず、複数の個所で超音波が検知可能に構成されていればよい。すなわち、複数の個所で音波を検知できれば同じ効果が得られるため、１個の音波検出器を生体表面上で２次元に走査しても良いし、１次元や２次元に超音波受信エレメントが配列したアレイ型を複数の場所に配置しても良い。２次元アレイ状の超音波検出器は、複数の位置で同時に超音波信号を検出できるため、スキャン方式の検出器よりも測定時間を短くできる利点がある。

なお、超音波検出器８から得られた電気信号が小さい場合は増幅器を用いて、信号強度を増幅することが好ましい。また、超音波検出器と測定対象である生体物質との間には、超音波の反射を抑えるための音響インピーダンスマッチング剤を使うことが望ましい。

（信号処理装置）
信号処理装置６は、超音波検出器８から出力される電気信号から画像構成を行い、これにより生体の光学特性値分布に関する情報（特に光吸収特性値分布に関する情報）を導くことができる。

本実施形態の信号処理装置６は、図１に示すように、タイムドメイン法による画像再構成処理を実行するタイムドメイン処理部６１と、フーリエドメイン法による画像再構成処理を実行するフーリエドメイン処理部６２とを備える。信号処理装置６は、超音波検出器８より得られた電気信号を元に、タイムドメイン処理部６１とフーリエドメイン処理部６２のそれぞれで画像を構成する。そして、画像演算部６３が、両処理部６１、６２から得られた画像を演算し、その結果を画像表示装置５に出力する。

信号処理装置６は、少なくとも、超音波の強さとその時間変化を記憶する記憶部と、この記憶部を参照して、超音波検出器８から得られた信号を光学特性値に関連したデータに変換する演算部とを具備する装置により実現可能である。例えば、データを解析できるコンピュータとソフトウエア、データを解析できるハードウエア、或いはその両方により、信号処理装置６を構成可能である。

（画像表示装置）
画像表示装置５としては、超音波信号を解析して得られた画像を表示できればよく、液晶ディスプレイなど様々なものを使用することができる。

なお、複数の波長の光を生体に照射する場合は、それぞれの波長に関して、生体内の光
学特性値分布、例えば、光吸収係数値分布に関係する情報を算出できる。それらの値と生体組織を構成する物質（グルコース、コラーゲン、酸化・還元ヘモグロビンなど）固有の波長依存性とを比較することによって、生体を構成する物質の濃度分布を画像化することも可能である。

このような実施形態に示された生体情報イメージング装置を用いることで、従来技術よりも、よりアーティファクトの少ない診断画像を獲得することが可能になり、信頼性の高い診断を行うことが可能となる。

（画像構成方法）
次に、図３を用いて、本実施形態の生体情報イメージング装置による画像構成方法について説明する。図３は、生体内の光学特性値分布である光吸収係数値分布に関係する情報を画像化する処理のフローチャートである。以下の（Ｓ１）から（Ｓ５）は、図３に示されたフローチャートの各処理Ｓ１〜Ｓ５にそれぞれ対応させたものである。この番号順に各処理を説明する。

（Ｓ１）光源４がパルス光２を生体１の表面へ照射する。

（Ｓ２）生体内の光吸収体９が、パルス光のエネルギーの一部を吸収し、超音波を発生する。この超音波信号を超音波検出器８で検出する。通常、球状の光吸収体により発生した超音波信号は、横軸に時間、縦軸に超音波の強度をとると、図２に示すようなＮ字型の波形を示す。

（Ｓ３）タイムドメイン処理部６１が、（Ｓ２）で得られた信号から、上記説明したタイムドメイン法により画像を再構成する。タイムドメイン法としては、上記説明した、ＤｎＳ法、ＣＢＰ法、ＨＴＡ法等が選択可能である。

（Ｓ４）同様にフーリエドメイン処理部６２が、（Ｓ２）で得られた信号から、上記説明したフーリエドメイン法により画像を再構成する。

（Ｓ５）次に、画像演算部６３が、タイムドメイン法によって得られた画像とフーリエドメイン法によって得られた画像の両者を演算することにより、最終的な診断画像を出力する。画像の演算としては、例えば、和（ＳＵＭ）、積（ＰＲＯＤＵＣＴ）などが考えられる。

図７Ａ〜図７Ｅは、画像の一例を示す図である。図７Ａは、生体内に存在する２つの円状音波源（光吸収体）を模式的に示す図である。図７Ｂは、タイムドメイン法としてＤｎＳ法を用いて構成された画像であり、図７Ｃは、フーリエドメイン法を用いて構成された画像である。図７Ｂ、図７Ｃから分かるように、タイムドメイン法で得られる画像とフーリエドメイン法で得られる画像とでは、発生するアーティファクトが異なる傾向がある。そこで、２つの画像を演算する（組み合わせる）ことにより効果的にアーティファクトを消去あるいは低減することができる。図７Ｄは、２つの画像の和演算の結果を示している。音波源に比べて、アーティファクトが目立たなくなっていることが分かる。また図７Ｅは、２つの画像の積演算の結果を示している。音波源の輪郭が明確になり、効果的にアーティファクトが消去されていることが分かる。なお、タイムドメイン法としてＣＢＰ法もしくはＨＴＡ法を用いた場合でも、同様の効果が得られる。ＣＢＰ法およびＨＴＡ法で発生するアーティファクトの性質はＤｎＳ法のものと類似しているからである。

本実施形態では、図３のフローチャートに基づいて、生体内の光エネルギー吸収密度分布、あるいは光吸収係数値に関連した特性値分布を求める例を説明したが、本発明はこの
フローチャートの処理に限定されるものではない。本発明の本質は、アーティファクトの性質の異なる複数の画像再構成手法を用いて超音波信号から複数の画像を生成し、その複数の画像を組み合わせることによりアーティファクトを低減することにある。すなわち、本実施形態では２つの画像から最終的な診断画像を構成したが、３つ以上の画像から診断画像を構成してもよい。また画像の演算としては、和と積に限らず、他の演算処理を適用してもよい。

（実施例）
本発明の一実施例を説明する。生体中の光吸収係数値の分布に関する情報を表す画像を求めることを目的とした生体情報イメージング装置の例について説明する。なお、ここでは生体を模擬したファントムをイメージングした実験例を示す。

光源としては１０６４ｎｍのナノ秒パルス光を発振できるＱスイッチＮｄ：ＹＡＧレーザーを用いた。パルスの幅は約５ナノ秒であり、繰り返し速度は１０Ｈｚである。例えば、１パルスのエネルギーは１２０ｍＪである。このパルス光を誘電体反射ミラーを用いて、生体の表面まで導く。なお、ここで用いたファントムはイントラリピッドを寒天で固めたものの中に、光吸収体として１ｍｍ直径のゴムチューブを埋め込んだものである。

光照射後に超音波検出器で測定された超音波信号は、コンピュータで構成された信号処理装置に送られる。信号処理装置は、タイムドメイン法とフーリエドメイン法とを用いて２つの画像を再構成し、２つの画像の積を診断画像として出力する。ここで、タイムドメイン法としては、ＤｎＳ法を用いた。

このような装置を用いて、ファントム全体の光吸収エネルギー密度値分布を画像化したところ、フーリエドメイン法を単独で使用した場合に比べ、よりアーティファクトの少ない光吸収体の画像を獲得することが可能になった。

図１は、本発明の実施形態における生体情報イメージング装置の構成例を説明する図である。 図２は、生体情報イメージング装置において検出される超音波信号の一例を示す図である。 図３は、生体情報イメージング装置による画像構成方法の一例を示すフローチャートである。 図４Ａ〜図４Ｄは、タイムドメイン法（ＤｎＳ法）による画像構成の概要を説明する図である。 図５Ａおよび図５Ｂは、タイムドメイン法（ＣＢＰ法）による画像構成の概要を説明する図である。 図６Ａ〜図６Ｃは、フーリエドメイン法による画像構成の概要を説明する図である。 図７Ａ〜図７Ｅは、本発明の効果を示す図である。

符号の説明

１：生体
２：パルス光
３：光伝播装置
４：光源
５：画像表示装置
６：信号処理装置
７：光吸収体から発生した超音波
８：超音波検出器
９：光吸収体
６１：タイムドメイン処理部
６２：フーリエドメイン処理部
６３：画像演算部
１０１：点音源の信号
１０２：音波源
１０３、１０４：アーティファクト

Claims (12)

1. 生体に光を照射する光源と、
   光吸収によって発生する超音波を受信し信号を出力する超音波検出器と、
   前記超音波検出器から出力される信号からタイムドメイン法により第１の画像を再構成するタイムドメイン処理部と、
   前記超音波検出器から出力される信号からフーリエドメイン法により第２の画像を再構成するフーリエドメイン処理部と、
   前記第１の画像と前記第２の画像とから診断画像を生成する画像演算部と、
   を備えることを特徴とする生体情報イメージング装置。
2. 前記タイムドメイン処理部は、前記フーリエドメイン処理部が前記第２の画像を再構成するために用いる信号よりも少ないデータ点数の信号から、前記第１の画像を再構成することを特徴とする請求項１に記載の生体情報イメージング装置。
3. 前記画像演算部は、前記第１の画像と前記第２の画像との和または積をとることにより、前記診断画像を生成することを特徴とする請求項１または２に記載の生体情報イメージング装置。
4. 前記タイムドメイン処理部は、ＤｎＳ法により前記第１の画像を再構成することを特徴とする請求項１〜３のうちいずれか１項に記載の生体情報イメージング装置。
5. 前記タイムドメイン処理部は、ＣＢＰ法により前記第１の画像を再構成することを特徴とする請求項１〜３のうちいずれか１項に記載の生体情報イメージング装置。
6. 前記タイムドメイン処理部は、ＨＴＡ法により前記第１の画像を再構成することを特徴とする請求項１〜３のうちいずれか１項に記載の生体情報イメージング装置。
7. 生体に光を照射し、光吸収によって発生する超音波を受信し、画像化する生体情報イメージング装置における、画像構成方法であって、
   前記超音波を受信して得られる信号からタイムドメイン法により第１の画像を再構成する工程と、
   前記超音波を受信して得られる信号からフーリエドメイン法により第２の画像を再構成する工程と、
   前記第１の画像と前記第２の画像とから診断画像を生成する工程と、
   を備えることを特徴とする画像構成方法。
8. 前記第１の画像を再構成する工程では、前記第２の画像を再構成するために用いられる信号よりも少ないデータ点数の信号から前記タイムドメイン法により前記第１の画像が再構成されることを特徴とする請求項７に記載の画像構成方法。
9. 前記診断画像を生成する工程では、前記第１の画像と前記第２の画像との和または積をとることにより、前記診断画像が生成されることを特徴とする請求項７または８に記載の画像構成方法。
10. 前記第１の画像を再構成する工程では、ＤｎＳ法により前記第１の画像が再構成されることを特徴とする請求項７〜９のうちいずれか１項に記載の画像構成方法。
11. 前記第１の画像を再構成する工程では、ＣＢＰ法により前記第１の画像が再構成されることを特徴とする請求項７〜９のうちいずれか１項に記載の画像構成方法。
12. 前記第１の画像を再構成する工程では、ＨＴＡ法により前記第１の画像が再構成されることを特徴とする請求項７〜９のうちいずれか１項に記載の画像構成方法。

Images