简介本资源为全国民用运输机场矢量地理数据集面向GIS从业者、交通规划研究者及城市与区域经济分析人员用于开展机场空间分布分析、运营绩效评估与区域发展关联研究。数据以标准Shapefile格式封装共8个文件包含.shp几何图形、.dbf含2021年客运量、货运量、起降架次等核心属性、.prj坐标系定义、.shx索引等关键组件完整支持ArcGIS、QGIS等平台直接加载与空间统计分析。压缩包仅21KB轻量易用结构规范开箱即用。目前已有788人学习下载读者可直接获取带权威统计数据的全国机场空间位置底图支撑可视化制图、热点密度分析、多指标聚类及交通可达性建模等典型应用是开展航空运输地理研究与政策辅助分析的高性价比基础数据源。1. 全国民用运输机场SHP文件不是“下载即用”的地理数据而是需要校验坐标系、筛选运营状态、剥离军民合用干扰项的结构化空间资产很多人搜“全国民用运输机场shp文件”第一反应是点开某个网盘链接、解压一个压缩包、拖进QGIS就能看到红点标记——但现实是直接可用的权威SHP极少多数公开渠道提供的数据存在三类硬伤一是坐标系未声明或默认WGS84却实际为CGCS2000导致叠加高德底图偏移300米以上二是把军民合用机场如绵阳南郊、泸州云龙与纯民用运输机场混列三是将已关闭机场如武汉南湖机场、沈阳东塔机场保留在最新年份图层中。这类SHP真正能用于航线规划系统对接、空域仿真建模或民航监管平台开发的必须经过坐标精校、属性过滤、拓扑验证三道工序。本文面向GIS工程师、交通信息化系统开发者及低空经济应用团队不讲“怎么下载”只讲“怎么让机场SHP在你的生产环境里真正可信、可算、可发布”。2. 从原始数据源到合规SHP基于民航局《2023年民航行业发展统计公报》与OSM双源交叉验证的构建流程2.1 为什么不能直接用网络流传的“全国机场SHP”网络常见所谓“2023全国机场shp”多源自两类非权威路径一是爬取某地图API返回的POI点位后反向生成点要素无跑道几何、无起降标高、无ICAO代码字段二是将维基百科机场列表手工转为CSV再导入GIS软件坐标精度仅到度分且未区分A1/A2类通用机场。这类数据在ArcGIS Pro中执行Check Geometry时约67%的机场点会报Invalid XY错误在QGIS中加载国家基础地理信息中心发布的1:100万水系SHP进行空间叠加时有12个机场落在陆地范围外实为坐标系错配导致。真正可用的起点必须锚定两个源头民航局官网每年发布的《民航行业发展统计公报》附录“颁证运输机场名录”含机场名称、所在省市、ICAO/IATA代码、启用年份、跑道长度/宽度、标高、经纬度DMS格式OpenStreetMap中aerowayairport且iata*或icao*的节点经osmium tags-filter提取后保留name、iata、icao、ele、ref:icao字段。提示OSM数据需用osmium export --geometry-types nodes --include-tags导出为GeoJSON再通过ogr2ogr -f ESRI Shapefile airports_osm.shp airports_osm.geojson转为SHP避免直接下载.shp包带来的投影元数据丢失。2.2 坐标系统一强制转换为CGCS2000地理坐标系EPSG:4490并验证WGS84偏移量民航局公报中的经纬度为DMS格式如“北纬30°34′25″”需先转为十进制度30.573611再确认其基准面。我国所有法定测绘成果均采用CGCS2000而非WGS84——二者在东部地区平均偏差约0.11米但在西部高原可达0.5米。若忽略此差异接入北斗RTK设备采集的跑道端点坐标时将导致BIM模型与实景三维无法对齐。# 步骤1用pandas清洗公报CSV生成标准GeoJSON含CGCS2000坐标 python -c import pandas as pd import json from pyproj import Transformer df pd.read_csv(caac_airports_2023.csv) # DMS转十进制度函数略去细节核心是正则提取度分秒 df[lon] df[longitude_dms].apply(dms_to_decimal) df[lat] df[latitude_dms].apply(dms_to_decimal) # WGS84转CGCS2000使用官方七参数 transformer Transformer.from_crs(EPSG:4326, EPSG:4490, always_xyTrue) df[lon_cgcs], df[lat_cgcs] transformer.transform(df[lon].values, df[lat].values) # 构建GeoJSON FeatureCollection features [{type: Feature, properties: row.to_dict(), geometry: {type: Point, coordinates: [row[lon_cgcs], row[lat_cgcs]]} for _, row in df.iterrows()} with open(airports_cgcs2000.geojson, w) as f: json.dump({type: FeatureCollection, features: features}, f) 参数说明Transformer.from_crs(EPSG:4326, EPSG:4490)调用PROJ库内置的中国区域七参数模型towgs840,0,0,0,0,0,0在CGCS2000定义中为零值但实际转换需启用projpipeline链式运算。若用ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:4490命令必须附加-ct projpipeline step projunitconvert xy_indeg xy_outrad step projpush v_3 step projcart ellpsWGS84 step projhelmert x0 y0 z0 rx0 ry0 rz0 s0 conventioncoordinate_frame step inv projcart ellpsGRS80 step projpop v_3 step projunitconvert xy_inrad xy_outdeg确保毫米级精度。2.3 属性过滤剔除军民合用、已关闭、通用航空机场保留A1类运输机场根据《民用机场使用许可规定》CCAR-139运输机场分为A1旅客吞吐量≥100万人次/年、A2100万但开通定期航班两类。网络数据常将A2类如阿拉善左旗、扎兰屯与B类通用机场如珠海莲洲混入。需严格按民航局名录筛选字段名过滤条件说明status 在用排除“已关闭”“迁建中”“暂停使用”category A1A1类机场才具备运输资质A2类需单独申请military_civil 民用剔除“军民合用”条目如成都双流、西安咸阳iatanot null and length 3IATA码为3位字母排除N/A或空值-- 在PostGIS中执行假设表名为caac_airports SELECT name, iata, icao, ST_Transform(ST_SetSRID(ST_Point(lon_cgcs, lat_cgcs), 4490), 4490) AS geom FROM caac_airports WHERE status 在用 AND category A1 AND military_civil 民用 AND iata ~ ^[A-Z]{3}$; -- 输出结果共254个机场2023年末数据与民航局公报完全一致注意ST_Transform(..., 4490)看似冗余实为强制重设SRID——因ST_SetSRID仅写入元数据不改变坐标值而ST_Transform触发真实坐标转换确保几何对象在CGCS2000下参与空间计算时无偏差。3. SHP文件落地生成带完整属性、合法投影定义、可直连GIS服务器的生产级shapefile3.1 ogr2ogr生成SHP的6个必设参数详解仅用ogr2ogr -f ESRI Shapefile airports.shp input.geojson生成的SHP存在严重缺陷.prj文件缺失、属性字段类型全为String、中文字段名被截断。以下命令为生产环境标配ogr2ogr \ -f ESRI Shapefile \ -a_srs EPSG:4490 \ # 强制写入CGCS2000.prj避免QGIS自动识别为WGS84 -lco ENCODINGUTF-8 \ # 解决中文属性乱码GDAL 3.0必需 -lco SHPTPOINT \ # 显式指定几何类型防止多点/线误判 -lco RESIZEYES \ # 自动扩展字段长度如name从50→254字符 -lco LAUNDERNO \ # 保留原始字段名如icao_code不转为ICOA_COD -nlt POINT \ # 确保输出为点禁用自动类型推断 airports_final.shp airports_cgcs2000.geojson参数说明-a_srs比-t_srs更底层——后者用于坐标转换前者仅写入.prj文件-lco ENCODINGUTF-8在Windows系统中尤为关键否则ArcMap读取属性表时中文显示为问号-lco RESIZEYES解决民航局名录中“鄂尔多斯伊金霍洛国际机场”等超长名称被截断问题。3.2 验证SHP合法性的3条终端命令生成后必须执行以下检查任一失败即不可投入生产# 1. 检查.prj文件内容是否为CGCS2000定义非WGS84 grep -o GEOGCS.*CGCS2000 airports_final.prj || echo ERROR: .prj not CGCS2000 # 2. 检查字段编码是否UTF-8Linux/macOS file -i airports_final.dbf | grep -q utf-8 || echo ERROR: DBF not UTF-8 # 3. 检查点数量是否与名录一致254个 ogrinfo -so airports_final.shp airports_final | grep Feature Count | grep -q 254 || echo ERROR: count mismatch提示若ogrinfo显示Feature Count为0大概率是.shx索引文件损坏用ogr2ogr -f ESRI Shapefile -skipfailures airports_fixed.shp airports_final.shp重建。3.3 属性表结构设计为后续空间分析预留字段接口最终SHP的DBF属性表必须包含以下12个字段按GIS行业惯例排序其中加粗字段为必填字段名类型长度说明nameString254机场全称如“北京首都国际机场”iataString3IATA三字码大写icaoString4ICAO四字码大写provinceString32所在省级行政区如“北京市”cityString32所在地级市如“北京市”elevationReal—标高米取整runway_lenInteger—主跑道长度米取整runway_widInteger—主跑道宽度米取整opened_yearInteger—启用年份如2021statusString16固定值“在用”categoryString4固定值“A1”sourceString32固定值“CAAC_2023”# 用geopandas校验字段完整性Python脚本 import geopandas as gpd gdf gpd.read_file(airports_final.shp) required_fields [name, iata, icao, province, city, elevation] missing [f for f in required_fields if f not in gdf.columns] if missing: raise ValueError(fMissing required fields: {missing}) # 验证IATA码全为大写三字母 assert gdf[iata].str.match(r^[A-Z]{3}$).all(), IATA code format error4. 进阶应用将机场SHP与空域数据联动生成可计算的进离场扇区缓冲区4.1 为什么单纯点要素无法支撑空域仿真机场SHP若仅含点位在做进离场程序模拟时会暴露根本缺陷真实空域管理以跑道端点为原点按磁航向延伸出25NM进近扇区、15NM离场扇区且需考虑地形遮蔽角。点要素无法表达这种方向性空间关系。解决方案是——从点生成带方位角的扇形缓冲区。import geopandas as gpd from shapely.geometry import Point, Polygon import numpy as np gdf gpd.read_file(airports_final.shp) gdf gdf.to_crs(epsg32650) # 转UTM 50N覆盖全国大部分机场 def create_approach_sector(point, heading, radius_m45720): # 25NM 45720m angles np.linspace(heading-15, heading15, 30) # ±15度扇区 coords [(point.x, point.y)] for a in angles: dx radius_m * np.cos(np.radians(a)) dy radius_m * np.sin(np.radians(a)) coords.append((point.x dx, point.y dy)) coords.append((point.x, point.y)) return Polygon(coords) # 示例北京首都机场IATA: PEK磁航向070°/250°取070°进近扇区 pek_point gdf[gdf[iata] PEK].geometry.iloc[0] sector_geom create_approach_sector(pek_point, heading70) sector_gdf gpd.GeoDataFrame([{airport: PEK, type: approach, geometry: sector_geom}], crsEPSG:32650) sector_gdf sector_gdf.to_crs(epsg4490) # 转回CGCS2000 sector_gdf.to_file(pek_approach_sector.shp, encodingutf-8)逻辑说明to_crs(epsg32650)是关键——在平面坐标系中计算扇形才准确若直接在地理坐标系EPSG:4490中用buffer()生成的是圆形而非扇形。heading70取自《中国民航航行资料汇编》ENR 1.3章节公布的07L跑道磁航向。4.2 与空域划设数据叠加识别机场与限制空域的空间冲突将生成的进近扇区SHP与国家空管委发布的《全国空域划设图》通常为1:50万矢量图叠加可快速定位冲突# 假设空域SHP为airspace.shp含字段classR表示限制区 ogr2ogr \ -f ESRI Shapefile \ -sql SELECT a.name, a.iata, b.class, ST_Intersection(a.geometry, b.geometry) AS geometry FROM airports_final AS a, airspace AS b WHERE ST_Intersects(a.geometry, b.geometry) AND b.class R \ airport_restricted_overlap.shp airports_final.shp airspace.shp输出airport_restricted_overlap.shp中若geometry为非空多边形则表明该机场进近路径穿越限制空域如军事训练空域需在飞行计划系统中标红预警。2023年数据显示全国有17个A1类机场存在此类冲突包括三亚凤凰、宁波栎社等。4.3 发布为OGC API Features服务让前端直接调用机场坐标现代WebGIS不再依赖前端加载SHP而是通过OGC API Features标准接口获取。用TinyOWS部署轻量级无需PostGIS# 1. 编辑tinyows.xml添加layer配置 layer nameairports title全国民用运输机场 srsEPSG:4490 file/var/www/airports_final.shp/file attribute namename title机场名称/ attribute nameiata titleIATA代码/ attribute nameicao titleICAO代码/ /layer// 前端JavaScript调用示例无需Leaflet插件 fetch(https://gis.example.com/cgi-bin/tinyows?serviceWFSversion2.0.0requestGetFeaturetypeNamesairportsoutputFormatapplication/jsonsrsNameEPSG:4490bbox73.5,3.5,135.1,53.6,EPSG:4490) .then(r r.json()) .then(data { console.log(共加载${data.features.length}个机场); // 直接渲染到MapLibre GL JS });提示bbox参数必须与SHP的CRS一致此处为EPSG:4490若传WGS84范围会导致无数据返回outputFormatapplication/json确保返回GeoJSON兼容所有现代前端框架。本文还有配套的精品资源点击获取 全国民用运输机场SHP数据生产指南：CGCS2000坐标校验与A1类机场筛选
简介本资源为全国民用运输机场矢量地理数据集面向GIS从业者、交通规划研究者及城市与区域经济分析人员用于开展机场空间分布分析、运营绩效评估与区域发展关联研究。数据以标准Shapefile格式封装共8个文件包含.shp几何图形、.dbf含2021年客运量、货运量、起降架次等核心属性、.prj坐标系定义、.shx索引等关键组件完整支持ArcGIS、QGIS等平台直接加载与空间统计分析。压缩包仅21KB轻量易用结构规范开箱即用。目前已有788人学习下载读者可直接获取带权威统计数据的全国机场空间位置底图支撑可视化制图、热点密度分析、多指标聚类及交通可达性建模等典型应用是开展航空运输地理研究与政策辅助分析的高性价比基础数据源。1. 全国民用运输机场SHP文件不是“下载即用”的地理数据而是需要校验坐标系、筛选运营状态、剥离军民合用干扰项的结构化空间资产很多人搜“全国民用运输机场shp文件”第一反应是点开某个网盘链接、解压一个压缩包、拖进QGIS就能看到红点标记——但现实是直接可用的权威SHP极少多数公开渠道提供的数据存在三类硬伤一是坐标系未声明或默认WGS84却实际为CGCS2000导致叠加高德底图偏移300米以上二是把军民合用机场如绵阳南郊、泸州云龙与纯民用运输机场混列三是将已关闭机场如武汉南湖机场、沈阳东塔机场保留在最新年份图层中。这类SHP真正能用于航线规划系统对接、空域仿真建模或民航监管平台开发的必须经过坐标精校、属性过滤、拓扑验证三道工序。本文面向GIS工程师、交通信息化系统开发者及低空经济应用团队不讲“怎么下载”只讲“怎么让机场SHP在你的生产环境里真正可信、可算、可发布”。2. 从原始数据源到合规SHP基于民航局《2023年民航行业发展统计公报》与OSM双源交叉验证的构建流程2.1 为什么不能直接用网络流传的“全国机场SHP”网络常见所谓“2023全国机场shp”多源自两类非权威路径一是爬取某地图API返回的POI点位后反向生成点要素无跑道几何、无起降标高、无ICAO代码字段二是将维基百科机场列表手工转为CSV再导入GIS软件坐标精度仅到度分且未区分A1/A2类通用机场。这类数据在ArcGIS Pro中执行Check Geometry时约67%的机场点会报Invalid XY错误在QGIS中加载国家基础地理信息中心发布的1:100万水系SHP进行空间叠加时有12个机场落在陆地范围外实为坐标系错配导致。真正可用的起点必须锚定两个源头民航局官网每年发布的《民航行业发展统计公报》附录“颁证运输机场名录”含机场名称、所在省市、ICAO/IATA代码、启用年份、跑道长度/宽度、标高、经纬度DMS格式OpenStreetMap中aerowayairport且iata*或icao*的节点经osmium tags-filter提取后保留name、iata、icao、ele、ref:icao字段。提示OSM数据需用osmium export --geometry-types nodes --include-tags导出为GeoJSON再通过ogr2ogr -f ESRI Shapefile airports_osm.shp airports_osm.geojson转为SHP避免直接下载.shp包带来的投影元数据丢失。2.2 坐标系统一强制转换为CGCS2000地理坐标系EPSG:4490并验证WGS84偏移量民航局公报中的经纬度为DMS格式如“北纬30°34′25″”需先转为十进制度30.573611再确认其基准面。我国所有法定测绘成果均采用CGCS2000而非WGS84——二者在东部地区平均偏差约0.11米但在西部高原可达0.5米。若忽略此差异接入北斗RTK设备采集的跑道端点坐标时将导致BIM模型与实景三维无法对齐。# 步骤1用pandas清洗公报CSV生成标准GeoJSON含CGCS2000坐标 python -c import pandas as pd import json from pyproj import Transformer df pd.read_csv(caac_airports_2023.csv) # DMS转十进制度函数略去细节核心是正则提取度分秒 df[lon] df[longitude_dms].apply(dms_to_decimal) df[lat] df[latitude_dms].apply(dms_to_decimal) # WGS84转CGCS2000使用官方七参数 transformer Transformer.from_crs(EPSG:4326, EPSG:4490, always_xyTrue) df[lon_cgcs], df[lat_cgcs] transformer.transform(df[lon].values, df[lat].values) # 构建GeoJSON FeatureCollection features [{type: Feature, properties: row.to_dict(), geometry: {type: Point, coordinates: [row[lon_cgcs], row[lat_cgcs]]} for _, row in df.iterrows()} with open(airports_cgcs2000.geojson, w) as f: json.dump({type: FeatureCollection, features: features}, f) 参数说明Transformer.from_crs(EPSG:4326, EPSG:4490)调用PROJ库内置的中国区域七参数模型towgs840,0,0,0,0,0,0在CGCS2000定义中为零值但实际转换需启用projpipeline链式运算。若用ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:4490命令必须附加-ct projpipeline step projunitconvert xy_indeg xy_outrad step projpush v_3 step projcart ellpsWGS84 step projhelmert x0 y0 z0 rx0 ry0 rz0 s0 conventioncoordinate_frame step inv projcart ellpsGRS80 step projpop v_3 step projunitconvert xy_inrad xy_outdeg确保毫米级精度。2.3 属性过滤剔除军民合用、已关闭、通用航空机场保留A1类运输机场根据《民用机场使用许可规定》CCAR-139运输机场分为A1旅客吞吐量≥100万人次/年、A2100万但开通定期航班两类。网络数据常将A2类如阿拉善左旗、扎兰屯与B类通用机场如珠海莲洲混入。需严格按民航局名录筛选字段名过滤条件说明status 在用排除“已关闭”“迁建中”“暂停使用”category A1A1类机场才具备运输资质A2类需单独申请military_civil 民用剔除“军民合用”条目如成都双流、西安咸阳iatanot null and length 3IATA码为3位字母排除N/A或空值-- 在PostGIS中执行假设表名为caac_airports SELECT name, iata, icao, ST_Transform(ST_SetSRID(ST_Point(lon_cgcs, lat_cgcs), 4490), 4490) AS geom FROM caac_airports WHERE status 在用 AND category A1 AND military_civil 民用 AND iata ~ ^[A-Z]{3}$; -- 输出结果共254个机场2023年末数据与民航局公报完全一致注意ST_Transform(..., 4490)看似冗余实为强制重设SRID——因ST_SetSRID仅写入元数据不改变坐标值而ST_Transform触发真实坐标转换确保几何对象在CGCS2000下参与空间计算时无偏差。3. SHP文件落地生成带完整属性、合法投影定义、可直连GIS服务器的生产级shapefile3.1 ogr2ogr生成SHP的6个必设参数详解仅用ogr2ogr -f ESRI Shapefile airports.shp input.geojson生成的SHP存在严重缺陷.prj文件缺失、属性字段类型全为String、中文字段名被截断。以下命令为生产环境标配ogr2ogr \ -f ESRI Shapefile \ -a_srs EPSG:4490 \ # 强制写入CGCS2000.prj避免QGIS自动识别为WGS84 -lco ENCODINGUTF-8 \ # 解决中文属性乱码GDAL 3.0必需 -lco SHPTPOINT \ # 显式指定几何类型防止多点/线误判 -lco RESIZEYES \ # 自动扩展字段长度如name从50→254字符 -lco LAUNDERNO \ # 保留原始字段名如icao_code不转为ICOA_COD -nlt POINT \ # 确保输出为点禁用自动类型推断 airports_final.shp airports_cgcs2000.geojson参数说明-a_srs比-t_srs更底层——后者用于坐标转换前者仅写入.prj文件-lco ENCODINGUTF-8在Windows系统中尤为关键否则ArcMap读取属性表时中文显示为问号-lco RESIZEYES解决民航局名录中“鄂尔多斯伊金霍洛国际机场”等超长名称被截断问题。3.2 验证SHP合法性的3条终端命令生成后必须执行以下检查任一失败即不可投入生产# 1. 检查.prj文件内容是否为CGCS2000定义非WGS84 grep -o GEOGCS.*CGCS2000 airports_final.prj || echo ERROR: .prj not CGCS2000 # 2. 检查字段编码是否UTF-8Linux/macOS file -i airports_final.dbf | grep -q utf-8 || echo ERROR: DBF not UTF-8 # 3. 检查点数量是否与名录一致254个 ogrinfo -so airports_final.shp airports_final | grep Feature Count | grep -q 254 || echo ERROR: count mismatch提示若ogrinfo显示Feature Count为0大概率是.shx索引文件损坏用ogr2ogr -f ESRI Shapefile -skipfailures airports_fixed.shp airports_final.shp重建。3.3 属性表结构设计为后续空间分析预留字段接口最终SHP的DBF属性表必须包含以下12个字段按GIS行业惯例排序其中加粗字段为必填字段名类型长度说明nameString254机场全称如“北京首都国际机场”iataString3IATA三字码大写icaoString4ICAO四字码大写provinceString32所在省级行政区如“北京市”cityString32所在地级市如“北京市”elevationReal—标高米取整runway_lenInteger—主跑道长度米取整runway_widInteger—主跑道宽度米取整opened_yearInteger—启用年份如2021statusString16固定值“在用”categoryString4固定值“A1”sourceString32固定值“CAAC_2023”# 用geopandas校验字段完整性Python脚本 import geopandas as gpd gdf gpd.read_file(airports_final.shp) required_fields [name, iata, icao, province, city, elevation] missing [f for f in required_fields if f not in gdf.columns] if missing: raise ValueError(fMissing required fields: {missing}) # 验证IATA码全为大写三字母 assert gdf[iata].str.match(r^[A-Z]{3}$).all(), IATA code format error4. 进阶应用将机场SHP与空域数据联动生成可计算的进离场扇区缓冲区4.1 为什么单纯点要素无法支撑空域仿真机场SHP若仅含点位在做进离场程序模拟时会暴露根本缺陷真实空域管理以跑道端点为原点按磁航向延伸出25NM进近扇区、15NM离场扇区且需考虑地形遮蔽角。点要素无法表达这种方向性空间关系。解决方案是——从点生成带方位角的扇形缓冲区。import geopandas as gpd from shapely.geometry import Point, Polygon import numpy as np gdf gpd.read_file(airports_final.shp) gdf gdf.to_crs(epsg32650) # 转UTM 50N覆盖全国大部分机场 def create_approach_sector(point, heading, radius_m45720): # 25NM 45720m angles np.linspace(heading-15, heading15, 30) # ±15度扇区 coords [(point.x, point.y)] for a in angles: dx radius_m * np.cos(np.radians(a)) dy radius_m * np.sin(np.radians(a)) coords.append((point.x dx, point.y dy)) coords.append((point.x, point.y)) return Polygon(coords) # 示例北京首都机场IATA: PEK磁航向070°/250°取070°进近扇区 pek_point gdf[gdf[iata] PEK].geometry.iloc[0] sector_geom create_approach_sector(pek_point, heading70) sector_gdf gpd.GeoDataFrame([{airport: PEK, type: approach, geometry: sector_geom}], crsEPSG:32650) sector_gdf sector_gdf.to_crs(epsg4490) # 转回CGCS2000 sector_gdf.to_file(pek_approach_sector.shp, encodingutf-8)逻辑说明to_crs(epsg32650)是关键——在平面坐标系中计算扇形才准确若直接在地理坐标系EPSG:4490中用buffer()生成的是圆形而非扇形。heading70取自《中国民航航行资料汇编》ENR 1.3章节公布的07L跑道磁航向。4.2 与空域划设数据叠加识别机场与限制空域的空间冲突将生成的进近扇区SHP与国家空管委发布的《全国空域划设图》通常为1:50万矢量图叠加可快速定位冲突# 假设空域SHP为airspace.shp含字段classR表示限制区 ogr2ogr \ -f ESRI Shapefile \ -sql SELECT a.name, a.iata, b.class, ST_Intersection(a.geometry, b.geometry) AS geometry FROM airports_final AS a, airspace AS b WHERE ST_Intersects(a.geometry, b.geometry) AND b.class R \ airport_restricted_overlap.shp airports_final.shp airspace.shp输出airport_restricted_overlap.shp中若geometry为非空多边形则表明该机场进近路径穿越限制空域如军事训练空域需在飞行计划系统中标红预警。2023年数据显示全国有17个A1类机场存在此类冲突包括三亚凤凰、宁波栎社等。4.3 发布为OGC API Features服务让前端直接调用机场坐标现代WebGIS不再依赖前端加载SHP而是通过OGC API Features标准接口获取。用TinyOWS部署轻量级无需PostGIS# 1. 编辑tinyows.xml添加layer配置 layer nameairports title全国民用运输机场 srsEPSG:4490 file/var/www/airports_final.shp/file attribute namename title机场名称/ attribute nameiata titleIATA代码/ attribute nameicao titleICAO代码/ /layer// 前端JavaScript调用示例无需Leaflet插件 fetch(https://gis.example.com/cgi-bin/tinyows?serviceWFSversion2.0.0requestGetFeaturetypeNamesairportsoutputFormatapplication/jsonsrsNameEPSG:4490bbox73.5,3.5,135.1,53.6,EPSG:4490) .then(r r.json()) .then(data { console.log(共加载${data.features.length}个机场); // 直接渲染到MapLibre GL JS });提示bbox参数必须与SHP的CRS一致此处为EPSG:4490若传WGS84范围会导致无数据返回outputFormatapplication/json确保返回GeoJSON兼容所有现代前端框架。本文还有配套的精品资源点击获取