title: Why is my geopandas installation missing sjoin? AttributeError: 'GeoDataFrame' object has no attribute 'sjoin'
body: <p>I need to use geopandas sjoin on a project. I have successfully created GeoDataFrame objects and used class methods readfile, crs and to_crs. However when it comes to using sjoin I am getting an error: AttributeError: 'GeoDataFrame' object has no attribute 'sjoin'.</p>
  <p>So I thought I would attempt to get sjoin working on a simple example taken directly from the official geopandas website: <a href="https://geopandas.org/en/stable/docs/reference/api/geopandas.GeoDataFrame.sjoin.html" rel="noreferrer">https://geopandas.org/en/stable/docs/reference/api/geopandas.GeoDataFrame.sjoin.html</a></p>
  <pre><code>import geopandas
  
  countries = geopandas.read_file(geopandas.datasets.get_path(&quot;naturalearth_lowres&quot;))
  cities = geopandas.read_file(geopandas.datasets.get_path(&quot;naturalearth_cities&quot;))
  print(countries.head())
  print(cities.head())
  
  cities_w_country_data = cities.sjoin(countries)
  print(cities_w_country_data.head())
  </code></pre>
  <p>Yet I am encountering the same problem. Please could someone explain this error and offer a solution? Thanks.</p>
  <p>The output is as follows:</p>
  <pre><code>C:\Users\mmm\AppData\Local\ESRI\conda\envs\arcgispro-py3_clone\lib\site-packages\geopandas\array.py:85: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.
    aout[:] = out
  C:\Users\mmm\AppData\Local\ESRI\conda\envs\arcgispro-py3_clone\lib\site-packages\geopandas\geodataframe.py:35: ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.
    out = from_shapely(data)
       pop_est      continent  ...  gdp_md_est                                           geometry
  0     920938        Oceania  ...      8374.0  MULTIPOLYGON (((180.00000 -16.06713, 180.00000...
  1   53950935         Africa  ...    150600.0  POLYGON ((33.90371 -0.95000, 34.07262 -1.05982...
  2     603253         Africa  ...       906.5  POLYGON ((-8.66559 27.65643, -8.66512 27.58948...
  3   35623680  North America  ...   1674000.0  MULTIPOLYGON (((-122.84000 49.00000, -122.9742...
  4  326625791  North America  ...  18560000.0  MULTIPOLYGON (((-122.84000 49.00000, -120.0000...
  
  [5 rows x 6 columns]
             name                   geometry
  0  Vatican City  POINT (12.45339 41.90328)
  1    San Marino  POINT (12.44177 43.93610)
  2         Vaduz   POINT (9.51667 47.13372)
  3    Luxembourg   POINT (6.13000 49.61166)
  4       Palikir  POINT (158.14997 6.91664)
  Traceback (most recent call last):
    File &quot;geopandas_sjoin_official_eg.py&quot;, line 8, in &lt;module&gt;
      cities_w_country_data = cities.sjoin(countries)
    File &quot;C:\Users\mmm\AppData\Local\ESRI\conda\envs\arcgispro-py3_clone\lib\site-packages\pandas\core\generic.py&quot;, line 5465, in __getattr__
      return object.__getattribute__(self, name)
  AttributeError: 'GeoDataFrame' object has no attribute 'sjoin'
  
  </code></pre>
  <p>Note:
  I am using Windows 10.
  I am running my code from within a cloned ArcGIS Pro Conda environment.
  I recently installed geopandas using the command &quot;conda install geopandas -y&quot;</p>