手把手教你用osmium-tool和tilemaker从osm.pbf提取中国铁路网数据并生成mbtiles

从OSM数据构建中国铁路地图的完整技术指南

地理信息系统开发者经常需要处理特定领域的地图数据，而铁路网络作为基础设施的核心组成部分，其数据提取与可视化具有广泛的应用价值。本文将详细介绍如何从OpenStreetMap的原始数据中精确提取中国铁路网络，并转换为轻量级的mbtiles格式，适用于各类WebGIS应用和移动端地图服务。

1. 环境准备与数据获取

在开始处理地理空间数据前，需要搭建合适的工作环境。推荐使用Ubuntu 20.04或更高版本作为基础系统，因其对地理信息工具链的支持最为完善。

必备工具安装：

sudo apt update sudo apt install -y git cmake make g++ libboost-dev libexpat1-dev zlib1g-dev \ libbz2-dev libsparsehash-dev libprotobuf-dev protobuf-compiler \ libshp-dev libsqlite3-dev

对于原始数据源，我们有两个主要选择：

• 全球完整数据集：planet-latest.osm.pbf（约60GB）
• 亚洲区域数据集：asia-latest.osm.pbf（约8GB）

下载命令示例：

wget https://download.geofabrik.de/asia-latest.osm.pbf

提示：如果只需要中国区域数据，建议先下载亚洲数据集以减少处理时间，后续再通过边界提取进一步缩小范围。

2. 精确提取中国区域数据

中国国境线复杂，需要精确的地理围栏定义。我们使用GeoJSON格式来描述提取范围：

{ "type": "FeatureCollection", "features": [{ "type": "Feature", "properties": {}, "geometry": { "type": "Polygon", "coordinates": [[ [73.08, 39.33], [135.35, 48.42], [134.99, 48.64], [125.53, 53.36], [119.48, 53.03], [87.90, 49.25], [79.63, 45.10], [73.08, 39.33] ]] } }] }

使用osmium-tool进行区域提取：

osmium extract -p china.geojson asia-latest.osm.pbf -o china.osm.pbf --overwrite

参数说明：

• -p：指定GeoJSON边界文件
• --overwrite：覆盖已存在的输出文件

3. 铁路网络数据过滤

OpenStreetMap中的铁路数据主要包含两类要素：

1. 铁路线路（railway=rail）
2. 火车站节点（railway=station）

执行精细化过滤：

osmium tags-filter china.osm.pbf \ wr/railway=rail \ wr/railway=station \ wr/railway=halt \ wr/railway=tram_stop \ wr/railway=subway_entrance \ -o china-railways.osm.pbf

关键标签说明：

4. 使用TileMaker生成MBTiles

MBTiles是一种基于SQLite的地图切片存储格式，特别适合Web地图应用。我们需要准备TileMaker的配置文件和处理脚本。

编译安装TileMaker：

git clone https://github.com/systemed/tilemaker.git cd tilemaker mkdir build cd build cmake .. make sudo make install

配置文件示例（config.json）：

{ "layers": { "railways": { "minzoom": 5, "maxzoom": 14, "simplify_below": 10, "simplify_level": 2.5 }, "stations": { "minzoom": 10, "maxzoom": 14 } } }

处理脚本关键部分（process.lua）：

function process_node(node) if node.tags.railway == "station" then way:layer("stations") way:attribute("name", node.tags.name) end end function process_way(way) if way.tags.railway == "rail" then way:layer("railways") way:attribute("gauge", way.tags.gauge or "standard") end end

生成最终MBTiles文件：

tilemaker --input china-railways.osm.pbf \ --output china-railways.mbtiles \ --config config.json \ --process process.lua \ --bbox 73.5,18.2,135.0,53.5

5. 性能优化与质量控制

处理大规模地理数据时，效率和质量控制至关重要。以下是几个关键优化点：

内存管理技巧：

• 使用--fsync参数降低I/O负载
• 设置--shard-stores进行并行处理
• 限制中间文件大小--temp-dir=/dev/shm

数据质量检查：

osmium fileinfo -e china-railways.osm.pbf

典型输出分析：

Header: Bounding boxes: (73.5,18.2,135.0,53.5) With history: no Options: generator=osmium/1.14.0 osmosis_replication_timestamp=2023-11-01T00:00:00Z Features: nodes: 1245678 ways: 234567 relations: 1234

常见问题解决方案：

问题：处理过程中内存不足 解决：添加交换空间或使用--compact模式

问题：生成的MBTiles文件过大 解决：调整--min-zoom和--max-zoom参数

问题：铁路线显示不连续 解决：检查原始OSM数据质量，可能需要手动修复

在实际项目中，我曾遇到处理中国西部铁路数据时节点密度不足的问题。通过调整osmium的--buffer参数和tilemaker的--simplify-level参数，最终获得了理想的视觉效果和性能平衡。