磁共振弥散张量成像方法和系统.pdf

上传人:xia****o6 文档编号:6630793 上传时间:2019-09-03 格式:PDF 页数:16 大小:853.04KB
返回 下载 相关 举报
摘要
申请专利号:

CN201010612326.1

申请日:

20101229

公开号:

CN102008307B

公开日:

20120725

当前法律状态:

有效性:

有效

法律详情:

IPC分类号:

A61B5/055

主分类号:

A61B5/055

申请人:

中国科学院深圳先进技术研究院

发明人:

吴垠,戴睿彬,刘新,郑海荣,邱本胜,邹超,张娜,谢国喜

地址:

518055 广东省深圳市南山区西丽大学城学苑大道1068号

优先权:

CN201010612326A

专利代理机构:

广州华进联合专利商标代理有限公司

代理人:

吴平

PDF下载: PDF下载
内容摘要

一种磁共振弥散张量成像方法,包括以下步骤:对k空间进行采集,并获取导航数据和测量数据;根据所述导航数据获取时间基函数和频率成分参数,以及根据测量数据获取空间基函数;根据所述时间基函数、频率成分参数和空间基函数计算获得k-t空间数据;根据所述k-t空间数据进行傅里叶反变换并获得重建图像。基于部分可分离函数技术,加快获取扫描数据,快速获取运动物体的成像信息,达到快速成像。

权利要求书

1.一种磁共振弥散张量成像方法,包括以下步骤:对k空间进行采集,并获取导航数据和测量数据;从K空间采集的数据S(k,t)基于可分离函数技术,表示为:所述L为频率成分参数、c(k)为空间基函数及为时间基函数;根据所述导航数据获取时间基函数和频率成分参数,以及根据测量数据获取空间基函数;所述导航数据确定时间基函数和频率成分参数,该测量数据确定空间基函数;所采集的点需要满足3个条件:1.T必须满足S(k,t)的时间奈奎斯特速率,T是重复时间,S(k,t)是导航数据;2.Δk必须满足S(k,t)的时间奈奎斯特速率,Δk是相邻相位编码线的k空间间距;3.从S(k,t)获取的采样帧数N必须大于或等于阶数L,S(k,t)是测量数据,帧数N表示不同时刻获得的图像的总数;根据所述时间基函数、频率成分参数和空间基函数计算获得k-t空间数据;根据所述k-t空间数据进行傅里叶反变换并获得重建图像。 2.根据权利要求1所述的磁共振弥散张量成像方法,其特征在于,在对k空间进行采集,并获取导航数据和测量数据的步骤之前还包括:心电触发获取采集信号;所述心电触发获取采集信号的步骤包括:采集心电图信号并确定心动周期;获取心脏动态成像,且根据所述心脏动态成像确定在心动周期内的延迟时间及准稳态期;在所述准稳态期内施加弥散梯度脉冲且进行数据采集。 3.根据权利要求1所述的磁共振弥散张量成像方法,其特征在于,在对k空间进行采集,并获取导航数据和测量数据的步骤之前还包括:呼吸导航获取采集信号;所述呼吸导航获取采集信号的步骤包括:发射二维脉冲信号;根据所述二维脉冲信号获取导航信号;根据所述导航信号获取呼吸运动幅度值;判断所述呼吸运动幅度值是否超出预设的导航窗,是,则对所述呼吸运动幅度值所对应的时间区间的采集信号进行摒弃,否,则所述呼吸运动幅度值所对应的时间区间的采集信号进行采集。 4.根据权利要求1所述的磁共振弥散张量成像方法,其特征在于,在对k空间进行采集,并获取导航数据和测量数据的步骤之前还包括:心电触发获取采集信号及呼吸导航获取采集信号的步骤;根据心电触发获取采集信号获取准稳态时间,根据所述导航获取采集信号获取在预设导航窗内的呼吸运动幅度值;在所述准稳态期内以及符合所述呼吸运动幅度值的时间区间进行信号采集。 5.一种磁共振弥散张量成像系统,其特征在于,包括:采集模块,用于对k空间进行采集,并获取导航数据和测量数据;从K空间采集的数据S(k,t)基于可分离函数技术,表示为:所述L为频率成分参数、c(k)为空间基函数及为时间基函数;提取模块,用于根据所述导航数据获取时间基函数和频率成分参数,以及根据测量数据获取空间基函数;所述导航数据确定时间基函数和频率成分参数,该测量数据确定空间基函数;所采集的点需要满足3个条件:1.T必须满足S(k,t)的时间奈奎斯特速率,T是重复时间,S(k,t)是导航数据;2.Δk必须满足S(k,t)的时间奈奎斯特速率,Δk是相邻相位编码线的k空间间距;3.从S(k,t)获取的采样帧数N必须大于或等于阶数L,S(k,t)是测量数据,帧数N表示不同时刻获得的图像的总数;计算模块,用于根据所述时间基函数、频率成分参数和空间基函数计算获得k-t空间数据;重建模块,根据所述k-t空间数据进行傅里叶反变换并获得重建图像。 6.根据权利要求5所述的磁共振弥散张量成像系统,其特征在于,还包括:与采集模块连接的心电触发模块,所述心电触发模块包括:心电采集单元,用于采集心电图信号并确定心动周期;获取单元,用于获取心脏动态成像,且根据所述心脏动态成像确定在心动周期内的延迟时间及准稳态期;心电执行单元,用于在准稳态期内施加弥散梯度脉冲且进行数据采集。 7.根据权利要求5所述的磁共振弥散张量成像系统,其特征在于,还包括:与采集模块连接的呼吸导航模块,所述呼吸导航模块包括:发射单元,用于发射二维脉冲信号;接收单元,用于根据所述二维脉冲信号获取导航信号;呼吸运动幅度值单元,用于根据所述导航信号获取呼吸运动幅度值;处理单元,用于判断所述呼吸运动幅度值是否超出预设的导航窗,是,则对所述呼吸运动幅度值所对应的时间区间的采集信号进行摒弃,否,则所述呼吸运动幅度值所对应的时间区间的采集信号进行采集。

说明书



【技术领域】

本发明涉及磁共振技术,特别是涉及一种磁共振弥散张量成像方法和系统。

【背景技术】

