JP2003180641A - Blood flow velocity measuring apparatus - Google Patents

Blood flow velocity measuring apparatus

Info

Publication number
JP2003180641A
JP2003180641A JP2001379858A JP2001379858A JP2003180641A JP 2003180641 A JP2003180641 A JP 2003180641A JP 2001379858 A JP2001379858 A JP 2001379858A JP 2001379858 A JP2001379858 A JP 2001379858A JP 2003180641 A JP2003180641 A JP 2003180641A
Authority
JP
Japan
Prior art keywords
pixel
pixels
image
images
blood flow
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.)
Granted
Application number
JP2001379858A
Other languages
Japanese (ja)
Other versions
JP3855024B2 (en
Inventor
Hitoshi Fujii
仁 藤居
Naoki Konishi
直樹 小西
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Topcon Corp
Original Assignee
Topcon Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Topcon Corp filed Critical Topcon Corp
Priority to JP2001379858A priority Critical patent/JP3855024B2/en
Publication of JP2003180641A publication Critical patent/JP2003180641A/en
Application granted granted Critical
Publication of JP3855024B2 publication Critical patent/JP3855024B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Abstract

<P>PROBLEM TO BE SOLVED: To provide a blood flow velocity measuring apparatus by which the same level measurement accuracy can be obtained as conventional apparatus using at least three images without requiring multiple images. <P>SOLUTION: The blood flow velocity measuring apparatus has a storage section storing in turn multiple images formed on a light-receiving part based on the reflecting signals from a body tissue by applying laser beam to the blood cells of the body tissue, and calculates the distribution condition of the blood cell velocity based on at least three images of the stored images. The calculation section detects the blood flow distribution condition in the body tissue by selecting a first pixel group comprising a plurality of pixels selected from the pixels of a predetermined address position and the pixels present adjacent to the above predetermined address position within one image, selecting a second and third pixel groups present respectively at the same address as the above first pixel group within other images, selecting predetermined two combined pixels from the above the first to third pixel groups, calculating the average of the total of the differences of the outputs of the selected two pixels or the inverse number, and carrying out the difference calculation on each predetermined address position pixel. <P>COPYRIGHT: (C)2003,JPO

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【発明の属する技術分野】本発明は、血球を有する生体
組織にレーザ光を照射し、その血球から反射されたスペ
ックル信号に基づき血流速度を測定するための血流速度
測定装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a blood flow velocity measuring device for irradiating a living tissue having blood cells with laser light and measuring a blood flow velocity based on a speckle signal reflected from the blood cells.

【0002】[0002]

【従来の技術】従来、被検眼の眼底等の生体組織の血球
にレーザ光を照射して、その血球からの反射光により形
成された画像を固体撮像装置(CCD)等のイメージセ
ンサー上に導き、この画像を連続的に所定時間間隔で多
数枚取り込んで記憶し、その記憶された多数の画像の中
から所定枚数の画像を選択し、各画像の各画素における
出力の時間変動量を積算した値、或いは、その逆数の値
を算出し、この値から血球の速度(血流速度)を算出す
る血流速度測定装置が知られている。
2. Description of the Related Art Conventionally, a blood cell of a living tissue such as a fundus of an eye to be examined is irradiated with laser light and an image formed by reflected light from the blood cell is guided onto an image sensor such as a solid-state imaging device (CCD). , This image is continuously captured and stored at a predetermined time interval, a predetermined number of images are selected from the stored many images, and the time variation amount of output at each pixel of each image is integrated. There is known a blood flow velocity measuring device that calculates a value or its reciprocal value and calculates the blood cell velocity (blood flow velocity) from this value.

【0003】この種の血流速度測定装置では、各画素の
出力変動量が血球の移動速度に対応するので、この算出
された各画素の出力変動量値、或いは、その逆数の値に
基づき、生体組織での血流速度分布を二次元画像として
モニター画面上に表示している。
In this type of blood flow velocity measuring device, since the output fluctuation amount of each pixel corresponds to the blood cell moving speed, based on the calculated output fluctuation amount value of each pixel or its reciprocal value, The blood flow velocity distribution in living tissue is displayed as a two-dimensional image on the monitor screen.

【0004】その各画素で演算される値、すなわち、血
流速度を示す演算値をSBR値とすると、例えば、下記
式が使用されている。
Assuming that the value calculated by each pixel, that is, the calculated value indicating the blood flow velocity is the SBR value, the following formula is used, for example.

【0005】[0005]

【数1】 [Equation 1]

【0006】この数1式において、SBRn,mは、図1
に示すように、例えば、1秒間にN枚撮影された一連の
時系列的な画像のうちのn,m番目の画素でのSBR値を
示し、In,m,tはn,m番目の画素における時間ti(i
=1、2、…、N)での画素出力を示すもので、〈 〉
は平均値を示すものである。
In this equation 1, SBR n, m is given by
As shown in, for example, the SBR value at the n-th and m-th pixels in a series of time-series images taken in one second is shown, and I n, m, t is the n-th and m-th pixels. The time t i (i
= 1, 2, ..., N), and <>
Indicates an average value.

【0007】その数1式の分子は、時間により変化する
n,m番目の画素の出力の所定時間内での平均値の二乗
を示す。その数1式の分母は、n,m番目の画素の各時
間での画素出力と、n,m番目の画素の所定時間内での
平均出力値との差分の絶対値の平均値の二乗を示し、こ
の分母がn,m番目の画素位置での出力変動の分散値に
対応した値を意味するものである。
The numerator of the equation (1) shows the square of the average value of the outputs of the nth and mth pixels which change with time within a predetermined time. The denominator of Equation 1 is the square of the average value of the absolute values of the differences between the pixel output of the nth and mth pixel at each time and the average output value of the nth and mth pixel within a predetermined time. This denominator means a value corresponding to the variance value of the output fluctuation at the n-th and m-th pixel positions.

【0008】このようにして、SBR値をn=1、m=
1番地の画素Xからn=n、m=m番地の画素Xまでに
渡って求め、図2に示すように、画面1にSBR値に対
応して色づけを行って血流速度を表示する。
In this way, the SBR values are n = 1 and m =
It is obtained from the pixel X at the address 1 to the pixel X at the address n = n, m = m, and the blood flow velocity is displayed on the screen 1 by coloring in accordance with the SBR value, as shown in FIG.

【0009】すなわち、スペックル変動が速いほどシャ
ッター時間内で積分され、出力変動のばらつきが減少す
る。従って、数1式のSBR値の分母が減少し、SBR
値が増大する。このSBR値の大小に対応させて色づけ
を行い、血流速度の変化を画面1に表示させている。
That is, the faster the speckle fluctuation is, the more the speckle fluctuation is integrated within the shutter time, and the fluctuation in the output fluctuation is reduced. Therefore, the denominator of the SBR value in Equation 1 decreases, and
The value increases. Coloring is performed according to the magnitude of this SBR value, and the change in blood flow velocity is displayed on the screen 1.

【0010】[0010]

【発明が解決しようとする課題】このように、この種の
血流速度測定装置では、血流速度の分布の時間的変化を
見るため、例えば、1秒間に500フレームの高速シャ
ッター間隔で画像を取り込み、その画像の中から例えば
N=64枚の一連の時系列的画像を取り出し、その64
枚の画像に基づいて、数1式により、血流速度に対応す
るSBR値を算出している。
As described above, in this type of blood flow velocity measuring apparatus, in order to see the temporal change in the blood flow velocity distribution, for example, an image is taken at a high-speed shutter interval of 500 frames per second. Captured, for example, a series of N = 64 time-series images is extracted from the image, and 64
The SBR value corresponding to the blood flow velocity is calculated by the equation 1 based on the images.

【0011】この数1式の演算方法によれば、高精度な
血流速度結果を得るためには少なくとも30枚程度の画
像が必要となり、高精度の血流速度を測定でき、かつ、
その血流速度の時間的変化の測定分解能を上げるには、
シャッター間隔が短い高速の連続画像取り込み装置が必
要となっている。
According to the calculation method of the equation (1), at least about 30 images are required to obtain a highly accurate blood flow velocity result, and the highly accurate blood flow velocity can be measured, and
To improve the measurement resolution of the temporal change in blood flow velocity,
There is a need for a high speed continuous image capture device with short shutter intervals.

【0012】本発明は、この従来の問題点を解決するこ
とを目的とし、血流速度演算のための画像としては多数
枚の画像を必要とせず、少なくとも3枚の画像を用いて
従来の血流速度測定装置と同等レベルの測定精度を得る
ことができる血流速度測定装置を提供することにある。
An object of the present invention is to solve this conventional problem and does not require a large number of images for blood flow velocity calculation, and at least three images can be used for the conventional blood flow calculation. It is an object of the present invention to provide a blood flow velocity measuring device that can obtain the same level of measurement accuracy as the flow velocity measuring device.

【0013】[0013]

【課題を解決するための手段】請求項1に記載の血流速
度測定装置は、生体組織の血球にレーザ光を照射するレ
ーザ光照射部と、前記生体組織からの反射光を検出する
多数の受光素子からなる受光部と、前記受光部からの信
号に基づき受光部上に形成された画像を順次複数枚記録
する記憶部と、前記記憶された画像の中の少なくとも3
枚の画像に基づき血球の速度の分布状態を演算する演算
部とからなり、前記演算部は、一つの画像の中の一つの
所定番地位置の画素及び該所定番地近傍に位置する画素
の中から選択された複数の画素からなる第1の画素群
と、その他の画像の中で前記第1の画素群とそれぞれ同
じ番地に位置する第2及び第3の画素群とを選択し、前
記第1ないし第3画素群から所定の2つの組み合わせ画
素を選択し、選択された2つの画素の出力の差分の総和
の平均値もしくはその逆数を演算し、その差分演算を各
所定番地の画素毎に行うことにより、生体組織での血流
分布状態を検出することを特徴とする。
A blood flow velocity measuring apparatus according to a first aspect of the present invention includes a laser light irradiating section for irradiating blood cells of a living tissue with a laser light, and a large number of detecting a reflected light from the living tissue. A light-receiving unit formed of a light-receiving element, a storage unit that sequentially records a plurality of images formed on the light-receiving unit based on a signal from the light-receiving unit, and at least three of the stored images.
And a calculation unit that calculates the distribution state of blood cell velocities based on a single image, wherein the calculation unit selects from among pixels at one predetermined address position and pixels located near the predetermined address in one image. The first pixel group including a plurality of selected pixels and the second and third pixel groups located at the same address as the first pixel group in other images are selected, and the first pixel group is selected. Or, a predetermined two combination pixels are selected from the third pixel group, the average value of the sums of the differences between the outputs of the two selected pixels or the reciprocal thereof is calculated, and the difference calculation is performed for each pixel at each predetermined address. Thus, the state of blood flow distribution in the living tissue is detected.

【0014】請求項2に記載の血流速度測定装置は、前
記受光部からの信号に基づき奇数走査線に基づく奇数フ
ィールド画像と偶数走査線の信号に基づく偶数フィール
ド画像とを1組として順次取り込み、記憶部に取り込ん
だ前記偶数フィールド画像と前記奇数フィールド画像と
を合成して合成フレーム画像として記憶し、前記演算部
は前記記憶された画像の中の少なくとも3枚の合成フレ
ーム画像に基づき血球の速度の分布状態を演算すること
を特徴とする。
A blood flow velocity measuring apparatus according to a second aspect of the invention sequentially captures an odd field image based on an odd scanning line based on a signal from the light receiving unit and an even field image based on a signal from an even scanning line as one set. , The even-numbered field image and the odd-numbered field image captured in the storage unit are combined and stored as a combined frame image, and the calculation unit calculates a blood cell based on at least three combined frame images in the stored image. It is characterized in that the velocity distribution state is calculated.

【0015】請求項3に記載の血流速度測定装置は、前
記演算部で演算に使用する画像は、受光素子上での像の
ずれ量を補正するために座標変換を行った画像であるこ
とを特徴とする。
In the blood flow velocity measuring apparatus according to the third aspect, the image used for the calculation in the calculation unit is an image which has been subjected to coordinate conversion in order to correct the shift amount of the image on the light receiving element. Is characterized by.

【0016】請求項4に記載の血流速度測定装置は、第
1ないし第3の画素群はさらに、それぞれ画素番地位置
が異なる2つの画素群に分けられ、該2つの画素群のそ
れぞれの画素群の中で所定の2つの組み合わせ画素を選
択することを特徴とする。
In the blood flow velocity measuring apparatus according to the fourth aspect, the first to third pixel groups are further divided into two pixel groups having different pixel address positions, and each pixel of the two pixel groups is divided. It is characterized in that two predetermined combination pixels are selected from the group.

【0017】[0017]

【発明の実施の形態】図19は血流速度測定装置の光学
系の概要を示すもので、1はレーザー光照射系、2は受
光系、Eは被検眼である。レーザー光照射系1のレーザ
ー光は、例えば、ハーフミラー3を介して被検眼Eの生
体組織としての例えば眼底Erに照射される。
DESCRIPTION OF THE PREFERRED EMBODIMENTS FIG. 19 shows an outline of an optical system of a blood flow velocity measuring device, 1 is a laser beam irradiation system, 2 is a light receiving system, and E is an eye to be examined. The laser light of the laser light irradiation system 1 is applied to the fundus Er, which is a living tissue of the eye E to be inspected, through the half mirror 3, for example.

【0018】受光系2は受光レンズ4、受光部としての
CCD(固体撮像素子)5、信号増幅回路6を有する。
眼底Erからのレーザー反射光は受光レンズ4によりC
CD5に生体組織像として結像される。CCD5はその
受光面上に多数の画素を有し、受光レンズ4により結像
された生体組織像を電気信号に変換し、フレーム蓄積方
式で信号電荷を読み出して映像信号として出力する。そ
の映像信号は、信号増幅回路6により増幅され、信号増
幅回路6により増幅された映像信号は利得制御等を行う
アナログ信号処理手段7に出力され、A/D変換器8によ
りデジタル信号に変換される。
The light receiving system 2 has a light receiving lens 4, a CCD (solid-state image sensor) 5 as a light receiving portion, and a signal amplification circuit 6.
The laser reflected light from the fundus Er is C by the light receiving lens 4.
An image of a living tissue is formed on the CD5. The CCD 5 has a large number of pixels on its light-receiving surface, converts the biological tissue image formed by the light-receiving lens 4 into an electrical signal, reads out the signal charges by the frame accumulation method, and outputs it as a video signal. The video signal is amplified by the signal amplification circuit 6, and the video signal amplified by the signal amplification circuit 6 is output to the analog signal processing means 7 for performing gain control and the like, and converted into a digital signal by the A / D converter 8. It

【0019】9はタイミングパルス発生器であり、10
は電子シャッター制御手段、11は固体撮像素子駆動手
段であり、タイミングパルス発生器9は電子シャッター
制御手段10と信号選択手段12とにタイミングパルス
を出力する。固体撮像素子駆動手段11はタイミングパ
ルスに基づき駆動される。
Reference numeral 9 is a timing pulse generator, and 10
Is an electronic shutter control means, 11 is a solid-state image sensor driving means, and the timing pulse generator 9 outputs a timing pulse to the electronic shutter control means 10 and the signal selection means 12. The solid-state image pickup device driving means 11 is driven based on the timing pulse.

【0020】信号選択手段12にはA/D変換器8によ
りA/D変換された映像信号としてのデジタル信号が入
力される。信号選択手段12は、タイミングパスル発生
器9からのタイミングパルスに基づいて、A/D変換器8
からの映像信号のうち、偶数フィールドの画像を第1画
像記録手段13へ送信し、奇数フィールドの画像を第2
画像記録手段14へ送信し、第1画像記録手段13には
偶数フィールドの画像が記憶され、第2画像記録手段1
4には奇数フィールドの画像が記録される。その信号選
択手段12、第1画像記録手段13、第2画像記録手段
14は、所定時間間隔で複数枚の画像を取り込む画像取
り込み部の一部として機能する。
A digital signal as a video signal A / D converted by the A / D converter 8 is input to the signal selection means 12. The signal selecting means 12 is based on the timing pulse from the timing pulse generator 9 and is based on the A / D converter 8
Of the video signal from the first image recording means 13 is transmitted to the first image recording means 13, and the image of the odd field is transmitted to the second image recording means 13.
The image of the even-numbered field is stored in the first image recording means 13 and transmitted to the image recording means 14, and the second image recording means 1
An image of an odd field is recorded in 4. The signal selecting unit 12, the first image recording unit 13, and the second image recording unit 14 function as a part of an image capturing unit that captures a plurality of images at predetermined time intervals.

【0021】その第1画像記録手段13、第2画像記録
手段14に記録された奇数フィールドの画像データと偶
数フィールドの画像データとは、画像合成手段15によ
り合成され、1/30秒間隔で撮影した1フレームの画
像データとされ、その1フレームの画像データは画像記
憶部としての画像記録器16に記憶される。
The image data of the odd field and the image data of the even field recorded in the first image recording means 13 and the second image recording means 14 are combined by the image combining means 15 and photographed at 1/30 second intervals. The image data of one frame is stored in the image recorder 16 as an image storage unit.

【0022】この画像記録器16に記憶された画像信号
は演算処理部17に入力され、演算処理部17は、後述
する演算処理を行う。なお、18はTVモニターであ
る。
The image signal stored in the image recording unit 16 is input to the arithmetic processing unit 17, and the arithmetic processing unit 17 performs arithmetic processing described later. In addition, 18 is a TV monitor.

【0023】演算処理部17は、画像記録器16に記憶
された画像から、例えば1/30秒間隔で撮影された3
枚の連続した画像を取り出し、その画像信号に基づき演
算を行うものである。
The arithmetic processing unit 17 takes three images, for example, taken at 1/30 second intervals from the image stored in the image recording device 16.
This is to take out a continuous image of one sheet and perform an operation based on the image signal.

【0024】演算処理部17は、眼底血流の測定演算を
行う前に、まず、3枚の画像の像のずれ量を補正する。
この補正は、3枚の画像の取り込みの間に、被検眼の固
視微動等の影響で像の位置ずれを生ずるおそれがある場
合に行う。その補正は、まず、3枚の画像のそれぞれの
画像間で、いわゆるテンプレート法で画像の相関を検出
し、その相関により、各画像の相互の位置ずれ量を計算
し、このずれ量に基づき、像のずれ量を考慮して座標変
換を行い、3枚の画像間のずれを補正した画像情報を作
成する(この画像間の像の位置ズレ量の演算について
は、例えば、特願2001−365419号参照)。な
お、この座標変換は、像のずれが無視できる場合には、
必要としない。
The calculation processing section 17 first corrects the image shift amount of the three images before performing the calculation calculation of the fundus blood flow.
This correction is performed when there is a possibility that an image position shift will occur due to the influence of the involuntary eye movement of the subject's eye during the capture of the three images. The correction is performed by first detecting the correlation of images among the three images by the so-called template method, calculating the mutual positional deviation amount of each image based on the correlation, and based on this deviation amount, Coordinate conversion is performed in consideration of the image shift amount, and image information in which the shift between the three images is corrected is created (for the calculation of the image positional shift amount between these images, for example, Japanese Patent Application No. 2001-365419). No.). In addition, this coordinate transformation, if the deviation of the image can be ignored,
do not need.

【0025】以下、この像の位置ずれを補正した3枚の
画像に基づく演算方法について述べる。
The calculation method based on the three images in which the positional deviation of the image has been corrected will be described below.

【0026】図3は画像記録器16に記憶されている時
系列的画像としての3枚の画像G1、G2、G3を示し
ている。その図3において、4×4の升目は着目すべき
画素群を含む画素群を示し、A1、A2、A3は番地
n、mの画素Xの出力値を示している。
FIG. 3 shows three images G1, G2, G3 as time series images stored in the image recorder 16. In FIG. 3, a 4 × 4 square represents a pixel group including a pixel group of interest, and A1, A2, and A3 represent output values of the pixel X at addresses n and m.

【0027】ここで、第1行目を奇数フィールドの走査
線Oに対応する画素群とすると、第2行目は偶数フィー
ルドの走査線eに対応する画素群であり、第3行目は奇
数フィールドの走査線Oに対応する画素群であり、第4
行目は偶数フィールドの走査線eに対応する画素群であ
る。
Here, assuming that the first row is a pixel group corresponding to the scanning lines O of the odd field, the second row is a pixel group corresponding to the scanning lines e of the even field, and the third row is an odd number. A pixel group corresponding to the scanning line O of the field,
The row is a pixel group corresponding to the scanning lines e of even fields.

【0028】従来の数1式に基づくSBR値を演算する
ものとすると、数1式の分母の値は、<│A1−<A>
│+│A2−<A>│+│A3−<A>│>2となり、
分母の差分の個数は3である。<A>は平均値であり、
<A>=(A1+A2+A3)/3である。従って、分
子の値は、<A>2である。
Assuming that the SBR value based on the conventional equation 1 is calculated, the value of the denominator of the equation 1 is <│A1- <A>.
│ + │A2- <A> │ + │A3- <A>│> 2 ,
The number of denominator differences is three. <A> is the average value,
<A> = (A1 + A2 + A3) / 3. Therefore, the numerator value is <A> 2 .

【0029】ゆえに、この場合のSBR値は、 SBR=<A>2/<│A1−<A>│+│A2−<A
>│+│A3−<A>│>2 となる。
Therefore, the SBR value in this case is SBR = <A> 2 / <│A1- <A> │ + │A2- <A
> │ + │A3- <A>│> 2 .

【0030】このように図3に示す3枚の画像のみを用
いて演算を行うものとすると、分母のサンプル個数が3
となり、そのサンプル個数が少ないため、3枚の画像の
みを用いて血流速度の演算を行うとするとは、到底安定
した高精度の測定を期待できないことになる。
If the calculation is performed using only the three images shown in FIG. 3, the denominator sample number is 3
Since the number of samples is small, if the blood flow velocity is calculated using only three images, it is impossible to expect stable and highly accurate measurement.

【0031】そこで、この発明では、以下に説明する演
算方法を用いることにした。 (実施例1)この実施例では、n、m番地の画素XのS
BR値を求めるのに、画素Xの近傍の出力値を演算に組
み入れることにした。
Therefore, in the present invention, the calculation method described below is used. (Embodiment 1) In this embodiment, S of pixel X at addresses n and m
In order to obtain the BR value, the output value near the pixel X is incorporated in the calculation.

【0032】図4、図5は時系列的に連続する一連の少
なくとも3枚の画像の一部分を拡大して示している。
4 and 5 are enlarged views showing a part of a series of at least three images which are continuous in time series.

【0033】ここで、各升目は画素を示し、ここでは、
4×4の画素が示されている。番地n、mの画素(一つ
の画像の中の一つの所定番地位置)XのSBR値に相当
する演算値を求める場合、図4に示すように、3枚の時
系列的画像についてそのn、m番地の画素Xの出力値A
i(i=1、2、3)とその画素Xの右斜め下の番地
(n+1)、(m+1)の画素の出力値Bi(i=1、
2、3)とからなる組み合わせと、図5に示すように、
3枚の時系列的画像について、その番地n、mの画素X
の右隣りの番地n、(m+1)の画素の出力値Di(i
=1、2、3)とその番地n、mの画素Xの直下の番地
(n+1)、mの画素の出力値Ci(i=1、2、3)
とからなる組み合わせとに分ける。
Here, each square represents a pixel, and here,
4x4 pixels are shown. When obtaining the calculated value corresponding to the SBR value of the pixel of address n, m (one predetermined address position in one image) X, as shown in FIG. Output value A of pixel X at address m
i (i = 1, 2, 3) and the output values B i (i = 1, 1) of the pixels at the addresses (n + 1) and (m + 1) at the lower right of the pixel X
2 and 3), and as shown in FIG.
For three time-series images, the pixel X of the address n, m
The output value D i (i of the pixel at the address n, (m + 1) on the right of
= 1, 2, 3) and the output value C i (i = 1, 2, 3) of the address (n + 1) and the pixel m immediately below the pixel X at the addresses n and m.
And a combination consisting of and.

