docs: def pysal_Join_Counts(self, **kwargs):
          """
          Compute join count statistics for GeoRaster
  
          Usage:
          geo.pysal_Join_Counts(permutations = 1000, rook=True)
  
          arguments passed to raster_weights() and pysal.Join_Counts
          See help(gr.raster_weights), help(pysal.Join_Counts) for options
          """
          if self.weights is None:
              self.raster_weights(**kwargs)
          rasterf = self.raster.flatten()
          rasterf = rasterf[rasterf.mask==False]
          self.Join_Counts = pysal.Join_Counts(rasterf, self.weights, **kwargs)