Subset and levelpaths optimizations

This is a minor release with primarily performance improvements and validations for compatibility with ArcPro.

• subset_nhdplus() now validates geometry and ensures all outputs are in NAD83
• subset_nhdplus() queries the NHDPlus database rather than loading then filtering
• get_levelpaths() rewritten for performance. Package now uses data.table.