Uses latitude and longitude coordinates to return the PLSS section geometry from the BLM PLSS web service.

`LL2PLSS(x, y, returnlevel = c("I", "S"))`

## Arguments

- x
longitude coordinates (WGS84)

- y
latitude coordinates (WGS84)

- returnlevel
'S' for "Section" or 'I' for "Intersection" (subsections)

## Value

`sf`

object with geometry and PLSS definition.

## Details

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.

## Note

This function requires the following packages: `httr`

, `jsonlite`

, and `sp`

.

## Author

D.E. Beaudette, Jay Skovlin, A.G. Brown