OGR几何关系与操作
关系
本文脚本需要导入ogr模块
from osgeo import ogr
Equals
两个几何的边界、内部和外部重合
a1=ogr.CreateGeometryFromWkt("POINT(1 1)")
a2=ogr.CreateGeometryFromWkt("POINT(1 1)")
print(a1.Equals(a2)) # True
Contains
几何包含另一个几何的内部和边界,并且边界不接触,适用于所有几何类型,并且a.Contains(b)==b.Within(a)
p=ogr.CreateGeometryFromWkt("POLYGON((0 0,0 2,2 2,2 0,0 0))")
l=ogr.CreateGeometryFromWkt("LINESTRING(0.5 0.5,1.5 1.5)")
a=ogr.CreateGeometryFromWkt("POINT(0.5 0.5)")
b=ogr.CreateGeometryFromWkt("POINT(1 1)")
c=ogr.CreateGeometryFromWkt("POINT(1.5 1.5)")
print(p.Contains(l)) # True