基于Python实现全局和局部双变量Moran指数计算

#calculate Local Bivariate Morans‘s I
import osarcpy
from arcpy import env
import numpy as np
from dbfread import DBF
import pandas as pd

env.overwriteOutput = True

def LBMI（shpNameqryfieldXfieldYkthresholdLBMIfield）:
    ‘‘‘qry is criteria to select data
        fieldX and fieldY are variables to calcualte
        k is the k-nearest neighbors
        threshold: the weight is zero if distance over threshold
        LBMIfield is local bivariate moran‘s I
        ‘‘‘
    #add idVariable
    if ‘ID‘ not in [n.name for n in arcpy.ListFields（shpName）]:
        arcpy.AddField_management（shpName“ID““LONG“）
    arcpy.CalculateField_management（shpName“ID“‘!FID!+1‘“PYTHON_9.3“）
    #select and copy data
    temp = ‘temp.shp‘
    arcpy.Copy_management（shpName“temp