【Python】ogr使用小记

• 图层的属性筛选

    # 通过属性进行要素筛选
    lyr.SetAttributeFilter(field_name + "=" + str(val))

    # 还原图层值筛选前的状态
    lyr.SetAttributeFilter(None)

• ogr中的geometry与shapely互相转换

test_poly_file = r"F:\Temp\temp_gis_data\testSplittingPolygon.shp"
ds = shp_driver.Open(test_poly_file, 0)
test_lyr = ds.GetLayer(0)
feature = test_lyr[0]
test_poly = loads(feature.geometry().ExportToWkt())

• ogr对geometry使用clone方法，当feature被销毁以后该geometry还能使用

geo = feat.geometry().Clone()

• 使用SQL语句查询地块的平均面积

def getAvgAreaByAttribute(ds, clause):
    lyr_name = ds.GetLayer(0).GetName()
    sql_query = f'''SELECT AVG(OGR_GEOM_AREA) as avg_area
             FROM {lyr_name}
             WHERE {clause}'''
    res_lyr = ds.ExecuteSQL(sql_query)
    avg_area = res_lyr[0].GetField('avg_area')
    ds.ReleaseResultSet(res_lyr)  # 清除查询结果
    return avg_area

• 关于SetSpatialAttribute()这个函数：只要与filter geometry相交的geometries都会被返回。（This function set the geometry to be used as a spatial filter when fetching features via the OGR_L_GetNextFeature() function. Only features that geometrically intersect the filter geometry will be returned.）
• 对于shapefile的datasource可以进行SQL查询，但是不支持CopyLayer方法。
• ogr的Buffer方法要与SetSpatialAttribute结合起来使用

如果你认为对你有用，或者想要与我有更多交流，关注我的微信公众号支持我一下吧！~