常用的一些地理数据处理(geoprocessing)计算函数:
多边形(Polygon):
- 交:
poly1.Intersection(poly2)
- 并:
poly1.Union(poly2)
- 差:
poly1.Difference(poly2)
- 补:
poly1.SymmetricDifference(poly2)
几何形状(Geometry):
- 对几何要素进行缓冲:
<geom>.Buffer(<distance>)
; - 两个geometry是否相等:
<geom1>.Equal(<geom2>)
; - 两个geometry之间的最短距离:
<geom1>.Distance(<geom2>)
; - 用方框边界框住这个几何形状,并返回四个角的点坐标(minx,maxx, miny, maxy):
<geom>.GetEnvelope()
。
空间分析的方法
建议采用GetArea()
获得整个多边形覆盖的面积。
geom.GetArea()`
不过此用法以“度”为量纲,已经失去它作为面积的意义了。若要较为准确的面积,可根据空间参考将度转换为米,再将面积计算为平方米。
判断两个对象的关系
Intersect
语句是判断两个要素是否相交:
poly2.Intersect(poly1)
返回0表示不相交,返回1表示相交。
Disjoint
是判断两个要素是否不相交:
poly2.Disjoint(poly1)
返回1表示不相交,返回0表示相交,跟Intersect
正好相反。
Touch
表示相邻:
poly2.Touches(poly1)
返回0表示不相邻,返回1表示相邻。
Crosses
表示横穿,一般是指一条线穿过一个多边形:
poly2.Crosses(line)
返回0表示不穿过,返回1表示穿过。
Within
表示包含,指一个要素完全被另一个要素包含:
ptB.Within(poly1)
返回0表示点在多边形外,返回1表示点在多边形内。
Contains
也表示包含,跟Within
正好相反:
poly1.Contains(ptB)
与Within
的主调对象与参数对换。
Overlaps
重叠,只有两个多边形之间才能算作overlap
poly2.Overlaps(poly3)
返回0表示不重叠,返回1表示重叠。