modl-master

所属分类:数值算法/人工智能
开发工具:Python
文件大小:10636KB
下载次数:7
上传日期:2019-06-04 17:07:25
上 传 者工大特仑苏
说明:  这段代码建议我们如何使用深度卷积神经网络(CNN)作为正则化器来解决优化问题。 此代码解决了以下优化问题: argmin_X_AX-B__U 2^2+X-DW(X)_^2_2 ` A`可以是任何度量运算符。这里我们考虑MRI中的平行成像问题 “A”运算符由欠采样掩模、FFT和线圈灵敏度图组成。 ` dw(x)`:它表示使用残差学习CNN的去噪器。
(This code suggest how we can use a deep convolutional neural netwrok (CNN) as a regularizer to solve an optimization problem. This code solves the following optimization problem: argmin_x ||Ax-b||_2^2 + ||x-Dw(x)||^2_2 `A` can be any measurement operator. Here we consider parallel imaging problem in MRI where the `A` operator consists of undersampling mask, FFT, and coil sensitivity maps. `Dw(x)`: it represents the denoiser using a residual learning CNN)

文件列表:
_config.yml (26, 2019-01-10)
demoImage.hdf5 (6238208, 2019-01-10)
model.py (6220, 2019-01-10)
outputImage.jpeg (127559, 2019-01-10)
savedModels (0, 2019-01-10)
savedModels\14Mar_1105pm (0, 2019-01-10)
savedModels\14Mar_1105pm\checkpoint (150, 2019-01-10)
savedModels\14Mar_1105pm\checkpointTst (73, 2019-01-10)
savedModels\14Mar_1105pm\events.out.tfevents.1521086760.argon-hm-p100-compute-3-34.hpc (9272154, 2019-01-10)
savedModels\14Mar_1105pm\model-100.data-00000-of-00001 (1366124, 2019-01-10)
savedModels\14Mar_1105pm\model-100.index (2423, 2019-01-10)
savedModels\14Mar_1105pm\model-100.meta (4913641, 2019-01-10)
savedModels\14Mar_1105pm\model-50.data-00000-of-00001 (1366124, 2019-01-10)
savedModels\14Mar_1105pm\model-50.index (2423, 2019-01-10)
savedModels\14Mar_1105pm\model-50.meta (4913641, 2019-01-10)
savedModels\14Mar_1105pm\model.data-00000-of-00001 (1366124, 2019-01-10)
savedModels\14Mar_1105pm\model.index (2423, 2019-01-10)
savedModels\14Mar_1105pm\model.meta (4913488, 2019-01-10)
savedModels\14Mar_1105pm\modelTst.data-00000-of-00001 (456748, 2019-01-10)
savedModels\14Mar_1105pm\modelTst.index (1072, 2019-01-10)
savedModels\14Mar_1105pm\modelTst.meta (646526, 2019-01-10)
supportingFunctions.py (5253, 2019-01-10)
trn.py (6586, 2019-01-10)
tstDemo.py (2910, 2019-01-10)

# MoDL MoDL: Model Based Deep Learning Architecture for Inverse Problems ### Reference paper: MoDL: Model Based Deep Learning Architecture for Inverse Problems by H.K. Aggarwal, M.P Mani, and Mathews Jacob in IEEE Transactions on Medical Imaging, 2018 Link: https://arxiv.org/abs/1712.02862 IEEE Xplore: https://ieeexplore.ieee.org/document/8434321/ #### What this code do: In the above paper, we propose a technique to combine the power of deep-learning with the model-based approaches. This code suggest how we can use a deep convolutional neural netwrok (CNN) as a regularizer to solve an optimization problem. This code solves the following optimization problem: argmin_x ||Ax-b||_2^2 + ||x-Dw(x)||^2_2 `A` can be any measurement operator. Here we consider parallel imaging problem in MRI where the `A` operator consists of undersampling mask, FFT, and coil sensitivity maps. `Dw(x)`: it represents the denoiser using a residual learning CNN. #### Dependencies We have tested the code in Anaconda python 2.7 and 3.6. The code should work with Tensorflow-1.7 onwards. The dataset is in the hdf5 format. You may require to install hdf5 library in python. In Anaconda you can give following command `conda install h5py` The training code requires tqdm library. It is a nice library that is helpful in tracking the training progress. It can be installed using: `conda install tqdm` In addition, matplotlib is required to visualize the output images. #### Dataset This git repository also includes a single image in the file `demoImage.hdf5`. The testing script `tstDemo.py` will use this image by default and does not require full data download for the testing purpose. We have released the parallel imaging dataset used in this paper. You can download the full dataset from the below link: **Download Link** : https://drive.google.com/file/d/1qp-l9kJbRfQU1W5wCjOQZi7I3T6jwA37/view?usp=sharing This dataset consist of parallel magnetic resonance imaging (MRI) brain data of five human subjects. Four of which are used during training of the model and fifth subject is used during testing. Above link contain fully sampled preprocessed data in numpy format for both training and testing. We also provide coil-sensitity-maps (CSM) pre-computed using E-SPIRIT algorithm. Total file size is 3 GB and contains following arrays: `trnOrg`: This is complex arrary of 256x232x360 containing 90 slices from each of the 4 training subjects. Each slice is of spatial dimension 256x232. This is the original fully sampled data. `trnCSM`: This is a complex array of 256x232x12x360 representing coil sensitity maps (csm). Here 12 represent number of coils. `trnMask`: This is the random undersampling mask to do 6-fold acceleration. We use different mask for different slices. `tstOrg`,`tstCSM`, `tstMask`: These are similar arrays for testing purpose. There are total 1*** testing images. The undersampling mask, for both training and testing cases, is shared corresponding to 6-fold acceleration case. #### How to run the code First, ensure that Tensorflow 1.7 or higher version is installed and working with GPU. Second, just clone or download this reporsitory. The `tstDemo.py` file should run without any changes in the code. On the command prompt `CD` to this `modl` directory i.e. the directory containig `tstDemo.py`. Then you can run the test code using the command: `$python tstDemo.py` from the command prompt. #### Output: ![alt text](https://raw.githubusercontent.com/hkaggarwal/modl/master/outputImage.jpeg) You will need the file `dataset.hdf5` to run the training code `trn.py`. You can download the dataset from the link provided above. Please ignore the future warning by python. #### Files description The folder `savedModels` contain the learned tensorflow model parameters. `tstDemo.py` will use it to read the model and run on the demo image in the file `demoImage.hdf5`. `supportingFunctions.py`: This file contain some supporting functions to calculate the time, PSNR, and read the dataset. `model.py`: This file contain the code for creating the residual learning CNN model as well as the algorithm for conjugate-gradient on complex data. `trn.py`: This is the training code `tstDemo.py`: This is the testing code #### Contact The code is provided to support reproducible research. If the code is giving syntax error in your particular python configuration or some files are missing then you may open an issue or directly email me at jnu.hemant@gmail.com

近期下载者

相关文件


收藏者