【0034】そして、その図4に示す時系列的画像につ
いて2つの画素の組み合わせの中の出力値から差分をと
ると共に、その図5に示す時系列的画像について2つ画
素の組み合わせの中の出力値から差分をとる。
Then, the difference is calculated from the output values in the combination of two pixels for the time series image shown in FIG. 4 and the output in the combination of two pixels for the time series image shown in FIG. Take the difference from the value.

【0035】すなわち、一つの所定番地位置(n、m番
地)の画素X及び所定番地位置近傍に位置する画素の中
から選択された複数の画素からなる第1画像G1の第1
の画素群と、その他の画像(第2画像G2、第3画像G
3)の第1画素群と同じ番地に位置する第2及び第3画
素群とを選択する。そして、第1画素群ないし第3画素
群から所定の二つの組み合わせ画素を選択する。
That is, the first image G1 consisting of a plurality of pixels selected from the pixel X at one predetermined address position (n, m addresses) and the pixels located near the predetermined address position.
Pixel group and other images (second image G2, third image G
3) Select the second pixel group and the third pixel group located at the same address as the first pixel group. Then, two predetermined combination pixels are selected from the first pixel group to the third pixel group.

【0036】その図4に示す2つの画素の組み合わせの
中の出力値から差分をとる組み合わせとしては、3枚の
時系列的画像の各画像について、演算に用いる画素の個
数が2個であるので、演算に用いる画素の全個数は6個
であり、全ての各画素の組み合わせについて、その差分
を演算すると、62個数(15個)の出力値の差分結果
が得られることになる。この場合、差分の演算結果は
(A1−B1)、(A1−A2)、(A1−B2)、
(A1−A3)、(A1−B3)、(B1−A2)、
(B1−B2)、(B1−A3)、(B1−B3)、
(A2−B2)、(A2−A3)、(A2−B3)、
(B2−A3)、(B2−B3)、(A3−B3)とな
る。
As a combination for obtaining the difference from the output value in the combination of the two pixels shown in FIG. 4, the number of pixels used for the calculation is two for each of the three time-series images. The total number of pixels used for the calculation is 6, and when the difference is calculated for all combinations of each pixel, 6 C 2 (15) difference values of the output values are obtained. In this case, the difference calculation results are (A1-B1), (A1-A2), (A1-B2),
(A1-A3), (A1-B3), (B1-A2),
(B1-B2), (B1-A3), (B1-B3),
(A2-B2), (A2-A3), (A2-B3),
(B2-A3), (B2-B3), and (A3-B3).

