JP2014000291A

JP2014000291A - 超音波診断装置

Info

Publication number: JP2014000291A
Application number: JP2012138554A
Authority: JP
Inventors: Eiji Kasahara; 英司笠原
Original assignee: Hitachi Aloka Medical Ltd
Current assignee: Hitachi Ltd
Priority date: 2012-06-20
Filing date: 2012-06-20
Publication date: 2014-01-09
Anticipated expiration: 2032-06-20
Also published as: JP6063154B2

Abstract

【課題】再構成処理により得られる画像の表示に係る改良技術を提供する。
【解決手段】再構成処理部２０は、前メモリ１４に記憶された複数周期に亘る複数の断層画像データの中から、互いに周期的に対応した複数の断層画像データを抽出し、それら互いに周期的に対応した複数の断層画像データから各再構成画像データを形成して、後メモリ２６に記憶する。マスク処理部３０は、各再構成画像データ内において対象物を含んだ注目領域以外の領域をマスク処理することにより注目領域を決定する。表示処理部４０は複数の時相に亘って複数の再構成画像データから次々に得られる注目領域の画像を一定の位置に配置した表示画像を形成し、その表示画像が表示部４２に表示される。
【選択図】図１

Description

本発明は、周期的な運動をする対象物の画像を表示する超音波診断装置に関する。

心臓などの周期的な運動をする対象物の三次元画像を形成する技術において再構成処理が知られている。例えば、特許文献１には、対象物に関する運動の複数の周期に亘って走査面を移動させつつ三次元空間内で複数の走査面を形成し、複数の走査面に対応した複数の画像の中から互いに周期的に対応した複数の画像を抽出し、抽出した複数の画像に基づいて対象物の表示画像を形成する旨の技術が記載されている。これにより、例えば周期内の複数の時相に亘って対象物が運動する様子を三次元的に示した表示画像を得ることができる。

ところが、具体的な処理内容次第で、再構成処理により形成される複数時相の表示画像内において、対象物の表示位置が全体的に動いてしまう場合がある。一方、対象物の表示位置を固定しようとすると表示画像の外枠が動いてしまう場合がある。

特開２０１０−２４６６３０号公報

上述した背景技術に鑑み、本願の発明者は、超音波における再構成処理について研究開発を重ねてきた。特に、再構成処理により得られる画像の表示に係る技術に注目した。

本発明は、その研究開発の過程において成されたものであり、その目的は、再構成処理により得られる画像の表示に係る改良技術を提供することにある。

上記目的にかなう好適な超音波診断装置は、周期的な運動をする対象物を含む三次元空間内において、当該運動の複数周期に亘って走査面を移動させつつ複数の走査面を形成するように超音波を送受する超音波送受部と、複数周期に亘って形成される複数の走査面に対応した複数の画像の中から、互いに周期的に対応した複数の画像を抽出して各再構成画像を形成することにより、周期内の互いに異なる複数の時相に対応した複数の再構成画像を得る再構成処理部と、各再構成画像内において前記対象物を基準とした注目領域を決定する領域決定部と、複数の時相に亘って複数の再構成画像から次々に得られる注目領域の画像を一定の位置に表示させる表示処理部と、を有することを特徴とする。

上記好適な超音波診断装置によれば、各再構成画像内において対象物を基準とした注目領域を決定し、複数の時相に亘って複数の再構成画像から次々に得られる注目領域の画像を一定の位置に表示させる旨の改良技術が提供される。例えば、対象物を基準としてその対象物の全体的な移動が抑制されるように注目領域を決定することにより、対象物の全体的な移動が抑制された表示画像を得ることができる。また、例えば、注目領域の大きさを一定にしておけば、注目領域に対応した表示画像の外枠が動くことも抑制できる。

望ましい具体例において、前記領域決定部は、各再構成画像内において前記対象物を含んだ注目領域以外の領域をマスク処理することにより注目領域を決定する、ことを特徴とする。

望ましい具体例において、前記領域決定部は、複数の再構成画像について注目領域の大きさを一定にして各再構成画像内から注目領域の画像を抽出する、ことを特徴とする。

