Uses latitude and longitude coordinates to return the PLSS section geometry from the BLM PLSS web service.
LL2PLSS(x, y, returnlevel = c("I", "S"))
sf
object with geometry and PLSS definition.
This function accepts geographic coordinates and returns the PLSS fabric geometry to the quarter-quarter section. returnlevel
defaults to 'I' which returns smallest intersected sectional aliquot geometry, 'S' will return the section geometry of the coordinates. See https://gis.blm.gov/arcgis/rest/services/Cadastral/BLM_Natl_PLSS_CadNSDI/MapServer for details.
This function requires the following packages: httr
, jsonlite
, and sp
.