要检查一个新对象是否与所有其他图层的几何体相交并输出所属图层，建议采用以下两种方案。由于 QGIS 的单个相交算法通常仅支持两个图层之间的运算，因此需要先进行数据整合或使用 PyQGIS 代码。
方案一：使用处理工具
这种方法通过将所有参考图层合并为一个临时图层，然后进行相交分析，从而在结果中保留原图层的属性信息。
第一步：合并所有参考图层
- 算法名称：合并矢量图层
- 算法 ID：native:mergevectorlayers
- 操作：选择所有需要检查的"其他图层"。合并后，QGIS 通常会生成一个包含原图层来源信息的字段（或者你可以预先为每个图层添加一个标识字段），这样在后续步骤中就能知道相交的对象来自哪个图层。
第二步：执行相交分析
- 算法名称：相交
- 算法 ID：native:intersection
- 操作：
- 输入图层：你的新对象图层。
- 叠加图层：上一步生成的"合并后的图层"。
- 结果：输出图层将仅包含相交的部分。通过查看结果图层的属性表，你可以直接看到该相交对象所属的原始图层名称。
方案二：使用 PyQGIS 实现
下面是一段可直接在 QGIS Python 控制台或作为独立脚本运行的 PyQGIS 代码。它会创建一个示例几何体作为"新对象"，然后遍历项目中所有矢量图层，检查该新对象是否与图层中的要素相交，并打印出相交的图层名称及要素 ID。
from qgis.core import (
QgsProject, QgsVectorLayer, QgsGeometry, QgsFeatureRequest,
QgsWkbTypes
)
from qgis.PyQt.QtCore import Qt, QVariant
from qgis.PyQt.QtGui import QColor, QFont, QBrush
# 获取当前项目
project = QgsProject.instance()
# ============================================================
# 1. 定义"新对象"的几何体
# ============================================================
# 示例：使用 WKT 创建一个多边形。
# 你可以将其替换为实际对象，例如：
# new_geom = feature.geometry() # 来自某个要素
new_geom = QgsGeometry.fromWkt(
"POLYGON((105.15 32.70, 105.25 32.70, 105.25 32.80, 105.15 32.80, 105.15 32.70))"
)
if new_geom.isNull() or new_geom.isEmpty():
print("新对象几何体无效，请检查输入。")
else:
print("开始检查新对象与项目中其他矢量图层的相交情况...")
print("-" * 60)
# 获取新对象的边界框，用于快速空间过滤
bbox = new_geom.boundingBox()
# ============================================================
# 2. 遍历项目中的所有地图图层
# ============================================================
for layer_id, layer in project.mapLayers().items():
# 仅处理矢量图层
if not isinstance(layer, QgsVectorLayer):
continue
# 跳过无效图层
if not layer.isValid():
print(f"图层 '{layer.name()}' 无效，已跳过。")
continue
# 跳过无几何信息的图层
if layer.geometryType() == QgsWkbTypes.NullGeometry:
print(f"图层 '{layer.name()}' 无几何信息，已跳过。")
continue
# ============================================================
# 3. 使用边界框快速过滤，再精确判断相交
# ============================================================
request = QgsFeatureRequest()
request.setFilterRect(bbox)
intersected_ids = []
for feature in layer.getFeatures(request):
feature_geom = feature.geometry()
if feature_geom.isNull() or feature_geom.isEmpty():
continue
# 精确相交检查
if new_geom.intersects(feature_geom):
intersected_ids.append(feature.id())
# ============================================================
# 4. 输出结果
# ============================================================
if intersected_ids:
print(f"图层 '{layer.name()}' 与新对象相交，相交要素 ID: {intersected_ids}")
else:
print(f"图层 '{layer.name()}' 与新对象不相交。")
print("-" * 60)
print("检查完成。")
使用说明
替换新对象：将 new_geom = QgsGeometry.fromWkt(...) 一行替换为你实际要检查的对象。例如：
source_layer = QgsProject.instance().mapLayersByName("你的图层名")[0]
new_geom = next(source_layer.getFeatures()).geometry()
更多 QGIS 使用问题，欢迎留言，转载须注明出处。