nhd_plus_query function

Select NHDplus features via polygon or circular buffer of coordinate pair