Get NEI point summary data

National Emissions Inventory (NEI) data is summarized as the sum of all point emissions within the buffer distance of each s2 geohash weighted by the inverse of the distance squared to each emission point.

Usage

get_nei_point_summary(
  x,
  year = c("2020", "2017"),
  pollutant_code = c("PM25-PRI", "EC", "OC", "SO4", "NO3", "PMFINE"),
  buffer = 1000
)

install_nei_point_data(year = c("2020", "2017"))

Arguments

x

a vector of s2 cell identifers (s2_cell object)

year

a character string that is the year of the NEI data

pollutant_code

the NEI pollutant to summarize

buffer

distance from s2 cell (in meters) to summarize data

Value

for get_nei_point_summary(), a numeric vector (the same length as x)

for get_nei_point_data(), a character string path to NEI point data RDS file

Details

The full NEI is conducted every three years, with the latest release being 2020.

The NEI file is downloaded, unzipped, and filtered to observations with a pollutant code of EC, OC, SO4, NO3, PMFINE, or PM25-PRI. Latitude and longitude are encoded as an s2 vector, column names are cleaned, and rows with missing values (including total emissions or emissions units) are excluded.

References

https://www.epa.gov/air-emissions-inventories/national-emissions-inventory-nei

https://www.epa.gov/air-emissions-inventories/2020-national-emissions-inventory-nei-technical-support-document-tsd

Examples

get_nei_point_summary(s2::as_s2_cell(c("8841b399ced97c47", "8841b38578834123")))
#> 8841b399ced97c47 8841b38578834123 
#>     6.302855e-05     1.770889e-05