【0037】同様にして、その図5に示す2つの画素の
組み合わせの中の出力値から差分をとる組み合わせとし
ては、3枚の時系列的画像の各画像について、演算に用
いる画素の個数が2個であるので、演算に用いる画素の
全個数は6個であり、全ての各画素の組み合わせについ
て、その差分を演算すると、62個数(15個)の出力
値の差分結果が得られることになる。
Similarly, as a combination for obtaining the difference from the output value in the combination of two pixels shown in FIG. 5, the number of pixels used for calculation is 2 for each of the three time-series images. Since the total number of pixels used in the calculation is 6, the difference between the output values of 6 C 2 (15) can be obtained by calculating the difference of all combinations of each pixel. become.

【0038】従って、演算に用いる全ての画素の個数は
12個であり、その出力の総和Sn, mは、(A1+A2
+A3)+(B1+B2+B3)+…+(D1+D2+
D3)である。この式は、番地n,mの画素Xの分子に
対応するもので、総和の平均値に係数2.5倍を掛けた
ものを分子の演算値Sn,m’とする。分母には、その3
0個の差分の総和の平均値を用いる。 Sn,m’=(A1+A2+A3)+…+(D1+D2+
D3)×2.5/12 ここで、係数2.5倍は、演算に用いた分母の画素の個
数と演算に用いた分子の画素の個数とを合わせたもので
ある。
Therefore, the total number of pixels used in the calculation is 12, and the total sum S n, m of the outputs is (A1 + A2
+ A3) + (B1 + B2 + B3) + ... + (D1 + D2 +
D3). This equation corresponds to the numerator of the pixel X at the address n, m, and the calculated sum Sn, m 'is obtained by multiplying the average value of the sums by a factor of 2.5. No. 3 in the denominator
The average value of the sum of 0 differences is used. S n, m '= (A1 + A2 + A3) + ... + (D1 + D2 +
D3) × 2.5 / 12 Here, the coefficient 2.5 times is the sum of the number of pixels of the denominator used for the calculation and the number of pixels of the numerator used for the calculation.