望ましい具体例において、前記再構成処理部は、互いに周期的に対応した複数の画像をそれらの配列順に応じた位置に対応付けて再構成空間内に配置して、それら複数の画像から各再構成画像を形成し、これにより、複数の時相に亘る複数の再構成画像内において前記対象物が全体的に移動するように形成され、前記領域決定部は、複数の再構成画像に亘って前記対象物の移動に追従するように注目領域を移動させて各再構成画像内から注目領域の画像を抽出する、ことを特徴とする。

望ましい具体例において、前記再構成処理部は、互いに周期的に対応した複数の画像の各々をそれに対応した走査面の位置に対応付けて再構成空間内に配置して、それら複数の画像から各再構成画像を形成し、これにより、複数の時相に亘る複数の再構成画像内において前記対象物が全体的に移動しないように形成され、前記領域決定部は、複数の再構成画像に亘って前記対象物を含むように注目領域の位置を一定にして各再構成画像内から注目領域の画像を抽出する、ことを特徴とする。

また、上記目的にかなう好適なプログラムは、周期的な運動をする対象物を含む三次元空間内において、当該運動の複数周期に亘って走査面を移動させつつ複数の走査面を形成するように超音波を送受して得られる画像を処理するプログラムであって、複数周期に亘って形成される複数の走査面の各々から得られる画像のうちの互いに周期的に対応した複数の画像に基づいて各再構成画像を形成することにより、周期内の互いに異なる複数の時相に対応した複数の再構成画像を得る再構成処理機能と、各再構成画像内において前記対象物を基準とした注目領域を決定する領域決定機能と、複数の時相に亘って複数の再構成画像から次々に得られる注目領域の画像を一定の位置に表示させる表示処理機能と、をコンピュータに実現させる、ことを特徴とする。

上記プログラムは、例えば、ディスクやメモリなどのコンピュータが読み取り可能な記憶媒体に記憶され、その記憶媒体を介してコンピュータに提供される。もちろん、インターネット等の電気通信回線を介して上記プログラムがコンピュータに提供されてもよい。

本発明により、再構成処理で得られる画像の表示に係る改良技術が提供される。例えば本発明の好適な態様によれば、対象物の全体的な移動が抑制された表示画像を得ることができ、また、表示画像の外枠が動くことを抑制できる。

本発明の実施において好適な超音波診断装置の全体構成図である。図１の超音波診断装置による三次元的な走査を説明するための図である。図１の再構成処理部２０における処理を説明するための図である。再構成処理から表示画像を得るまでの具体例１を示す図である。再構成処理から表示画像を得るまでの具体例２を示す図である。

図１は、本発明の実施において好適な超音波診断装置（本超音波診断装置）の全体構成図である。プローブ１０は、診断の対象物を含む三次元空間内において超音波を送受波する。プローブ１０は、超音波を送受する複数の振動素子を備えており、複数の振動素子が送受信部１２によって送信制御されて送信ビームが形成される。また、複数の振動素子が対象物から反射された超音波を受波し、これにより得られた信号が送受信部１２へ出力され、送受信部１２が受信ビームを形成する。

プローブ１０は、超音波ビーム（送信ビームと受信ビーム）を三次元空間内において走査して立体的にエコーデータを収集する３Ｄプローブである。例えば、一次元的に配列された複数の振動素子（１Ｄアレイ振動子）によって電子的に形成される走査面を機械的に動かすことにより超音波ビームが三次元的に走査される。また、二次元的に配列された複数の振動素子（２Ｄアレイ振動子）を電子的に制御して超音波ビームを三次元的に走査してもよい。

送受信部１２は、プローブ１０が備える複数の振動素子の各々に対応した送信信号を供給することにより超音波の送信ビームを形成する。また、送受信部１２は、プローブ１０が備える複数の振動素子の各々から得られる受信信号に対して整相加算処理などを施すことにより超音波の受信ビームを形成し、受信ビームに沿って得られるエコーデータを出力する。

図２は、本超音波診断装置による三次元的な走査を説明するための図である。本超音波診断装置において、診断の対象物は例えば心臓などの周期的に運動する組織である。そして、例えば心臓などの周期的な運動をする対象物を含む三次元空間内において、対象物の運動の複数周期に亘って走査面Ｓを移動させつつ、例えば数百から数千枚程度の走査面Ｓを形成するように超音波が送受される。

