Question

我正在尝试创建一个新字典（intersections），该字典包含名为zones的字典中多边形的相交区域。

我使用combinations查找zones的所有可能的唯一组合。然后，我使用Shapely的.intersects()测试if区域相交。如果他们这样做，我希望将它们的几何图形保存在变量int_geometry中，然后存储在字典中（使用：Shapely的.intersection()）。

我知道有四个交集，因为此代码返回了它们：

for a, b in combinations(zones.values(), 2):
  a_geom = a['location']
  b_geom = b['location']
  if a_geom.intersects(b_geom) == True:
    print a_geom.intersection(b_geom)

但是，如果像下面的代码那样替换if语句之后的内容，它将开始覆盖自身。

intersections = {}    

for a, b in combinations(zones.values(), 2):
  a_geom = a['location']  
  b_geom = b['location']  
  if a_geom.intersects(b_geom) == True:
    int_geometry = a_geom.intersection(b_geom)
    int_area = round(int_geometry.area,2)
    int_perimeter = round(int_geometry.length,2)
    intersections = {
      'geometry' : int_geometry,
      'attributes' : {
        'area' : int_area,
        'perimeter' : int_perimeter,
      }
    }    

pprint(intersections)

关于相似问题，有多个主题，尽管我找不到答案。我知道我在这里忽略了很明显的事情，但是我无法察觉。有人可以向我解释我在做什么错吗？

Answer 1

index中的

int_index = "int_index " + str(index + 1)将始终与第一个循环后的值相同，其值是恒定的。因此，数据将被覆盖。

...

# PART 2 - ANALYSE THE DATA
intersections = {}
index = 0

for a, b in combinations(zones.values(), 2):
  a_geom = a['location']
  b_geom = b['location']
  if a_geom.intersects(b_geom) == True:
    int_index = "int_index " + str(index + 1)
    int_geometry = a_geom.intersection(b_geom)
    int_area = round((a_geom.intersection(b_geom).area),2)
    int_perimeter = round((a_geom.intersection(b_geom).length),2)
    intersections[int_index] = {
      'geometry' : int_geometry,
      'attributes' : {
        'area' : int_area,
        'perimeter' : int_perimeter,
      }
    }
    index += 1

pprint(intersections)

防止Python For Loop覆盖字典

1 个答案: