[R-sig-Geo] zonal statistics as table

[Forrest Stevens]

Hi Joseph.. With datasets that large, in my experience you're looking at
spatial subsetting to get something manageable and I've given up on zonal()
as it's just too slow for my needs. Rather, I would do your own feature to
raster conversion on the polygons based on the raster layer so cells
overlap exactly, then for your subsets convert the data to data.table
objects and run your aggregations using that very fast package.

When it comes to zonal statistics R is just a bit cumbersome for anything
larger than the smallest datasets.  I'd be happy to hear anyone else's
thoughts and be corrected if there's a more efficient way to do it!

Sincrerly,
Forrest

--
Forrest R. Stevens
Ph.D. Candidate, QSE3 IGERT Fellow
Department of Geography
Land Use and Environmental Change Institute
University of Florida
www.clas.ufl.edu/users/forrest

On Sat, Sep 27, 2014 at 3:41 PM, Joseph Bechara <jbechara at lri-lb.org> wrote:

>
> Dear all,
>   I'm looking for a R code that do the same work of "zonal statistics as
> table" of ArcIS.I found that the function zonal can do it but the issue is
> that I have large inputs raster and polygons.I'll appreciate if somebody ca
> help me.
> polygons layers contains 9 millions rowsraster is 30 Gb
> Regards,
>
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

	[[alternative HTML version deleted]]