【0039】このような演算を、n=1、m=1番地か
らn=n、m=m番地まで行って、SBR値に相当する
演算値を得るものである。
The above operation is performed from the address n = 1, m = 1 to the address n = n, m = m to obtain the operation value corresponding to the SBR value.

【0040】なお、分母と分子とを逆にした逆数をSB
R値に基づく演算値として用いても良い。
The reciprocal of the inverse of the denominator and the numerator is SB
It may be used as a calculated value based on the R value.

【0041】ここでは、奇数フィールドに対応する画素
と偶数フィールドに対応する画素とで1/60秒ずれて
いる。
Here, the pixels corresponding to the odd fields and the pixels corresponding to the even fields are displaced by 1/60 seconds.

【0042】このような演算方法を用いると、血流速度
演算のための多数枚の画像を必要とせず、少なくとも3
枚の画像を用いて従来と同等の速度で血流速度を示すS
BR値に相当する演算値を求めることができる。 (実施例2)図7は図6に示す番地n、mの画素XのS
BR値に相当する演算値を求めるのに、図6に示す番地
n、mの画素Xの近傍の番地n、(m+1)の画素、番
地(n+1)、mの画素、番地(n+1)、(m+1)
の画素の出力値Ai〜Di(i=1、2、3)を演算に用
いることにしたものであり、この12個の選択された画
素のうちの2つを組み合わせてその差分を取ることにし
たものである。
When such a calculation method is used, a large number of images for calculating the blood flow velocity are not required, and at least 3
S that shows blood flow velocity at the same speed as the conventional one using one image
A calculated value corresponding to the BR value can be obtained. (Embodiment 2) FIG. 7 shows S of the pixel X at the addresses n and m shown in FIG.
To obtain the calculated value corresponding to the BR value, the addresses n, (m + 1) pixel, address (n + 1), m pixel, address (n + 1), ( m + 1)
The output values A i to D i (i = 1, 2, 3) of these pixels are used for the calculation, and two of the twelve selected pixels are combined to obtain the difference. It was decided.

【0043】演算に用いる画素の個数は12個であり、
差分を取るのに2個の画素を用いるので、差分の個数は
122(=66)個である。出力の総和Sn,mは、実施例
1と同じである。
The number of pixels used for the calculation is 12,
Since two pixels are used to take the difference, the number of differences is
It is 12 C 2 (= 66). The total output S n, m is the same as that in the first embodiment.

【0044】この実施例2では、分子の演算に用いる画
素の個数が12個であり、差分をとるのに用いた個数が
66個であるので、出力の総和Sn,mの平均値に掛ける
係数は、66/12=5.5である。
In the second embodiment, the number of pixels used for calculating the numerator is 12, and the number used for obtaining the difference is 66. Therefore, the average value of the total sum S n, m of the outputs is multiplied. The coefficient is 66/12 = 5.5.

【0045】この実施例2による場合も実施例1と同様
の効果を奏する。 (実施例3)図9、図10は図8に示す番地n、mの画
素XのSBR値に相当する演算値を求めるのに、図9に
示す(n+1)、m番地の画素の出力値Ci(i=1、
2、3)と(n+1)、(m+1)番地の画素の出力値
i(i=1、2、3)とからなる組み合わせと、図1
0に示す(n+1)、(m+1)番地の画素の出力値B
i(i=1、2、3)とn、(m+1)番地の画素の出
力値Di(i=1、2、3)とからなる組み合わせとを
用いて、出力値の差分と総和Sn,mとを求めたものであ
る。
In the case of the second embodiment, the same effect as that of the first embodiment can be obtained. (Embodiment 3) FIGS. 9 and 10 show the output value of the pixel at the address (n + 1), m shown in FIG. 9 for obtaining the calculated value corresponding to the SBR value of the pixel X at the address n, m shown in FIG. C i (i = 1,
2 and 3) and output values B i (i = 1, 2, 3) of the pixels at addresses (n + 1) and (m + 1), and FIG.
Output value B of the pixel at addresses (n + 1) and (m + 1) shown in 0
Using the combination of i (i = 1, 2, 3) and the output value D i (i = 1, 2, 3) of the pixel at address n, (m + 1), the difference between the output values and the sum S n , m and are obtained.

【0046】図9に示す画素を選択して、画素の出力値
の差分を求める組み合わせは、6個の画素から2個の組
み合わせを選択するので、62=15通り、図10に示
す画素を選択して、画素の出力値の差分を求める組み合
わせも、62=15通りであるので、分母の演算に用い
る画素の個数は30個である。
The pixel shown in FIG. 9 is selected and the combination for obtaining the difference between the output values of the pixels is selected from two combinations of 6 pixels. Therefore, 6 C 2 = 15 combinations, and the pixel shown in FIG. There are also 6 C 2 = 15 combinations for selecting the difference between the output values of the pixels by selecting, and therefore the number of pixels used for the denominator calculation is 30.

【0047】一方、総和Sn,mに用いる画素の個数は9
個であり、総和Sn,mは、 Sn,m=C1+B1+C2+B2+C3+B3+D1+
D2+D3 なお、係数は3.3である。
On the other hand, the number of pixels used for the sum S n, m is 9
And the total sum S n, m is S n, m = C1 + B1 + C2 + B2 + C3 + B3 + D1 +
D2 + D3 The coefficient is 3.3.

【0048】この実施例3の場合にも、実施例1と同様
の効果を奏する。 (実施例4)図12、図13は図11に示す番地n、m
の画素Xの値を求めるのに、図12に示す番地n、mの
画素の出力値Ai(i=1、2、3)と番地n、(n+
1)の画素の出力値Ci(i=1、2、3)とからなる
組み合わせと、図13に示す番地n、(m+1)の画素
の出力値Di(i=1、2、3)と番地(n+1)、
(m+1)の画素の出力値Bi(i=1、2、3)とか
らなる組み合わせとに分け、図12に示す6個の画素の
中から2つを選択してその出力値の差分を取ると共に、
図13に示す6個の画素の中から2つを選択してその出
力値の差分を取ることにしたものである。
Also in the case of the third embodiment, the same effect as that of the first embodiment is obtained. (Embodiment 4) FIGS. 12 and 13 show addresses n and m shown in FIG.
12 to obtain the value of the pixel X, the output values A i (i = 1, 2, 3) and the addresses n, (n +
1) the output value C i (i = 1, 2, 3) of the pixel and the output value D i (i = 1, 2, 3) of the pixel of address n, (m + 1) shown in FIG. And address (n + 1),
The output value B i (i = 1, 2, 3) of the (m + 1) pixel is divided into two combinations, and two of the six pixels shown in FIG. Along with taking
Two of the six pixels shown in FIG. 13 are selected and the difference between their output values is determined.

【0049】このものによる場合も、実施例3と同様に
分母の演算に用いる画素の個数は30個であり、総和S
n,mの演算に用いる画素の個数は12個であり、係数は
2.5である。
Also in this case, the number of pixels used for the calculation of the denominator is 30 as in the third embodiment, and the total sum S
The number of pixels used for the calculation of n and m is 12, and the coefficient is 2.5.

【0050】この実施例4の場合にも、実施例1と同様
の効果を奏する。 (実施例5)図15は図14に示す画素XのSBR値に
相当する演算値を演算するのに番地n、mの画素の出力
値Ai(i=1、2、3)、番地n、(m+2)の画素
の出力値Fi(i=1、2、3)、番地(n+1)、
(m+1)の画素の出力値Bi(i=1、2、3)、番
地(n+2)、mの画素の出力値Ei(i=1、2、
3)、番地(n+2)、(m+2)の画素の出力値Gi
(i=1、2、3)の5個の画素の出力値の中から2つ
を選択して、その差分を求めることにしたものである。
総和Sn,mは、 Sn,m=A1+A2+A3+B1+B2+B3+E1+
E2+E3+F1+F2+F3+G1+G2+G3 総和Sn,mの演算に用いる個数は15個であり、15個
の中から2個を選んで差分を取る組み合わせは、152
であるので、その分母の差分に用いる個数の組み合わせ
は105通りであり、従って、総和Sn,mに掛ける係数
は、7.0である。 (実施例6)図17は図16に示す画素XのSBR値に
相当する演算値を求めるのに番地n、(m+1)の画素
の出力値Di(i=1、2、3)、番地(n+1)、m
の画素の出力値Ci(i=1、2、3)、番地(n+
1)、(m+1)の画素の出力値Bi(i=1、2、
3)、番地(n+1)、(m+2)の画素の出力値Ji
(i=1、2、3)、番地(n+2)、(m+1)の画
素の出力値Hi(i=1、2、3)の5個の画素の出力
値の中から2つを選択して、その差分を求めることにし
たものである。総和Sn,mは、 Sn,m=B1+B2+B3+C1+C2+C3+D1+
D2+D3+H1+H2+H3+J1+J2+J3 総和の演算に用いる個数は15個であり、15個の中か
ら2個を選んで差分を取る組み合わせは、152である
ので、その分母の差分に用いる個数の組み合わせは10
5通りであり、従って、総和Sn,mに掛ける係数は、
7.0である。 (実施例7)図18に示すように、時系列的に連続する
複数枚の画像G1〜GNの中から時系列的に連続する3枚
の画像G1、G2、G3を取り出し、実施例1ないし実施
例6のいずれかの演算結果を用いて番地n,mの画素の
値をn=1、m=1からn=n、m=mまで演算し、そ
の演算結果に基づき血流速度分布を求める。
Also in the case of the fourth embodiment, the same effect as that of the first embodiment can be obtained. (Embodiment 5) FIG. 15 shows the output values A i (i = 1, 2, 3) and n of the pixels at addresses n and m for calculating the calculated value corresponding to the SBR value of the pixel X shown in FIG. , (M + 2) pixel output values F i (i = 1, 2, 3), address (n + 1),
Output value B i (i = 1, 2, 3) of pixel (m + 1), address (n + 2), output value E i of pixel m (i = 1, 2,
3), output values G i of pixels at addresses (n + 2) and (m + 2)
Two of the output values of the five pixels (i = 1, 2, 3) are selected, and the difference between them is determined.
The sum S n, m is S n, m = A1 + A2 + A3 + B1 + B2 + B3 + E1 +
E2 + E3 + F1 + F2 + F3 + G1 + G2 + G3 The number of the sum S n, m used for the calculation is 15, and the combination of selecting 2 out of 15 and taking the difference is 15 C 2
Therefore, there are 105 combinations of the numbers used for the difference of the denominator, and therefore, the coefficient by which the sum S n, m is multiplied is 7.0. (Embodiment 6) FIG. 17 shows output values D i (i = 1, 2, 3) of an address n and a pixel of (m + 1) for obtaining an operation value corresponding to the SBR value of the pixel X shown in FIG. (N + 1), m
Pixel output values C i (i = 1, 2, 3), address (n +
1), output values B i (i = 1, 2,
3), output values J i of pixels at addresses (n + 1) and (m + 2)
Two of the output values H i (i = 1, 2, 3) of the five pixels at (i = 1, 2, 3), the address (n + 2), and (m + 1) are selected. Then, the difference is determined. The total sum S n, m is S n, m = B1 + B2 + B3 + C1 + C2 + C3 + D1 +
D2 + D3 + H1 + H2 + H3 + J1 + J2 + J3 The number used to calculate the sum is 15, and the combination that selects two from 15 and takes the difference is 15 C 2 , so the combination of the number used for the difference of the denominator is 10
There are five ways, and therefore the coefficient by which the sum S n, m is multiplied is
It is 7.0. (Embodiment 7) As shown in FIG. 18, three time-sequential images G 1 , G 2 , and G 3 are taken out from a plurality of time-sequential images G 1 to G N. The value of the pixel at the address n, m is calculated from n = 1, m = 1 to n = n, m = m using the calculation result of any of the first to sixth embodiments, and based on the calculation result. Obtain the blood flow velocity distribution.

