We develop a method for estimating the functional surface of a regression coefficient that varies over a complex spatial domain with irregular boundaries, peninsulas and interior holes. The method is motivated by, and applied to, data on housing markets, where the central object of inference is...