QGIS 3.32实战5分钟搞定OpenStreetMap城市路网数据下载与导入在数字化城市研究与规划领域获取精准的路网数据往往成为项目启动的第一道门槛。传统方法需要手动划定经纬度范围不仅效率低下还会引入大量冗余数据。现在借助QGIS 3.32与Overpass API的黄金组合即使是GIS新手也能在咖啡冷却前完成城市路网的精准采集与可视化。1. Overpass API数据获取实战1.1 中文城市ID智能查询OpenStreetMap的Overpass API如同城市数据的搜索引擎但直接使用英文查询中国城市常会遇到匹配偏差。这里有个专业技巧通过name:zh参数锁定中文行政区划。打开Overpass Turbo在线查询工具https://overpass-turbo.eu/在左侧编辑区输入以下代码模板osm-script query typerelation has-kv kboundary vadministrative/ has-kv kname:zh v武汉市/ /query print/ /osm-script执行查询后结果面板会显示类似这样的关键信息relation id3797255 tag kname:zh v武汉市/ tag kboundary vadministrative/ /relation注意国内部分新城区可能出现多个关联ID建议选择admin_level6的市级行政区划关系。若查询超时可添加timeout25参数延长等待时间。1.2 路网数据精准下载获取城市ID后需要将其转换为Overpass特有的area引用值。这个转换公式让很多初学者困惑area_ref 3600000000 relation_id以武汉市为例ID3797255其area_ref就是3603797255。现在使用这个魔法数字构建路网查询语句osm-script timeout1800 union area-query ref3603797255/ recurse typenode-relation intorels/ recurse typenode-way/ recurse typeway-relation/ /union union item/ recurse typeway-node/ /union print modebody/ /osm-script这段代码的精妙之处在于area-query限定数据获取范围三重recurse确保获取完整拓扑关系timeout1800避免大数据量城市超时中断点击运行后通过导出按钮下载OSM格式数据。建议选择map.osm命名方便后续识别。2. QGIS 3.32数据处理技巧2.1 OSM数据智能转换新版QGIS 3.32内置的OSM数据处理器堪称黑科技。无需安装插件三步完成格式转换菜单导航图层 → 数据源管理器 → 矢量 → 协议类型选择OSM XML文件加载点击...按钮选择下载的map.osm文件图层选择勾选lines类型道路数据主要存储于此常见报错解决方案若遇到CRS未定义警告右键图层 → 属性 → 源 → 将CRS强制设置为EPSG:43262.2 路网数据深度优化原始OSM数据包含大量冗余属性需要针对性过滤# 在QGIS Python控制台执行以下脚本优化数据 layer iface.activeLayer() request QgsFeatureRequest() request.setFilterExpression(highway IS NOT NULL) # 只保留道路要素 optimized_layer layer.materialize(request) QgsProject.instance().addMapLayer(optimized_layer)通过属性表筛选器可以进一步精炼数据道路类型筛选条件适用场景主干道highway primary交通流量分析次干道highway secondary社区规划人行道highway footway步行系统研究3. 可视化与样式设计3.1 智能分级渲染QGIS 3.32的规则式渲染功能可以让路网自动分级显示右键图层 → 属性 → 符号化选择规则式渲染器添加以下典型规则highway motorway → 红色6px实线 highway IN (primary,trunk) → 橙色4px实线 highway secondary → 黄色3px虚线3.2 动态标签显示道路名称智能标注的关键设置大小调整数据定义覆盖 → 使用length字段动态控制字号避让规则设置道路走向为标注方向参考优先级给name:zh字段更高权重# 动态标签大小表达式示例 CASE WHEN highway motorway THEN 12 WHEN highway primary THEN 10 ELSE 8 END4. 高级应用场景4.1 路网拓扑检查利用QGIS的拓扑检查器Processing Toolbox → Vector geometry → Topology Checker可快速发现数据问题悬挂节点Dangling nodes重复线段Duplicate geometries自相交道路Self intersections4.2 网络分析准备进行可达性分析前需要确保路网拓扑正确使用几何工具 → 多部分转单部分处理复杂要素执行网络分析 → 服务区分析前务必设置正确的方向字段对于单向道路需要添加oneway字段过滤# 构建网络分析图 import networkx as nx G nx.Graph() for feature in layer.getFeatures(): if feature[highway] not in [footway,cycleway]: G.add_edge(feature[from_node], feature[to_node], weightfeature[length]/feature[maxspeed])在实际项目中这套方法帮助我快速完成了深圳前海片区15平方公里路网数据的采集与分析。最耗时的步骤其实是数据清洗建议在Overpass查询阶段就通过更精细的过滤条件减少后续工作量。例如添加[highway!~footway|cycleway|path]可以提前排除非机动车道数据。 QGIS 3.32实战：5分钟搞定OpenStreetMap城市路网数据下载与导入
QGIS 3.32实战5分钟搞定OpenStreetMap城市路网数据下载与导入在数字化城市研究与规划领域获取精准的路网数据往往成为项目启动的第一道门槛。传统方法需要手动划定经纬度范围不仅效率低下还会引入大量冗余数据。现在借助QGIS 3.32与Overpass API的黄金组合即使是GIS新手也能在咖啡冷却前完成城市路网的精准采集与可视化。1. Overpass API数据获取实战1.1 中文城市ID智能查询OpenStreetMap的Overpass API如同城市数据的搜索引擎但直接使用英文查询中国城市常会遇到匹配偏差。这里有个专业技巧通过name:zh参数锁定中文行政区划。打开Overpass Turbo在线查询工具https://overpass-turbo.eu/在左侧编辑区输入以下代码模板osm-script query typerelation has-kv kboundary vadministrative/ has-kv kname:zh v武汉市/ /query print/ /osm-script执行查询后结果面板会显示类似这样的关键信息relation id3797255 tag kname:zh v武汉市/ tag kboundary vadministrative/ /relation注意国内部分新城区可能出现多个关联ID建议选择admin_level6的市级行政区划关系。若查询超时可添加timeout25参数延长等待时间。1.2 路网数据精准下载获取城市ID后需要将其转换为Overpass特有的area引用值。这个转换公式让很多初学者困惑area_ref 3600000000 relation_id以武汉市为例ID3797255其area_ref就是3603797255。现在使用这个魔法数字构建路网查询语句osm-script timeout1800 union area-query ref3603797255/ recurse typenode-relation intorels/ recurse typenode-way/ recurse typeway-relation/ /union union item/ recurse typeway-node/ /union print modebody/ /osm-script这段代码的精妙之处在于area-query限定数据获取范围三重recurse确保获取完整拓扑关系timeout1800避免大数据量城市超时中断点击运行后通过导出按钮下载OSM格式数据。建议选择map.osm命名方便后续识别。2. QGIS 3.32数据处理技巧2.1 OSM数据智能转换新版QGIS 3.32内置的OSM数据处理器堪称黑科技。无需安装插件三步完成格式转换菜单导航图层 → 数据源管理器 → 矢量 → 协议类型选择OSM XML文件加载点击...按钮选择下载的map.osm文件图层选择勾选lines类型道路数据主要存储于此常见报错解决方案若遇到CRS未定义警告右键图层 → 属性 → 源 → 将CRS强制设置为EPSG:43262.2 路网数据深度优化原始OSM数据包含大量冗余属性需要针对性过滤# 在QGIS Python控制台执行以下脚本优化数据 layer iface.activeLayer() request QgsFeatureRequest() request.setFilterExpression(highway IS NOT NULL) # 只保留道路要素 optimized_layer layer.materialize(request) QgsProject.instance().addMapLayer(optimized_layer)通过属性表筛选器可以进一步精炼数据道路类型筛选条件适用场景主干道highway primary交通流量分析次干道highway secondary社区规划人行道highway footway步行系统研究3. 可视化与样式设计3.1 智能分级渲染QGIS 3.32的规则式渲染功能可以让路网自动分级显示右键图层 → 属性 → 符号化选择规则式渲染器添加以下典型规则highway motorway → 红色6px实线 highway IN (primary,trunk) → 橙色4px实线 highway secondary → 黄色3px虚线3.2 动态标签显示道路名称智能标注的关键设置大小调整数据定义覆盖 → 使用length字段动态控制字号避让规则设置道路走向为标注方向参考优先级给name:zh字段更高权重# 动态标签大小表达式示例 CASE WHEN highway motorway THEN 12 WHEN highway primary THEN 10 ELSE 8 END4. 高级应用场景4.1 路网拓扑检查利用QGIS的拓扑检查器Processing Toolbox → Vector geometry → Topology Checker可快速发现数据问题悬挂节点Dangling nodes重复线段Duplicate geometries自相交道路Self intersections4.2 网络分析准备进行可达性分析前需要确保路网拓扑正确使用几何工具 → 多部分转单部分处理复杂要素执行网络分析 → 服务区分析前务必设置正确的方向字段对于单向道路需要添加oneway字段过滤# 构建网络分析图 import networkx as nx G nx.Graph() for feature in layer.getFeatures(): if feature[highway] not in [footway,cycleway]: G.add_edge(feature[from_node], feature[to_node], weightfeature[length]/feature[maxspeed])在实际项目中这套方法帮助我快速完成了深圳前海片区15平方公里路网数据的采集与分析。最耗时的步骤其实是数据清洗建议在Overpass查询阶段就通过更精细的过滤条件减少后续工作量。例如添加[highway!~footway|cycleway|path]可以提前排除非机动车道数据。