【0051】次いで、画像を1枚ずらして時系列的に連
続する3枚の画像G2、G3、G4を取り出し、次の血流
速度分布を求める。このように、血流速度分布の演算に
用いる画像を1枚ずつずらしながら3枚の画像に基づき
血流分布速度を求めると、時間的に変化する血流速度分
布が求められる。
Next, one image is shifted and three images G 2 , G 3 , and G 4 which are continuous in time series are taken out, and the next blood flow velocity distribution is obtained. As described above, when the blood flow distribution velocity is calculated based on the three images while shifting the images used for calculating the blood flow velocity distribution one by one, the blood flow velocity distribution that changes with time is obtained.

【0052】これによって、血流速度分布の時間的変化
を求めることができ、血流速度分布の変化を動画像とし
て画面1に表示させることができる。
With this, it is possible to obtain a temporal change of the blood flow velocity distribution and display the change of the blood flow velocity distribution on the screen 1 as a moving image.

【0053】[0053]

【発明の効果】本発明は、以上説明したように構成した
ので、従来の血流速度測定装置のように血流速度演算の
ための多数枚の画像を必要とせず、少なくとも3枚の画
像があれば充分な精度で血流速度を演算できるものであ
り、シャッター間隔が短い高速の連続画像取り込み装置
を必要とせず、例えば、1秒間60フレームの一般的な
画像取り込み装置であっても、充分高い測定精度を得る
ことができるという効果を奏する。
Since the present invention is configured as described above, it is not necessary to provide a large number of images for blood flow velocity calculation as in the conventional blood flow velocity measurement device, and at least three images can be displayed. If it can calculate the blood flow velocity with sufficient accuracy, it does not require a high-speed continuous image capturing device with a short shutter interval. For example, even a general image capturing device of 60 frames per second is sufficient. An effect that high measurement accuracy can be obtained is obtained.

【図面の簡単な説明】[Brief description of drawings]

【図1】 N枚の一連の画像のうちの番地n、mの画素
のSBR値の演算を、各画像の番地n、mの各画素の出
力値に基づき行うことを説明するための説明図である。
FIG. 1 is an explanatory diagram for explaining that an SBR value of pixels at addresses n and m in a series of N images is calculated based on output values of pixels at addresses n and m of each image. Is.

【図2】 図1に示す一連の画像を用いて求められたS
BR値に基づき画面に表示される血流速度分布の一例を
説明するための説明図である。
FIG. 2 is an S obtained using the series of images shown in FIG.
It is an explanatory view for explaining an example of a blood flow velocity distribution displayed on a screen based on a BR value.

【図3】 一連の時系列的画像のうちの時系列的に連続
する3枚の画像のうちのn、m番地の画素のSBR値に
相当する演算値を求めるのに用いる近傍画素群の説明図
である。
FIG. 3 is a description of a neighboring pixel group used to obtain a calculated value corresponding to an SBR value of pixels at addresses n and m of three images that are consecutive in time series among a series of time series images. It is a figure.

【図4】 本発明の実施例1に示すSBR値に相当する
演算値を求めるために二つの群に分けてその一方の群に
用いる画素の説明図であって、図3に示す画素群の4行
4列の画素を拡大して示す図である。
FIG. 4 is an explanatory diagram of pixels which are divided into two groups and used for one of the groups to obtain an operation value corresponding to the SBR value shown in the first embodiment of the present invention. It is a figure which expands and shows the pixel of 4 rows 4 columns.

【図5】 本発明の実施例1に示すSBR値に相当する
演算値を求めるために二つの群に分けてその他方の群に
用いる画素の説明図であって、図3に示す画素群の4行
4列の画素を拡大して示す図である。
5 is an explanatory diagram of pixels that are used in the other group by dividing them into two groups in order to obtain a calculated value corresponding to the SBR value shown in Example 1 of the present invention. It is a figure which expands and shows the pixel of 4 rows 4 columns.

【図6】 本発明の実施例2に示すSBR値に相当する
演算値を求める画素Xの説明図であって、図3に示す画
素群の4行4列の画素を拡大して示す図である。
FIG. 6 is an explanatory diagram of a pixel X for obtaining a calculated value corresponding to an SBR value according to the second embodiment of the present invention, which is an enlarged view of a pixel in row 4, column 4 of the pixel group shown in FIG. 3; is there.

【図7】 本発明の実施例2に示すSBR値に相当する
演算値を求めるのに使用する画素Xの近傍の説明図であ
る。
FIG. 7 is an explanatory diagram in the vicinity of a pixel X used for obtaining a calculated value corresponding to an SBR value shown in Example 2 of the present invention.

【図8】 本発明の実施例3に示すSBR値に相当する
演算値を求める画素Xの説明図であって、図3に示す画
素群の4行4列の画素を拡大して示す図である。
FIG. 8 is an explanatory diagram of a pixel X for obtaining a calculated value corresponding to an SBR value according to the third embodiment of the present invention, which is an enlarged view of a pixel in row 4 and column 4 of the pixel group shown in FIG. is there.

【図9】 本発明の実施例3に示すSBR値に相当する
演算値を求めるために二つの群に分けてその一方の群に
用いる画素の説明図である。
FIG. 9 is an explanatory diagram of pixels which are divided into two groups and used in one of the groups to obtain a calculated value corresponding to an SBR value shown in Example 3 of the present invention.

【図10】 本発明の実施例3に示すSBR値に相当す
る演算値を求めるために二つの群に分けてその他方の群
に用いる画素の説明図である。
FIG. 10 is an explanatory diagram of pixels which are divided into two groups and are used in the other group in order to obtain a calculated value corresponding to the SBR value shown in Example 3 of the present invention.

【図11】 本発明の実施例4に示すSBR値に相当す
る演算値を求める画素Xの説明図であって、図3に示す
画素群の4行4列の画素を拡大して示す図である。
FIG. 11 is an explanatory diagram of a pixel X for obtaining a calculated value corresponding to an SBR value according to the fourth embodiment of the present invention, which is an enlarged view of a pixel in row 4, column 4 of the pixel group shown in FIG. 3; is there.

