欧洲大黄蜂分布点源数据栅格化遇阻，求生成.tif及复刻文章地图方法

先排查常见报错根源 #

• 确认点数据的复合概率字段是纯数值类型，无空值、字符串或异常值。打开属性表检查，若有空值，要么删除对应点要素，要么用字段计算器填充0这类合理值。
• 检查坐标系：点数据与输出栅格必须使用同一投影（推荐欧洲常用的EPSG:3035，避免地理坐标系的分辨率单位混乱）。
• 分辨率设置：若报错为内存溢出，说明分辨率设得过小（比如10m级），改成1km或5km级分辨率重试。

QGIS可视化栅格化步骤（匹配文章图2效果） #

1. 预处理点数据

   • 打开点图层属性表，筛选出复合概率字段的有效记录，删除无效行。
   • 右键点图层 → Export → Save Features As，保存为GeoPackage（确保坐标系正确）。

2. 运行Rasterize工具

   • 打开Processing Toolbox → 搜索Rasterize (Vector to Raster)
   • 参数设置：
     • Input layer：选择预处理后的点数据
     • Field to use for a burn-in value：选择复合概率对应的字段（比如comp_prob，以你数据里的实际字段名为准）
     • Output raster size units：选Map units，设置Width/Horizontal resolution和Height/Vertical resolution为1000（对应1km分辨率，匹配文章图2精度）
     • Output extent：选择Calculate from layer → 选中你的点图层，自动获取范围
     • Output raster：指定保存路径，格式选GeoTIFF
   • 点击Run，若仍报错，查看Processing Log的具体错误信息（比如字段名错误就修正字段名，内存不足就调大分辨率）

3. 匹配文章图2的可视化风格

   • 生成栅格后，打开图层属性 → Symbology → 选择Singleband pseudocolor
   • 配色方案选渐变（比如从浅蓝到深红，对应概率从低到高），调整色阶范围为0-1（复合概率的常规区间）
   • 添加欧洲边界矢量图层作为底图，与文章图2的范围对齐

GDAL命令行替代方案（更稳定） #

如果GUI工具持续报错，试试用GDAL命令行：

gdal_rasterize -a 复合概率字段名 -tr 1000 1000 -te 点数据最小x 点数据最小y 点数据最大x 点数据最大y -ot Float32 -of GTiff 输入点数据.gpkg 输出栅格.tif

• -a：指定用来生成栅格值的目标字段
• -tr：设置分辨率（1000即1km）
• -te：替换成你点数据的实际范围坐标（可从QGIS图层属性的Information面板获取）
• -ot Float32：保证概率值的精度

内容的提问来源于stack exchange，提问作者kar-g