Python地学分析 — 创建矢量数据点、线、面 04

欢迎关注博主的微信公众号：“智能遥感”。

该公众号将为您奉上Python地学分析、爬虫、数据分析、Web开发、机器学习、深度学习等热门源代码。

Python的小伙伴们，咱们又见面啦！

本节给大家带来得知识点是通过OGR创建点、线、面。

本人的GitHub代码资料主页（持续更新中，多给Star，多Fork）：

https://github.com/xbr2017

~~~~~~~~~~~~~~~~~~~~~~~~~~

编程环境：

操作系统：windows

Python版本：2.7

IDE版本：PyCharm 2018.2.4专业版

~~~~~~~~~~~~~~~~~~~~~~~~~~

创建点的代码：

# _*_ coding: utf-8 _*_
__author__ = 'xbr'
__date__ = '2018/11/4 23:17'

from osgeo import ogr
import matplotlib.pyplot as plt

from ospybook.vectorplotter import VectorPlotter

point = ogr.Geometry(ogr.wkbPoint)  # 构建几何类型:点
point.AddPoint(59.5, 11.5)          # 创建点01
x, y = point.GetX(), point.GetY()   # Python的任性赋值方式
# 调用VectorPlotter类
vp = VectorPlotter(True)
vp.plot(point, 'bo')      # 画出蓝色圆点
point.AddPoint(59.5, 13)  # 在点01基础上添加点02
vp.plot(point, 'rs')      # 画出红色方点

plt.show()   # 少了这句话则图像不显示

生成的结果图：

创建线的代码：

# _*_ coding: utf-8 _*_
__author__ = 'xbr'
__date__ = '2018/11/4 23:25'

from osgeo import ogr
import matplotlib.pyplot as plt

from ospybook.vectorplotter import VectorPlotter

line = ogr.Geometry(ogr.wkbLineString)  # 构建几何类型:线
line.AddPoint(54, 37)      # 添加点01
line.AddPoint(62, 35.5)    # 添加点02
line.AddPoint(70.5, 38)    # 添加点03
line.AddPoint(74.5, 41.5)  # 添加点04
# 调用VectorPlotter类
vp = VectorPlotter(True)
vp.plot(line, 'r-')

plt.show()   # 少了这句话则图像不显示

生成的结果图：

创建多边形的代码：

# _*_ coding: utf-8 _*_
__author__ = 'xbr'
__date__ = '2018/11/4 23:33'

from osgeo import ogr
import matplotlib.pyplot as plt

from ospybook.vectorplotter import VectorPlotter

ring = ogr.Geometry(ogr.wkbLinearRing)  # 构建几何类型:线
ring.AddPoint(58, 38.5)  # 添加点01
ring.AddPoint(53, 6)     # 添加点02
ring.AddPoint(99.5, 19)  # 添加点03
ring.AddPoint(73, 42)    # 添加点04
yard = ogr.Geometry(ogr.wkbPolygon)  # 构建几何类型:多边形
yard.AddGeometry(ring)
yard.CloseRings()
# 调用VectorPlotter类
vp = VectorPlotter(True)
vp.plot(yard, fill=False, edgecolor='blue')
ring = yard.GetGeometryRef(0)
for i in range(ring.GetPointCount()):
    ring.SetPoint(i, ring.GetX(i) - 5, ring.GetY(i))
vp.plot(yard, fill=False, ec='red', linestyle='dashed')
plt.show()   # 少了这句话则图像不显示

生成的结果：

注：需要解释一下的就是面与多边形的关系，本节中第三个代码生成的是多边形，但是本节标题中提到的是面，其实面就是填充后的多边形 ^_^。

最后，希望能够推荐给您的好友，让更多人来学习Python，感受Python的魔力与强大。