在发展中国家,心肌梗塞是最大的死亡原因之一。在我国,近年来心肌梗 塞的发病率呈明显上升趋势,占心血管疾病的一半以上,成为卫生保健和卫生 资源的沉重负担。心肌结构的改变是心肌梗塞引发心脏衰竭的最主要原因之一。 因此,对于心肌结构的研究,不仅能够从微观结构角度了解心脏功能的机理, 也可以从细胞水平上的变化发现病症,并为诊断和治疗提供相应的依据。

目前,传统的成像技术,其获取数据的速度慢,而且对于运动十分敏感, 导致带来的图像伪影十分严重。

【发明内容】

基于此,有必要提供一种快速的磁共振弥散张量成像方法。

另外,还有必要提供一种快速的磁共振弥散张量成像系统。

一种磁共振弥散张量成像方法,包括以下步骤:对k空间进行采集,并获 取导航数据和测量数据;从K空间采集的数据S(k,t)基于可分离函数技术,表示 为:所述L为频率成分参数、cl(k)为空间基函数及为时间 基函数;根据所述导航数据获取时间基函数和频率成分参数,以及根据测量数 据获取空间基函数;所述导航数据确定时间基函数和频率成分参数,该测量数 据确定空间基函数;所采集的点需要满足3个条件:1.TR必须满足Snav(k,t)的时 间奈奎斯特速率,TR是重复时间,Snav(k,t)是导航数据;2.Δky必须满足Simg(k,t) 的时间奈奎斯特速率,Δky是相邻相位编码线的k空间间距;3.从Simg(k,t)获取 的采样帧数N必须大于或等于阶数L,Simg(k,t)是测量数据,帧数N表示不同时 刻获得的图像的总数;根据所述时间基函数、频率成分参数和空间基函数计算 获得k-t空间数据;根据所述k-t空间数据进行傅里叶反变换并获得重建图像。

优选地,在对k空间进行采集,并获取导航数据和测量数据的步骤之前还 包括:心电触发获取采集信号;所述心电触发获取采集信号的步骤包括:采集 心电图信号并确定心动周期;获取心脏动态成像,且根据所述心脏动态成像确 定在心动周期内的延迟时间及准稳态期;在所述准稳态期内施加弥散梯度脉冲 且进行数据采集。

优选地,在对k空间进行采集,并获取导航数据和测量数据的步骤之前还 包括:呼吸导航获取采集信号;所述呼吸导航获取采集信号的步骤包括:发射 二维脉冲信号;根据所述二维脉冲信号获取导航信号;根据所述导航信号获取 呼吸运动幅度值;判断所述呼吸运动幅度值是否超出预设的导航窗,是,则对 所述呼吸运动幅度值所对应的时间区间的采集信号进行摒弃,否,则所述呼吸 运动幅度值所对应的时间区间的采集信号进行采集。

优选地,在对k空间进行采集,并获取导航数据和测量数据的步骤之前还 包括:心电触发获取采集信号及呼吸导航获取采集信号的步骤;根据心电触发 获取采集信号获取准稳态时间,根据所述导航获取采集信号获取在预设导航窗 内的呼吸运动幅度值;在所述准稳态期内以及符合所述呼吸运动幅度值的时间 区间进行信号采集。

还有必要提供一种磁共振弥散张量成像系统,包括:采集模块,用于对k 空间进行采集,并获取导航数据和测量数据;从K空间采集的数据S(k,t)基于可 分离函数技术,表示为:所述L为频率成分参数、cl(k)为空间 基函数及为时间基函数;提取模块,用于根据所述导航数据获取时间基函数 和频率成分参数,以及根据测量数据获取空间基函数;所述导航数据确定时间 基函数和频率成分参数,该测量数据确定空间基函数;所采集的点需要满足3 个条件:1.TR必须满足Snaw(k,t)的时间奈奎斯特速率,TR是重复时间,Snaw(k,t)是 导航数据;2.Δky必须满足Simg(k,t)的时间奈奎斯特速率,Δky是相邻相位编码线 的k空间间距;3.从Simg(k,t)获取的采样帧数N必须大于或等于阶数L,Simg(k,t) 是测量数据,帧数N表示不同时刻获得的图像的总数;计算模块,用于根据所 述时间基函数、频率成分参数和空间基函数计算获得k-t空间数据;重建模块, 根据所述k-t空间数据进行傅里叶反变换并获得重建图像。

优选地,还包括:与采集模块连接的心电触发模块,所述心电触发模块包 括:心电采集单元,用于采集心电图信号并确定心动周期;获取单元,用于获 取心脏动态成像,且根据所述心脏动态成像确定在心动周期内的延迟时间及准 稳态期;心电执行单元,用于在准稳态期内施加弥散梯度脉冲且进行数据采集。

优选地,还包括:与采集模块连接的呼吸导航模块,所述呼吸导航模块包 括:发射单元,用于发射二维脉冲信号;接收单元,用于根据所述二维脉冲信 号获取导航信号;呼吸运动幅度值单元,用于根据所述导航信号获取呼吸运动 幅度值;处理单元,用于判断所述呼吸运动幅度值是否超出预设的导航窗,是, 则对所述呼吸运动幅度值所对应的时间区间的采集信号进行摒弃,否,则所述 呼吸运动幅度值所对应的时间区间的采集信号进行采集。

采用上述磁共振弥散张量成像方法及系统,基于部分可分离函数技术,加 快获取扫描数据,快速获取运动物体的成像信息,达到快速成像。

【附图说明】

图1为磁共振弥散张量成像方法的流程图;

图2为一实施例的基于部分可分离函数技术的数据采集示意图;

图3为心电触发获取采集信号的流程图;

图4为一实施例的心电触发获取采集信号的心电触发成像原理图;

图5为一实施例的心电触发获取采集信号的一个心动周期内心脏解剖结构 “收缩-舒张-收缩”示意图;

图6为一实施例的心电触发获取采集信号的同一区域的心脏剖面信号连接 图;

图7为一实施例的呼吸导航获取采集信号流程图;

图8为一实施例的呼吸导航获取采集信号中放置于隔肌上的发射二维脉冲 信号的示意图;

图9为一实施例的同时采用心电触发获取采集信号及呼吸导航获取采集信 号的方法流程图;

图10为一实施例的同时采用心电触发及呼吸导航获取采集信号的示意图;

图11为磁共振弥散张量成像装置的原理框图;

图12为磁共振弥散张量成像装置中采用心电触发模块的原理框图;

