CN105769117B - 一种眼底光学全场微血管造影成像装置及方法 - Google Patents
一种眼底光学全场微血管造影成像装置及方法 Download PDFInfo
- Publication number
- CN105769117B CN105769117B CN201610017830.4A CN201610017830A CN105769117B CN 105769117 B CN105769117 B CN 105769117B CN 201610017830 A CN201610017830 A CN 201610017830A CN 105769117 B CN105769117 B CN 105769117B
- Authority
- CN
- China
- Prior art keywords
- signal
- mrow
- speckle
- frequency
- imaging
- 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.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
- A61B3/1241—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes specially adapted for observation of ocular blood flow, e.g. by fluorescein angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0071—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by measuring fluorescence emission
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/026—Measuring blood flow
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Hematology (AREA)
- Pathology (AREA)
- Ophthalmology & Optometry (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Vascular Medicine (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明所述一种眼底光学全场微血管造影成像装置及方法，包括眼底照明光路、眼底成像光路、同步触发源和数据处理系统，其中眼底照明光路和眼底成像光路由同步触发源同步触发；眼底照明光路包括依次连接的光源控制器、激光光源、扩束镜、分光镜和接目物镜，所述眼底成像光路包括依次连接的接目物镜、分光镜、成像目镜、镜头和COMS相机。本发明使光源以与相机采样帧率相同的频率发出激光，在保证实现眼底微血管成像的同时减小激光功率使之达到人眼所能承受的范围。本发明可以实现眼底光学全场微血管造影成像，该方法无需造影剂，具有较高的分辨率、信噪比和探测深度，可以重建出眼底的毛细血管。
Description
技术领域
本发明涉及血管造影成像技术领域，特别是涉及一种眼底光学全场微血管造影成像装置及方法。
背景技术
血流成像在医学临床诊断与病理基础研究中起到了重要作用，其中微血管病变是糖尿病伴有的诸多并发症中较为普遍的一种，它会使多脏器发生病变，包括糖尿病性肾病、糖尿病性神经系统病变和糖尿病性视网膜病变。临床医学中，对于糖尿病性视网膜微血管病变的诊断通常采用眼底荧光血管造影的方法实现，该方法需要将荧光素钠光学造影剂注射到患者眼中，可能造成人体的不良反应。
应用激光散斑原理进行全场微血管造影成像是目前世界各地的研究热点，包括激光多普勒频移和激光散斑衬比成像等方法。激光多普勒频移方法已经实现了皮肤烧伤修复以及手指血流灌注的检测，但是该方法中散斑多普勒频移光子被高速面阵相机捕捉，需要较高的采样速度，因而会降低像素分辨率，故空间分辨率较低。激光散斑衬比成像方法以衬比值作为成像参数，已经实现了活体皮肤、大脑皮层等组织的微血管造影成像，可作为视网膜成像的手段，但其理论是基于散斑的时间与空间自相关特性，在一定环境下散斑衬比反比于血液流速的平方根，其抑制背景噪声的能力不强，成像质量有待改善。同时，上述方法以激光作为光源，应用于眼底成像时需使用激光直接照射视网膜，因此对于激光功率要求十分严格，若功率过大，可能对人体造成损伤；若功率过小，则会损失成像质量。
发明内容
本发明克服了现有技术中的缺点，提供了一种眼底光学全场微血管造影成像装置及方法，可实现对眼底光学全场微血管造影成像，该方法不需要造影剂，对人体无害，具有较高的空间分辨率。
为了解决上述技术问题，本发明是通过以下技术方案实现的：
一种眼底光学全场微血管造影成像装置，包括眼底照明光路(1)、眼底成像光路(2)、同步触发源(10)和数据处理系统(20)，其中眼底照明光路(1)和眼底成像光路(2)由同步触发源(10)同步触发；所述眼底照明光路(1)包括依次连接的光源控制器(11)、激光光源(12)、扩束镜(13)、分光镜(16)和接目物镜(15)，所述眼底成像光路(2)包括依次连接的接目物镜(15)、分光镜(16)、成像目镜(17)、镜头(18)和COMS相机(19)。
一种眼底光学全场微血管造影成像的方法，光源控制器(11)和所述COMS相机(19)受同步触发源(10)发出的同步触发信号控制，一旦光源控制器(11)接收到同步触发信号，则触发激光光源(12)发出激光；一旦COMS相机19接收到同步触发信号，则受到触发进行采集，扩束镜(13)用于对激光扩束，靠近人眼的接目物镜(15)和靠近COMS相机(19)的成像目镜(17)同作为摄影物镜，用于控制成像焦距，在眼底照明光路(1)中，分光镜(16)用于反射激光到达眼底；在眼底成像光路(2)中，分光镜(16)用于透射由眼底反射的散斑信号，以便相机采集图像。
其包括以下步骤：
步骤1、设置COMS相机的镜头参数，选择短曝光的相机工作模式；
设置相机的镜头参数，旨在保证红细胞、散斑、像素单元的成像面积接近相等，当红细胞进、出像素单元的成像范围时，局域的浓度涨落产生干涉强度变化的动态散斑。
选择短曝光的相机工作模式，旨在保证相机的曝光时间足够短，以捕捉频率较高的动态散斑信号的瞬时状态，而非散斑信号强度的变化积分值。
步骤2、同步触发源以频率F0(100/s≤F0≤300/s)发出n(n为大于或等于10的自然数)次同步触发信号，分别作用于眼底照明光路(1)和眼底成像光路(2)；
①眼底照明光路(1)中光源控制器接收到所述频率为F0、n次发出的同步触发信号，n次触发激光光源以频率F0出射激光，所述n次出射的激光经扩束后由分光镜反射，产生的n次反射光以所述频率F0通过接目物镜照射人眼；
②眼底成像光路(2)中，由人眼以所述频率F0、n次反射的激光散斑信号经由分光镜透射到成像目镜，COMS相机接收到所述同步触发源以频率F0发出的n个同步触发信号，n次连续采集以所述频率F0透过成像目镜的n张散斑图像，每张散斑图像上有m(m为面阵相机的像素单元个数)个散斑光强点，所述F0同时为采样帧率；
步骤3、数据处理系统中，由计算机对所述相机n次采集到的m个原始时间序列信号进行FFT变换以获得m个散斑频域信号，设定滤波频率u0，对于第k(1≤k≤m)个所述散斑频域信号，将背景噪声和集中在零频以及低频范围内的静态散斑频域信号与集中在高频范围内的动态散斑频域信号分离，计算所述第k个点(坐标为(x,y))的调制深度MD(x,y)，以全部m个点的调制深度MD(x,y)作为成像参数进行成像。
所述同步触发源的触发频率F0在100/s至300/s的范围内，旨在保证相机处于长采样的工作模式，因此保证相机拥有足够的像素单元，当干涉光信号的周期远小于采样时间时，相机的像素单元数量足够用来记录散斑信号，得到具有高分辨率的图像。
所述步骤2中，利用COMS相机采集样本的后向散射光，从而进行散斑图像的重建。
所述步骤2中，COMS相机捕捉的散斑强度信号是不同干涉项频率的强度总和，对于COMS相机第t(1≤t≤n)次采集到的所述m个原始时间序列信号，第k(1≤k≤m)个像素单元(坐标为(x,y))散斑强度信号表示为：
其中，I0表示静态散斑的信号强度，该强度始终保持不变，对应于背景信号。表示移动散射物质之间以及静止散射物质与移动散射物质之间相互干涉产生的干涉信号强度，其中Ii是第i个干涉项的强度，fi是相对应的时间调制频率，∏表示方脉冲函数，定义为∏T/2[t-t0]＝H[t-t0+T/2]-H[t-t0-T/2]，用来描述信号的采集过程，t0表示某一时间点，T为总采集时间。
所述步骤3中，通过进行快速傅里叶变换将权利要求4中所述相机n次采集到的m个原始时间序列信号进行快速傅里叶变换将相机接收的时域散斑强度信号IP沿着时间序列转变为频域散斑强度信号，其中低频范围主要为静态散斑信号，高频范围主要为动态散斑信号，相机第t次采集到的所述信号的散斑频谱图可以表示为：
FTt→u[IP(t)]＝iS[u]+id[u±fi]
其中，iS[u]表示静态散斑信号频谱，即以零频为中心的小范围内的频谱，其表达式为iS[u]＝I0sinc(Tu)；id[u±fi]表示动态散斑频谱，即频域中除去静态散斑信号频谱的剩余部分，其表达式为
坐标为(x,y)的像素点位置重建图像的时域的调制深度表示为：
其中，<id(x,y,u)>是坐标为(x,y)的像素点位置对应信号经滤波后的动态散斑信号对于时间的平均值，<is(x,y,u)>是坐标为(x,y)的像素点位置对应信号经滤波后的静态散斑信号对于时间的平均值。
坐标为(x,y)的像素点位置重建图像的频域的调制深度表示为：
其中，id(x,y,u)是坐标为(x,y)的像素点位置对应信号经滤波后的动态散斑信号；is(x,y,u)是坐标为(x,y)的像素点位置对应信号经滤波后的静态散斑信号。
所述步骤3中，设置合适的滤波频率u0，以分开静态散斑信号和动态散斑信号。所述静态散斑信号为背景信号，反映了生物的背景组织；所述动态散斑信号为血流信号。
所述步骤3中，调制深度MD(x,y)为成像参数，可以分为时域的调制深度和频域的调制深度。所述时域调制深度定义为经滤波后动态散斑信号对于时间的平均值与经滤波后静态散斑信号对于时间的平均值的比值；所述频域调制深度定义为直接在频域进行计算的调制深度，即动态散斑数据之和与静态散斑数据之和的比值。以频域的调制深度作为成像参量，无需对于频域数据再进行逆傅立叶变换，可以缩减数据处理的时间，同时降低对于计算机内存大小的要求。
与现有技术相比，本发明的有益效果是：
本发明所述一种眼底光学全场微血管造影成像装置及方法，可以用于眼底光学全场微血管造影成像，以激光作为光源，采用短曝光和长采样的相机工作模式捕捉样品的后向散射光，以调制深度为成像参数可以清晰的还原毛细血管图像，具有信噪比高、空分辨率高、探测深度高等优点，除了对于微血管进行全场造影成像，该方法还可以实时监测血流的动态过程。
本发明提供的方法和装置无需造影剂，本发明使用同步触发装置控制相机和光源同步运行，可以在保证成像质量的前提下减小激光功率，达到医用标准，无论用于临床检测还是生物实验，对于人体和小动物都无伤害。
本发明使光源以与相机采样帧率相同的频率发出激光，在保证实现眼底微血管成像的同时减小激光功率使之达到人眼所能承受的范围。本发明可以实现眼底光学全场微血管造影成像，该方法无需造影剂，具有较高的分辨率、信噪比和探测深度，可以重建出眼底的毛细血管。
附图说明
附图用来提供对本发明的进一步理解，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制，在附图中：
图1为本发明所述一种眼底光学全场微血管造影成像装置原理图。
图2为本发明一种眼底光学全场微血管造影成像方法流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
本发明中以激光作为光源，以调制深度作为成像参量，可实现具有较高空间分辨率的眼底光学全场微血管造影成像。为使本发明更加清晰，以下参照附图，结合实施例，对本发明作进一步详细说明。
本发明提供了一种眼底光学全场微血管造影成像装置及方法，以调制深度作为成像参量，可实现具有较高空间分辨率的全场微血管造影成像。本实施例提供一种对于人眼底光学全场微血管造影成像的方法及装置，图1是本发明的装置示意图，本实施例采用图1所示装置，包括：同步触发源10、眼底照明光路1、眼底成像光路2和数据处理系统20，其中眼底照明光路1和眼底成像光路2由同步触发源10同步触发。所述眼底照明光路包括光源控制器11、激光光源12、扩束镜13、分光镜16和接目物镜15，所述眼底成像光路包括接目物镜15、分光镜16、成像目镜17、镜头18和COMS相机19。
所述光源控制器受同步触发源10发出的同步触发信号控制，一旦光源控制器11接收到所述同步触发信号，则触发激光光源12发出激光，一旦所述COMS相机19接收到所述同步触发信号，则受到触发进行采集。扩束镜13用于对激光扩束，靠近人眼14的接目物镜15和靠近相机的成像目镜17同作为摄影物镜，用于控制成像焦距，从而使本发明所述方法适用于各类情况的人眼14(包括可能存在光焦度不正的病眼)。接目物镜15和分光镜16既是眼底照明光路中的一部分，也是眼底成像光路的一部分，在眼底照明光路中，分光镜16用于反射激光到达人眼14；在眼底成像光路中，分光镜16用于透射由人眼14反射的散斑信号，以便相机采集图像。
具体地，图2为本发明所述方法的流程图，依照图2流程，设置相机的镜头参数，选择短曝光的相机工作模式。
所述同步触发源以频率F0(100/s≤F0≤300/s)发出n(n为大于或等于10的自然数)次同步触发信号，分别作用于眼底照明光路1和眼底成像光路2。
所述眼底照明光路1中光源控制器接收到所述频率为F0、n次发出的同步触发信号，n次触发激光光源以频率F0出射激光，所述n次出射的激光经扩束后由分光镜反射，产生的n次反射光以所述频率F0通过接目物镜照射人眼。
所述眼底成像光路2中，由人眼以所述频率F0、n次反射的激光散斑信号经由分光镜透射到成像目镜，COMS相机接收到所述同步触发源以频率F0发出的n个同步触发信号，n次连续采集以所述频率F0透过成像目镜的n张散斑图像，每张散斑图像上有m(m为面阵COMS相机的像素单元个数)个散斑光强点，所述F0同时为采样帧率。
所述数据处理系统中，对于相机第t(1≤t≤n)次采集到的所述m个原始时间序列信号，第k(1≤k≤m)个像素单元(坐标为(x,y))散斑强度信号表示为：
其中，I0表示静态散斑的信号强度，该强度始终保持不变，对应于背景信号。表示移动散射物质之间以及静止散射物质与移动散射物质之间相互干涉产生的干涉信号强度，其中Ii是第i个干涉项的强度，fi是相对应的时间调制频率，∏表示方脉冲函数，定义为∏T/2[t-t0]＝H[t-t0+T/2]-H[t-t0-T/2]，用来描述信号的采集过程，t0表示某一时间点，T为总采集时间。
对所述相机n次采集到的m个原始时间序列信号进行快速傅里叶变换将相机接收的时域散斑强度信号IP沿着时间序列转变为频域散斑强度信号，其中低频范围主要为静态散斑信号，高频范围主要为动态散斑信号，相机第t次采集到的所述信号的散斑频谱图可以表示为：
FTt→u[IP(t)]＝iS[u]+id[u±fi]
其中，iS[u]表示静态散斑信号频谱，即以零频为中心的小范围内的频谱，其表达式为iS[u]＝I0 sin c(Tu)；id[u±fi]表示动态散斑频谱，即频域中除去静态散斑信号频谱的剩余部分，其表达式为
对于第k个所述散斑频域信号，通过设置适合的滤波范围u0，可以将反映组织背景信息的静态散斑信号和反映血液流动的动态散斑信号分离，计算所述第k个点的调制深度MD(x,y)，以所述m个点的调制深度MD(x,y)作为成像参数进行成像。通过上述方法可以提取出血液流动的信息，即可实现全场微血管造影成像。
所述调制深度MD(x,y)，其在时域的表达式为
其中，MD(x,y)是坐标为(x,y)的像素点位置重建图像的调制深度值，<id(x,y,u)>是坐标为(x,y)的像素点位置对应信号经滤波后的动态散斑信号对于时间的平均值，<is(x,y,u)>是坐标为(x,y)的像素点位置对应信号经滤波后的静态散斑信号对于时间的平均值。
进一步地，所述MD(x,y)作为成像参数，直接在频域进行计算对成像结果没有影响，故可使用频域的调制深度作为参数进行成像，其表达式为
其中MD(x,y)是坐标为(x,y)的像素点位置重建图像的频域调制深度值，定义为动态散斑数据之和与静态散斑数据之和的比值；id(x,y,u)是坐标为(x,y)的像素点位置对应信号经滤波后的动态散斑信号；is(x,y,u)是坐标为(x,y)的像素点位置对应信号经滤波后的静态散斑信号。以频域的调制深度作为成像参量，即对于频域数据不再进行逆傅立叶变换，可以缩减数据处理的时间，同时降低对于计算机内存大小的要求。
本领域普通技术人员可以理解：上述实施例所述步骤中部分或全部可通过程序控制相关硬件完成，可将所述程序储存于计算机可读取介质中，该程序执行包括本发明所述方法以及上述实施例的步骤。前述储存介质包括RAM、ROM、光盘、磁盘等各种可用于储存程序代码的储存介质。
最后应说明的是：以上仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，但是凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种眼底光学全场微血管造影成像的方法，其特征在于：包括以下步骤：
步骤1、设置COMS相机的镜头参数，选择短曝光的相机工作模式；
步骤2、同步触发源以频率F0发出n次同步触发信号，其中，100/s≤F0≤300/s，n为大于或等于10的自然数，分别作用于眼底照明光路(1)和眼底成像光路(2)；
①眼底照明光路(1)中光源控制器接收到所述频率为F0、n次发出的同步触发信号，n次触发激光光源以频率F0出射激光，所述n次出射的激光经扩束后由分光镜反射，产生的n次反射光以所述频率F0通过接目物镜照射人眼；
②眼底成像光路(2)中，由人眼以所述频率F0、n次反射的激光散斑信号经由分光镜透射到成像目镜，COMS相机接收到所述同步触发源以频率F0发出的n个同步触发信号，n次连续采集以所述频率F0透过成像目镜的n张散斑图像，每张散斑图像上有m个散斑光强点，所述F0同时为采样帧率，m为面阵相机的像素单元个数；
步骤3、数据处理系统中，由计算机对所述相机第n次采集到的m个散斑光强点的原始时间序列信号进行FFT变换以获得m个散斑频域信号，设定滤波频率u0，对于第k个所述散斑频域信号，其中1≤k≤m，将背景噪声和集中在零频以及低频范围内的静态散斑频域信号与集中在高频范围内的动态散斑频域信号分离，计算所述第k个点的调制深度MD(x,y)，第k个点的坐标为(x,y)，以全部m个点的调制深度MD(x,y)作为成像参数进行成像；
利用这个方法的成像装置包括眼底照明光路(1)、眼底成像光路(2)、同步触发源(10)和数据处理系统(20)，其中眼底照明光路(1)和眼底成像光路(2)由同步触发源(10)同步触发；所述眼底照明光路(1)包括依次连接的光源控制器(11)、激光光源(12)、扩束镜(13)、分光镜(16)和接目物镜(15)，所述眼底成像光路(2)包括依次连接的接目物镜(15)、分光镜(16)、成像目镜(17)、镜头(18)和COMS相机(19)；
光源控制器(11)和COMS相机(19)受同步触发源(10)发出的同步触发信号控制，一旦光源控制器(11)接收到同步触发信号，则触发激光光源(12)发出激光；一旦COMS相机(19)接收到同步触发信号，则受到触发进行采集，扩束镜(13)用于对激光扩束，靠近人眼的接目物镜(15)和靠近COMS相机(19)的成像目镜(17)同作为摄影物镜，用于控制成像焦距，在眼底照明光路(1)中，分光镜(16)用于反射激光到达眼底；在眼底成像光路(2)中，分光镜(16)用于透射由眼底反射的散斑信号，以便相机采集图像；
所述步骤2中，COMS相机捕捉的散斑强度信号是不同干涉项频率的强度总和，对于COMS相机第t次采集到的所述m个散斑光强点的原始时间序列信号，其中，1≤t≤n，第k个像素单元的散斑强度信号表示为，其中，1≤k≤m，像素单元坐标为(x,y)，
其中，I0表示静态散斑的信号强度，该强度始终保持不变，对应于背景信号，表示移动散射物质之间以及静止散射物质与移动散射物质之间相互干涉产生的干涉信号强度，其中Ii是第i个干涉项的强度，fi是相对应的时间调制频率，∏表示方脉冲函数，定义为∏T/2[t-t0]＝H[t-t0+T/2]-H[t-t0-T/2]，用来描述信号的采集过程，t0表示某一时间点，T为总采集时间。
2.根据权利要求1所述一种眼底光学全场微血管造影成像的方法，其特征在于：所述步骤2中，利用COMS相机采集样本的后向散射光，从而进行散斑图像的重建。
3.根据权利要求1所述一种眼底光学全场微血管造影成像的方法，其特征在于：所述步骤3中，通过所述相机第n次采集到的m个原始时间序列信号进行快速傅里叶变换，将相机接收的时域散斑强度信号IP沿着时间序列转变为频域散斑强度信号，其中低频范围主要为静态散斑信号，高频范围主要为动态散斑信号，相机第t次采集到的所述信号的散斑频谱图可以表示为：
FTt→u[IP(t)]＝iS[u]+id[u±fi]，
其中，iS[u]表示静态散斑信号频谱，即以零频为中心的小范围内的频谱，其表达式为iS[u]＝I0sinc(Tu)；id[u±fi]表示动态散斑频谱，即频域中除去静态散斑信号频谱的剩余部分，其表达式为
4.根据权利要求1所述一种眼底光学全场微血管造影成像的方法，其特征在于：坐标为(x,y)的像素点位置重建图像的时域的调制深度表示为：
<mrow>
<mi>M</mi>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo><</mo>
<msub>
<mi>i</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
</mrow>
<mrow>
<mo><</mo>
<msub>
<mi>i</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中，<id(x,y,u)>是坐标为(x,y)的像素点位置对应信号经滤波后的动态散斑信号对于时间的平均值，<is(x,y,u)>是坐标为(x,y)的像素点位置对应信号经滤波后的静态散斑信号对于时间的平均值。
5.根据权利要求1所述一种眼底光学全场微血管造影成像的方法，其特征在于：坐标为(x,y)的像素点位置重建图像的频域的调制深度表示为：
<mrow>
<mi>M</mi>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
<mrow>
<mi>f</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msub>
<mi>i</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>u</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>u</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
</msubsup>
<msub>
<mi>i</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>u</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>u</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中，id(x,y,u)是坐标为(x,y)的像素点位置对应信号经滤波后的动态散斑信号；is(x,y,u)是坐标为(x,y)的像素点位置对应信号经滤波后的静态散斑信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610017830.4A CN105769117B (zh) | 2016-01-11 | 2016-01-11 | 一种眼底光学全场微血管造影成像装置及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610017830.4A CN105769117B (zh) | 2016-01-11 | 2016-01-11 | 一种眼底光学全场微血管造影成像装置及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105769117A CN105769117A (zh) | 2016-07-20 |
| CN105769117B true CN105769117B (zh) | 2017-08-29 |
Family
ID=56403123
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610017830.4A Active CN105769117B (zh) | 2016-01-11 | 2016-01-11 | 一种眼底光学全场微血管造影成像装置及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105769117B (zh) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106419889B (zh) * | 2016-10-31 | 2023-12-05 | 佛山科学技术学院 | 一种基于片层光的血流三维成像装置及方法 |
| CN106691410B (zh) * | 2016-12-13 | 2024-01-16 | 佛山科学技术学院 | 一种脉搏及血红细胞浓度监测仪及方法 |
| CN106580328B (zh) * | 2016-12-22 | 2024-01-12 | 佛山科学技术学院 | 基于血红细胞浓度测量的烧伤程度与面积检测系统及方法 |
| CN107468234B (zh) * | 2017-08-28 | 2023-08-08 | 佛山科学技术学院 | 一种血管血流投影层析三维成像系统及成像方法 |
| EP4293484A3 (en) * | 2017-09-28 | 2024-03-13 | Apple Inc. | Method and device for eye tracking using event camera data |
| CN107907518A (zh) * | 2017-12-12 | 2018-04-13 | 佛山科学技术学院 | 一种在高散射介质中的动态荧光成像装置及方法 |
| CN107907519A (zh) * | 2017-12-12 | 2018-04-13 | 佛山科学技术学院 | 一种荧光分子信号的成像方法及其装置 |
| CN109124615B (zh) * | 2018-09-06 | 2023-07-11 | 佛山科学技术学院 | 一种可选区高动态激光散斑血流成像装置及方法 |
| CN109124686A (zh) * | 2018-09-11 | 2019-01-04 | 中山大学 | 一种血液流量测量设备及方法 |
| CN109394181A (zh) * | 2018-12-05 | 2019-03-01 | 吉林大学 | 一种脑部功能区域定位系统、方法以及可移动设备 |
| CN109431458A (zh) * | 2018-12-21 | 2019-03-08 | 合肥奥比斯科技有限公司 | 多光谱光源及眼底成像系统 |
| CN111067479B (zh) * | 2019-12-31 | 2021-08-13 | 西安电子科技大学 | 一种眼底成像装置及成像方法 |
| CN111528793B (zh) * | 2020-05-12 | 2022-06-21 | 宁波蓝明信息科技有限公司 | 一种视网膜血管光学造影成像系统及方法 |
| CN112229822B (zh) * | 2020-08-25 | 2022-05-13 | 西安电子科技大学 | 对流动液体内多目标的反射式单帧散射成像装置及方法 |
| CN112712521B (zh) * | 2021-01-18 | 2023-12-12 | 佛山科学技术学院 | 基于全局梯度搜索的眼底视盘自动定位方法及其存储介质 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6215957B1 (en) * | 1999-12-03 | 2001-04-10 | Peter Arthur Van Houten | Synchronizer for fundus camera |
| CN1498593A (zh) * | 2002-11-11 | 2004-05-26 | ������������ʽ���� | 眼科摄影装置 |
| CN201624637U (zh) * | 2009-08-21 | 2010-11-10 | 上海市第一人民医院 | 可调谐光源结合共焦激光同步眼底血管造影系统 |
| CN101926640A (zh) * | 2009-06-22 | 2010-12-29 | 中国科学院光电技术研究所 | 活体人眼视网膜动态成像仪 |
| CN102008288A (zh) * | 2010-12-17 | 2011-04-13 | 中国科学院光电技术研究所 | 一种线扫描共焦检眼镜的系统和方法 |
| CN102429636A (zh) * | 2011-09-30 | 2012-05-02 | 中国科学院长春光学精密机械与物理研究所 | 大视场液晶自适应光学眼底成像的系统 |
| CN105286779A (zh) * | 2015-10-30 | 2016-02-03 | 温州医科大学 | 一种在体视网膜血流动力学的成像与绝对流速测量方法 |
-
2016
- 2016-01-11 CN CN201610017830.4A patent/CN105769117B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6215957B1 (en) * | 1999-12-03 | 2001-04-10 | Peter Arthur Van Houten | Synchronizer for fundus camera |
| CN1498593A (zh) * | 2002-11-11 | 2004-05-26 | ������������ʽ���� | 眼科摄影装置 |
| CN101926640A (zh) * | 2009-06-22 | 2010-12-29 | 中国科学院光电技术研究所 | 活体人眼视网膜动态成像仪 |
| CN201624637U (zh) * | 2009-08-21 | 2010-11-10 | 上海市第一人民医院 | 可调谐光源结合共焦激光同步眼底血管造影系统 |
| CN102008288A (zh) * | 2010-12-17 | 2011-04-13 | 中国科学院光电技术研究所 | 一种线扫描共焦检眼镜的系统和方法 |
| CN102429636A (zh) * | 2011-09-30 | 2012-05-02 | 中国科学院长春光学精密机械与物理研究所 | 大视场液晶自适应光学眼底成像的系统 |
| CN105286779A (zh) * | 2015-10-30 | 2016-02-03 | 温州医科大学 | 一种在体视网膜血流动力学的成像与绝对流速测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105769117A (zh) | 2016-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105769117B (zh) | 一种眼底光学全场微血管造影成像装置及方法 | |
| JP4621496B2 (ja) | ラインスキャン検眼鏡 | |
| CN106166058B (zh) | 一种应用于光学相干断层扫描血管成像方法及oct系统 | |
| JP2013059551A (ja) | 眼底観察装置 | |
| CN103745257B (zh) | 一种基于图像识别的锥细胞密度的计算方法 | |
| CN105231990B (zh) | 基于oct三维成像分析角膜生物力学性能的装置及方法 | |
| CA3132140C (en) | Intrinsic contrast optical cross-correlated wavelet angiography | |
| JP6923392B2 (ja) | 眼科装置 | |
| EP3801208B1 (en) | Optical coherence imager and method | |
| US10909683B2 (en) | Methods and systems for enhancing optical image quality | |
| CN105249933A (zh) | 光声分子三维成像仪 | |
| Choi et al. | Volumetric cutaneous microangiography of human skin in vivo by VCSEL swept-source optical coherence tomography | |
| CN111067479B (zh) | 一种眼底成像装置及成像方法 | |
| JP2003532461A (ja) | 網膜像を用いた非観血的測定方法 | |
| CN206261606U (zh) | 一种体外鼻咽腔体静脉血氧饱和度成像装置 | |
| Saldan et al. | Efficiency of optical-electronic systems: methods application for the analysis of structural changes in the process of eye grounds diagnosis | |
| Zhang et al. | A low-cost and smartphone-based laser speckle contrast imager for blood flow | |
| Samanta et al. | Ocular imaging of Carassius auratus eyes using photoacoustic microscopy | |
| Gologorsky et al. | Principles of Ocular Imaging | |
| Bracco et al. | BIE-PInCS: Brain injury evaluation with pupillometer based on infrared camera system | |
| CN115944270B (zh) | 超分辨率视网膜成像方法和系统 | |
| Talib et al. | A General Overview of Retinal Lesion Detecting Devices | |
| WO2023042577A1 (ja) | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム | |
| Wei et al. | 75-degree high-resolution wide-field OCT angiography enabled by self-tracking method | |
| WO2023238729A1 (ja) | 眼科装置、眼科装置を制御する方法、プログラム、及び記録媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210510 Address after: No. 422-21, 4th floor, block C, F2 complex building, No. 21, B District, Sanshui center science and Technology Industrial Zone, Foshan City, Guangdong Province, 528100 Patentee after: Foshan Linguistic Technology Co.,Ltd. Address before: 528000 No. 18, Jiangwan Road, Chancheng District, Guangdong, Foshan Patentee before: FOSHAN University |
|
| TR01 | Transfer of patent right |