CN113112563B - 一种优化区域知识先验的稀疏角cb-xlct成像方法 - Google Patents

一种优化区域知识先验的稀疏角cb-xlct成像方法 Download PDF

Info

Publication number
CN113112563B
CN113112563B CN202110429169.9A CN202110429169A CN113112563B CN 113112563 B CN113112563 B CN 113112563B CN 202110429169 A CN202110429169 A CN 202110429169A CN 113112563 B CN113112563 B CN 113112563B
Authority
CN
China
Prior art keywords
sparse
xlct
model
angle
convex
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110429169.9A
Other languages
English (en)
Other versions
CN113112563A (zh
Inventor
张海波
寇姣姣
海琳琦
黄启辰
刘英杰
耿国华
周明全
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NORTHWEST UNIVERSITY
Original Assignee
NORTHWEST UNIVERSITY
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NORTHWEST UNIVERSITY filed Critical NORTHWEST UNIVERSITY
Priority to CN202110429169.9A priority Critical patent/CN113112563B/zh
Publication of CN113112563A publication Critical patent/CN113112563A/zh
Application granted granted Critical
Publication of CN113112563B publication Critical patent/CN113112563B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开一种优化区域知识先验的稀疏角CB‑XLCT成像方法,包括:步骤1,根据稀疏角CB‑XLCT成像不适定逆问题构建稀疏非凸Lp,0<p<1正则模型;步骤2,将稀疏非凸Lp正则模型转化为重加权凸L1正则模型;步骤3,设计迭代重加权裂分增广拉格朗日收缩算法求解;步骤4,结合经典非凸p值,求解稀疏解集对应的可行区域集;步骤5,将可行区域集优化为稳定的知识先验进而指导最终逆问题的求解;步骤6,根据步骤5中的求解结果显示重建后的稀疏角CB‑XLCT图像;本发明不仅有效解决了稀疏角CB‑XLCT成像逆问题,而且具有良好的可扩展性。

Description

一种优化区域知识先验的稀疏角CB-XLCT成像方法
技术领域
本发明医学成像相关技术领域,具体涉及一种优化区域知识先验的稀疏角 CB-XLCT成像方法。
背景技术
X射线发光断层成像(XLCT)做为一种新兴的三维分子影像模态,已成为医学成 像领域的研究热点。相对于传统的生物发光断层成像(BLT)和荧光分子断层成像 (FMT),XLCT消除了FMT自体荧光的干扰,弥补了BLT成像深度有限的不足。
根据CT射线放源方式,传统XLCT技术主要分为笔束型XLCT(PB-XLCT)、扇形束XLCT(FB-XLCT)以及锥束型XLCT(CB-XLCT)。尽管PB-XLC和FB-XLCT具有成像分辨率高的 优点,但是系统耗时过长,不利于生物的快速成像,鉴于此,CB-XLCT借助锥束型射 线扫描时间短和X剂量利用率高的优点,为今后在体外快速成像提供了可行性。
类似于BLT和FMT光学断层成像,CB-XLCT成像问题在数学上属于逆问题研究, 已有相关方法可有效降低逆问题的病态性,随后,稀疏角CB-XLCT技术相继出现,然 而由于稀疏角度数据的约束,成像逆问题病态性相比传统CB-XLCT的逆问题病态性更 为严重,使得传统CB-XLCT成像方法不能有效应用于稀疏角CB-XLCT重建中。
发明内容
为克服现有技术的不足,本发明目的在于提供一种优化区域知识先验的稀疏角CB-XLCT成像方法,可有效降低逆问题的病态性,准确重建靶标。
为实现上述目的,本发明采用的技术方案是:
一种优化区域知识先验的稀疏角CB-XLCT成像方法,具体包括以下步骤;
步骤1:基于稀疏角CB-XLCT成像模型理论,构建稀疏非凸Lp,0<p<1正则模型;
所述稀疏非凸Lp正则模型的计算公式表示为:
其中,τ为正则化参数,||ρ||p=(Σii|p)1/p为Lp,0<p<1范数,ρ表示待求解的纳米靶标,J表示采集的生物体表光信号,A为系统矩阵;
步骤2:将步骤1中所述的稀疏非凸Lp正则模型转化为重加权凸L1正则模型进行 迭代求解;
步骤3:将步骤2中所述的重加权凸L1正则模型迭代求解过程中经过迭代重加权裂分增广拉格朗日收缩算法进行求解;
步骤4:将步骤3得到的求解结果结合经典非凸p值,求解出稀疏解集对应的可 行区域集;
步骤5:将步骤4中得到的可行区域集优化为稳定的知识先验指导逆问题的求解,准确重建靶标;
步骤6:根据步骤5中重建后的靶标显示稀疏角CB-XLCT图像。
进一步地,所述步骤1中CB-XLCT成像模型的计算公式为:J=Aρ;
其中,ρ表示待求解的纳米靶标,J表示采集的生物体表光信号。
进一步地,所述步骤2中的迭代过程为:
上式中,ρk+1为第k+1次迭代得到的逆问题解,和ρi分别为相应逆问题解的分量,ε为避免所在分母为零的非零常数,可利用下式自动选取:/>其中,i0=m/[4log(n/m)],n代表系统矩阵A的列数,m代表系统矩阵A的行数;
转化为:
其中/>表示权重系数;
基于变量裂分思想,子问题可转化为:
其中/>
进一步地,所述步骤3中迭代重加权裂分增广拉格朗日收缩算法求解过程为:
采用经典的交替方向乘子法交替求解 中两个裂分变量:
其中,dk+1=dk-(vk+1k+1);
v为裂分变量,dt为与拉格朗日乘子相关的迭代向量,dk+1和dk为与拉格朗日乘 子相关的迭代向量。
进一步地,所述步骤4中求解出稀疏解集对应的可行区域集的过程为:对于常用按照间隔/>进行不同的p值选取并分别带入到 />进行模型求解,接着对于稀疏解确定相应可行区 域,得到稀疏解集对应的可行区域集合,进一步对得到的可行区域集合进行求交优化。
进一步地,所述步骤5中将可行区域集的优化结果视为知识先验,定义含有0或 1的列向量算子
将上式代入式进行逆问题规模约简并
优化求解空间:
其中,表示分别从式/>中的A和ρ中去除列向量算子/>中非零元操作,如果/>中的节点处在PSR中则置1, 否则置0;对于式/>采用L0正则模型或L1正则模型求解。
与现有技术相比,本发明的有益效果是:
1)本发明通过构建稀疏非凸Lp,0<p<1模型,利用迭代重加权裂分增广拉格朗 日收缩算法,进一步结合经典非凸p值提出一种鲁邦稳定的可行区域提取方法,进而 作为优化的知识先验指导靶标的准确重建,有效地解决了稀疏角CB-XLCT成像严重不 适定的逆问题。
2)本发明采用的成像方法通过与L0正则模型和L1正则模型各自经典的算法相结合,具有良好的扩展性。
附图说明
图1是本发明的成像算法框架图;
图2(a)是仿真实验采用的数字鼠躯干模型;
图2(b)是数字鼠躯干模型在稀疏4角度下的前向仿真结果图一;
图2(c)是数字鼠躯干模型在全24角度下的前向仿真结果图二;
图3(a)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图一;
图3(b)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图二;
图3(c)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图三;
图3(d)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图四;
图3(e)是本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图一;
图3(f)是本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图二;
图3(g)是本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图三;
图3(h)是本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在稀疏4角度投影下重建结果对比图四;
图4(a)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在全24角度投影下重建结果对比图一;
图4(b)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在全24角度投影下重建结果对比图二;
图4(c)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在全24角度投影下重建结果对比图三;
图4(d)是本发明分别结合L1正则模型的两种经典算法和L0正则模型的两种经 典算法,在全24角度投影下重建结果对比图四;
图4(e)本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经典 算法,在全24角度投影下重建结果对比图(一);
图4(f)本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经典 算法,在全24角度投影下重建结果对比图(二);
图4(g)本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经典 算法,在全24角度投影下重建结果对比图(三);
图4(h)本发明直接采用L1正则模型的两种经典算法和L0正则模型的两种经典 算法,在全24角度投影下重建结果对比图(四);
图5(a)是本发明在稀疏4角度下靶标深度水平(50mm,58mm,66mm,74mm,82mm,90mm)的鲁棒性测试重建结果曲线图;
图5(b)是本发明在噪声水平(5%,10%,15%,20%,25%,30%,35%,40%)的鲁棒性测试重建 结果曲线图;
图6(a)是本发明真实物理仿体实验表面光信号分布图;
图6(b)是本发明稀疏4角度重建X-Y方向重建结果与CT切片叠加展示图一;
图6(c)是本发明稀疏4角度重建X-Z方向重建结果与CT切片叠加展示图二;
图6(d)是本发明稀疏4角度重建Y-Z方向重建结果与CT切片叠加展示图三。
具体实施方式
下面结合具体实施例对本发明作进一步详细描述,但不作为对本发明的限定。
如图1所示,本发明一种优化区域知识先验的稀疏角CB-XLCT成像方法,具体包 括以下步骤;
步骤1:CB-XLCT成像模型计算公式表示为:J=Aρ (1);
式(1)中,ρ表示待求解的纳米靶标,J表示采集的生物体表光信号,A为系 统矩阵。
生物组织经过散射、折射、吸收等过程,最终达到生物体表的光信息非常有限, 而待求解的纳米靶标ρ尾数通常很高,无法直接求解,此外,模拟早期肿瘤的纳米靶 标在生物体内具有明显的稀疏分布特性,因此,基于稀疏表示理论,构建稀疏非凸Lp (0<p<1)正则模型;如图2所示,其中2(a)表示仿真实验采用的数字鼠躯干模型, 2(b)表示数字鼠躯干模型在稀疏4角度下的前向仿真结果,2(c)表示数字鼠躯干 模型在全24角度下的前向仿真结果,稀疏非凸Lp(0<p<1)正则模型的计算公式表示 为:
式(2)中,τ为正则化参数,||ρ||p=(Σii|p)1/p为Lp(0<p<1)范数。
步骤2:上式(2)中目标函数的非凸特点,直接求解非常困难,基于迭代重加 权思想,将求解过程转化为重加权凸L1正则模型进行迭代求解,如图3所示,其中3 (a)至3(d)分别结合L1正则模型的两种经典算法和L0正则模型的两种经典算法, 在稀疏4角度投影下重建结果对比图,图3(e)至图3(h)是不结合本发明,分别 直接采用L1正则模型的两种经典算法和L0正则模型的两种经典算法,在稀疏4角度 投影下重建结果对比图;如图(4)所示,其中4(a)至4(d)分别结合L1正则模型 的两种经典算法和L0正则模型的两种经典算法,在全24角度投影下重建结果对比图, 图4(e)至图4(h)是不结合本发明,分别直接采用L1正则模型的两种经典算法和 L0正则模型的两种经典算法,在全24角度投影下重建结果对比图,迭代求解过程表 示如下:
式(3)中,ρk+1为第k+1次迭代得到的逆问题解,ρi k和ρi分别为相应逆问题解 的分量,ε为避免所在分母为零的非零常数,可利用下式自动选取: (4);式(4)中,i0=m/[4log(n/m)],n代表系统矩阵A的列数,m代表系统矩阵A的行 数。
将式(3)转化为:
式(5)中表示权重系数。
基于变量裂分思想,将上式(5)的子问题可转化为:
式(6)中/>
步骤3:通过迭代重加权裂分增广拉格朗日收缩算法求解逆问题;
采用经典的交替方向乘子法交替求解式(7)中两个裂分变量:
式(8)中dk+1=dk-(vk+1k+1);
v为裂分变量,dt为与拉格朗日乘子相关的迭代向量,dk+1和dk为与拉格朗日乘 子相关的迭代向量。
步骤4:结合经典非凸p值,求解稀疏解集对应的可行区域集,并进一步优化为 稳定的知识先验指导逆问题求解;
对于常用按照间隔/>进行不同的p值选取并分别带入到式(5)进行模型求解,接着对于稀疏解确定相应可行区域,得到稀疏解集对应的可行区域集合,进一 步对得到的可行区域集合进行求交优化。
步骤5:将可行区域集的优化结构视为知识先验指导,定义含有0或1的列向量 算子作为知识先验指导逆问题求解:
将上式(9)代入式(5)进行逆问题规模约简并优化求解空间:
式(10)中表示分别从式(5)的A和ρ中去除列向量算子中非零元操作,如果/>中的节点处在PSR中则置1,否则置0。对于式(10),可采用 常见的L0正则模型或L1正则模型来有效求解。
步骤6:根据步骤5中重建后的靶标显示稀疏角CB-XLCT图像,如图5所示,其 中5(a)表示在稀疏4角度下靶标深度水平(50mm,58mm,66mm,74mm,82mm,90mm)的鲁 棒性测试重建结果曲线图,5(b)表示在噪声水平(5%,10%,15%,20%,25%,30%,35%,40%)的鲁棒性测试重建结果曲线图;如图6所示,其中6(a)表示真实物理仿体实验表面 光信号分布图,6(b)表示稀疏4角度重建X-Y方向重建结果与CT切片叠加展示,6 (c)表示稀疏4角度重建X-Z方向重建结果与CT切片叠加展示,6(d)表示稀疏4 角度重建Y-Z方向重建结果与CT切片叠加展示。
综上所述,本发明通过构建稀疏非凸Lp(0<p<1)模型,利用迭代重加权裂分增 广拉格朗日收缩算法,进一步结合经典非凸p值提出一种鲁邦稳定的可行区域提取方 法,进而作为优化的知识先验指导靶标的准确重建,有效地解决了稀疏角CB-XLCT成 像严重不适定的逆问题。

