code: def test_sjoin_empty_geometries(self):
      # https://github.com/geopandas/geopandas/issues/944
      empty = GeoDataFrame(
          geometry=[GeometryCollection()] * 3, crs=self.crs
      )
      df = sjoin(
          pd.concat([self.pointdf, empty]), self.polydf, how="left"
      )
      assert df.shape == (24, 8)
      df2 = sjoin(
          self.pointdf, pd.concat([self.polydf, empty]), how="left"
      )
      assert df2.shape == (21, 8)
  
imports: import math
  from typing import Sequence
  import numpy
  import pandas
  import shapely
  from shapely.geometry import Point, Polygon, GeometryCollection
  import geopandas
  import geopandas._compat
  from geopandas import GeoDataFrame, GeoSeries, read_file, sjoin, sjoin_nearest
  from geopandas.testing import assert_geodataframe_equal, assert_geoseries_equal
  from pandas.testing import assert_frame_equal, assert_series_equal
  import pytest
label_desc: Use image node_id 13 for calling the TestSpatialJoinNYBB obj's underlying member method code with example usage: obj.test_sjoin_empty_geometries() without return types