Monte Carlo (MC) method is a statistical method for simulating photon propagation in media in the optical molecular imaging field. However, obtaining an accurate result using the method is quite time-consuming, especially because the boundary of the media is complex. A voxel classification method is proposed to reduce the computation cost. All the voxels generated by dividing the media are classified into three types (outside, boundary, and inside) according to the position of the voxel. The cla