发明名称 一种快速激光散斑血流成像系统和方法
摘要 一种快速激光散斑血流成像系统和方法,所述系统包括激光器、激光扩束器、平面镜、毛玻璃、步进电机、被测对象、CCD相机和计算机。其中,激光器发射出的激光通过激光扩束器、平面镜、毛玻璃后,照射在被测对象表面,被测对象所产生的后向散射光在CCD相机表面形成散斑,CCD相机将散斑的成像图像送至计算机进行处理,计算机遍历第p帧图像中所有像素点,计算得到所有像素点的散斑图像模糊度,根据模糊度计算对应的相对血流速度、建立相对血流速度分布图。相对于现有的激光散斑血流成像方法和系统,本发明的方法和系统具有运算量小,实时性好的优点。
申请公布号 CN103300841B 申请公布日期 2014.11.05
申请号 CN201310234441.3 申请日期 2013.06.13
申请人 上海理工大学 发明人 杨晖;杨海马;孔平;郑刚;于小强;宋磊磊
分类号 A61B5/0265(2006.01)I 主分类号 A61B5/0265(2006.01)I
代理机构 上海脱颖律师事务所 31259 代理人 脱颖
主权项 一种快速激光散斑血流成像系统,包括:激光器(1)、激光扩束器(2)、平面镜(3)、毛玻璃(4)、步进电机(5)、被测对象(6)、图像采集设备(7)、计算机(8),其特征在于,激光器发射出的激光通过激光扩束器进行扩束,经过平面镜反射后,再通过由步进电机带动的匀速旋转的毛玻璃散射后,照射在被测对象表面,被测对象对入射来的激光产生后向散射,所产生的后向散射光在图像采集设备表面形成散斑,图像采集设备对所形成的散斑进行成像,并将所成像的图像送至计算机进行处理,计算机遍历第p帧图像中所有像素点,计算得到所有像素点的散斑图像模糊度,然后计算所有像素点对应的相对血流速度,并根据像素点的相对血流速度建立相对血流速度分布图;所述计算机按公式(1)‑(3)计算第p帧图像上横纵坐标分别为x,y的像素点的散斑图像模糊度K<sub>B</sub>(x,y,p):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>K</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mi>REF</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>I</mi><mi>SD</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000537456270000011.GIF" wi="1139" he="181" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>I</mi><mi>REF</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>q</mi><msup><mrow><mo>(</mo><mn>2</mn><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mi>p</mi><mo>-</mo><mi>q</mi></mrow><mi>p</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mi>y</mi><mo>-</mo><mi>i</mi></mrow><mrow><mi>y</mi><mo>+</mo><mi>i</mi></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mi>x</mi><mo>-</mo><mi>i</mi></mrow><mrow><mi>x</mi><mo>+</mo><mi>i</mi></mrow></munderover><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000537456270000012.GIF" wi="1345" he="195" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>I</mi><mi>SD</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mi>p</mi><mo>-</mo><mi>q</mi></mrow><mi>p</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mi>y</mi><mo>-</mo><mi>i</mi></mrow><mrow><mi>y</mi><mo>+</mo><mi>i</mi></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mi>x</mi><mo>-</mo><mi>i</mi></mrow><mrow><mi>x</mi><mo>+</mo><mi>i</mi></mrow></munderover><mo>|</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>I</mi><mi>REF</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000537456270000013.GIF" wi="1356" he="190" /></maths>其中,I(x,y,p)表示图像采集设备连续采集到的第p帧图像上x,y点的像素灰度值,I<sub>REF</sub>(x,y,p)表示第p帧图像上x,y点的平均像素灰度值,I<sub>SD</sub>(x,y,p)表示第p帧图像上x,y点的统计像素灰度值与平均像素灰度值的差,i*i表示空间统计窗口的大小,q表示当前帧p之前的q帧图像;计算机遍历图像中所有的像素点,得到所有像素点的散斑图像模糊度;所述计算机按照公式V(x,y,p)=C·K<sub>B</sub>(x,y,p)计算所有像素点对应的相对血流速度V(x,y,p),其中,C为归一化系数。
地址 200093 上海市杨浦区军工路516号