创建点的代码：
# _*_ 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() # 少了这句话则图像不显示