arcpy Thiessen&arcgis泰森多边形

Thiessen和Voronoi是一个意思。

https://blog.csdn.net/elongn/article/details/102689955

arcengine voronoi

https://max.book118.com/html/2018/0827/8076132105001121.shtm

arcgis泰森多边形

https://jingyan.baidu.com/article/09ea3ede2a8b8781aede39ae.html

前提：需要先根据雨量站生成泰森多边形。并且把小流域根据attribute进行split。接下来：

1. Clip

 

各个一维小流域的泰森多边形

2. 添加area属性。（批量）

3. 属性计算器（计算面积）

ArcPy.Calculator_Area

https://blog.csdn.net/gognzixiaobai666/article/details/114065829

4. 添加Percent属性。

5. 计算百分比

 

import arcpy
import numpy
from arcpy import da
def GetFieldValueList(inTable,inField):
     value_list=[]
     rows=arcpy.da.SearchCursor(inTable,inField)
     for row in rows:
         value_list.append(row[0])
     del row
     del rows
     return value_list
def main(fc,c):
    if c=="Sum":
        valueList=GetFieldValueList(fc,"Area")
        return sum(valueList)

!Area! /main("623023AFD1600603_clip","Sum")

ERROR 00539: SyntaxEror EOL while scanning string literal (<expression>, line 1) ：放到记事本里重新复制

[若报错打不开文件，需要先在arcmap中加载xxx_clip.shp系列文件]

6. 导出Excel

7. 对Excel进行编辑，删除不必要的名字部分，删除不必要的属性

只留下STCD和Percent即可。