发明名称 一种光声显微成像的自适应图像重建方法
摘要 本发明涉及一种针对生物组织脉管系统进行的自适应合成孔径图象重建的光声显微成像方法,首先利用光声显微成像系统获取初始的脉管系统光声成像数据;接着对原始三维成像数据进行图象预处理运算,获取成像分支脉管的特征信息;针对成像分支脉管的特征信息,对原始三维光声显微成像进行自适应合成孔径运算;最后获取经过自适应合成孔径运算后的脉管图像。本发明使离焦情况下各个不同走向的分支脉管都获得最佳的侧向空间分辨率,从而使离焦情况下系统成像的侧向分辨率达到焦点处的侧向分辨率,提高分支脉管系统成像的准确性。为脉管系统量化信息如脉管长度、脉管直径、脉管密度等形态结构参数的获取提供精确保证。
申请公布号 CN102579073B 申请公布日期 2013.09.11
申请号 CN201110392620.0 申请日期 2011.12.01
申请人 华中科技大学 发明人 骆清铭;龚辉;邓子林;杨孝全
分类号 A61B8/00(2006.01)I;A61B5/00(2006.01)I 主分类号 A61B8/00(2006.01)I
代理机构 武汉开元知识产权代理有限公司 42104 代理人 唐正玉
主权项 1.一种光声显微成像的自适应图像重建方法,其特征在于包括以下步骤:(1)通过光声显微成像系统扫描头的二维扫描移动获取原始生物组织脉管系统的三维光声显微成像数据;步骤(1)中准备成像的生物组织脉管系统吸收的ns脉冲光依次通过光纤、聚焦透镜、圆锥透镜和锥台镜后聚焦照射在吸收体上,生物组织脉管系统吸收光产生超声信号,且产生的光声信号是柱状波,光声显微成像系统中用于拾取超声信号的超声探头是高频球面聚焦超声探头,且超声探头的声焦点远小于共聚焦光的光焦点,成像过程中聚焦的光路部分和超声探头进行二维扫描获取生物组织脉管系统的光声信号;(2)对原始三维成像数据进行图象预处理运算,获取成像的生物组织脉管系统的分支脉管的骨架走向的特征信息;步骤(2)对原始三维成像数据进行图象预处理运算,获取成像的生物组织脉管系统的分支脉管的骨架走向的特征信息,具体步骤为:当超声探头沿x、y方向进行二维扫描时,对两个扫描方向的孔径进行合成;将二维合成孔径运算后的光声显微成像数据沿着A线方向进行最大值投影,获得扫描平面内的血管投影成像数据;接着调整最大值投影图为8位灰度图,并适当调整图象的灰度范围,然后利用矩阵实验室(Matlab)语言编写的滤波程序对图象进行滤波,接着重新将滤波后的图象转化为8位灰度图,最后利用Matlab图象处理工具箱中的骨架特征中的提取函数来提取图象中分支脉管的特征信息;针对不同分支的脉管进行一阶线性拟合,获得不同分支脉管在二维扫描平面内与扫描方向的夹角即走向角θ;(3)针对成像的生物组织脉管系统的分支脉管的特征信息,对原始三维光声显微成像数据进行自适应合成孔径运算;(4)获取经过自适应合成孔径运算后的生物组织脉管系统的血管图像。2、根据权利要求1所述的一种光声显微成像的自适应图像重建方法,其特征在于:步骤(3)针对成像的生物组织脉管系统的分支脉管的特征信息,对原始三维光声显微成像数据进行自适应合成孔径运算,具体步骤为:针对不同分支脉管的走向角θ,在垂直脉管走向的方向将平面内一点(x<sub>i</sub>,y<sub>j</sub>)到立体角范围内的虚拟探测器的超声信号进行叠加,可以用下式描述:<img file="987286DEST_PATH_IMAGE002.GIF" wi="252" he="45" /><img file="655159DEST_PATH_IMAGE004.GIF" wi="113" he="27" /><img file="60994DEST_PATH_IMAGE006.GIF" wi="120" he="27" />其中表示实际接收的空间中某点的初始光声信号强度,<img file="DEST_PATH_IMAGE007.GIF" wi="32" he="30" />表示光声信号从合成点到不同虚拟探测器位置时所经历的时间,<img file="708401DEST_PATH_IMAGE008.GIF" wi="101" he="26" />表示空间中某点经过自适应合成孔径算法后的光声信号大小,N为沿合成方向的最大合成线数,上面式中的[]表示不大于某个数的最大整数,k是沿合成孔径方向的单位象素长度;<img file="273506DEST_PATH_IMAGE010.GIF" wi="288" he="96" />为了进一步的提高重建图象的分辨率和信噪比,提出适合于自适应合成孔径算法的相干权重因子:其中各个参数的意义如上,这样对每个点进行自适应合成孔径算法后的信号强度可以表示成<img file="DEST_PATH_IMAGE011.GIF" wi="324" he="26" />。
地址 430074 湖北省武汉市洪山区珞瑜路华中科技大学光电国家实验室F206室