つまり、図２に示すように、深さｒ方向に沿って形成される超音波ビームがビーム角度θ方向に走査されて走査面Ｓが形成され、さらに、スキャンφ方向に走査面Ｓを移動させつつ複数の走査面Ｓが形成される。例えば、電子的な走査により走査面Ｓが形成され、機械的に走査により走査面Ｓをスキャンφ方向に移動する。もちろん、電子的な走査により走査面Ｓをスキャンθ方向に移動させてもよい。対象物が胎児の心臓であれば、例えば約８秒程度で約２０心拍を含む期間に亘って、走査面Ｓがスキャンφ方向にゆっくりと移動して形成される。

図１に戻り、複数の走査面が形成されると、各走査面ごとにその走査面に関する断層画像データが収集され、複数の走査面に対応した複数の断層画像データが次々に前メモリ１４に記憶される。なお、各断層画像データは、ｒθ座標系（図２）で前メモリ１４に記憶されてもよいし、ｘｙ直交座標系に変換してから前メモリ１４に記憶されてもよい。

再構成処理部２０は、前メモリ１４に記憶された複数周期に亘る複数の断層画像データの中から、互いに周期的に対応した複数の断層画像データを抽出し、それら互いに周期的に対応した複数の断層画像データから各再構成画像データを形成して、後メモリ２６に記憶する。

図３は、再構成処理部２０における処理を説明するための図である。図３には、前メモリ１４（図１）に記憶されるデータと後メモリ２６（図１）に記憶されるデータの対応関係が示されている。図３において、「断層画像φｎ（ｎ＝１，２，３，・・・，６０）」は、スキャンφ方向（図２）の座標φｎに対応した断層画像データを意味している。

前メモリ１４には、φ方向に沿って次々に形成される複数の走査面に対応した複数の断層画像データが形成された順に記憶されている。つまり、前メモリ１４には、座標φｎの小さい方から、断層画像φ１，断層画像φ２，・・・，断層画像φ６０，・・・の順に複数の断層画像データが記憶されている。

再構成処理においては、代表的な時相の断層画像を基準にして、画像の並び替えが行われる。対象組織が心臓の場合には、例えば拡張末期が基準とされる。つまり、前メモリ１４に記憶された複数の断層画像のうち、拡張末期に対応した断層画像が基準画像として抽出される。図３において、断層画像φ１，断層画像φ１５，・・・，断層画像φ５１が複数の基準画像である。再構成処理部２０は、互いに周期的に対応した複数の断層画像データとして、まず、基準画像である断層画像φ１，断層画像φ１５，・・・，断層画像φ５１を抽出する。そして、抽出された断層画像φ１，断層画像φ１５，・・・，断層画像φ５１が一つのデータブロックとなって後メモリ２６内に記憶される。

次に、再構成処理部２０は、互いに周期的に対応した複数の断層画像データとして、複数の基準画像の各々に対してφ方向の正側に隣接する複数の断層画像を抽出する。つまり断層画像φ２，断層画像φ１６，・・・，断層画像φ５２が抽出され、これらが一つのデータブロックとなって後メモリ２６内に記憶される。

さらに、再構成処理部２０は、断層画像φ２，断層画像φ１６，・・・，断層画像φ５２の各々に対してφ方向の正側に隣接する複数の断層画像を抽出する。こうして、複数の基準画像の各々を起点として、互いに周期的に対応した複数の断層画像のデータブロックからなる再構成画像データが次々に抽出されて、後メモリ２６内に記憶される。

図１に戻り、マスク処理部３０は、各再構成画像データ内において対象物を基準とした注目領域を決定する。マスク処理部３０は、各再構成画像データ内において対象物を含んだ注目領域以外の領域をマスク処理することにより注目領域を決定する。また、表示処理部４０は、複数の時相に亘って複数の再構成画像データから次々に得られる注目領域の画像を一定の位置に配置した表示画像を形成し、その表示画像が表示部４２に表示される。なお、本超音波診断装置内の各部は、制御部５０により制御される。

そこで、本超音波診断装置における再構成処理から表示画像が得られるまでの処理について詳述する。その処理には以下に説明する２つの具体例がある。

図４は、再構成処理から表示画像を得るまでの具体例１を示す図である。具体例１の再構成処理においては、互いに周期的に対応した複数の断層画像データが、それらの配列順に応じた位置に対応付けて再構成空間内に配置され、それら複数の断層画像データから各時相の再構成画像データが形成される。

