冬小麦抽穗前单位面积穗数( PNPA )的快速准确估测对于评估产量潜力和调控作物生长以提高最终产量至关重要。由于忽略了光谱饱和和背景影响，现有方法利用抽穗前的遥感数据估算PNPA的精度较低，本研究提出了一种基于无人机多光谱影像的光谱-纹理PNPA敏感指数( spectral-textural PNPA sensitive index，SPSI )，以降低冬小麦抽穗前的光谱饱和度，提高PNPA的估算精度，结合背景物质对纹理指数( TI )估算PNPA的影响，并将最优光谱指数( SI )与TI相结合，构建了综合指数SPSI。随后，通过与其他指标( SI和TIs)的比较来评估SPSI的性能。结果表明，在8种纹理特征的所有指数中，除了TI_{[ HOM ]}，TI_{[ ENT ]}和TI_{[ SEM ]}外，绿像素TI比全像素TI具有更好的性能，与DATT_{[ 850，730，675]}、TI_{NDRE [ MEA ]}和NDTI_{COR[ 850，730 ]}相比，由公式DATT_{[ 850，730，675]}+ NDTI_{COR[ 850，730 ]}计算的SPSI在任何数据集中的任何日期都具有最高的总体精度。对于组装了两个实验数据集的统一模型，Rv²的SPSI值增加了0.11~0.23，RMSE和RRMSE在各日期均比次优指数降低了16.43 %~38.79 %。这些发现表明SPSI在降低光谱饱和的价值，并且利用高分辨率卫星影像可以更好估计PNPA的巨大潜力。

图1 不同数据集在抽穗前个别日期的PNPA与SIs线性关系的决定系数(R² )值，无价值区域的值均小于0.00。在单独的( VE和SE数据集)和合并的数据集中显示了抽穗前不同日期的PNPA~SI关系的R²热图。对于任何一个数据集，9个传统SI中的4个较好的SI在所有日期上与PNPA的关系最密切，即NDRE、CI_red-edge、NEI和DATT_{[ 850，730，675]}，此外DATT_{[ 850，730，675]}在合并数据集上获得了最高的R²值。

图2 不同日期冬小麦PNPA与DATT_{[ 850，730，675]} ( A & C )和TI_NDRE[MEA]( B & D )的关系，黑色线条表示所有数据点的拟合线条。展示了PNPA与DATT_{［850，730，675]}和TINDRE_［MEA]在3月14日和抽穗前部分个体和混合数据集的关系，无论3月14日还是抽穗前期，当应用于合并数据集DATT_{［850，730，675]}时，各组个体数据集之间均存在明显的分层现象，尤其是高值范围( 图2A和C )，此外，从个别数据集的3月14日部分估算PNPA时，TI_NDRE［MEA] ( R²=0.66 ( VE )和0.77 ( SE )，图2B )的R²值与DATT_{［850，730，675]} ( R²= 0.68 ( VE )和0.76 ( SE )，图2A)的R²值相近。然而，TI_NDRE［MEA]在应用于混合数据集的3月14日或预抽穗部分时，表现出比DATT_{［850，730，675]} ( R² > 0.45，2A & C)更明显的分层( R²< 0.26，图2B & D)

图3 R²为PNPA与所有可能波段组合( A-D )的NDTICOR之间的线性关系，以及PNPA与NDTI_{COR[ 850，730 ]}( E )和NDTICOR_{[ 730，850 ]}( F )之间的关系，灰色网状区域代表无效组合，白色区域代表R²小于0.45的波段组合

图4 冬小麦PNPA与DATT_{[ 850，730，675]}、NDTI_{COR[ 850，730 ]}和SPSI在个别日期( A-I )或整个孕穗期( J-L )的关系，A-C为所有样本，D-F为DATT_{[ 850，730，675]}分别在0.35-0.45 ( D )、0.38-0.48 ( E )和0.41-0.51 ( F )范围内的子集，G-I是对于PNPA样本(点数)分别与D-F中相同的子集，利用DATT_{[ 850，730，675]}和NDTI_{COR[ 850，730 ]}推导SPSI

图5 5种情景下冬小麦PNPA与DATT_{[ 850，730，675]} ( A-E )和SPSI ( F-J )的关系。2021年3月14日，SCE # 1代表场景# 1，其中N水平和播期一致，但V9 ( A & F )的密度不同；SCE # 2代表情景# 2，其中N水平、播期、密度和叶型一致，但品种不同( B & G )；SCE # 3代表情景# 3，其中N水平、播期和叶型一致，但密度和品种不同( C & H )；SCE # 4代表情景# 4，其中N水平、播期和密度一致，但叶型和品种不同( D & I )；SCE # 5代表情景# 5，其中N水平和播期一致，但密度、叶型和品种不同( E & J)

图6 不同实验条件下PNPA与DATT_{[ 850，730，675]} ( A-D )和SPSI ( E-H )的关系，A ~ H分别为氮水平( A & E)、种植密度( B & F)、播期( C & G)和叶型( D & H)的数据分析，黑色线条表示所有数据点的拟合线条