BrainWeb: Simulated Brain Database使用（原图以及GroundTruth切片可视化）和矢状冠状横断面

---

注：此篇文章参考脑图像的数据预处理
作者：凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/
由于我在做实验时需要脑部横断面，矢状，冠状的原始图片，以及ground truth，
重新梳理BrainWeb数据集的使用。

---

一、BrainWeb数据集简介？

BrainWeb: Simulated Brain Database
http://brainweb.bic.mni.mcgill.ca/brainweb/
数据集选自McGill大学Montreal神经所大脑成像中心的Brain Web反震脑部MR图像数据库。该数据库包含基于两种解剖模型的模拟脑MRI数据：正常和多发性硬化（MS），并且提供了人工合成三种模态(T1-, T2-, PD-)下的三维脑MR图像。图像中包含不同的扫描厚度、噪声以及偏移场，最重要的是且该数据库对于所有脑MR图像均提供了标准分割结果。
用户可对所需图片自定义进行下载：

1. 模态 T1,T2,PD

2. 协议 icmb

3. 脑图像名字：正常

4. 切片厚度 1mm,3mm,5mm,7mm,9mm

5. 噪声水平 pn0=0%;1%,3%,5%,7%,9%

6. 灰度不均匀水平 rf0=0%;20%,40%

二、使用步骤

1.文件下载

原始图片以ms病变脑数据库为例，模态= T1，协议= ICBM，幻像名称=正常，切片厚度= 1mm，噪声= 0％，INU = 0％：

groundtruth同样以以ms病变脑数据库为例

下载好之后就可以处理数据了

2.数据处理

得到第90切片

function g = readrawb(filename, num)
%函数readrawb(filename, num)中的第一个参数filename是欲读取的rawb文件的文件名，第二个参数num就是第多少张。
fid = fopen(filename);
%连续读取181*217*181个数据，这时候temp是一个长度为181*217*181的向量。
%先将rawb中的所有数据传递给temp数组
temp = fread(fid, 181 * 217 * 181);
%然后把它变成了一个181*217*181的数组
images = reshape(temp, 181 , 217, 181);  
%获取图片大小
 [xrange, yrange, zrange] = size(images)