図４（１）には、時相０、時相Ｎ／２、時相Ｎ（Ｎは自然数）の各再構成画像データが示されている。図４（１）に示す各再構成画像データは、ｒθφ座標系（図２）におけるｒφ面の画像に対応している。

例えば、図３において後メモリ２６に記憶された断層画像φ１，断層画像φ１５，・・・，断層画像φ５１が、再構成空間内においてその空間の端（先頭位置）から順にφ方向に沿って配置されて、図４（１）の時相０の再構成画像データが形成される。また、例えば、図３の後メモリ２６に記憶された断層画像φ１０，断層画像φ２４，・・・，断層画像φ６０が、再構成空間内においてその空間の端（先頭位置）から順にφ方向に沿って配置されて、図４（１）の時相Ｎの再構成画像データが形成される。

この例において、時相０の再構成画像データ内で端（先頭位置）に配置される断層画像データは、断層画像φ１（図３）となり、時相Ｎの再構成画像データ内で同じく端（先頭位置）に配置される断層画像データは、断層画像φ１０（図３）となる。つまり、スキャンφ方向において互いに異なる位置にある断層画像φ１と断層画像φ１０が、再構成画像データ内において同じ位置に配置される。

このように、時相０と時相Ｎの再構成画像データは、互いに異なる位置（ずれた位置）から得られる断層画像によって形成されるため、時相０と時相Ｎの再構成画像データ内において、対象物である例えば心臓が異なる位置に配置される。この関係は時相０と時相Ｎとの間に限らず、時相が異なると心臓が全体的にずれた位置に配置されることになる。その結果、図４（１）に示す例のように、複数の時相に亘る複数の再構成画像データ内において、対象物である心臓が全体的に移動するように形成される。

そこで、図４（２）に示すマスク処理では、複数の再構成画像データに亘ってマスク領域を移動させることにより、対象物の移動に追従するように注目領域を移動させて各再構成画像データ内から注目領域の画像を抽出する。

図４（２）に示す時相０、時相Ｎ／２、時相Ｎの各再構成画像データ内において、斜線で示される領域がマスク領域であり、そのマスク領域以外の心臓を含んだ画像領域が注目領域である。注目領域の大きさ（面積）は、複数の時相に亘って一定とする。したがってマスク領域の大きさも、複数の時相に亘って一定となる。

図４（１）の再構成画像データにおいて互いに隣接する時相間では、図２のスキャン方向φに沿って断層画像一枚分のずれがある。つまり、図２において互いに隣接する断層画像間の角度Δφだけのずれがある。そこで、図４（２）に示すマスク処理では、各時相ごとに、スキャン方向φに沿って負方向に（図の左方向に）マスク領域をΔφずつ移動させる。

つまり、時相０においてスキャン方向φの負側（左側）を占めているマスク領域が、その後の各時相において、Δφずつ負方向に移動し、その移動によりはみ出した分だけスキャン方向φの正側（右側）にマスク領域が現れる。そして、中間の時相である時相Ｎ／２において正側と負側（右側と左側）のマスク領域の面積が等しくなる。さらに、時相が進んで、最終の時相である時相Ｎになると、スキャン方向φの正側（右側）に全てのマスク領域が移動する。

時相０から時相Ｎまでの間に各時相ごとにマスク領域をΔφだけ移動させて、図４（２）に示すマスク処理を実現するためには、マスク領域の大きさ（面積）をＮ×Δφとすればよい。また、再構成画像データの大きさ、例えば図４（１）に示す各再構成画像データの面積から、マスク領域の面積を差し引くことにより、注目領域の面積が得られる。複数の時相に亘って、マスク領域の面積Ｎ×Δφが一定であれば、注目領域の面積も一定となる。

こうして、図４（２）に示すように、複数の時相に亘ってマスク領域を移動させることにより、対象物である心臓の全体的な移動に追従するように注目領域が移動され、各時相ごとにその再構成画像データ内から注目領域の画像が抽出される。なお、各時相ごとに、θ方向（図２参照）に沿って並ぶ全てのｒφ面において同じマスク領域が適用される。

各時相ごとに注目領域が決定されると、図４（３）に示す表示処理において、各時相ごとに、注目領域の画像が表示領域の中央に配置されるように、再構成画像データを回転移動する。そして、表示領域の中央に配置された注目領域の画像のみが、図４（４）に示す表示画像として表示される。

