automated brain image analysis has a huge potential to objectively help in the diagnosis and followup of many neurological diseases . to perform such analysis tasks ,
one of the first image processing operations is the delimitation of the area of interest . for brain image analysis ,
this operation has received many different names such as brain extraction , skull stripping , or intracranial cavity masking . in each case , the aim is to isolate the brain or intracranial tissues ( depending on area definition ) from the raw image . the accurate estimation of the intracranial volume plays crucial role to obtain robust and reliable normalized measurements of brain structures .
the importance of this operation is reflected by the large number of methods proposed over the past decade [ 213 ] .
many of these methods are based on the modeling of brain intensities ( normally using t1 weighted images due to their excellent contrast for brain tissues ) combined with a set of morphological operations [ 3 , 5 , 12 ] or atlas priors .
for example , the bet ( brain extraction tool ) software from the fsl image processing library   is one of the most used techniques probably due to its accuracy , ease of use , and low computational load .
other techniques like 3dintracranial , hybrid watershed algorithm ( hwa ) , or brain surface extractor ( bse )   have been also widely used .
intracranial cavity extraction can also be obtained indirectly as part of a full modeling of brain intensities using a parametric model such as that done in statistical parametric mapping ( spm )   or vbm8 ( http:/dbm.neuro.uni - jena.de / vbm/ ) software packages . over the last decade
, methods have been proposed to automatically measure the intracranial cavity volume ( icv ) by using nonlinear registration atlas - based approaches [ 15 , 16 ] .
more recent works of special interest for the brain extraction problem are methods like maps   and beast .
maps uses multiple nonlinear registrations followed by a voxel - wise label fusion while beast uses a single linear registration in combination with nonlocal patch - based label fusion . both techniques scored well on the loni segmentation validation engine ( sve )   comparison for brain extraction ( see http://sve.bmap.ucla.edu/archive/ )
although maps has a much larger computational load compared to beast . in this paper we present an extension of the beast methodology where we aim to improve the accuracy while reducing the computational load .
the main contributions of the proposed method are threefold : first , the use of a new pipeline for the multiatlas library construction for improved normalization between template library subjects ; second , the use of a new bilateral patch similarity measure to better estimate pattern similarities ; and finally , a blockwise labeling approach that enables significant savings in computational cost and imposing at the same time a regularization constraint that increases the method 's accuracy .
since the method proposed in this paper is based on the use of a library of prelabeled cases to perform the segmentation process , we will first describe the template library construction and then present the proposed method .
a library of manually labeled templates was constructed using subjects from different publically available datasets . to include as large age range as possible ,
thirty normal subjects ( age range : 2475 years ) were randomly selected from the open access ixi dataset ( http://www.brain-development.org/ ) .
this dataset contains images from nearly 600 healthy subjects from several hospitals in london ( uk ) .
both 1.5  t ( 7 cases ) and 3  t ( 23 cases ) images were included in this dataset .
3  t images were acquired on a philips intera 3  t scanner ( tr = 9.6  ms , te = 4.6  ms , flip angle = 8 , slice thickness = 1.2  mm , volume size = 256  256  150 , voxel dimensions = 0.94  0.94  1.2  mm ) .
1.5  t images were acquired on a philips gyroscan 1.5  t scanner ( tr = 9.8  ms , te = 4.6  ms , flip angle = 8 , slice thickness = 1.2  mm , volume size = 256  256  150 , voxel dimensions = 0.94  0.94  1.2  mm ) . 

nine patients with alzheimer 's disease ( age range = 7580 years , mmse = 23.7  3.5 , cdr = 1.1  0.4 ) scanned using a 1.5  t general electric signa hdx mri scanner ( general electric , milwaukee , wi ) were randomly selected .
this dataset consisted of high resolution t1-weighted sagittal 3d mp - rage images ( tr = 8.6  ms , te = 3.8  ms , ti = 1000  ms , flip angle = 8 , slice thickness = 1.2  mm , matrix size = 256  256 , voxel dimensions = 0.938  0.938  1.2  mm ) .
these images were downloaded from the brain segmentation testing protocol   website ( https://sites.google.com/site/brainseg/ ) while they belong originally to the open access oasis dataset ( http://www.oasis-brains.org/ ) . 

ten infant datasets were also downloaded from the brain segmentation testing protocol   website ( https://sites.google.com/site/brainseg/ ) .
these data were originally collected by gousias et al .   and are also available at http://www.brain-development.org/ ( this dataset is property of the imperial college of science technology & medicine and has been used after accepting the license agreement ) .
the selected 10 cases are from the full sample of 32 two - year - old infants born prematurely ( age = 24.8  2.4 months ) .
sagittal t1 weighted volumes were acquired from each subject ( 1.0  t phillips hpq scanner , tr = 23  ms , te = 6  ms , slice thickness = 1.6  mm , matrix size = 256  256 , voxel dimensions = 1.04  1.04  1.6  mm resliced to isotropic 1.04  mm ) .
downloaded images from the different websites consisted of raw images with no preprocessing and no intracranial cavity masks were supplied with these data . to generate the template library ,
all images in the database were denoised using the spatially adaptive nonlocal means ( sanlm ) filter   to enhance the image quality .
the sanlm filter can deal with spatially varying noise levels across the image without the need of explicitly estimating the local noise level which makes it ideal to process data with either stationary or spatially varying noise ( as in the case of parallel imaging ) in a fully automatic manner . to further improve the image quality ,
the n4 method is an incremental improvement of the n3 method   that has been implemented in the itk toolbox   and has proven to be more efficient and robust . in order to perform the segmentation process , templates and the subject to be segmented
therefore , a spatial normalization based on a linear registration to the montreal neurological institute ( mni 152 ) space was performed using ants software .
the resulting images in the mni space have a size of 181  217  181 voxels with 1  mm voxel resolution .
as the proposed method is based on the estimation of image similarities using intensity - derived measures , every image in the library must be intensity normalized in order to make the intensity distributions comparable among them .
we use a tissue - derived approach to force mean intensities of white matter ( wm ) , grey matter ( gm ) , and cerebrospinal fluid ( csf ) to be as similar as possible across subjects of the library in a similar manner as done by ltjnen et al . .
for this purpose , mean values of csf , gm , and wm tissues were estimated using the trimmed mean segmentation ( tms ) method   which robustly estimates the mean values of the different tissues by excluding partial volume voxels from the estimation jointly with the use of an unbiased robust mean estimator .
such estimation was performed using only voxels within the standard brain mask area of mni 152 template to minimize the inclusion of external tissues .
finally , a piecewise linear intensity mapping [ 25 , 27 ] was applied ensuring that wm had an average intensity of 250 , gm of 150 , and csf of 50 ( see figure 1 ) . as commented previously ,
there is no standard definition of what should be included in brain or intracranial masks ( it all depends on what you are looking for ) . in beast , the mask definition included the following tissues : all cerebral and cerebellar white matters , all cerebral and cerebellar gray matters , csf in ventricles ( lateral , 3rd , and 4th ) and the cerebellar cistern , csf in deep sulci and along the surface of the brain and brain stem , the brainstem ( pons , medulla),internal brain blood vessels . in the present work we extended that definition by including all external csf ( thus covering total csf of ic ) and therefore selecting most of the intracranial cavity volume .
we have not included other intracranial tissues in our mask definition such as dura , exterior blood vessels , or veins because they are normally of no interest for brain analysis .
this mask definition has been traditionally used to estimate the total intracranial volume ( tiv ) in many methods such rbm , spm8 , or vbm8 methods to normalize brain tissue volumes [ 29 , 30 ] as it is expected to be nearly constant in each subject during the adult lifespan .
all cerebral and cerebellar white matters , all cerebral and cerebellar gray matters , csf in ventricles ( lateral , 3rd , and 4th ) and the cerebellar cistern , csf in deep sulci and along the surface of the brain and brain stem , the brainstem ( pons , medulla ) , internal brain blood vessels .
to generate the template masks we followed a similar approach as described in beast paper since full manual labeling was too time consuming and error prone as discussed in eskildsen et al . .
all template images in the library were automatically segmented using beast software to have an initial mask .
conditional mask dilation ( only over csf voxels ) was applied to include external csf not already included in the beast mask .
finally , all the images were manually corrected by an expert on brain anatomy using the itk - snap software   to remove segmentation errors . in figure 2
we show an example of our mask definition compared to beast definition for a patient with alzheimer 's disease .
to further increase the number of available priors on the library , all the cases were flipped along the midsagittal plane using the symmetric properties of the human brain yielding a total number of 98 labeled templates ( original and flipped ) as done in beast paper .
compared to beast template library creation , the main differences are the use of a denoising method to improve data quality , the use of a different registration method ( ants instead of animal ) , and the application of different intensity normalization method .
the scheme of the template library construction pipeline is summarized in figure 3 . while the beast technique was designed to improve downstream analysis such as the assessment of cortical thickness
, our proposed method has extended the mask definition to include extracerebral spinal fluid as it can be interesting to obtain normalized brain and tissue specific volumes in many neurological diseases such as alzheimer or parkinson .
we will refer our proposed method as nice ( nonlocal intracranial cavity extraction ) . since the method proposed in this paper is an evolution of the beast brain masking method
, we refer the reader to the original paper for the detailed method overview . here
, we summarize the nice method and present the main improvements introduced to increase the method performance . to segment a new case , it must be first preprocessed using the proposed normalization pipeline ( see section 2.1 and figure 2 ) so that the new case is spatially aligned with the template library and to ensure that it has the same intensity characteristics . in the classical nonlocal means label fusion technique proposed by coup et al .
, for each voxel xi from the new image to be segmented the method estimates the final label by performing a weighted label fusion v(xi ) of all surrounding samples inside the search area vi from n subjects selected from the library : 
 ( 1)v(xi)=s=1njviw(xi , xs , j)ls , js=1njviw(xi , xs , j ) , 
 where ls , j is the label from the voxel xsj at the position j in the template library subject s and w(xi , xs , j ) is the weight calculated by patch comparison which is computed depending on the similarity of the surrounding patch for xi and for xs , j this weight is estimated as follows : 
 ( 2)w(xi , xs , j)={exp||p(xi)p(xs , j)||22/h2,if  ss < th0,otherwise , 
 where p(xi ) is the patch around the voxel xi , p(xs , j ) is the patch around the voxel xj in the templates , and ||||2 is the normalized l2-norm ( normalized by the number of elements ) calculated by the distance between each pair of voxels from both patches p(xi ) and p(xs , j ) and modulated by h parameter . if ss ( structural similarity index   between patches ) is less than a threshold th , then w is not computed thus avoiding unneeded computations . the structural similarity index ss is calculated as follows : 
 ( 3)ss=2is , ji2+s , j22is , ji2+s , j2 , 
 where  and  are the mean and standard deviation of the patches surrounding xi and xs , j at location j of the template s. finally , the final label l(xi ) is computed as 
 ( 4)l(xi)={1v(xi)0.50v(xi)<0.5 . 
 in this paper
first , we make use of the fact that all the images are registered to a common space and therefore a locality principle can be used , assuming that samples that are spatially closer are likely to be more similar in their labels . however , this locality principle is limited by residual anatomical variability and registration errors in the template library space . therefore , we redefined the similarity weight to take into account not only intensity similarity but also spatial patch proximity : 
 ( 5)w(xi , xs , j)={exp[(||xixj||2/d)+(||p(xi)p(xs , j)||22/h ) ] ,  if  ss <
th0 ,  otherwise , 
 where xi and xj are the coordinates of patch centers and d is normalization constant . we set d = 8  mm experimentally which curiously coincides with the typical gaussian blurring kernel size normally used on voxel based morphometry ( vbm ) to deal with registration error and subject anatomical variability .
this approach shares some similarities to the bilateral filter proposed by tomasi and manduchi   for image denoising .
we experimentally set the threshold th to 0.97 instead of 0.95 as used in beast ( this difference can be explained due to the use of filtered data and a different intensity normalization method ) . also a comment about h parameter of ( 5 ) is required since it plays a major role in the weight computation process . in   this value was set to 
 ( 6)h(xi)=argminxs , l||p(xi)p(xj , s)||2+ , 
 where  is a small constant to ensure numerical stability . in    was set to 1 but we found experimentally that a value of 0.1 produced better results in the proposed method possibly due to the improved intensity normalization .
classical nonlocal label fusion works in a voxelwise manner which sometimes results in a lack of regularization on the final labels .
given that we wish to segment a continuous anatomical structure , some level of regularization can be used as a constraint achieved by a blockwise vote scheme , similar to the one used by rousseau et al .  
this bilateral blockwise vote is computed as follows : 
 ( 7)v(b(xi))=s=1njviw(xi , xs , j)l(b(ys , j))k=0ms=1njviw(xi , xs , j ) , 
 where b(xi ) is a 3d region which is labeled at the same time .
finally , the vote for the voxel xi is obtained in an overcomplete manner by averaging over all blocks containing xi and the label l(xi ) is decided as in ( 4 ) . with overlapping blocks ,
it is worth noting that the distance between adjacent block centers can be increased to be equal or higher than 2 voxels .
therefore , we can obtain important accelerating factors compared to the voxelwise version of the algorithm ( e.g. , for a distance equal to 2 voxels in all three directions a speedup factor of 2 = 8 can be obtained ) .
the described approach is used within the multiresolution framework as described in the beast paper .
these datasets can be classified two categories : ( a ) those that were used to measure the accuracy of the different methods compared and ( b ) those used to measure their reproducibility . 

, we used the template library dataset by using a leave - one - out ( loo ) cross validation .
each of the 49 ( nonflipped ) library images was processed with the remaining images as priors ( after removing the current case and its flipped version ) .
independent validation dataset . to avoid any factor associated to our ic mask definition that could bias the comparison of the compared methods
therefore , we performed a validation using an independent dataset available in the online segmentation validation engine ( sve ) .
this dataset consists of 40 t1w mri scans and its associated manual labels ( 20 males and 20 females ; age range 1940 ) .
this high - resolution 3d spoiled gradient echo ( spgr ) mri volume was acquired on a ge 1.5  t system as 124 contiguous 1.5  mm coronal brain slices ( tr range 10.0  ms12.5  ms ; te range 4.22  ms4.5  ms ; fov 220  mm or 200  mm ; flip angle 20 ) with in - plane voxel resolution of 0.86  mm ( 38 subjects ) or 0.78  mm ( 2 subjects ) .
although the accuracy of a method is very important , another important feature is its reproducibility .
indeed , the capability to detect changes induced by the pathology in a consistent manner is a key aspect . to measure the reproducibility of the different compared methods
, we used the reproducibility dataset of the brain segmentation testing protocol website ( https://sites.google.com/site/brainseg/ ) .
this dataset consists of a test - retest set of 20 subjects scanned twice in the same scanner and sequence ( sss ) and another set of 36 subjects scanned twice on different scanner and different magnetic field strength ( dsdf ) ( 1.5 and 3 tesla ) . 

sss dataset . to measure the reproducibility of the different methods compared on the same subjects and using the same mri scanner
, we used a subset of the oasis ( www.oasis-brains.org ) dataset consisting in 20 subjects ( age = 23.4  4.0 years , 8 females ) who were scanned using the same pulse sequence two times ( 1.5  t siemens vision scanner , tr = 9.7  ms , te = 4  ms , ti = 20  ms , flip angle = 10 , slice thickness = 1.25  mm , matrix size = 256  256 , voxel dimensions = 1  1  1.25  mm resliced to 1  1  1  mm , averages = 1 ) . 
 dsdf dataset . to determine the consistency of the segmentations when different mri scanners and different magnetic field strength were used ,
36 adult subjects were scanned using two mri scanners ( 1.5  t and 3.0  t general electric signa hdx scanner ) , mean interscan interval between 1.5  t and 3  t scanner = 6.7  4.2 days ) . to study the impact of the method parameters , an exhaustive search of the optimum values
each one of the 49 subjects in the library was processed using the remaining cases of the library as priors and the resulting segmentation was compared to the manual labeling . to measure segmentation accuracy ,
method parameters such as patch size and search area were set as in beast method while an exhaustive search for the optimal number of templates n used for the segmentation process was carried out ( see figure 4 ) .
this search demonstrated that the segmentation accuracy stabilizes around n = [ 2030 ] range which is in good agreement with previous results from beast .
we decided to use n = 30 as default value given the reduced computational cost of the proposed method .
another parameter of our proposed block - based approach is the spacing between adjacent blocks which jointly with patch size defines the degree of overlap between blocks .
we observed experimentally that the optimal value for that parameter was 2 voxels in all 3 dimensions since the resulting accuracy was virtually the same from full overlap ( 1 voxel spacing ) while computation time was greatly reduced .
both beast and vbm8 methods were selected because of their public accessibility and because they are among the highest ranking methods on the online segmentation validation engine website   ( http://sve.bmap.ucla.edu/archivel/ ) . to ensure a fair comparison all three methods
, we used the same preprocessing pipeline with the exception of the intensity normalization step ( i.e. , using ants registration to ensure the same image space and the same homogenization and filtering to ensure the same image quality ) . in this way , only the labeling process was evaluated eliminating other sources of variability . both nice and beast
were run with the same number of preselected templates ( n = 30 ) to ensure a fair comparison .
we used release 435 of vbm8 , which was the latest version at the time of writing . to compare the segmentation results of the different methods , several quantitative metrics
in table 1 , the average dice coefficient , sensitivity , and specificity for all 49 cases of loo dataset using the different methods compared are provided .
results for all the cases together and separated by dataset subtype are provided ( alzheimer 's disease ( ad ) , normal infants ( infant ) , and normal adult subjects ( adult ) ) .
as can be noticed , nice method obtained the best results in all the situations .  
table 2 shows the statistical significance of these differences ( paired t - test ) .
intracranial cavity volume is normally used to normalize brain tissue volumes to provide a tissue measure independent of head size .
therefore , the ability of the compared methods to provide an accurate icv estimation has to be assessed . to this end
, volume estimations using the different compared methods were obtained and compared to gold standard manual volumes .  
figure 5 shows the automatic versus manual volume correlation for all the compared methods and dataset used . as can be noticed , the nice method had highest overall correlation ( 0.976 ) while beast and vbm8 had 0.923 and 0.778 , respectively . in figure 6
, a visual comparison of the segmentation results of three examples belonging to the three different subject populations can be performed . to perform an independent validation of the compared methods ,
the sve web service allows the comparison of results with hand - corrected brain masks .
as done in beast and maps papers , we used the brain masks provided by the sve website which included all the internal ventricular csf and some external sulcal csf .
although this mask definition slightly differs from our mask definition ( not all csf was included ) , this does not represent a problem for the method 's comparison since all the methods shared the same references .
validation of nice using the sve test dataset resulted in a mean dice of 0.9819  0.0024 ( see http://sve.bmap.ucla.edu/archivel/ ) . at the time of writing , this result was the best ( p < 0.01 ) of all the methods published on the website .
beast had a dsc of 0.9781  0.0047 and vbm8 obtained a dsc of 0.9760  0.0025 .
sensitivity and specificity results are also included in table 3 . a visual representation of false positive and false negative as supplied by the website is presented at figure 7 . in table 4 ,
the average icv differences for the different methods and datasets is provided . as can be noticed , nice method obtained the most reproducible results in all situations . for the sss dataset experiment ( test - retest ) ,
nice significantly improved beast method while these differences were not significant for vbm8 method . for the dsdf dataset experiment ,
volume differences were higher than in the previous experiment . in this case , nice was found to yield significantly improved estimates ( p < 0.05 ) compared to the two other methods .
nice method took around 4 minutes ( nice was implemented as a multithreaded mex c file ) , beast method took around 25 minutes ( we have to note that no multithreading optimizations were used here ) , and vbm8 took around 8 minutes on average ( in this time it was also included the different tissue segmentations ) .
all the experiments were performed using matlab 2009b 64 bits ( mathworks , inc . ) on a desktop pc with an intel core i7 with 16  gb ram running windows 7 .
however , it is worth noting that if we reduce the number of selected templates to 10 cases , we can reduce the processing time to less than 1 minute with only a small reduction of the segmentation accuracy ( 0.9911 to 0.9901 in the loo accuracy experiment ) .
we have presented a new method for intracranial cavity extraction that outperforms related state - of - the - art methods and a previously proposed method ( beast ) by our group both in terms of accuracy and reduced computational load .
in addition , we demonstrated that the new proposed method is more robust in terms of measurement reproducibility .
this last point is of special interest since in many cases we are not only interested in the specific brain volume at one time point but in its evolution in a longitudinal study .
in addition , vbm8 was found to be almost as reproducible as nice but at the expense of introducing larger systematic errors on the segmentations .
the high level of reproducibility of vbm8 may be explained by the fact that it uses a single template and thus a more deterministic pipeline is applied .
also the fact that it operates at 1.5  mm resolution introduces a blurring effect which increases the method reproducibility at the expense of the accuracy .
the increased reproducibility / accuracy of our proposed method may have a significant impact on the brain image analysis methods by increasing their sensitivity to detect subtle changes produced by the disease .
while the advantage of nice in segmentation accuracy of 0.9911 versus 0.9762 for vbm8 may appear small when compared over the three datasets evaluated , it is statistically significant and corresponds to more than a 2-fold reduction in error , from 2.38% to 0.89% . in a large volume such as the intracranial cavity ( 1500  cc )
, this reduction in error can represent a volume of approximately 20  cc , a nonnegligible amount .
the improvement over beast is smaller ( 35% ) but still statistically significant . when evaluated on the sve dataset , the nice yields a dice overall of 0.9819 , while beast and vbm8 yield 0.9781 and 0.9760 , corresponding to 20% and 32% less error on average , respectively .
the improved results of nice over beast can be understood thanks to improvements on two parts of the proposed method .
first , improvements on template library construction such as the improved intensity normalization yield more coherent and better defined priors .
one limitation of the first part of our validation is in the use of manually corrected masks that may induce a favourable bias toward beast and nice .
however , after the conditional dilation and manual correction steps , almost all edge voxels were modified , thus minimizing any bias .
second , the blockwise and new bilateral label fusion scheme results in more regular and accurate segmentations .
the advantages of using a 3d blockwise approach in comparison to the previously used voxelwise are twofold : first , the fact that we label together the whole block imposes an intrinsic regularization which forces connected voxels to have similar labels and second if a space between block centers is used , a significant speed - up factor can be obtained in comparison with the voxelwise version .
finally , the new similarity measure using spatial distance weighting takes into account a locality principle that favors the contribution of closer patches by assuming that after linear registration similar structures are close in a similar manner as done for the well - known bilateral filter for image denoising .
it is also worth noting that the segmentation accuracy depends on the preexistence of similar local patterns within the library . in our method
, we do not need to have totally similar templates to the case to be segmented within the library since it is able to find locally similar patterns from different templates in the library .
however , it is also true that if some specific pattern is not present in the library , it will not be correctly identified and therefore the resulting segmentation will be incorrect .
this risk is normally reduced when using nonlinear registrations at the expense of a much higher computational load and the introduction of interpolation artifacts in both images and associated labels .
however , this issue can be solved more efficiently ( mainly in terms of computational cost ) by increasing library size with uncommon cases and their associated corrected manual labels making it unnecessary to perform costly nonlinear registrations ( but making necessary the manual label correction of new library cases ) . we experimentally found that increasing the size of the library just using the symmetric versions of the original library improved the segmentation results as previously reported in the beast paper .
finally , it is also possible to construct disease specific libraries ( as done for templates in spm ) maximizing the likelihood to find suitable matches for the segmentation process or to improve templates preselection by adding extra information such as age or sex which could help to find optimal matches ( especially useful when the library size will grow ) .
we have presented an improved method to perform intracranial cavity extraction that has been shown to be fast , robust , and accurate .
the improvements proposed have been shown to increase segmentation quality and reduce the computational load at the same time ( the proposed method is able to work in a reasonable time of approximately 4 minutes ) .
we plan to make the nice pipeline publicly accessible through a web interface in the near future so everybody can benefit from its use independently of their location and local computational resources
. finally , the usefulness of the proposed approach to provide accurate icv based normalization brain tissue measurements has to be addressed on future clinical studies .
 automatic and accurate methods to estimate normalized regional brain volumes from mri data are valuable tools which may help to obtain an objective diagnosis and followup of many neurological diseases . to estimate such regional brain volumes , the intracranial cavity volume ( icv ) 
 is often used for normalization . 
 however , the high variability of brain shape and size due to normal intersubject variability , normal changes occurring over the lifespan , and abnormal changes due to disease makes the icv estimation problem challenging . in this paper 
 , we present a new approach to perform icv extraction based on the use of a library of prelabeled brain images to capture the large variability of brain shapes . to this end , an improved nonlocal label fusion scheme based on beast technique 
 is proposed to increase the accuracy of the icv estimation . 
 the proposed method is compared with recent state - of - the - art methods and the results demonstrate an improved performance both in terms of accuracy and reproducibility while maintaining a reduced computational burden . 
1. Introduction
2. Materials and Methods
3. Experiments
4. Results
5. Discussion
6. Conclusion