图13为磁共振弥散张量成像装置中采用呼吸导航模块的原理框图。

【具体实施方式】

对于运动物体的磁共振成像过程中,图像是在其空间位置和时间t的函数。 若磁共振的扫描速度足够快,在物体的空间位置还来不及发生很大变化的瞬间 就可以采集重建图像所需要的全部数据,则可近似认为信号扫描的瞬间物体处 于静止状态。但是,目前的磁共振设备无法达到如此高的扫描速度。

基于此,很难在相对静止的瞬间采集足够的数据进行图像重建。此时所采 集的磁共振信号与自旋质子密度实际是空间位置和时间t的函数,其中接收信 号S(k,t)和所需要的图像函数ρ(r,t)之间的关系为:

S ( k , t ) = ∫ - ∞ + ∞ ρ ( r , t ) e - i 2 πk · r dt ]]>

请参阅附图1和2,一种磁共振弥散张量成像方法,包括以下步骤:

S10:对k空间进行采集,并获取导航数据和测量数据。

具体地,经过一段时间对k空间的信号进行采集,此时空间中某一平面数 据将包含来自运动物体各个不同部位的信息。

对该采集的数据基于部分可分离函数技术,即认为图像函数ρ(r,t)的空间变 化和时间变化是L阶可分离的(阶数与成像物体的运动模式相关,理论上阶数 越大,对运动的描述越精确),利用部分可分离函数的性质,S(k,t)可以表示为 空间和时间两个独立变量函数cl(k)和

所述L为频率成分参数、cl(k)为空间基函数及为时间基函数。

S20:根据导航数据获取时间基函数和频率成分参数,以及根据测量数据获 取空间基函数。

具体地,为了获得频率成分参数L、空间基函数cl(k)及时间基函数只 需要采集两组k空间数据即可。即对k空间进行采集,获取高时间、低空间分 辨率导航数据(Navigator data,空心圆圈Snav(k,t),见图2),获取高时间、低空 间分辨率测量数据(Measurement data,实心圆点Simg(k,t),见图2)。该导航数据 确定和L,该测量数据确定cl(k)。

而且,所采集的点需要满足3个条件:

1.TR必须满足Snav(k,t)的时间奈奎斯特速率,TR是重复时间,Snav(k,t)是导航 数据;

2.Δky必须满足Simg(k,t)的时间奈奎斯特速率,Δky是相邻相位编码线的k空 间间距;

3.从Simg(k,t)获取的采样帧数N必须大于或等于阶数L,Simg(k,t)是测量数据, 帧数N表示不同时刻获得的图像的总数。

S30:根据时间基函数、频率成分参数和空间基函数计算获得k-t空间数据。 根据获得的时间基函数、频率成分参数和空间基函数后,则可以计算出扩大的 k-t空间数据(Synthetic data,实心叉,见图2)。

S40:根据所述k-t空间数据进行傅里叶反变换并获得重建图像。

基于部分可分离函数技术,加快获取扫描数据,快速获取运动物体的成像 信息,达到快速成像。

在一实施例中,结合附图3~6,在步骤S10之前还包括心电触发获取采集信 号的步骤A。在弥散梯度脉冲施加期间,物体中的水分子在弥散梯度方向上的 布朗运动将积累净相位,在成像时与静止的分子相比信号强度减弱。可见,运 动的物体在相同的位置上施加弥散梯度脉冲,才能够保证静止的分子净相位累 计为0。

具体地,心电触发获取采集信号的步骤A包括:

A1:采集心电图信号并确定心动周期。具体地,采集心电图(ECG)R波 为信号测量的参考点,并确定心脏的心动周期。

A2:获取心脏动态成像,且根据所述心脏动态成像确定在心动周期内的延 迟时间及准稳态期。具体地,通过施加一种磁共振脉冲序列,例如平衡式稳态 自由进动梯度回波(Balance-FFE)序列,根据该序列快速获得心脏动态成像。然 后根据心脏动态成像确定在心动周期内的延迟时间及准稳态期。

对于延迟时间(Delay time)D选取是否能够达到预计的效果,即保证运动 物体的相对位移值为零,现选取一具体实施例进行详细说明。

按上述方法,首先通过施加一种磁共振脉冲序列获得心脏动态成像。然后, 在一个心动周期内捕获20个或以上心脏相位对应的心脏解剖结构,包括从“舒 张态-收缩态-舒张态”的周而复始的变化,根据心内膜和心外膜的轮廓的位移判 断心脏运动情况(见附图5)。结合附图6,首先,放置一条跨越左心室中心点 的采值线(涵盖心肌和左心室腔),则可提取该线所覆盖区域的信号强度,从而 获得相应的信号轮廓。然后,将20个心脏相位上同一区域的信号连接起来,可 以得出在一个心动周期内心肌及左心室腔的变化情况。不难看出,在白色线段 标注范围内,心肌和左心室腔基本没有明显的变化,表明在该时间段内心脏运 动的变化不大。可见在该时间段,心脏跳动是比较舒缓的,因此可以据此确定 延迟时间(Delay time)D。在延迟时间D后,且在一个心动周期内,确定为准 稳态期。

A3:在准稳态期内施加弥散梯度脉冲且进行数据采集。具体地,在该准稳 态期内发射弥散张量成像(DTI)的弥散梯度脉冲,并进行数据采集。以至于所 获取的采集数据更加精确,更少受到心脏运动的影响。

采用该方案,通过对心电图信号以及心脏动态成像确定心动周期及准稳态 期,进而在准稳态期内施加弥散梯度脉冲且进行数据采集,准确获得成像数据, 减少运动伪影。

在另一实施例中,结合附图7~8,在步骤S10之前还包括呼吸导航获取采集 信号的步骤B。采用呼吸导航技术,主要是用于检测自由呼吸下隔面位置的变化, 根据隔面的位置信息确定采集信号。

具体地,呼吸导航获取采集信号的步骤B包括:

B1:发射二维脉冲信号。具体地,在隔肌之上(见图8)发射二维脉冲信号。

B2:根据二维脉冲信号获取导航信号。具体地,根据二维脉冲信号获取导 航笔形波束的导航信号。

B3:根据导航信号获取呼吸运动幅度值。具体地,在整个扫描过程中回波 信号代表了该时间段内隔肌的运动情况,膈肌的运动同时反映了呼吸运动,故 可以获得呼吸运动幅度值。

