基于连通域优化的超记忆效应范围非侵入式散射成像方法与流程

本发明涉及计算机成像学领域，尤其涉及一种基于连通域优化的超记忆效应范围非侵入式散射成像方法。

背景技术：

散射介质在自然界中广泛存在，散射介质具有不均匀的折射率和反射率分布，光通过散射介质后原有信息被严重干扰，阻碍了传统光学系统对象信息的直接分析。散射介质作用下的成像技术在生物医学成像、遥感等领域都有很好的应用前景。

现有的散射成像方法主要分为两类：基于弹道光或基于记忆效应理论。其中，基于弹道光的方法包括光学相干断层扫描(oct)、光子累积成像、时域光子统计分布和双光子显微镜等。这些方法需要大型设备作为支持，同时由于弹道光的数量随散射介质散射作用的增强呈指数式衰减，因此这一类方法只适用于散射效应较弱的场景。基于记忆效应的技术包含波前整形、散斑扫描采集系统、单次散斑相关成像系统和psf反卷积方法，其中波前整形和反卷积方法通常需要系统校准或引导星等，由于其侵入式的成像方法，应用场景受到了较大的限制。散斑扫描采集系统和单次散斑相关成像系统虽然可以实现强散射介质的非侵入式成像，但成像范围都受到记忆效应范围的限制，无法对大视角物体成像。

以上背景技术内容的公开仅用于辅助理解本发明的构思及技术方案，其并不必然属于本专利申请的现有技术，在没有明确的证据表明上述内容在本专利申请的申请日已经公开的情况下，上述背景技术不应当用于评价本申请的新颖性和创造性。

技术实现要素：

为解决上述技术问题，本发明提出一种基于连通域优化的超记忆效应范围非侵入式散射成像方法，能够实现大视角范围的多个物体散射成像。

为了达到上述目的，本发明采用以下技术方案：

本发明公开了一种基于连通域优化的超记忆效应范围非侵入式散射成像方法，包括以下步骤：

a1：采集待成像的样本的散斑信息，其中待成像的样本包括尺寸在记忆效应范围内但距离超出记忆效应范围的连续的第一物体和第二物体；

a2：计算步骤a1采集的散斑信息的自相关图像，并根据自相关图像得到自相关图形；

a3：根据步骤a2的自相关图形得到第一物体或第二物体的重建结果，并利用连通域优化第一物体或第二物体的重建结果；

a4：计算经步骤a3优化后的第一物体或第二物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第二物体或第一物体的初始自相关；

a5：根据前一步骤得到的第二物体或第一物体的自相关相应得到第二物体或第一物体的重建结果，并利用连通域优化第二物体或第一物体的重建结果，然后计算优化后的第二物体或第一物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第一物体或第二物体的自相关；

a6：根据前一步骤得到的第一物体或第二物体的自相关相应得到第一物体或第二物体的重建结果，并利用连通域优化第一物体或第二物体的重建结果，然后计算优化后的第一物体或第二物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第二物体或第一物体的自相关；

a7：循环进行步骤a5和a6直到预定的循环次数，再利用相位恢复算法和计算得到的第一物体的自相关和第二物体的自相关进行空域重构，实现非侵入式成像。

优选地，步骤a1中具体包括：使用基于单次散斑相关成像系统并设置被散射介质遮挡的尺寸在记忆效应范围之内但距离超出记忆效应范围的第一物体和第二物体作为待成像的样本，利用非相关光照射采集待成像的样本的散斑信息i为：

i＝o1*psf1+o2*psf2

其中，o1、o2分别代表第一物体和第二物体，psf1和psf2分别代表第一物体和第二物体对应散射介质部分的点扩散函数，*代表卷积运算。

优选地，步骤a2具体包括：计算步骤a1采集的散斑信息的自相关图像为：

i★i＝(o1*psf1+o2*psf2)★(o1*psf1+o2*psf2)

