Sven, It looks like myPointGrid_utm has more points that rows in the associated dataframe. Could you check that? For example, with nrow(coordinates(myPointGrid_utm)) nrow(myPointGrid_utm at data) Best, Virgilio