히트 맵을 형성하기 위해 여러 매끈한 다각형을 쌓기

Question

겹치는 매끄러운 다각형이 있습니다. 각 모양은 야생에서 개체의 특정 관찰을 나타냅니다. 폴리곤을 함께 결합하여 일종의 누적 관측 (히트 맵?)을 만들고 싶습니다. 노동 조합 만이 아닙니다 : 나는 그것들을 함께 임계 값을 설정하고 실제로 대상이 어디 있는지 더 잘 예측할 수있는 방식으로 결합하고 싶습니다. 매끈한 다각형을 "래스터 화"하는 것이 최선의 방법은 무엇입니까?

Answer 1

하나가 아마도 한 번에 하나의 다각형을 추가하여 진행할 수 많은 폴리곤이 기여하는 방법들을 각각 분리 된 형태의 목록을 유지하고 기억

import copy 
from itertools import groupby 
from random import randint, seed 
from shapely.geometry import Polygon, box 
from shapely.ops import cascaded_union 

polygons = [ 
    Polygon([(0, 0), (2, 0), (2, 2), (0, 2), (0, 0)]), 
    Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]), 
    Polygon([(1, 1), (2, 1), (2, 2), (1, 2), (1, 1)]) 
] 

def check_shape(s): 
    return (s.geom_type in ['Polygon', 'MultiPolygon'] and not s.is_empty) 

shapes = [] 
for p in polygons: 

    polygon = copy.deepcopy(p) 
    new_shapes = [] 
    for shape_cnt, shape in shapes: 

     new_shapes.extend([ 
      (shape_cnt, shape.difference(polygon)), 
      (shape_cnt+1, shape.intersection(polygon)) 
     ]) 

     polygon = polygon.difference(shape) 
    new_shapes.append((1, polygon)) 

    shapes = list(filter(lambda s: check_shape(s[1]), new_shapes)) 

for p in polygons: 
    print(p) 

for cnt, g in groupby(sorted(shapes, key = lambda s: s[0]), key = lambda s: s[0]): 
    print(cnt, cascaded_union(list(map(lambda s: s[1], g)))) 

를이는 생산 :

POLYGON ((0 0, 2 0, 2 2, 0 2, 0 0)) 
POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)) 
POLYGON ((1 1, 2 1, 2 2, 1 2, 1 1)) 
1 MULTIPOLYGON (((2 1, 2 0, 1 0, 1 1, 2 1)), ((0 1, 0 2, 1 2, 1 1, 0 1))) 
2 MULTIPOLYGON (((1 0, 0 0, 0 1, 1 1, 1 0)), ((1 2, 2 2, 2 1, 1 1, 1 2)))