code: def geo_contains(left, right):
      
      op = ops.GeoContains(left, right)
      return op.to_expr()
docs: Check if the first geometry contains the second one
  
      Parameters
      ----------
      left : geometry
      right : geometry
  
      Returns
      -------
      contains : bool scalar