Voronoi tessellation has been used widely to approximate and model various cellular structures and stochastic patterns appearing in nature. In this work, we present an extended version of the Voronoi tessellation method that partitions the space with certain constraints commonly encountered in either experimental measurements or theoretical models, such as cell volume or size distribution. The new Voronoi method is implemented using an inverse Monte Carlo method. We calculate the topological and statistical properties of tessellated Voronoi cells in several model systems with cell volumes obeying lognormal and bimodal distributions. We also compare the results with those obtained using the conventional PoissonVoronoi method. We observed systematic changes in the topological properties as well as deviations from some established topological relations as the parameters in the constraint were varied. The application of this constrained Voronoi method in microstructure modelling and characterization in poly- and nano-crystalline materials is briefly discussed.