pythongdal教程之：过滤器，简单的空间分析，函数和模块

Layer对象有一个方法叫SetAttributeFilter(<where_clause>)可以将Layer中符合某一条件的Feature过滤出来。设定了Filter之后就可以用GetNextFeature()方法依次取出符合条件的Feature了。SetAttributeFilter(None)可以清楚一个Filter。例如

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)

就是把主调对象和参数换一下啦

Overlaps重叠，好像只有两个多边形之间才能overlap

poly2.Overlaps(poly3)

返回0表示不重叠，返回1表示重叠

下面看看简单的地理数据处理geoprocessing

多边形的：

交：poly3.Intersection(poly2)

并：poly3.Union(poly2)

差：poly3.Difference(poly2)

补：poly3.SymmetricDifference(poly2)

geometry的：

<geom>.Buffer(<distance>) 给geometry加buffer，就是把点线变成多边形，变粗了

<geom1>.Equal(<geom2>) 两个geometry相等吗？

<geom1>.Distance(<geom2>) 返回两个geometry之间的最短距离

<geom>.GetEnvelope() 信封，有意思，其实就是用一个方框框住这个几何形状，返回四个角的坐标(minx, maxx, miny, maxy)

python的函数function，异常exception和模块module可以从任何一本python教材上找到，在此不多赘述

以上就是python gdal教程之：过滤器，简单的空间分析，函数和模块的内容，更多相关内容请关注PHP中文网（HdhCmsTestgxlcms测试数据）！

查看更多关于pythongdal教程之：过滤器，简单的空间分析，函数和模块的详细内容...