JP2000504411A

JP2000504411A - 流体が通過するシステムをモニタする装置

Info

Publication number: JP2000504411A
Application number: JP9526230A
Authority: JP
Inventors: デガニ，ハダッサ
Original assignee: イエダリサーチアンドデベロツプメントカンパニーリミテツド
Priority date: 1996-01-18
Filing date: 1997-01-21
Publication date: 2000-04-11
Also published as: IL116810A0; EP0877917A1; WO1997026507A1; AU721771B2; CN1134747C; CN1213430A; CA2242599A1; CA2242599C; EP2339303A3; AU1703397A; EP0877917A4; EP2339303A2

Abstract

(57)【要約】流体が通過するシステムをモニタする方法および装置において、この流体は、このシステムにおける時間空間での変化によって特徴付けられている。このシステム中の予じめ選択された場所をモニタ（１４）して、システムイベントに相関した、２つまたはそれ以上の時刻（１２）においてデータを収集する。このデータは、システムパラメータを表わし、このシステムパラメータは、システムウォッシュインおよびシステムウォッシュアウトビィフェイビア（１５）に関連した少なくとも２つの変数の関数として、時間と共に変化する。各ピクセルまたはボックセルを有する演算結果に基いてキャリブレーションマップを作成（２０）する。このピクセル、またはボックセルは、ウォッシュアウトビィフェイビア（１９）を表示する色相を表わすと共に、ウォッシュインビィフェイビア（１６）を表示する色強度を表わしている。満足できるマップが得られた場合には、収集したデータをこのマップを基準として処理して、色相および色強度と相関性を有する各空間ユニットを保有した予じめ選ばれた場所のイメージを得る。ソフトウエアおよびデータ処理システムを設けて、このキャリブレーションマップを生成する。このキャリブレーションマップおよび予じめ選ばれた場所のイメージは、また、新規な手段である。

Description

