JP2003164431A - 血流速度測定装置 - Google Patents
血流速度測定装置Info
- Publication number
- JP2003164431A JP2003164431A JP2001365419A JP2001365419A JP2003164431A JP 2003164431 A JP2003164431 A JP 2003164431A JP 2001365419 A JP2001365419 A JP 2001365419A JP 2001365419 A JP2001365419 A JP 2001365419A JP 2003164431 A JP2003164431 A JP 2003164431A
- Authority
- JP
- Japan
- Prior art keywords
- image
- blood flow
- images
- flow velocity
- pixels
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Eye Examination Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
(57)【要約】
【課題】 血流速度を演算するための複数枚の画像のず
れを補正し、高精度で鮮明な血流速度分布画像を得るこ
とができる血流速度測定装置を提供する。 【解決手段】 本発明の血流速度測定装置は、生体組織
の血球にレーザ光を照射するレーザ光照射系１と、生体
組織からの反射光を検出する多数の画素からなる受光部
５を有する受光系２と、受光部５からの信号に基づき所
定時間間隔で連続的に複数の画像を取り込む画像取り込
み部１２と、複数の画像を記憶する画像記憶部１６と、
記憶された複数画像の対応する各画素の出力信号の時間
的変化から生体組織内の血流速度を演算する演算部１７
とからなる血流速度測定装置であり、前記演算部は前記
複数画像の相互の位置ずれ情報を検出する位置ずれ検出
部を有し、位置ずれ情報に基づき記憶された複数画像の
位置ずれを補正する。
れを補正し、高精度で鮮明な血流速度分布画像を得るこ
とができる血流速度測定装置を提供する。 【解決手段】 本発明の血流速度測定装置は、生体組織
の血球にレーザ光を照射するレーザ光照射系１と、生体
組織からの反射光を検出する多数の画素からなる受光部
５を有する受光系２と、受光部５からの信号に基づき所
定時間間隔で連続的に複数の画像を取り込む画像取り込
み部１２と、複数の画像を記憶する画像記憶部１６と、
記憶された複数画像の対応する各画素の出力信号の時間
的変化から生体組織内の血流速度を演算する演算部１７
とからなる血流速度測定装置であり、前記演算部は前記
複数画像の相互の位置ずれ情報を検出する位置ずれ検出
部を有し、位置ずれ情報に基づき記憶された複数画像の
位置ずれを補正する。
Description
【０００１】
【発明の属する技術分野】本発明は、血球を有する生体
組織にレーザ光を照射し、その血球から反射されたスペ
ックル信号に基づき血流速度を測定するための血流速度
測定装置に関する。
組織にレーザ光を照射し、その血球から反射されたスペ
ックル信号に基づき血流速度を測定するための血流速度
測定装置に関する。
【０００２】
【従来の技術】従来、被検眼の眼底等の生体組織の血球
にレーザ光を照射して、その血球からの反射光により形
成された画像を固体撮像装置（ＣＣＤ）等のイメージセ
ンサー上に導き、この画像を連続的に所定時間間隔で多
数枚取り込み・記憶し、その記憶された多数の画像の中
から所定枚数の画像を選択し、各画像の各画素における
出力の時間変動量を積算した値を算出し、この値から血
球の速度（血流速度）を算出する血流速度測定装置が知
られている。
にレーザ光を照射して、その血球からの反射光により形
成された画像を固体撮像装置（ＣＣＤ）等のイメージセ
ンサー上に導き、この画像を連続的に所定時間間隔で多
数枚取り込み・記憶し、その記憶された多数の画像の中
から所定枚数の画像を選択し、各画像の各画素における
出力の時間変動量を積算した値を算出し、この値から血
球の速度（血流速度）を算出する血流速度測定装置が知
られている。
【０００３】この種の血流速度測定装置では、各画素の
出力変動量が血球の移動速度に対応するので、この算出
された各画素の出力変動量値に基づき、生体組織での血
流速度分布を二次元画像としてモニター画面上に表示し
ている。
出力変動量が血球の移動速度に対応するので、この算出
された各画素の出力変動量値に基づき、生体組織での血
流速度分布を二次元画像としてモニター画面上に表示し
ている。
【０００４】
【発明が解決しようとする課題】しかしながら、この種
の血流速度測定装置では、所定数枚の一連の画像につい
て演算を行っているが、これらを取り込む間に測定対象
画像の対応する画素の位置関係がずれ、それらを演算し
た結果、精細な血流速度分布画像を得られないという欠
点を有している。
の血流速度測定装置では、所定数枚の一連の画像につい
て演算を行っているが、これらを取り込む間に測定対象
画像の対応する画素の位置関係がずれ、それらを演算し
た結果、精細な血流速度分布画像を得られないという欠
点を有している。
【０００５】特に、眼底像をCCD上に導く場合、固視微
動或いは眼球の運動の影響により、CCD上に形成される
画像は常にその位置が変動し、鮮明な血流速度分布画像
を得ることが非常に困難となっている。
動或いは眼球の運動の影響により、CCD上に形成される
画像は常にその位置が変動し、鮮明な血流速度分布画像
を得ることが非常に困難となっている。
【０００６】本発明は、この従来の問題点を解決するこ
とを目的とし、血流速度を演算するための複数枚の画像
のずれを補正し、高精度で鮮明な血流速度分布画像を得
ることができる血流速度測定装置を提供することを目的
とする。
とを目的とし、血流速度を演算するための複数枚の画像
のずれを補正し、高精度で鮮明な血流速度分布画像を得
ることができる血流速度測定装置を提供することを目的
とする。
【０００７】
【課題を解決するための手段】請求項１に記載の血流速
度測定装置は、生体組織の血球にレーザ光を照射するレ
ーザ光照射系と、前記生体組織からの反射光を検出する
多数の画素からなる受光部を有する受光系と、前記受光
部からの信号に基づき所定時間間隔で連続的に複数の画
像を取り込む画像取り込み部と、前記複数の画像を記憶
する画像記憶部と、該記憶された複数画像の対応する各
画素の出力信号の時間的変化から生体組織内の血流速度
を演算する演算部とからなる血流速度測定装置におい
て、前記演算部は前記複数画像の相互の位置ずれ情報を
検出する位置ずれ検出部を有し、前記位置ずれ情報に基
づき記憶された複数画像の位置ずれを補正することを特
徴とする。
度測定装置は、生体組織の血球にレーザ光を照射するレ
ーザ光照射系と、前記生体組織からの反射光を検出する
多数の画素からなる受光部を有する受光系と、前記受光
部からの信号に基づき所定時間間隔で連続的に複数の画
像を取り込む画像取り込み部と、前記複数の画像を記憶
する画像記憶部と、該記憶された複数画像の対応する各
画素の出力信号の時間的変化から生体組織内の血流速度
を演算する演算部とからなる血流速度測定装置におい
て、前記演算部は前記複数画像の相互の位置ずれ情報を
検出する位置ずれ検出部を有し、前記位置ずれ情報に基
づき記憶された複数画像の位置ずれを補正することを特
徴とする。
【０００８】請求項２に記載の血流速度測定装置は、前
記位置ずれ検出部が、各画素出力を該画素の近傍の複数
の画素の出力に基づき演算した値とした暫定マッチング
用画像を作成することを特徴とする。
記位置ずれ検出部が、各画素出力を該画素の近傍の複数
の画素の出力に基づき演算した値とした暫定マッチング
用画像を作成することを特徴とする。
【０００９】請求項３に記載の血流速度測定装置は、前
記位置ずれ検出部が、２つの前記暫定マッチング用画像
のマッチングをとることにより各画像の相互の位置ずれ
を検出することを特徴とする。
記位置ずれ検出部が、２つの前記暫定マッチング用画像
のマッチングをとることにより各画像の相互の位置ずれ
を検出することを特徴とする。
【００１０】請求項４に記載の血流速度測定装置は、前
記画像のマッチングが、一方の画像を他方の画像上でと
びとびの所定画素毎に所定ステップ数移動させて各移動
ステップで求めた相関値から第１画像移動量を算出する
第１測定ステップと、前記第１測定ステップで算出され
た第１画像移動量を基準として１画素毎に所定ステップ
数移動させて各移動ステップで求めた相関値から第２画
像移動量を算出する第２測定ステップとからなり、前記
第１画像移動量と第２画像移動量とから全体の画像移動
量を算出することを特徴とする。
記画像のマッチングが、一方の画像を他方の画像上でと
びとびの所定画素毎に所定ステップ数移動させて各移動
ステップで求めた相関値から第１画像移動量を算出する
第１測定ステップと、前記第１測定ステップで算出され
た第１画像移動量を基準として１画素毎に所定ステップ
数移動させて各移動ステップで求めた相関値から第２画
像移動量を算出する第２測定ステップとからなり、前記
第１画像移動量と第２画像移動量とから全体の画像移動
量を算出することを特徴とする。
【００１１】請求項５に記載の血流速度測定装置は、前
記第１測定ステップが、移動画素数を順次少なくする複
数の測定ステップからなることを特徴とする。
記第１測定ステップが、移動画素数を順次少なくする複
数の測定ステップからなることを特徴とする。
【００１２】請求項６に記載の血流速度測定装置は、前
記第１測定ステップは、所定間隔離間した各画素間で相
関をとることを特徴とする。
記第１測定ステップは、所定間隔離間した各画素間で相
関をとることを特徴とする。
【００１３】請求項７に記載の血流速度測定装置は、前
記受光部が固体撮像素子であり、前記画像取り込み部
は、奇数走査線に基づく奇数フィールド画像と偶数走査
線の信号に基づく偶数フィールド画像とを１組として取
り込み、前記画像記憶部は前記奇数フィールド画像と偶
数フィールド画像との組み合わせからなる合成フレーム
画像として記憶することを特徴とする。
記受光部が固体撮像素子であり、前記画像取り込み部
は、奇数走査線に基づく奇数フィールド画像と偶数走査
線の信号に基づく偶数フィールド画像とを１組として取
り込み、前記画像記憶部は前記奇数フィールド画像と偶
数フィールド画像との組み合わせからなる合成フレーム
画像として記憶することを特徴とする。
【００１４】
【発明の実施の形態】図１は血流測定装置の光学系の概
要を示すもので、１はレーザー光照射系、２は受光系、
Ｅは被検眼である。レーザー光照射系１のレーザー光
は、例えば、ハーフミラー３を介して被検眼Ｅの生体組
織としての例えば眼底Ｅｒに照射される。
要を示すもので、１はレーザー光照射系、２は受光系、
Ｅは被検眼である。レーザー光照射系１のレーザー光
は、例えば、ハーフミラー３を介して被検眼Ｅの生体組
織としての例えば眼底Ｅｒに照射される。
【００１５】受光系２は受光レンズ４、受光部としての
ＣＣＤ（固体撮像素子）５、増幅回路６を有する。眼底
Ｅｒからのレーザー反射光は受光レンズ４によりＣＣＤ
５に生体組織像として結像される。ＣＣＤ５はその受光
面上に多数の画素を有し、受光レンズ４により結像され
た生体組織像を電気信号に変換し、フレーム蓄積方式で
信号電荷を読み出して映像信号として出力する。その映
像信号は、信号増幅回路６により増幅され、信号増幅回
路６により増幅された映像信号は利得制御等を行うアナ
ログ信号処理手段７に出力され、A/D変換器８によりデ
ジタル信号に変換される。
ＣＣＤ（固体撮像素子）５、増幅回路６を有する。眼底
Ｅｒからのレーザー反射光は受光レンズ４によりＣＣＤ
５に生体組織像として結像される。ＣＣＤ５はその受光
面上に多数の画素を有し、受光レンズ４により結像され
た生体組織像を電気信号に変換し、フレーム蓄積方式で
信号電荷を読み出して映像信号として出力する。その映
像信号は、信号増幅回路６により増幅され、信号増幅回
路６により増幅された映像信号は利得制御等を行うアナ
ログ信号処理手段７に出力され、A/D変換器８によりデ
ジタル信号に変換される。
【００１６】９はタイミングパルス発生器であり、１０
は電子シャッター制御手段、１１は固体撮像素子駆動手
段であり、タイミングパルス発生器９は電子シャッター
制御手段１０と信号選択手段１２とにタイミングパルス
を出力する。固体撮像素子駆動手段１１はタイミングパ
ルスに基づき駆動される。
は電子シャッター制御手段、１１は固体撮像素子駆動手
段であり、タイミングパルス発生器９は電子シャッター
制御手段１０と信号選択手段１２とにタイミングパルス
を出力する。固体撮像素子駆動手段１１はタイミングパ
ルスに基づき駆動される。
【００１７】信号選択手段１２にはＡ／Ｄ変換器８によ
りＡ／Ｄ変換された映像信号としてのデジタル信号が入
力され、タイミングパスル発生器９からのタイミングパ
ルスに基づいて、信号選択手段１２は、A/D変換器８か
らの映像信号のうち、偶数フィールドの画像を第１画像
記録手段１３へ送信し、奇数フィールドの画像を第２画
像記録手段１４へ送信し、第１画像記録手段１３には偶
数フィールドの画像が記憶され、第２画像記録手段１４
には奇数フィールドの画像が記録される。その信号選択
手段１２、第１画像記録手段１３、第２画像記録手段１
４は、所定時間間隔で複数枚の画像を取り込む画像取り
込み部の一部として機能する。
りＡ／Ｄ変換された映像信号としてのデジタル信号が入
力され、タイミングパスル発生器９からのタイミングパ
ルスに基づいて、信号選択手段１２は、A/D変換器８か
らの映像信号のうち、偶数フィールドの画像を第１画像
記録手段１３へ送信し、奇数フィールドの画像を第２画
像記録手段１４へ送信し、第１画像記録手段１３には偶
数フィールドの画像が記憶され、第２画像記録手段１４
には奇数フィールドの画像が記録される。その信号選択
手段１２、第１画像記録手段１３、第２画像記録手段１
４は、所定時間間隔で複数枚の画像を取り込む画像取り
込み部の一部として機能する。
【００１８】その第１画像記録手段１３、第２画像記録
手段１４に記録された奇数フィールドの画像データと偶
数フィールドの画像データとは、画像合成手段１５によ
り合成され、１／３０秒間隔で撮影した１フレームの画
像データとされ、その１フレームの画像データは画像記
憶部としての画像記録器１６に記憶される。
手段１４に記録された奇数フィールドの画像データと偶
数フィールドの画像データとは、画像合成手段１５によ
り合成され、１／３０秒間隔で撮影した１フレームの画
像データとされ、その１フレームの画像データは画像記
憶部としての画像記録器１６に記憶される。
【００１９】この画像記憶器１６に記憶された画像信号
は演算処理部１７に入力され、演算処理部１７は、後述
する演算処理を行う。なお、１８はＴＶモニターであ
る。
は演算処理部１７に入力され、演算処理部１７は、後述
する演算処理を行う。なお、１８はＴＶモニターであ
る。
【００２０】ところで、複数枚の画像を多数枚撮影する
場合、ＣＣＤ受光面上での生体組織像の位置がずれる。
これは、被検眼等に固視微動等があるからである。すな
わち、例えば、図２（ａ）に示すように、画像Ｇ０が時
刻ｔ０において撮影された画像データであるとして、図
２（ｂ）で示す画像Ｇ１が時刻ｔ１において撮影された
画像データとするとき、時刻ｔ０で撮影された画像Ｇ０
の点ｇ１と時刻ｔ１において撮影された画像Ｇ１の点ｇ
１’とが対応関係にあるものであるが、ＣＣＤ５の受光
面の画素上での位置がずれている。なお、図２におい
て、破線は、血管の輪郭線ｇ２を示す。なお、時刻ｔ０
と時刻ｔ１との時間間隔Δは３０分の１秒とする。
場合、ＣＣＤ受光面上での生体組織像の位置がずれる。
これは、被検眼等に固視微動等があるからである。すな
わち、例えば、図２（ａ）に示すように、画像Ｇ０が時
刻ｔ０において撮影された画像データであるとして、図
２（ｂ）で示す画像Ｇ１が時刻ｔ１において撮影された
画像データとするとき、時刻ｔ０で撮影された画像Ｇ０
の点ｇ１と時刻ｔ１において撮影された画像Ｇ１の点ｇ
１’とが対応関係にあるものであるが、ＣＣＤ５の受光
面の画素上での位置がずれている。なお、図２におい
て、破線は、血管の輪郭線ｇ２を示す。なお、時刻ｔ０
と時刻ｔ１との時間間隔Δは３０分の１秒とする。
【００２１】通常、テンプレート法により、二つの画像
を重ねて、一方の画像を基準にして他方の画像を順次ず
らしてゆき、各画素毎にその差分を求め、その差分の総
和を算出して、これらを演算に使用した総画素数で割っ
た値の最小値を画像が一致した状態として求めるもので
あるが、画像データにはスペックルパターン（スペック
ル信号）が混在しているので、各画像の輪郭線ｇ２がは
っきりせず、画像記憶器１６に記憶されている画像デー
タをそのまま用いたとしても、ずれ量Δｘ、Δｙを求め
ることができない。
を重ねて、一方の画像を基準にして他方の画像を順次ず
らしてゆき、各画素毎にその差分を求め、その差分の総
和を算出して、これらを演算に使用した総画素数で割っ
た値の最小値を画像が一致した状態として求めるもので
あるが、画像データにはスペックルパターン（スペック
ル信号）が混在しているので、各画像の輪郭線ｇ２がは
っきりせず、画像記憶器１６に記憶されている画像デー
タをそのまま用いたとしても、ずれ量Δｘ、Δｙを求め
ることができない。
【００２２】従って、本発明では、まず、演算処理部１
７は、画像の輪郭線ｇ２の像をはっきりさせて、画像マ
ッチングを容易にするために暫定マッチング画像を作成
する。
７は、画像の輪郭線ｇ２の像をはっきりさせて、画像マ
ッチングを容易にするために暫定マッチング画像を作成
する。
【００２３】図３は、例えばその合成画像Ｇ０の一部の
領域を示すもので、Ｏは奇数フィールドの画像に相当す
る走査線、Eは偶数フィールドの画像に相当する走査線
である。この合成画像Ｇ０から奇数フィールドにおける
暫定マッチング用画像と、偶数フィールドにおける暫定
マッチング#用画像とをそれぞれ演算する。この演算は
演算処理部１７の位置ずれ検出部によって行う。
領域を示すもので、Ｏは奇数フィールドの画像に相当す
る走査線、Eは偶数フィールドの画像に相当する走査線
である。この合成画像Ｇ０から奇数フィールドにおける
暫定マッチング用画像と、偶数フィールドにおける暫定
マッチング#用画像とをそれぞれ演算する。この演算は
演算処理部１７の位置ずれ検出部によって行う。
【００２４】いま、偶数フィールドの暫定マッチング用
画像を演算により求める場合について説明する。
画像を演算により求める場合について説明する。
【００２５】図３に示すｐｉｊの画素について演算を行
う場合には、画素ｐｉｊを中心として縦方向と横方向と
について一画素おきの画素ｐｉｊ１〜ｐｉｊ８の８個の
出力値（画素値）に基づいて演算を行う。なお、画素ｐ
ｉｊの（）内の文字はその画素の出力値（グレイスケー
ル値）を表す。
う場合には、画素ｐｉｊを中心として縦方向と横方向と
について一画素おきの画素ｐｉｊ１〜ｐｉｊ８の８個の
出力値（画素値）に基づいて演算を行う。なお、画素ｐ
ｉｊの（）内の文字はその画素の出力値（グレイスケー
ル値）を表す。
【００２６】すなわち、縦方向については、奇数の走査
線の画素を飛ばして一つおきの画素ｐｉｊ１、ｐｉｊ
２、ｐｉｊ３、ｐｉｊ４の出力値を用い、横方向につい
ては、一画素を飛ばして一つおきの画素ｐｉｊ５、ｐｉ
ｊ６、ｐｉｊ７、ｐｉｊ８の出力値を用いて演算を行
う。
線の画素を飛ばして一つおきの画素ｐｉｊ１、ｐｉｊ
２、ｐｉｊ３、ｐｉｊ４の出力値を用い、横方向につい
ては、一画素を飛ばして一つおきの画素ｐｉｊ５、ｐｉ
ｊ６、ｐｉｊ７、ｐｉｊ８の出力値を用いて演算を行
う。
【００２７】その演算式は以下の通りである。
【００２８】ｘ＝（ｘ＋ａ）＋（ｘ＋ｂ）＋………＋
（ｘ＋ｈ）／｜ｘ−ａ｜＋｜ｘ−ｂ｜＋…＋｜ｘ−ｈ｜ このｘ値が暫定マッチング用画像のｐｉｊ画像の出力と
なる。
（ｘ＋ｈ）／｜ｘ−ａ｜＋｜ｘ−ｂ｜＋…＋｜ｘ−ｈ｜ このｘ値が暫定マッチング用画像のｐｉｊ画像の出力と
なる。
【００２９】この式に基づいて、偶数フィールドについ
て、ｉ＝０、ｊ＝０の第１番目の画素からｉ＝ｎ、ｊ＝
ｍの最後の画素までについて行う。
て、ｉ＝０、ｊ＝０の第１番目の画素からｉ＝ｎ、ｊ＝
ｍの最後の画素までについて行う。
【００３０】奇数フィールドの暫定マッチング用画像を
演算により求める場合についても、同様に上記式を用い
る。
演算により求める場合についても、同様に上記式を用い
る。
【００３１】このようにして、時刻ｔ０において撮影さ
れた画像に対応する暫定マッチング用の合成画像が求め
られる。同様に、時刻ｔ１において撮影された画像に対
応する暫定マッチング用の合成画像が求められ、画像記
録部１６に記憶された残りの複数枚の合成画像について
も、同様の処理によって暫定マッチング用画像が求めら
れる。
れた画像に対応する暫定マッチング用の合成画像が求め
られる。同様に、時刻ｔ１において撮影された画像に対
応する暫定マッチング用の合成画像が求められ、画像記
録部１６に記憶された残りの複数枚の合成画像について
も、同様の処理によって暫定マッチング用画像が求めら
れる。
【００３２】この一連の暫定マッチング用画像の作成処
理によって、スペックルパターンが大略除去されて、輪
郭線ｇ２が浮き上がった暫定マッチング用画像が得られ
る。
理によって、スペックルパターンが大略除去されて、輪
郭線ｇ２が浮き上がった暫定マッチング用画像が得られ
る。
【００３３】すなわち、位置ずれ検出部は、その各画素
の出力を、その各画素近傍の出力に基づいて補正し、こ
れによって、輪郭線が強調される。
の出力を、その各画素近傍の出力に基づいて補正し、こ
れによって、輪郭線が強調される。
【００３４】次に、これらの暫定マッチング用画像に基
づいて、所定時間間隔毎の複数枚の画像の移動量を測定
する方法について述べる。
づいて、所定時間間隔毎の複数枚の画像の移動量を測定
する方法について述べる。
【００３５】図４は、時間間隔Δ(1/30秒間隔で)での２
枚の暫定マッチング用画像Ｇ０’と暫定マッチング用画
像Ｇ１’とを示している。暫定マッチング用画像Ｇ０’
の各画素の番地と暫定マッチング用画像Ｇ１’の番地と
が一致しているときを基準位置とする。
枚の暫定マッチング用画像Ｇ０’と暫定マッチング用画
像Ｇ１’とを示している。暫定マッチング用画像Ｇ０’
の各画素の番地と暫定マッチング用画像Ｇ１’の番地と
が一致しているときを基準位置とする。
【００３６】ここでは、この２つの暫定マッチング用画
像の相関が最も高い位置を求め、暫定マッチング用画像
の移動量（ずれ量）を求めるテンプレート手法を採用す
る。このテンプレート手法では、一方の画像を他方の画
像上で所定ステップ毎に移動させ、それぞれの位置での
画像の相関度を求める。
像の相関が最も高い位置を求め、暫定マッチング用画像
の移動量（ずれ量）を求めるテンプレート手法を採用す
る。このテンプレート手法では、一方の画像を他方の画
像上で所定ステップ毎に移動させ、それぞれの位置での
画像の相関度を求める。
【００３７】この相関度を検出する際、例えば、移動さ
せた位置での２つの画像の対応する各画素の出力の差の
積算値を求め、この積算に用いた画素の個数で割り算し
たものを用いる。この積算値が小さいほど、相関度が高
いと判定される。ここで精度の高いマッチングを行う場
合、１画素分づつ移動させて画像の全ての各画素での出
力の差の絶対値の積算値を算出して１画素での精度での
マッチングを行うことが必要であるが、この演算を行う
と極めて膨大な演算を必要とするため、本実施例では、
下記の順序でマッチングを行い、演算量を少なくしてい
る。 (１) 出力の差の積算値を求めるための画素は、図５に
拡大して示すように、４画素刻みとし、４画素刻みで演
算に用いる画素を選択する。ここでは、基準位置での演
算に用いる画素が斜線で示されている。例えば、画素の
番地（ｇ０，０）を演算に用いるとすると、ｘ方向につ
いて次ぎに演算に用いる画素の番地は（ｇ４，０）であ
り、ｙ方向について次ぎに演算に用いる画素の番地は
（ｇ０，４）であり、一般に演算に用いる画素の番地は
（ｇｉ＋４，ｊ＋４）である。
せた位置での２つの画像の対応する各画素の出力の差の
積算値を求め、この積算に用いた画素の個数で割り算し
たものを用いる。この積算値が小さいほど、相関度が高
いと判定される。ここで精度の高いマッチングを行う場
合、１画素分づつ移動させて画像の全ての各画素での出
力の差の絶対値の積算値を算出して１画素での精度での
マッチングを行うことが必要であるが、この演算を行う
と極めて膨大な演算を必要とするため、本実施例では、
下記の順序でマッチングを行い、演算量を少なくしてい
る。 (１) 出力の差の積算値を求めるための画素は、図５に
拡大して示すように、４画素刻みとし、４画素刻みで演
算に用いる画素を選択する。ここでは、基準位置での演
算に用いる画素が斜線で示されている。例えば、画素の
番地（ｇ０，０）を演算に用いるとすると、ｘ方向につ
いて次ぎに演算に用いる画素の番地は（ｇ４，０）であ
り、ｙ方向について次ぎに演算に用いる画素の番地は
（ｇ０，４）であり、一般に演算に用いる画素の番地は
（ｇｉ＋４，ｊ＋４）である。
【００３８】また、画像の移動ステップ幅Ｗは、図６に
示すように、基準位置Ｏ０を基準に８画素刻みとする。
画像の移動のステップ数は基準位置Ｏ0を基準としてプ
ラスマイナス１２ステップとする。例えば、図６（ａ）
に示すように、暫定マッチング用画像Ｇ１’を右横方向
（＋ｘ方向）に８画素刻みに移動させることを１２回繰
り返し、次いで、基準位置Ｏ0に戻って暫定マッチング
用画像Ｇ１’を左横方向（−ｘ方向）に８画素刻みに移
動させることを１２回繰り返した後、図６（ｂ）に示す
ように、暫定マッチング用画像Ｇ１’を下縦方向（＋ｙ
方向）に８画素刻み移動させる。次いで、暫定マッチン
グ用画像Ｇ１’を右横方向（＋ｘ方向）に８画素刻みに
移動させることを１２回繰り返し、基準位置Ｏ0に戻っ
て暫定マッチング用画像Ｇ１’を左横方向（−ｘ方向）
に８画素刻みに移動させることを１２回繰り返す。この
暫定マッチング用画像Ｇ１’の８画素刻みの移動を＋ｙ
方向に１２回、−ｙ方向に１２回繰り返すと、総計６２
５回（＝２５×２５）の移動処理が行われる。
示すように、基準位置Ｏ０を基準に８画素刻みとする。
画像の移動のステップ数は基準位置Ｏ0を基準としてプ
ラスマイナス１２ステップとする。例えば、図６（ａ）
に示すように、暫定マッチング用画像Ｇ１’を右横方向
（＋ｘ方向）に８画素刻みに移動させることを１２回繰
り返し、次いで、基準位置Ｏ0に戻って暫定マッチング
用画像Ｇ１’を左横方向（−ｘ方向）に８画素刻みに移
動させることを１２回繰り返した後、図６（ｂ）に示す
ように、暫定マッチング用画像Ｇ１’を下縦方向（＋ｙ
方向）に８画素刻み移動させる。次いで、暫定マッチン
グ用画像Ｇ１’を右横方向（＋ｘ方向）に８画素刻みに
移動させることを１２回繰り返し、基準位置Ｏ0に戻っ
て暫定マッチング用画像Ｇ１’を左横方向（−ｘ方向）
に８画素刻みに移動させることを１２回繰り返す。この
暫定マッチング用画像Ｇ１’の８画素刻みの移動を＋ｙ
方向に１２回、−ｙ方向に１２回繰り返すと、総計６２
５回（＝２５×２５）の移動処理が行われる。
【００３９】そして、暫定マッチング用画像Ｇ１’の各
移動ステップ位置で、４画素刻みのそれぞれの各画素の
出力の差の積算値を計算し、その積算値を演算に用いた
画素数で割り算する。そして、その積算値の割り算値の
うち最小(最も相関度が高い)ものの方（小さい方）から
５点の移動位置を選んで、その5点の平均移動位置を算
出する。この平均移動位置をΔｘ0，Δｙ0とする。これ
によって、大まかな移動量（ズレ量）を把握できる。
移動ステップ位置で、４画素刻みのそれぞれの各画素の
出力の差の積算値を計算し、その積算値を演算に用いた
画素数で割り算する。そして、その積算値の割り算値の
うち最小(最も相関度が高い)ものの方（小さい方）から
５点の移動位置を選んで、その5点の平均移動位置を算
出する。この平均移動位置をΔｘ0，Δｙ0とする。これ
によって、大まかな移動量（ズレ量）を把握できる。
【００４０】暫定マッチング用画像Ｇ０’、Ｇ１’を使
用し、画像同士の対応関係をつけやすくなっているの
で、相関量の検出が容易となっている。
用し、画像同士の対応関係をつけやすくなっているの
で、相関量の検出が容易となっている。
【００４１】この平均移動位置Δｘ0,Δｙ0を基準位置
とする。 (２) 次いで、出力の差の積算値を求めるための画素
を、図７に拡大して示すように、３画素刻みとし、３画
素刻みで演算に用いる画素を選択する。その図７におい
て、ハッチングは相関の演算に用いる画素を示す。
とする。 (２) 次いで、出力の差の積算値を求めるための画素
を、図７に拡大して示すように、３画素刻みとし、３画
素刻みで演算に用いる画素を選択する。その図７におい
て、ハッチングは相関の演算に用いる画素を示す。
【００４２】画像の移動ステップ幅Ｗは、図８に示すよ
うに６画素刻みとする。移動のステップ数は項目(1)で
算出した平均移動位置（Δｘ，Δｙ）を基準位置として
プラスマイナス１４ステップとする。
うに６画素刻みとする。移動のステップ数は項目(1)で
算出した平均移動位置（Δｘ，Δｙ）を基準位置として
プラスマイナス１４ステップとする。
【００４３】項目(１)で述べたと同様に、この移動ステ
ップ幅Ｗでプラスマイナス４ステップ数(総計８１ステ
ップ)で移動させ、各移動位置で、３画素刻みの各の各
画素の出力の差の積算値を計算し、その積算値を演算に
用いた画素数で割り算する、その積算値の割り算値の最
小(最も相関度が高い)ものの方（小さい方）から５点の
移動位置を選んで、その5点の平均移動位置を算出す
る。この平均移動位置をΔｘ１，Δｙ２とする。これに
よって、項目（１）で粗く求めたズレ量よりもより精度
の良好なズレ量を求めることができる。この平均移動位
置を基準位置（Δｘ１，Δｙ２）とする。 (３) 出力の差の積算値を求めるための画素は、すべて
の画素として、画像の移動ステップ幅Ｗは、3画素刻み
とする。
ップ幅Ｗでプラスマイナス４ステップ数(総計８１ステ
ップ)で移動させ、各移動位置で、３画素刻みの各の各
画素の出力の差の積算値を計算し、その積算値を演算に
用いた画素数で割り算する、その積算値の割り算値の最
小(最も相関度が高い)ものの方（小さい方）から５点の
移動位置を選んで、その5点の平均移動位置を算出す
る。この平均移動位置をΔｘ１，Δｙ２とする。これに
よって、項目（１）で粗く求めたズレ量よりもより精度
の良好なズレ量を求めることができる。この平均移動位
置を基準位置（Δｘ１，Δｙ２）とする。 (３) 出力の差の積算値を求めるための画素は、すべて
の画素として、画像の移動ステップ幅Ｗは、3画素刻み
とする。
【００４４】この移動ステップ幅Ｗで、同様にプラスマ
イナス４ステップ数(総計８１ステップ)で移動させ、各
移動位置で、項目(２)で算出した平均移動位置（Δｘ
１，Δｙ１）を基準位置として全各画素の出力の差の積
算値を計算し、その積算値を演算に用いた画素数で割り
算する。その積算値の割り算値のうち最小(最も相関度
が高い)ものの方（小さい）方から５点の移動位置を選
んで、その5点の平均移動位置を算出する。この平均移
動位置をΔｘ３，Δｙ３とする。
イナス４ステップ数(総計８１ステップ)で移動させ、各
移動位置で、項目(２)で算出した平均移動位置（Δｘ
１，Δｙ１）を基準位置として全各画素の出力の差の積
算値を計算し、その積算値を演算に用いた画素数で割り
算する。その積算値の割り算値のうち最小(最も相関度
が高い)ものの方（小さい）方から５点の移動位置を選
んで、その5点の平均移動位置を算出する。この平均移
動位置をΔｘ３，Δｙ３とする。
【００４５】項目（１）ないし項目（３）の処理を、第
１画像移動量を算出する第１測定ステップといい、第１
測定ステップでは飛び飛びに所定ステップ数暫定マッチ
ング用画像を移動させて、相関量（画像のずれ量）が求
められる。 (４) 出力の差の積算値を求めるための画素は、全画素
として、画像の移動ステップ幅Ｗは１画素刻みとする。
１画像移動量を算出する第１測定ステップといい、第１
測定ステップでは飛び飛びに所定ステップ数暫定マッチ
ング用画像を移動させて、相関量（画像のずれ量）が求
められる。 (４) 出力の差の積算値を求めるための画素は、全画素
として、画像の移動ステップ幅Ｗは１画素刻みとする。
【００４６】この移動ステップ幅Ｗで同様にプラスマイ
ナス６ステップ数(総計１６９ステップ)で移動させ、各
移動位置で、すべての各画素の出力の差の積算値を計算
し、その積算値を演算に用いた画素数で割り算する。そ
の積算値のうち最小値(最も相関度が高い)を選んで移動
量を算出する。この移動位置をΔｘ４，Δｙ４とする。
ナス６ステップ数(総計１６９ステップ)で移動させ、各
移動位置で、すべての各画素の出力の差の積算値を計算
し、その積算値を演算に用いた画素数で割り算する。そ
の積算値のうち最小値(最も相関度が高い)を選んで移動
量を算出する。この移動位置をΔｘ４，Δｙ４とする。
【００４７】この項目（４）の処理を、第２画像移動量
を算出する第２測定ステップといい、第２測定ステップ
では、１画素毎に暫定マッチング用画像Ｇ１’を移動さ
せて相関量（画像のずれ量）が求められる。
を算出する第２測定ステップといい、第２測定ステップ
では、１画素毎に暫定マッチング用画像Ｇ１’を移動さ
せて相関量（画像のずれ量）が求められる。
【００４８】このように、漸次サーチ範囲を狭くしなが
ら画像の移動量（ズレ量）を順次算出し、最終的には、
精密な画像移動量ΔＸ=Δｘ１+Δｘ２+Δｘ３+Δｘ４，
ΔＹ=Δｙ１+Δｙ２+Δｙ３+Δｙ４が算出される。
ら画像の移動量（ズレ量）を順次算出し、最終的には、
精密な画像移動量ΔＸ=Δｘ１+Δｘ２+Δｘ３+Δｘ４，
ΔＹ=Δｙ１+Δｙ２+Δｙ３+Δｙ４が算出される。
【００４９】なお、Δｘ１〜Δｘ４、Δｙ１〜Δｙ４は
正負の符号を含むものとする。
正負の符号を含むものとする。
【００５０】以上の例では、暫定マッチング用画像の２
つの画像についてその移動量（ズレ量）を算出する方法
について述べたが、同様の手法により、下記記載する画
像に基づき移動量を求めることができる。
つの画像についてその移動量（ズレ量）を算出する方法
について述べたが、同様の手法により、下記記載する画
像に基づき移動量を求めることができる。
【００５１】最初に偶数フィールド画像の暫定マッチン
グ用画像の移動量（ズレ量）を求める。 最初の偶数フィールドでの暫定マッチング用画像Ｅ１
を基準として、次の偶数フィールドの暫定マッチング用
画像E2の移動量（ズレ量）を算出する。 偶数フィールドでの暫定マッチング用画像E2を基準と
して、次の偶数フィールドにおける暫定マッチング用画
像E3の移動量（ズレ量）を算出する。 以下、同様に移動量（ズレ量）の算出を行い、この
算出結果から、暫定マッチング用画像Ｅ１を基準とした
各Ｅ２、…、Eｎの画像の移動量（ズレ量）を算出す
る。
グ用画像の移動量（ズレ量）を求める。 最初の偶数フィールドでの暫定マッチング用画像Ｅ１
を基準として、次の偶数フィールドの暫定マッチング用
画像E2の移動量（ズレ量）を算出する。 偶数フィールドでの暫定マッチング用画像E2を基準と
して、次の偶数フィールドにおける暫定マッチング用画
像E3の移動量（ズレ量）を算出する。 以下、同様に移動量（ズレ量）の算出を行い、この
算出結果から、暫定マッチング用画像Ｅ１を基準とした
各Ｅ２、…、Eｎの画像の移動量（ズレ量）を算出す
る。
【００５２】次に、奇数フィールド画像の暫定マッチン
グ用画像の移動量（ズレ量）を算出する。 最初の偶数フィールドでの暫定画像E1を基準として、
最初の奇数フィールドの暫定マッヒ゜ンク゛画像Ｏ１の移動量
（ズレ量）を算出する。 次いで、奇数フィールドの暫定マッチング用画像Ｏ１
を基準として次の奇数フィールドでの暫定マッチング用
画像Ｏ2の移動量（ズレ量）を算出する。 次いで、奇数フィールドでの暫定マッチング用画像Ｏ
2を基準として、次の奇数フィールドでの暫定マッチン
グ用画像Ｏ3の移動量（ズレ量）を算出する。
グ用画像の移動量（ズレ量）を算出する。 最初の偶数フィールドでの暫定画像E1を基準として、
最初の奇数フィールドの暫定マッヒ゜ンク゛画像Ｏ１の移動量
（ズレ量）を算出する。 次いで、奇数フィールドの暫定マッチング用画像Ｏ１
を基準として次の奇数フィールドでの暫定マッチング用
画像Ｏ2の移動量（ズレ量）を算出する。 次いで、奇数フィールドでの暫定マッチング用画像Ｏ
2を基準として、次の奇数フィールドでの暫定マッチン
グ用画像Ｏ3の移動量（ズレ量）を算出する。
【００５３】以下、同様に奇数フィールドでの移動量
（ズレ量）の算出を行い、この算出結果から、最初の偶
数フィールド画像Ｅ０を基準とする各奇数フィールドで
の暫定マッチング用画像Ｏ１、…、Ｏｎの画像の移動量
（ズレ量）を算出する。
（ズレ量）の算出を行い、この算出結果から、最初の偶
数フィールド画像Ｅ０を基準とする各奇数フィールドで
の暫定マッチング用画像Ｏ１、…、Ｏｎの画像の移動量
（ズレ量）を算出する。
【００５４】この算出により、１つの基準画像、ここで
は、偶数フィールドでの暫定マッチング用画像Ｅ１を基
準として、奇数フィールド、偶数フィールド毎に記憶さ
れた全ての暫定マッチング用画像の移動量（ズレ量）が
算出される。
は、偶数フィールドでの暫定マッチング用画像Ｅ１を基
準として、奇数フィールド、偶数フィールド毎に記憶さ
れた全ての暫定マッチング用画像の移動量（ズレ量）が
算出される。
【００５５】この暫定マッチング用画像を用いて算出さ
れた移動量（ズレ量）に基づき、画像記録器１６に記憶
されている１つの画像を基準として、残りの各複数画像
の位置関係が座標変換され、像の移動量（ズレ量）が補
正された複数画像を得ることができる。すなわち、画像
記録器１６に記憶されている複数の画像の移動量（相対
的ズレ量）が補正される。
れた移動量（ズレ量）に基づき、画像記録器１６に記憶
されている１つの画像を基準として、残りの各複数画像
の位置関係が座標変換され、像の移動量（ズレ量）が補
正された複数画像を得ることができる。すなわち、画像
記録器１６に記憶されている複数の画像の移動量（相対
的ズレ量）が補正される。
【００５６】この像の移動量（ズレ量）が補正されてか
つ座標変換された複数の画像情報に基づき、血流速度
は、下記式で演算される。
つ座標変換された複数の画像情報に基づき、血流速度
は、下記式で演算される。
【００５７】すなわち、１秒間に例えば６０フレーム分
の速度で連続的に画像記録器１６に連続的に複数個のフ
レームが記憶されているものとする。記憶された画像情
報に基づき固体撮像素子５の（ｍ、ｎ）番目の画素につ
いて、ｋ回目の出力をＩk（ｍ、ｎ）とし、ｋ＋１回目
の出力をＩk+1（ｍ、ｎ）とするとき、下記数１式に基
づいて、各画素におけるスペックルの変動率を積算した
ＡＤ値（平均差異）を演算する。
の速度で連続的に画像記録器１６に連続的に複数個のフ
レームが記憶されているものとする。記憶された画像情
報に基づき固体撮像素子５の（ｍ、ｎ）番目の画素につ
いて、ｋ回目の出力をＩk（ｍ、ｎ）とし、ｋ＋１回目
の出力をＩk+1（ｍ、ｎ）とするとき、下記数１式に基
づいて、各画素におけるスペックルの変動率を積算した
ＡＤ値（平均差異）を演算する。
【００５８】
【数１】
【００５９】この演算式において、分母は各画素の出力
値をノーマライズ化するもので、これにより、血管部の
反射率の差に起因する要因は除去され、ＡＤ値はスペッ
クルの変動量、すなわち、血流速度の関数となる。この
ＡＤ値は血流の速度を示す。
値をノーマライズ化するもので、これにより、血管部の
反射率の差に起因する要因は除去され、ＡＤ値はスペッ
クルの変動量、すなわち、血流速度の関数となる。この
ＡＤ値は血流の速度を示す。
【００６０】このＡＤ値の演算は各画素毎に行い、この
演算で求めた各画素でのＡＤ値に基づいて、二次元的に
血流速度が例えばＴＶモニター１８に表示される。
演算で求めた各画素でのＡＤ値に基づいて、二次元的に
血流速度が例えばＴＶモニター１８に表示される。
【００６１】なお、この演算式については、特開平８−
１１２２６２号により公知である。
１１２２６２号により公知である。
【００６２】
【発明の効果】本発明は、以上説明したように構成した
ので、血流速度の演算に用いる複数枚の画像の位置ズレ
を補正し、高精度で鮮明な血流速度分布画像を得ること
ができるという効果を奏する。
ので、血流速度の演算に用いる複数枚の画像の位置ズレ
を補正し、高精度で鮮明な血流速度分布画像を得ること
ができるという効果を奏する。
【図１】 本発明の血流速度測定装置の要部構成を示す
図である。
図である。
【図２】 図１に示す画像記録部に記録された画像の一
例を示す説明図であって、（ａ）はある時刻において撮
影された眼底像を示し、（ｂ）は他の時刻において撮影
された眼底像を示し、（ａ）に示す画像と（ｂ）に示す
画像とでその位置関係がずれていることを示している。
例を示す説明図であって、（ａ）はある時刻において撮
影された眼底像を示し、（ｂ）は他の時刻において撮影
された眼底像を示し、（ａ）に示す画像と（ｂ）に示す
画像とでその位置関係がずれていることを示している。
【図３】 図２（ａ）に示す画像に基づき暫定マッチン
グ用画像Ｇ０’、Ｇ１’を作成するときに用いる画素の
説明図である。
グ用画像Ｇ０’、Ｇ１’を作成するときに用いる画素の
説明図である。
【図４】 図３に示す画素を用いて作成された暫定マッ
チング用画像Ｇ０’、Ｇ１’とを基準位置で重ね合わせ
た状態を示す図である。
チング用画像Ｇ０’、Ｇ１’とを基準位置で重ね合わせ
た状態を示す図である。
【図５】 暫定マッチング用画像に基づきずれ量を算出
するのに使用する画素の一例を示す説明図であって、ず
れ量を算出する画素を４画素毎にした例を示している。
するのに使用する画素の一例を示す説明図であって、ず
れ量を算出する画素を４画素毎にした例を示している。
【図６】 暫定マッチング用画像に基づきずれ量を算出
する場合の画像のずらし方の説明図で、（ａ）は暫定マ
ッチング用画像Ｇ１’を暫定マッチング用画像Ｇ０’に
対して右に８画素毎にずらしている状態を示し、（ｂ）
は暫定マッチング用画像Ｇ１’を暫定マッチング用画像
Ｇ０’に対して下に８画素ずらした後、右に８画素毎ず
らしている状態を示している。
する場合の画像のずらし方の説明図で、（ａ）は暫定マ
ッチング用画像Ｇ１’を暫定マッチング用画像Ｇ０’に
対して右に８画素毎にずらしている状態を示し、（ｂ）
は暫定マッチング用画像Ｇ１’を暫定マッチング用画像
Ｇ０’に対して下に８画素ずらした後、右に８画素毎ず
らしている状態を示している。
【図７】 暫定マッチング用画像に基づきずれ量を算出
するのに使用する画素の一例を示す説明図であって、ず
れ量を算出する画素を３画素毎にした例を示している。
するのに使用する画素の一例を示す説明図であって、ず
れ量を算出する画素を３画素毎にした例を示している。
【図８】 暫定マッチング用画像に基づきずれ量を算出
する場合の画像のずらし方の説明図で、先に演算で求め
た平均ずれ量を基準位置として、暫定マッチング用画像
Ｇ１’を暫定マッチング用画像Ｇ０’に対して右に６画
素毎にずらしている状態を示す説明図である。
する場合の画像のずらし方の説明図で、先に演算で求め
た平均ずれ量を基準位置として、暫定マッチング用画像
Ｇ１’を暫定マッチング用画像Ｇ０’に対して右に６画
素毎にずらしている状態を示す説明図である。
１…レーザー照射系
２…受光系
５…ＣＣＤ（受光部）
１２…信号選択手段（画像取り込み部）
１３…第１画像記録手段（画像取り込み部）
１４…第２画像記録手段（画像取り込み部）
１５…画像合成手段（画像取り込み部）
１６…画像記録器（画像記憶部）
１７…演算部
─────────────────────────────────────────────────────
フロントページの続き
(51)Int.Cl.7 識別記号 ＦＩ テーマコート゛(参考）
Ｇ０１Ｐ 5/20 Ａ６１Ｂ 5/02 ３４０Ｄ
Ｇ０６Ｔ 1/00 ２９０ 3/12 Ｅ
7/20 3/10 Ｒ
(72)発明者 藤居 仁
福岡県宗像市日の里２丁目26番８
(72)発明者 高良 一弘
福岡県福岡市早良区南庄４丁目12−201号
Ｆターム(参考） 2G059 AA05 BB12 BB13 CC16 EE02
FF01 GG01 JJ11 KK04 MM01
MM03 MM05 MM09 MM10 PP04
4C017 AA11 AB07 AC28 BC01 BC11
CC03 FF30
5B057 AA07 BA02 CA12 CA16 DB02
DC02 DC32
5L096 BA06 BA13 CA02 GA08 HA04
Claims (7)
- 【請求項１】 生体組織の血球にレーザ光を照射するレ
ーザ光照射系と、前記生体組織からの反射光を検出する
多数の画素からなる受光部を有する受光系と、前記受光
部からの信号に基づき所定時間間隔で連続的に複数の画
像を取り込む画像取り込み部と、前記複数の画像を記憶
する画像記憶部と、該記憶された複数画像の対応する各
画素の出力信号の時間的変化から生体組織内の血流速度
を演算する演算部とからなる血流速度測定装置におい
て、 前記演算部は前記複数画像の相互の位置ずれ情報を検出
する位置ずれ検出部を有し、前記位置ずれ情報に基づき
記憶された複数画像の位置ずれを補正することを特徴と
する血流速度測定装置。 - 【請求項２】 前記位置ずれ検出部は、各画素出力を該
画素の近傍の複数の画素の出力に基づき演算した値とし
た暫定マッチング用画像を作成する請求項１に記載の血
流速度測定装置。 - 【請求項３】 前記位置ずれ検出部は、２つの前記暫定
マッチング用画像のマッチングをとることにより各画像
の相互の位置ずれを検出する請求項２に記載の血流速度
測定装置。 - 【請求項４】 前記画像のマッチングは、一方の画像を
他方の画像上でとびとびの所定画素毎に所定ステップ数
移動させて各移動ステップで求めた相関値から第１画像
移動量を算出する第１測定ステップと、前記第１測定ス
テップで算出された第１画像移動量を基準として１画素
毎に所定ステップ数移動させて各移動ステップで求めた
相関値から第２画像移動量を算出する第２測定ステップ
とからなり、前記第１画像移動量と第２画像移動量とか
ら全体の画像移動量を算出する請求項３に記載の血流速
度測定装置。 - 【請求項５】 前記第１測定ステップは、移動画素数を
順次少なくする複数の測定ステップからなる請求項４に
記載の血流速度測定装置。 - 【請求項６】 前記第１測定ステップは、所定間隔離間
した各画素間で相関をとる請求項４に記載の血流速度測
定装置。 - 【請求項７】 前記受光部は固体撮像素子であり、前記
画像取り込み部は、奇数走査線に基づく奇数フィールド
画像と偶数走査線の信号に基づく偶数フィールド画像と
を１組として取り込み、前記画像記憶部は前記奇数フィ
ールド画像と偶数フィールド画像との組み合わせからな
る合成フレーム画像として記憶する請求項１に記載の血
流速度測定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001365419A JP2003164431A (ja) | 2001-11-30 | 2001-11-30 | 血流速度測定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001365419A JP2003164431A (ja) | 2001-11-30 | 2001-11-30 | 血流速度測定装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2003164431A true JP2003164431A (ja) | 2003-06-10 |
Family
ID=19175445
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001365419A Pending JP2003164431A (ja) | 2001-11-30 | 2001-11-30 | 血流速度測定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2003164431A (ja) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008017961A (ja) * | 2006-07-12 | 2008-01-31 | Image One Co Ltd | 血流速度の測定方法及び装置 |
| US7597443B2 (en) | 2004-10-28 | 2009-10-06 | Kyushu Institute Of Technology | Wide viewing angle ocular fundus blood flow imaging device |
| WO2010084927A1 (ja) * | 2009-01-22 | 2010-07-29 | 日本電気株式会社 | 画像処理装置、生体認証装置、画像処理方法及び記録媒体 |
| US7817256B2 (en) | 2006-02-22 | 2010-10-19 | Kyushu Institute Of Technology | Personal authentication method and personal authentication device utilizing finger-tip blood flow measurement by laser light |
| WO2010131550A1 (ja) | 2009-05-13 | 2010-11-18 | 国立大学法人九州工業大学 | 血流画像診断装置 |
| US7844083B2 (en) | 2004-06-18 | 2010-11-30 | Kyushu Institute Of Technology | Method for acquiring personal identification data, personal identification method, apparatus for acquiring personal identification data, and personal identification apparatus |
| US8285003B2 (en) | 2006-06-07 | 2012-10-09 | Kyushu Institute Of Technology | Personal authentication method and personal authentication device utilizing ocular fundus blood flow measurement by laser light |
| US8462995B2 (en) | 2008-04-06 | 2013-06-11 | Kyushu Institute Of Technology | Method of authenticating individual and apparatus of authenticating individual used for the same |
| US8494228B2 (en) | 2008-04-03 | 2013-07-23 | Kyushu Institute Of Technology | Personal authentication method using subcutaneous bloodstream measurement and personal authentication device |
| WO2014175154A1 (ja) | 2013-04-23 | 2014-10-30 | ソフトケア有限会社 | 血流画像診断装置及び診断方法 |
| WO2017175470A1 (ja) * | 2016-04-05 | 2017-10-12 | ソニー株式会社 | スペックル測定装置およびスペックル測定方法 |
| WO2018003139A1 (ja) | 2016-06-28 | 2018-01-04 | ソフトケア有限会社 | 血流動態画像化診断装置及び診断方法 |
| JP2020099591A (ja) * | 2018-12-25 | 2020-07-02 | 花王株式会社 | 皮膚毛細血管の撮影及び解析方法 |
-
2001
- 2001-11-30 JP JP2001365419A patent/JP2003164431A/ja active Pending
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7844083B2 (en) | 2004-06-18 | 2010-11-30 | Kyushu Institute Of Technology | Method for acquiring personal identification data, personal identification method, apparatus for acquiring personal identification data, and personal identification apparatus |
| US7597443B2 (en) | 2004-10-28 | 2009-10-06 | Kyushu Institute Of Technology | Wide viewing angle ocular fundus blood flow imaging device |
| US7817256B2 (en) | 2006-02-22 | 2010-10-19 | Kyushu Institute Of Technology | Personal authentication method and personal authentication device utilizing finger-tip blood flow measurement by laser light |
| US8285003B2 (en) | 2006-06-07 | 2012-10-09 | Kyushu Institute Of Technology | Personal authentication method and personal authentication device utilizing ocular fundus blood flow measurement by laser light |
| JP2008017961A (ja) * | 2006-07-12 | 2008-01-31 | Image One Co Ltd | 血流速度の測定方法及び装置 |
| US8494228B2 (en) | 2008-04-03 | 2013-07-23 | Kyushu Institute Of Technology | Personal authentication method using subcutaneous bloodstream measurement and personal authentication device |
| US8462995B2 (en) | 2008-04-06 | 2013-06-11 | Kyushu Institute Of Technology | Method of authenticating individual and apparatus of authenticating individual used for the same |
| US9544146B2 (en) | 2009-01-22 | 2017-01-10 | Nec Corporation | Image processing apparatus, biometric authentication apparatus, image processing method and recording medium |
| WO2010084927A1 (ja) * | 2009-01-22 | 2010-07-29 | 日本電気株式会社 | 画像処理装置、生体認証装置、画像処理方法及び記録媒体 |
| JP2010170348A (ja) * | 2009-01-22 | 2010-08-05 | Nec Corp | 画像処理装置、生体認証装置、画像処理方法及びプログラム |
| WO2010131550A1 (ja) | 2009-05-13 | 2010-11-18 | 国立大学法人九州工業大学 | 血流画像診断装置 |
| US9028421B2 (en) | 2009-05-13 | 2015-05-12 | Kyushu Institute Of Technology | Blood flow image diagnosing device |
| WO2014175154A1 (ja) | 2013-04-23 | 2014-10-30 | ソフトケア有限会社 | 血流画像診断装置及び診断方法 |
| US10098592B2 (en) | 2013-04-23 | 2018-10-16 | Softcare Co., Ltd. | Blood flow image diagnosing device and method |
| WO2017175470A1 (ja) * | 2016-04-05 | 2017-10-12 | ソニー株式会社 | スペックル測定装置およびスペックル測定方法 |
| JPWO2017175470A1 (ja) * | 2016-04-05 | 2019-02-14 | ソニー株式会社 | スペックル測定装置およびスペックル測定方法 |
| EP3440997A4 (en) * | 2016-04-05 | 2019-04-17 | Sony Corporation | SPECKLEMESSVORRICHTUNG AND SPECKLEMESS PROCEDURE |
| US10801831B2 (en) | 2016-04-05 | 2020-10-13 | Sony Corporation | Speckle measurement apparatus and speckle measurement method |
| WO2018003139A1 (ja) | 2016-06-28 | 2018-01-04 | ソフトケア有限会社 | 血流動態画像化診断装置及び診断方法 |
| US11330995B2 (en) | 2016-06-28 | 2022-05-17 | Softcare Co., Ltd. | Apparatus and method for imaging and analyzing hemodynamics |
| JP2020099591A (ja) * | 2018-12-25 | 2020-07-02 | 花王株式会社 | 皮膚毛細血管の撮影及び解析方法 |
| JP7164090B2 (ja) | 2018-12-25 | 2022-11-01 | 花王株式会社 | 皮膚毛細血管の撮影及び解析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lindner et al. | Lateral and depth calibration of PMD-distance sensors | |
| US6445814B2 (en) | Three-dimensional information processing apparatus and method | |
| US5832106A (en) | Method for camera calibration of range imaging system by use of neural network | |
| JP2003164431A (ja) | 血流速度測定装置 | |
| US6950550B1 (en) | Tracing technique and recording media of object motion | |
| US20060098846A1 (en) | Movement analysis apparatus | |
| CN113543718B (zh) | 用于确定超声探头的包括前后方向性的运动的装置和方法 | |
| US9183634B2 (en) | Image processing apparatus and image processing method | |
| EP3417919B1 (en) | Transformation matrix deriving device, position estimation apparatus, transformation matrix deriving method, and position estimation method | |
| JP3855024B2 (ja) | 血流速度測定装置 | |
| JP3637226B2 (ja) | 動き検出方法、動き検出装置及び記録媒体 | |
| CN101909165A (zh) | 一种基于混合测度的视频数据宽景成像方法 | |
| D'Apuzzo | Surface measurement and tracking of human body parts from multi station video sequences | |
| US11972544B2 (en) | Method and apparatus for optical coherence tomography angiography | |
| JP7297705B2 (ja) | 処理装置、処理方法、学習装置およびプログラム | |
| JP3534551B2 (ja) | 動き検出装置 | |
| JP2001338280A (ja) | ３次元空間情報入力装置 | |
| JP2005141655A (ja) | ３次元モデリング装置及び３次元モデリング方法 | |
| JP7281660B2 (ja) | 変位計測装置及び変位計測方法 | |
| JPH0443204B2 (ja) | ||
| JP2007235769A (ja) | 動きベクトル検出方法及び装置 | |
| JP2961140B2 (ja) | 画像処理方法 | |
| JPS63240832A (ja) | 肺野部輪郭線を求める画像計測方法 | |
| CN117647263B (en) | Nonlinear optimization-based single photon camera visual inertial odometer method and system | |
| JP3077773B2 (ja) | 三次元画像表示装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20051019 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20051025 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20060314 |