To model the true shape of MRI brain images, automatically classified Ti-weighted 3D MRI images (gray matter, white matter, cerebrospinal fluid, scalp/bone and background) are utilized for simulation of grayscale data and imaging artifacts. For each class, Gaussian distribution of grayscale values is assumed, and mean and variance are computed from grayscale images. A random generator fills up the class images with Gauss-distributed grayscale values. Since grayscale values of neighboring voxels are not correlated, a Gaussian low-pass filtering is done, preserving class region borders. To simulate anatomical variability, a Gaussian distribution in space with user-defined mean and variance can be added at any user-defined position. Several imaging artifacts can be added: 1 . to simulate partial volume effects, every voxel is averaged with neighboring voxels if they have a different class label; 2. a linear or quadratic bias field can be added with user-defined strength and orientation; 3. additional background noise can be added; and 4. artifacts left over after spoiling can be simulated by adding a band with increasing/decreasing grayscale values. With this method, realistic-looking simulated MRI images can be produced to test classification and segmentation algorithms regarding accuracy and robustness even in the presence of artifacts.Evaluation of classification or segmentation results and validation of automatic algorithms is difficult in the case of nonsimulated MRI images because the true borders of tissue types are unknown. Therefore, evaluation is often done by comparing automatically obtained results with results supplied by a physician. This has the disadvantage of subjectivity and inter-or intra-observer variability. Furthermore, 3D brain images have very complex region borderlines. A complete interactive delineation is practically impossible.In this approach, Ti-weighted 3D-FLASH MRI grayscale images are utilized to generate a software phantom for evaluation and validation purposes. The images are classified into gray matter (GM), white matter (WM), cerebrospinal fluid (CSF), scalpfbone (SB) and background (BG) using a simple thresholding method. To model the true shape of MRI brain images, these automatically classified 3D MIRI images are utilized for simulation of grayscale data and imaging artifacts.The main topics of this paper are the detailed descriptions of the methods developed for classification and simulation of Ti-weighted 3D-MRI images and possible artifacts. Examples of applications of this software phantom are briefly described.