【発明の詳細な説明】流体が通過するシステムをモニタする装置発明の分野本発明は、時間空間で、システムをモニタする装置および方法に関し、特に、以下のような新規な装置に関するものである。即ち、ユニークなＭＲＩマシーンや、かんがい（灌漑）システムの効率をテストする新規なかんがい装置や、暖房および冷房システムの効率を決定する新規なテスト装置で、流体が通過するシステムをテストまたは制御する装置であり、ここで、システムの流体は、部分的に消費するか、または再生を必要としている。従来例の説明現在、流体が通過するシステムをモニタしたり、テストしたり、測定する種々の装置が知られている。また、この流体中の物質の一部分は、このシステムを横切る時に消費したり、再生を必要としたりする。例えば、今日、ＭＲＩマシーンを用いて、トレーサ造影剤を注入したり、注入しないでイメージを生成している。慣例上、このＭＲＩマシーンを制御して、個別の時間間隔で、一連のイメージを採取しており、次に、これらイメージを動的に分析して、出力を得ている。例えば、胸部腫瘍におけるコントラスト（造影剤）強調による動的検査により、信号強度における変化率は、胸部塊りの識別用の重要なパラメータであり、これによって薬物動態学が発展している。しかし乍ら、腫瘍の不均一性の結果として、コントラスト強調における時間変化の局部的な大きな変動が存在している。その結果として、記録および分析ステップの両者における高い空間解像度を維持することは、極めて重要なことである。胸部の標準的な医学的なＭＲＩでは、高い空間解像度を達成すると共に、高い時間解像度を維持することは、困難である。以前から実行されてきた動的検査の大部分においては、この強調は、高い時間解像度（空間解像度を犠牲にして）によるものであり、これは細胞内空間における平衡をモニタすると共に、組織の細胞外空間に、初期に分散するものであった。従って、標準的なＭＲＩマシーンでは、出力結果が、時として、想像外のものとなってしまう。前述したことは、他のシステムの特性でもある。これら他のシステム中では、流体が通過するか、または、この流体の成分が、一部分消費するか、または、再生を要求するようになり、例えば、かんがいシステム、暖房および冷房システム等が相当する。発明の要旨従って、本発明の目的は、流体が通過する、時間空間的にシステムをモニタする装置、またはマシーンおよびそれの方法を提供するものであり、これらによって、システム偏差またはシステム効率に関連した、更に相像可能な結果が得られる。また、本発明は、時間空間におけるシステムをモニタする装置に関する。このシステムとしては、物理的、化学的、生物学的、生理学的、環境的、医学的、または、他のシステムが部分的、または全体的に考案でき、このシステムは、所定方法において、時間空間的に徐々に変化している。本発明の装置は、１次、２次、３次またはそれ以上の高次元を基準にして機能できるものである。この装置によって選択される空間解像度のタイプや範囲、時刻数、およびその時間間隔は、このシステムに依存するものであると共に、２つの時刻に対する低い限度値と共に変化することができる。例えば、放射線検査の時間依存データの処理に利用でき、例えば、ＭＲＩ、超音波ホログラフィ、Ｘ線断層診断、従来のＸ線または核医学診断、予後、および腫瘍の治療や、更に、他の病理学的異常診断に利用できる。また、用水かんがいの環境データのモニタおよびコントロールにも利用できる。また、データを分析して、パイプ内の漏洩エリアの決定が行なえる。更にまた、食物、化粧や他の産業で得られるデータを分析することもでき、これらデータには混合物および溶液が包含されており、それらの均一性を決定できる。また、暖房および／または冷房システムの効率を評価できる。新規なアプローチを駆使して本発明による処理方法において、時間と共に、空間に亘って徐々に変化する現象には多数の現象が存在しており、このことを、本明細書では、ウォッシュイン（ｗａｓｈ−ｉｎ）およびウォッシュアウト（ｗａｃｈ−ｏｕｔ）ビィフェイビア（振舞）と称するものとする。これらウォッシュインおよびウォッシュアウトの用語は、一方向に変化すること（ウォッシュイン）および逆方向に変化すること（ウォッシュアウト）をシンボル的に記述したものであり、これらは、実際に逆のものでないが、変化を伴うあらゆる通路である。特に、流体または流体成分が消費するか、または再生要求が生じるシステムにおける流体の流れは、本発明においては、ウォッシュインおよびウォッシュアウトと記述してある。あらゆるウォッシュイン／ウォッシュアウト状況の下では、“ｎ”が１からあらゆる整数２，３，４等に変化する場合には、ｎ個の時刻を基準にして、ｎ個のウォッシュアウトパターンを記述でき、“ｍ”が２から、プロセスの時間変化における特別な時刻のあらゆる整数値に変化する場合にも記述できる。このウォッシュアウトの定義は、厳密なものでなく、ウォッシュインは、ウォッシュアウトにもなれると共に、逆にもなることが可能である。各種のシステムに対して、本発明による装置によれば、このシステムの時刻および他の最適なパラメータを設定する手段によって、このシステムをモニタしたり、コントロールしたり、調整したりする手段が得られる。このパラメータ等の設定には、物理的モデルに基いた新規なキャリブレーションマップが利用され、このマップによって、近似な手法または正確な手法で、時間を伴う変化について記述する。また、これらキャリブレーションマップによって、この装置製品の一つとして得られた最終色相／色強度を、定量的に解釈することもできる。この新規な装置の１つの特定な利用法としては、ガンの特別な診断を行なう製品を得るために、造影剤強調されたＭＲＩデータ用の装置が存在する。造影剤投与の開始時間は、時刻ｔ０であり、続く、２つの造影剤時刻ｔ１とｔ２とを利用する。これら造影剤時刻は、造影剤（コントラスト）強調の運動のモデルに基いてキャリブレーションマップを作成することによって選択される。この造影剤強調によって、ウォッシュイン／ウォッシュアウトレートを２つの病理学的パラメータに関連付けている。即ち、微小血管透過率×表面積（短かく、微小血管透過率と称し、文字Ｋで表わす）、およびνで表わされた細胞外ボリュームの有理数である。このキャリブレーションマップは、この装置によって形成される。図面の簡単な説明第１図は、時間空間において、システムをモニタするための、本発明による新規な装置のブロックダイヤグラム；第２図は、データ収集用パラメータを設定するための、第１図の装置の選択手段の詳細なブロックダイヤグラム；第３図は、第１図の装置の選択手段の着色機能Ｃを、詳細に表現するブロックダイヤグラム；第４ａ図および第４ｂ図は、第１図の装置の選択手段による、データ収集用のパラメータ設定を説明するフローチャート；第５ａ図および第５ｂ図は、ＭＲＩイメージ収集用の新規なＭＲＩ装置として使用する新規な装置の選択手段の選択を説明するフローチャート；第６図は、第１図装置における、システムに関するデータ収集の制御および実行する装置の制御／実行手段のブロックダイヤグラム；第７図は、上記装置をＭＲＩに適用した時のブロックダイヤグラム、および、被検体に関するＭＲＩデータ収集を制御／実行する装置の制御／実行手段の詳細を示すブロックダイヤグラム；第８図は、第１図装置における、収集したデータを処理する処理手段の詳細を表わすブロックダイヤグラム；第９図は、第１図装置における、二次元で収集したデータを処理する処理手段の詳細を表わすフローチャート；第１０図は、第１図の装置における、ＭＲＩ用に、被検体から収集したイメージを処理する処理手段の詳細を示すフローチャート；第１１図は、第１図の装置における、処理したデータを分析する分析手段を詳細に表わすブロックダイヤグラム；第１２図は、第１図の装置における、３ＴＰイメージを分析するＭＲＩに使用される分析手段を詳細に示すブロックダイヤグラム；第１３図は、３つの時刻ｔ０、ｔ１、ｔ２およびデータ値Ｉ（ｔ０）、Ｉ（ｔ１）、Ｉ(ｔ２)に対する、ウォッシュアウト／スローウォッシュインおよびスローウォッシュアウト／ファーストウォッシュインのパターンを表わすグラフ；第１４図は、第１３図と同様に、モデレートウォッシュアウト／スローウォッシュインおよびモデレートウォッシュアウト／ファーストウォッシュインのパターンを表わすグラフ；第１５図は、第１３図と同様に、ファーストウォッシュアウト／スローウォッシュインおよびファーストウォッシュアウト／ファーストウォッシュインのパターンを表わすグラフ；第１６図は、ＭＲＩ用の新規な装置に利用される、代表的なキャリブレーションマップ（較正図）を、表わす線図。好適実施例の詳細な説明第１図には、本発明を具現化したモニタ装置のブロックダイヤグラムが図示されており、このモニタ装置によって、システムを時間空間でモニタする。この装置には、選択手段１が設けられており、この選択手段１によって、データ収集のための時刻を設定し、データ収集用に利用されるこの装置の各部におけるパラメータを設定し、更に、処理したデータの分析用の着色（ｃｏｌｏｒｉｚｅｄ）済みキャリブレーションマップ（較正図）を処理する。これら処理データは、独特のウォッシュアウト／ウォッシュインビフェイビア（振舞）に対するこのシステムの変数に関連しており、これらによって色相および色強度を利用して、上述したビフェイビアをコーデング（符号化する）。この選択手段を記憶手段２に接続する。この記憶手段２には、選ばれた回数、パラメータ、およびキャリブレーションマップが記憶されている。この記憶手段２を、制御／実行手段３に接続する。この制御／実行手段３は、選択手段１によって決定されたデータ収集の回数およびパラメータを制御すると共に、このデータ収集を実行する。このデータ収集実行手段３は、既知なもので、この装置で利用される。また、制御／実行手段３を、記憶手段４に接続し、この記憶手段４には、時間空間で収集したデータが記憶されている。この記憶手段４を、処理手段５に接続し、この処理手段５によって、新規でユニークなウォッシュアウト／ウォッシュインスキームに従って、これら記憶／収集データを処理する。このスキームによって、色相および色強度が利用される。また、このスキームは、この処理手段５にプログラムされている。また、この処理手段５を、記憶手段２にも接続して、１で設定した、その入力パラメータおよび値の一部として受信する。この処理手段５を、これら処理データをストアする記憶手段６に接続する。この記憶手段６を、処理データを分析する分析手段７に接続する。この分析手段７を、また、記憶手段２に接続して、ストアされているキャリブレーションマップに基いて、処理データを分析する。これら全部の記憶手段、一部の記憶手段、または、全くこれら記憶手段を用いないで、これらすべての手段を、順次作動されることが可能となる。記憶手段の代わり、ブロック１〜３、ブロック３〜５、またはブロック５〜７間に、ダイレクト出力−入力リンクを設けることもできる。この装置の各部は、全く同時に機能する必要もないし、同一位置において機能する必要もない。また、第２図には、第１図のブロック１を示された、選択手段の詳細なブロックダイヤグラムが示されている。この選択手段の機能は、データ収集の時刻、およびデータ収集のパラメータを選択するためのものであり、且つ、新規なキャリブレーションマップ（較正図）を生成する機能である。このマップは、モニタすべき特定のシステムの変数に対して、色相／色強度でコーディングしたウォッシュアウト／ウォッシュイン特性に関連している。このプログラムの入力としては、１０、１１、１２および１３が存在し、これらには、以下のものが包含されている。即ち、モニタ中のシステムの変数に対する最小値および最大値（ブロック１０）；データ収集用に用いられる一定パラメータの入力（ブロック１１）；可能できるデータ収集用のパラメータ（ブロック１２）；および、３の時刻である。これら時刻は、繰返し手順の開始点として、このシステムの一般的な知識から直感的想像として、当初より選択されている（ブロック１３）。ブロック１０の入力によって、このシステムにおいて、これら変数によって、規定されたシステムスペースをブロック１４中のグリッドに分割することができる。１個の変数に対しては、各グリッドポイントは、単位長さを有する。２個の変数に対しては、各グリッドポイントによってピクセルを規定する。３個の変数に対しては、各グリッドポイントは、ボリューム（ボックセル）を規定する。ブロック１４で発生されたグリッドおよびブロック１１、１２、１３の入力によって、ブロック１５において、各グリッドポイントに対して、パラメータに対する値が決定されるように作用し、このパラメータは、予じめ選ばれた回数Ｉ（ｔ０）、Ｉ（ｔ１）、およびＩ（ｔ２）のそれぞれで測定される。この決定、即ち演算では、モニタ中の特定のシステム特有の方程式が利用され、この式によって、このシステムの変数の各組に対して、測定すべきパラメータの値における時間変化を、正確または実質的に概算することができる。このような方程式は、特定のシステムに対する分野では周知である。ブロック１５において演算されたＩ（ｔ０）とＩ（ｔ１）との値を利用して、ブロック１６において、強度関数を決定、即ち演算する。この強度関数は、ウォッシュインの初期レートのウォッシュイン特性を表わす。次に、ブロック１７において、従来のループによって、最大強度関数を有するグリッドポイントを決定すると共に、次に、これを最大強度としてブロック１８に出力する。ブロック１８において、残りのグリッドポイントにおける他の強度関数の全部を、この最大強度を参照し乍ら正規化する。次に、ブロック１９において、この装置によって、各グリッドポイントを、第３図に示したプログラムに基いて着色する。第３図において、ウォッシュアウト／色相のパターンを決定する関数Ｃを、各グリッドポイントに対して、Ｉ（ｔ１）およびＩ（ｔ２）に関して、ブロック２２と２３とにおいて、演算、即ち、トパターン／色相を選択する。第１３図、第１４図および第１５図は、色相でコーディングされたウォッシュアウトパターンを選択する原理および、色の強度でコーディングされた初期レート（見掛け上の初期レート）におけるウォッシュを選択する原理を、系統的に説明する。第１のパターン（第１３図）は、Ｉ（ｔ１）＜Ｉ（ｔ２）によって規定されたスローウォッシュアウトパターンであり、このパターンは、ンを、赤色でコーディング（符号化）する。（これは、ウォッシュインの見掛け上の初期レートの測定値である）が、スローウォッシュアウトパターンに対して高い値を有する場合には、この赤色は、ファーストウォッシュイン（第１３図）を反映して明るくなる。この強度関数が、スローウォッシュアウトパターンに対して、低い値を有する場合には、この赤色は、スローウォッシュイン（第１３図）を反映して、暗くなる。ォッシュアウトパターン（第１４図）であり、のパターンを緑色でコーディングする。上述の赤色の場合のように、この強度関数が、このモデレートパターンに対して、高い値を有する場合に、この緑色は、明るく（第１４図）なる。しかし乍ら、この強度関数が低い値を有する場合には、この緑色は暗くなる（第１４図）、また、第３パターン（第１５図）は、Ｉ（ｔ１）＞Ｉ（ｔ２）によって規定されるファーストウォッシュアウトパターば、このパターンを青色でコーディングする。再度、この強度関数が高い場合に、即ち、ウォッシュイン見掛けの初期レートが速い場合には、この青色は、明るくなる（第１５図）。しかし乍ら、この強度関数が低い場合には、この青色は暗くなる（第１５図）。その結果、色相によってコーディングされたウォッシュアウトの各パターンに対して、ウォッシュインレートを規定することができ、このレートは、色強度を規定する強度関数に関連するものである。異なったウォッシュインレート間の分離は、選定した色強度の範囲に依存するものである。再び、第２図に戻り、ブロック２０では、選択されたｔ０、ｔ１、ｔ２および選択されたシステムパラメータに対する最終出力において、この装置によって、色相および色強度が各グリッドポイントに割当てられる。この各グリッドポイントによって、このシステムの変数の値を規定する。これら着色されると共に、強度正規化されたグリッドポイントを、集合的に、“キャリブレーションマップ” と称する。後述するように、選択されたｔ０、ｔ１、ｔ２および選択されたシステムパラメータに対する合成キャリブレーションマップが満足できるものであれば、このプログラムは終了し、更に、この出力をブロック２（第１図）中のストレージに送信するか、またはブロック３（第１図）に直接送信する。このキャリブレーションマップが満足できない場合には、この装置は、例えば、ｔ１とｔ２とを増大させることによって、３つの時刻の適当な調整を行なうと共に、ブロック１３に入力する。また、システムのパラメータに対して新しい値をインクレメントするか、および／または測定に対して新しい値（予じめ選択して）をインクレメントして、更に、ブロック１２に入力することも可能である。また、これら時刻およびパラメータの両者を変化させることもできる。ブロック１０および１１の元の入力を利用し乍ら、このプログラムにおける各ステップを繰返えすと共に、ブロック１０〜２０から実行する。この巡回（第２図中のブロック２１）は、予じめ選択された回数およびパラメータの一組に対して、満足なキャリブレーションマップが得られるまで、繰返えすことができる。ブロック１６（第２図）の強度関数の定義およびブロック２２と２３（第３図）のウォッシュアウト／色相関数のパターンの定義において、Ｉ（ｔ１）＞Ｉ（ｔ０）であるものと仮定している。Ｉ（ｔ１）が負であるが、｜Ｉ（ｔ１）｜の絶対値がＩ（ｔ０）より大きい場合には、これら定義はこれら絶対値に対して支持される。ブロック２３（第３図）における着色関数に対する選択は、Ｃ＜−０．１、−０．１≦Ｃ＜０．１、およびＣ≧０．１の特定な例に限定されるものでなく、他の有理数（例えば、０．０５または０．２）を利用して、このシステムに依存して、Ｃを規定することが可能である。第４図には、データ収集用のパラメータの設定およびキャリブレーションマップの生成のフローチャート（ステップまたはブロック）が図示されている。このフローチャートは、モニタされたパラメータ用に書かれたもので、このパラメータは、このシステムの２つの変数（文字“Ｋ”および“ν”で割当てられている）の関数として、時間と共に変化する。Ｋとνとの二次元グリッドにおける各グリッドポイントに対して、Ｋの０．０１ユニットおよびνの０．０１ユニットの次元のピクセルが、ステップ、またはブロック３８、４０、５３、５５で規定されている。このプログラムは、ブロック３０から開始し、時刻ｔ０、ｔ１、ｔ２；システムおよび測定パラメータ；ならびにそれらの最小値および最大値間のＫの範囲およびνの範囲の入力が、ブロック３１で得られる。このプログラムは、ブロック３２において、ピクセル（Ｋｍｉｎ、νｍｉｎ）から開始して、モニタ中のシステムに相関している近似式または正確な式を利用し、ブロック３３でＩ（ｔ０）、Ｉ（ｔ１）、Ｉ（ｔ２）を演算する。この式は、既知なように、時間Ｉ（ｔ）でモニタしているパラメータが、どのように、値Ｋおよびν、ならびに他のシステムパラメータに依存するかを予測するものである。これら演算した、即ち決定したＩ（ｔ１）およびＩ（ｔ０）を用いて、ブロック３４で示したように、各ピクセル強度（Ｋ、ν）に対して演算する。これらＩ（ｔ１）およびＩ（ｔ０）は、ウォッシュイン初期レートを表わしている。次に、このプログラムでは、ピクセルが最大強度を有すかどうか（ブロック３５〜４０）を探索すると共に、これをすべてのピクセルに対して、ブロック３３に戻るループモードで進める。また、再度、ステップ３４からブロック４０まで進み、これら処理は、最大値Ｋおよび最大値νを有するピクセルに達するまで行なう。このループ中において、最大強度を有するピクセルが識別されると共に、すべてのピクセル（Ｋν ）に対して、強度を演算する。ここで、ピクセル（Ｋｍｉｎ、νｍｉｎ）ブロック４３から開始して、各ピクセルに対して、正規化した強度を演算し、ブロック４４で、この最大強度に関連して正規化する。最大強度を有するピクセルに、強度Ｎに対する最大値を割当てる。このＮとしては、システムの要求に応じて、１、２、３または、例えば、８、６４、２５６（コンピュータ数字）等の他の数を使用することが可能となる。次に、ピクセル（Ｋｍｉｎ、νｍｉｎ）からスタートして、各ピクセルに対して、ウォッシュアウトパターンを演算し、この演算は、ピクセル（Ｋｍａｘ、νｍａｘ）に到達するまで実行される。また、ブロック４５〜５４までに図示したように、色相で、各パターンをコーディングする。ここで、すべてのピクセルは、色相および正規化された色強度を有するようになる。このことによって、出力において、選択したｔ０、ｔ１、ｔ２および、システムと測定パラメータに対して、ブロック５６において、Ｋ、νのキャリブレーションマップが生成される。ここで、このキャリブレーションマップが満足できない場合には、例えば、過度に、１つの色相に傾斜する場合、新しい時刻、新しいシステム又は測定パラメータ値、また、これらすべてが一方向に調整されて、このキャリブレーションマップを補正すると共に、更に、このマップを満足できるように、色分布の基準からバランスさせる。ここで、これら新しい入力を利用して、満足できるキャリブレーションマップが得られるまで、このフローチャートにおける全ステップを実行する。このプログラムによって、選択された時刻およびシステムパラメータが設定される。また、何がこのような満足できるキャリブレーションマップかは、このシステムに依存すると共に、当業者にとっては明らかなものである。大部分のシステムにおいては、満足できるマップとは、ピクセルの１／３を赤色、１／３を緑色、および１／３を青色を有するものである。第５図には、トレーサ変調型ＭＲＩ、別称“コントラスト強調型ＭＲＩ用のパラメータを設定するためのフローチャートの特別な例が図示されている。ＭＲＩ信号は、モニター中のパラメータであり、このパラメータは、別名“造影剤”として知られているトレーサを投与した結果として、時間と共に変化している。ステップ６における入力パラメータには、３つの時刻ｔ０、ｔ１、ｔ２が包含されている。これら３つの時刻は、初期値として、ステップ８７において、経験によって得られ、これら時刻は、最適の３つの時間が得られるまで、巡回処理される。他の入力としては、トレーサ線量およびＭＲＩパラメータであり、これらによって、ＭＲＩ信号を、何如にして記録するかを規定する。これらトレーサ線量およびＭＲＩパラメータの両者を一定値とすることもでき、また、ステップ８７において、巡回処理によって最適化することもできる。また、この入力には、薬物動態学パラメータが包まれており、これらパラメータによって、血液中のトレーサの時間変化が規定される。また、この入力には、２つの変数Ｋ、およびν用の最大／最小値が含まれており、これら値によって、このシステムの病理形態学、即ち、被検体を規定するものである。また、変数Ｋによって、微小血管透過率を規定するものであり、この透過率により、トレーサを漏出させるための血管の容量を概算する。また、変数νによって、細胞外の体積の有理数を規定し、これによって、組織中の自由空間の量を概算する。このフローチャート内のステップ６２〜８５は、第４図のステップ３２〜５４に追従する。このフローチャートにおいて、最大強度を、ステップ７４に割当てて、２５６の値を有するようにする。ステップ８６の出力は、２つの変数Ｋとνとのキャリブレーションマップから構成される。これら変数は、時刻と他の入力との特定の組合せに対して、Ｋｍｉｎ、 νｍｉｎから、Ｋｍａｘ、νｍａｘまでの間の変数である。特定のＫとνとの値を有するこのマップにおける各ピクセルは、色相と色強度とを有する。満足可能なキャリブレーションマップは、色、または色と色強度との或る分布に達することによって規定される。例えば、満足可能なマップとしては、Ｋ−ν平面を分割するマップ、または、これら３色間のボリュームを、ほぼ均等の３つの領域に分割するマップとすることができる。換言すれば、このキャリブレーションマップにおけるピクセルの１／３が赤、１／３が緑、および１／３が青である。第１６図には、本発明による、代表的なキャリブレーションマップが図示されている。このマップは、変数、微小血管透過率Ｋおよび細胞外ボリュームνの有理数に対する、公知であるコントラスト強調式に基いて生成されたものである。このマップは、３個、またはそれ以上の選ばれた時刻に対する予期されたパターン（例えば、色および色強度）として構成されている。更に、具体的には、第１６図には、２次元図面中に、３次元表示したものが利用されている。第１６図は、グラディエントエコーおよび０．０８ｍｍｏｌ／ｋｇのトレーサ線量を駆使して、時刻ｔ０＝０、ｔ１＝４、ｔ２＝１２ｍｉｎ．で、人体に対する時刻を表わしている。これら３つの時刻は、例えば、繊維腺腫と悪性腫瘍とを区別するために選択されている。これら等曲線は、各パターン中の同一強度、例えば、ウォッシュインの同じ初期レートの領域を表わす。１つの次元は、関心値（例えば、０．００ｍｉｎ^-1から０．３ｍｉｎ^-1）の間の微小血管透過率Ｋである。第２の次元は、．１〜１間で変化する細胞外ボリュームνの有理数で、第３の次元は、正規化された強度であり、これは、実際には、各色の強度（例えば、０〜２５６強度内の何れかの値である）。このキャリブレーションマップによって、最適なプリセット時刻ｔ０、ｔ１、ｔ２ならびに例えば線量等の他のパラメータを決定するように作用する。また、他の異なるシステムに対しては、他の異なる最適な時刻を選ぶことができる。また、このキャリブレーションマップによって、この装置の処理手段の出力を解釈するように作用する。この装置としては、後述するように、ＭＲＩの特定例として、３ＴＰイメージがある。また、これら３つの時刻を、互いに近接するように選択した場合には、このキャリブレーションマップによって、スローウォッシュアウトビフェイビアのみが表わされる。即ち、上述した例では、赤色パターンが支配的である。他方、最後の時刻を、かなり遅れて選んだ場合には、このキャリブレーションマップは、ファーストウォッシュアウト、即ち、青によって支配される。ある特定システムに対する好適な３つの時刻は、多く表わされてる方法において、３色全部が分布するような値を選ぶ、即ち、キャリブレーションマップでは、この領域の１／３が、各色によって占有されるようになる。第５図に示したフローチャートに類似したフローチャートが、他の特定のシステムに対して存在する。最終出力において、ｔ０、ｔ１、ｔ２および他の入力の最適な組合せに対する変数の満足できるキャリブレーションマップが、常時得られる。時刻およびパラメータが、ブロック１（第１図）において設定されると共に、制御／実行手段３によって利用されることによって、データ収集を制御および実行する（第１図）。また、第６図には、データ収集の制御／実行のための詳細なブロックダイヤグラムが表われている。ブロック９０および９１の入力が、ストレージ２（第１図）から直接得られるか、また、ブロック１（第１図）から直接得られる。ブロック９２、９５、および９７における測定が、各システムにおける特定の手段によって実行される。このような手段は、既知なものである。トレーサの投与の回数、および測定の回数を制御することによって、入力のｔ０、ｔ１、ｔ２が、正確に制御されるようになる。このシステム中へのブロック９４でのトレーサの注入または投与は、あらゆる既知な手法によって実行できる。トレーサ投与のタイミングは、時刻ｔ０で開始するように固定されており、通常、このｔ０に近い時刻で終了することが望ましいが、時刻ｔ１以前にする必要は無い。ブロック９２、９５、９７で収集したデータを処理手段５（第１図）に、ブロック３（第１図）から直接、またはストレージ４（第１図）から転送する。第７図には、データ収集を制御／実行する制御／実行手段のブロックダイヤグラムの特別な例が示されており、この手段は、変形したＭＲＩ装置の一部分または、トレーサ変調型ＭＲＩ用のマシーンとして作用する。第７図のブロックダイヤグラムに示した制御／実行手段は、第１図のブロック１における選択手段によって設定された時刻ｔ０、ｔ１、ｔ２をブロック１２０の入力として受信する。これら時刻ｔ０、ｔ１、ｔ２は、第４図のフローチャートにおいて、この特別な例に対して、詳細されている。他の入力は、ＭＲＩパラメータに関連すると共に、トレーサの線量や、このトレーサの薬物動態的パラメータに関連しており、このトレーサを、被検体の血液中に注入する。これらＭＲＩパラメータおよびトレーサの線量を、第１図のブロック１の選択手段によって設定すると共に、これらＭＲＩパラメータおよびトレーサ線量について、第４図のフローチャート中でトレーサ変調型ＭＲＩについて詳細に説明している。次に、第７図のブロック１２１において、装置は、このシステム中の規定された領域、または規定されたボリュームのイメージ、即ち、人体を、既知手段によって記録する。記録用パラメータは、第１図における選択手段１によって設定されると共に、第７図のブロック１２に入力されるパラメータである。次に、このイメージをブロック１２２にストアすると共に、別名、造影剤と称されるトレーサを、既知の手法によって、ブロック１２３において、人体中に、時刻ｔ０で、投与する。この投与は、ｔ０において開始するようにタイミング調整され、この投与期間は、好適には（しかし、必然的でなく）、時刻ｔ０に近接すると共に、時間ｔ１における記録動作の前に終了する。トレーサ投与の後、この装置は、時刻ｔ１において、ブロック１２４において、同一領域またはボリュームのイメージを記録する。このイメージは、ブロック１２１における同一ＭＲＩパラメータを利用して、同一人体において、このブロック１２１において、記録されたイメージのようなものである。このイメージを、ブロック１２５においてストアする。次に、時刻ｔ２において、この装置は、ブロック１２１におけるような同一のＭＲＩパラメータを利用し乍ら、同一人体中の同一領域、または同一ボリュームの他のイメージを記録し、次に、このイメージを、ブロック１２７のストレージにストアする。また、第８図は、第１図のブロック５に示した処理手段の更に、詳細なブロックダイヤグラムを表わす。第８図のブロック１００で示した入力は、第１図のブロック１でセットすると共に、第１図のブロック３で用いた時刻ｔ０、ｔ１、ｔ２である。第８図のブロック１０１に示した他の入力には、第１図のブロック３で収集したデータが包含されている。これらデータを、第１図のブロック３から、または第１図のブロック４のストレージから直接、転送することができる。各空間ユニット（単位）に対して、第８図のブロック１０１におけるデータを、３つの異なった時間で得られた３つのデータ値によって表わす。まず、各空間ユニットに対する第１データ値を、Ｉ（ｔ０）と称し、このデータ値は、トレーサの投与の前に得られる。トレーサおよび／または第３を、時刻ｔ０に投入する。従って、この第１データ値を、時刻ｔ０以前であるが、このｔ０に可能な限り近接した時刻において測定する。各空間ユニットに対する第２データ値が、時刻ｔ１で得られ、これをＩ（ｔ１）と称する。各空間ユニットに対して、第３データ値が、時刻ｔ２で得られると共に、このデータ値をＩ（ｔ２）と称す。第８図のブロック１０１におけるもう１つの入力は、最大強度である。この最大強度は、第２図のブロック２０で示したような選択手段の出力の一部である。これらデータＩ（ｔ０）、Ｉ（ｔ１）および時刻ｔ０、ｔ１から、次に、正規化された強度が、各空間ユニットに対して、ブロック１０２で決定されている。この正規化された強度は、以下の方程式によって与えられる：ここで、Ｉ（ｔ０）、Ｉ（ｔ１）および最大強度が、ブロック１０１に入力されている。“Ｎ”は整数であり、例えば、１、２、３または８、６４、２５６（コンピータ数字）とすることができる。システムの同一モニタ動作上では、第８図のブロック１０２におけるＮは、第４図のブロック４４におけるＮと等しいものである。次に、ブロック１０４において、この装置によって、第３図に示すプログラムに従って、各空間ユニットを着色する。また、ブロック１０６の最終出力は、色相および色強度で割当てられた各空間値を有するようになる。色相は、ウォッシュアウトパターンを表わすと共に、色の強度は、ウォッシュインの初期レートを表わす。ブロック１０６における着色された出力を、第１図のストレージ６に送給するか、または、第１図のブロック７における分析手段へ直接送給する。第８図のブロック１０６における出力中の、各空間ユニットでの各色および色強度に対して、キャリブレーションマップ中の少なくとも１個のグリッドポイント、または少なくとも１つの場所における同一色および同一色強度が存在する。このキャリブレーションマップには、このシステムの変数の規定値が設けられている。従って、処理手段５によって得られた、各空間ユニットにおける、色／色強度を、このキャリブレーションマップで決定した変数の規定値に関連付ることができる。第１図に示した装置は、正規化された強度と一緒に機能する。この正規化は、すでに詳述したように、キャリブレーションマップを得るためのプロセスにおいて、第１図のブロック１における選択手段によって実行されるものである。この選択手段によって、最大強度が規定されると共に、この最大強度の値を利用して、この最大強度に対して測定したすべての強度を正規化する。この代りに、好適なものではないが、最大強度を参照し乍ら選択手段によって実行した正規化を、選択手段および処理手段によって独立して最大強度を選択することにより、この処理手段によって実行した正規化を分離することもできる。しかし乍ら、このことによって、第８図のブロック１０６に示すと共に、第９図のブロック１２６、および第１０図のブロック１６６に示したデータ処理の最終出力と、キャリブレーションマップとの間の相関性が弱められてしまう。第９図には、収集したデータを処理する装置のステップまたはブロックのフローチャートが示されている。このフローチャートは、２次元で記憶したデータを取扱っているので、各データポイントはピクセルである。第９図における第１ステップ１１０は開始である。このステップに続いて、ステップ１１１の入力となる。この入力は、第１図のブロック１における選択手段によってセットされた選ばれた時刻ｔ０、ｔ１、ｔ２から成っており、これら時刻は、第１図のブロック３における制御／実行手段によって利用されている。空間ユニット（単位）“ｎ ”および“ｍ”によって、グリッド中のピクセルを規定する。第１ピクセルに対して、ｎ＝１およびｍ＝１で、従って、ｎは１からｎまで進み、ｍは１からｍまで、１のステップで進む。もう１つの入力は、各ピクセルの収集したデータＩ（ｔ０）、Ｉ（ｔ１）およびＩ（ｔ２）である。この入力は、第１図のブロック３における制御／実行手段、または第１図のブロック４のストレージより直接、得ることができる。この入力には、また、最大強度が含まれており、この最大強度は、第２図のブロック２０に示した選択手段の出力または、第４図のフローチャートのステップ５６で得られる。次に、第９図のステップ１１２と１１３とにおいて、この装置によって、第１ピクセルｎ＝ｍ＝１が選択されると共に、ステップ１１４において、このピクセル中の正規化された強度を決定する。この強度の最大強度を参考にした正規化によって規定される。ここで前述したようにＮ＝整数≧１である。この正規化された強度は、ウォッシュインビィフェイバに対する初期レートの測定値である。ステップ１１４のＮは、システムの同じモニタリング動作中では、第４図におけるステップ４４のＮと等しいものである。次に、このピクセル（ｎ＝ｍ＝１）を、ステップ１１５〜１２０において、第３図のブロックダイヤグラムに従った色関数Ｃおよび第４図のステップ４５〜５０に従って着色する。このことによって、このピクセルのウォッシュアウトパターン／色相を決定する。可能なパターン（３つのパターン）のいずれも発生しない場合には、例えば、Ｉ（ｔ１）＝０の場合には、このピクセルはステップ１２１で黒に色付けされる。ステップ１２２、１２３、１２４、１２５では、次のピクセルが選ばれると共に、ステップ１２０まで続行するステップ１１４へのループによって、この次のピクセルに対して、正規化された強度およびパターン／色相を決定する。このループは、すべてのピクセル（ｎ×ｍ）に対して、繰返えされる。第９図のステップ１２６の出力において、各ピクセルは、色によってコーディングされたウォッシュアウトパターンおよび、色強度によってコーディングされたウォッシュイン初期レートを有する。この出力は、ステップ１２７に進んで、この第９図のフローチャートを終了する。第１０図には、トレーサ変調型ＭＲＩにおいて収集したデータを処理する装置のフローチャートの特定の例を示す。このフローチャートは、ステップ１５０でスタートする。次のステップ１５１では、選択した時刻ｔ０、ｔ１、ｔ２の入力が行われる。これら時刻は、第５図のフローチャートで説明した装置におけるトレーサ変調型ＭＲＩ用に示した、この装置の選択手段によって設定される。またこのステップは、第５図の出力ステップ８６の一部分である。これら３つの時刻は、また、第１０図のステップ１５１における人体の同一モニタリング動作におけるデータ収集において、制御／実行手段によって利用される。また、最大強度が、第４図の出力ステップ８６から得られると共に、このシステムの同じモニタリング動作、即ち、人体から得られる。また、この入力には、第７図のトーサ変調型ＭＲＩで示した制御／実行手段によって収集されたデータが包含されている。これらデータがイメージとして記録される。同一領域、またはボリュームの３つの記録されたイメージ、即ち、３組の記録イメージが存在している。第１のイメージ、またはイメージの１組が、時刻ｔ１で記録されると共に、第３のイメージ、または、１組のイメージが時刻ｔ２で記録される。このイメージにおける各ピクセルまたはボックセルは、ＭＲＩ信号強度を有し、この信号強度は、トレーサの投与後、時間と共に変化する。トレーサ投与以前であるが、この投与時間に近接した時間に記録した第１イメージにおいて、この強度を、Ｉ_image1(n,m)と称する。また、ｎ＝１およびｍ＝１であるピクセルを、強度Ｉ_image1(1,1)を有する第１ピクセルと称する。ｎ×ｍ個のピクセルが存在し、ここでは、ｎは、１からｎまで進み、ｍは１からｍまで進む。時刻ｔ１において記録された第２イメージにおける各ピクセル（ｎ、ｍ）における強度を、Ｉ_image2(n,m)と称する。同様に、時刻ｔ２において記録した第３イメージの各ピクセルにおける強度を、Ｉ_image3(n,m)と称する。第１０図におけるステップ１５２、１５３において、ｎ＝１およびｍ＝１から開始して、これらピクセルが割当てられる。また、第１０のステップ１５４において、正規化された色強度は、ピクセル（ｎ、ｍ）に対して、Ｉ_image1(n,m)およびＩ_image2(n,m) におけるピクセル（ｎ、ｍ）の値から、以下の式に従って決定される：すべてのピクセルに対して、すべてのピクセルが処理されるまで、ステップ１６２、１６３、１６４および１６５中のループによって、正規化された強度が決定される。次に、ステップ１５５〜１６０において、各ピクセル（ｎ、ｍ）に対して、ウォッシュアウトパターン／色が、値Ｉ_image2(n,m)およびＩ_image3(n,m)から、色関数Ｃを利用して決定される。この色関数Ｃは、この特定な例に対して、以下のように規定される：ならびに、第３図のブロック２３において示したように、および第１０図のブロック１５５〜１５９において示したように、色が選択される。ステップ１５５〜１６０までにおいて、着色されずに残っているピクセルは、ステップ１６１において、Ｉ_image2(n,m) ＝０が黒に割当てられる場合である。着色ステップは、ステップ１６２、１６３、１６４、１６５によってループに形成されて、すべてのピクセルを包含する。ステップ１６６における出力は、ｎ×ｍイメージを表わし、このイメージにおいて、各ピクセル（ｎ、ｍ）は、規定されたウォッシュアウトパターン／色相、および規定されたウォッシュインレート／色強度を有する。この着色されたイメージは、３時刻イメージ、短縮して、３ＴＰイメージと称する。色相／色強度は、キャリブレーションマップに相関すると共に、微小血管透過率Ｋおよび細胞外ボリュームνの有理数の値および、第５図のステップ８６の出力のキャリブレーションマップ中の被検体の２つの変数に照して、解釈される。また、第１１図には、第１図の分析手段７の構造を表わすブロックダイヤグラムである。ブロック２００における入力は、キャリブレーションマップから成り、このキャリブレーションマップは、選択手段１によって得られると共に、第２図のブロック２０中の出力または、第４図のブロック５６中の出力に含まれている。第１１図のブロック２０１の他の入力は、第９図のブロック１２６と、第８図のブロック１０６に示した処理手段の出力である。これらブロック２００中の入力と、ブロック２０１中の入力とは、対象となるシステムの同一のモニタリング動作から得られる。ブロック２０２の分析は、例えば、ヒストグラムを利用して、色の分布および色強度の分析から構成される。また、この分析は、ブロック２０１の入力の各空間ユニットに対するキャリブレーションマップと、色／色強度との間の相関性から成る。この相関性によって、ブロック２０１の入力における各空間ユニットに対するキャリブレーションマップの変数の値を概算する。最終的に、ブロック２０３の装置によって、この分布分析およびキャリブレーションマップ付きの相関性が出力される。第１２図には、トレーサ変調型ＭＲＩで得られた３ＴＰイメージの分析用の手段を包含した装置の一部分の特定例を示す。ブロック２１０の入力は、３ＴＰイメージである。この３ＴＰイメージのブロック２１２での分析は、色分布および色強度分布の分析から成り、この分析としては、例えば、何個のピクセルを、或る強度で赤に着色するかを決定すると共に、他の色に対して、同様の決定を行なうかである。また、強度分布および色分布の独立した分析を行なうことができる。ブロック２１１における装置の一部分によって、第５図のトレーサ変調型ＭＲＩに対して図示したように、選択手段によって得られるキャリブレーションマップを入力する手段が得られる。このキャリブレーションマップは、同じｔ０、ｔ１、ｔ２用のもので、および３ＴＰイメージとして、他の割当てられたパラメータ用のものである。キャリブレーションマップの２つの変数Ｋおよびνに関連して、３ＴＰイメージの分析が、第１２図のブロック２１３に表わされている。これにより、３ＴＰイメージの各ピクセル中の色相／色強度を、キャリブレーションマップ中の色相／色強度によって決定された、薬物動態学パラメータ値、微小血管透過率Ｋ、および細胞外ボリュームνの有理数に相関させる。従って、この分析が、撮影された領域またはボリューム中の２つの薬物動態学変数の分布に照らして実行される。最終的に、ブロック２１４において、この装置によって、ストアまたは出力中に、ディジタル信号をストアするか、または、モニタのようなディスプレイデバイス上で表示するか、または、プリンタに送給して、１つまたは、複数の３ＴＰイメージのカラープリントを行なうと共に、対応する相関的なキャリブレーションマップを得るようにする。他の特定な例としては、かんがい（灌漑）システム用のコントロール／モニタ装置がある。通常使われる最新の技術で、有効な手法で大規模領域を、かんがい（用水）する技術の１つとして、ドリップイリゲータ（ｄｒｉｐｉｒｒｉｇａｔｏｒ）によるものがある。そのようなシステムを構築する場合には、以下のような、ニーズによって調整可能なパラメータが存在する。即ち、パイプの寸法、滴下（ドリップ）孔のサイズおよび長さ、およびかんがいの圧力およびタイミングのパラメータ。これらパラメータの全体によって、単位面積当りの水の滴下速度が、決定され、本例では、文字Ｋによって割当てられている。しかし乍ら、このかんがいの効率を決定するもう１つのパラメータは、本例において、文字νで割当てられている、地中における見掛けの水分散定数である。この分散率、または定数は、地中の土の物理的および化学的特性に依存するもので、この地中を、水が通過する。換言すれば、例えば砂のような軽い土を保有する領域では、見掛けの分散定数は大きなものとなり、他方、重い土を保有する領域では、この定数は小さなものとなる。従って、この分散率は、かんがいすべきフィールドによって変化する。本発明の装置によれば、Ｋおよびνを見積ることができ、これによって、かんがい効率を最適化できる。予備条件として、土の単位重量当りの水量を測定する必要がある。水定数を決定する方法には、数種類の方法がある。一例としては、堀った、直後の土の量を正確に測る。土を完全に乾燥させると共に、次に、再度、その重量を測定する。この重量における損失は、このサンプルの水の重さに等価なものとなる。フィールドのサイズおよび到達すべき精度に依存して変化する空間解像度を有するサンプルを採用することができる。例えば、極めて大きなフィールド（数十〜数百kｍ² ）に対しては、このフィールドを１ｋｍ² 単位に分割することが当然のことである。従って、本発明の新規な装置（特に処理手段）によって構成されたかんがいの最終イメージにおける各ピクセルによって、１ｋｍ² 当りのビィフェイビアが反映される。土のサンプルは、この１ｋｍ² の単位面積中の中央、または他の規定された位置から採取する必要がある。このサンプルのサイズを変化させることは可能であるが、約１ｇより小さくすることが可能である。また、このサンプルを採取する深さは、ニーズに従って変化させられる。例えば、かんがいすべき植物の成長は、根のレベルの水の量に依存する場合には、このサンプルを、このレベルから採取する必要がある。また、この装置を、３Ｄで利用することができ、且つ、サンプルを同一エリアから異なった深さで採取できる。複数のサンプルは、ほぼ同一場所（このサンプルのサイズは、全体の単位面積よりかなり小さいものである）から採取すべきである。本発明の装置によれば、以下の動作が行われる。テストされるエリアには、テスト準備ができた滴下システムが設けられると共に、サンプリングの位置が割当てられている。このかんがいシステムの動作直前の時間において、これら割当てられたすべての位置からサンプルが採取される。次に、時刻ｔ０において、このかんがいシステムの動作が、プリセット時間だけ開始し、時刻ｔ１の前に終了する。この時刻は、キャリブレーションマップによって決定される。かんがい動作が停止した後に、時刻ｔ１において、各割当て位置から、測定用に、第２サンプルが採取される。最終的に、プリセットされた時刻ｔ２において、各割当て位置から第３のサンプルが、測定用に採取される。次に、水の量がこれらサンプル中で決定される。時刻ｔ０、ｔ１、およびｔ２において、各位置において、採用された３つのサンプル中の水分内容は、ウォッシュインレートおよびウォッシュアウトパターンに従って変化する。また、この水分内容は、本発明による新規な装置および方法によって開発されたような色相／色強度コードによって特徴付けられている。ウォッシュインおよびウォッシュアウトビィフェイビアは、単位面積当りの水滴下レート＝Ｋおよび地中の、見掛け上の水分散定数＝νに依存する。同じνに対しては、初期レートは、Ｋと一緒に増大し、同一Ｋに対しては、この初期レートは、νと一緒に増大する。このウォッシュアウトパターンは、Ｋおよびνにも依存する。例えばνがＫより大きい場合には、ファーストウォッシュアウトパターン（青に割当てられた）が支配的になる。また、Ｋおよびνが同じ程度の場合には、緑色のコーディングモデレートウォッシュアウトが支配的になる。σがＫより小さい場合には、スローウォッシュアウトパターン（赤色）が支配的になる。全位置（各位置は、ピクセルよって記述されている）上のＫ、νの分布は、このキャリブレーションマップを参照し乍ら、これらピクセルの着色および色強度（各々は、色相および色強度を有している）を相関させることによって決定される。このキャリブレーションマップによって、時刻ｔ０、ｔ１、ｔ２に対して、Ｋ −ν平面上で、ウォッシュイン強度関数およびウォッシュアウトパターンが、以下のように与えられる。即ち、時刻ｔ１とｔ２とにおいて、累積された水分量の各Ｋ、ν対に対する演算が、既知のモデルに基いて実行される。このＫの範囲が、対応のかんがいシステムの最高レベルまで０から選ばれる一方、νの範囲が、０から、純水の分散定数、又は、既知ならば、かんがい処理すべきフィールドにおける最大の水分散定数までから選ばれる。Ｋとνとの範囲以内において、最良の分解能を得るために、時刻ｔ０、ｔ１、ｔ２を以下のように選ぶ、即ち、Ｋ−ν 平面を、３つのウォッシュアウトパターン／色間で、３つのほぼ同一面積に分解する。これら測定が一旦行われると共に、本発明による新規な装置が、その初期機能の下で実行されると、この装置を調整して、Ｋ（このＫは、かんがいシステムに依存するので、可変である）を以下の方法によって変化させることが可能となる。即ち、ピクセルによって表わされた各規定領域におけるかんがいが最高効率で行われるように変化させる。例えば、モデレートウォッシュアウトパターン（緑色）は、所定期間において、一定の水分量を得るために好適なものである。Ｋを増大させることによって、同じνに対して、青領域から緑領域まで移動できる。最終調整を達成できるようにするため、Ｋを変化させ、次に、対象フィールド中の位置および同じ３つの時刻で、この装置を作動させることができる。かんがいの色／色強度コーディングされたマップの他の利用としては、植物の種類を調整することによる植物プランの準備が存在する。即ち、土の種類およびかんがいシステムによって規定されるかんがいの数量に対する植物の密度を調整することである。また、本発明の新規な装置を利用することによって、空調プランをテストすると共に変更するために利用でき、例えば、加熱または冷却、更には、両方の機能に利用できる。以下に、加熱プロセスを制御することおよび／またはモニタリングするための新規な装置について詳述する。複数の室または規定された複数空間において、家屋、工場、オフィスビル、ショッピング街、または集合住宅を有する全体構造に対して建築した空調システム（暖房、冷房、またはそれらの両方）に対しては、各規定のエリアを、効率的に空調すると共に、その結果として、実物そっくりにすることを確認できるようにこの空調システムを設計することが重要なことである。或るエリアを過剰に加熱する一方、他のエリアを過剰に冷却することが可能となる。そのような構造体用の新規な装置による調整は、単位面積当りの熱量および本例では、文字Ｋで割当てられた単位時間当りの熱量を変更することによって実行できる。熱容量を決定する他の変数としては、周囲環境への単位面積当りの熱伝導率がある。これは、本例では、文字νで割当てられた隔離における不完全性によるものである。この加熱に対する評価は、各室、または規定空間内のあらゆる位置に配置された熱電対または温度計を使用して温度を測定することによって行われる。各温度計位置によって、対象システム（構造物）の２Ｄまたは３Ｄプランにおけるピクセル位置が規定される。このシステムがオフされている、ｔ０の前の時刻において、すべて位置におけるすべての温度測定デバイスの読取りが行われ、温度Ｔ０が決定される。時刻ｔ０において、空調システムが、プリセット時間だけオンとなる。このプリセット時間は、時刻ｔ１の前で、近いものである。この時刻ｔ１において、システムがオフされてしまった後で、すべてのデバイスの２回目の読取りが行なわれて、温度Ｔ１が決定される。次に、時刻ｔ２で終了する第２プリセット時間の後に、第３回目の読取りが行なわれて、各位置の温度Ｔ２が決定される。時刻ｔ０、ｔ１、ｔ２間での温度変化、即ち、各測定位置における温度Ｔ０、Ｔ１、Ｔ２間の差は、室に放出された、時間（Ｋ）／面積当りの空調量、および拡散（ν）によって壁部を介して外部への放出される熱損失または熱吸収量に依存するものである。後者のパラメータνは、測定が行なわれる位置と、環境との間の熱流方向に依存して、負または、正の値とすることができる。これらの変化は、空気流のウォッシュインパターンとウォッシュアウトパターンによって記述できる。このウォッシュアウトのパターンは、色相で記述される。時刻ｔ１における温度Ｔ１が、時刻ｔ２における温度Ｔ２より高い場合には、ウォッシュアウトプロセスは、ファースト（ｆａｓｔ）に規定されて、青が割当てられる。Ｔ１がＴ２より低い場合には、このウォッシュアウトプロセスは、スローとして規定されて、赤が割当てられる。この温度が、同一時刻Ｔ１＝Ｔ２（装置によって決められた範囲内において）で維持される場合には、ウォッシュアウトプロセスは、モデレート（ｍｏｄｅｒａｔｅ）として規定され、緑色が割当てられる。これら色の輝度は、れる初期レートに依存するようになる。また、この初期レートは、Ｋとνとに依存する。同じνに対して、この初期レートは、Ｋと共に増大する。また、同じＫに対して、この初期レートは、νを増加させ乍ら、減少すると共に、νを負の値に減少させ乍ら、増大する。また、その場所が十分に隔離されていない共に、その周囲がより冷たい場合には、νは大きくなると共に、温度変化によって、ファーストウォッシュアウトのパターンが得られる。また、この隔離が良好でない場合（νが小さく、０に近似する場合）、温度変化によって、モデレートウォッシュアウトのパターンが得られる。周囲が暖かく、室が十分に隔離されていない場合には、νは、大きな負の値を有し、温度変化によってスローウォッシュアウトパターンが得られる。本発明による新規な色相／色強度のコンセプトを利用して、十分に空調、例えば暖房されていない場所を識別できると共に、十分に、隔離されていない場所、および過剰空調されている、例えば、過剰暖房されている場所を識別できる。空調が暖房の場合には、この温度は、熱量によって決定される。この熱は、温度測定器、例えば温度計が設置された場所に達するものであり、および周囲からの熱損失または熱利得の結果として、この場所に残存する熱量によって決定される。室のサイズおよび熱源の分布によって規定された所定ケースにおいて、加熱時間中、熱流は、比較的速くなると共に、この室での平衡が急速に達成されるものと仮定できる。従って、この温度は、この加熱期間中に発注された総熱量に依存するものである。従って、Ｋは、０とこの加熱システムの最大容量との間で存在している。また、周囲からまた、周囲への熱流は、同じ絶縁、即ち、同じνであるが、逆の符号によって決定される。従って、νは、−νから＋νまでの間で、例えば、無絶縁での値によって決められる実際の値と共に、存在する。時刻ｔ０、ｔ１、ｔ２は、既知の方程式に基いて形成されたキャリブレーションマップを利用して選択され、この方程式は、前述したように、Ｋまで、時間と共に、変化する温度に関連している。特定な例として、加熱について説明しているが、この新規な装置は、同じプログラムと、冷却を評価する手段と共に利用でき、この場合には、ウォッシュアウトパターンおよびウォッシュインパターンを決定する温度変化に対する絶対値を利用している。本発明の装置には、電子的、光学的、または両方で作動するコンピュータシステムが設けられており、このコンピュータシステムには、メモリ、ＣＰＵ、ディスプレイ、デバイスイベント信号を発生する入力デバイス、ならびに、これらと協動するソフトウエアが設けられている。このソフトウエア（二進数、または関連したフォーム）は、コンピュータ解読可能なプログラムコードを有するコンピュータ利用可能な媒体を有する。このプログラムコードには、種々のフローチャートを実行するプログラムロジックおよび上述したブロックダイヤグラムを実行するロジックが包含されている。これらコンピュータの詳細については、周知であると共に、当業者であれば、これらフローチャート／ブロックダイヤグラムを実行できる十分な経験および知識を有しているので、上述した本発明を十分且つ完全に理解／認識するのに不必要且つ過剰なので、特別なハードウェアの詳述は省略するものとする。また、当業者であれば、本発明の方法および装置を、記載の詳細な説明／技術思想より実現できる。以上、本発明について詳述したが、種々の変形および変更を、本発明の技術思想から逸脱しないで実行できることは、当業者にとって明白であり、これら変形については、添付の請求範囲から規定されるものである。

