Compared with the functionally-based implicit models, the level set models are usually defined by discretely sampling the implicit function values over the volume grid to represent the shapes, which are suitable for scientific computing and engineering. However, the voxel-based level set methods usually lose some features and details to represent the solids due to the finite grid resolutions. In this paper, we propose a more accurate voxel-based implicit modeling method based on discrete level sets. Compared with the traditional level set methods, more than one level set function can be defined in an implicit model with our approach by classifying the cells of the volume grid. The local Boolean compositions conducted over the compound cells are used to depict the interactions of different implicit surfaces represented by corresponding level set functions, and the sharp and thin details caused by Boolean operations are preserved owing to the definition of the compound cells. We show in our experiments that our method is robust to provide the more detailed implicit description of the shapes with less loss of the features than the conventional voxel-based level set methods. INDEX TERMS Compound cells, discrete level sets, local Boolean composition, voxel-based implicit modeling.