【図12】 本発明の実施例4に示すSBR値に相当す
る演算値を求めるために二つの群に分けてその一方の群
に用いる画素の説明図である。
FIG. 12 is an explanatory diagram of pixels which are divided into two groups and used in one of the groups to obtain a calculated value corresponding to an SBR value shown in Example 4 of the present invention.

【図13】 本発明の実施例4に示すSBR値に相当す
る演算値を求めるために二つの群に分けてその他方の群
に用いる画素の説明図である。
FIG. 13 is an explanatory diagram of pixels which are divided into two groups and are used in the other group in order to obtain a calculated value corresponding to the SBR value shown in Example 4 of the present invention.

【図14】 本発明の実施例5に示すSBR値に相当す
る演算値を求める画素Xの説明図であって、図3に示す
画素群の4行4列の画素を拡大して示す図である。
FIG. 14 is an explanatory diagram of a pixel X for obtaining a calculation value corresponding to an SBR value according to the fifth embodiment of the present invention, which is an enlarged view of a pixel in row 4 and column 4 of the pixel group shown in FIG. is there.

【図15】 本発明の実施例5に示すSBR値に相当す
る演算値を求めるのに使用する画素Xの近傍の説明図で
ある。
FIG. 15 is an explanatory diagram of the vicinity of a pixel X used to obtain a calculated value corresponding to an SBR value shown in Example 5 of the present invention.

【図16】 本発明の実施例6に示すSBR値に相当す
る演算値を求める画素Xの説明図であって、図3に示す
画素群の4行4列の画素を拡大して示す図である。
FIG. 16 is an explanatory diagram of a pixel X for obtaining a calculated value corresponding to an SBR value according to the sixth embodiment of the present invention, which is an enlarged view of a pixel in row 4 and column 4 of the pixel group shown in FIG. is there.

【図17】 本発明の実施例6に示すSBR値に相当す
る演算値を求めるのに使用する画素Xの近傍の説明図で
ある。
FIG. 17 is an explanatory diagram in the vicinity of a pixel X used to obtain a calculated value corresponding to an SBR value shown in Example 6 of the present invention.

【図18】 本発明の実施例7に示す時系列的画像群の
説明図である。
FIG. 18 is an explanatory diagram of a time-series image group according to the seventh embodiment of the present invention.

【図19】 本発明の血流速度測定装置の要部構成を示
す図である。
FIG. 19 is a diagram showing a main configuration of a blood flow velocity measuring device of the present invention.

【符号の説明】[Explanation of symbols]

1…レーザー光照射系(レーザー光照射部) 5…CCD(受光部) 16…画像記録器(記憶部) 17…演算処理部 1 ... Laser light irradiation system (laser light irradiation part) 5 ... CCD (light receiving part) 16 ... Image recorder (storage unit) 17 ... Arithmetic processing unit

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 FI テーマコート゛(参考) G06T 7/20 A61B 5/02 340D (72)発明者 藤居 仁 福岡県宗像市日の里2丁目26番8 (72)発明者 小西 直樹 福岡県嘉穂郡庄内町仁保232−15 Fターム(参考) 2G059 AA05 BB12 CC16 EE02 GG01 JJ22 JJ23 KK04 MM01 MM09 MM10 4C017 AA11 AC28 BC11 5B057 AA10 BA02 CA08 CA12 CA16 CD01 DA20 DB02 DB09 DC32 5L096 AA06 BA06 BA13 CA04 FA25─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 7 identification code FI theme code (reference) G06T 7/20 A61B 5/02 340D (72) Inventor Hitoshi Fujii 2-26, Hisato, Munakata-shi, Fukuoka 8 (72) Inventor Naoki Konishi, 232-15 F term, Niho, Shonai-cho, Kaho-gun, Fukuoka Prefecture (reference) 2G059 AA05 BB12 CC16 EE02 GG01 JJ22 JJ23 KK04 MM01 MM09 MM10 4C017 AA11 AC28 BC11 5B057 AA10 BA02 CA08 CA20 DB16 DB01 CD01 DA01 5L096 AA06 BA06 BA13 CA04 FA25

Claims (4)

【特許請求の範囲】[Claims] 【請求項1】 生体組織の血球にレーザ光を照射するレ
ーザ光照射部と、前記生体組織からの反射光を検出する
多数の受光素子からなる受光部と、前記受光部からの信
号に基づき受光部上に形成された画像を順次複数枚記録
する記憶部と、前記記憶された画像の中の少なくとも3
枚の画像に基づき血球の速度の分布状態を演算する演算
部とからなり、 前記演算部は、一つの画像の中の一つの所定番地位置の
画素及び該所定番地近傍に位置する画素の中から選択さ
れた複数の画素からなる第1の画素群と、その他の画像
の中で前記第1の画素群とそれぞれ同じ番地に位置する
第2及び第3の画素群とを選択し、前記第1ないし第3
画素群から所定の2つの組み合わせ画素を選択し、選択
された2つの画素の出力の差分の総和の平均値もしくは
その逆数を演算し、その差分演算を各所定番地の画素毎
に行うことにより、生体組織での血流分布状態を検出す
ることを特徴とする血流速度測定装置。
1. A laser light irradiating section for irradiating blood cells of living tissue with laser light, a light receiving section composed of a plurality of light receiving elements for detecting reflected light from the living tissue, and light reception based on a signal from the light receiving section. A storage unit for sequentially recording a plurality of images formed on a set, and at least 3 of the stored images
Comprising a calculation unit for calculating the distribution state of the blood cell velocity based on a single image, the calculation unit, from the pixel at one predetermined address position in one image and the pixel located near the predetermined address The first pixel group including a plurality of selected pixels and the second and third pixel groups located at the same address as the first pixel group in other images are selected, and the first pixel group is selected. Or third
By selecting two predetermined combination pixels from the pixel group, calculating the average value of the sum of the differences between the outputs of the two selected pixels or the reciprocal thereof, and performing the difference calculation for each pixel at each predetermined address, A blood flow velocity measuring device characterized by detecting a blood flow distribution state in a biological tissue.
【請求項2】前記受光部からの信号に基づき奇数走査線
に基づく奇数フィールド画像と偶数走査線の信号に基づ
く偶数フィールド画像とを1組として順次取り込み、記
憶部に取り込んだ前記偶数フィールド画像と前記奇数フ
ィールド画像とを合成して合成フレーム画像として記憶
し、前記演算部は前記記憶された画像の中の少なくとも
3枚の合成フレーム画像に基づき血球の速度の分布状態
を演算する請求項1記載の血流速度測定装置。
2. An even field image sequentially captured as a set of an odd field image based on an odd scanning line and an even field image based on a signal of an even scanning line based on a signal from the light receiving unit, and the even field image captured in a storage unit. 2. The odd-numbered field image is combined and stored as a combined frame image, and the calculation unit calculates the distribution state of blood cell velocities based on at least three combined frame images in the stored images. Blood flow velocity measuring device.
【請求項3】 前記演算部で演算に使用する画像は、受
光素子上での像のずれ量を補正するために座標変換を行
った画像であることを特徴とする請求項1記載の血流速
度測定装置。
3. The blood flow according to claim 1, wherein the image used in the calculation by the calculation unit is an image that has been subjected to coordinate conversion in order to correct the image shift amount on the light receiving element. Speed measuring device.
【請求項4】 第1ないし第3の画素群はさらに、それ
ぞれ画素番地位置が異なる2つの画素群に分けられ、該
2つの画素群のそれぞれの画素群の中で所定の2つの組
み合わせ画素を選択する請求項1記載の血流速度測定装
置。
4. The first to third pixel groups are further divided into two pixel groups having different pixel address positions, and two predetermined combination pixels are included in each of the two pixel groups. The blood flow velocity measuring device according to claim 1, which is selected.
JP2001379858A 2001-12-13 2001-12-13 Blood flow velocity measuring device Expired - Lifetime JP3855024B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001379858A JP3855024B2 (en) 2001-12-13 2001-12-13 Blood flow velocity measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001379858A JP3855024B2 (en) 2001-12-13 2001-12-13 Blood flow velocity measuring device