Claims (6)

1.一种优化区域知识先验的稀疏角CB-XLCT成像方法,其特征在于,具体包括以下步骤;
步骤1:基于稀疏角CB-XLCT成像模型理论,构建稀疏非凸Lp,0<p<1正则模型;
所述稀疏非凸Lp正则模型的计算公式表示为:
其中,τ为正则化参数,||ρ||p=(Σii|p)1/p为Lp,0<p<1范数,ρ表示待求解的纳米靶标,J表示采集的生物体表光信号,A为系统矩阵;
步骤2:将步骤1中所述的稀疏非凸Lp正则模型转化为重加权凸L1正则模型进行迭代求解;
步骤3:将步骤2中所述的重加权凸L1正则模型迭代求解过程经过迭代重加权裂分增广拉格朗日收缩算法进行求解;
步骤4:将步骤3得到的求解结果结合经典非凸p值,求解出稀疏解集对应的可行区域集;
步骤5:将步骤4中得到的可行区域集优化为稳定的知识先验指导逆问题的求解,准确重建靶标;
步骤6:根据步骤5中重建后的靶标显示稀疏角CB-XLCT图像。
2.根据权利要求1所述的一种优化区域知识先验的稀疏角CB-XLCT成像方法,其特征在于:所述步骤1中CB-XLCT成像模型的计算公式为:J=Aρ;
其中,ρ表示待求解的纳米靶标,J表示采集的生物体表光信号。
3.根据权利要求1所述的一种优化区域知识先验的稀疏角CB-XLCT成像方法,其特征在于:所述步骤2中的迭代过程为:
上式中,ρk+1为第k+1次迭代得到的逆问题解,和ρi分别为相应逆问题解的分量,ε为避免所在分母为零的非零常数,可利用下式自动选取:/>其中,i0=m/[4log(n/m)],n代表系统矩阵A的列数,m代表系统矩阵A的行数;
转化为:
其中/>表示权重系数;
基于变量裂分思想,子问题可转化为:
其中/>
4.根据权利要求3所述的一种优化区域知识先验的稀疏角CB-XLCT成像方法,其特征在于:所述步骤3中迭代重加权裂分增广拉格朗日收缩算法求解过程为:
采用经典的交替方向乘子法交替求解
中两个裂分变量:
其中,dk+1=dk-(vk+1k+1);
v为裂分变量,dt为与拉格朗日乘子相关的迭代向量,dk+1和dk为与拉格朗日乘子相关的迭代向量。
5.根据权利要求4所述的一种优化区域知识先验的稀疏角CB-XLCT成像方法,其特征在于:所述步骤4中求解出稀疏解集对应的可行区域集的过程为:对于常用按照间隔/>进行不同的p值选取并分别带入到/>进行模型求解,接着对于稀疏解确定相应可行区域,得到稀疏解集对应的可行区域集合,进一步对得到的可行区域集合进行求交优化。
6.根据权利要求5所述的一种优化区域知识先验的稀疏角CB-XLCT成像方法,其特征在于:所述步骤5中将可行区域集的优化结果视为知识先验,定义含有0或1的列向量算子
将上式代入式进行逆问题规模约简并优化求解空间:
其中, 表示分别从式/>中的A和ρ中去除列向量算子/>中非零元操作,如果/>中的节点处在PSR中则置1,否则置0;对于式/>采用L0正则模型或L1正则模型求解。
CN202110429169.9A 2021-04-21 2021-04-21 一种优化区域知识先验的稀疏角cb-xlct成像方法 Active CN113112563B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110429169.9A CN113112563B (zh) 2021-04-21 2021-04-21 一种优化区域知识先验的稀疏角cb-xlct成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110429169.9A CN113112563B (zh) 2021-04-21 2021-04-21 一种优化区域知识先验的稀疏角cb-xlct成像方法

