压缩采样光声显微成像方法及装置的制造方法
【技术领域】
[0001] 本发明涉及图像处理领域，具体涉及一种压缩采样光声显微成像方法及装置。
【背景技术】
[0002] 光声显微镜成像技术（PhotoacousticMicroscopy，PAM)由于其能够对深层组织 进行高分辨率成像受到广泛地关注，并被用于很多生物医学研究，如乳腺癌检测。
[0003]目前的光声显微镜成像技术一般是对成像区域逐点扫描实现的。为了获取更高分 辨率的图像，就需要更密集的采样，这将大大增加PAM信号的采集和处理时间。目前为了实 现该目的，大多是通过改进成像系统或成像原理来提高成像分辨率，这无疑会增加系统的 成本。
【发明内容】
[0004] 有鉴于此，本发明实施例提供一种压缩采样光声显微成像方法及装置，在不增加 系统成本的同时，能够实现高分辨率的光声显微图像的快速采集。
[0005] 为此目的，本发明提出一种压缩采样光声显微成像方法，包括：
[0006] 根据预设的待采样目标的感兴趣区域的采样比、待采样目标的背景区域的采样 比、水平方向采样点数和垂直方向采样点数，并基于边膨胀图理论生成压缩采样模板；
[0007] 利用所述压缩采样模板获取光声显微压缩采样数据矩阵；
[0008] 利用低秩矩阵填充方法对所述光声显微压缩采样数据矩阵进行恢复，得到光声显 微图像。
[0009] 另一方面，本发明提出一种压缩采样光声显微成像装置，包括：
[0010] 压缩采样模板生成单元，用于根据预设的待采样目标的感兴趣区域的采样比、待 采样目标的背景区域的采样比、水平方向采样点数和垂直方向采样点数，并基于边膨胀图 理论产生压缩米样模板；
[0011] 压缩采样数据获取单元，用于利用所述压缩采样模板获取光声显微压缩采样数据 矩阵；
[0012] 光声显微图像生成单元，用于利用低秩矩阵填充方法对所述光声显微压缩采样数 据矩阵进行恢复，得到光声显微图像。
[0013] 本发明实施例所述的压缩采样光声显微成像方法及装置，基于边膨胀图理论产生 不同采样比压缩采样模板，利用该压缩采样模板获取光声显微压缩采样数据矩阵，直接在 该压缩采样数据矩阵的基础上，采用低秩矩阵填充方法获得高分辨率的光声显微图像，不 需要改进成像系统或成像原理，因而相较于现有技术，能够在不增加系统成本的基础上，实 现光声显微图像的快速采集。
【附图说明】
[0014]图1为本发明压缩采样光声显微成像方法一实施例的流程示意图；
[0015] 图2为一个传统光声显微图像；
[0016]图3为本发明得到的光声显微图像（图3与图2的待成像目标相同）；
[0017] 图4为本发明压缩采样光声显微成像装置一实施例的方框结构示意图。
【具体实施方式】
[0018] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明 一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有 做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0019] 如图1所不，本实施例公开一种压缩米样光声显微成像方法，包括：
[0020]S1、根据预设的待采样目标的感兴趣区域的采样比、待采样目标的背景区域的采 样比、针对目标传统采样模式下的水平方向采样点数和垂直方向采样点数，并基于边膨胀 图理论产生压缩采样模板；
[0021]S2、利用所述压缩采样模板获取光声显微压缩采样数据矩阵；
[0022]S3、利用低秩矩阵填充方法对所述光声显微压缩采样数据矩阵进行恢复，得到光 声显微图像。
[0023] 本发明实施例所述的压缩采样光声显微成像方法，基于边膨胀图理论产生不同采 样比压缩采样模板，利用该压缩采样模板获取光声显微压缩采样数据矩阵，直接在该压缩 采样数据矩阵的基础上，采用低秩矩阵填充方法获得高分辨率的光声显微图像，不需要改 进成像系统或成像原理，因而相较于现有技术，能够在不增加系统成本的基础上，实现光声 显微图像的快速采集。
[0024] 可选地，在本发明压缩采样光声显微成像方法的另一实施例中，所述利用所述压 缩采样模板获取光声显微压缩采样数据矩阵，包括：
[0025] 对于每一个待采样位置，判断所述压缩采样模板在该位置处的值是否为1，若为 1，则采集该位置处待成像目标的光声数据，并记录到光声显微压缩采样数据矩阵的相应位 置。
[0026] 本发明实施例中，所述压缩采样模板在每一个待采样位置处的值为1或者0,对于 任一待采样位置，若所述压缩采样模板在该位置处的值为1，则采集光声数据，并记录到光 声显微压缩采样数据矩阵的相应位置，否则，不进行采集，最后生成的光声显微压缩采样数 据矩阵的值包含对待成像目标的光声数据进行选择采集的结果。
[0027] 可选地，在本发明压缩采样光声显微成像方法的另一实施例中，所述利用低秩矩 阵填充方法对所述光声显微压缩采样数据矩阵进行恢复，得到光声显微图像，包括：
[0028] 基于低秩矩阵填充构建高分辨率光声显微图像恢复模型 其中，τ为常数，| | | |，为矩阵的核范数，1￡为所述光声显微压缩采样数据矩阵，
Μ、,为矩阵Μ第i行第j列位置处的值，E(i，j)为所述压缩采样 模板第i行第j列位置处的值，所述恢复模型的输入为IE，输出为所述光声显微图像的灰度 值矩阵I;
[0029] 利用凸优化方法求解所述恢复模型，得到所述光声显微图像。
[0030] 本发明基于边膨胀图设计的压缩采样模板，直接利用压缩采样数据矩阵，使用低 秩矩阵填充的方法获取高分辨率的光声显微图像，能够实现PAM数据的快速采集。
[0031] 可选地，在本发明压缩采样光声显微成像方法的另一实施例中，所述利用凸优化 方法求解所述恢复模型，得到所述光声显微图像，包括如下步骤：
[0032] S321、初始化Y= 0，k= &，其中，矩阵Y的行数为所述水平方向采样点数，矩阵Y 的列数为所述垂直方向采样点数，&为预设的整数；
[0033] S322、对Y做奇异值分解，得到矩阵U、Σ和V，其中，UXVT=SVD(Y);
[0034]S323、计算对角矩阵f，其丨Σ"-r,0}，为矩阵Σ第i行第i列位 置处的值，良为矩阵f第i行第i列位置处的值；
[0035] S324、计算矩阵I，其中，/ =
[0036] 3325、更新￥的值为￥+5(?^1〇^^(1))，其中，6为常数；
[0037] S326、更新k的值为k-Ι;
[0038] S327、判断k是否为0,若k为0,则根据所述矩阵I生成所述光声显微图像，否则， 重复执行步骤S322至S327直至k为0,并根据所述矩阵I生成所述光声显微图像（对于同 一待成像目标，如图2所示为一个传统光声显微图像，图3所示为本发明得到的光声显微图 像）。
[0039] 可选地，在本发明压缩米样光声显微成像方法的另一实施例中，δ=N2| |IE| |2,N2 为预设的整数。
[0040] 本发明所述基于边膨胀图和低秩矩阵填充的压缩采样光声显微成像方法，该方法 包括以下步骤：
[0041] 步骤一：根据边膨胀图理论产生不同采样比的压缩采样模板E;步骤11 :初始化以 下参数：感兴趣区域采样比、其他背景区域采样比α2、全采样条件下水平方向采样点数 Μ和垂直方向采样点数Ν。步骤12 :根据边膨胀图理论产生压缩采样模板Ε，其中Ε为Μ行 Ν列的0，1矩阵。
[0042] 步骤二：利用压缩采样模板获取光声显微压缩采样数据矩阵ΙΕ;
[0043] 步骤21 :给定采样的初始位置，并初始化压缩采样数据矩阵^为0。
[0044] 步骤22 :若该位置处Ε的值为1，则采集光声数据，并记录到ΙΕ的对应位置；若为 0,则不进行采集。
[0045] 步骤23 :移动到下一个采样点处，重复步骤22,直到处理完所有的Μ*Ν个采样位 置。
[0046] 步骤三：利用低秩矩阵填充方法对压缩采样数据^进行恢复，得到高分辨率的光 声显微图像?。
[0047] 步骤31 :构建基于低秩矩阵填充的高分辨率光声显微图像恢复模型
[0049] 其中III|$为矩阵的核范数，即矩阵的奇异值之和，
[0050] 步骤32 :利用凸优化方法求解上述模型，得到矩阵I。
[0051 ] 步骤321 :给定τ，δ，k参数值，并初始化Y= 0 ;
[0052] 步骤322 :对Y做奇异值分解，UΣVT=SVD⑴；
[0053] 步骤323 :对分解得到的奇异值做阈值处理，?ν/ =max%.,. -Μ)；·;
[0054] 步骤324 :用更新的奇异值得到/ = ;
[0055] 步骤 325:更新Y，Y-Y+δ(PE (IE) -PE ⑴）；
[0056] 步骤 326 :更新k，k=k_l;
[0057] 步骤326 :若k= 0,进行步骤327 ;否则回到步