CN110279393A - 微血管检测装置和方法 - Google Patents

微血管检测装置和方法 Download PDF

Info

Publication number
CN110279393A
CN110279393A CN201810222949.4A CN201810222949A CN110279393A CN 110279393 A CN110279393 A CN 110279393A CN 201810222949 A CN201810222949 A CN 201810222949A CN 110279393 A CN110279393 A CN 110279393A
Authority
CN
China
Prior art keywords
capilary
point
value
digitized video
mark point
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
CN201810222949.4A
Other languages
English (en)
Other versions
CN110279393B (zh
Inventor
林孟癸
侯军纬
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.)
Young Fast Optoelectronics Co Ltd
Original Assignee
Young Fast Optoelectronics Co Ltd
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 Young Fast Optoelectronics Co Ltd filed Critical Young Fast Optoelectronics Co Ltd
Priority to CN201810222949.4A priority Critical patent/CN110279393B/zh
Publication of CN110279393A publication Critical patent/CN110279393A/zh
Application granted granted Critical
Publication of CN110279393B publication Critical patent/CN110279393B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0062Arrangements for scanning
    • A61B5/0064Body surface scanning
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0077Devices for viewing the surface of the body, e.g. camera, magnifying lens
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1075Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1079Measuring physical dimensions, e.g. size of the entire body or parts thereof using optical or photographic means

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Hematology (AREA)
  • Image Analysis (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明公开一种微血管检测装置,待测一手指放在一指槽中,经由一手指皮下组织中的至少一微血管影像,检测所述微血管的血流流速以及管径,包含有:一计算机,具有一显示器以及一处理器;一感光耦合组件,电性讯号链接该计算机;以及一显微镜镜头,经由显微镜镜头撷取所述微血管影像,该微血管影像由该感光耦合组件形成复数帧数字影像,其中时间连续的复数帧该数字影像,经由该处理器显示于该显示器。

Description

微血管检测装置和方法
技术领域
本发明涉及一种人体血管的检测装置和方法,更具体的说,是指一种检测微血管的血流流速以及管径的微血管检测装置和方法。
背景技术
现有技术的微血管检测,如中国台湾专利号为I246910的专利,该发明提出一种直接实时检测微血管血流流速的方法与微循环功能的评估装置。包括利用红外光雷射血管显微摄影仪的影像动画,以指针选取影像中特定微血管分支,沿着微血管纵向标定分析范围,透过连续的动画处理,可以绘制实时红血球位移影像。侦测实时血球位移影像的斜率变化可以分析红血球的位移速度变化及加速度。综合此部位各个微血管流速,进行统计分析,可以观察微血管群的特性差异。
上述常用的红外光雷射血管显微摄影仪,用灰阶判定的方式容易有误差;再者,计算路径方式单一化,无其他计算路径方式变化来做为参考;此外,目前用于微血管检测的装置大都昂贵的专用仪器,无法普及应于一般人的居家护理使用。
有鉴于此,本发明人乃潜心研思、设计组制,期能提供一种廉价且可简便操作以便快速检测微血管的装置和方法,其以简易低价的显微摄像装置搭配一般家庭计算机使用即可快速地自我检测微血管的血流流速以及管径,让使用者随时自我简易检测评估血液循环状态,以便注意维护身体健康。
发明内容
本发明的主要目的,在于提供一种利用白血球定位以及像素运算,来达到检测微血管的血流流速以及管径的目的。
为达上述目的,本发明的一实施例为一种微血管检测装置,经由一手指皮下组织中的至少一微血管影像,检测所述微血管的血流流速以及管径,包含有:一计算机,具有一显示器以及一处理器;一感光耦合组件,电性讯号链接该计算机;以及一显微镜镜头,经由该显微镜镜头撷取该微血管影像,该微血管影像由该感光耦合组件形成复数帧数字影像,其中时间连续的复数帧该数字影像,经由该处理器显示于该显示器。
所述检测装置在一实施例中,该处理器标定复数帧该数字影像,对应该微血管中一白血球的时间连续的标示点,包括一起点标示点以及一终点标示点,该起点标示点以及该终点标示点的时间差为第1时间差,该处理器计算加总连续标示点的第1路径长,第1路径长除以第1时间差的一血流流速值,于该显示器显示该血流流速值。
所述检测装置在一实施例中,该处理器标定复数帧该数字影像,对应该微血管中一白血球的时间连续的标示点,包括一起点标示点以及一终点标示点,该处理器具有一工字型框架的模块,包含有45度标点、90度标点、以及135度标点的路经搜寻到一最大边缘位置,该处理器计算一血管中心标示点,由该血管中心标示点,经由该工字型框架的模块,依次标示出至少一计算标示点,直到该终点标示点,该起点标示点以及该终点标示点的时间差为第2时间差,该处理器计算加总连续的该计算标示点的第2路径长,除以第2时间差的一血流流速值,于该显示器显示该血流流速值。
所述检测装置在一实施例中,该处理器扫描复数帧该数字影像成灰阶讯号,由纵轴灰阶讯号加总值最大值,标定该微血管管径的两边缘端点,该处理器由横轴该微血管的两边缘端点的相应像素值,计算该微血管的一管径值,于该显示器显示该微血管的该管径值。
本发明的另一实施例为一种微血管检测方法,分解该数字影像作为量测一微血管的一血流流速值,其中,检测步骤,包含有:点选该数字影像中一白血球起始位置;找出时间连续的该数字影像中同一个该白血球的位置,并点选该白血球的位置,搜寻到路径,则标示路径并计算该血流流速值;以及搜寻不到路径,则显示错误,并重新回到点选该数字影像中该白血球起始位置的步骤。
在一实施例中,分解该数字影像作为量测一微血管的一管径值,其中,检测步骤,包含有:点选该数字影像中该微血管内的任一位置;找出该微血管内的两边缘端点位置,并量测该微血管的管径;以及手动微调量测位置,并显示该微血管的该管径值于该显示器。
所述检测方法在一实施例中,分解该数字影像作为量测一微血管的一血流流速值,其中,同一个该白血球在时间连续的该数字影像中的位置,被标定为一起点标示点以及一终点标示点、以及至少一追踪点,计算所有该追踪点的一平均位置点,各追踪点和平均位置的角度,依角度进行排序,得到该白血球流经的顺序,该起点标示点以及该终点标示点的时间差为第1时间差,依序加总相邻距离,得到流经该微血管的第1路径长,第1路径长除以第1时间差的一血流流速值;以及于该显示器显示该血流流速值。
所述检测方法在一实施例中,分解该数字影像作为量测一微血管的一血流流速值,其中,检测步骤,更包含有:点选该数字影像中该白血球的一起点标示点以及一终点标示点;使用一工字型框架,以45度标点、90度标点、以及135度标点的路经搜寻到一最大边缘位置;计算一血管中心标示点;由该血管中心标示点,经由该工字型框架依次标示出至少一计算标示点,直到该终点标示点;该起点标示点以及该终点标示点的时间差为第2时间差,计算加总连续的该计算标示点的第2路径长,除以第2时间差的一血流流速值;以及于该显示器显示该血流流速值。
所述检测方法在一实施例中,找出该微血管内的两边缘端点位置,其步骤,包含:由扫描该数字影像成灰阶讯号,形成纵轴灰阶讯号值,横轴像素值;纵轴灰阶讯号加总值最大值,标定该微血管管径的两边缘端点;该微血管管径的两边缘端点的相应横轴像素值,计算该微血管的该管径值;以及于该显示器显示该微血管的该管径值。
本“发明内容”是以简化形式介绍一些选定概念,在下文的“具体实施方式”中将进一步对其进行描述。本“发明内容”并非意欲辨识申请专利的关键特征或基本特征,亦非意欲用于限制申请专利的保护范围。
附图说明
图1为本发明装置示意图。
图2为本发明起点标示点示意图。
图3为本发明终点标示点示意图。
图4为本发明以白血球解析路径示意图。
图5为本发明以沿血管边缘解析路径示意图。
图6为本发明量测管径示意图。
图7为本发明微血管血流流速检测方法步骤图。
图8为本发明以白血球析路径的分解画面示意图。
图9为本发明微血管血流流速检测方法步骤图。
图10为本发明由平均位置的角度,计算流经路径长示意图。
图11为本发明以沿血管边缘解析路径的分解画面示意图。
图12为本发明以血管边缘解析路径的初始搜寻路径。
图13A为本发明工字型框架135度搜寻示意图。
图13B为本发明工字型框架90度搜寻示意图。
图13C为本发明工字型框架45度搜寻示意图。
图14为本发明最大边缘位置中心定位示意图。
图15为本发明以沿血管边缘解析路径的前进标示示意图。
图16A为本发明工字型框架135度搜寻示意图。
图16B为本发明工字型框架90度搜寻示意图。
图16C为本发明工字型框架45度搜寻示意图。
图17为本发明最大边缘位置中心定位示意图。
图18为本发明以血管边缘解析路径的轨迹图。
图19为本发明以沿血管边缘解析路径血流流速计算示意图。
图20为本发明管径量测示意图。
符号说明:
11 计算机
12 显示器
13 感光耦合组件
14 显微镜镜头
15 指槽
16 手指
17 处理器
21 微血管
22 白血球
24 第1标线
25 第2标线
81 起点标示点
82 终点标示点
85 135度标点
86 90度标点
87 45度标点
88 工字型框架
89 最大边缘位置
90 血管中心标示点
91 计算标示点
92 范围外白血球
101 平均位置点
102 第1追踪点
103 第2追踪点
104 第3追踪点
105 第4追踪点
106 第5追踪点
107 第6追踪点
110 第1分布图
111 第1曲线图
120 第2分布图
121 第2曲线图
具体实施方式
图1所示,为本发明装置示意图。在一实施例中,本发明提供一种微血管检测装置,经由一手指16皮下组织中的至少一微血管21影像,检测该微血管21的血流流速以及管径,本发明装置包含有:计算机11、感光耦合组件13(charge-coupled device,CCD)、以及显微镜镜头14。其中该计算机11,具有一显示器12以及一处理器17;该感光耦合组件13,电性讯号链接该计算机11,例如使用通讯传输接口的USB接口、IEEE 1394接口、Ethernet接口、以及CVBS加上影像撷取卡接口等。以及该显微镜镜头14具有放大显微影像功能,经由该显微镜镜头14撷取该微血管21影像,该微血管21影像由该感光耦合组件13形成复数帧数字影像,其中时间连续的复数帧该数字影像,经由该处理器17显示于该显示器12。
图2和图3所示,该处理器17标定复数帧该数字影像,对应该微血管21中一白血球22的时间连续的标示点,包括一起点标示点81以及一终点标示点82,该起点标示点81以及该终点标示点82的时间差为第1时间差,该处理器17计算加总连续标示点的第1路径长,第1路径长除以第1时间差的一血流流速值,于该显示器12显示该血流流速值。
图4和图5所示为本发明两种血流流速值的显示结构,图4所示为本发明白血球22解析路径,找出若干分解图中被选取血管的白血球22,并分析其流动路径。以及图5所示为本发明沿微血管21边缘解析路径。从微血管21内一点选位置出发,沿微血管21边缘,找到另一点选位置。由流动路径长度,及总共有的分解数字影像的图数,即可得到总时间,以便可计算出血流流速。
图6所示,是量测管径示意图。为本发明一实施例,经由处理器17计算,显示在显示器12的功能显示画面,包括显示微血管21的管径值,以及微调左(left)的第1标线24的边境值,微调右(right)的第2标线25的边境值,以及使用鼠标点选第1标线24、以及第2标线25的位置后,处理器17计算微血管21的管径值,图6显示画面,与图7中步骤S71、S72、以及S73的共享该显示器12的功能。显示在显示器12的功能显示画面的下端,是图形功能接口的选项。
图7所示,在另一实施中,是本发明微血管21血流流速检测方法步骤图。其中,检测步骤,包含有:步骤S1,待测一手指16放在一指槽15中;步骤S2,使一感光耦合组件13取得一数字影像;步骤S3,经传输将该数字影像显示在一计算机11的一显示器12;步骤S4,依该显示器12的该数字影像调整该手指16的位置及一显微镜镜头14的焦距;以及步骤S5,撷取该数字影像并依连续时间分解该数字影像。
图7所示,由步骤S6,功能选择,来选择分解该数字影像作为量测一微血管21的一血流流速值,其中,检测步骤,包含有:步骤S81,点选该数字影像中一白血球22起始位置;步骤S82,找出时间连续的该数字影像中同一个该白血球22的位置,并点选该白血球22的位置;步骤S85,搜寻到路径,则标示路径并计算该血流流速值;以及步骤S84,搜寻不到路径,则显示错误,并重新回到点选该数字影像中该白血球22起始位置的步骤。
图7所示,由步骤S6,功能选择,来选择分解该数字影像作为量测一微血管21的一管径值,其中,检测步骤,包含有:步骤S71,点选该数字影像中该微血管21内的任一位置;步骤S72,找出该微血管21内的两边缘端点位置,并量测该微血管21的管径;以及步骤S73,手动微调量测位置,并显示该微血管21的该管径值于该显示器12。
图7至图10所示,分解该数字影像作为量测一微血管21的一血流流速值,其中,同一个该白血球22在时间连续的该数字影像中的位置,被标定为一起点标示点81以及一终点标示点82、以及至少一追踪点,图10所示,第1追踪点102、第2追踪点103、到第6追踪点107;计算所有该追踪点的一平均位置点101,各追踪点和平均位置的角度,依角度进行排序;得到该白血球22流经的顺序,该起点标示点81以及该终点标示点82的时间差为第1时间差;依序加总相邻距离,得到流经该微血管21的第1路径长;第1路径长除以第1时间差的一血流流速值;以及于该显示器12显示该血流流速值。
图9中,同一白血球22的选取,以影像分析各分解图间的灰阶差异,找出其间所有白血球22位置,并去除选取范围外白血球92。
图8中,使用者选取分解画面的白血球22起点,并在后续几张分解画面中选取白血球22终点,依计算机11效能,录像压缩后,可以从数字图像文件案得知图框率(FrameRate),以下以每秒25张的图框率(Frame Rate),举例:每张为1/25秒,分解图的数字影像有9张,头尾间隔8/25秒。图10中,假设上述路径长计算得96像素(pixel),而显微镜镜头14撷取的数字影像,像素比为1.164594μm/pixel.,第2路径长是96像素(pixel)乘以1.164594μm/pixel=111.801024μm,显微镜镜头14撷取的数字影像,每秒有25张数字影像,由此算出该起点标示点81以及该终点标示点82的时间差为第1时间差,因此第2路径长除以第1时间差,可以计算血流流速,[96pixels/(8/25sec)]X 1.164594μm/pixel=349μm/sec≒0.35mm/sec。
图11至图9中,是本发明一实施例工字型框架88计算血流流速值。本发明分解该数字影像作为量测一微血管21的一血流流速值,其中,检测步骤,更包含有:图11所示,点选该数字影像中该白血球22的一起点标示点81以及一终点标示点82;图12至图13C所示,使用一工字型框架88,以45度标点87、90度标点86、以及135度标点85的路经,图14所示,搜寻到一最大边缘位置89;图14所示,计算一血管中心标示点90;由该血管中心标示点90,图15至图17所示,经由该工字型框架88依次标示出至少一计算标示点91,直到该终点标示点82;该起点标示点81以及该终点标示点82的时间差为第2时间差,图18所示,计算加总连续的该计算标示点91的第2路径长,除以第2时间差的一血流流速值;以及于该显示器12显示该血流流速值。
图12所示,是由起点标示点81或终点标示点82较低者搜寻,初始的搜寻方向是向上的90度方向。本发明,内定每一步的搜寻角度是+/-45度。所以下一步的候选位置为135、90,45度,如图13A至图13C所示,三点位置。由左至右分别是往135、90,45度前进的位置。工字型框架88为搜寻范围,从血管外部往中心方向计算该角度的最大边缘量,其中,最大边缘量定义是依点选位置沿着工字型框架88的工字型左右两侧的箭头方向,从血管外部往中心方向各扫描若干像素(pixels)的范围,扫描方法先产生邻近像素间灰阶值差异,并将工字型框架88的工字型中轴方向的灰阶值加总,结果得到工字型框架88的工字型左右侧的曲线图,找出左右侧曲线中的坡峰,即差异最大处,为该角度的左侧的最大边缘量(Lmax),以及右侧的最大边缘量(Rmax)。
左右两侧最大边缘量的总和(Lmax+Rmax),是此角度的总边缘量。如图13A实施例中135度的右侧的最大边缘量(Rmax)=181,左侧的最大边缘量(Lmax)=29,因此,135度的总边缘量181+29=210。如图13B实施例中90度的右侧的最大边缘量(Rmax)=354,左侧的最大边缘量(Lmax)=672,90度的总边缘量672+354=1026。如图13C实施例中45度的右侧的最大边缘量(Rmax)=229,左侧的最大边缘量(Lmax)=243,因此,45度的总边缘量229+243=472。由上述计算,90度总边缘量大于45度总边缘量,大于135度总边缘量,所以最大总边缘量是90度总边缘量的值,即1026。由于在此阶段的最大总边缘量是90度的值,即1026,所以此阶段的方向往最大总边缘量,即90度前进。
工字型框架88两侧各有一个最大边缘量,两侧最大边缘量的总和,是此角度的总边缘量。假设此阶段的最大总边缘量是90度,所以此阶段的方向往90度前进,如图14所示。在图14中,90度标点86为90度方向的候选位置。算法会根据此最大边缘量的位置于两侧的位置,即,最大边缘位置89,将候选坐标,在图14中,90度标点86往血管中心方向调整到血管中心标示点90。最后,此一血管中心标示点90即是此一阶段中,所找出的下一开启搜索点的位置,而下一步的候选位置,是以血管中心标示点90为开启搜索点,继续向135、90,45度方向搜索。
如图15所示。以90度方向一步步往前计算,可以以血管中心标示点90为开启搜索点,逐次找出最大总边缘量的方向,前进的计算标示点91位置。如图16A至图16C所示,以起始搜索的计算标示点91位置,重新上述工字型框架88的搜寻,搜寻方向是向上的90度方向。本发明的算法,内定每一步的可旋转角度是+/-45度。所以图16A至图16C所示的下一步的候选位置为135、90,45度,如图16A至图16C所示的三个位置,包括:45度标点87、90度标点86、以及135度标点85的位置。由左至右分别是往135、90,45度前进的位置。
以工字型框架88为搜寻范围,从血管外部往中心方向计算该角度的最大边缘量。其中,最大边缘量定义是依点选位置沿着工字型框架88的工字型左右两侧的箭头方向,从血管外部往中心方向各扫描若干像素(pixels)的范围,扫描方法先产生邻近像素间灰阶值差异,并将工字型框架88的工字型中轴方向的灰阶值加总,结果得到工字型框架88的工字型左右侧的曲线图,找出左右侧曲线中的坡峰,即差异最大处,为该角度的左侧的最大边缘量(Lmax),以及右侧的最大边缘量(Rmax)。
左右两侧最大边缘量的总和(Lmax+Rmax),是此角度的总边缘量。如图16A实施例中135度的右侧的最大边缘量(Rmax)=473,左侧的最大边缘量(Lmax)=29,因此,135度的总边缘量473+29=502。如图16B实施例中90度的右侧的最大边缘量(Rmax)=701,左侧的最大边缘量(Lmax)=540,90度的总边缘量701+540=1241。如图16C实施例中45度的右侧的最大边缘量(Rmax)=758,左侧的最大边缘量(Lmax)=979,因此,45度的总边缘量758+979=1737。由上述计算,45度总边缘量大于90度总边缘量,90度总边缘量大于135度总边缘量,所以最大总边缘量是45度总边缘量的值,即1737。由于在此阶段的最大总边缘量是45度的值,即1737,所以此阶段的方向往最大总边缘量,即45度前进。
如图16A至图16C所示,工字型框架88为搜寻范围,沿着双箭头方向,从血管外部往双箭头中心方向计算该角度的边缘量。工字型框架88两侧各有一个最大边缘位置89,两侧最大边缘位置89的总和,即是此角度的总边缘量。假设此阶段的最大总边缘量是45度,所以此阶段的方向往45度前进,如图17所示。
在图17中,45度方向的候选位置是45度标点87。假设最大边缘位置89在两侧形成,算法会根据此两侧的最大边缘位置89,将候选坐标的45度标点87往血管中心方向调整到血管中心标示点90。最后血管中心标示点90即是此一阶段所找出的下一点起始搜索点的位置,而下一步的候选位置,是以血管中心标示点90为开启搜索点,继续向135、90,45度方向搜索。
在图18中,是依照上述方法,一步一步往下一个计算出来的位置推进,并且调整行进方向,直到进入终点区域的终点标示点82,即找出此条微血管21的轨迹。
图11中,使用者选取分解画面的白血球22起点,并在后续几张分解画面中选取白血球22终点,依计算机11效能,录像压缩后,可以从数字图像文件案得知图框率(FrameRate),以下以每秒25张的图框率(Frame Rate),举例:每张为1/25秒,分解图的数字影像有9张,头尾间隔8/25秒。图18以及图19中,假设上述路径长计算得86像素(pixel),而显微镜镜头14撷取的数字影像,像素比为1.164594μm/pixel.,第2路径长是86像素(pixel)乘以1.164594μm/pixel=100.155084μm,显微镜镜头14撷取的数字影像,每秒有25张数字影像,由此算出该起点标示点81以及该终点标示点82的时间差为第2时间差,因此第2路径长除以第2时间差,可以计算血流流速,86pixels/(8/25sec)乘以1.164594μm/pixel=294μm/s=0.29mm/s。
图20中,找出该微血管21内的两边缘端点位置,其步骤,包含:由扫描该数字影像成灰阶讯号,形成纵轴灰阶讯号值,横轴像素值;纵轴灰阶讯号加总值最大值,标定该微血管21管径的两边缘端点;该微血管21管径的两边缘端点的相应横轴像素值,计算该微血管21的该管径值;以及于该显示器12显示该微血管21的该管径值。
图20中,依点选位置,左右各扫描16像素(pixel),正方型范围,如果找不到自动再向左右扩大扫描,扫描方法先产生像素间灰阶值差异,结果得到如图20中第1分布图110、以及第2分布图120的像素间灰阶值差异分布图,再将垂直方向的灰阶值加总,结果得到如图20中的第1曲线图111以及第2曲线图121的曲线图。找出曲线中的坡峰,如图中的差异最大值,即为微血管21边缘,依显微镜镜头14,撷取的数字影像,像素比为1.164594μm/pixel,图20中左侧坡峰,在第7像素(pixel)处,右侧在第10像素(pixel)处,因此共间隔16像素(pixel)。计算出微血管21管径,为16像素(pixel)乘以1.164594μm/pixel等于18.6μm。
虽然本发明已以实施例公开如上,然其并非用以限定本发明,任何熟悉此技术者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰,因此本发明的保护范围以权利要求书为准。

Claims (10)

1.一种微血管检测装置,经由一手指皮下组织中的至少一微血管影像,检测所述微血管的血流流速以及管径,其特征在于,包含有:
一计算机,具有一显示器以及一处理器;
一感光耦合组件,电性讯号链接该计算机;以及
一显微镜镜头,经由该显微镜镜头撷取该微血管影像,该微血管影像由该感光耦合组件形成复数帧数字影像,其中时间连续的复数帧该数字影像,经由该处理器显示于该显示器。
2.如权利要求1所述的微血管检测装置,其特征在于,其中,该处理器标定复数帧该数字影像,对应该微血管中一白血球的时间连续的标示点,包括一起点标示点以及一终点标示点,该起点标示点以及该终点标示点的时间差为第1时间差,该处理器计算加总连续标示点的第1路径长,第1路径长除以第1时间差的一血流流速值,于该显示器显示该血流流速值。
3.如权利要求1所述的微血管检测装置,其特征在于,其中,该处理器标定复数帧该数字影像,对应该微血管中一白血球的时间连续的标示点,包括一起点标示点以及一终点标示点,该处理器具有一工字型框架的一搜寻模块,包含有一45度标点、一90度标点以及一135度标点的路经搜寻到一最大边缘位置,该处理器计算一血管中心标示点,由该血管中心标示点,经由该工字型框架的该搜寻模块,依次标示出至少一计算标示点,直到该终点标示点,该起点标示点以及该终点标示点的时间差为第2时间差,该处理器计算加总连续的该计算标示点的第2路径长,除以第2时间差的一血流流速值,于该显示器显示该血流流速值,其中,经由该搜寻模块计算一最大边缘位置的搜寻,是使用一最大边缘量,该最大边缘量依该起点标示点以及该计算标示点的点选位置,沿着该工字型框架的工字型左右两侧的方向,从该微血管外部往该工字型框架的中心方向各扫描复数个像素,扫描产生邻近像素间灰阶值差异,将该工字型框架的工字型中轴方向的灰阶值加总,得到左侧中差异最大处为相应角度的左侧的该最大边缘量,以及得到右侧中差异最大处为相应角度的右侧的该最大边缘量,其中,同一相应角度的左右两侧最大边缘量的总和,分别是该45度标点、该90度标点以及该135度标点的相应角度的一总边缘量,该搜寻模块将该45度标点、该90度标点、以及该135度标点的相应角度分别的该总边缘量中,选出一相应阶段中最大该总边缘量,该相应阶段的下一阶段搜寻的前进方向,是往该相应阶段的最大该总边缘量的相应方向前进。
4.如权利要求1所述的微血管检测装置,其特征在于,其中,该处理器扫描复数帧该数字影像成灰阶讯号,由纵轴灰阶讯号加总值最大值,标定该微血管管径的两边缘端点,该处理器由横轴该微血管的两边缘端点的相应像素值,计算该微血管的一管径值,于该显示器显示该微血管的该管径值。
5.一种使用权利要求1所述的微血管检测装置的微血管检测方法,其特征在于,其中,检测步骤,包含有:
待测一手指放在一指槽中;
使一感光耦合组件取得一数字影像;
经传输将该数字影像显示在一计算机的一显示器;
依该显示器的该数字影像调整该手指的位置及一显微镜镜头的焦距;以及
撷取该数字影像并依连续时间分解该数字影像。
6.如权利要求5所述的微血管检测方法,分解该数字影像作为量测一微血管的一血流流速值,其特征在于,其中,检测步骤,包含有:
点选该数字影像中一白血球起始位置;
找出时间连续的该数字影像中同一个该白血球的位置,并点选该白血球的位置;
搜寻到路径,则标示路径并计算该血流流速值;以及
搜寻不到路径,则显示错误,并重新回到点选该数字影像中该白血球起始位置的步骤。
7.如权利要求5所述的微血管检测方法,分解该数字影像作为量测一微血管的一管径值,其特征在于,其中,检测步骤,包含有:
点选该数字影像中该微血管内的任一位置;
找出该微血管内的两边缘端点位置,并量测该微血管的管径;以及
手动微调量测位置,并显示该微血管的该管径值于该显示器。
8.如权利要求6所述的微血管检测方法,分解该数字影像作为量测一微血管的一血流流速值,其特征在于,其中,
同一个该白血球在时间连续的该数字影像中的位置,被标定为一起点标示点以及一终点标示点、以及至少一追踪点;
计算所有该追踪点的一平均位置点,各追踪点和平均位置的角度,依角度进行排序;
得到该白血球流经的顺序,该起点标示点以及该终点标示点的时间差为第1时间差;
依序加总相邻距离,得到流经该微血管的第1路径长;
第1路径长除以第1时间差的一血流流速值;以及
于该显示器显示该血流流速值。
9.如权利要求6所述的微血管检测方法,分解该数字影像作为量测一微血管的一血流流速值,其特征在于,其中,检测步骤,更包含有:
点选该数字影像中该白血球的一起点标示点以及一终点标示点;
使用一工字型框架,以45度标点、90度标点、以及135度标点的路经搜寻到一最大边缘位置;
计算一血管中心标示点;
由该血管中心标示点,经由该工字型框架依次标示出至少一计算标示点,直到该终点标示点;
该起点标示点以及该终点标示点的时间差为第2时间差,计算加总连续的该计算标示点的第2路径长,除以第2时间差的一血流流速值;以及
于该显示器显示该血流流速值;
其中,该最大边缘位置的搜寻步骤,包括:
使用一最大边缘量,该最大边缘量依该起点标示点以及该计算标示点的点选位置,沿着该工字型框架的工字型左右两侧的方向,从该微血管外部往该工字型框架的中心方向各扫描复数个像素;
扫描产生邻近像素间灰阶值差异;
将该工字型框架的工字型中轴方向的灰阶值加总;
得到左侧中差异最大处为相应角度的左侧的该最大边缘量,以及得到右侧中差异最大处为相应角度的右侧的该最大边缘量;
将同一相应角度的左右两侧最大边缘量的总和,分别是该45度标点、该90度标点以及该135度标点的相应角度的一总边缘量;
将该45度标点、该90度标点以及该135度标点的相应角度分别的该总边缘量中,选出一相应阶段中最大该总边缘量;以及
使该相应阶段的下一阶段搜寻的前进方向,是往该相应阶段的最大该总边缘量的相应方向前进。
10.如权利要求7所述的微血管检测方法,其特征在于,其中,找出该微血管内的两边缘端点位置,其步骤,包含:
由扫描该数字影像成灰阶讯号,形成纵轴灰阶讯号值,横轴像素值;
纵轴灰阶讯号加总值最大值,标定该微血管管径的两边缘端点;
该微血管管径的两边缘端点的相应横轴像素值,计算该微血管的该管径值;以及
于该显示器显示该微血管的该管径值。
CN201810222949.4A 2018-03-19 2018-03-19 微血管检测装置和方法 Active CN110279393B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810222949.4A CN110279393B (zh) 2018-03-19 2018-03-19 微血管检测装置和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810222949.4A CN110279393B (zh) 2018-03-19 2018-03-19 微血管检测装置和方法

Publications (2)

Publication Number Publication Date
CN110279393A true CN110279393A (zh) 2019-09-27
CN110279393B CN110279393B (zh) 2023-04-28

Family

ID=68000787

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810222949.4A Active CN110279393B (zh) 2018-03-19 2018-03-19 微血管检测装置和方法

Country Status (1)

Country Link
CN (1) CN110279393B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111956200A (zh) * 2020-09-01 2020-11-20 徐州利华电子科技发展有限公司 微循环高速血流测算***和测算方法
CN117752298A (zh) * 2023-12-29 2024-03-26 山东大学 基于octa的毛细血管血流流速测量方法及***

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN85201718U (zh) * 1985-04-28 1986-05-28 清华大学 微循环多参数综合测试仪
JPH0765154A (ja) * 1993-08-31 1995-03-10 Toshiba Corp 血管像の定量解析装置及びその定量解析方法
TW200618774A (en) * 2004-12-01 2006-06-16 Yio-Wha Shau Real time inspection system and inspection method for micro-circulation
JP2008017961A (ja) * 2006-07-12 2008-01-31 Image One Co Ltd 血流速度の測定方法及び装置
JP2010075354A (ja) * 2008-09-25 2010-04-08 Takenori Uchiyama 毛細血管血流測定装置、毛細血管血流測定方法およびプログラム
CN101779967A (zh) * 2009-01-21 2010-07-21 重庆医科大学超声影像学研究所 应用超声微泡(微球)造影剂直接测量血流速度
CN101984916A (zh) * 2010-11-17 2011-03-16 哈尔滨工程大学 基于数字图像处理技术的血管管径测量方法
TW201309264A (zh) * 2011-08-26 2013-03-01 Univ Nat Sun Yat Sen 自動化甲襞微循環分析系統

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN85201718U (zh) * 1985-04-28 1986-05-28 清华大学 微循环多参数综合测试仪
JPH0765154A (ja) * 1993-08-31 1995-03-10 Toshiba Corp 血管像の定量解析装置及びその定量解析方法
TW200618774A (en) * 2004-12-01 2006-06-16 Yio-Wha Shau Real time inspection system and inspection method for micro-circulation
JP2008017961A (ja) * 2006-07-12 2008-01-31 Image One Co Ltd 血流速度の測定方法及び装置
JP2010075354A (ja) * 2008-09-25 2010-04-08 Takenori Uchiyama 毛細血管血流測定装置、毛細血管血流測定方法およびプログラム
CN101779967A (zh) * 2009-01-21 2010-07-21 重庆医科大学超声影像学研究所 应用超声微泡(微球)造影剂直接测量血流速度
CN101984916A (zh) * 2010-11-17 2011-03-16 哈尔滨工程大学 基于数字图像处理技术的血管管径测量方法
TW201309264A (zh) * 2011-08-26 2013-03-01 Univ Nat Sun Yat Sen 自動化甲襞微循環分析系統

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111956200A (zh) * 2020-09-01 2020-11-20 徐州利华电子科技发展有限公司 微循环高速血流测算***和测算方法
CN117752298A (zh) * 2023-12-29 2024-03-26 山东大学 基于octa的毛细血管血流流速测量方法及***
CN117752298B (zh) * 2023-12-29 2024-05-24 山东大学 基于octa的毛细血管血流流速测量方法及***

Also Published As

Publication number Publication date
CN110279393B (zh) 2023-04-28

Similar Documents

Publication Publication Date Title
CN102436652B (zh) 一种多源遥感图像自动配准方法
CN103646250B (zh) 基于距离图像头肩特征的行人监测方法及装置
CN113920097B (zh) 一种基于多源图像的电力设备状态检测方法及***
CN102341810B (zh) 用于小瓶和容器上的帽体的存在和类型的自动化探测的方法和设备
CN116721107B (zh) 一种电缆生产质量智能监控***
WO2022262054A1 (zh) 图像处理方法、装置、设备及存储介质
US20230098480A1 (en) Apparatus and method for detecting intersection edges
CN103090796B (zh) 火箭偏摆、沉降的测量***及方法
CN110279393A (zh) 微血管检测装置和方法
CN107092905B (zh) 一种电力巡检机器人的待识别仪器定位方法
CN107123126A (zh) 一种人流运动场景热度估计方法
CN107067430A (zh) 一种基于特征点聚类的麦田作物行检测方法
TWI796317B (zh) 微血管檢測裝置和方法
TWM563247U (zh) 微血管檢測裝置
CN103337080A (zh) 基于梯度方向Hausdorff 距离的红外图像和可见光图像配准技术
CN112446304A (zh) 一种火焰检测方法及其***
CN112132823A (zh) 一种基于人工智能的光伏电池图像局部光照不均处理方法
CN103776376A (zh) 一种激光定位灯十字线图像检测方法及装置
RU2538340C1 (ru) Способ совмещения изображений, полученных с помощью разнодиапазонных фотодатчиков
CN107240107B (zh) 一种基于图像检索的显著性检测元评估方法
CN113534296B (zh) 基于神经网络沙尘天气预报强度误差的测算方法及装置
CN112255141B (zh) 一种热成像气体监测***
CN116030373A (zh) 鱼道智能监测***
CN106996934A (zh) 多相机视觉检测***及检测方法
TWI689285B (zh) 臉部對稱性檢測方法及其系統

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant