【新书推荐】【2020】基于MATLAB®的并行MRI正则化图像重建

在这里插入图片描述

在加速并行磁共振成像(pMRI)中,正则化成为重建过程中不可或缺的一部分,因为需要以简约模型的形式利用最具鉴别能力的信息来生成高质量图像,同时减少噪声和伪影。

Regularization becomes an integral part of the reconstruction process in accelerated parallel magnetic resonance imaging (pMRI) due to the need for utilizing the most discriminative information in the form of parsimonious models to generate high quality images with reduced noise and artifacts.

除了提供各种pMRI重建方法的详细概述和实现细节之外,基于MATLAB的并行MRI正则化图像重建将pMRI正则化图像重建解释为一种有效控制两种特定类型误差信号之间平衡的方法,既可以提高缺失样本估计的精度,也可以加快估计过程。

Apart from providing a detailed overview and implementation details of various pMRI reconstruction methods, Regularized image reconstruction in parallel MRI with MATLAB examples interprets regularized image reconstruction in pMRI as a means to effectively control the balance between two specific types of error signals to either improve the accuracy in estimation of missing samples, or speed up the estimation process.

第一种类型对应于所获得的和其估计值之间的建模误差。

The first type corresponds to the modeling error between acquired and their estimated values.

第二类是由于自标定方法中k空间值的摄动或基于压缩感知的重建模型中的稀疏近似引起的。

The second type arises due to the perturbation of k-space values in autocalibration methods or sparse approximation in the compressed sensing based reconstruction model.

本书特点:

提供了在每种重建类型中优化正则化参数的详细信息。

比较了各种pMRI重建的正则化方法。

包括使用临床获得的数据进行案例研究讨论。

为每种重建类型提供了MATLAB代码。

包含自适应正则化以优化速度和精度的方法描述。

Features:

Provides details for optimizing regularization parameters in each type of reconstruction.

Presents comparison of regularization approaches for each type of pMRI reconstruction.

Includes discussion of case studies using clinically acquired data.

MATLAB codes are provided for each reconstruction type.

Contains method-wise description of adapting regularization to optimize speed and accuracy.

This book serves as a reference material for researchers and students involved in development of pMRI reconstruction methods. Industry practitioners concerned with how to apply regularization in pMRI reconstruction will find this book most useful.

更多精彩文章请关注公众号:在这里插入图片描述

TV正则化是一种用于图像重建的技术,它能够有效地处理图像中的噪声和模糊,产生更清晰的图像。下面是一个使用MATLAB实现TV正则化图像重建的例子: ```matlab % 读入原始图像 im = imread('lena.png'); im = im2double(im); % 添加高斯噪声 sigma = 0.1; noise = sigma * randn(size(im)); im_noisy = im + noise; % 设置TV正则化参数 lambda = 0.1; max_iter = 500; tol = 1e-6; % 进行TV正则化图像重建 im_recon = TV(im_noisy, lambda, max_iter, tol); % 显示结果 figure; subplot(1,3,1); imshow(im); title('原始图像'); subplot(1,3,2); imshow(im_noisy); title('带噪声图像'); subplot(1,3,3); imshow(im_recon); title('重建图像'); % TV正则化函数 function [im_recon] = TV(im_noisy, lambda, max_iter, tol) [h, w] = size(im_noisy); dx = [-1 1; 0 0]; dy = [-1 0; 1 0]; D = [kron(speye(h), sparse(dx)); kron(sparse(dy), speye(w))]; Dt = D'; x = im_noisy(:); y = x; L = max(eigs(Dt*D)); tau = 1/(1.01*L); sigma = 1/(1.01*L); for i = 1:max_iter x_old = x; y = y + sigma*D*x; x = soft_threshold(y, lambda*sigma); x = x./max(1, abs(x)); if norm(x-x_old)/norm(x) < tol break; end end im_recon = reshape(x, [h,w]); end % 软阈值函数 function [x] = soft_threshold(y, lambda) x = sign(y).*max(abs(y)-lambda, 0); end ``` 在这个例子中,我们首先读入了一张原始图像,然后添加了高斯噪声。接着,我们设置了TV正则化的参数,包括正则化系数lambda、最大迭代次数max_iter和收敛精度tol。最后,我们调用了TV函数进行图像重建,得到了去噪后的图像。 TV函数中的核心部分是一个迭代过程,它通过求解一个带约束的优化问题来实现图像重建。具体来说,我们定义了一个代价函数,它由两部分组成:图像平滑度和与观测数据的拟合度。使用迭代算法,我们不断优化这个代价函数,直到满足收敛精度为止。 最后,我们还定义了一个soft_threshold函数,用于对图像进行软阈值处理,以实现稀疏性约束。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值