Greetings! I'm using the extract command (raster package) with a raster layer (x) and a spatial polygons dataframe (y). I'm using it to extract the sum of raster values by each spatial polygon. However, I'm getting NAs as a result for some polygon IDs, which I know isn't true because that doesn't happen in ArcGIS. Does anyone experienced the same issue? Thanks! Joao
Raster extract by polygon generating NAs
3 messages · João Carreiras, Ben Tupper
Hi, It's hard to know without any reproducible code, but you will want to pay close attention to the value of the na.rm argument to raster::extract(Raster,SpatialPolygons) See ?extract for all the details. Cheers, Ben
On May 16, 2018, at 2:27 PM, Jo?o Carreiras <jmbcarreiras at gmail.com> wrote: Greetings! I'm using the extract command (raster package) with a raster layer (x) and a spatial polygons dataframe (y). I'm using it to extract the sum of raster values by each spatial polygon. However, I'm getting NAs as a result for some polygon IDs, which I know isn't true because that doesn't happen in ArcGIS. Does anyone experienced the same issue? Thanks! Joao
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Ben Tupper Bigelow Laboratory for Ocean Sciences 60 Bigelow Drive, P.O. Box 380 East Boothbay, Maine 04544 http://www.bigelow.org Tick Forecasting: https://eco.bigelow.org/
Thank you Ben, setting na.rm = T did the job. Best wishes Joao
On 16 May 2018 at 19:36, Ben Tupper <btupper at bigelow.org> wrote:
Hi, It's hard to know without any reproducible code, but you will want to pay close attention to the value of the na.rm argument to raster::extract(Raster,SpatialPolygons) See ?extract for all the details. Cheers, Ben On May 16, 2018, at 2:27 PM, Jo?o Carreiras <jmbcarreiras at gmail.com> wrote: Greetings! I'm using the extract command (raster package) with a raster layer (x) and a spatial polygons dataframe (y). I'm using it to extract the sum of raster values by each spatial polygon. However, I'm getting NAs as a result for some polygon IDs, which I know isn't true because that doesn't happen in ArcGIS. Does anyone experienced the same issue? Thanks! Joao
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo Ben Tupper Bigelow Laboratory for Ocean Sciences 60 Bigelow Drive, P.O. Box 380 East Boothbay, Maine 04544 http://www.bigelow.org Tick Forecasting: https://eco.bigelow.org/