arcgispython空间插值_基于python的空间距离权重计算——arcgis中的arcpy

arcgis出了10后&＃xff0c;python代替vb成为官方脚本语言&＃xff0c;python中的gis库也是有很多&＃xff0c;但因为本人最先接触的是arcgis所以就直接学习他自带的库arcpy了&＃xff0c;主要以空间权重矩阵作为一个契机顺便学习arcpy。

代码及数据资料&＃xff1a;https://github.com/shikanon/WeightDistance

因为还是不太习惯arcpy中的表操作&＃xff0c;所以借助了pandas进行表操作&＃xff0c;构建一个表转换&＃xff1a;

#ListFields包含field类的数组

fields&＃61;arcpy.ListFields(path)

def GetTable():

&＃39;&＃39;&＃39;将arcpy表单变为pandas表单&＃xff0c;还是喜欢pandas些~&＃39;&＃39;&＃39;

table&＃61;[]

fieldname&＃61;[field.name for field in fields]

#游标集合&＃xff0c;用for 循环一次后没办法循环第二次!一个游标实例只能循环一次

data&＃61;arcpy.SearchCursor(path)

for row in data:

#Shape字段中的要数是一个几何类

r&＃61;[]

for field in fields:

r.append(row.getValue(field.name))

table.append(r)

return pd.DataFrame(table,columns&＃61;fieldname)

因为求欧几里距离在arcpy 中polyline直接有length属性所以只需要构建一个线类&＃xff0c;然后直接求长度就行了

def getdistance(point1,point2,sparef):

&＃39;&＃39;&＃39;求两点的距离&＃39;&＃39;&＃39;

#两点构造一条线

l&＃61;arcpy.Polyline(arcpy.Array([point1,point2]),sparef)

#求线长

return l.getLength()

在求长度时候需要考虑到shp的投影坐标&＃xff0c;因此需要先得到一个spatialReference

#空间坐标投影参数

#随便抽取一个shape出来求得该图层的空间投影

sparef&＃61;data[&＃39;Shape&＃39;][0].spatialReference

print sparef.name

#计算质心centre

#Shape字段下是一个geometry object对象,其下有中心点和面积等工具

#貌似由centroid抽出来的点没有空间投影属性,

#所以data[&＃39;centre&＃39;][0].spatialReference将会报错

data[&＃39;centre&＃39;]&＃61;[d.centroid for d in data[&＃39;Shape&＃39;]]

#计算欧几里权重矩阵

weight&＃61;[]

for i in data[&＃39;centre&＃39;]:

row&＃61;[]

for j in data[&＃39;centre&＃39;]:

row.append(getdistance(i,j,sparef))

weight.append(row)

weight&＃61;pd.DataFrame(weight,index&＃61;data[&＃39;FID&＃39;],columns&＃61;data[&＃39;FID&＃39;])

最后导出csv:

#输出csv

weight.to_csv(path_or_buf&＃61;&＃39;weight.csv&＃39;,sep&＃61;&＃39;,&＃39;)

小小练习完成&＃xff0c;查看一下结果&＃xff0c;图为广东省市级单元的空间权重矩阵&＃xff1a;

贴上完整代码&＃xff1a;

# -*- coding: cp936 -*-

import arcpy

import numpy as np

import pandas as pd

path&＃61;&＃39;data/gd.shp&＃39;

#ListFields包含field类的数组

fields&＃61;arcpy.ListFields(path)

def GetTable():

&＃39;&＃39;&＃39;将arcpy表单变为pandas表单&＃xff0c;还是喜欢pandas些~&＃39;&＃39;&＃39;

table&＃61;[]

fieldname&＃61;[field.name for field in fields]

#游标集合&＃xff0c;用for 循环一次后没办法循环第二次!一个游标实例只能循环一次

data&＃61;arcpy.SearchCursor(path)

for row in data:

#Shape字段中的要数是一个几何类

r&＃61;[]

for field in fields:

r.append(row.getValue(field.name))

table.append(r)

return pd.DataFrame(table,columns&＃61;fieldname)

def getdistance(point1,point2,sparef):

&＃39;&＃39;&＃39;求两点的距离&＃39;&＃39;&＃39;

#两点构造一条线

l&＃61;arcpy.Polyline(arcpy.Array([point1,point2]),sparef)

#求线长

return l.getLength()

data&＃61;GetTable()

#空间坐标投影参数

#随便抽取一个shape出来求得该图层的空间投影

sparef&＃61;data[&＃39;Shape&＃39;][0].spatialReference

print sparef.name

#计算质心centre

#Shape字段下是一个geometry object对象,其下有中心点和面积等工具

#貌似由centroid抽出来的点没有空间投影属性,

#所以data[&＃39;centre&＃39;][0].spatialReference将会报错

data[&＃39;centre&＃39;]&＃61;[d.centroid for d in data[&＃39;Shape&＃39;]]

#计算欧几里权重矩阵

weight&＃61;[]

for i in data[&＃39;centre&＃39;]:

row&＃61;[]

for j in data[&＃39;centre&＃39;]:

row.append(getdistance(i,j,sparef))

weight.append(row)

weight&＃61;pd.DataFrame(weight,index&＃61;data[&＃39;FID&＃39;],columns&＃61;data[&＃39;FID&＃39;])

#输出csv

weight.to_csv(path_or_buf&＃61;&＃39;weight.csv&＃39;,sep&＃61;&＃39;,&＃39;)

http://my.oschina.net/Kanonpy/

------------窈----------------窈---------------窈------------窈------------窈-------------窈------------窈-------