coords <-lead$coords
poly <- bbox(as.matrix(coords))
h<- # CHOOSE A SUITABLE VALUE FOR THE GRID-SPACING
grid.pred <- expand.grid(seq(poly[1,1],poly[1,2],h),
                         seq(poly[2,1],poly[2,2],h))