def pysal_Join_Counts(self, **kwargs):
if self.weights is None:<EOL><INDENT>self.raster_weights(**kwargs)<EOL><DEDENT>rasterf = self.raster.flatten()<EOL>rasterf = rasterf[rasterf.mask==False]<EOL>self.Join_Counts = pysal.Join_Counts(rasterf, self.weights, **kwargs)<EOL>
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