Nonparametric maximum likelihood estimation of features in spatial point processes using Vorono tessellation