我最近开始处理具有SRTM数据的项目,并已使用phyghtmap
提取了pbf文件。
首先,我要获取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
最后使用phyghtmap --max-nodes-per-tile=0 -s 10 -0 --pbf warp-90.tif
结果是pbf
个文件的列表。当我使用osm2pgsql
将它们加载到PostGIS中时,它们非常好。但我想将它们合并以固定导入。
我尝试了所有主要的解决方案:
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
它们都不起作用,结果只是结果文件中合并数据的一小部分。
我做错什么了吗?
注意:如果我用--append
加载所有pbf,则可以使用,但是在世界上很小的一部分需要花一些时间。
答案 0 :(得分:2)
我发现了问题。我没有在脚本中设置--start-node-id
和--start-way-id
,所以我所有的pbf
都使用相同的ID范围。现在,我分配了唯一的ID,它的工作方式就像一个超级魅力:)