一种用于光学分辨率光声显微镜下活体肿瘤血管新生图像的血管提取方法
【专利摘要】本发明提供一种用于光学分辨率光声显微镜(OR?PAM)下活体肿瘤血管新生图像的血管提取方法，属于图像处理领域。本发明中首先使用了多种增强算法结合对肿瘤血管新生图像进行预处理，有效的增强了细小血管的对比度，最大程度的降低了噪声的影响。然后采用Otsu算法来自动选取阈值进行二值化图像，有效的降低错分概率。最后使用两步快速行进法去测定血管的中心线，从而使所计算的中心线能够更准确地定位在血管内，而且针对传统的水平集血管提取算法速度太慢的问题，使用快速行进算法进行血管提取的算法运行速度有了很大的提升。实验结果表明，该方法可以更有效的提取出肿瘤血管新生图像的血管信息，包括细小血管的信息也提取效果良好。
【专利说明】
-种用于光学分辨率光声显微镜下活体肿瘤血管新生图像的 血管提取方法
技术领域
[0001] 本发明设及一种用于光学分辨率光声显微镜(OR-PAM)下活体肿瘤血管新生图像 的血管提取方法，属于图像处理技术领域。
【背景技术】
[0002] 血管新生是指从已有的毛细血管或毛细血管后静脉发展而形成新的血管。肿瘤血 管新生是一个极其复杂的过程，一般包括血管内皮基质降解、内皮细胞移行、内皮细胞增 殖、内皮细胞管道化分支形成血管环和形成新的基底膜等步骤。肿瘤的生长和转移是一个 依赖于血管的过程，当肿瘤体积增长到一定程度时，维持其生长就需要依靠新生血管。大量 的研究表明，良性肿瘤血管新生稀少，血管生长缓慢，而大多数恶性肿瘤的血管生长密集且 生长迅速。因此，血管新生在肿瘤的发展转移过程中起到重要作用，抑制运一过程将能明显 阻止肿瘤组织的发展和扩散转移。因此研究血管新生的结构变化在肿瘤等疾病的预防与诊 治中具有非常重要的意义，而从肿瘤血管新生的图像中精确的提取出血管又是血管量化和 可视化诊断的重要条件，因此如何精确的提取血管的研究成为了肿瘤血管新生图像方向的 研究热点。
[0003] 目前超声造影是检测新生血管的重要手段之一，但是常规二维超声不能完全显示 新生血管的立体化与空间化的分布情况，作为新一代的无损医学成像技术，光声成像可W 无标记地对单个细胞成像、可W对血管形态的高分辨成像、对不同组织的成份进行解析和 对血液参数高特异性的功能检测。光声成像技术具有对比性强、灵敏度高、成像深度深的优 点，能够提高成像系统的灵敏度，从而增加探测早期癌症的可能性。
[0004] 目前光声成像的主要研究分支有光声断层成像(PAT)、光声显微成像(PAM)、光声 内窥成像(IVPAI)。光声断层成像清晰地探测到活体小鼠脑血管分布，根据血容量、血流、血 氧等参数反映了脑功能信息。光声成像技术为肿瘤的早期诊断与治疗监控提供了强大的技 术支持。随着光声显微镜的出现，光声成像发展到了一个新的阶段。光声显微镜将横向分辨 率提高了一个数量级达到了 45皿。利用光声显微成像技术不仅可W获得高分辨率黑色素瘤 的实体和周围的微血管的形态结构图像，还可W得到活体动物的血氧饱和度信息。光学分 辨率的光声显微镜(OR-PAM)可W轻而易举地对黑色素瘤细胞和血红细胞进行单细胞成像。
[0005] 目前已经有许多相关的血管提取算法，但是对于血管新生图像来说，细小的血管 比较多，又存在受到噪声污染、血管图像对比度低等原因，提取效果并不理想，不能较好地 提取出细小血管，所W如何在灰度不均匀的血管新生图像中准确的提取出血管信息，是急 需解决的问题。
【发明内容】
[0006] 为解决传统血管提取算法对于肿瘤血管新生图像不能较好地提取细小血管的问 题，实现准确提取血管新生图像中的血管信息，提高血管提取的鲁棒性，本发明提出一种融 合最大类间方差法(OtSU)与快速行进法思想，提取光声显微镜下肿瘤血管新生图像血管信 息的方法。本方法首先使用高频强调滤波等多种图像增强算法对肿瘤血管新生图像进行增 强预处理，有效的提高细小血管的对比度，再结合化SU与两步快速行进法对血管信息进行 提取。高频强调滤波算法在高频滤波器上加上一个偏移量，再将滤波器乘W-个大于1的常 数，通过常量乘数突出高频部分，也增加了低频部分的幅度，只要能保证偏移量与乘数相比 较小，就能保证低频增强的影响弱于高频增强的影响，可W有效的增强图像的对比度，并且 有效的解决了频域空间中高频滤波器增强后的图像，边缘和亮度急速变化的区域得到了增 强，但图像却失去了大部分原图像所呈现出的背景色调的问题。Otsu方法的基本思想是:设 使用某一个阔值将灰度图像根据灰度大小，分成目标部分和背景部分两类，在运两类的类 内方差最小和类间方差最大的时候，得到的阔值是最优的二值化阔值，Otsu方法计算简单， 不受图像亮度和对比度的影响，可W有效的降低错分概率。快速行进法基于当前零水平集 对像素进行分类，并按照一定机制对水平集进行演化。本发明中采用二步快速行进算法提 取血管中屯、线，可W更有效的提取出肿瘤血管新生图像中的细小血管。
[0007]为实现上述发明目的，本发明采用下述技术方案：
[000引一种用于光学分辨率光声显微镜OR-PAM下活体肿瘤血管新生图像的血管提取方 法，该方法首先使用图像增强算法对光学分辨率光声显微镜活体肿瘤血管新生图像进行预 处理，然后利用最大类间方差法Otsu提取血管新生图像的二值化图像，再利用两步快速行 进法提取血管新生图像的血管中屯、线，该方法包括W下步骤：
[0009]步骤1:使用高频强调滤波去增强血管新生图像中的细小血管的对比度；
[001日]步骤2:计算图像的Hessian矩阵；
[0011 ] 步骤3:计算化SSian矩阵的特征值；
[0012] 步骤4:对特征图进行梯度增强得到一个增强的特征映射；
[0013] 步骤5:使用化SU算法选取最佳阔值T，把血管图像提取为二值化图像；
[0014] 步骤6:基于所提取的二值化血管图像，使用两步快速行进法提取出血管的中屯、 线。
[0015] 优选的，在步骤1中采用高频强调滤波和直方图均衡法相结合得到增强后的图像。
[0016] 优选的，在步骤2使用一个尺度S，计算图像在尺度S的化SSian矩阵，通过尺度S来 控制内核的大小。
[0017] 优选的，在步骤3计算化SSian矩阵的特征值，其表示沿着对应于各特征向量的主 方向的曲率。
[0018] 优选的，在步骤4中利用梯度域图像增强方法去进一步增强血管新生图像中的细 小血管的对比度。
[0019] 本发明中，使用多种增强算法结合可W有效的提高细小血管的对比度，有效的地 降低了噪声的影响，Otsu算法使用的是聚类的思想，把图像的灰度数按灰度级分成2个部 分，使得两个部分之间的灰度值差异最大，每个部分之间的灰度差异最小，通过方差的计算 来寻找一个合适的灰度级别来划分，不仅计算简单，不受图像亮度和对比度的影响，而且可 W将错分概率降低。使用两步快速行进法去测定中屯、线，可W使所计算的中屯、线能够更准 确地定位在血管内，针对传统的水平集血管提取算法速度太慢的问题，快速行进算法进行 血管提取算法运行速度提升了许多。实验结果表明，使用该方法，可W成功地从原始图像中 提取出大多数的血管，包括细小血管的提取效果也良好，并且所提取血管的连续性和强度 的均匀性都优于其他传统的提取方法。
【附图说明】
[0020] 图1是本发明所述的用于OR-PAM活体肿瘤血管新生图像的血管提取方法的流程 图；
[0021] 图2是选取的老鼠活体的OR-PM成像的某个原始图像片；
[0022] 图3是对图像片二值化后的图像；
[0023] 图4是对图像片提取的血管中屯、线效果图；
[0024] 图5是图像片提取的最终血管骨架图。
【具体实施方式】
[0025] 本发明提出一种用于光学分辨率光声显微镜OR-PAM下活体肿瘤血管新生图像的 血管提取方法。首先使用多种增强算法结合可W有效的提高细小血管的对比度，有效地降 低了噪声的影响，然后使用Otsu算法对增强后的图像进行二值化图像，采用Otsu算法来自 动选取阔值进行二值化是图像分割中阔值选取的最佳算法，不仅计算简单，不受图像亮度 和对比度的影响，而且可W将错分概率降至最小。最后使用两步快速行进法测定中屯、线，使 所计算的中屯、线能够更准确地定位在血管内，针对传统的水平集血管提取算法速度太慢的 问题，快速行进算法进行血管提取算法运行速度有了很大的提升。图1中明确表示了本发明 所述方法的整体流程，具体实施步骤如下：
[0026] 步骤1:使用高频强调滤波去增强血管新生图像中的细小血管的对比度，高频强调 滤波定义如下：
[0027] Hhfe(u, V) =a+bHhp(u, V)
[002引其中Hhp (u，v)是一个高通算子，U和V代表二维空间频率分量，a指偏移量，b指乘数；
[0029] 步骤2:计算图像在尺度S的化SSian矩阵，定义如下：
[0030]
[0031] 其中I(x，y，s)为Kx, y)卷积，S用来控制高斯平滑的内核的大小，X和y代表空间坐 标；
[0032] 步骤3:计算化SSian矩阵的特征值，用来表示沿着对应于各特征向量的主方向的 曲率，为了进一步提高微血管和抑制噪声，定义特征映射f在尺度S为：
[0033]
[0034] 其中|Ai| < |A2|，e为一个很小的数，Al和A2分别为Hessian矩阵的特征值；
[0035] 步骤4:在经过上面的步骤之后，已经可W在特征图中显示出大多数血管，但毛细 血管地区仍然有些暗淡，所W在运一步对图像的梯度场进行了直方图均衡化，使图像的低 频区域的细节能够在梯度域得到增强，再利用最小二乘算法重建梯度域图像，即可得到梯 度域增强后的图像；
[0036] 步骤5:使用化SU算法选取最佳阔值T，把血管图像提取为二值化图像。Otsu算法步 骤如下：
[0037] (1)计算输入图像的归一化直方图，使用914 = 〇，1，2，...丄-1表示该直方图的各 个分量，L代表图像中可能出现的灰度级总数；
[003引（2)使用公； ，对于4 = 0，1，2，...，心1，计算累积和口1化）；
[0039] (3)使用公； ，对于4 = 0，1，2，...，心1，计算累积均值111化）；
[0040] (4)使用公； 算全局灰度均值niG;
[004U (5)使用公;
，对于k = 0，l，2，. ..，1^-1，计算类间方差
[0042] (6)得到Os化阔值即使得好別巧最大的k值，如果最大值不唯一，则用相应检测到 的各个最大值的平均值来作为k%
[0043] (7)在k = k郑计算
得到可分性度量rf。
[0044] 步骤6:基于所提取的二值化血管图像，使用两步快速行进法提取出血管的中屯、 线。首先将血管网络分离为各个子图像，使每一个子图像中仅包括一组相互连接的血管树。 第一次使用快速行进算法在每一个子图像内，血管内的点被用作种子在初始速度均匀地扩 大，直到它们的轮廓表面达到接近血管的边界。第二次使用快速行进算法，在血管树的一个 分支提取中屯、线，使用回溯算法进行连接，使用迭代的方法，直到所有的血管树的分支的中 屯、线都被提取，最后提取出一个血管树的所有分支的中屯、线，形成运个血管树的骨架。
[0045] 由于在本方法的实验中采集的小白鼠活体数据为=维体数据，应用本方法进行处 理时运行时间过长，而且结果图像太大，看不出细节效果，所W运里只列出选取的某个图像 片的处理效果，图2为选取的老鼠活体的光学分辨率光声显微镜OR-PAM下成像的某个原始 图像片；
[0046] 针对图2所选的图像片应用化SU算法对其进行二值化，图3是对图像片二值化后的 图像效果；
[0047] 针对图3所示的二值化后的图像应用快速行进算法进行血管提取，在提取前先将 血管网络分离成只包含一组相互连接的子血管树，先对每个子血管树进行提取，最后将各 个子血管树的提取连接起来，图4是对图像片提取的血管中屯、线效果图，图5是图像片提取 的最终血管骨架图。
[0048] W上公开的仅为本发明的具体实施例，根据本发明提供的思想，本领域的技术人 员能思及的变化，都应落入本发明的保护范围内。
【主权项】
1. 一种用于光学分辨率光声显微镜OR-ΡΑΜ下活体肿瘤血管新生图像的血管提取方法， 该方法首先使用图像增强算法对光学分辨率光声显微镜下活体肿瘤血管新生图像进行预 处理，然后利用最大类间方差法Otsu提取血管新生图像的二值化图像，再利用两步快速行 进法提取血管新生图像的血管中心线，该方法包括以下步骤： 步骤1:使用高频强调滤波去增强血管新生图像中的细小血管的对比度； 步骤2:计算图像的Hessian矩阵； 步骤3:计算He s s i an矩阵的特征值； 步骤4:对特征图进行梯度增强得到一个增强的特征映射； 步骤5:使用Otsu算法选取最佳阈值T，把血管图像提取为二值化图像； 步骤6:基于所提取的二值化血管图像，使用两步快速行进法提取出血管的中心线。2. 根据权利要求1所述的用于光学分辨率光声显微镜OR-ΡΑΜ下活体肿瘤血管新生图像 的血管提取方法，其特征在于:在步骤1中采用高频强调滤波和直方图均衡相结合得到增强 后的图像。3. 根据权利要求1所述的用于光学分辨率光声显微镜OR-ΡΑΜ下活体肿瘤血管新生图像 的血管提取方法，其特征在于:在步骤2使用一个尺度s，计算图像在尺度s的Hessian矩阵， 通过尺度s来控制内核的大小。4. 根据权利要求1所述的用于光学分辨率光声显微镜OR-ΡΑΜ下活体肿瘤血管新生图像 的血管提取方法，其特征在于:在步骤3计算Hessian矩阵的特征值，其表示沿着对应于各特 征向量的主方向的曲率。5. 根据权利要求1所述的用于光学分辨率光声显微镜OR-ΡΑΜ下活体肿瘤血管新生图像 的血管提取方法，其特征在于:在步骤4中利用梯度域图像增强方法去进一步增强血管新生 图像中的细小血管的对比度。
【文档编号】G06T7/00GK106023165SQ201610312546
【公开日】2016年10月12日
【申请日】2016年5月12日
【发明人】朱荣, 孟静, 龚小竞, 刘成波
【申请人】曲阜师范大学