# fetch map unit keys from the bounding-box:
#
# +------------- (-120.8, 37.8)
# | |
# | |
# (-120.95, 37.7) --------------+
b <- c(-120.9,37.7,-120.8,37.8)
m <- MUKEYS_by_ll_bbox(b)
Run the code above in your browser using DataLab