Uses lattitude and longitude coordinates to return the PLSS section geometry from the BLM PLSS web service.
LL2PLSS(x, y, returnlevel = "I")longitude coordinates
lattitude coordinates
'S' for "Section" or 'I' for "Intersection" (subsections)
list of of PLSS codes and coordinates.
This function takes xy coordinates and returns the PLSS section geometry to the quarter-quarter section. returnlevel options are defaulted 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.