1.
说明:本文把“大豆苗带宽度”理解为空间频谱上的最高显著空间频率(如行间、株间、叶片尺度变化),并把“采样频率”理解为空间采样点密度(单位:点/米或采样间距m)。
小分段:为什么重要——若采样频率不足,会产生混叠,丢失小尺度变异;过高则浪费资源。
2.
传感器:高分辨率RGB相机或多光谱相机;可选光谱辐射计。定位:RTK-GNSS或差分GPS。平台:UAV(无人机)或田间手持滑轨。
小分段:计算与记录工具——笔记本、GIS软件(QGIS/ArcGIS)、MATLAB或Python(numpy、scipy、rasterio)。
3.
步骤1:实地观察并测量典型最小结构尺寸。例如:行距通常为0.38 m,单株冠幅0.1–0.3 m。记录最小感兴趣结构Lmin(米)。
小分段:经验转换——空间带宽(最高空间频率 fmax)≈ 1/(2*Lmin)。若感兴趣结构0.25 m,则 fmax≈2 cycles/m。
4.
奈奎斯特:采样频率fs(samples/m)应 ≥ 2 * fmax。换算为采样间距Δx(m)则 Δx ≤ 1/fs。
小分段:示例计算——若fmax=2 cycles/m,需fs≥4 samples/m,Δx≤0.25 m。即每0.25m一个样点可保证不混叠。
5.
原则:采样带宽(测带宽度)至少覆盖若干行——建议横向宽度≥3×行距,以包含行内与行间变异。对于0.38 m行距,带宽≥1.2 m。
小分段:若使用遥感影像,应保证横向地面分辨率(GSD)小于或等于Δx,并且航带重叠率确保无盲区。
6.
步骤1:计算所需像素大小(GSD)= 高度H(m) * 像元尺寸(mm) / (镜头焦距(mm) * 1000)。确保GSD ≤ Δx。
步骤2:计算沿轴向采样间距 = 飞行速度V (m/s) / 相机触发频率R (Hz)。保证该间距 ≤ Δx。
7.
实例:目标Δx=0.25m,使用相机触发频率R=5Hz,求最大速度Vmax=Δx * R = 1.25 m/s。若想飞2 m/s,则R需≥8Hz。
小分段:若相机频率不能满足,可降低飞行高度(降低GSD)或改用高频触发和滚动快门同步方案。
8.
采样前:如使用扫描器或相机,尽量用光学或机动方式进行低通(如稍降低镜头分辨率、加小光圈或使用低通滤镜)。
采样后:数字上用低通滤波(例如高斯平滑、移动平均或频域滤波)去除高于奈奎斯特带宽的成分,再下采样以防混叠。
9.
步骤A:现场勘测并标定坐标系,放置基准点(RTK定位)。记录行向与行距。
步骤B:设定目标最小结构Lmin,计算fmax、fs和Δx;按第6节匹配相机频率与飞行/移动速度。
步骤C:完成试飞/试采样,下载影像,检查GSD与重叠;若不满足,调整高度或速度重新采集。
10.
处理步骤:影像拼接(orthomosaic)、几何校正、辐射校正、提取行间/株间剖面并对剖面做FFT(快速傅里叶变换)估算实际空间带宽。
小分段:质量控制——检查能否在频域看到预期fmax峰值,若存在高频噪声超标,回到反混叠步骤重新处理或补采。
11.
建议:先做小范围试点获取Lmin与频谱,再据此确定全田采样密度与航次。预算要覆盖RTK、传感器频率升级及数据存储/计算资源。
小分段:重复与复测——每块试验地至少三次重复采样以评估变异与系统误差。
12.
问:如果田间最小感兴趣结构不确定,该如何设定采样频率?

答:先做高分辨率试点(尽可能小的GSD,例如0.05–0.1m)并对样片做频谱分析,估算fmax;然后按奈奎斯特乘以安全系数1.5–2反推全场采样频率,能有效避免低估带宽。
13.
问:手持或地面滑轨采样与UAV采样在带宽需求上有何不同?
答:本质相同,区别在于控制采样间距的方式:地面平台可更精确控制位置以达到小Δx,但覆盖效率低;UAV覆盖快但需严格匹配触发频率与速度并注意振动带来的高频噪声,需加滤波。
14.
问:如何在有限资源下在保证不混叠的前提下降低采样成本?
答:采取策略:1) 先做频谱测定决定必要的fs;2) 在关键区域高密度采样,在其他区域用较低密度并做插值;3) 使用低通滤波减少必需的fs;4) 合理设计航线与重叠率,减少重复覆盖浪费。