空间计算

地理数据处理(geoprocessing)计算函数:

多边形(Polygon):

交:poly3.Intersection(poly2)

并:poly3.Union(poly2)

差:poly3.Difference(poly2)

补:poly3.SymmetricDifference(poly2)

几何形状(Geometry):

给geometry加buffer,效果就是把点、线变成多边形,线条加粗: <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表示重叠。