対象物である心臓の全体的な動きに追従するように注目領域が決定されるため、図４（４）に示す表示画像では、複数の時相に亘って心臓の全体的な位置が一定とされる。また注目領域の大きさも一定とされているため、図４（４）に示す表示画像では、表示される画像の枠（つまり注目領域の外延）も固定された表示画像が得られる。

なお、表示画像は、例えばＸＹＺ直交座標系で形成されてもよい。例えば、図２に示したｒθ座標系の各断層画像データが、ｘｙ直交座標系に変換されてから図１の前メモリ１４に記憶され、さらに、図４（２）に示すマスク処理と並行して、ｙφ座標系（ｙはｒ方向に対応）からＹＺ直交座標系への変換が行われる。これにより、先に直交座標系に変換された各断層画像データのｘ方向をそのままＸ方向とすれば、ＸＹＺ直交座標系で表示画像が得られることになる。

図５は、再構成処理から表示画像を得るまでの具体例２を示す図である。具体例２の再構成処理においては、互いに周期的に対応した複数の断層画像データの各々が、その断層画像データが得られた走査面の位置に対応付けて再構成空間内に配置され、それら複数の断層画像データから各時相の再構成画像データが形成される。

図５（１）には、時相０、時相Ｎ／２、時相Ｎ（Ｎは自然数）の各再構成画像データが示されている。図５（１）に示す各再構成画像データは、ｒθφ座標系（図２）におけるｒφ面の画像に対応している。

例えば、図３において後メモリ２６に記憶された断層画像φ１，断層画像φ１５，・・・，断層画像φ５１が、再構成空間内において、それぞれ座標φ１，φ１５，・・・，φ５１に配置されて、図５（１）の時相０の再構成画像データが形成される。また、例えば図３の後メモリ２６に記憶された断層画像φ１０，断層画像φ２４，・・・，断層画像φ６０が、再構成空間内において、それぞれ座標φ１０，φ２４，・・・，φ６０に配置されて、図５（１）の時相Ｎの再構成画像データが形成される。

図５の例では、再構成画像データ内において走査面の位置に応じて断層画像データが配置されるため、心臓の画像も再構成処理前の当初の位置に配置される。したがって、再構成処理前において心臓が全体的に移動していなければ、複数の時相に亘る複数の再構成画像データ内においても心臓は全体的に移動することはない。

ところが、図５の例においては、時相０の再構成画像データの端に配置される断層画像データは、断層画像φ１（図３）となり、時相Ｎの再構成画像データの端に配置される断層画像データは、断層画像φ１０（図３）となる。そして、断層画像φ１と断層画像φ１０は、スキャンφ方向において互いに異なる位置にあるため、再構成画像データ内においても互いに異なる位置に配置される。つまり、時相０の再構成画像データの端の位置と、時相Ｎの再構成画像データの端の位置とが、互いに異なることになる。

このように、時相０と時相Ｎの再構成画像データは、互いに画像の端の位置がずれてしまう。この関係は時相０と時相Ｎとの間に限らず、時相が異なると再構成画像データの端の位置がずれることになる。その結果、図５（１）に示す例のように、複数の時相に亘る複数の再構成画像データにおいて、再構成画像データの端の位置、つまり再構成画像データの外枠がずれてしまう。

そこで、図５（２）に示すマスク処理では、複数の再構成画像データに亘って、外枠のずれてしまう部分にマスク領域を設定し、対象物を含む注目領域の位置を中央に固定して各再構成画像データ内から注目領域の画像を抽出する。

図５（２）に示す時相０、時相Ｎ／２、時相Ｎの各再構成画像データ内において、斜線で示される領域がマスク領域であり、そのマスク領域以外の心臓を含んだ画像領域が注目領域である。注目領域の大きさ（面積）は、複数の時相に亘って一定とする。したがってマスク領域の大きさも、複数の時相に亘って一定となる。また、注目領域の位置も複数の時相に亘って一定とする。なお、各時相ごとに、θ方向（図２参照）に沿って並ぶ全てのｒφ面において同じマスク領域が適用される。

各時相ごとに注目領域が決定されると、表示処理においては、注目領域の画像が表示領域の中央に配置されるようにそのまま維持され、そして、表示領域の中央に配置された注目領域の画像のみが、図５（３）に示す表示画像として表示される。

