Ternary Interpolation

data(Feldspar)
ggtern(Feldspar,aes(x=Or,y=An,z=Ab)) +
geom_interpolate_tern(aes(value=P.Gpa,color=..level..),binwidth=20) +
geom_point()