Start by looking at which watersheds we need polygon shapes for.

Watersheds with paired fish-habitat data:

          Watershed  n percent
             Entiat 61   18.6%
           John Day 75   22.9%
              Lemhi 34   10.4%
              Minam 20    6.1%
  South Fork Salmon 30    9.1%
 Upper Grande Ronde 86   26.2%
          Wenatchee 22    6.7%

Watersheds with CHaMP habitat data

          Watershed  n percent
             Entiat 89   16.8%
           John Day 80   15.1%
              Lemhi 30    5.7%
             Methow 45    8.5%
              Minam 20    3.8%
  South Fork Salmon 43    8.1%
           Tucannon 47    8.9%
 Upper Grande Ronde 87   16.4%
          Wenatchee 45    8.5%
        Yankee Fork 44    8.3%

Read in some HUC boundaries.

NA

Test a plot

Looks good. Let’s save the polygon sf to use later.

Now let’s make a map that we can import into the manuscript

LS0tCnRpdGxlOiAiQ3JlYXRlIENIYU1QIG1hcCIKb3V0cHV0OiBodG1sX25vdGVib29rCmVkaXRvcl9vcHRpb25zOiAKICBjaHVua19vdXRwdXRfdHlwZTogaW5saW5lCi0tLQoKYGBge3Igc2V0dXAsIGVjaG89RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0UsIHJlc3VsdHM9J2hpZGUnfQpsaWJyYXJ5KGtuaXRyKQprbml0cjo6b3B0c19jaHVuayRzZXQoCiAgY29sbGFwc2UgPSBUUlVFLAogIGNvbW1lbnQgPSAiIz4iCikKa25pdHI6Om9wdHNfY2h1bmskc2V0KGVjaG89RkFMU0UsIHdhcm5pbmc9RkFMU0UsIG1lc3NhZ2U9RkFMU0UpCgojIGxvYWQgcGFja2FnZXMgZm9yIGFuYWx5c2lzCmxpYnJhcnkodGlkeXZlcnNlKQpsaWJyYXJ5KGphbml0b3IpCmxpYnJhcnkobWFncml0dHIpCmxpYnJhcnkoc2YpCmxpYnJhcnkoZ2dyZXBlbCkKbGlicmFyeShnZ3NwYXRpYWwpCgojIHNldCBkZWZhdWx0IHRoZW1lIGZvciBnZ3Bsb3QKdGhlbWVfc2V0KHRoZW1lX2J3KCkpCiMgc2V0d2QoJ2FuYWx5c2lzL1Jfc2NyaXB0cycpCmBgYAoKU3RhcnQgYnkgbG9va2luZyBhdCB3aGljaCB3YXRlcnNoZWRzIHdlIG5lZWQgcG9seWdvbiBzaGFwZXMgZm9yLgoKYGBge3J9CmZpc2hfaGFiID0gcmVhZF9yZHMoIi4uL2RhdGEvZGVyaXZlZF9kYXRhL2Zpc2hfaGFiLnJkcyIpCm1vZGVsX3N2eV9kZiA9IHJlYWRfcmRzKCIuLi9kYXRhL2Rlcml2ZWRfZGF0YS9tb2RlbF9zdnlfZGYucmRzIikKYGBgCgpXYXRlcnNoZWRzIHdpdGggcGFpcmVkIGZpc2gtaGFiaXRhdCBkYXRhOgoKYGBge3J9CnRhYnlsKGZpc2hfaGFiLAogICAgICBXYXRlcnNoZWQpICU+JQogIGFkb3JuX3BjdF9mb3JtYXR0aW5nKCkKYGBgCgpXYXRlcnNoZWRzIHdpdGggQ0hhTVAgaGFiaXRhdCBkYXRhCgpgYGB7cn0KbW9kZWxfc3Z5X2RmJGRhdGFbWzFdXSAlPiUKICB0YWJ5bChXYXRlcnNoZWQpICU+JQogIGFkb3JuX3BjdF9mb3JtYXR0aW5nKCkKYGBgCgpSZWFkIGluIHNvbWUgSFVDIGJvdW5kYXJpZXMuCgpgYGB7cn0KaHVjMyA9IHJlYWRfc2YoJ34vT25lRHJpdmUgLSBNZXJjayBTaGFycCAmIERvaG1lLCBDb3JwL0RhdGEvV2F0ZXJzaGVkQm91bmRhcmllcy9XQkRIVTYuc2hwJykgJT4lCiAgc2VsZWN0KFNUQVRFUzpOQU1FKSAlPiUKICBzdF90cmFuc2Zvcm0oY3JzID0gNTA3MCkKaHVjNCA9IHJlYWRfc2YoJ34vT25lRHJpdmUgLSBNZXJjayBTaGFycCAmIERvaG1lLCBDb3JwL0RhdGEvV2F0ZXJzaGVkQm91bmRhcmllcy9XQkRIVTguc2hwJykgJT4lCiAgc2VsZWN0KFNUQVRFUzpOQU1FKSAlPiUKICBzdF90cmFuc2Zvcm0oc3RfY3JzKGh1YzMpKQpodWM1ID0gcmVhZF9zZignfi9PbmVEcml2ZSAtIE1lcmNrIFNoYXJwICYgRG9obWUsIENvcnAvRGF0YS9XYXRlcnNoZWRCb3VuZGFyaWVzL1dCREhVMTAuc2hwJykgJT4lCiAgc2VsZWN0KFNUQVRFUzpOQU1FKSAlPiUKICBzdF90cmFuc2Zvcm0oc3RfY3JzKGh1YzMpKQoKIyBmaWx0ZXIoaHVjNCwKIyAgICAgICAgZ3JlcGwoJ01ldGhvdycsIE5BTUUpKSAlPiUKIyAgIHNlbGVjdChOQU1FKSAlPiUgcGxvdCgpCgpjaGFtcF93dHNkcyA9IGh1YzMgJT4lCiAgZmlsdGVyKE5BTUUgJWluJSBjKCdKb2huIERheScpKSAlPiUKICBzZWxlY3QoTkFNRSkgJT4lCiAgcmJpbmQoaHVjNCAlPiUKICAgICAgICAgIGZpbHRlcihOQU1FICVpbiUgYygnTGVtaGknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICdNZXRob3cnLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICdXZW5hdGNoZWUnLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICdTb3V0aCBGb3JrIFNhbG1vbicsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ1VwcGVyIEdyYW5kZSBSb25kZScpKSAlPiUKICAgICAgICAgIHNlbGVjdChOQU1FKSkgJT4lCiAgcmJpbmQoaHVjNSAlPiUKICAgICAgICAgIGZpbHRlcihOQU1FICVpbiUgYygnTWluYW0gUml2ZXInLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICdTZWNlc2ggUml2ZXInLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICdZYW5rZWUgRm9yaycpKSAlPiUKICAgICAgICAgIHNlbGVjdChOQU1FKSAlPiUKICAgICAgICAgIG11dGF0ZShOQU1FID0gc3RyX3JlbW92ZShOQU1FLCAnIFJpdmVyJCcpKSkgJT4lCiAgcmJpbmQoaHVjNSAlPiUKICAgICAgICAgIGZpbHRlcihOQU1FICVpbiUgYygnRW50aWF0IFJpdmVyJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnTWFkIFJpdmVyJykpICU+JQogICAgICAgICAgc3RfdW5pb24oKSAlPiUKICAgICAgICAgIHN0X3NmKCkgJT4lCiAgICAgICAgICBtdXRhdGUoTkFNRSA9ICdFbnRpYXQnKSkgJT4lCiAgcmJpbmQoaHVjNSAlPiUKICAgICAgICAgIGZpbHRlcihncmVwbCgnVHVjYW5ub24nLCBOQU1FKSkgJT4lCiAgICAgICAgICBzdF91bmlvbigpICU+JQogICAgICAgICAgc3Rfc2YoKSAlPiUKICAgICAgICAgIG11dGF0ZShOQU1FID0gJ1R1Y2Fubm9uJykpICU+JQogIGFycmFuZ2UoTkFNRSkgJT4lCiAgbXV0YXRlKGlkID0gMTpuKCksCiAgICAgICAgIGhhYl9kYXRhID0gVCwKICAgICAgICAgZmlzaF9kYXRhID0gaWZfZWxzZShOQU1FICVpbiUgYygnU2VjZXNoJywgYXMuY2hhcmFjdGVyKHVuaXF1ZShmaXNoX2hhYiRXYXRlcnNoZWQpKSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgVCwgRiksCiAgICAgICAgIGZpc2hfZGF0YSA9IGlmX2Vsc2UoTkFNRSA9PSAnU291dGggRm9yayBTYWxtb24nLCBGLCBmaXNoX2RhdGEpKQoKYGBgCgpUZXN0IGEgcGxvdAoKYGBge3J9CmdncGxvdChkYXRhID0gY2hhbXBfd3RzZHMpICsKICBnZW9tX3NmKCkgKwogIGdlb21fc2YoZGF0YSA9IGNoYW1wX3d0c2RzICU+JQogICAgICAgICAgICBmaWx0ZXIoZmlzaF9kYXRhKSwKICAgICAgICAgIGZpbGwgPSAnYmxhY2snKSArCiAgZ2VvbV9zZl9sYWJlbChhZXMobGFiZWwgPSBpZCksCiAgICAgICAgICAgICAgICBwb3NpdGlvbiA9IHBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMjAwKSkgKwogICAgICAgICAgICAgICAgIyBwb3NpdGlvbiA9IHBvc2l0aW9uX251ZGdlKHggPSA1MDAwLAogICAgICAgICAgICAgICAgIyAgICAgICAgICAgICAgICAgICAgICAgICAgIHkgPSAtNTAwMCkpICsKICBhbm5vdGF0aW9uX25vcnRoX2Fycm93KGxvY2F0aW9uID0gJ3RyJywKICAgICAgICAgICAgICAgICAgICAgICAgIHdoaWNoX25vcnRoID0gJ3RydWUnLAogICAgICAgICAgICAgICAgICAgICAgICAgc3R5bGUgPSBub3J0aF9hcnJvd19mYW5jeV9vcmllbnRlZXJpbmcpICsKICB0aGVtZShheGlzLnRpdGxlID0gZWxlbWVudF9ibGFuaygpKQogIApgYGAKCkxvb2tzIGdvb2QuIExldCdzIHNhdmUgdGhlIHBvbHlnb24gc2YgdG8gdXNlIGxhdGVyLgoKYGBge3J9CndyaXRlX3JkcyhjaGFtcF93dHNkcywKICAgICAgICAgIHBhdGggPSAnLi4vZGF0YS9kZXJpdmVkX2RhdGEvY2hhbXBfd3Rkcy5yZHMnKQpgYGAKCk5vdyBsZXQncyBtYWtlIGEgbWFwIHRoYXQgd2UgY2FuIGltcG9ydCBpbnRvIHRoZSBtYW51c2NyaXB0CgpgYGB7cn0KY2hhbXBfd3RzZHMgPSByZWFkX3JkcygnLi4vZGF0YS9kZXJpdmVkX2RhdGEvY2hhbXBfd3Rkcy5yZHMnKQpyaXZlcl9saW5lcyA9IFFSRmNhcGFjaXR5OjpnZXRfZmxvd2xpbmVzKHN0cmVhbV9vcmRlciA9IDYsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY2hhbXBfd3RzZHMgJT4lCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzdF91bmlvbigpICU+JQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc3RfbWFrZV9ncmlkKG4gPSAxKSkgJT4lCiAgICAgICAgICAgIGZpbHRlcighZ3JlcGwoJ14xNzAxJywgcmVhY2hjb2RlKSwKICAgICAgICAgICAgICAgICAgICFncmVwbCgnXjE3MScsIHJlYWNoY29kZSkpCgojIHJpdmVyX2xpbmVzICU+JQojICAgZmlsdGVyKCFnbmlzX25hbWUgJWluJSBjKCdEZXNjaHV0ZXMgUml2ZXInLAojICAgICAgICAgICAgICAgICAgICAgICAgICAgICdDcm9va2VkIFJpdmVyJywKIyAgICAgICAgICAgICAgICAgICAgICAgICAgICAnUGF5ZXR0ZSBSaXZlcicsCiMgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ05vcnRoIEZvcmsgUGF5ZXR0ZSBSaXZlcicpKSAlPiUKIyAgIGdncGxvdCgpICsKIyAgIGdlb21fc2YoYWVzKGNvbG9yID0gZ25pc19uYW1lKSkKCnd0c2RfbWFwID0gZ2dwbG90KGRhdGEgPSBjaGFtcF93dHNkcykgKwogIGdlb21fc2YoZGF0YSA9IHJpdmVyX2xpbmVzICU+JQogICAgICAgICAgICBmaWx0ZXIoIWduaXNfbmFtZSAlaW4lIGMoJ0Rlc2NodXRlcyBSaXZlcicsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnQ3Jvb2tlZCBSaXZlcicsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnUGF5ZXR0ZSBSaXZlcicsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnTm9ydGggRm9yayBQYXlldHRlIFJpdmVyJykpLAogICAgICAgICAgY29sb3IgPSAnZGFya2JsdWUnKSArCiAgZ2VvbV9zZigpICsKICBnZW9tX3NmKGRhdGEgPSBjaGFtcF93dHNkcyAlPiUKICAgICAgICAgICAgZmlsdGVyKGZpc2hfZGF0YSksCiAgICAgICAgICBmaWxsID0gJ2JsYWNrJykgKwogIGdlb21fc2ZfbGFiZWwoYWVzKGxhYmVsID0gaWQpLAogICAgICAgICAgICAgICAgcG9zaXRpb24gPSBwb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDIwMCkpICsKICBhbm5vdGF0aW9uX25vcnRoX2Fycm93KGxvY2F0aW9uID0gJ3RyJywKICAgICAgICAgICAgICAgICAgICAgICAgIHdoaWNoX25vcnRoID0gJ3RydWUnLAogICAgICAgICAgICAgICAgICAgICAgICAgc3R5bGUgPSBub3J0aF9hcnJvd19mYW5jeV9vcmllbnRlZXJpbmcpICsKICB0aGVtZShheGlzLnRpdGxlID0gZWxlbWVudF9ibGFuaygpKQoKIyBzYXZlIGEgY291cGxlIHZlcnNpb25zCmdnc2F2ZSgnLi4vZmlndXJlcy9DaGFtcFd0ZHMuc3ZnJywKICAgICAgIHd0c2RfbWFwLAogICAgICAgZHBpID0gNjAwLAogICAgICAgaGVpZ2h0ID0gNiwKICAgICAgIHdpZHRoID0gNikKZ2dzYXZlKCcuLi9maWd1cmVzL0NoYW1wV3Rkcy5wbmcnLAogICAgICAgd3RzZF9tYXAsCiAgICAgICBkcGkgPSA2MDAsCiAgICAgICBoZWlnaHQgPSA2LAogICAgICAgd2lkdGggPSA2KQpnZ3NhdmUoJy4uL2ZpZ3VyZXMvQ2hhbXBXdGRzLnRpZmYnLAogICAgICAgd3RzZF9tYXAsCiAgICAgICBkcGkgPSA2MDAsCiAgICAgICBoZWlnaHQgPSA2LAogICAgICAgd2lkdGggPSA2KQpnZ3NhdmUoJy4uL2ZpZ3VyZXMvQ2hhbXBXdGRzLmpwZWcnLAogICAgICAgd3RzZF9tYXAsCiAgICAgICBkcGkgPSA2MDAsCiAgICAgICBoZWlnaHQgPSA2LAogICAgICAgd2lkdGggPSA2KQpgYGAKCg==