Publications (2)

Publication Number Publication Date
CN113112563A CN113112563A (zh) 2021-07-13
CN113112563B true CN113112563B (zh) 2023-10-27

Family

ID=76719029

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110429169.9A Active CN113112563B (zh) 2021-04-21 2021-04-21 一种优化区域知识先验的稀疏角cb-xlct成像方法

Country Status (1)

Country Link
CN (1) CN113112563B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201410528D0 (en) * 2013-06-14 2014-07-30 Nuctech Co Ltd CT imaging methods and systems
CN107016653A (zh) * 2017-03-29 2017-08-04 中国人民解放军信息工程大学 基于总曲率联合总变分的ct图像稀疏角度重建方法及装置
CN107133930A (zh) * 2017-04-30 2017-09-05 天津大学 基于低秩矩阵重建和稀疏表示的行列缺失图像填充方法
WO2018099321A1 (zh) * 2016-11-30 2018-06-07 华南理工大学 一种基于广义树稀疏的权重核范数磁共振成像重建方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201410528D0 (en) * 2013-06-14 2014-07-30 Nuctech Co Ltd CT imaging methods and systems
WO2018099321A1 (zh) * 2016-11-30 2018-06-07 华南理工大学 一种基于广义树稀疏的权重核范数磁共振成像重建方法
CN107016653A (zh) * 2017-03-29 2017-08-04 中国人民解放军信息工程大学 基于总曲率联合总变分的ct图像稀疏角度重建方法及装置
CN107133930A (zh) * 2017-04-30 2017-09-05 天津大学 基于低秩矩阵重建和稀疏表示的行列缺失图像填充方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
具有约束的稀疏正则图像重建模型及其在CT成像中的应用;唐玉超;陈宝;朱传喜;;南昌大学学报(工科版)(第03期);全文 *
实现稀疏角度下的精确CT重建:利用ADMM-LP算法求解非凸模型;宋洁;陈平;潘晋孝;;中国组织工程研究(第31期);全文 *

