Spatial analyses involving binning require that every bin have the same area, but this is impossible using a rectangular grid laid over the Earth or over any projection of the Earth. Discrete global grids use hexagons, triangles, and diamonds to overcome this issue, overlaying the Earth with equally-sized bins. This package provides utilities for working with discrete global grids, along with utilities to aid in plotting such data.
Spatial Analysis Done Right
You want to do spatial statistics, and it's going to involve binning.
Binning with a rectangular grid introduces messy distortions. At the macro-scale using a rectangular grid does things like making Greenland bigger than the United States and Antarctica the largest continent.
But this kind of distortion is present no matter what the resolution is; in fact, it shows up whenever you project a sphere onto a plane.
What you want are bins of equal size, regardless of where they are on the globe, regardless of their resolution.
dggridR solves this problem.
dggridR builds discrete global grids which partition the surface of the Earth into hexagonal, triangular, or diamond cells, all of which have the same size. (There are some minor caveats which are detailed in the vignettes.)
(Naturally, you can use much smaller cells than those shown in the image above.)
This package includes everything you need to make spatial binning great again.
Many details and examples are included in the vignette.
dggridR is available from CRAN via:
If you want your code to be as up-to-date as possible, you can install it using:
library(devtools) #Use `install.packages('devtools')` if need be install_github('r-barnes/dggridR', vignette=TRUE)
Your analysis could be as easy as this:
library(dggridR)library(dplyr)dggs <- dgconstruct(spacing=1000, metric=FALSE, resround='down')#Load included test data setdata(dgquakes)#Get the corresponding grid cells for each earthquake epicenter (lat-long pair)dgquakes$cell <- dgtransform(dggs,dgquakes$lat,dgquakes$lon)#Get the number of earthquakes in each equally-sized cellquakecounts <- dgquakes %>% group_by(cell) %>% summarise(count=n())
In R, typing
will bring up many examples.
Many different grid resolutions are available for many different grids. The following chart shows the number of cells, their area, and statistics regarding the spacing of their center nodes for the ISEA3H grid type.
|Res||Number of Cells||Cell Area (km^2)||Min||Max||Mean||Std|
The code in the 'src' directory is based off of DGGRIDv6.2b by Kevin Sahr.
However, Richard Barnes has made some significant alterations. These include:
The package relies on several libraries, as noted in the Licensing section below.
This package uses the following libraries:
clipper: The clipper library has been released under the Boostv1 license by Angus Johnson
dggrid: Kevin Sahr has released dggrid as a "public domain software program"
proj4lib: Gerald Evenden has released this code into the public domain. More recent versions of the library, not used here, are under an MIT-style license.
shapelib: Frank Warmerdam has released shapelib under an MIT-style license. There is also an LGPL licensing option which I have chosen not to use.
This package, and all code and documentation not otherwise mentioned above
(essentially anything outside the
src/ directory of this package) are released
under the MIT (Expat) license, as stated in the
LICENSE file. The
file exists for use with CRAN.
This package should operate in the manner described here, in the package's
main documentation, and in Kevin Sahr's
dggrid documentation. Unfortunately,
none of us are paid enough to make absolutely, doggone certain that that's the
case. That said, if you find bugs or are seeking enhancements, we want to hear
Please cite this package as:
long double output formatting issue that gives bad results when
code is compiled using MingW on certain Windows systems.
Adds a debugging capability to the R code to make it easier to get diagnostics from users.
wrapcellsargmuent to grid generation functions: makes plotting easier
savegridargument to make it easier to output grids for other purposes
First official release