Кажется немного слабым по объему данных, но добьется успеха, используя комбинацию функции интерполяции Акимы и функции построения базового контура:
yield <- read.table(textConnection("Lat Lon mxgYield
38.96152878 -123.5578995 7.925971605
39.24150085 -123.6392975 19.6720449
39.52162933 -123.7214966 19.777747
39.80192947 -123.8044968 12.23218451
40.36297989 -123.9729004 15.00058428
40.64371872 -124.058403 1.443492322
38.74369812 -123.1184998 1.938840925"), header=TRUE)
require(akima)
ak.yld <- with( yield, interp(x=Lat, y=Lon, z=mxgYield))
contour(ak.yld)
(выглядит лучше на моем экране, но вы все равно должны увидеть контурную диаграмму для региона с данными.)