我是新手，很抱歉，如果这很明显/已经有人问过并回答了，但我找不到任何东西。
我有两个shapefile：1.英国某县的行政边界层，称为LSOA边界，其中有500个小区域2.洪水区。
理想情况下，我想找出洪水泛滥区内的哪个小LSOA区域≥50％，最后对500个LSOA区域中的每一个以是/否或1/0结尾。
但是我不知道该怎么做。我以为可以合并两个shapefile，但是它们之间没有共同的属性。然后，我想我可以使用“按位置连接属性”功能，该功能有效，并向我显示了哪些LSOA处于洪泛区，但是几乎所有这些都存在（参见图2）。
我认为这是一个SQL问题，但我不知道。我是QGIS的新手，从未使用过PostgreSQL。
任何帮助将不胜感激。我可以提供您可爱的人需要帮助的任何信息。
Answers:
使用QGIS中包含的地理处理工具，这是一项相对简单的任务。
计算您的LSOA区域的面积。
将洪泛区图层合并为单个多部分要素。
Vector > Geometry Tools > Singleparts to Multipart。将LSOA区域层与多部分洪水区域层相交。
Vector > Geoprocessing Tools > Intersect。结果层将是LSOA区域的一部分（具有LSOA区域层的属性），该部分与洪泛区层重叠。要计算洪泛区内每个LSOA区域的比例：
使用两层共享的唯一ID，将原始LSOA层连接到相交的层。
将连接的图层导出为新的shapefile。
删除重复的属性。
等等！
如果没有步骤2，将为每个LSOA功能的每个不同的洪泛区功能创建一个单独的功能。如果您仅对每个LSOA区域的总覆盖率感兴趣，那么这可能不是您想要的。如果要区分河流洪水/潮汐洪水/洪灾洪水（且洪水区数据支持该洪水洪水），则可以将“部分”转换为“多部分”，将“ TYPE”字段指定为“唯一ID”字段。
Vector > Research Tools > Vector Grid... Output grid as polygons）具有相同程度的矢量网格，然后将网格与泛洪层相交。然后在第3步中使用输出而不是Flood层。我猜测该层为空是因为它崩溃了。
您可以使用spatialite和一些空间SQL函数。
Select t1.geometry, t1.ID, area(t1.geometry), area(t2.geometry) ...... (anything you need to have in the table results)
(area(intersection(t1.geometry,t2.geometry))) as "Commun_AREA"
, ("Commun_AREA"*100/(area(t1.geometry))) as "Percent_AREA"
From lsoa as t1, flood_zone as t2
Where Intersects( t1.geometry,t2.geometry ) = 1
这似乎比提交的答案要容易得多。我会亲自使用一个简单的python脚本：
floodName = "the layer name here"
boundryName = "the layer name here"
fieldName = "the name of the field to contain the output 1/0"
minCoverage = 0.5 # the minimum amount of area covered to write 1
updateMap = [] # this will store values to be written
# get layers
floodLayer = QgsMapLayerRegistry.instance().mapLayersByName(floodName)[0]
boundryLayer = QgsMapLayerRegistry.instance().mapLayersByName(boundryName)[0]
fieldIndex = boundryLayer.dataProvider().fieldNameIndex(fieldName)
# iterate through boundries
for b in boundryLayer.getFeatures():
# get only flood features that intersect with this feature's bounding box
# this will make the script go way faster than it would otherwise
request = QgsFeatureRequest().setFilterRect(b.geometry().boundingBox())
floodGeom = geometry()
floodFeat = QgsFeature()
iter = floodLayer.getFeatures(request)
iter.nextFeature(feat)
while iter.nextFeature(feat):
floodGeom = floodGeom.combine(feat.geometry())
intersectGeom = b.geometry().intersection(feat.geometry())
if intersectGeom.area() > minCoverage * b.geometry().area():
updateMap[b.id()] = {fieldIndex : 1}
else:
updateMap[b.id()] = {fieldIndex : 0}
boundryLayer.dataProvider().changeAttributeValues(updateMap)
这只会评估与每个边界图层的边界框相交的洪水多边形，因此应该相当快地运行，然后只更新现有图层中的一个字段（而不是制作一个新的新图层并复制旧值的复杂操作然后删除）