mri去偏场代码

matlab实现的一个MRI图像去偏场(bias correction)代码。

warning off

addpath .\biascorrector
clc

% load data

name=‘brain.img‘;
V  = spm_vol（name）;
ima=（spm_read_vols（V））; 
ima=double（ima）;
ima=ima（::30）‘;
res=abs（[V.mat（11）V.mat（22）V.mat（33）]）;

% filter

[fimaB]=BiasCorrector3Dv1（ima14res）;
  
colormap（gray）
clf
subplot（221）imagesc（ima）;
subplot（222）imagesc（fima）;
subplot（223）imagesc（B）;
h1=histogram（ima）;
h2=histogram（fima）;
subplot（224）plot（h1）hold onplot（h2‘r‘）;
legend（‘Original‘‘Corrected‘）


%************************************************

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        6044  2007-06-28 10:41  biascorrector\biascorrector3Dv1.m
     文件        1683  2007-03-21 21:03  biascorrector\cgradientmodule.c
     文件        6656  2007-03-21 21:03  biascorrector\cgradientmodule.dll
     文件        1636  2007-03-21 12:40  biascorrector\cvolfilter.c
     文件        6656  2007-03-21 12:41  biascorrector\cvolfilter.dll
     文件         386  2007-06-27 21:27  biascorrector\histogram.m
     文件          43  2007-07-16 17:02  biascorrector\mexall.m
     文件      227328  2003-05-10 03:00  biascorrector\spm_bsplinc.dll
     文件       20480  2003-05-10 03:00  biascorrector\spm_bsplins.dll
     文件         471  2002-08-16 15:51  biascorrector\spm_close_vol.m
     文件       11791  2007-05-03 10:58  biascorrector\spm_create_vol.m
     文件         457  2004-10-09 15:46  biascorrector\spm_flip_analyze_images.m
     文件        2049  2002-08-14 03:00  biascorrector\spm_matrix.m
     文件        7760  2005-07-21 12:54  biascorrector\spm_platform.m
     文件        5294  2005-01-11 01:00  biascorrector\spm_read_hdr.m
     文件        2224  2006-02-05 19:08  biascorrector\spm_read_vols.m
     文件      224768  2003-05-10 03:00  biascorrector\spm_slice_vol.dll
     文件        3039  1999-04-27 03:00  biascorrector\spm_type.m
     文件        3510  2005-01-11 01:00  biascorrector\spm_vol.m
     文件        5715  2005-01-11 01:00  biascorrector\spm_vol_ana.m
     文件        4567  2005-01-11 15:56  biascorrector\spm_write_plane.m
     文件        2393  2003-02-26 19:58  biascorrector\spm_write_vol.m
     目录           0  2007-07-16 17:03  biascorrector\
     文件        1663  2007-07-16 17:06  Readme.txt
     文件         493  2007-06-27 23:44  demo3D.m
     文件         558  2007-06-27 21:39  demo2D.m
     文件     3932160  2007-06-27 16:14  brain.img
     文件         348  2007-06-27 16:14  brain.hdr
     文件      198656  2009-08-19 17:35  biascorrector\spm_bsplinc.mexw32
     文件       29184  2009-08-19 17:35  biascorrector\spm_bsplins.mexw32
     文件      190976  2009-08-19 17:35  biascorrector\spm_slice_vol.mexw32
............此处省略0个文件信息