再構成処理において走査面の位置が考慮されているため、図５（３）に示す表示画像では、対象物である心臓の位置も再構成処理前の当初の位置となる。また、表示領域の中央に配置された注目領域の画像のみが表示されるため、図５（３）に示す表示画像では、表示される画像の枠（つまり注目領域の外延）も固定された表示画像が得られる。

なお、表示画像は、例えばＸＹＺ直交座標系で形成されてもよい。例えば、図２に示したｒθ座標系の各断層画像データが、ｘｙ直交座標系に変換されてから図１の前メモリ１４に記憶され、さらに、図５（２）に示すマスク処理と並行して、ｙφ座標系（ｙはｒ方向に対応）からＹＺ直交座標系への変換が行われる。これにより、先に直交座標系に変換された各断層画像データのｘ方向をそのままＸ方向とすれば、ＸＹＺ直交座標系で表示画像が得られることになる。

以上、本発明の好適な実施形態である超音波診断装置について説明したが、例えば、図４と図５を利用して説明した再構成処理から表示画像を得るまでの処理の一部または全てに対応したプログラムにより、図１に示した再構成処理部２０とマスク処理部３０と表示処理部４０の機能の一部または全てをコンピュータで実現し、そのコンピュータを超音波画像処理装置として機能させてもよい。

なお、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。

１０プローブ、１２送受信部、２０再構成処理部、３０マスク処理部、４０表示処理部。

Claims

周期的な運動をする対象物を含む三次元空間内において、当該運動の複数周期に亘って走査面を移動させつつ複数の走査面を形成するように超音波を送受する超音波送受部と、
複数周期に亘って形成される複数の走査面に対応した複数の画像の中から、互いに周期的に対応した複数の画像を抽出して各再構成画像を形成することにより、周期内の互いに異なる複数の時相に対応した複数の再構成画像を得る再構成処理部と、
各再構成画像内において前記対象物を基準とした注目領域を決定する領域決定部と、
複数の時相に亘って複数の再構成画像から次々に得られる注目領域の画像を一定の位置に表示させる表示処理部と、
を有する、
ことを特徴とする超音波診断装置。
請求項１に記載の超音波診断装置において、
前記領域決定部は、各再構成画像内において前記対象物を含んだ注目領域以外の領域をマスク処理することにより注目領域を決定する、
ことを特徴とする超音波診断装置。
請求項１または２に記載の超音波診断装置において、
前記領域決定部は、複数の再構成画像について注目領域の大きさを一定にして各再構成画像内から注目領域の画像を抽出する、
ことを特徴とする超音波診断装置。
請求項１から３のいずれか１項に記載の超音波診断装置において、
前記再構成処理部は、互いに周期的に対応した複数の画像をそれらの配列順に応じた位置に対応付けて再構成空間内に配置して、それら複数の画像から各再構成画像を形成し、これにより、複数の時相に亘る複数の再構成画像内において前記対象物が全体的に移動するように形成され、
前記領域決定部は、複数の再構成画像に亘って前記対象物の移動に追従するように注目領域を移動させて各再構成画像内から注目領域の画像を抽出する、
ことを特徴とする超音波診断装置。
請求項１から３のいずれか１項に記載の超音波診断装置において、
前記再構成処理部は、互いに周期的に対応した複数の画像の各々をそれに対応した走査面の位置に対応付けて再構成空間内に配置して、それら複数の画像から各再構成画像を形成し、これにより、複数の時相に亘る複数の再構成画像内において前記対象物が全体的に移動しないように形成され、
前記領域決定部は、複数の再構成画像に亘って前記対象物を含むように注目領域の位置を一定にして各再構成画像内から注目領域の画像を抽出する、
ことを特徴とする超音波診断装置。
周期的な運動をする対象物を含む三次元空間内において、当該運動の複数周期に亘って走査面を移動させつつ複数の走査面を形成するように超音波を送受して得られる画像を処理するプログラムであって、
複数周期に亘って形成される複数の走査面の各々から得られる画像のうちの互いに周期的に対応した複数の画像に基づいて各再構成画像を形成することにより、周期内の互いに異なる複数の時相に対応した複数の再構成画像を得る再構成処理機能と、
各再構成画像内において前記対象物を基準とした注目領域を決定する領域決定機能と、
複数の時相に亘って複数の再構成画像から次々に得られる注目領域の画像を一定の位置に表示させる表示処理機能と、
をコンピュータに実現させる、
ことを特徴とするプログラム。