如何将GRIB2文件转换为指定格式的GeoJSON？

GRIB2属于栅格气象数据，要转换成你需要的带LineString几何的GeoJSON，核心是先提取等值线，再调整属性与格式，用GDAL/OGR工具链即可完成，步骤如下：

1. 从GRIB2提取等值线（生成中间矢量文件） #

使用gdal_contour工具从GRIB2栅格中提取对应变量的等值线，以目标格式中的sigwaveht（有效波高）为例，假设要提取数值为1.0的等值线：

gdal_contour -a sigwaveht -fl 1.0 input.grib2 contour_temp.shp

• -a sigwaveht：给生成的等值线要素添加属性字段sigwaveht，值为指定的等值线数值
• -fl 1.0：指定要提取的等值线数值，若需多个值用空格分隔（如-fl 1.0 2.0）
• 输出的contour_temp.shp为临时Shapefile，包含等值线的LineString几何与sigwaveht属性

2. 转换并调整为目标GeoJSON格式 #

用你熟悉的ogr2ogr工具，将临时Shapefile转换为目标格式，同时自动添加自增的ID属性：

ogr2ogr -f GeoJSON -sql "SELECT *, @ogr_fid AS ID FROM contour_temp" output.geojson contour_temp.shp

• -f GeoJSON：指定输出格式为GeoJSON
• -sql：通过SQL语句给每个要素添加ID字段，值为OGR内部的要素ID（从0开始自增，匹配目标格式要求）
• output.geojson即为最终符合要求的文件

3. 一步完成（跳过临时文件） #

若不想生成中间Shapefile，可借助GDAL虚拟文件系统实现工具间直接传参：

gdal_contour -a sigwaveht -fl 1.0 input.grib2 /vsistdout/ | ogr2ogr -f GeoJSON -sql "SELECT *, @ogr_fid AS ID FROM /vsistdin/" output.geojson /vsistdin/

• /vsistdout/和/vsistdin/是GDAL虚拟文件系统，实现数据无临时文件传递

注意事项 #

• 确认GDAL支持GRIB2：执行gdalinfo --formats查看是否包含GRIB格式，若未包含需重新编译GDAL时添加GRIB驱动（通常默认自带）
• 多变量GRIB2需指定波段：先用gdalinfo input.grib2查看波段对应的变量，再用-b <波段号>参数指定处理的波段（如gdal_contour -b 1 ...）

内容的提问来源于stack exchange，提问作者jms1980