def pysal_Moran(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.Moran = pysal.Moran(rasterf, self.weights, **kwargs)<EOL>
Compute Moran's I measure of global spatial autocorrelation for GeoRaster

Usage:
geo.pysal_Moran(permutations = 1000, rook=True)

arguments passed to raster_weights() and pysal.Moran
See help(gr.raster_weights), help(pysal.Moran) for options