Extreme value inference deals with fitting the generalized extreme value distribution model and the generalized Pareto distribution model, which are recently combined to give a single model, namely a two-dimensional non-homogeneous Poisson exceedance point process model. In this paper, we extend the two-dimensional non-homogeneous Poisson process model to include non-stationary effect or dependence on covariates and then derive the likelihood for the extended model.