B4:判断呼吸运动幅度值是否超出预设的导航窗,是,则对呼吸运动幅度 值所对应的时间区间的采集信号进行摒弃,否,则所述呼吸运动幅度值所对应 的时间区间的采集信号进行采集。具体地,预设一导航窗,当呼吸运动的幅值 在导航窗的范围内,该范围内所对应的时间区间则认为此数据是有效且采集该 数据。当呼吸运动的幅值在导航窗的范围外,该范围内所对应的时间区间则认 为该呼吸运动幅度较大,应该摒弃在该时间所采集的数据。

采用该实施例的方案,所采集的数据为运动物体(例如在体心脏)较为平 稳缓和处,即该呼吸运动较为平缓处,则所获取的数据准确,在很大程度上去 除运动伪影的干扰。

在其它实施例,请参阅附图9~10,同时采用步骤A和步骤B的方法,即在 步骤S10之前还包括步骤C,包括:

C1:根据心电触发获取采集信号获取准稳态时间,根据所述呼吸导航获取 采集信号获取在预设导航窗内的呼吸运动幅度值。

C2:在准稳态期内以及符合所述呼吸运动幅度值所对应的时间区间施加弥 散梯度脉冲且进行数据采集。

结合附图10,采用该实施例方案,获取了运动较为平缓处以及呼吸幅度值 较为平稳处的信号,则所采集的数据更为准确,在很大程度上去除运动伪影的 干扰。

请参阅附图11,基于上述的磁共振弥散张量成像方法还提供一种磁共振弥 散张量成像系统。

一种磁共振弥散张量成像系统,包括:

采集模块10,用于对k空间进行采集,并获取导航数据和测量数据。

提取模块20,用于根据所述导航数据获取时间基函数和频率成分参数,以 及根据测量数据获取空间基函数。

计算模块30,用于根据所述时间基函数、频率成分参数和空间基函数计算 获得k-t空间数据。

重建模块40,根据所述k-t空间数据进行傅里叶反变换并获得重建图像。

请参阅附图12,在一实施例中,该磁共振弥散张量成像系统还包括:与采 集模块连接的心电触发模块50,该心电触发模块50包括:

心电采集单元51,用于采集心电图信号并确定心动周期。

获取单元52,用于获取心脏动态成像,且根据所述心脏动态成像确定在心 动周期内的延迟时间及准稳态期。

心电执行单元53,用于在准稳态期内施加弥散梯度脉冲且进行数据采集。

请参阅附图13,在一实施例中,该磁共振弥散张量成像系统还包括:与采 集模块连接的呼吸导航模块60,该呼吸导航模块60包括:

发射单元61,用于发射二维脉冲信号。

接收单元62,用于根据所述二维脉冲信号获取导航信号。

呼吸运动幅度值单元63,用于根据所述导航信号获取呼吸运动幅度值。

处理单元64,用于判断所述呼吸运动幅度值是否超出预设的导航窗,是, 则对所述呼吸运动幅度值所对应的时间区间的采集信号进行摒弃,否,则所述 呼吸运动幅度值所对应的时间区间的采集信号进行采集。

基于新型的磁共振弥散张量成像系统,加快获取扫描数据,快速获取运动 物体的成像信息。

以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细, 但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域 的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和 改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附 权利要求为准。

磁共振弥散张量成像方法和系统.pdf_第1页
第1页 / 共16页
磁共振弥散张量成像方法和系统.pdf_第2页
第2页 / 共16页
磁共振弥散张量成像方法和系统.pdf_第3页
第3页 / 共16页
点击查看更多>>
资源描述

《磁共振弥散张量成像方法和系统.pdf》由会员分享,可在线阅读,更多相关《磁共振弥散张量成像方法和系统.pdf(16页珍藏版)》请在专利查询网上搜索。

1、(10)授权公告号 CN 102008307 B (45)授权公告日 2012.07.25 CN 102008307 B *CN102008307B* (21)申请号 201010612326.1 (22)申请日 2010.12.29 A61B 5/055(2006.01) (73)专利权人 中国科学院深圳先进技术研究院 地址 518055 广东省深圳市南山区西丽大学 城学苑大道 1068 号 (72)发明人 吴垠 戴睿彬 刘新 郑海荣 邱本胜 邹超 张娜 谢国喜 (74)专利代理机构 广州华进联合专利商标代理 有限公司 44224 代理人 吴平 CN 101384215 A,2009.03.。

2、11, 全文 . US 20090267603 A1,2009.10.29, 全文 . CN 101666865 A,2010.03.10, 全文 . CN 101283288 A,2008.10.08, 全文 . US 20070164740 A,2007.07.19, 全文 . WO 2004059335 A1,2004.07.15, 全文 . 郭睿, 等 . 磁共振弥散张量成像在中枢神经 系统的应用研究进展 .医学影像学杂志 .2009, 第 19 卷 ( 第 6 期 ),762-765. Alexander Karaus,et.alBlack-Blood Imaging of the 。

3、Human Heart Using Rapid Stimulated Echo Acquisition Mode (STEAM) MRI.Journal of Magnetic Resonance Imaging .2007,1666-1671. (54) 发明名称 磁共振弥散张量成像方法和系统 (57) 摘要 一种磁共振弥散张量成像方法, 包括以下步 骤 : 对 k 空间进行采集, 并获取导航数据和测量数 据 ; 根据所述导航数据获取时间基函数和频率成 分参数, 以及根据测量数据获取空间基函数 ; 根 据所述时间基函数、 频率成分参数和空间基函数 计算获得 k-t 空间数据 ; 根据所述 k。

4、-t 空间数据 进行傅里叶反变换并获得重建图像。基于部分可 分离函数技术, 加快获取扫描数据, 快速获取运动 物体的成像信息, 达到快速成像。 (51)Int.Cl. (56)对比文件 审查员 沈研研 权利要求书 2 页 说明书 6 页 附图 7 页 (19)中华人民共和国国家知识产权局 (12)发明专利 权利要求书 2 页 说明书 6 页 附图 7 页 1/2 页 2 1. 一种磁共振弥散张量成像方法, 包括以下步骤 : 对 k 空间进行采集, 并获取导航数据和测量数据 ; 从 K 空间采集的数据 S(k, t) 基于可分离函数技术, 表示为 : 所述 L 为频率成分参数、 cl(k) 为空。

