Collectives™ on Stack Overflow
Find centralized, trusted content and collaborate around the technologies you use most.
Learn more about Collectives
Teams
Q&A for work
Connect and share knowledge within a single location that is structured and easy to search.
Learn more about Teams
I'm using the function
rasterToPoints()
from the package
raster
. By default, this function omits the
NA
values.
raster <- raster(matrix(c(1,2,NA, 1,2,3, 1,2,3), nrow = 3, ncol = 3 ))
raster_points <- rasterToPoints(raster, na.rm = FALSE)
Therefore, I put a raster with 200 cells and I get a data.frame with only 150 rows. How can I obtain the point with empty cells in my final data.frame?
You could combine xyFromCell
, values
, and cbind
. This returns the central coordinates of each cell (when you initialized the raster, it assumes a 1 x 1 extent unless defined otherwise) and its value, including NAs
, in matrix form:
cbind(xyFromCell(raster, 1:ncell(raster)), values(raster))
Thanks for contributing an answer to Stack Overflow!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.