Mapping Bombus species diversity in the US

I'm working on a museum display as part of the outreach activities for our NSF Understanding the Rules of Life mountain epibeenomics grant and decided to illustrate why mountains are especially interesting for bumbles by making a map overlaying species richness with elevation. I just used some fairly heavily post-filtered data downloaded from GBIF for this, so could probably be easily improved on with some more rigorous filtering (although the richness blobs correlate pretty well with Paul Williams' maps, so they are pretty good I think). I determined a gridded richness estimate using R with mostly the sf and raster packages and made the pretty map with rayshader. It was a lot of work to wrangle the data and get all the layers in the correct formats, but came out pretty well I think.