Sometimes it happens, sometimes it doesn’t while using read.csv, with that said, you can avoid the problem entirely by running
data.omitted <-na.omit(data). That has cleared up that issue.
But do you happen to know what is the algorithm the contour inter/extrapolation is based on?
So far with fitting trend surfaces, idw, kriging and akima I haven’t been able to replicate it.