5、间基函数及为时间基函数 ; 根据所述导航数据获取时间基函数和频率成分参数, 以及根据测量数据获取空间基函 数 ; 所述导航数据确定时间基函数和频率成分参数, 该测量数据确定空间基函数 ; 所采集 的点需要满足 3 个条件 : 1.TR必须满足 Snav(k, t) 的时间奈奎斯特速率, TR是重复时间, Snav(k, t) 是导航数据 ; 2.ky必须满足 Simg(k, t) 的时间奈奎斯特速率, ky是相邻相位编码线的 k 空间间 距 ; 3. 从 Simg(k, t) 获取的采样帧数 N 必须大于或等于阶数 L, Simg(k, t) 是测量数据, 帧数 N 表示不同时刻获得的图像的总。

6、数 ; 根据所述时间基函数、 频率成分参数和空间基函数计算获得 k-t 空间数据 ; 根据所述 k-t 空间数据进行傅里叶反变换并获得重建图像。 2. 根据权利要求 1 所述的磁共振弥散张量成像方法, 其特征在于, 在对 k 空间进行采 集, 并获取导航数据和测量数据的步骤之前还包括 : 心电触发获取采集信号 ; 所述心电触发获取采集信号的步骤包括 : 采集心电图信号并确定心动周期 ; 获取心脏动态成像, 且根据所述心脏动态成像确定在心动周期内的延迟时间及准稳态 期 ; 在所述准稳态期内施加弥散梯度脉冲且进行数据采集。 3. 根据权利要求 1 所述的磁共振弥散张量成像方法, 其特征在于, 在对。

7、 k 空间进行采 集, 并获取导航数据和测量数据的步骤之前还包括 : 呼吸导航获取采集信号 ; 所述呼吸导航获取采集信号的步骤包括 : 发射二维脉冲信号 ; 根据所述二维脉冲信号获取导航信号 ; 根据所述导航信号获取呼吸运动幅度值 ; 判断所述呼吸运动幅度值是否超出预设的导航窗, 是, 则对所述呼吸运动幅度值所对 应的时间区间的采集信号进行摒弃, 否, 则所述呼吸运动幅度值所对应的时间区间的采集 信号进行采集。 4. 根据权利要求 1 所述的磁共振弥散张量成像方法, 其特征在于, 在对 k 空间进行采 集, 并获取导航数据和测量数据的步骤之前还包括 : 心电触发获取采集信号及呼吸导航获取采集信。

8、号的步骤 ; 根据心电触发获取采集信号获取准稳态时间, 根据所述导航获取采集 权 利 要 求 书 CN 102008307 B 2 2/2 页 3 信号获取在预设导航窗内的呼吸运动幅度值 ; 在所述准稳态期内以及符合所述呼吸运动幅度值的时间区间进行信号采集。 5. 一种磁共振弥散张量成像系统, 其特征在于, 包括 : 采集模块, 用于对 k 空间进行采集, 并获取导航数据和测量数据 ; 从 K 空间采集的数据 S(k, t) 基于可分离函数技术, 表示为 : 所述L为频率成分参数、 cl(k)为空间基函数及为时间基函数 ; 提取模块, 用于根据 所述导航数据获取时间基函数和频率成分参数, 以及。

9、根据测量数据获取空间基函数 ; 所述导航数据确定时间基函数和频率成分参数, 该测量数据确定空间基函数 ; 所采集 的点需要满足 3 个条件 : 1.TR必须满足 Snav(k, t) 的时间奈奎斯特速率, TR是重复时间, Snav(k, t) 是导航数据 ; 2.ky必须满足 Simg(k, t) 的时间奈奎斯特速率, ky是相邻相位编码线的 k 空间间 距 ; 3. 从 Simg(k, t) 获取的采样帧数 N 必须大于或等于阶数 L, Simg(k, t) 是测量数据, 帧数 N 表示不同时刻获得的图像的总数 ; 计算模块, 用于根据所述时间基函数、 频率成分参数和空间基函数计算获得 k。

10、-t 空间 数据 ; 重建模块, 根据所述 k-t 空间数据进行傅里叶反变换并获得重建图像。 6. 根据权利要求 5 所述的磁共振弥散张量成像系统, 其特征在于, 还包括 : 与采集模块 连接的心电触发模块, 所述心电触发模块包括 : 心电采集单元, 用于采集心电图信号并确定心动周期 ; 获取单元, 用于获取心脏动态成像, 且根据所述心脏动态成像确定在心动周期内的延 迟时间及准稳态期 ; 心电执行单元, 用于在准稳态期内施加弥散梯度脉冲且进行数据采集。 7. 根据权利要求 5 所述的磁共振弥散张量成像系统, 其特征在于, 还包括 : 与采集模块 连接的呼吸导航模块, 所述呼吸导航模块包括 : 。

11、发射单元, 用于发射二维脉冲信号 ; 接收单元, 用于根据所述二维脉冲信号获取导航信号 ; 呼吸运动幅度值单元, 用于根据所述导航信号获取呼吸运动幅度值 ; 处理单元, 用于判断所述呼吸运动幅度值是否超出预设的导航窗, 是, 则对所述呼吸运 动幅度值所对应的时间区间的采集信号进行摒弃, 否, 则所述呼吸运动幅度值所对应的时 间区间的采集信号进行采集。 权 利 要 求 书 CN 102008307 B 3 1/6 页 4 磁共振弥散张量成像方法和系统 【技术领域】 0001 本发明涉及磁共振技术, 特别是涉及一种磁共振弥散张量成像方法和系统。 【背景技术】 0002 在发展中国家, 心肌梗塞是最。

12、大的死亡原因之一。 在我国, 近年来心肌梗塞的发病 率呈明显上升趋势, 占心血管疾病的一半以上, 成为卫生保健和卫生资源的沉重负担。 心肌 结构的改变是心肌梗塞引发心脏衰竭的最主要原因之一。 因此, 对于心肌结构的研究, 不仅 能够从微观结构角度了解心脏功能的机理, 也可以从细胞水平上的变化发现病症, 并为诊 断和治疗提供相应的依据。 0003 目前, 传统的成像技术, 其获取数据的速度慢, 而且对于运动十分敏感, 导致带来 的图像伪影十分严重。 【发明内容】 0004 基于此, 有必要提供一种快速的磁共振弥散张量成像方法。 0005 另外, 还有必要提供一种快速的磁共振弥散张量成像系统。 0。

