(57)【要約】電気的パラメータイメージング装置および方法は、イメージングされる対象を囲む電極アレイ上で電荷分布パターンを取得することを含む。さらに、対象の、外面境界、すなわち、輪郭が、光ビームおよび関連する光センサのアレイによって測定される。輪郭測定を使用して、取得された電荷分布パターンから、電気的パラメータ画像を計算するのに必要とされる対象の幾何形状の最初の評価値が提供される。
2006505352.tif 000002
【特許請求の範囲】
【請求項１】
対象の電気的特徴を示す画像を作成する方法であって、
ａ）前記対象の表面と電気接続するセンサエレメントのアレイにより前記対象の表面に電圧を印加するステップと、
ｂ）前記印加された電圧から生ずるそれぞれのセンサエレメントにおける表面電荷を測定するステップと、
ｃ）表面電荷測定値を変換するステップと、
ｄ）前記対象の表面輪郭を示す情報を使用して前記対象の前記電気的特徴を評価するステップと、
ｅ）前記変換された表面電荷測定値および前記対象の評価された電気的特徴から、前記対象の全体にわたる電気的特徴値を計算するステップと、
ｆ）前記計算された電気的特徴値から前記対象の前記電気的特徴の画像を作成するステップとを含む対象の電気的特徴を示す画像を作成する方法。
【請求項２】
前記ステップｃ）がフーリエ変換である請求項１に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項３】
前記ステップｄ）が、
ｉ）前記対象の表面輪郭を測定すること、および、
ｉｉ）前記表面輪郭測定値を使用して、前記対象の物理的モデルを作成することを含む請求項１に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項４】
前記ステップｄ）が、
ｉｉｉ）前記変換された表面電荷測定値を使用して、前記対象の異常の存在を判断することをさらに含む請求項３に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項５】
前記ステップｅ）が、
ｉ）前記対象の前記評価された電気的特徴によって作成される評価された表面電荷パターンを計算すること、
ｉｉ）前記評価された表面電荷パターンをフーリエ変換すること、
ｉｉｉ）前記フーリエ変換された表面電荷パターンを前記フーリエ変換された表面電荷測定値と比較すること、および、
ｉｖ）前記対象の前記評価された電気的特徴を変更し、前記フーリエ変換され評価された表面電荷パターンと前記フーリエ変換された表面電荷測定値の間の差が所定の量より小さくなるまで、前記ステップｉ）、ｉｉ）、ｉｉｉ）を反復することを含む請求項２に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項６】
前記ステップｄ）が、
ｉ）前記対象の表面輪郭を測定すること、および、
ｉｉ）前記表面輪郭測定値および評価された電気的特徴の対応する複数のセットを使用して、前記対象の複数の物理的モデルを作成することを含む請求項５に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項７】
前記ステップｅ）が、
ｖ）前記評価された電気的特徴のセットのそれぞれについて、前記ステップｉ）、ｉｉ）、ｉｉｉ）、およびｉｖ）を反復すること、および、
ｖｉ）前記フーリエ変換され評価された表面電荷パターンと前記フーリエ変換された表面電荷測定値の間の最小の差を生ずる前記評価された電気的特徴を選択することをさらに含む請求項６に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項８】
前記対象がヒトの乳房であり、ステップｄ）が、
ｉ）前記乳房の輪郭を測定することを含む請求項１に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項９】
前記乳房の輪郭が、光ビームのアレイを確立し、前記乳房を前記光ビームのアレイ内に配置し、どの光ビームが前記乳房によって遮断されるかを判定することによって測定される請求項８に記載の対象の電気的特徴を示す画像を作成する方法。
【請求項１０】
対象から電気的特性データを取得する装置であって、
前記対象を収容する空洞を画定する容器と、
前記容器によって支持され、前記対象の周りに配設される金属電極のアレイと、
前記容器によって支持され、前記空洞を通って延びる光ビームを生成するように動作可能である発光デバイスのアレイと、
前記容器によって支持され、前記光ビームのうちのそれぞれ１つの光ビームを受け取るように配置された光センサのアレイと、
前記光センサによって生成された信号に応答する手段であって、それによって、前記空洞内に設置された前記対象の輪郭を求める、応答する手段と、
前記金属電極のアレイに電圧を印加する手段と、
それぞれの金属電極において生成される、結果として得られる電荷を測定する手段とを備える対象から電気的特性データを取得する装置。
【請求項１１】
既知の電気的特性を有する流体が、前記容器内に含まれ、前記対象の表面および前記金属電極と接触する請求項１０に記載の対象から電気的特性データを取得する装置。
【請求項１２】
前記流体が、前記対象の前記電気的特性に一致する生理食塩水(saline  water  solution)である請求項１１に記載の対象から電気的特性データを取得する装置。
【請求項１３】
該装置の温度をほぼ一定温度に維持する温度コントローラを含む請求項１０に記載の対象から電気的特性データを取得する装置。
【請求項１４】
前記流体をほぼ一定温度に維持する温度コントローラを含む請求項１１に記載の対象から電気的特性データを取得する装置。
【請求項１５】
前記対象がヒト被験者の乳房であり、前記容器が、支持構造に取り付けられ、それによって、前記ヒト被験者が前記支持構造によって支持されると、前記乳房が前記空洞内で下に延びる請求項１０に記載の対象から電気的特性データを取得する装置。
【発明の詳細な説明】
【関連出願の相互参照】
【０００１】
本出願は、２００２年１１月７日に出願され、「METHOD  AND  APPARATUS  FOR  PRODUCING  AN  ELECTRICAL  PROPERTY  IMAGE  OF  SUBSTANTIALLY  HOMOGENEOUS  OBJECTS  CONTAINING  INHOMOGENEITIES」と題する米国仮特許出願第６０／４２４，５６８号の利益を主張する。
【発明の背景】
【０００２】
本発明は電気的イメージング技術、さらに詳しくは、腫瘍などの異質部分を除いて組成が実質的に均一である、ヒトの乳房などの対象の正確な電気的特性値により、高解像度画像を作成するための装置および方法に関する。
【０００３】
新しい医療イメージングモダリティに対する需要は、既存のイメージングモダリティを使用しては現在識別可能ではない組織特徴を識別する必要によって高まっている。肺がんに次いで、乳がんは、依然として女性にとって最も致命的ながんであり、米国国立がん研究所(National  Cancer  Institute)によれば、２００１年に約４０，２００人の女性の命を奪った。２００１年に１９２，０００の新たな乳がんの事例が存在した。米国の約２千８百万人の女性が、毎年、乳がんのスクリーニングを受けている。
【０００４】
スクリーニング段階では、高い割合の乳がんが検出されない。調査が示すところでは、乳がんの２０％〜５０％が、スクリーニング段階で検出されないままとなる。早期検出に対する動機は大きく、早期段階で検出される乳がんは、１１，０００ドルの平均治療費および約９６％の５年生存率を有するが、末期乳がんは、治療するのに平均で１４０，０００ドルかかり、５年生存率は２０％に低下する。医療専門家は、がん性組織を判定するために費用のかかる生検に頼ることが多い。これらの手技は、速くも、患者に優しくもない。がん性腫瘍の放射線治療は、がん性細胞の完全な破壊を確実に行うために、腫瘍領域の全体にわたって広範囲に、かつ、過剰に適用される。明らかに、乳がんを検出し、放射線治療処置中におけるがん性細胞破壊をリアルタイムで追跡するための、より優れたイメージング技術の必要性が存在する。
【０００５】
Ｘ線マンモグラフィは、乳がん検出のための好ましいモダリティである。デジタルシステムの開発、および、ｘ線フィルムを走査することによって、医師が疑わしい病巣を識別するのを助けるコンピュータ支援診断（ＣＡＤ）の使用に伴い、マンモグラフィシステムの販売の大幅な増加が期待される。しかしながら、先に述べたように、多数のがんは、ｘ線マンモグラフィを使用して検出されず、ｘ線被爆を低減するために、乳房圧迫技法が使用され、検査を苦痛なものにしている。
【０００６】
疑わしい病巣が見つかった後、標準的な処置は、生検を実施することである。外科生検は、悪性の可能性が高い疑わしい病巣について推奨されるが、穿刺吸引(fine-needle  aspiration)細胞診（ＦＮＡＣ）およびコア生検は、費用がかからず、有効な代替である可能性がある。ＦＮＡＣとコア生検は共に、外科生検の数を減らし、患者の不安をなくし、処置のコストを低減するのに役立ってきた。しかしながら、コア生検は、浸潤(invasive)がんを示すことにしばしば失敗してきた。また、ＦＮＡＣとコア生検は共に、悪性細胞を標的から離れたところに変位させる可能性があり、誤診をもたらす。
【０００７】
米国がん協会(American  Cancer  Society)によれば、乳房生検の約８０％は良性である。このため、テラヘルツパルスイメージング(ＴＰＩ)、赤外を含む熱および光学イメージング技法、蛍光および電気インピーダンスイメージングを含む、新しい侵襲性の少ない技術が開発されてきた。ほとんどの場合、これらの技術は、コンピュータ断層撮影法、磁気共鳴イメージング、陽電子放出断層撮影法、超音波、および、ＰＥＴ−ＣＴなどのハイブリッドシステムを含む、従来のイメージングモダリティに次ぐものとして追求されている。
【０００８】
正常細胞に対するがん性細胞の生化学的特性は、３つの因子、すなわち、ナトリウム、カリウム、および他のイオンの細胞内含量の増加、水の細胞内含量の増加、および、細胞膜の電気化学的特性の著しい差を特徴とする。ナトリウム、カリウム、および他のイオンの細胞内濃度の増加によって、細胞内の導電率が高くなる。同様に脂肪質の細胞が、がん性細胞を取り囲む時に、水が脂肪より良好な導電体であるために、水含量の増加によって導電率が高くなる。そして、最後に、がん性細胞の細胞膜における生化学的差によって、誘電率が高くなる。
【０００９】
乳がんの研究では、３つの別個の組織の分類、すなわち、腫瘍バルク(tumor  bulk)、浸潤性縁部(infiltrating  margins)、および遠方（正常）組織(distant(normal)tissue)が説明されている。病巣の中心は、腫瘍バルクと呼ばれ、高い割合のコラーゲン、弾性線維、および多くの腫瘍細胞を特徴とする。腫瘍細胞がわずかであること、および大部分が、影響を受けていない乳房組織内で正常に分布したコラーゲンと脂肪であることが浸潤性縁部を特徴付ける。最後に、遠方組織（病巣から２ｃｍ以上離れている）は、正常組織とみなされる。
【００１０】
がん性組織の特徴付けは、２つのグループ、すなわち、 部位内(in  situ：非浸潤性)病巣および浸潤性病巣に分けられる。部位内病巣は、病巣がそこから発生した上皮組織内に閉じ込められたままになる腫瘍である。腫瘍は、基底膜を超えず(cross)、したがって、腫瘍および健康な組織は、同じ性質である（上皮性）。そのため、部位内病巣の電気インピーダンスは、巨視的な導電率（ナトリウムおよび水の増加によって影響を受ける）および誘電率（細胞膜の電気化学的性質の差によって影響を受ける)に影響を与えることになる悪性細胞の量によって決まる。
【００１１】
対照的に、浸潤性病巣は、基底膜を通過する腫瘍である。悪性組織は、正常組織と異なる性質を有する（上皮対脂肪）。上皮組織は、ぎっしり詰まり、密である。脂肪組織は、ほとんどがトリグリセリドである大細胞で構成される。これらの構造上の差は、以下の影響を与える。第１に、正常組織は細胞密度が低い。第２に、正常組織の細胞液は、上皮細胞ほどには豊富ではない。一般に、上皮細胞の半径は、脂肪細胞の半径より小さく、このことから、本発明者らは、がん性細胞の半径は正常細胞の半径より小さいと結論付ける。がん性細胞対正常細胞の分数容積(fractional  volume)に与える影響は、がん性細胞の分数容積が、正常細胞より大きいことである。その理由は、上皮の固体群(population)が、正常な脂肪細胞より大きいためである。最後に、がん性細胞の細胞内導電率は、正常細胞の細胞内導電率より大きいことに、本発明者らは留意する。同様に、（正常細胞とがん性細胞の間の大きな空間のために）細胞外の流体が豊富なため、細胞外導電率は大きい。そのため、浸潤した組織の導電率は、正常組織より大きいことになる。
【００１２】
１９５０年代以来、一部の研究者は、生物学的組織の電気的特性を測定し、表にしてきた。ヒトの組織の電気的特性（導電率および誘電率）は、周波数依存性（分散)を示す。Ｄ．Ｃ．〜１ＧＨｚの範囲の周波数において、３つの分散領域（α、β、およびγ）が存在する。組織におけるこれらの分散は、細胞数、細胞形状、および細胞の向き、ならびに、組織の化学組成（すなわち、間隙性空間および細胞質の組成およびイオン濃度)によって決まる。
【００１３】
種々の調査が示すところでは、生物学的組織の固有抵抗値は多くの理由で変わる。たとえば、がん性腫瘍は、周囲の健康な組織に比べて、２桁大きい（１００倍の）導電率および誘電率を有する。医療治療の適用によって、組織の電気的特性の変化も生ずる。放射線によって治療された筋肉組織の場合、組織インピーダンスに対する大きな変化が報告されている。大きな変化は、ハイパーサーミア治療中の、低周波での骨格筋の電気インピーダンスにおいて起こり、電気的特性のこの変化は、細胞壊死の始まりを予示する。
【００１４】
電気インピーダンス断層撮影法（ＥＩＴ）は、対象内のインピーダンス分布をマップするプロセスである。このマップは通常、電流の印加およびその対象の境界に沿う電位差の測定から作成される。３つの種類のＥＩＴシステム、すなわち、電流注入デバイス、印加式電位デバイス、および誘導デバイスが存在する。HendersonおよびWebsterは、インピーダンス分布の全体のマップを作成した、インピーダンスカメラとして知られるデバイスを初めて導入した。Ｓｈｅｆｆｉｅｌｄのシステムおよびそれを具体化したものは、第１世代のＥＩＴシステムである。８０年代の後期に、LiおよびKrugerは、誘導電流デバイスについて報告している。こうしたシステムでは、コイルの組み合わせ体が、試験下の対象の周りに設置される。コイル内の変化する電流によって、変動する磁界が生成され、磁界は次に、試験下の対象に電流を誘導する。他の駆動方法と同様に、電極を、対象の境界上に設置して、境界に沿う電位降下を測定する。
【００１５】
そのような電気的特性イメージング技術は、しばしば、「インピーダンス断層撮影法」といわれる。最も慣用的な電気的特性イメージング技法は、１）複数の電極またはセンサを測定すべき試料に直接付着させる（医療適用では、試料はヒトの身体である）および、２）電流が順次各電極を通じて試料に注入され、引き続いての電圧が測定されることを前提としている。したがって、これらの慣用的イメージング技法は「定電流／測定電圧」スキームを行う。
【００１６】
そのような慣用的電気的特性イメージング技法から離れて、本発明者らのうちの一人は、米国特許第４，４９３，０３９号に開示されているように、測定すべき対象の外部のアレイにセンサを配列した。さらに、試料のイメージングの間に、電流を測定しつつ、交流電圧を固定された振幅にて印加した。対象とセンサアレイの間の空間をインピーダンス整合媒体で満たすことによって、係属中の特許出願ＷＯ９９／１２４７０に記載されているように、この手法はさらに改良された。加えて、測定された表面電荷に基づく内部電荷分布を計算するための２つの技法が記載されており、これは、スケールファクター法および反復法と言われる。反復法およびスケールファクター法は、Ｘ線ＣＴシステムのような関連イメージングシステムから得られた内部構造の幾何学形状を初めに判断することを必要とする。また、反復法は、各領域の電気的特性を最初に仮定することも要し、その後、仮定の有効性をチェックするために、境界における予測された電流の順方向計算を行う。このプロセスを仮定が測定された値に近い境界電流を生じるまで反復する。スケールファクター法は、「ルックアップ」テーブルまたは神経ネットアルゴリズムを作成し、これにより、モデル計算の莫大なデータ組を用い、試料の電気的特性または内部を外的に測定されたパラメータと相関させることが可能となる。モデルの制限、そしてデータ組のサイズを合理的に維持するために結果を外挿する必要性のため、スケールファクター法は精度が制限されるが、試料のおおよその電気的特性を前もって知っておく必要はない。事実、スケールファクター計算の結果は、反復法の初期評価値となり得る。両法はコンピュータにより集中的に行われる。
【発明の概要】
【００１７】
本発明は、過剰な計算時間を必要とすることなく、また、ＣＴまたはＭＲＩ画像に基づく初期評価を必要とすることなく、対象の電気的特性の正確な画像を作成する装置および方法を提供することによって、従来の電気的パラメータイメージング技法に関連する問題を解決する。
【００１８】
より詳しくは、ほぼ均一な構造を有する、ヒトの乳房などの対象をイメージングする時に、既知の電圧を印加した時に境界に蓄積した電荷に加えて、対象の外側境界の場所が検出される。対象の外側境界の形状が分かれば、対象の全体にわたる特性の初期評価を使用して、電気的特性のイメージング手順を使用することができる。腫瘍などの異常は、センサ上で、蓄積された電荷に対する特別な寄与として見られ、これらの特別な電荷の寄与を使用して、内部の異常の位置が特定され、そのサイズおよび電気的特性を測定することができる。
【００１９】
本発明の一般的な目的は、個別の断層撮影画像を必要とすることなく、電気的特性画像を作成することが可能である、電気的特性イメージング方法およびシステムを提供することである。対象の周辺境界の場所は、電気的特性測定センサと交互に配置した光源および光学センサのアレイを使用して検出される。こうして、対象境界は、位置が特定され、電気的特性センサによって行われた電荷測定と共に自動的に登録される。
【００２０】
本発明の、上述の、また、他の目的および利点は、以下の説明から明らかになるであろう。説明において、説明の一部を形成し、本発明の好ましい実施形態が例として示される、添付図面が参照される。しかしながら、こうした実施形態は、必ずしも本発明の全範囲を示しておらず、したがって、本発明の範囲を解釈するために、添付の特許請求の範囲および本明細書が参照される。
【好ましい実施態様の説明】
【００２１】
本発明のイメージング技法の基礎となる数学的理論を図２〜図４を参照して説明する。図２は、領域Ｆ  ２０６が試料２０４と表面２０２の間のスペースであるように、表面２０２によって囲まれた、試料２０４および内部領域Ｆ  ２０６を含む閉鎖体積空間（closed  volume  space）２００の平面図である。試料２０４は、サブ領域Ａ  ２０８、サブ領域Ｂ  ２１０、サブ領域Ｃ  ２１２、サブ領域Ｄ  ２１４およびサブ領域Ｅ  ２１６と便宜的に表された複数の結合されたサブ領域を含む。各サブ領域２０８〜２１６は、ヒト被験者の異なる組織のような異なる材料から成るものであってもよい。
【００２２】
或る特定の周波数（ｆ）における電磁場を閉鎖体積空間スペース２００中の試料２０４に印加すると、例えば各サブ領域２０８〜２１６の誘電率および伝導率の電気的特性に違いがある場合の、試料２０４の各サブ領域２０８〜２１６の間の境界のように、電気的特性が変化する場所でのみに全電荷は生じる。これらの全電荷は、今度は、閉鎖体積空間２００の表面の全電荷の再分布を誘導する。これらの誘導された電荷分布は、（独立して自由に動く）見掛電荷ならびに閉鎖体積空間２００の表面２０２に位置する分極電荷双方に由来すると評価される。表面２０２の電荷は、全（見掛＋分極）電荷でもあり、ここに、表面２０２上の或る一点の全電荷は大文字「Ｑ」で示され、他方、閉鎖体積空間２００の内部の一点の全電荷は小文字「ｑ」で示される。全電荷Ｑの測定は、電荷Ｑの実際の測定、または電荷Ｑの経時的変化の率である電流の小さな増分ｌに由来する電荷Ｑのいずれかを含むことができることに注意することが重要である。
【００２３】
表面２０２上の一点における全電荷Ｑ、および内部の一点における全電荷ｑは電磁気学理論を介して関連付けることができる。時間により変化する電場が電気媒体に適用されると、それらは媒体中に電流を誘導する。これらの電流は、今度は、ファラデーの法則に従い、かけられた電場に加えて、電場を誘導することができる時間により変化する磁場を生じる。電場に対するこの余分な寄与は低周波数においてはわずかであり、無視することができる。本発明者らはこのいわゆる「準静的」近似を用いる。静電気学の基本的理論によれば、内部全電荷ｑおよび表面２０２上の全電荷Ｑは一意的に関連する。
【００２４】
図３は、表面２０２上の一点３０２における全電荷Ｑと内部の一点における全電荷ｑの間の、グリーン関数を介して関連付けられる関係を示す閉鎖体積空間２００の平面図である。具体的には、グリーン関数は、表面２０２上の点ｊにおける全電荷Ｑを点ｉにおける内部全電荷ｑと関係づける。
【００２５】
【数１】
2006505352.tif 000003
【００２６】
この関係は、試料２０４の内部サブ領域２０８〜２１６の電気的特性についての所望の情報を提供する。図３は、この議論で用いる座標系および関連幾何図形的配列のいくつかを示す。場の点３０４、源の点３０６および表面の点３０２についての座標系で用いる表記は、各々、
【数２】
2006505352.tif 000004
、
【数３】
2006505352.tif 000005
プライム
【数４】
2006505352.tif 000006
、および
【数５】
2006505352.tif 000007
ダブルプライム
【数６】
2006505352.tif 000008
である。源の点３０６における試料２０４の内部の全電荷ｑを、表面の点３０２における全電荷Ｑと関連付けることによって、試料２０４の内部の増強された画像を発生させることができる。したがって、電場が測定される位置は場の点３０４である。
【００２７】
本方法のイメージング技法は、従来の電気的特性イメージング技法とはかなり異なる。図４は、そのような従来のイメージング技法によって測定される閉鎖体積空間２００の平面図である。試料の電気的特性が集中回路エレメントのネットワークによって表される。そのような方法では、電流は、閉鎖体積空間２００の表面２０２の、既知の場所、例えばＰ１  ４０２に注入され、既知の場所、例えばＰ２  ４０４で引き出される。次いで、周囲のセンサ上の電圧を測定し、集中回路インピーダンスを電流−電圧測定値の組から計算する。これに対し、本発明の技法では、閉鎖体積空間２００における電場の波様性質、および測定容量を完全に記載することが可能となり、測定すべき試料２０４のサブ領域２０８〜２１６を特徴付けるのに用いられる集中回路エレメントネットワークまたは同等回路の構造に関して、いずれの具体的仮定も要しない。
電磁波理論のマックスウェルの方程式をすぐ上に記した問題に適用すると、式１Ａが得られる。
【００２８】
【数７】
2006505352.tif 000009
【００２９】
ここで、
σ＝伝導率
εｒ＝比誘電率
ε０＝自由空間の誘電率
Ф＝ポテンシャルである。
加えて、電磁気理論の標準的な結果は、ポアソン方程式として知られる、ポテンシャル（Ф）、および全電荷密度ρの間の関係となり、即ち、式１Ｂである。
【００３０】
【数８】
2006505352.tif 000010
【００３１】
ここで、ρＴｏｔａｌは体積全電荷密度である。電場Ｅは以下の方程式から得られる。
【００３２】
【数９】
2006505352.tif 000011
【００３３】
式１Ａおよび１Ｂからは、スカラーポテンシャルファイ（Ф）、重要な電荷密度を全電荷、すなわち、見掛電荷＋分極電荷に関連付けることが示される。
【００３４】
電気的特性をイメージングするための他の方法は、測定値から直接各領域の誘電率および伝導率を計算することを試みる。全内部電荷を中間ステップとして計算する。伝導率または誘電率を直接的に求めるよりもむしろ電荷を求めることの１つの利点は、電気的映像を全体的に支配する内部電荷が、実質的に、対象内の非連続部分において存在する境界のみに出現し、かくして、計算を要する数値がかなり少ないことである。以下の式２は、このことを示しており伝導率の勾配および誘電率の勾配は全電荷密度に寄与している。したがって、全電荷は、伝導率および誘電率が距離に従って変化する率に依存する。
【００３５】
【数１０】
2006505352.tif 000012
【００３６】
電磁気学的理論における標準的な定理は一意性の定理である。準静的場合についての一意性の定理は、もしポテンシャルまたはその通常の導関数が閉鎖体積を囲む表面で既知であれば、場の点３０４におけるポテンシャルは一意的に決定することができることを述べる。ポテンシャルまたはポテンシャルの通常の導関数の双方が知られている必要はないことに注意することが重要である。事実、もし双方が知られていれば、問題の解は重複するであろう。境界表面の或る部分での既知ポテンシャルおよび他の部分での通常の導関数に関する問題を定義することができるが、以下の式（３）は、表面２０２上でのポテンシャルが既知である単純な場合を考える。これは、ジリクレ境界条件として知られている。
式３はグリーンの関数を用いてのポワソンの方程式（式２）に対する解である。
【００３７】
【数１１】
2006505352.tif 000013
【００３８】
ここで、ＧＤがジリクレグリーン関数であり、ｄτは体積の要素であり、Ｓは体積τを囲う表面の要素である。
【００３９】
式３は、まさに一意性定理が予測するごとく、内部の全電荷ｑおよび表面２０２上のポテンシャルによって決定される場の点３０４におけるポテンシャルである。解は、標準的な処理である幾何関数、グリーンの関数により得られる。試料２０４が存在すると、全電荷ｑ密度の体積積分および表面２０４の表面積分が存在する。もし表面の同一ポテンシャル分布を、試料が存在しない場合で考えれば、電荷密度はゼロであるが、表面積分は依然として同一である。表面項（式３における２番目の積分）は試料２０４を挿入しても不変である。なぜならば、表面２０２上では電圧は前もって決定された値に設定され、試料２０４の挿入前後でその値に維持されるからである。このため、２つの項を差し引くと、残りの式はグリーンの関数（これは、測定センサのアレイの一定の形状に対して既知数である）および電荷密度のみを含む。したがって、センサの間に試料２０４が挿入された場合と、試料２０４が挿入されない場合との間のポテンシャルの差を用いるのが便利である。このポテンシャルの差は、ポテンシャルの差の通常の導関数を取ることによって表面２０２の電荷と関連付けることができ、これは、電場の垂直の成分を生じさせる。というのは、ガウスの法則により、導電性表面近くの場の垂直成分はその表面の面積当たりの電荷に直接比例するからである。次に、連続体モデルから離散電荷にわたる和に変更し、下記式（４）は、添字「ｊ」で表される表面２０２における電荷Ｑｊが、「ｊ」および「ｉ」双方を含む行列（matrix）の要素によって、添字「ｉ」で表される内部の電荷ｑｉに関連付けられることを示すものであって、結合行列要素は、単に、グリーンの関数の通常の導関数である。
【００４０】
【数１２】
2006505352.tif 000014
【００４１】
以下の式５は、「ｊ」における一連の式が、ベクトルとしての表面２０２の電荷を含む行列式にて一緒に記載して、グループ化できることを示すもので、ベクトルの各項は全電荷のうちの１つである。表面２０２の電荷については、大文字「Ｑ」が用いられ、それらは、各項が内部の全電荷のうちの１つである小文字「ｑ」を用いた同様のベクトルに関連付けられる。
【００４２】
【数１３】
2006505352.tif 000015
【００４３】
この一連の式は内部電荷「ｑ」を与えるべく反転され、行列それ自体は逆行列を有する。さらに詳しくは、もし逆行列が存在するならば、グリーンの関数の導関数行列にその逆行列を掛け合わせると、式（６）に示すように単位行列が得られる。
【００４４】
【数１４】
2006505352.tif 000016
【００４５】
以下でより詳細に述べるように、イメージングされる被験者は、測定アレイ内に設置され、測定アレイは、被験者を通して電界
【数１５】
2006505352.tif 000017
を確立するように、所望の周波数で、かつ、１５ボルトｒｍｓ以下の正弦波電圧を、被験者の表面に印加することを可能にする。印加されたこの電界によって生ずる表面電荷Ｑｊが測定される。表面電荷測定は、印加される電界を異なる方向に向けて反復され、また、１０ＫＨｚ〜１０ＭＨｚの異なる周波数で反復されることができる。
【００４６】
以下の式７は、被験者の表面の全電荷Ｑｊについての式において現れる二重和をいかにして一つの和にするかを示す。一旦、（まさにサイン関数である）直交関数の完全な組においてグリーンの関数を展開すると、結果は、式７におけるサイン関数内に現れるパラメータ「Ｌ」にわたる和、及び式４に現れる電荷ｑｉにわたる和でもある。与えられた値「Ｌ」についての適当なサイン関数を掛け合わせ、測定アレイの一面を合計する結果、「Ｌ」にわたる和が排除され、それにより、残ったただ１つの項が得られる。この結果は、サインおよびコサインの関数の直交特性のため起こる。測定アレイの対応する反対側の測定の結果を付け加え、式８に示すように、フーリエ変換（サイン変換）についての「Ｌ」の値を与えるための式が得られることによって、精度はさらに改良することができる。
【００４７】
【数１６】
2006505352.tif 000018
【００４８】
手順が比較的簡単である。「Ｌ」の各値では、１つの式を作成することができ、その各々は「ｉ」によって表される電荷にわたる和を含み、その行列の要素は式８に示される。
【００４９】
下記の繰り返しプロセスを使用することにより、対象の内部における電荷ｑの正確な表示を作成することが可能である。一旦これがなされれば、式（３）で提示した既知のグリーンの関数の解を使用して対象の内部のいずれの場所におけるポテンシャルについても解を出すことに、内部電荷分布画像を用いることができる。一旦、対象の内部の各場所のポテンシャルが知られれば、式（１Ｃ）を用いてそれらのポテンシャルから電場を容易にもたらすことができる。したがって、内部の１つの地点からもう１つの地点に進むにつれ、これは、次いで、対象の内部における各点での電気的特性の見積もりをもたらすものであるが、電場の変化を得ることができる。
【００５０】
内部の測定された電荷を取り、対象の内部のポテンシャルおよび／または電場を生じさせるのにおける最終ステップは、電荷ｑをグリーン関数の解の式（３）に挿入することによって達成することができる。内部静電ポテンシャルに対する式には２つの項、Ａ）上記にてまさに計算した電荷ｑに関連する体積の積分；およびＢ）表面に設定された既知のポテンシャルに関連する表面にわたる積分、があることに留意する。最初の積分は、上記にて計算された電荷を用いてまさに得られる。第２の積分はまさに容易に得られる。なぜならば、表面センサにおいて設定されたポテンシャルは実験者によって決定され、知られているからである。したがって、全てが知られており、単純なプラグイン操作および既知の式によって計算して、式（３）に従って内部ポテンシャルを得ることができる。
【００５１】
一旦内部ポテンシャルが知られれば、いずれの箇所の電場も、式（１Ｃ）において与えられるように各方向のポテンシャルの変化の率から得ることができる。一旦電場が内部のどこかで知られれば、この対象内の各領域の電気的特性（σおよびεｒ）は、対象内の各境界を横切ってこの電場の垂直成分の変化から計算することができる。したがって、対象を囲む媒体の既知の電気的特性で開始することにより、対象における隣接領域の電気的特性を計算することができる。領域（ｎ）から領域（ｎ＋１）の各境界において、電場の正常成分の比率は、以下のように、領域（ｎ）および（ｎ＋１）における電気的特性に関連する。
【００５２】
【数１７】
2006505352.tif 000019
【００５３】
σおよびεrの値は被験者を囲む媒体について知られているので、式（１２）を用いると、次の領域（ｎ＝２）についての対応する値が得られるであろう。領域２および３の間の境界につき式を再び適用し、領域３についての値が得られ、対象全体の電気的特性が得られるまでこれらを続ける。たとえば、乳房をイメージングする時、周囲の生理食塩水媒体の電気的特性が既知であり、乳房の任意の異常領域の電気的特性を計算するプロセスが適用される。組織の性質および悪性腫瘍の存在を示すのは、これらの異常な領域の電気的特性である。
【００５４】
表面電荷データおよび被験者輪郭データを取得し、被験者の電気的特徴を示す画像をそこから作成するシステムが図５に示される。システムは、図７Ａおよび７Ｂにより詳細に示され、以下で詳細に説明される測定アレイ支持構造５００を含む。支持構造５００は、被験者の電気的特性にできる限り厳密に一致する、既知の電気的特性の生理食塩水で満たされた容器を形成する４つの垂直側面および底面を有する。イメージングされる被験者は、開いた上面５０２を通して挿入される。乳房をイメージングするのに使用されると、支持構造５００は、患者テーブルの開口の下に取り付けられ、乳房は、容器内に垂れ下がるように位置合わせされる。
【００５５】
システムは、図１により詳細に示され、以下で述べるコンピュータコントローラ５０４によって制御される。システムは、インピーダンスアナライザ５０６を動作させて、電圧ドライバ５０８を通して電荷測定アレイの個別のエレメントに電荷を印加し、これらのエレメントのそれぞれにおいて、結果として得られる電荷Ｑを測定する。インピーダンスアナライザ５０６は、「１２６０  インピーダンス／ゲイン位相アナライザ（１２６０  Ｉｍｐｅｄａｎｃｅ／Ｇａｉｎ  Ｐｈａｓｅ  Ａｎａｌｙｚｅｒ）」という商標名の下でSolartron  Analyticalから販売されている。インピーダンスアナライザ５０６は、コンピュータコントローラ５０４上で実行されるその「Ｚプロット」ソフトウェアを使用して動作する。
【００５６】
電圧ドライバおよび電荷センサは図６に詳細に示される。演算増幅器５１０は、その入力端子５１２と電荷測定アレイエレメントに接続される一対の出力５１４の間で単位ゲインを有するインバータとして動作する。直列接続された出力抵抗ＲＳの両端の電圧降下は、アナライザ５０６への出力として役立ち、出力５１４が接続される電荷測定アレイエレメントにおける、結果として得られる表面電荷Ｑｊを計算するのに使用される。
【００５７】
コンピュータコントローラ５０４はまた、走査回路５２０を動作させて、イメージングされている被験者の表面輪郭を指示するデータを取得する。走査回路５２０は、ＬＥＤドライバ回路５２２を順次イネーブルして、支持構造５００の２つの側面に沿って配設されるＬＥＤアレイのＬＥＤに電圧供給する。同時に、支持構造５００の２つの対向する側面に沿って配設される光センサアレイ５２６の、対応するイネーブルされた光センサから、信号が５２４で入力される。ＬＥＤドライバ５２２および対応する光センサ５２６の対を順次イネーブルすることによって、支持構造５００内の被験者の輪郭を求めることができる。走査回路５２０は、支持構造５００内の被験者によって遮断された光経路を指示するマップをコントローラ５０４に提供する。コントローラ５０４によって実施されるプログラムは、このマップデータを補間して、イメージングされている被験者の外部表面が画定される。
【００５８】
測定精度を維持するために、測定アレイ支持構造５００内の生理食塩水の温度を制御することが必要である。これは、温度コントローラ５０５によって達成され、温度コントローラ５０５は、温度センサ（図示せず）から受け取る信号に応答して支持部５００内の加熱エレメント（図示せず）を動作させ、温度センサは同様に支持部５００内にある。好ましくは、温度は、患者の快適さのために、体温に維持される。
【００５９】
特に図７Ａおよび７Ｂを参照すると、測定アレイ支持構造５００は、容器の４つの面の全てに配設された金属エレメント５５０の２Ｄアレイを含む。これらのエレメント５５０は、対応する電圧ドライバ５０８の出力５１４に接続される正方形金属電極である。エレメントは、被験者５５４を囲む生理食塩水媒体５５２と電気接触する。これらのエレメント５５０に印加される電圧は、容器内で、かつ、被験者５５４全体を通して電界Ｅを確立し、被験者５５４の電気的特性によって決まる表面電荷Ｑｊを蓄積する。好ましい実施形態２２５では、エレメント５５０は、４つの面のそれぞれに配設され、塩化銀コーティングを有する銀で構成される。
【００６０】
発光ダイオード５５８の２次元アレイは、支持構造５００の２つの隣接する壁に配設され、光センサ５６０の２次元アレイは、反対の２つの壁に配設される。それぞれのＬＥＤ５５８は、反対の壁の上に配設される対応する光センサ５６０と対をなして、容器を通る一意の光経路が画定される。上述したように、各対を順次イネーブルすることによって、容器を通る光経路はすべて、系統的に調べられて、どの光経路が被験者５５４によって遮断されたかが求められる。より正確に輪郭を測定するために、支持構造５００全体が回転して、複数の異なる視野から輪郭データが取得される。
【００６１】
特に図１を参照すると、コンピュータコントローラシステムは、記憶システム２３の一部を形成するメモリ２２に記憶されたプログラム命令を実行するプロセッサ２０を含む。プロセッサ２０は、Ｍｉｃｒｏｓｏｆｔ  Ｃｏｒｐｏｒａｔｉｏｎ  Ｗｉｎｄｏｗｓオペレーティングシステムの１つによって動作するように設計された市販デバイスである。プロセッサ２０は、全ての外部メモリ２２を操作するために、システム統合および統合メモリ管理回路機構(circuitry)を容易にするための、内部メモリおよびＩ／Ｏ制御部を含む。プロセッサ２０はまた、３２ビットＰＣＩバス２４との直接インタフェースを提供するＰＣＩバスドライバを含む。
【００６２】
ＰＣＩバス２４は、プロセッサ２０と複数の周辺コントローラカード間で３２ビットのデータを転送する業界規格バスである。これらは、ＣＤ  ＲＯＭドライブ２８およびディスクドライブ３０に対して高速データ転送を提供するＰＣＩ  ＥＩＤＥコントローラ２６を含む。グラフィックコントローラ３４は、標準ＶＧＡ接続部３６を通してＰＣＩバス２４をＣＲＴモニタ１２に連結させ、キーボードおよびマウスコントローラ３８は、キーボードおよびマウス１４を通して手入力されるデータを受け取る。
【００６３】
ＰＣＩバス２４はまた、インピーダンスアナライザ・インタフェースカード４０および輪郭測定インタフェースカード４２に接続する。インタフェースカード４０は、手順のデータ取得段階中に、インピーダンスアナライザ５０６に対してデータを結合させる。プロセッサ２０によって実行されるプログラムは、インピーダンスアナライザ５０６を制御して、電荷測定アレイに電圧が印加され、結果として得られる表面電荷Ｑｊを示すデータが入力される。インタフェースカード４２は、発光ダイオード（ＬＥＤ）のアレイを駆動する走査回路５２０に接続し、光センサの対応するアレイから応答信号を受け取る。プロセッサ２０によって実行されるプログラムは、インタフェースカード４２を通して走査回路５２０を制御して、支持構造５００内に設置された対象の輪郭を測定し、イメージングされている被験者の表面境界を指示するデータをインタフェースカード４２を通して入力する。
【００６４】
特に図８を参照すると、手順は、画像取得段階および画像再構成段階からなる。プロセスブロック６００で指示するように、画像取得段階の第１ステップは、被験者を所定場所に置かないで表面電荷データＱｊを取得することである。この「空の」表面電荷データは、再構成段階の間に必要とされ、上述したように、測定アレイ５００に対して選択された周波数の電圧を印加することによって取得される。有限時間間隔にわたって蓄積する、結果として得られる表面電荷Ｑｊが入力される。システムは、印加されるＥフィールドの全ての可能な向きで電荷データが採取されるように、決定ブロック６０１でループバックする。
【００６５】
被験者は、次に、プロセスブロック６０２で指示するように、測定アレイ支持構造５００内に挿入され、被験者の輪郭は、プロセスブロック６０４で指示するように測定される。上述したように、これは、被験者によって遮断される光経路を指示するマップの取得を含む。この入力データは、補間されて、被験者の外部表面の場所および輪郭が画定される。
【００６６】
次に、表面電荷データＱｊが、規定の周波数および規定のＥフィールドの向きで取得されるループに入る。表面電荷データＱｊは、電荷測定エレメント５５０に規定の周波数で電圧を印加し、それぞれのエレメント５５０で蓄積する電荷Ｑｊを読み取ることによって、プロセスブロック６０６で取得される。測定は、それぞれの規定周波数で反復される。決定ブロック６０８で指示するように、システムは、次に、これらの測定を他のＥフィールドの向きで反復するようにループバックする。プロセスブロック６１０で指示するように、電荷測定エレメント５５０に印加される電圧振幅は、被験者内で生成される電界Ｅの方向を再制御する(reorient)ために変えられる。３Ｄ情報を提供するために、異なるＥフィールドの向きの少なくとも３つの取得が必要とされ、最終画像のＳＮＲを改善するために、さらなる取得が得られる場合がある。最後のＥフィールドの向きについて表面電荷データが取得されたと決定ブロック６０８にて判断されると、プロセスブロック６１２で指示するように、画像再構成を始めることができる。
【００６７】
特に図９を参照すると、画像再構成は、プロセスブロック６２０で指示するように、取得された表面電荷データＱｊの空間フーリエ変換（ＦＴｍｅａｓ）を計算することによって始まる。プロセスブロック６２２で指示される次のステップは、被験者の、物理的な場所、形状、およびサイズ、ならびに、被験者内の任意の異常についての第１評価値を生成することである。この第１の物理的評価値（モデルｋ）は、取得した被験者表面輪郭データを使用し、「正常な」組織の電気的特性値を、被験者内部の全ての場所に割り当てることによって生成される。腫瘍などの異常は、被験者上の測定された電荷分布Ｑｊを、全体にわたる正常組織に基づいて計算された期待される表面電荷分布と比較することによって評価される。これらの期待される表面電荷分布と測定された電荷分布の差は、被験者の普通なら均一な電気的特性内において異常を特定するのに必要とされる情報を提供する。
【００６８】
プロセスブロック６２４で指示される次のステップは、被験者のこの最初の物理的モデル（モデルｋ）の可能な変動範囲を予め計算することである。これらの変動は通常、評価される異常の、サイズ、形状、および場所の変動を含む。以下で述べるように、これらのｋｍａｘ個の物理的モデルのそれぞれを使用して、最適な電気的特性画像が見つかるであろう。
【００６９】
次に、測定された電荷（ＦＴｍｅａｓ）に最もよく「合う」電気的特性画像が作成される反復プロセスが始められる。最初に、プロセスブロック６２６で指示するように、異常を含む、被験者組織の電気的値（αｅｓｔ）が、物理的モデルｋおよび測定された電荷Ｑｊを使用して評価される。次に、プロセスブロック６２８にて、物理的モデルｋおよび評価された電気的値αｅｓｔを使用して、期待される表面電荷Ｑｅｓｔが計算される。プロセスブロック６３０にて、評価される電荷Ｑｅｓｔの空間フーリエ変換（ＦＴｅｓｔ）が計算され、プロセスブロック６３２にて、これが、測定された電荷ＦＴｍｅａｓと比較される。差（δＬ）がプリセットされた閾値より小さいと決定ブロック６３４にて判断される場合、プロセスブロック６３６で指示するように、この物理的モデルｋについての評価された電気的値（αｅｓｔ）が、差の値（δＬ）と共に記憶される。そうでない場合、システムは、評価された電気的値（αｅｓｔ）をプロセスブロック６３８で指示される量δαだけ変えた後、ループバックし、プロセスを反復する。
この反復プロセスは、以下のように表現される。
【００７０】
【数１８】
2006505352.tif 000020
【００７１】
ＦＴｍｅａｓとＦＴｅｓｔの差を使用して、
【数１９】
2006505352.tif 000021
と書く、次元Ｎｍａｘ×Ｌｍａｘの数値計算された行列によって与えられるＦＴｅｓｔの変化のテイラー級数展開の第１項を使用して、電気的特定ベクトルδαに対する補正が計算される。ここで、Ｎｍａｘは、段落００５０及び００５１で述べた異なる領域の数であり、Ｌｍａｘは、フーリエ項の最大数である。特異値分解技法を利用する、通常、優決定の式のセットによって式（１０）を反転することによってδαが得られる。
【００７２】
決定ブロック６４０において決定されたように、物理的モデルｋのｋｍａｘ個の変動の全てが上述した方法で処理されて、それぞれのモデルについて評価された電気的値（αｅｓｔ）のセットが作成される。プロセスブロック６４２で指示するように、最も小さい差の値（δＬ）をもたらした、記憶された電気的値（αｅｓｔ）のセットが記憶された値の中で特定され、プロセスブロック６４４で指示するように、この情報を使用して、最終画像が作成される。
【００７３】
個別のイメージングシステムを使用する従来の技法に勝る本発明の利点は、比較的費用がからない光学アレイを使用して、物理的モデルｋを評価するのに必要とされる情報を取得することができることである。したがって、このシステムは、より望ましい乳がんスクリーニングデバイスである。しかしながら、本発明によって作成される画像が、可能性のある悪性疾患を暴露する場合、取得された表面電荷データＱｊを、高解像度ＭＲＩまたはｘ線ＣＴまたは超音波画像から作成された物理的モデルと組み合わせて使用して、電気的値のより改良された画像を作成することができる。
【図面の簡単な説明】
【００７４】
【図１】本発明を実施するのに有用な例示的なコンピュータシステムを示すブロック図である。
【図２】空間内の閉鎖体積の平面図である。
【図３】測定された外面全電荷Ｑｊと内部全電荷ｑｉの関係を示す、空間内の閉鎖体積の平面図である。
【図４】従来の電気的特性イメージング技法によって測定される空間内の閉鎖体積の平面図である。
【図５】本発明を使用する、電気的特性イメージングシステムの好ましい実施形態のブロック図である。
【図６】図５のシステムの一部を形成する電圧駆動式回路の回路図である。
【図７】図５のシステムの一部を形成する測定アレイ支持部の概略図である。
【図８】図５のコンピュータコントローラによって実施されるデータ取得プログラムのフローチャートである。
【図９】図５のコンピュータコントローラによって実施される画像再構成プログラムのフローチャートである。
【図１】
2006505352.tif 000022
【図２】
2006505352.tif 000023
【図３】
2006505352.tif 000024
【図４】
2006505352.tif 000025
【図５】
2006505352.tif 000026
【図６】
2006505352.tif 000027
【図７Ａ】
2006505352.tif 000028
【図７Ｂ】
2006505352.tif 000029
【図８】
2006505352.tif 000030
【図９】
2006505352.tif 000031
【国際調査報告】
2006505352.tif 000032
2006505352.tif 000033
2006505352.tif 000034