H3 六边形呈现交换纬度，长 kepler.gl

Question

我想绘制 H3 六边形。来自奥地利。

下载解压 https://biogeo.ucdavis.edu/data/gadm3.6/gpkg/gadm36_AUT_gpkg.zip

完整代码可在 https://gist.github.com/geoHeil/b5b74887e20e4b659d4bb693a700a402 生成生成六边形，如：

size = 7
hexagons = pd.DataFrame(h3.polyfill(geoJson, size), columns=['hexagons'])
hexagons.head()

8752e5b80ffffff
8752ee6c1ffffff

注意 h3 期望 epsg:4326 并稍后再次生成相同的投影 (https://github.com/uber/h3/issues/121)

这给出了一个类似于以下内容的文件：

现在，当移动到 https://kepler.gl/ 并上传数据时，我看到发生了三件奇怪的事情

1. 来自 WKT 线串的多边形被扭曲。这表明使用了错误的投影。但是尝试转换为受支持的 https://github.com/keplergl/kepler.gl/blob/6b380ac6db94e10fed0a76f5e78ef7e55406df21/docs/user-guides/b-kepler-gl-workflow/a-add-data-to-the-map.md Webmercator 并没有修复它

2. 当手动添加六边形层时，它在也门渲染（基于 H3 地址。这看起来很奇怪。这可能是开普勒演示中的错误吗？ 。这看起来真的很奇怪，因为几何图形是使用以下六边形生成的：h3_to_geo_boundary

3. 六边形质心未填充。现在，当使用 h3_to_geo 转换为六边形质心并将数据添加回 ha HexBin 层时，并不是所有的六边形都被填充。但这很奇怪，因为最初所有六边形都可用（参见 1 和 2）。

notice how in (3) the hexbin hexagons are projected correctly as hexagons and not distorted.

Answer 1

我认为这里发生了一些事情：

• 假设您正在使用 h3-py 的 master 分支，signature of polyfill 是 polyfill(geo_json, res, geo_json_conformant=False)。您需要将 geo_json_conformant=True 添加到 polyfill 调用中，否则多边形中的坐标将被解释为 lat,lng 而不是 lng,lat。这可能是您问题的根源。

• 我不是开普勒专家，但我相信 HexBin 层使用生成的笛卡尔 north/south 对齐的六角网格，这就是为什么它们看起来 "correct" 在屏幕上。 H3 六边形具有低失真，但它们确实有一些形状和面积失真，并且它们永远不会 north/south 对齐。当您使用墨卡托投影显示它们时，就像在开普勒中一样，作为投影的函数，它们将有更多的变形，尤其是朝向两极。然而，这里的主要失真问题可能是由于切换 lat,lng - h3_to_multi_polygon 函数 也 需要一个额外的布尔参数来输出符合 GeoJSON 的坐标。

• 我相信 Kepler 也支持 H3 六边形层，所以一种选择是将原始点输入 Kepler，让 Kepler 对 H3 索引进行聚合。

Answer 2

kepler 在所有六边形上使用实例渲染，它假定所有 h3 六边形彼此相对靠近。它使用您当前的地图中心来计算六边形变形并将其应用于所有六边形（实例渲染）。不完美但显着提高了性能。因为计算每个六边形的变形成本太高了。