───────────────────────────────────────────────────── フロントページの続き (81)指定国ＥＰ(ＡＴ，ＢＥ，ＣＨ，ＤＥ，ＤＫ，ＥＳ，ＦＩ，ＦＲ，ＧＢ，ＧＲ，ＩＥ，ＩＴ，ＬＵ，ＭＣ，ＮＬ，ＰＴ，ＳＥ)，ＯＡ(ＢＦ，ＢＪ，ＣＦ，ＣＧ，ＣＩ，ＣＭ，ＧＡ，ＧＮ，ＭＬ，ＭＲ，ＮＥ，ＳＮ，ＴＤ，ＴＧ)，ＡＰ(ＫＥ，ＬＳ，ＭＷ，ＳＤ，ＳＺ，ＵＧ)，ＵＡ(ＡＭ，ＡＺ，ＢＹ，ＫＧ，ＫＺ，ＭＤ，ＲＵ，ＴＪ，ＴＭ)，ＡＬ，ＡＭ，ＡＴ，ＡＵ，ＡＺ，ＢＡ，ＢＢ，ＢＧ，ＢＲ，ＢＹ，ＣＡ，ＣＨ，ＣＮ，ＣＵ，ＣＺ，ＤＥ，ＤＫ，ＥＥ，ＥＳ，ＦＩ，ＧＢ，ＧＥ，ＨＵ，ＩＬ，ＩＳ，ＪＰ，ＫＥ，ＫＧ，ＫＰ，ＫＲ，ＫＺ，ＬＣ，ＬＫ，ＬＲ，ＬＳ，ＬＴ，ＬＵ，ＬＶ，ＭＤ，ＭＧ，ＭＫ，ＭＮ，ＭＷ，ＭＸ，ＮＯ，ＮＺ，ＰＬ，ＰＴ，ＲＯ，ＲＵ，ＳＤ，ＳＥ，ＳＧ，ＳＩ，ＳＫ，ＴＪ，ＴＭ，ＴＲ，ＴＴ，ＵＡ，ＵＧ，ＵＳ，ＵＺ，ＶＮ【要約の続き】ータ処理システムを設けて、このキャリブレーションマップを生成する。このキャリブレーションマップおよび予じめ選ばれた場所のイメージは、また、新規な手段である。

