空间分析：2-4.Python生成泰森多边形

最新推荐文章于 2025-05-24 21:10:33 发布

才华横溢吴道简

最新推荐文章于 2025-05-24 21:10:33 发布

阅读量4.3k

点赞数 1

CC 4.0 BY-SA版权

分类专栏：空间分析

本文链接：https://blog.csdn.net/sinat_41310868/article/details/109611006

这篇博客详细介绍了如何利用Python3的scipy、matplotlib、numpy和shapely库生成泰森多边形和德洛内三角网。代码注释清晰，读者可以通过阅读代码直接理解实现过程。示例中使用matplotlib展示图形，并提到可以将生成的数据导入QGIS进行进一步处理。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

代码注释很详细，如何生成泰森多边形和德洛内三角网。

用的是Python3，使用了scipy和matplotlib、numpy、shapely，引用matplotlib主要是为了看一下效果。

注释很详细，看代码就行。

from scipy.spatial import Voronoi,voronoi_plot_2d
import matplotlib.pyplot as plt
import numpy as np
from shapely.geometry import Polygon, LineString, Point, MultiPoint
from shapely.ops import polygonize
# 要生成德劳内三角网和泰森多边形的点文件
f = open(r'point.txt','r',encoding='utf-8')
flines = f.readlines()
# 所有的点
points = []
for line in flines:
    # line是这样的：1 Point (116.32600952681723072 39.87965125233392882)
    linelist = line.strip('\n').split('\t')
    pointstr = linelist[1].replace('Point (','').replace(')','').split(' ')
    points.append([float(pointstr[0]),float(pointstr[1])])
# 记录点个数，德劳内三角形用
pointlength = len(points)
# 取凸包，确保所有点都能获取一个多边形
# radius单位为度，0.0