Linear interpolation and extrapolation on 2D grid

Linear interpolation and extrapolation on 2D grid

I have a mesh of points on 2D irregular grid. How do I interpolate values
that fall on the boundary and extrapolate the points that fall outside.
The data I have for the points p that I want to extrapolate are
coordinates for the two closest vertices and function values at those
vertices.