Publications (2)

Publication Number Publication Date
JP2003180641A true JP2003180641A (en) 2003-07-02
JP3855024B2 JP3855024B2 (en) 2006-12-06

Family

ID=27591197

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001379858A Expired - Lifetime JP3855024B2 (en) 2001-12-13 2001-12-13 Blood flow velocity measuring device

Country Status (1)

Country Link
JP (1) JP3855024B2 (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006061328A (en) * 2004-08-26 2006-03-09 Kowa Co Ophthalmologic apparatus
US7597443B2 (en) 2004-10-28 2009-10-06 Kyushu Institute Of Technology Wide viewing angle ocular fundus blood flow imaging device
JP2010110371A (en) * 2008-11-04 2010-05-20 Topcon Corp Apparatus and method for measuring eye, and program
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 (en) 2009-05-13 2010-11-18 国立大学法人九州工業大学 Blood flow image diagnosing device
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
US8733933B2 (en) 2010-12-02 2014-05-27 Canon Kabushiki Kaisha Ophthalmologic apparatus and blood flow velocity calculation method
US8770752B2 (en) 2010-12-02 2014-07-08 Canon Kabushiki Kaisha Ophthalmic apparatus, ophthalmic system, processing apparatus, and blood flow velocity calculation method
WO2014175154A1 (en) 2013-04-23 2014-10-30 ソフトケア有限会社 Blood flow image diagnosis device and diagnosis method
CN104968511A (en) * 2014-01-31 2015-10-07 横滨橡胶株式会社 Stud pin and pneumatic tire
EP3213668A2 (en) 2016-03-03 2017-09-06 Nidek Co., Ltd. Ophthalmic image processing apparatus and ophthalmic image processing program
WO2018003139A1 (en) 2016-06-28 2018-01-04 ソフトケア有限会社 Blood flow dynamic imaging diagnosis device and diagnosis method
WO2018203429A1 (en) * 2017-05-02 2018-11-08 ソニー株式会社 Flow velocity measurement method, flow velocity measurement device, and program
JP2019195586A (en) * 2018-05-11 2019-11-14 キヤノン株式会社 Image processing device, image processing method, and program

Cited By (28)

* Cited by examiner, † Cited by third party
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
JP2006061328A (en) * 2004-08-26 2006-03-09 Kowa Co Ophthalmologic apparatus
JP4585814B2 (en) * 2004-08-26 2010-11-24 興和株式会社 Ophthalmic equipment
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
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
JP2010110371A (en) * 2008-11-04 2010-05-20 Topcon Corp Apparatus and method for measuring eye, and program
WO2010131550A1 (en) 2009-05-13 2010-11-18 国立大学法人九州工業大学 Blood flow image diagnosing device
US9028421B2 (en) 2009-05-13 2015-05-12 Kyushu Institute Of Technology Blood flow image diagnosing device
KR101477084B1 (en) * 2010-12-02 2014-12-29 캐논 가부시끼가이샤 Ophthalmic apparatus, ophthalmic system, processing apparatus, and blood flow velocity calculation method
US8770752B2 (en) 2010-12-02 2014-07-08 Canon Kabushiki Kaisha Ophthalmic apparatus, ophthalmic system, processing apparatus, and blood flow velocity calculation method
US8733933B2 (en) 2010-12-02 2014-05-27 Canon Kabushiki Kaisha Ophthalmologic apparatus and blood flow velocity calculation method
KR101447010B1 (en) 2010-12-02 2014-10-06 캐논 가부시끼가이샤 Ophthalmologic apparatus and blood flow velocity calculation method
US10098592B2 (en) 2013-04-23 2018-10-16 Softcare Co., Ltd. Blood flow image diagnosing device and method
WO2014175154A1 (en) 2013-04-23 2014-10-30 ソフトケア有限会社 Blood flow image diagnosis device and diagnosis method
CN104968511A (en) * 2014-01-31 2015-10-07 横滨橡胶株式会社 Stud pin and pneumatic tire
CN104968511B (en) * 2014-01-31 2018-06-05 横滨橡胶株式会社 Anti-skid stud and pneumatic tire
EP3213668A2 (en) 2016-03-03 2017-09-06 Nidek Co., Ltd. Ophthalmic image processing apparatus and ophthalmic image processing program
WO2018003139A1 (en) 2016-06-28 2018-01-04 ソフトケア有限会社 Blood flow dynamic imaging diagnosis device and diagnosis method
US11330995B2 (en) 2016-06-28 2022-05-17 Softcare Co., Ltd. Apparatus and method for imaging and analyzing hemodynamics
WO2018203429A1 (en) * 2017-05-02 2018-11-08 ソニー株式会社 Flow velocity measurement method, flow velocity measurement device, and program
JPWO2018203429A1 (en) * 2017-05-02 2020-05-14 ソニー株式会社 Flow velocity measuring method, flow velocity measuring device and program
JP6999656B2 (en) 2017-05-02 2022-01-18 ソニーグループ株式会社 Flow velocity measuring method, flow velocity measuring device and program
US11231436B2 (en) 2017-05-02 2022-01-25 Sony Corporation Flow rate measuring method, flow rate measuring device, and program
JP2019195586A (en) * 2018-05-11 2019-11-14 キヤノン株式会社 Image processing device, image processing method, and program
JP7281872B2 (en) 2018-05-11 2023-05-26 キヤノン株式会社 Image processing device, image processing method and program

Also Published As

Publication number Publication date
JP3855024B2 (en) 2006-12-06

Similar Documents

Publication Publication Date Title
JP3855024B2 (en) Blood flow velocity measuring device
AU608807B2 (en) Apparatus for monitoring bloodstream
JP5294340B2 (en) Ultrasonic diagnostic equipment
EP2075722B1 (en) Ultrasound diagnostic device and method of displaying ultrasound images
JP2003164431A (en) Blood flow velocity measuring apparatus
CN100413468C (en) Ultrasonic diagnostic apparatus
US20200301121A1 (en) Method for high-resolution scanning microscopy
JPH0528133B2 (en)
JP5176081B2 (en) X-ray imaging method and X-ray imaging apparatus
JP4492107B2 (en) Diagnostic imaging equipment
JPH01124437A (en) Blood flow monitor apparatus
JP2006023494A (en) Method for obtaining microscope observation image and microscope system
JPH0243486B2 (en)
JP6763086B2 (en) Methods and devices for improving frame rates in infrared sensor systems
JPS63240832A (en) Image measuring method for calculating contour line of lung part
US7436929B2 (en) Radiographic system
JP3077773B2 (en) 3D image display
JP3643398B2 (en) Arrival order and time determination device
JP3699088B2 (en) Arrival order and time determination device
JPH0599623A (en) Displacement measuring apparatus
JP3858094B2 (en) Image capturing method and apparatus for reproducible phenomenon
JPH11151234A (en) X-ray imaging instrument
JPH07151516A (en) Speckle displacement gauge
JPH0627383A (en) Method and device for picture display
JP3087132B2 (en) Digital video conversion processing system

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

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051201

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20051201

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20051201

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20060227

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20060227

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20060228

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060530

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060613

R155 Notification before disposition of declining of application

Free format text: JAPANESE INTERMEDIATE CODE: R155

R150 Certificate of patent or registration of utility model

Ref document number: 3855024

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

EXPY Cancellation because of completion of term