13、006 一种磁共振弥散张量成像方法, 包括以下步骤 : 对 k 空间进行采集, 并获取 导航数据和测量数据 ; 从 K 空间采集的数据 S(k, t) 基于可分离函数技术, 表示为 : 所述L为频率成分参数、 cl(k)为空间基函数及为时间基函数 ; 根据 所述导航数据获取时间基函数和频率成分参数, 以及根据测量数据获取空间基函数 ; 所述 导航数据确定时间基函数和频率成分参数, 该测量数据确定空间基函数 ; 所采集的点需要 满足 3 个条件 : 1.TR必须满足 Snav(k, t) 的时间奈奎斯特速率, TR是重复时间, Snav(k, t) 是 导航数据 ; 2.ky必须满足 Simg(。

14、k, t) 的时间奈奎斯特速率, ky是相邻相位编码线的 k 空 间间距 ; 3. 从 Simg(k, t) 获取的采样帧数 N 必须大于或等于阶数 L, Simg(k, t) 是测量数据, 帧数 N 表示不同时刻获得的图像的总数 ; 根据所述时间基函数、 频率成分参数和空间基函 数计算获得 k-t 空间数据 ; 根据所述 k-t 空间数据进行傅里叶反变换并获得重建图像。 0007 优选地, 在对 k 空间进行采集, 并获取导航数据和测量数据的步骤之前还包括 : 心 电触发获取采集信号 ; 所述心电触发获取采集信号的步骤包括 : 采集心电图信号并确定心 动周期 ; 获取心脏动态成像, 且根据所。

15、述心脏动态成像确定在心动周期内的延迟时间及准 稳态期 ; 在所述准稳态期内施加弥散梯度脉冲且进行数据采集。 0008 优选地, 在对 k 空间进行采集, 并获取导航数据和测量数据的步骤之前还包括 : 呼 吸导航获取采集信号 ; 所述呼吸导航获取采集信号的步骤包括 : 发射二维脉冲信号 ; 根据 所述二维脉冲信号获取导航信号 ; 根据所述导航信号获取呼吸运动幅度值 ; 判断所述呼吸 运动幅度值是否超出预设的导航窗, 是, 则对所述呼吸运动幅度值所对应的时间区间的采 集信号进行摒弃, 否, 则所述呼吸运动幅度值所对应的时间区间的采集信号进行采集。 0009 优选地, 在对 k 空间进行采集, 并获。

16、取导航数据和测量数据的步骤之前还包括 : 心 说 明 书 CN 102008307 B 4 2/6 页 5 电触发获取采集信号及呼吸导航获取采集信号的步骤 ; 根据心电触发获取采集信号获取准 稳态时间, 根据所述导航获取采集信号获取在预设导航窗内的呼吸运动幅度值 ; 在所述准 稳态期内以及符合所述呼吸运动幅度值的时间区间进行信号采集。 0010 还有必要提供一种磁共振弥散张量成像系统, 包括 : 采集模块, 用于对 k 空间进行 采集, 并获取导航数据和测量数据 ; 从 K 空间采集的数据 S(k, t) 基于可分离函数技术, 表 示为 :所述 L 为频率成分参数、 cl(k) 为空间基函数及。

17、为时间基函 数 ; 提取模块, 用于根据所述导航数据获取时间基函数和频率成分参数, 以及根据测量数据 获取空间基函数 ; 所述导航数据确定时间基函数和频率成分参数, 该测量数据确定空间基 函数 ; 所采集的点需要满足 3 个条件 : 1.TR必须满足 Snaw(k, t) 的时间奈奎斯特速率, TR是 重复时间, Snaw(k, t) 是导航数据 ; 2.ky必须满足 Simg(k, t) 的时间奈奎斯特速率, ky是 相邻相位编码线的 k 空间间距 ; 3. 从 Simg(k, t) 获取的采样帧数 N 必须大于或等于阶数 L, Simg(k, t) 是测量数据, 帧数 N 表示不同时刻获得。

18、的图像的总数 ; 计算模块, 用于根据所述时 间基函数、 频率成分参数和空间基函数计算获得 k-t 空间数据 ; 重建模块, 根据所述 k-t 空 间数据进行傅里叶反变换并获得重建图像。 0011 优选地, 还包括 : 与采集模块连接的心电触发模块, 所述心电触发模块包括 : 心电 采集单元, 用于采集心电图信号并确定心动周期 ; 获取单元, 用于获取心脏动态成像, 且根 据所述心脏动态成像确定在心动周期内的延迟时间及准稳态期 ; 心电执行单元, 用于在准 稳态期内施加弥散梯度脉冲且进行数据采集。 0012 优选地, 还包括 : 与采集模块连接的呼吸导航模块, 所述呼吸导航模块包括 : 发射 。

19、单元, 用于发射二维脉冲信号 ; 接收单元, 用于根据所述二维脉冲信号获取导航信号 ; 呼吸 运动幅度值单元, 用于根据所述导航信号获取呼吸运动幅度值 ; 处理单元, 用于判断所述呼 吸运动幅度值是否超出预设的导航窗, 是, 则对所述呼吸运动幅度值所对应的时间区间的 采集信号进行摒弃, 否, 则所述呼吸运动幅度值所对应的时间区间的采集信号进行采集。 0013 采用上述磁共振弥散张量成像方法及系统, 基于部分可分离函数技术, 加快获取 扫描数据, 快速获取运动物体的成像信息, 达到快速成像。 【附图说明】 0014 图 1 为磁共振弥散张量成像方法的流程图 ; 0015 图 2 为一实施例的基于。

20、部分可分离函数技术的数据采集示意图 ; 0016 图 3 为心电触发获取采集信号的流程图 ; 0017 图 4 为一实施例的心电触发获取采集信号的心电触发成像原理图 ; 0018 图 5 为一实施例的心电触发获取采集信号的一个心动周期内心脏解剖结构 “收 缩 - 舒张 - 收缩” 示意图 ; 0019 图 6 为一实施例的心电触发获取采集信号的同一区域的心脏剖面信号连接图 ; 0020 图 7 为一实施例的呼吸导航获取采集信号流程图 ; 0021 图 8 为一实施例的呼吸导航获取采集信号中放置于隔肌上的发射二维脉冲信号 的示意图 ; 0022 图 9 为一实施例的同时采用心电触发获取采集信号及。

