code: def test_sjoin_predicate(self):
      # points within polygons
      df = sjoin(
          self.pointdf, self.polydf, how="left", predicate="within"
      )
      assert df.shape == (21, 8)
      assert df.loc[1]["BoroName"] == "Staten Island"
  
      # points contain polygons? never happens so we should have nulls
      df = sjoin(
          self.pointdf, self.polydf, how="left", predicate="contains"
      )
      assert df.shape == (21, 8)
      assert np.isnan(df.loc[1]["Shape_Area"])
  
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 6 for calling the TestSpatialJoinNYBB obj's underlying member method code with example usage: obj.test_sjoin_predicate() without return types