Uses latitude and longitude coordinates to return the PLSS section geometry from the BLM PLSS web service.
LL2PLSS(x, y, returnlevel = c("I", "S"))longitude coordinates (WGS84)
latitude coordinates (WGS84)
'S' for "Section" or 'I' for "Intersection" (subsections)
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.