Advertisement
Guest User

Untitled

a guest
Apr 24th, 2017
48
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.25 KB | None | 0 0
  1. grid <- raster(extent(one.ym))
  2. res(grid) <- res(one.ym)[1]
  3. proj4string(grid)<-proj4string(one.ym)
  4. gridpolygon <- rasterToPolygons(grid)
  5. writeOGR(obj = gridpolygon,dsn="/home/user/",layer="mygrid_modis.shp", driver = "ESRI Shapefile", overwrite_layer = T
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement