Gaussian process (GP) model has been proven effective in modeling photometric redshift as function of flux. However, such model suffers from a number of issues: poor predictibility at extreme redshifts, inability to account for arbitrary uncertainty and proper management of outliers. Instead, a GP mixture model based on spatial dirichlet process solves a number of issues. Motivating by this example, we develop a novel non-parametric spatial modeling technique based on Dirichlet Process - Gaussian Process (DPGP), which models data distribution as mixture of normals, allowing the mixting proportion to change over input space. This modeling approach uses more than one GP to model the data by discovering latent cluster distribution using a Dirichlet Mixture model. We show illustrations of DPGP model using a simulated and a real example from redshift estimation problem.