使用Python和QgsFeatureRequest()方法在QGIS中进行空间查询

QgsFeatureRequest()方法是在QGIS中进行空间查询的一种常用方法。它可以用来请求满足特定条件的要素。以下是使用Python和QgsFeatureRequest()方法在QGIS中进行空间查询的示例：

首先，需要导入所需的模块和类：

from qgis.core import QgsProject, QgsPointXY, QgsRectangle, QgsFeatureRequest

接下来，加载QGIS项目文件：

project = QgsProject.instance()
project.read('path_to_project_file.qgs')

然后，获取要查询的图层：

layer = project.mapLayersByName('layer_name')[0]

进行空间查询之前，可以选择性地设置查询条件。以下是一些常用的查询条件示例：

1. 矩形查询：

rectangle = QgsRectangle(xMin, yMin, xMax, yMax)
request = QgsFeatureRequest().setFilterRect(rectangle)

2. 距离查询：

point = QgsPointXY(x, y)
request = QgsFeatureRequest().setFilterExpression("distance($geometry, geomFromWkt('POINT(x y)')) < distance_threshold")

3. 属性查询：

request = QgsFeatureRequest().setFilterExpression("attribute_name = 'value'")

通过QgsFeatureRequest()方法设置的各种条件可以组合使用。

接下来，可以使用QgsFeatureRequest()方法执行空间查询并获取结果：

features = layer.getFeatures(request)
for feature in features:
    # 处理查询结果，例如打印要素属性值
    print(feature.attribute('attribute_name'))

最后，需要记得在完成查询后清理资源：

del request

完整的示例代码如下：

from qgis.core import QgsProject, QgsPointXY, QgsRectangle, QgsFeatureRequest

project = QgsProject.instance()
project.read('path_to_project_file.qgs')

layer = project.mapLayersByName('layer_name')[0]

# 矩形查询示例
rectangle = QgsRectangle(xMin, yMin, xMax, yMax)
request = QgsFeatureRequest().setFilterRect(rectangle)

# 距离查询示例
point = QgsPointXY(x, y)
request = QgsFeatureRequest().setFilterExpression("distance($geometry, geomFromWkt('POINT(x y)')) < distance_threshold")

# 属性查询示例
request = QgsFeatureRequest().setFilterExpression("attribute_name = 'value'")

features = layer.getFeatures(request)
for feature in features:
    print(feature.attribute('attribute_name'))

del request

使用Python和QgsFeatureRequest()方法在QGIS中进行空间查询非常简单。只需要设置相应的查询条件，然后执行查询并处理查询结果即可。