#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