Also Published As

Publication number Publication date
CN113112563A (zh) 2021-07-13

Similar Documents

Publication Publication Date Title
CN111492406B (zh) 训练机器学习算法的方法、图像处理系统和图像重建方法
Kuanar et al. Low dose abdominal CT image reconstruction: an unsupervised learning based approach
US9613441B2 (en) Joint reconstruction of electron density images
CN102488493B (zh) 小动物活体多模态分子成像系统和成像方法
CN103514629B (zh) 用于迭代重建的方法和设备
Li et al. DenseX-net: an end-to-end model for lymphoma segmentation in whole-body PET/CT images
Zhang et al. UHR-DeepFMT: ultra-high spatial resolution reconstruction of fluorescence molecular tomography based on 3-D fusion dual-sampling deep neural network
Fu et al. A deep learning reconstruction framework for differential phase-contrast computed tomography with incomplete data
Zhang et al. PET image reconstruction using a cascading back-projection neural network
CN111915733A (zh) 基于LeNet网络的三维锥束X射线发光断层成像方法
He et al. Downsampled imaging geometric modeling for accurate CT reconstruction via deep learning
US11350895B2 (en) System and method for spectral computed tomography using single polychromatic x-ray spectrum acquisition
Ohnishi et al. Three-dimensional vessel segmentation in whole-tissue and whole-block imaging using a deep neural network: Proof-of-concept study
US20240135603A1 (en) Metal Artifact Reduction Algorithm for CT-Guided Interventional Procedures
Zhang et al. Self-training strategy based on finite element method for adaptive bioluminescence tomography reconstruction
Yang et al. Accelerating vasculature imaging in tumor using mesoscopic fluorescence molecular tomography via a hybrid reconstruction strategy
Guo et al. Spectral2Spectral: Image-spectral similarity assisted deep spectral CT reconstruction without reference
CN113112563B (zh) 一种优化区域知识先验的稀疏角cb-xlct成像方法
Germer et al. Limited-angle tomography reconstruction via deep end-to-end learning on synthetic data
CN116712038A (zh) 基于螺旋交错稀疏采样的多光谱光声断层成像系统及方法
Wu et al. Deep learning-based low-dose tomography reconstruction with hybrid-dose measurements
Gao et al. Self-Absorption Correction in X-Ray Fluorescence-Computed Tomography With Deep Convolutional Neural Network
CN115423892A (zh) 一种基于最大期望网络的无衰减校正pet重建方法
CN113288188A (zh) 基于分组注意力残差网络的锥束x射线发光断层成像方法
Ren et al. High-resolution tomographic reconstruction of optical absorbance through scattering media using neural fields

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant