modis-sinu-grid1.R

## the projected extent
prj <- "+proj=sinu +a=6371007.181 +b=6371007.181"
ex <- extent(rgdal::project(cbind(c(-180, -180, 180, 180), c(90, 0, 0, -90)), prj))
#the MODIS tiles
mtiles <- raster(ex, ncols = 35, nrows = 17, crs = prj)

## for real use you wouldn't actually populate the grid here
## raster's cell logic can be used without ever filling the data
mtiles[] <- seq(ncell(mtiles))
#plot(setValues(mtiles, seq(ncell(mtiles))))


library(sf)
pmtiles <- st_as_sf(as(mtiles, "SpatialPolygonsDataFrame"))
pmtiles$row <- rowFromCell(mtiles, seq(ncell(mtiles)))
pmtiles$col <- colFromCell(mtiles, seq(ncell(mtiles)))
plot(pmtiles[c("row", "col")], asp = "")

plot(st_geometry(pmtiles), asp = "")
data("wrld_simpl", package = "maptools")
plot(spTransform(wrld_simpl, prj), add = TRUE)