Claims

【特許請求の範囲】 1. 流体がその中を流れるシステムをモニタするに当り、この流体が、システム中の時間空間における変化によって特徴付けられるモニタ装置において；（ａ）前記システム内の予じめ選択された場所をモニタして、システムパラメータを表わすデータを収集するモニタ手段と、このシステムパラメータは、システムウォッシュインおよびウォッシュアウトビィフェイビアに関連した、少なくとも１つの変数の関数として時間と共に変化し；（ｂ）前記モニタ手段に応答して、前記収集したデータを受信すると共に、このデータを処理する第１データ処理手段と；（ｃ）前記モニタ手段を制御して、データ収集のパラメータを設定すると共に、前記システムパラメータと、システムイベントに相関した複数の時刻とを包含するシステムパラメータの入力を受入れる選択手段と；（ｄ）第２のデータ処理手段と；この第２データ処理手段には：（１）各々が最小から最大まで変化する、少なくとも２つの変数によって規定されたスペースをグリッドに分割するグリッド手段と；（２）各グリッド位置に対して、各時刻でモニタされている前記パラメータのウォッシュインを表わす演算した値／強度関数を決定する第１決定手段と；（３）最大強度のグリッド位置を決定すると共に、前記最大値を参照し乍らすべてのグリッド位置を正規化する第１正規化手段と；（４）ウォッシュアウトビィフェイビアを表わす色関数に基いて、すべてのグリッド位置を複数色の一色で着色する第１着色手段と；（５）すべての時刻において、すべてのグリッド位置を合成物にアレンジして、色相／色強度に相関した各グリッド位置を有する前記グリッドのキャリブレーションマップを展開するアレンジ手段と；（６）このキャリブレーションマップを分析して、このキャリブレーションマップの属性が、予じめ選択された基準に合致するかどうか決定し、不一致の場合、複数の時刻、またはシステムパラメータを実行すると共に、結果としてのキャリブレーションマップが上記基準を満すまで巡回する分析手段とが設けられており；（ｅ）前記第２データ処理手段において決定されたような前記複数の時刻および前記強度関数の最大値を、更に入力として受入れる前記第１データ処理手段とを具備し；この第１データ処理手段には：（１）前記システム中の予じめ選択された場所を、空間単位に空間的に分割する分割手段と；（２）各空間ユニットに対して、各時刻で収集した前記データに基いて、ウォッシュインビフェイバーを表わす値／強度関数を決定する第２決定手段と；（３）前記第１正規手段から得られた最大強度に基いて、各空間単位の強度を正規化する第２正規化手段と；（４）前記第１着色手段のウォッシュアウトビィフェイビアを表わす色関数に基いて、すべての空間単位を複数色の１色で着色する第２着色手段と；（５）すべての時刻に対して、すべての空間単位を、色相／色強度に相関した各空間単位と共に、予じめ選択した空間の合成物をアレンジするアレンジ手段とが設けられたことを特徴とするモニタ装置。 2. 前記第１着色手段によって、前記色関数が、スローウォッシュアウトビィフェイビアを表わすＸ以上の場合に、１色が得られ；このＸより小さく、且つ、モデレートウォッシュアウトビィフェイビアの表わすＹ以上の場合に、第２色が得られ；ファーストウォッシュアウトビィフェイビアを表わすＹより小さい場合には、第３色が得られ、ここで、前記ＸおよびＹは、Ｘが正であり、且つＹが負であると共に、１より少ない小数であることを特徴とする請求項１記載の装置。 3. 前記Ｘが０．１に等しい前記Ｙが０．１に等しいことを特徴とする請求項２記載の装置。 4. ３つの時刻として、ｔ０、ｔ１、ｔ３を利用し、このｔ０は前記システムイベントの開始の時刻であり；ｔ１は、このｔ０の後の時刻であり、且つ、ｔ２はｔ１以降の時刻であることを特徴とする請求項２記載の装置。 5. 前記システムイベントは、トレーサ材料を、このシステム中に流れる前記流体中に注入することであることを特徴とする請求項１記載の装置。 6. 前記装置がＭＲＩシステムであることを特徴とする請求項５記載の装置。 7. ディスプレイ手段を設けて、前記最終キャリブレーションマップおよび、前記予じめ選択された場所の合成物を表示したことを特徴とする請求項１記載の装置。 8. ストレージ手段を設けて、前記最終キャリブレーションマップおよび前記予じめ選択した場所の合成物をストアしたことを特徴とする請求項１記載の装置。 9. プリンタ手段を設けて、前記最終キャリブレーションマップおよび前記予じめ選択した場所の合成を印刷したことを特徴とする請求項１記載の装置。 10. 流体がその中を流れるシステムをモニタするに当り、この流体がシステム中の時間空間における変化によって特徴付けられるモニタ装置において；（ａ）第１出力を有し、前記システム内の予じめ選択された場所をモニタして、システムパラメータを表わすデータを収集するモニタ手段と；このシステムパラメータは、システムウォッシュインおよびウォッシュアウトビィフェイビアに関連した、少なくとも１つの変数の関数として時間と共に変化し；（ｂ）第２出力と第１入力とを有し、前記第１出力から収集したデータを受信すると共に、これを処理する第１処理手段と；（ｃ）第３出力を有し、前記モニタ手段を制御し、データ収集のパラメータを設定し、および入力を前記第１入力に与える選択手段と；この選択手段には：（ｉ）前記システムパラメータと、システムイベントに相関した、複数の予じめ選択した時刻とを有するシステムパラメータの入力を受信する入力手段と；（ii）第２処理手段とが設けられ；この第２処理手段には：（１）各々が最小から最大まで変化する、少なくとも２つの変数によって規定されたスペースをグリッドに分割するグリッド手段と；（２）各グリッド位置に対して、各時刻でモニタされるパラメータの強度の演算された値を決定する第１決定手段と；（３）前記グリッド位置の最大強度を決定すると共に、この最大強度を参照し乍ら、すべてのグリッド位置を正規化する第１正規化手段と；（４）ウォッシュアウトビィフェイビアを表わす色関数に基いて、すべてのグリッド位置を複数色の一色で着色する第１着色手段と；（５）すべての時刻において、すべてのグリッド位置を合成物にアレンジして、色相／色強度に相関した各グリッド位置を有する前記グリッドのキャリブレーションマップを展開するアレンジ手段と；（６）このキャリブレーションマップを分析して、このキャリブレーションマップの色分布が予じめ決められた条件を満たすかどうかを決定し；および満たさない場合には、前記複数の時刻の１つおよびシステムパラメータを実行し；前記最終キャリブレーションマップが前記予じめ決められた条件を満たすまで巡回のための入力手段に入力する分析手段とが設けられ；（ｄ）前記第２処理手段中で決められたような前記複数の時刻と前記強度関数の最大値を、別の入力として受信する前記第１処理手段とを具備し；この第２処理手段には：（ｉ）前記システム中の前記予じめ選択された場所を、空間単位に空間的に分割する分割手段と；（ii）各時刻において収集したデータに基いて、各空間単位について、値／強度を決定する第２決定手段と； (iii)前記第１正規化手段から得られた最大強度に基いて各空間単位の強度を正規化する第２正規化手段と；（iv）前記第１着色手段で利用したウォッシュアウトビィフェイビアの相関性に基いて、すべての空間単位を、複数色の一色で着色する第２着色手段と；（ｖ）すべての時刻で、すべての空間単位をアレンジして、色相／色強度に相関した各空間単位と一緒に、前記予じめ選択した場所の合成物を構成するアレンジ手段とを有したことを特徴とするモニタ装置。 11. メモリと、デバイスイベント信号を発生する入力デバイスと、ディスプレイとを有するコンピュータと共に利用するソフトウエアにおいて、このソフトウエアには、コンピュータによって読取り可能なプログラムコードを有するコンピュータ読取り可能な媒体が設けられ；このプログラムコードには：システムの予じめ選択されたパラメータに相関した２つの変数によって規定されたスペースを、グリッドに分割する第１プログラムロジックと；このシステム中を流体が通過すると共に、このシステムのパラメータは、このシステムにおける時間空間での変化によって特徴付けられており、各変数は、最小値から最大値まで変化し；前記デバイスイベント信号に応答して、予じめ決められた強度関数に基いて、複数の時刻の各々の時刻における各グリッドの強度を決定する第２プログラムロジックと；この第２プログラムロジックに応答して、予じめ決められた色関数に基いてすべてのグリッド位置を着色する第３プログラムロジックと；およびこの第３プログラムロジックに応答して、色相／色強度に相関する各グリッド位置と共に、このグリッドの合成キャリブレーションマップを生成する第４プログラムロジックとが設けられたことを特徴とするソフトウエア。 12. 更に、前記第４プログラムロジックに応答して、前記キャリブレーションマップを、前記ディスプレイ上に表示する第５プログラムロジックを設けたことを特徴とする請求項１１記載のソフトウエア。 13. 更に、前記第４プログラムロジックに応答して、前記キャリブレーションマップを分析して、このマップが予じめ決められた条件を満たすかどうかを決定する第５プログラムロジックを設けたことを特徴とする請求項１１記載のソフトウエア。 14. 前記第５プログラムロジックに応答して、前記条件を満たすかどうかを決定すると共に、満たさない場合には、前記プログラムロジックを巡回させることによって、この条件を満たすまで前記キャリブレーションマップを調整する第６プログラムロジックを設けたことを特徴とする請求項１３記載のソフトウエア。 15. モニタされているシステムと一緒に利用されるキャリブレーションマップを生成するに当り、このシステム中に流体が流れ、この流体がこのモニタ中のシステムにおける時間空間での変化によって特徴付けられ、このモニタリング動作が予じめ決められた位置において起るようなデータ処理システムにおいて、（ａ）データを処理するコンピュータ処理手段と；（ｂ）データを記憶媒体にストアするストレージ手段と；（ｃ）この記憶媒体を初期化する手段と；（ｄ）各々が最小値から最大値まで変化する２つの変数によって規定される空間をグリッドに分割する処理手段と；これら変数は、前記モニタされたシステムの予じめ選択されたパラメータに相関しており；（ｅ）各グリッド位置に対して、複数の時刻の各々において、演算した値および強度関数を決定する処理手段と；（ｆ）前記グリッド位置の強度関数の最大値を決定すると共に、この最大値を参照し乍らすべてのグリッド位置を正規化する処理手段と；（ｇ）モニタ中のシステムビィフェイビアに相関した色関数に基いて、すべてのグリッド位置を、複数色の１色で着色する処理手段と；（ｈ）すべての時刻について、すべてのグリッド位置を合成物にアレンジして、色相／色強度に相関した各グリッド位置を有するグリッドのキャリブレーションマップを生成する処理手段と；および（ｉ）このキャリブレーションマップを分析して、このキャリブレーションマップの属性が予じめ決められた条件を満たすかどうかを決定すると共に、満たさない場合には、複数の時刻、または変数を実行すると共に、最終的なキャリブレーションマップの属性がこの予じめ決められた条件を満たすまで巡回する処理手段とを具備したことを特徴とするデータ処理システム。 16. 前記属性は色分布であることを特徴とする請求項１５記載のデータ処理システム。 17. 前記条件は、前記キャリブレーションマップの色の着色に実質的に同一であることを特徴とする請求項１６記載のデータ処理システム。 18. モニタ中のシステムと一緒に用いられるキャリブレーションマップを作成するに当り、このシステム中には流体が流れており、これがこのシステムにおける時間空間における変化によって特徴付けられており、予じめ決められた位置でモニタリング動作が行われるデータ処理方法において、以下のステップから成り；（ａ）データを処理し；（ｂ）データを記憶媒体にストアし；（ｃ）この記憶媒体を初期化し；（ｄ）各々が最小値から最大値まで変化する２つの変数によって規定されるスペースをグリッドに分割処理し、これら変数が前記モニタ中のシステムの予じめ選択されたパラメータに相関しており；（ｅ）各グリッド位置において、複数の時刻の各時刻におけるこのパラメータの演算した値／強度関数を決定するように処理し；（ｆ）これらグリッド位置の強度関数の最大値を決定すると共に、この最大値を参照し乍ら、すべてのグリッド位置を正規化するように処理し；（ｇ）このモニタ中のシステムビィフェイビアに相関した色関数に基いて、すべてのグリッド位置を複数色の一色で着色処理し；（ｈ）すべての時刻に対して、すべてのグリッド位置を合成物にアレンジすることによって、色相／色強度に相関した各グリッド位置を有するグリッドのキャリブレーションマップを生成処理し；および（ｉ）このキャリブレーションマップを分析して、このマップの属性が予じめ決められた条件に合致するかどうか決定処理し、合致しなければ、複数の時刻または変数を実行し、最終的なキャリブレーションマップの属性が上記条件に合致するまで巡回するように処理するステップを具備したことを特徴とするデータ処理方法。 19. 流体が通過するシステムをモニタするに当り、これがシステム中の時間空間における変化として特徴付けられるモニタ方法において、以下のステップより成り；（ａ）システムパラメータを表わすデータを収集するために、このシステムにおける予じめ選択された位置をモニタし、このシステムパラメータは、ウォッシュイン／ウォッシュアウトビィフェイビアに関連した少なくとも２つの変数の関数として時間と共に変化し；（ｂ）収集したデータを処理し；（ｃ）データ収集のパラメータを設定すると共に、システムイベントに相関した複数の時刻を選択することによって、上記モニタリング動作を制御し；（ｄ）各々が最小値から最大値まで変化する、少なくとも２つの変数によって規定される空間をグリッドに分割し；（ｅ）各グリッド位置において、各時刻において、モニタされるパラメータの演算された値／強度を決定し；（ｆ）これらグリッド位置の最大強度を決定すると共に、この最大強度を参照し乍ら、すべてのグリッド位置を正規化し；（ｇ）ウォッシュアウトビィフェイビアに相関した色関数に基いて、すべてのグリッド位置を、複数色の一色によって着色し；（ｈ）すべての時刻について、すべてのグリッド位置を、合成物にアレンジして、色相／色強度に相関した各グリッド位置を有するグリッドのキャリブレーションマップを生成し；（ｉ）このキャリブレーションマップを分析して、このマップの属性が予じめ決められた条件に合致するかどうかを決定すると共に、合致しない場合には、複数の時刻またはシステムパラメータを実行し、最終的なキャリブレーションマップの属性が予じめ選択された条件に合致するまで巡回するようにし；（ｊ）前記ステップ（ｆ）で決定された強度関数の最大値および複数の時刻を、前記ステップ（ｂ）の処理ステップに入力し；このステップ（ｂ）には、更に：（ｉ）このシステムの予じめ選択された場所を空間ユニットに、空間的に分割するステップと；（ii）各空間ユニットに対して、各時刻について収集したデータに基いて、値／強度関数を決定するステップと； (iii)前記ステップ（ｆ）から得られた最大強度に基いて、各空間ユニットの強度を正規化するステップと；（iv）前記ステップ（ｇ）の色関数に基いて、すべての空間ユニットを、複数色の一色で着色するステップと；（ｖ）すべての時刻について、すべての空間ユニットを、上記予じめ選択された場所の合成物にアレンジするステップが設けられ、この場所には、色相／色強度に相関した各空間ユニットが設定されていることを特徴とするモニタ方法。 20. 前記属性を、色分布としたことを特徴とする請求項１９記載の方法。 21. 前記条件は、実質的に着色であることを特徴とする請求項２０記載の方法。 22. 前記システムイベントは、前記予じめ選択された位置のシステム上流における前記流体中に、トレーサを注入することであることを特徴とする請求項１９記載の方法。 23. ２つの時刻を用いるようにしたことを特徴とする請求項１９記載の方法。 24. ３つの時刻ｔ０、ｔ１、ｔ２を用い、ｔ０は、注入時間に近い時刻であり、ｔ１は、システムウォッシュインビィフェイビアを表わすｔ０より予じめ決められた時間だけ遅い時刻であり、ｔ２は、システムウォッシュアウトビィフェイビアを表わすｔ１より、予じめ決められた時間だけ遅い時刻であることを特徴とする請求項２２記載の方法。２５．流体が通過するシステムをモニタする時に利用され、この流体が、システムパラメータにおける、時間空間での変化によって特徴付けられ；システムイベント以後の２つの時間間隔におけるシステムウォッシュインおよびシステムウォッシュアウトビィフェイビアに関連した２つの変数の関数として変化するようなキャリブレーションマップにおいて、このマップによって、二次元または三次元で、最小値から最大値まで変化する、これら２つの変数のイメージを描画し；このイメージの個別エレメントは、システムウォッシュアウトビィフェイビアを表わす複数色の１色の色相および、システムウォッシュインビィフェイビアを表わす色強度を有することを特徴とするキャリブレーションマップ。 26. 前記イメージは、モニタ上に表示されることを特徴とする請求項２５記載のキャリブレーションマップ。 27. 前記イメージを、コンピュータ読取り可能な媒体上にディジタル的に符号化したことを特徴とする請求項２５記載のキャリブレーションマップ。 28. 前記イメージを、印刷可能媒体上に印刷したことを特徴とする請求項２５記載のキャリブレーションマップ。 29. ３色を用いたことを特徴とする請求項２５記載のキャリブレーションマップ。 30. 前記色は、赤、青、緑であることを特徴とする請求項２９記載のキャリブレーションマップ。 31. 流体が通過するシステムにおける位置を、二次元または三次元で描画するに当り、これが、システムパラメータにおける変化によって、この位置において特徴付けられるイメージにおいて、システムイベントの後の第１時間間隔におけるシステムウォッシュインビィフェイビアおよび、第２時間間隔におけるシステムウォッシュアウトビィフェイビアに関連した２つの変数の関数として変化し；このイメージの個別エレメントは、システムウォッシュアウトビィフェイビアを表わす複数色の１色の色相およびシステムウォッシュインビィフェイビアを表わす色強度を有したことを特徴とするイメージ。 32. 前記イメージを、モニタ上に表示したことを特徴とする請求項３１記載のイメージ。 33. 前記イメージを、コンピュータ読取り可能媒体上に、ディジタル的に符号化したことを特徴とする請求項３１記載のイメージ。 34. 前記イメージを、印刷可能な媒体上に印刷したことを特徴とする請求項３１記載のイメージ。 35. 少なくとも３色を利用したことを特徴とする請求項３１記載のイメージ。 36. 前記色は、赤、青、緑であることを特徴とする請求項３５記載のイメージ。