vb_test_head

 vb_test_head
 test script for make 3-shell surface from one-shell Curry surface

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 % vb_test_head
0002 % test script for make 3-shell surface from one-shell Curry surface
0003 %
0004 % Copyright (C) 2011, ATR All Rights Reserved.
0005 % License : New BSD License(see VBMEG_LICENSE.txt)
0006 
0007 clear all
0008 
0009 proj_root = [getenv('MATHOME') '/Retino_TY'];
0010 MRI_ID    = 'Retino_TY';
0011 
0012 % Brain & image file for plot
0013 brainfile   = [proj_root '/' MRI_ID '.brain.mat'];
0014 analyzefile = [proj_root '/3d.hdr'];
0015 
0016 % Output smooth surface
0017 headsave{1} = [proj_root '/test_head1.head.mat'];
0018 headsave{2} = [proj_root '/test_head2.head.mat'];
0019 headsave{3} = [proj_root '/test_head3.head.mat'];
0020 
0021 % Brain & image file for plot
0022 Para.brain_file   = brainfile;
0023 Para.analyze_file = analyzefile;
0024 
0025 % --- バネモデル平滑化パラメタ
0026 Para.tangent_rate   = 0.3; % バネ強度
0027 Para.mask_ratio     = 0.5; % マスク強度力係数
0028 Para.mask_threshold = 0.3; % マスク強度閾値
0029 
0030 Para.vstep     = 1  ; % イメージ処理ボクセルサイズ [mm]
0031 Para.Nvertex = 3000; % 出力頂点数
0032 
0033 Para.Nloop = 200;     % バネモデル平滑化繰り返し回数
0034 Para.Nlast = 10;
0035 Para.plot_mode  = 1; % No plot for vb_job_smooth_head
0036 
0037 % --- Skull
0038 Para.Radius = [3 3];% dilasion (expansion)
0039 vb_job_smooth_head(headsave{1}, headsave{2}, Para);
0040 
0041 % --- Scalp
0042 Para.Radius = [2 2];% dilasion (expansion)
0043 vb_job_smooth_head(headsave{2}, headsave{3}, Para);
0044 
0045 % --- Check 3 surface on MRI
0046 vb_plot_slice_head(headsave,analyzefile,brainfile);
0047