CN112581385A - 基于多先验约束扩散峰度成像张量估计方法、介质和设备 - Google Patents
基于多先验约束扩散峰度成像张量估计方法、介质和设备 Download PDFInfo
- Publication number
- CN112581385A CN112581385A CN202011371326.7A CN202011371326A CN112581385A CN 112581385 A CN112581385 A CN 112581385A CN 202011371326 A CN202011371326 A CN 202011371326A CN 112581385 A CN112581385 A CN 112581385A
- Authority
- CN
- China
- Prior art keywords
- tensor
- dki
- diffusion
- model
- voxel
- 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.)
- Granted
Links
- 238000009792 diffusion process Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000003384 imaging method Methods 0.000 title claims abstract description 27
- 238000012937 correction Methods 0.000 claims abstract description 13
- 238000011524 similarity measure Methods 0.000 claims abstract description 8
- 238000013139 quantization Methods 0.000 claims abstract description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 230000035945 sensitivity Effects 0.000 claims description 6
- 238000003745 diagnosis Methods 0.000 abstract description 2
- 201000010099 disease Diseases 0.000 abstract description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 6
- 230000006872 improvement Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000002598 diffusion tensor imaging Methods 0.000 description 5
- 210000004556 brain Anatomy 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Databases & Information Systems (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Pathology (AREA)
- Data Mining & Analysis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
Abstract
基于多先验约束扩散峰度成像张量估计方法、介质和设备,包括步骤:采用有约束的加权线性最小二乘拟合方法得到DKI参数的初始估计值;根据图像背景计算高斯噪声的标准差;根据初始估计值构建参数集、或根据初始估计值和高斯噪声的标准差构建参数集;根据基于体素的非局部结构相似性测度模型,计算每个体素搜索窗中的权重;建立DKI扩散信号在非中心卡方分布下的一阶矩噪声校正模型;建立DKI张量的物理约束模型;由权重、物理约束模型、一阶矩噪声校正模型、基于局部全变分的约束项构建DKI张量估计目标函数模型;以参数集作为输入,计算DKI张量场,并根据DKI张量场计算得到DKI量化参数图。本发明准确得到DKI量化参数图,为临床疾病的量化诊断提供可靠信息。
Description
技术领域
本发明属于扩散磁共振参数成像技术领域,具体地说属于一种可以准确估计扩散峰度成像(diffusion kurtosis imaging,DKI)张量场的方法。
背景技术
扩散磁共振成像(diffusion magnetic resonance imaging,dMRI)技术可通过检测活体组织内水分子的微观扩散运动状态来反映机体组织结构的生理、病理等特点。扩散张量成像(diffusion tensor imaging,DTI)是一种较为传统的dMRI技术,通过在不同扩散编码方向上采集到的图像进行张量成像,从而获得临床相关的量化指标。DTI理论基础是水分子成高斯扩散,但是组织内水分子受微环境影响,呈现非高斯扩散特性。扩散峰度成像(diffusionkurtosisimaging,DKI)是基于DTI技术的扩展,采用非高斯扩散模型,更符合组织内水分子的实际扩散特性,从而更适合描述组织微观结构的变化。与传统DTI技术相比,DKI需要采集更高b值(2000~3000s/mm2)的扩散加权(diffusion-weighted,DW)图像。其中b值为扩散敏感系数,表示扩散加权程度,b值越高,信号衰减越快,因此高b值的DW图像受噪声影响更为严重。严重的噪声会导致DKI张量估计不准确,而具有临床诊断价值的DKI参数均是由DKI张量直接计算得到的,因此不准确的DKI张量直接影响了后续DKI参数的准确性,从而限制了DKI技术在临床上的推广。
目前临床常采用多次采集平均技术来提高DW图像的信噪比,但该方法增加了采集时间,提高了采集成本。为了降低噪声的影响,目前主要有两种后处理技术,一种是先对图像进行去噪,然后根据降噪图像来估计DKI张量。然而任何一种去噪算法都不能完美恢复图像,图像去噪过程中或多或少会引入一些误差,该误差会传播到后续的张量估计当中,导致误差累积;另一种方法是在张量估计的同时加入基于先验信息的正则化项,从而直接从噪声图像得到平滑的张量场。现有的基于先验约束的张量估计方法并未考虑磁共振图像噪声特性,会导致估计得到的张量产生偏差,尤其是对于DKI模型中的峰度张量。此外,现有的基于先验约束的张量估计方法并没有同时考虑多种约束,并且主要是用来解决DTI成像技术。由于DKI中高b值的DW图像受噪声污染特别严重,此类方法直接应用到DKI上,并不能产生理想的参数估计效果,尤其是峰度张量。
发明内容
为了克服上述技术缺陷,本发明提供了一种基于多先验约束扩散峰度成像张量估计方法、介质和设备,其能直接从原始噪声DW图像拟合得到准确、无偏的张量场,进而得到准确的DKI量化参数图。
为了解决上述问题,本发明按以下技术方案予以实现的:
一种基于多先验约束的扩散峰度成像张量估计方法,包括步骤:
采用有约束的加权线性最小二乘拟合方法得到DKI参数的初始估计值;
根据图像背景计算高斯噪声的标准差;
根据所述初始估计值构建参数集、或根据所述初始估计值和所述高斯噪声的标准差构建参数集;
根据基于体素的非局部结构相似性测度模型,计算每个体素搜索窗中的权重;
建立DKI扩散信号在非中心卡方分布下的一阶矩噪声校正模型;
建立DKI张量的物理约束模型;
由所述权重、所述物理约束模型、所述一阶矩噪声校正模型、基于局部全变分的约束项构建DKI张量估计目标函数模型;
以所述参数集作为输入,计算得到DKI张量场,并根据所述DKI张量场计算得到DKI量化参数图。
作为本发明的进一步改进,所述根据基于体素的非局部结构相似性测度模型,计算每个体素搜索窗中的权重的步骤,包括步骤:
利用扩散加权图像在扩散编码方向上的相关性,将同一位置不同扩散编码方向的像素视为一个体素,采用基于体素的非局部结构相似性测度模型计算得到所述权重。
作为本发明的进一步改进,所述建立DKI扩散信号在非中心卡方分布下的一阶矩噪声校正模型的步骤,包括步骤:
根据采集数据的线圈通道个数、DKI信号模型、所述高斯噪声标准差构建所述一阶矩噪声校正模型。
作为本发明的进一步改进,所述DKI信号模型由体素中的非扩散加权信号、梯度编码方向的分量、扩散敏感系数、体素对应的扩散张量和峰度张量的各个元素、体素的平均扩散系数构建而成。
作为本发明的进一步改进,所述建立DKI张量的物理约束模型的步骤,包括步骤:
根据梯度编码方向的分量、扩散敏感系数、体素对应的扩散张量和峰度张量的各个元素、体素的平均扩散系数构建所述DKI张量的物理约束模型。
作为本发明的进一步改进,所述计算得到DKI张量场的步骤,具体为:
采用预条件优化算法计算得到所述DKI张量场。
作为本发明的进一步改进,所述建立的DKI张量估计模型还包括步骤:
采用基于局部的全变分方法对所述参数集中的每个参数进行空间平滑约束。
与现有技术相比,本发明具有以下有益效果:本发明综合考虑了DKI的以下特性:①、DW图像噪声服从非中心卡方分布,而非高斯分布;②、DW图像在空间上具有冗余性,在扩散编码方向上具有相关性;③、参数图在空间上是连续缓慢变化的;④、估计得到的DKI张量需要保证其DKI模型的物理意义,基于此,本发明针对以上几种已知的先验特性,引入了一阶矩噪声校正模型、基于结构相似性测度以及局部全变分的平滑约束项,以及保证DKI模型物理意义的约束项,并将以上先验信息纳入到DKI张量估计目标函数模型中,从而直接从原始噪声DW图像拟合得到准确、无偏的张量场,进而到准确的DKI量化参数图,为临床疾病的量化诊断提供可靠信息。
此外,本发明还提供了一种计算机可读存储介质,所述可读存储介质中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由处理器加载并执行以实现上述的基于多先验约束的扩散峰度成像张量估计方法。
本发明还提供了一种计算机设备,所述计算机设备包括处理器和存储器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由所述处理器加载并执行以实现上述的基于多先验约束的扩散峰度成像张量估计方法。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明,其中:
图1基于多先验约束的扩散峰度成像张量估计方法示意图;
图2是实施例1脑部DKI仿真数据的实验方法对比结果,仿真数据的噪声是空间不变的;(a)扩散张量的相关参数图以及误差图;(b)峰度张量的相关参数图以及误差图。右下角数字是对应参数图的RMSE;
图3是实施例2脑部DKI仿真数据的实验方法对比结果,仿真数据的噪声是空间变化的;(a)扩散张量的相关参数图以及误差图;(b)峰度张量的相关参数图以及误差图。右下角数字是对应参数图的RMSE;
图4是实施例3脑部DKI真实数据的实验方法对比结果。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
实施例1
本实施例公开了一种基于多先验约束的扩散峰度成像张量估计方法,首先,采用HCP数据(https://www.humanconnectome.org)进行数据仿真,仿真数据包含1个b=0s/mm2,45个b=1000s/mm2和45个b=2000s/mm2,仿真的噪声是空间不变的,σg=0.05,通道个数L=8,
如图1和图2所示,该基于多先验约束的扩散峰度成像张量估计方法包括步骤:
S2、根据图像背景计算高斯噪声的标准差σg。
S3、根据所述初始估计值构建参数集,将噪声水平σg当作已知参数,不做优化处理,即参数集Θ={S0,ΘD,ΘK}。
S4、利用扩散加权图像在扩散编码方向上的相关性,将同一位置不同扩散编码方向的像素视为一个体素,根据基于体素的非局部结构相似性测度模型,计算每个体素搜索窗中的权重;其中,基于体素的非局部结构相似性测度模型为:
式(1)中,P(yi)和P(yj)表示以体素i和j为中心的DKI结构块,大小分别为M×M×N。M×M为邻域块大小,N为所有扩散加权图像和非扩散加权图像的总数。Vi为以体素i为中心的搜索窗,h大小决定了参数图的平滑度,本实施例中,参数h=0.2σg,搜索窗大小为11×11,块大小为5×5×91,91为所有扩散编码方向和非扩散编码方向的总数;
S5、根据采集数据的线圈通道个数L、DKI信号模型Sdki(Θi)、高斯噪声的标准差σg建立DKI扩散信号在非中心卡方分布下的一阶矩噪声校正模型f(Θi);
式(2)中,1F1为合流超几何函数,Θi∈RL为参数集Θ∈RM×N×L在体素i位置的所有参数,其中L为每个体素所对应的参数个数。对于空间不变的噪声,σg作为已知量,L=22,即每个体素包含6个描述扩散张量的参数,15个描述扩散峰度张量的参数,以及一个描述非扩散加权信号的参数。DKI信号模型Sdki(Θi)由体素i中的非扩散加权信号S0i、梯度编码方向的三个分量g1,g2,g3、扩散敏感系数b、体素i对应的扩散张量Di∈R3×3和峰度张量Wi∈R3 ×3×3×3的各个元素体素的平均扩散系数MDi构建而成,即:
S6、建立DKI张量的物理约束模型R(ΘD,ΘK);具体的,物理约束模型根据梯度编码方向的三个分量g1,g2,g3、扩散敏感系数b、体素i对应的扩散张量Di∈R3×3和峰度张量Wi∈R3 ×3×3×3的各个元素体素i的平均扩散系数MDi构建而成,即:
式(4)中,c为很大的一个常数,从而保证当所有体素的扩散张量和峰度张量均满足DKI模型的物理意义时,R(ΘD,ΘK)→0,否则,R(ΘD,ΘK)→∞。
S7、由所述权重wij、所述物理约束模型R(ΘD,ΘK)、所述一阶矩噪声校正模型f(Θi)、基于非局部全变分的约束项TV(·)构建DKI张量估计目标函数模型,公式如下:
式(5)中,yj为一个向量,表示体素j点的所有扩散加权和非扩散加权信号的灰度值,ΘD∈RM×N×6表示扩散张量场D对应的参数集合,ΘK∈RM×N×15表示峰度张量场K对应的参数集合。Θq∈RM×N表示第q个参数的参数图,α和β表示各自约束项的权重,表示梯度算子,||·||1表示1-范数。
S8、以参数集Θ作为输入,采用适合大规模参数计算的l-BFGS优化算法,计算得到DKI张量场,并根据所述DKI张量场计算得到DKI量化参数图,其包括:扩散张量场(diffusion tensor field,DT field)和峰度张量场(kurtosis tensor field,KTfield),并根据DKI张量场计算得到DKI量化参数图,包括由DT field得到的MD,AD,RD和FA参数图,以及由KT field得到的MK,AK和RK参数图。
实施例2
如图3所示,本实施例公开了另一种基于多先验约束的扩散峰度成像张量估计方法,其与实施例1的区别在于:本实施例考虑了空间变化的噪声,并将高斯噪声的标准差σg作为未知参数进行优化,此时,Θi∈RL为参数集Θ∈RM×N×L在体素i位置的所有参数,L=23,由于σg与其他参数{S0,ΘD,ΘK}的最优步长的尺度不同,因此采用预条件l-BFGS优化算法计算得到DKI张量场,并根据DKI张量场计算得到DKI量化参数图。预条件优化算法具体为:当把σg也当做未知参数进行优化时,将关于σg的偏微分导数的尺度设为0.01。
实施例3
如图4所示,本实施例公开了另一种基于多先验约束的扩散峰度成像张量估计方法,其与实施例1的区别在于:本实施例采用了真实数据,该数据来自健康志愿者的大脑扩散加权数据,包括一个b=0s/mm2;30个b=1000s/mm2和30个b=2000s/mm2,采集序列为double SE-EPI,部分傅立叶采样(在相位编码方向采集5/8条K空间线)。
实施例4
本实施例公开了一种计算机可读存储介质,该可读存储介质中存储有至少一条指令、至少一段程序、代码集或指令集,至少一条指令、至少一段程序、代码集或指令集由处理器加载并执行以实现实施例1-实施例3的基于多先验约束的扩散峰度成像张量估计方法。
实施例5
一种计算机设备,其特征在于,该计算机设备包括处理器和存储器,存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由处理器加载并执行以实现实施例1-实施例3的基于多先验约束的扩散峰度成像张量估计方法。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,故凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (9)
1.基于多先验约束的扩散峰度成像张量估计方法,其特征在于,包括步骤:
采用有约束的加权线性最小二乘拟合方法得到DKI参数的初始估计值;
根据图像背景计算高斯噪声的标准差;
根据所述初始估计值构建参数集、或根据所述初始估计值和所述高斯噪声的标准差构建参数集;
根据基于体素的非局部结构相似性测度模型,计算每个体素搜索窗中的权重;
建立DKI扩散信号在非中心卡方分布下的一阶矩噪声校正模型;
建立DKI张量的物理约束模型;
由所述权重、所述物理约束模型、所述一阶矩噪声校正模型、基于局部全变分的约束项构建DKI张量估计目标函数模型;
以所述参数集作为输入,计算得到DKI张量场,并根据所述DKI张量场计算得到DKI量化参数图。
2.根据权利要求1所述的基于多先验约束的扩散峰度成像张量估计方法,其特征在于,所述根据基于体素的非局部结构相似性测度模型,计算每个体素搜索窗中的权重的步骤,包括步骤:
利用扩散加权图像在扩散编码方向上的相关性,将同一位置不同扩散编码方向的像素视为一个体素,采用基于体素的非局部结构相似性测度模型计算得到所述权重。
3.根据权利要求1所述的基于多先验约束的扩散峰度成像张量估计方法,其特征在于,所述建立DKI扩散信号在非中心卡方分布下的一阶矩噪声校正模型的步骤,包括步骤:
根据采集数据的线圈通道个数、DKI信号模型、所述高斯噪声标准差构建所述一阶矩噪声校正模型。
4.根据权利要求3所述的基于多先验约束的扩散峰度成像张量估计方法,其特征在于,所述DKI信号模型由体素中的非扩散加权信号、梯度编码方向的分量、扩散敏感系数、体素对应的扩散张量和峰度张量的各个元素、体素的平均扩散系数构建而成。
5.根据权利要求1所述的基于多先验约束的扩散峰度成像张量估计方法,其特征在于,所述建立DKI张量的物理约束模型的步骤,包括步骤:
根据梯度编码方向的分量、扩散敏感系数、体素对应的扩散张量和峰度张量的各个元素、体素的平均扩散系数构建所述DKI张量的物理约束模型。
6.根据权利要求1所述的基于多先验约束的扩散峰度成像张量估计方法,其特征在于,所述计算得到DKI张量场的步骤,具体为:
采用预条件优化算法计算得到所述DKI张量场。
7.根据权利要求1所述的基于多先验约束的扩散峰度成像张量估计方法,其特征在于,所述建立的DKI张量估计模型还包括步骤:
采用基于局部的全变分方法对所述参数集中的每个参数进行空间平滑约束。
8.一种计算机可读存储介质,其特征在于,所述可读存储介质中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由处理器加载并执行以实现如权利要求1至7任一所述的基于多先验约束的扩散峰度成像张量估计方法。
9.一种计算机设备,其特征在于,所述计算机设备包括处理器和存储器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由所述处理器加载并执行以实现如权利要求1至7任一所述的基于多先验约束的扩散峰度成像张量估计方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011371326.7A CN112581385B (zh) | 2020-11-30 | 2020-11-30 | 基于多先验约束扩散峰度成像张量估计方法、介质和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011371326.7A CN112581385B (zh) | 2020-11-30 | 2020-11-30 | 基于多先验约束扩散峰度成像张量估计方法、介质和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112581385A true CN112581385A (zh) | 2021-03-30 |
CN112581385B CN112581385B (zh) | 2023-07-07 |
Family
ID=75126374
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011371326.7A Active CN112581385B (zh) | 2020-11-30 | 2020-11-30 | 基于多先验约束扩散峰度成像张量估计方法、介质和设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112581385B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113270168A (zh) * | 2021-05-19 | 2021-08-17 | 中科芯未来微电子科技成都有限公司 | 一种提高医学图像处理能力的方法及系统 |
CN117310581A (zh) * | 2023-10-11 | 2023-12-29 | 无锡鸣石峻致医疗科技有限公司 | 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104574298A (zh) * | 2014-12-25 | 2015-04-29 | 天津大学 | 一种基于互信息的多b值扩散权重图像的降噪方法 |
US20150362575A1 (en) * | 2013-02-01 | 2015-12-17 | Ucl Business Plc | Apparatus and method for correcting susceptibility artefacts in a magnetic resonance image |
CN106997034A (zh) * | 2017-04-25 | 2017-08-01 | 清华大学 | 基于以高斯模型为实例整合重建的磁共振扩散成像方法 |
WO2018001759A1 (en) * | 2016-06-28 | 2018-01-04 | Koninklijke Philips N.V. | Diffusion weighted mr imaging using multi-shot epi with motion detection and modified sense reconstruction |
CN108171690A (zh) * | 2017-12-22 | 2018-06-15 | 渤海大学 | 基于结构先验信息的人体心脏左心室扩散张量估计方法 |
-
2020
- 2020-11-30 CN CN202011371326.7A patent/CN112581385B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150362575A1 (en) * | 2013-02-01 | 2015-12-17 | Ucl Business Plc | Apparatus and method for correcting susceptibility artefacts in a magnetic resonance image |
CN104574298A (zh) * | 2014-12-25 | 2015-04-29 | 天津大学 | 一种基于互信息的多b值扩散权重图像的降噪方法 |
WO2018001759A1 (en) * | 2016-06-28 | 2018-01-04 | Koninklijke Philips N.V. | Diffusion weighted mr imaging using multi-shot epi with motion detection and modified sense reconstruction |
CN106997034A (zh) * | 2017-04-25 | 2017-08-01 | 清华大学 | 基于以高斯模型为实例整合重建的磁共振扩散成像方法 |
CN108171690A (zh) * | 2017-12-22 | 2018-06-15 | 渤海大学 | 基于结构先验信息的人体心脏左心室扩散张量估计方法 |
Non-Patent Citations (1)
Title |
---|
吴锡;何晋;王玉;谢明元;周激流;: "改进Wiener滤波弥散加权磁共振图像Rician噪声复原", 中国生物医学工程学报, no. 02, pages 10 - 15 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113270168A (zh) * | 2021-05-19 | 2021-08-17 | 中科芯未来微电子科技成都有限公司 | 一种提高医学图像处理能力的方法及系统 |
CN113270168B (zh) * | 2021-05-19 | 2024-01-02 | 中科芯未来微电子科技成都有限公司 | 一种提高医学图像处理能力的方法及系统 |
CN117310581A (zh) * | 2023-10-11 | 2023-12-29 | 无锡鸣石峻致医疗科技有限公司 | 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 |
CN117310581B (zh) * | 2023-10-11 | 2024-05-10 | 安徽峻德医疗科技有限公司 | 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112581385B (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Barbieri et al. | Deep learning how to fit an intravoxel incoherent motion model to diffusion‐weighted MRI | |
Bollmann et al. | DeepQSM-using deep learning to solve the dipole inversion for quantitative susceptibility mapping | |
Friman et al. | A Bayesian approach for stochastic white matter tractography | |
US9430854B2 (en) | System and method for model consistency constrained medical image reconstruction | |
CN107240125B (zh) | 一种扩散加权成像方法 | |
US11346912B2 (en) | Systems and methods of generating robust phase images in magnetic resonance images | |
CN110969614B (zh) | 基于三维卷积神经网络的脑龄预测方法及系统 | |
Tabelow et al. | Local estimation of the noise level in MRI using structural adaptation | |
CN112991483B (zh) | 一种非局部低秩约束的自校准并行磁共振成像重构方法 | |
Lorenzi et al. | 4D registration of serial brain’s MR images: a robust measure of changes applied to Alzheimer’s disease | |
CN114450599A (zh) | 麦克斯韦并行成像 | |
CN112581385B (zh) | 基于多先验约束扩散峰度成像张量估计方法、介质和设备 | |
CN111445546A (zh) | 图像重建方法、装置、电子设备、存储介质 | |
van der Heide et al. | High‐resolution in vivo MR‐STAT using a matrix‐free and parallelized reconstruction algorithm | |
Zhou et al. | Evaluation of non-local means based denoising filters for diffusion kurtosis imaging using a new phantom | |
Mouches et al. | Unifying brain age prediction and age-conditioned template generation with a deterministic autoencoder | |
CN109471053B (zh) | 一种基于双约束的介电特性迭代成像方法 | |
Aja-Fernández et al. | Validation of deep learning techniques for quality augmentation in diffusion MRI for clinical studies | |
Aja-Fernández et al. | Scalar diffusion-MRI measures invariant to acquisition parameters: a first step towards imaging biomarkers | |
Kleineisel et al. | Real‐time cardiac MRI using an undersampled spiral k‐space trajectory and a reconstruction based on a variational network | |
Lee et al. | Registration and quantification network (RQnet) for IVIM‐DKI analysis in MRI | |
WO2016090686A1 (zh) | 磁共振弥散张量去噪方法和装置以及计算机程序产品 | |
Christiaens et al. | Modeling fiber orientations using diffusion MRI | |
Fadnavis et al. | Nuq: A noise metric for diffusion mri via uncertainty discrepancy quantification | |
Thai et al. | Using Deep Convolutional Neural Network for Mouse Brain Segmentation in DT-MRI |
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 |