＝(o1★o1)+(o2★o2)+c

其中，★指自相关运算，c为背景项；

根据自相关图像得到自相关图形d为：d＝i★i-c。

优选地，步骤a3中根据步骤a2的自相关图形得到第一物体或第二物体的重建结果具体包括：对步骤a2的自相关图形使用多次相位恢复算法，得到恢复误差最小的第一物体或第二物体的重建结果。

优选地，步骤a3、a5和a6中利用连通域优化第一物体的重建结果或者利用连通域优化第二物体的重建结果分别具体包括：

将物体的重建结果sk移动到中心，然后将物体的重建结果sk按连通域分为n个区域，n代表重建结果连通域数目，且为不小于2的正整数，并按照连通域的大小降序排列为：ε1，ε2……εn；归一化物体的重建结果sk；得到利用连通域优化后的物体的重建结果s′k为：

其中，k为1时，sk代表第一物体的重建结果，s′k代表利用连通域优化后的第一物体的重建结果；k为2时，sk代表第二物体的重建结果，s′k代表利用连通域优化后的第二物体的重建结果；averagesk代表重建结果中所有大于0.1的像素点的均值，ε＝{(x,y),|(x,y)∈ε1，ε2}。

优选地，步骤a4中具体包括：计算经步骤a3优化后的第一物体或第二物体的重建结果的自相关并归一化，并根据自相关图形减去优化后的第一物体或第二物体的重建结果的自相关进行归一化后的一半，相应计算得到第二物体或第一物体的初始自相关。

优选地，步骤a5和a6中根据前一步骤得到的第一物体的自相关得到第一物体的重建结果具体包括：对第一物体的自相关使用多次相位恢复算法，得到恢复误差最小的结果即为第一物体的重建结果；

步骤a5和a6中根据前一步骤得到的第二物体的自相关得到第二物体的重建结果具体包括：对第二物体的自相关使用多次相位恢复算法，得到恢复误差最小的结果即为第二物体的重建结果。

优选地，步骤a5和a6中计算优化后的第一物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第二物体的自相关具体包括：计算优化后的第一物体的重建结果的自相关并归一化，并根据自相关图形减去优化后的第一物体的重建结果的自相关进行归一化后的一半，相应计算得到第二物体的初始自相关；

步骤a5和a6中计算优化后的第二物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第一物体的自相关具体包括：计算优化后的第二物体的重建结果的自相关并归一化，并根据自相关图形减去优化后的第二物体的重建结果的自相关进行归一化后的一半，相应计算得到第一物体的初始自相关。

优选地，步骤a7具体包括：循环进行步骤a5和a6直到预定的循环次数后第一物体和第二物体的自相关被正确分开，然后利用相位恢复算法和计算得到的第一物体的自相关和第二物体的自相关进行空域重构，实现非侵入式成像。

与现有技术相比，本发明的有益效果在于：本发明公开的基于连通域优化的超记忆效应范围非侵入式散射成像方法，采集的散斑信息的自相关近似等于每个单个对象的散斑自相关的线性叠加，并基于线性减法的自相关分离部分来分离混合信号，同时基于连通域优化的方法保证迭代优化的方向收敛，最终实现准确地分离和恢复自相关完全混叠的两个对象，实现两个距离超出记忆效应范围的物体的非侵入式成像，从而实现大视角范围的多个物体散射成像。本发明的方法，只需采集一次样本，就可以同时实现两个距离超出记忆效应范围的物体非侵入式成像，并且除了要求每个物体只有单个连通域，两个物体透光面积和透光率相近外不需要任何其他先验信息，极大地提升了散斑相关方法的成像范围。

附图说明

图1是本发明优选实施例的基于连通域优化的超记忆效应范围非侵入式散射成像方法的流程图；

图2是本发明另一优选实施例的基于连通域优化的超记忆效应范围非侵入式散射成像方法的流程图。