21、呼吸导航获取采集信号的 说 明 书 CN 102008307 B 5 3/6 页 6 方法流程图 ; 0023 图 10 为一实施例的同时采用心电触发及呼吸导航获取采集信号的示意图 ; 0024 图 11 为磁共振弥散张量成像装置的原理框图 ; 0025 图 12 为磁共振弥散张量成像装置中采用心电触发模块的原理框图 ; 0026 图 13 为磁共振弥散张量成像装置中采用呼吸导航模块的原理框图。 【具体实施方式】 0027 对于运动物体的磁共振成像过程中, 图像是在其空间位置 和时间 t 的函数。若磁 共振的扫描速度足够快, 在物体的空间位置还来不及发生很大变化的瞬间就可以采集重建 图像所需要。

22、的全部数据, 则可近似认为信号扫描的瞬间物体处于静止状态。 但是, 目前的磁 共振设备无法达到如此高的扫描速度。 0028 基于此, 很难在相对静止的瞬间采集足够的数据进行图像重建。此时所采集的磁 共振信号与自旋质子密度实际是空间位置 和时间 t 的函数, 其中接收信号 S(k, t) 和所需 要的图像函数 (r, t) 之间的关系为 : 0029 0030 请参阅附图 1 和 2, 一种磁共振弥散张量成像方法, 包括以下步骤 : 0031 S10 : 对 k 空间进行采集, 并获取导航数据和测量数据。 0032 具体地, 经过一段时间对 k 空间的信号进行采集, 此时空间中某一平面数据将包 。

23、含来自运动物体各个不同部位的信息。 0033 对该采集的数据基于部分可分离函数技术, 即认为图像函数 (r, t) 的空间变化 和时间变化是L阶可分离的(阶数与成像物体的运动模式相关, 理论上阶数越大, 对运动的 描述越精确 ), 利用部分可分离函数的性质, S(k, t) 可以表示为空间和时间两个独立变量 函数 cl(k) 和 0034 0035 所述 L 为频率成分参数、 cl(k) 为空间基函数及为时间基函数。 0036 S20 : 根据导航数据获取时间基函数和频率成分参数, 以及根据测量数据获取空间 基函数。 0037 具体地, 为了获得频率成分参数 L、 空间基函数 cl(k) 及时。

24、间基函数只需 要采集两组 k 空间数据即可。即对 k 空间进行采集, 获取高时间、 低空间分辨率导航数 据 (Navigator data, 空心圆圈 Snav(k, t), 见图 2), 获取高时间、 低空间分辨率测量数据 (Measurement data, 实心圆点 Simg(k, t), 见图 2)。该导航数据确定和 L, 该测量数据确 定 cl(k)。 0038 而且, 所采集的点需要满足 3 个条件 : 0039 1.TR必须满足 Snav(k, t) 的时间奈奎斯特速率, TR是重复时间, Snav(k, t) 是导航数 据 ; 0040 2.ky必须满足 Simg(k, t) 。

25、的时间奈奎斯特速率, ky是相邻相位编码线的 k 空间 间距 ; 说 明 书 CN 102008307 B 6 4/6 页 7 0041 3. 从 Simg(k, t) 获取的采样帧数 N 必须大于或等于阶数 L, Simg(k, t) 是测量数据, 帧数 N 表示不同时刻获得的图像的总数。 0042 S30 : 根据时间基函数、 频率成分参数和空间基函数计算获得 k-t 空间数据。根 据获得的时间基函数、 频率成分参数和空间基函数后, 则可以计算出扩大的 k-t 空间数据 (Synthetic data, 实心叉, 见图 2)。 0043 S40 : 根据所述 k-t 空间数据进行傅里叶反变。

26、换并获得重建图像。 0044 基于部分可分离函数技术, 加快获取扫描数据, 快速获取运动物体的成像信息, 达 到快速成像。 0045 在一实施例中, 结合附图 3 6, 在步骤 S10 之前还包括心电触发获取采集信号的 步骤A。 在弥散梯度脉冲施加期间, 物体中的水分子在弥散梯度方向上的布朗运动将积累净 相位, 在成像时与静止的分子相比信号强度减弱。 可见, 运动的物体在相同的位置上施加弥 散梯度脉冲, 才能够保证静止的分子净相位累计为 0。 0046 具体地, 心电触发获取采集信号的步骤 A 包括 : 0047 A1 : 采集心电图信号并确定心动周期。具体地, 采集心电图 (ECG)R 波为。

27、信号测量 的参考点, 并确定心脏的心动周期。 0048 A2 : 获取心脏动态成像, 且根据所述心脏动态成像确定在心动周期内的延迟时间 及准稳态期。具体地, 通过施加一种磁共振脉冲序列, 例如平衡式稳态自由进动梯度回波 (Balance-FFE)序列, 根据该序列快速获得心脏动态成像。 然后根据心脏动态成像确定在心 动周期内的延迟时间及准稳态期。 0049 对于延迟时间 (Delay time)D 选取是否能够达到预计的效果, 即保证运动物体的 相对位移值为零, 现选取一具体实施例进行详细说明。 0050 按上述方法, 首先通过施加一种磁共振脉冲序列获得心脏动态成像。 然后, 在一个 心动周期。

28、内捕获 20 个或以上心脏相位对应的心脏解剖结构, 包括从 “舒张态 - 收缩态 - 舒 张态” 的周而复始的变化, 根据心内膜和心外膜的轮廓的位移判断心脏运动情况 ( 见附图 5)。结合附图 6, 首先, 放置一条跨越左心室中心点的采值线 ( 涵盖心肌和左心室腔 ), 则可 提取该线所覆盖区域的信号强度, 从而获得相应的信号轮廓。然后, 将 20 个心脏相位上同 一区域的信号连接起来, 可以得出在一个心动周期内心肌及左心室腔的变化情况。不难看 出, 在白色线段标注范围内, 心肌和左心室腔基本没有明显的变化, 表明在该时间段内心脏 运动的变化不大。可见在该时间段, 心脏跳动是比较舒缓的, 因此。

29、可以据此确定延迟时间 (Delay time)D。在延迟时间 D 后, 且在一个心动周期内, 确定为准稳态期。 0051 A3 : 在准稳态期内施加弥散梯度脉冲且进行数据采集。 具体地, 在该准稳态期内发 射弥散张量成像 (DTI) 的弥散梯度脉冲, 并进行数据采集。以至于所获取的采集数据更加 精确, 更少受到心脏运动的影响。 0052 采用该方案, 通过对心电图信号以及心脏动态成像确定心动周期及准稳态期, 进 而在准稳态期内施加弥散梯度脉冲且进行数据采集, 准确获得成像数据, 减少运动伪影。 0053 在另一实施例中, 结合附图 7 8, 在步骤 S10 之前还包括呼吸导航获取采集信号 的步。

30、骤 B。采用呼吸导航技术, 主要是用于检测自由呼吸下隔面位置的变化, 根据隔面的位 置信息确定采集信号。 0054 具体地, 呼吸导航获取采集信号的步骤 B 包括 : 说 明 书 CN 102008307 B 7 5/6 页 8 0055 B1 : 发射二维脉冲信号。具体地, 在隔肌之上 ( 见图 8) 发射二维脉冲信号。 0056 B2 : 根据二维脉冲信号获取导航信号。 具体地, 根据二维脉冲信号获取导航笔形波 束的导航信号。 0057 B3 : 根据导航信号获取呼吸运动幅度值。 具体地, 在整个扫描过程中回波信号代表 了该时间段内隔肌的运动情况, 膈肌的运动同时反映了呼吸运动, 故可以获。

31、得呼吸运动幅 度值。 0058 B4 : 判断呼吸运动幅度值是否超出预设的导航窗, 是, 则对呼吸运动幅度值所对应 的时间区间的采集信号进行摒弃, 否, 则所述呼吸运动幅度值所对应的时间区间的采集信 号进行采集。具体地, 预设一导航窗, 当呼吸运动的幅值在导航窗的范围内, 该范围内所对 应的时间区间则认为此数据是有效且采集该数据。当呼吸运动的幅值在导航窗的范围外, 该范围内所对应的时间区间则认为该呼吸运动幅度较大, 应该摒弃在该时间所采集的数 据。 0059 采用该实施例的方案, 所采集的数据为运动物体 ( 例如在体心脏 ) 较为平稳缓和 处, 即该呼吸运动较为平缓处, 则所获取的数据准确, 。

32、在很大程度上去除运动伪影的干扰。 0060 在其它实施例, 请参阅附图910, 同时采用步骤A和步骤B的方法, 即在步骤S10 之前还包括步骤 C, 包括 : 0061 C1 : 根据心电触发获取采集信号获取准稳态时间, 根据所述呼吸导航获取采集信 号获取在预设导航窗内的呼吸运动幅度值。 0062 C2 : 在准稳态期内以及符合所述呼吸运动幅度值所对应的时间区间施加弥散梯度 脉冲且进行数据采集。 0063 结合附图 10, 采用该实施例方案, 获取了运动较为平缓处以及呼吸幅度值较为平 稳处的信号, 则所采集的数据更为准确, 在很大程度上去除运动伪影的干扰。 0064 请参阅附图 11, 基于上。

33、述的磁共振弥散张量成像方法还提供一种磁共振弥散张量 成像系统。 0065 一种磁共振弥散张量成像系统, 包括 : 0066 采集模块 10, 用于对 k 空间进行采集, 并获取导航数据和测量数据。 0067 提取模块 20, 用于根据所述导航数据获取时间基函数和频率成分参数, 以及根据 测量数据获取空间基函数。 0068 计算模块 30, 用于根据所述时间基函数、 频率成分参数和空间基函数计算获得 k-t 空间数据。 0069 重建模块 40, 根据所述 k-t 空间数据进行傅里叶反变换并获得重建图像。 0070 请参阅附图 12, 在一实施例中, 该磁共振弥散张量成像系统还包括 : 与采集模。

34、块 连接的心电触发模块 50, 该心电触发模块 50 包括 : 0071 心电采集单元 51, 用于采集心电图信号并确定心动周期。 0072 获取单元 52, 用于获取心脏动态成像, 且根据所述心脏动态成像确定在心动周期 内的延迟时间及准稳态期。 0073 心电执行单元 53, 用于在准稳态期内施加弥散梯度脉冲且进行数据采集。 0074 请参阅附图 13, 在一实施例中, 该磁共振弥散张量成像系统还包括 : 与采集模块 连接的呼吸导航模块 60, 该呼吸导航模块 60 包括 : 说 明 书 CN 102008307 B 8 6/6 页 9 0075 发射单元 61, 用于发射二维脉冲信号。 0。

35、076 接收单元 62, 用于根据所述二维脉冲信号获取导航信号。 0077 呼吸运动幅度值单元 63, 用于根据所述导航信号获取呼吸运动幅度值。 0078 处理单元 64, 用于判断所述呼吸运动幅度值是否超出预设的导航窗, 是, 则对所述 呼吸运动幅度值所对应的时间区间的采集信号进行摒弃, 否, 则所述呼吸运动幅度值所对 应的时间区间的采集信号进行采集。 0079 基于新型的磁共振弥散张量成像系统, 加快获取扫描数据, 快速获取运动物体的 成像信息。 0080 以上所述实施例仅表达了本发明的几种实施方式, 其描述较为具体和详细, 但并 不能因此而理解为对本发明专利范围的限制。应当指出的是, 对。

36、于本领域的普通技术人员 来说, 在不脱离本发明构思的前提下, 还可以做出若干变形和改进, 这些都属于本发明的保 护范围。因此, 本发明专利的保护范围应以所附权利要求为准。 说 明 书 CN 102008307 B 9 1/7 页 10 图 1 说 明 书 附 图 CN 102008307 B 10 2/7 页 11 图 2 图 3 说 明 书 附 图 CN 102008307 B 11 3/7 页 12 图 4 图 5 图 6 说 明 书 附 图 CN 102008307 B 12 4/7 页 13 图 7 图 8 说 明 书 附 图 CN 102008307 B 13 5/7 页 14 图 9 图 10 说 明 书 附 图 CN 102008307 B 14 6/7 页 15 图 11 图 12 说 明 书 附 图 CN 102008307 B 15 7/7 页 16 图 13 说 明 书 附 图 CN 102008307 B 16 。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 人类生活必需 > 医学或兽医学;卫生学


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1