def input_permutations(poly, res=5): g = [poly] g = chain_toggle_map(drop_last, g) g = chain_toggle_map(reverse, g) for p in g: hexes = h3.polyfill_polygon(p[0], res=res, holes=p[1:]) yield hexes
def test_polyfill_polygon_holes(): outer, hole1, hole2 = get_us_box_coords() assert 7063 == len( h3.polyfill_polygon(outer, 5) ) for res in 1, 2, 3, 4, 5: hexes_all = h3.polyfill_polygon(outer, res) hexes_holes = h3.polyfill_polygon(outer, res, [hole1, hole2]) hexes_1 = h3.polyfill_polygon(hole1, res) hexes_2 = h3.polyfill_polygon(hole2, res) assert len(hexes_all) == len(hexes_holes) + len(hexes_1) + len(hexes_2) assert hexes_all == set.union(hexes_holes, hexes_1, hexes_2)
def test_h3_set_to_multi_polygon(): h = '8928308280fffff' hexes = h3.k_ring(h, 1) mpoly = h3.h3_set_to_multi_polygon(hexes) out = h3.polyfill_polygon(mpoly[0][0], 9, holes=None, lnglat_order=False) assert out == hexes
def test_2_polys(): h = '8928308280fffff' hexes = h3.hex_ring(h, 2) hexes = hexes | {h} # hexes should be a center hex, and the 2-ring around it # (with the 1-ring being absent) out = [ h3.polyfill_polygon(poly[0], 9, holes=poly[1:], lnglat_order=False) for poly in h3.h3_set_to_multi_polygon(hexes, geo_json=False) ] assert set.union(*out) == hexes
def test_polyfill_polygon(): # lat/lngs for State of Maine maine = [ (45.137451890638886, -67.13734351262877), (44.8097, -66.96466), (44.3252, -68.03252), (43.98, -69.06), (43.68405, -70.11617), (43.090083319667144, -70.64573401557249), (43.08003225358635, -70.75102474636725), (43.21973948828747, -70.79761105007827), (43.36789581966826, -70.98176001655037), (43.46633942318431, -70.94416541205806), (45.3052400000002, -71.08482), (45.46022288673396, -70.6600225491012), (45.914794623389355, -70.30495378282376), (46.69317088478567, -70.00014034695016), (47.44777598732787, -69.23708614772835), (47.184794623394396, -68.90478084987546), (47.35462921812177, -68.23430497910454), (47.066248887716995, -67.79035274928509), (45.702585354182816, -67.79141211614706), (45.137451890638886, -67.13734351262877) ] # a very rough hexagonal approximation to the State of Maine expected = { '832b13fffffffff', '832b18fffffffff', '832b1afffffffff', '832b1efffffffff', '832ba9fffffffff', '832badfffffffff' } out = h3.polyfill_polygon(maine, 3) assert out == expected
def test_compact(): # lat/lngs for State of Maine maine = [ (45.137451890638886, -67.13734351262877), (44.8097, -66.96466), (44.3252, -68.03252), (43.98, -69.06), (43.68405, -70.11617), (43.090083319667144, -70.64573401557249), (43.08003225358635, -70.75102474636725), (43.21973948828747, -70.79761105007827), (43.36789581966826, -70.98176001655037), (43.46633942318431, -70.94416541205806), (45.3052400000002, -71.08482), (45.46022288673396, -70.6600225491012), (45.914794623389355, -70.30495378282376), (46.69317088478567, -70.00014034695016), (47.44777598732787, -69.23708614772835), (47.184794623394396, -68.90478084987546), (47.35462921812177, -68.23430497910454), (47.066248887716995, -67.79035274928509), (45.702585354182816, -67.79141211614706), (45.137451890638886, -67.13734351262877) ] res = 5 h_uncomp = h3.polyfill_polygon(maine, res) h_comp = h3.compact(h_uncomp) expected = {'852b114ffffffff', '852b189bfffffff', '852b1163fffffff', '842ba9bffffffff', '842bad3ffffffff', '852ba9cffffffff', '842badbffffffff', '852b1e8bfffffff', '852a346ffffffff', '842b1e3ffffffff', '852b116ffffffff', '842b185ffffffff', '852b1bdbfffffff', '852bad47fffffff', '852ba9c3fffffff', '852b106bfffffff', '852a30d3fffffff', '842b1edffffffff', '852b12a7fffffff', '852b1027fffffff', '842baddffffffff', '852a349bfffffff', '852b1227fffffff', '852a3473fffffff', '852b117bfffffff', '842ba99ffffffff', '852a341bfffffff', '852ba9d3fffffff', '852b1067fffffff', '852a3463fffffff', '852baca7fffffff', '852b116bfffffff', '852b1c6bfffffff', '852a3493fffffff', '852ba9dbfffffff', '852b180bfffffff', '842bad7ffffffff', '852b1063fffffff', '842ba93ffffffff', '852a3693fffffff', '852ba977fffffff', '852b1e9bfffffff', '852bad53fffffff', '852b100ffffffff', '852b102bfffffff', '852a3413fffffff', '852ba8b7fffffff', '852bad43fffffff', '852b1c6ffffffff', '852a340bfffffff', '852b103bfffffff', '852b1813fffffff', '852b12affffffff', '842a34dffffffff', '852b1873fffffff', '852b106ffffffff', '852b115bfffffff', '852baca3fffffff', '852b114bfffffff', '852b1143fffffff', '852a348bfffffff', '852a30d7fffffff', '852b181bfffffff', '842a345ffffffff', '852b1e8ffffffff', '852b1883fffffff', '852b1147fffffff', '852a3483fffffff', '852b12a3fffffff', '852a346bfffffff', '852ba9d7fffffff', '842b18dffffffff', '852b188bfffffff', '852a36a7fffffff', '852bacb3fffffff', '852b187bfffffff', '852bacb7fffffff', '842b1ebffffffff', '842b1e5ffffffff', '852ba8a7fffffff', '842bad9ffffffff', '852a36b7fffffff', '852a347bfffffff', '832b13fffffffff', '852ba9c7fffffff', '832b1afffffffff', '842ba91ffffffff', '852bad57fffffff', '852ba8affffffff', '852b1803fffffff', '842b1e7ffffffff', '852bad4ffffffff', '852b102ffffffff', '852b1077fffffff', '852b1237fffffff', '852b1153fffffff', '852a3697fffffff', '852a36b3fffffff', '842bad1ffffffff', '842b1e1ffffffff', '852b186bfffffff', '852b1023fffffff'} # noqa assert h_comp == expected return h_uncomp, h_comp, res
def test_polyfill_polygon_order(): lnglat, _, _ = get_us_box_coords(order='lnglat') out = h3.polyfill_polygon(lnglat, 5, lnglat_order=True) assert len(out) == 7063