具体实施方式

下面对照附图并结合优选的实施方式对本发明作进一步说明。

如图1所示，本发明优选实施例公开了一种基于连通域优化的超记忆效应范围非侵入式散射成像方法，包括以下步骤：

a1：采集待成像的样本的散斑信息，其中待成像的样本包括两个尺寸在记忆效应范围内的连续物体；

具体地，使用基于单次散斑相关成像系统并设置两个被散射介质遮挡的大小在记忆效应范围之内但距离超出记忆效应范围的物体作为待成像的样本，利用非相关光照射采集样本的散斑信息散斑可表示为：

i＝o1*psf1+o2*psf2

其中，o1、o2分别代表第一物体和第二物体，psf1和psf2分别代表第一物体和第二物体对应散射介质部分的点扩散函数，*代表卷积运算。

a2：计算步骤a1采集的散斑信息的自相关图像，并根据自相关图像得到自相关图形；

具体地，将采集的散斑信息除以其低频滤波图像包络从而锐化散斑图像并减少噪音，然后计算拟合后的散斑信息的自相关图像，再减去背景项(即散斑自相关图像强度的最小值)得到自相关图形；

其中拟合后的散斑信息的自相关图像的计算公式为：

i★i＝(o1*psf1+o2*psf2)★(o1*psf1+o2*psf2)

＝(o1★o1)+(o2★o2)+c

其中，i是指散斑信息，★指自相关运算，c为背景项。

从上式可以看出散斑信息的自相关图像为两物体自相关之和加上背景项c，为了去除背景项c，减去i★i的最小值(即为背景项c)并归一化即得到自相关图形d＝i★i-c。

a3：对自相关图形使用多次相位恢复算法，得到恢复误差最小的重建结果，并利用连通域优化该重建结果；

其中，连通域优化包括：

(1)将物体的重建结果sk移动到中心；

(2)将物体的重建结果sk按连通域分为多个区域，并按照连通域的大小降序排列为：ε1，ε2……εn(物体的重建结果由多个连通域组成，本实施例中按连通域理论拆分为n个连通域，n代表重建结果连通域数目，且为不小于2的正整数，并按照连通域的面积大小排序，最大的ε1，其次为ε2，以此类推)；

(3)令ε＝{(x,y),|(x,y)∈ε1，ε2}；

(4)归一化sk；

(5)令

其中，sk代表物体的重建结果，s′k代表利用连通域优化后的物体的重建结果；具体地：k为1时，sk代表第一物体的重建结果，s′k代表利用连通域优化后的第一物体的重建结果；k为2时，sk代表第二物体的重建结果，s′k代表利用连通域优化后的第二物体的重建结果；averagesk代表重建结果中所有大于0.1的像素点的均值。

a4：计算优化后的第一物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第二物体的初始自相关；

具体地，计算s′1的自相关，并做归一化处理得到c1，用自相关图形d减去强度减半后的c1得到c2即：

(1)c1＝s′1★s′1；

(2)归一化c1；

(3)c2＝s2★s2＝d-0.5·c1。

a5：根据前一步骤得到的第二物体的自相关相应得到第二物体的重建结果，并利用连通域优化第二物体的重建结果，然后计算优化后的第二物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到新的第一物体的自相关；

其中步骤a5中的根据前一步骤具体是指：在第一次执行步骤a5时，该前一步骤是指步骤a4，也即前一步骤得到的第二物体的自相关是指步骤a4得到的第二物体的初始自相关；在执行步骤a7的循环步骤时，也即第二次或者第二次之后执行步骤a5时，该前一步骤是指前一次循环的步骤a6，也即前一步骤得到的第二物体的自相关是指前一次循环的步骤a6得到的第二物体的自相关。

具体地，对第二物体的自相关c2使用多次相位恢复算法得到恢复误差最小的重建结果s2；利用连通域对重建结果s2进行优化得到优化后的重建结果s′2，计算优化后的重建结果s′2的自相关，并做归一化处理得到c2，用自相关图形d减去强度减半的c2得到c1即：

