Problem: given surface orientation at each image cell, compute depth over some reference plane.
Recall: 
, 
integrate to recover z.
Estimates of p,q are noisy 
use least-squares method to
fit a surface that will minimize 
.
Euler equation:
Discretize, iterate, apply boundary conditions.