1.一种激光反射层析成像欠采样的重建方法,其特征在于,所述重建方法包含以下步骤:步骤一、激光器发射窄激光脉冲照射目标,并以等投影角度间隔围绕目标进行激光扫描,获取完备角度下投影角度欠采样的多视角目标原始投影数据;步骤二、对完备角度下投影角度欠采样的多视角目标原始投影数据进行反投影中心对准处理;步骤三、利用局部哈达玛矩阵,对经反投影中心对准处理后的投影数据进行稀疏观测采样;步骤四、对稀疏观测采样后的投影数据依次进行联合代数重建和全变差最小化调整的迭代处理;步骤四的具体实现过程为:步骤401、设待重建的激光反射层析图像的初始值 k为迭代次数;步骤402、将经过投影中心对准处理后的投影数据赋值给投影矩阵P;并利用快速网格遍历法求解出系统矩阵W;步骤403、根据待重建激光反射层析图像的初始值以及投影矩阵P和系统矩阵W,按如下公式,计算得到 其中, 为第k次SART重建中由P中投影数据第S个采样点的投影值所得到的待重建图像的像素值; 为第k-1次SART重建的由P中投影数据第S个采样点的投影值所得到的待重建图像的像素值;P s 为投影矩阵P中投影数据第S个采样点的投影值;W s 为系统矩阵W的第S行行向量; 为W s 的转置向量;s=1,2,…,M,M为目标原始投影数据离散采样的点数;λ为松弛因子;步骤404、判断 是否大于等于0,如是,则令 进入步骤405;如否,则令F(k)=0,进入步骤405;步骤405、初始化 取 步骤406、计算增量因子 其中, 为待重建的激光反射层析图像的像素值初值;步骤407、计算全变差梯度 及负梯度方向 其中, f i,j 为待重建图像的第i行第j列像素值,f i-1,j 为待重建图像的第i-1行第j列像素值,f i,j-1 为待重建图像第i行第j-1列像素值,i=1,2,…N 1 ,j=1,2,…N 2 ,N 1 和N 2 分别为待重建图像的行数和例数; 为第n次TV最小化时所得到的待重建图像的像素值,n=1,2,…,m,m为TV最小化过程的最大迭代次数;步骤408、按照如下公式沿TV负梯度迭代修正图像:其中, 为第n次TV最小化时所得到的待重建图像的像素值; 为第n-1次TV时的第k-1次SART迭代的由P中投影数据第s个采样点的投影值所得到的待重建图像的像素值;α为调节因子;步骤409、判断n是否小于m,如是,则令n=n+1,返回步骤408;如否,再进入步骤410;步骤410、判断 是否小于等于阈值,如是,则令 结束;如否,则进入步骤411;步骤411、计算 和 步骤412、当 则令 进入步骤403;反之,则令k=k+1,进入步骤408。
2.根据权利要求1所述的重建方法,其特征在于,步骤一的具体实现过程为:步骤11、激光器发射窄激光脉冲照射目标,并以等投影角度间隔围绕目标物体进行“一步一停”式的360°范围内激光扫描;步骤12、探测接收目标反射的激光回波脉冲信号,并对回波信号进行离散采样,从而得到完备角度下投影角度欠采样的目标原始投影数据。
3.根据权利要求2所述的重建方法,其特征在于,所述窄激光脉冲宽度为100ps;所述探测接收带宽为4GHz,离散采样速率为20GS/s。
4.根据权利要求1或2所述的重建方法,其特征在于,步骤一中,所述投影角度间隔的取值范围为6°≤Δφ≤12°。
5.一种存储介质,其特征在于,所述存储介质存储有计算机执行指令;所述计算机执行指令执行时,实现权利要求1~4中任意一项所述的重建方法。
6.一种激光反射层析成像欠采样的重建系统,其特征在于,所述重建系统包括权利要求5所述的存储介质。