(1)第一次执行步骤a5时：将第二物体的初始自相关c2作为输入，使用多次相位恢复算法，取其中误差最小的结果作为本次的第二物体的重建结果s2；第二次或者第二次之后执行步骤a5时：将上一循环的步骤a6中的第二物体的自相关c2作为输入，使用多次相位恢复算法，取其中误差最小的结果作为本次的第二物体的重建结果s2；

(2)使用步骤a3中的连通域优化方法优化重建结果s2得到s′2；

(3)c2＝s′2★s′2；

(4)归一化c2；

(5)c1＝s1★s1＝d-0.5·c2。

a6：根据前一步骤得到的第一物体的自相关相应得到第一物体的重建结果，并利用连通域优化第一物体的重建结果，然后计算优化后的第一物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第二物体的自相关；

其中步骤a6中的根据前一步骤具体是指：相应循环次数的步骤a5。

具体地，对第一物体的自相关c1使用多次自相关算法得到恢复误差最小的重建结果s1，利用连通域对重建结果s1进行优化得到优化后的重建结果s′1，计算优化后的重建结果s′1的自相关，并作归一化处理得到c1，用自相关图形d减去强度减半的c1得到c2即：

(1)新的第一物体的自相关c1作为输入，使用多次相位恢复算法，取其中误差最小的结果作为本次的第一物体的重建结果s1；

(2)使用步骤a3中的连通域优化方法优化重建结果s1得到s′1；

(3)c1＝s′1★s′1；

(4)归一化c1；

(5)c2＝s2★s2＝d-0.5·c1。

a7：循环进行步骤a5和a6，直到达到设置的循环次数，分别对输出的两个物体的自相关使用相位恢复算法实现两个物体的成像。

具体地，多次循环步骤a5和步骤a6后，第一物体的重建结果s1和第二物体的重建结果s2分别越来越接近真实的第一物体与第二物体，在经过一定的循环次数后两个物体的自相关被正确分开，最后利用相位恢复算法实现两个距离超出记忆效应范围的物体的非侵入式成像。

如图2所示，在另一优选实施例中，步骤a4中也可为：计算优化后的第二物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第一物体的初始自相关；相应的步骤a5为：根据前一步骤得到的第一物体的自相关相应得到第一物体的重建结果，并利用连通域优化第一物体的重建结果，然后计算优化后的第一物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到新的第二物体的自相关；相应的步骤a6为：根据前一步骤得到的第二物体的自相关相应得到第二物体的重建结果，并利用连通域优化第二物体的重建结果，然后计算优化后的第二物体的重建结果的自相关并归一化，并结合自相关图形，相应计算得到第一物体的自相关。

本发明的优选实施例提出的基于连通域优化的超记忆效应范围非侵入式散射成像方法，当每个对象落入记忆效应范围内并且距离超出记忆效应范围时，采集的散斑的自相关可以简化为每个单个对象散斑自相关的线性叠加；并基于线性减法的自相关分离部分来分离混合信号，基于相位恢复算法恢复自相关丢失的相位信息，基于连通域优化的方法以保证迭代优化的方向收敛，最终准确地分离和恢复自相关完全混叠的两个对象，实现两个距离超出记忆效应范围的物体非侵入式成像，从而实现大视角范围的多个物体散射成像。通过本发明优选实施例的方法，只需采集一次样本，就可以同时实现两个距离超出记忆效应范围的物体非侵入式成像，并且除了要求每个物体只有单个连通域(也即每个物体空间连通)，两个物体透光面积和透光率相近(也即两个物体强度相当)外不需要任何其他先验信息，极大地提升了散斑相关方法的成像范围。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干等同替代或明显变型，而且性能或用途相同，都应当视为属于本发明的保护范围。