1

我最近开始使用 SRTM 数据开展一个项目，并使用phyghtmap.

首先，我正在获取hgt文件，将它们转换为tif使用以下命令：gdal_fillnodata.py data.hgt data.tif

然后我用它们扭曲它们gdalwarp -co BIGTIFF=YES -co TILED=YES -co COMPRESS=LZW -co PREDICTOR=2 -t_srs "+proj=merc +ellps=sphere +R=6378137 +a=6378137 +units=m" -r bilinear -tr 90 90 data.tif warp-90.tif

最后创建 pbf 文件phyghtmap --max-nodes-per-tile=0 -s 10 -0 --pbf warp-90.tif

结果是pbf文件列表。当我将它们加载到 PostGIS 中时，它们非常好osm2pgsql。但我想合并它们以加强导入。

我已经尝试了所有主要的解决方案：

  • osmium merge *.pbf -o merged.pbf

  • 转换pbf为o5m然后osmconvert64 *.o5m -o=merge.o5m再转换回pbf

  • 两两合并osmosis --read-pbf lon4.00_5.00lat44.00_45.00_local-source.pbf --read-pbf lon5.00_6.00lat44.00_45.00_local-source.osm.pbf --merge --write-pbf osmo_merge.osm.pbf

它们都不起作用，结果只是合并到结果文件中的数据的一小部分。

难道我做错了什么？

注意：如果我用它加载所有 pbf，--append它就可以工作，但是对于世界上很小的一部分来说，它需要很长时间。

4

1 回答 1

2

我发现了这个问题。我没有在我的脚本中设置--start-node-idand所以我所有的人都使用相同的 id 范围。现在我正在分配唯一的 ID，它就像一个魅力:)--start-way-idpbf

于 2018-12-17T08:39:36.627 回答