GDAL gdal_rasterize 怎么用：像元对齐、烧写规则与矢栅一致性检查
把地块、道路或掩膜面转成栅格时，最容易出现的不是命令报错，而是边界错半个像元、类别被覆盖、输出范围与分析栅格对不上。这些问题会在后续统计或机器学习里被放大。gdal_rasterize 的关键是先定义目标网格，再决定要烧什么值、怎样处理重叠要素。本文给出从参考栅格到结果核验的完整判断路径。
目标栅格网格必须先固定
分辨率、范围、原点、像元对齐和 CRS 共同决定每个像元位置。只指定 -tr 而不对齐参考栅格，会让看似相同分辨率的结果产生半像元偏移。
烧写值反映业务语义
固定 burn 值适合二值掩膜，-a 字段适合分类或连续属性；多个要素重叠时，后写覆盖前写并不一定符合业务，需先定义优先级或预处理。
ALL_TOUCHED 改变面积口径
默认按像元中心判定，ALL_TOUCHED 会包含所有接触像元。两者没有绝对对错，但面积统计、训练标签和边界小对象会明显不同，不能混用。
栅格化不修复输入问题
无效面、错误 CRS、重复地块或空属性不会因为变成像元而消失。输入质量应在矢量阶段先检查。
可执行实操流程
- 选定分析用参考栅格，读取其 CRS、extent、pixel size 和左上角原点；将这些参数写进处理配置。
- 检查输入矢量的 CRS、无效几何、烧写字段空值和重叠关系；重叠有优先级时先排序或生成无重叠结果。
- 使用 -te、-tr、-tap 和 -a 或 -burn 明确输出网格与值规则；先对一个小范围执行并可视化。
- 将输出与参考栅格叠加，检查四角、细长地物和边界处是否像元对齐；统计各类别像元数与预期面积量级。
- 把命令、GDAL 版本、输入数据 hash、范围和 NoData 写入日志；数据更新后复跑同样统计。
gdal_rasterize -a landuse -te xmin ymin xmax ymax
-tr 10 10 -tap -a_nodata 0 -ot Byte
parcels.gpkg landuse_10m.tif项目避坑与质量检查
不少人将 -at 当作“更精确”。事实上它会让窄条要素更容易占像元，也会系统性增大面积。若分类结果要与既有中心点判定产品比较，就不应擅自切换 -at；先在边界样本上量化差异，并将选择依据写进方法说明。
| 检查阶段 | 必须保留的证据 | 异常处理 |
|---|---|---|
| 输入 | 参考栅格 CRS、原点、范围与像元大小 | 隔离异常样本，不直接覆盖源数据 |
| 处理 | 烧写字段空值、重叠策略与 NoData | 回到参数、单位与筛选条件逐项复现 |
| 交付 | 边界抽样、类别像元数和对齐截图 | 用独立样本或第二个环境复核 |
让流程能被下一位同事复跑
留存的不只是输出 TIFF，还应包括参考网格定义、命令行、输入版本和 ALL_TOUCHED 选择理由。相同地块换一套网格就可能得到不同统计，复现必须能还原这个前提。
把源数据批次、坐标参考、关键参数、异常清单和前后统计放在同一份处理记录中。这样数据更新时，团队判断的是结果差异来自哪里，而不是重新猜测上一次做过什么。
抽样复核与交接记录
若成果将参与多期变化检测，应把每一期都锁定到同一参考网格，不能让单期按自身范围自动扩展。除比较类别面积外，还应在固定小窗口统计像元值交叉表；某一类突然大量落在相邻类别，常比总面积变化更早暴露坐标偏移或字段映射错误。
FAQ
-tap 有什么作用？
它将输出范围扩展到像元网格边界，帮助与同分辨率网格对齐；仍需与参考栅格原点和范围一起验证。
重叠面会输出哪个值？
通常后绘制的要素覆盖前者。应先明确业务优先级，不能依赖文件顺序。
为何栅格面积与矢量面积不同？
像元中心/接触规则、分辨率和边界效应都会造成差异。应使用同一规则比较，并报告误差。
总结
高质量栅格化从来不是一条命令，而是一份可追溯的网格契约。对齐、烧写规则和边界抽检到位，矢栅转换才能支撑后续分析。