title: st_intersects errors for st_crs(x) == st_crs(y) is not TRUE
body: <pre><code>library(tidyverse) 
  library(tigris) 
  library(sf)
  santacruz &lt;- tracts("CA", "Santa Cruz")
  coords_sf &lt;- locations %&gt;% st_as_sf(coords = c("Longitude", "Latitude"), crs=4269) 
  </code></pre>
  
  <p>This should have the same CRS, but when I try</p>
  
  <pre><code>st_intersects(coords_sf, santacruz)
  </code></pre>
  
  <p>I get</p>
  
  <blockquote>
    <p>Error: st_crs(x) == st_crs(y) is not TRUE</p>
  </blockquote>
  
  <p>I then tried </p>
  
  <pre><code>st_set_crs(santacruz, 4269)
  st_set_crs(coords_sf, 4269)
  st_transform(santacruz, 4269)
  st_transform(coords_sf, 4269)
  </code></pre>
  
  <p>and it doesn't work.  I also tried </p>
  
  <pre><code>st_transform(santacruz, crs = "+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")
  
  st_transform(coords_sf, crs = "+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")
  </code></pre>
  
  <p>No matter what I try with setting CRS and transforming it when I try</p>
  
  <pre><code>st_intersects(coords_sf, santacruz)
  </code></pre>
  
  <p>I get</p>
  
  <blockquote>
    <p>Error: st_crs(x) == st_crs(y) is not TRUE</p>
  </blockquote>
  
  <p>At this point I can't tell if it is something wrong with setting CRS or transforming or the <code>st_intersects</code> function. Thanks,</p>