基于无人机高光谱遥感的水稻氮营养诊断方法
许童羽 白驹驰 郭忠辉 金忠煜 于丰华
(1.沈阳农业大学信息与电气工程学院,沈阳 110866;
2.辽宁省农业信息化工程技术研究中心,沈阳 110866)
氮素是对作物生长、发育,以及最终产量影响最明显的元素,其含量变化会对光合作用、蛋白质合成以及碳氮代谢产生影响。氮素的缺乏会抑制作物地上部分和根系的生长,限制繁殖器官的形成和发育,并显著影响作物最终产量以及品质[1-2]。因此,快速、精准、大面积地对水稻田间氮素需求情况进行诊断,并依据诊断结果实现精准施肥,是实现水稻田间精准管理和保证水稻产量的重要手段。传统的水稻氮营养诊断标准主要基于临界氮浓度曲线(Critical N concentration curve)。GREENWOOD等[3]最早提出了临界氮浓度的概念,即作物达到最大生物量所需的最小氮浓度。随后LEMAIRE等[4]提出了基于临界氮浓度理论的氮营养指数(Nitrogen nutrition index,NNI)。氮营养指数作为作物氮营养诊断的重要指标,能够定量描述作物的氮营养丰缺状况[5-6],为作物氮亏缺量的确定提供了理论依据。吕茹洁等[7]基于临界氮浓度理论,计算出超级杂交稻和常规稻的临界氮浓度曲线,并依据曲线确定了二者的氮亏缺量与适宜施肥量,为水稻氮亏缺量的确定提供了理论支撑。
然而,基于氮营养指数计算氮亏缺量需要通过田间采样获取数据,成本较高、测量周期长,且结果具有滞后性[8],难以对实际农业生产做出指导。近年来,随着高光谱遥感技术的发展,利用无人机高光谱遥感技术获取水稻氮营养状况信息成为精准农业领域的重要发展方向。高光谱遥感技术能够通过分析作物光谱数据来获得其生长信息,具有快速、无损、准确等优点[9]。对此,国内外已在相关领域取得一定的研究进展[10-11]。这些研究大多数集中在利用高光谱对水稻氮含量进行反演,而单纯的氮含量信息难以反映水稻的氮营养丰缺状况,因此,结合无人机高光谱遥感技术与临界氮浓度曲线理论,通过光谱反演氮营养指数的形式对水稻进行氮营养诊断成为当前的研究热点[12-13]。ZHA等[14]通过无人机遥感获取水稻冠层光谱数据,构建多种植被指数并结合机器学习算法对水稻地上生物量、氮吸收量与氮营养指数进行反演建模,结果显示随机森林算法显著提高了氮营养指数反演精度;
QIU等[15]提取水稻RGB光谱信息并构建植被指数,最后比较多种机器学习方法反演水稻氮营养指数的精度,认为随机森林算法在各生育期表现最佳。
前人研究大都以对水稻进行氮营养诊断为主,而单纯的氮营养诊断只能判断水稻的生长状况,需要进一步获取水稻氮亏缺量数据才能以此对精准施肥做出定量指导。因此,本文拟以水稻无人机高光谱数据与氮亏缺量为研究对象,以水稻各个时期氮亏缺量约等于0的冠层光谱为标准,对水稻高光谱数据进行比值、差值与归一化差值变换,然后利用竞争性自适应重加权采样法处理原始光谱与变换光谱,筛选出效果最佳的光谱特征波段,并建立基于多元线性回归、极限学习机与蝙蝠算法优化极限学习机的水稻氮亏缺量反演模型,最后比较模型反演效果,确定基于高光谱数据的水稻氮亏缺量最佳反演方法,以期为快速获取水稻氮亏缺量提供研究方向,也为基于田间水稻氮营养状况的精准施肥提供定量指导。
1.1 实验设计
实验于2021年6—9月在辽宁省鞍山市海城市耿庄镇沈阳农业大学精准农业航空科研基地(北纬40°58′45.39″,东经122°43′47.01″)进行,为避免阴雨多云天气对遥感数据采集造成误差,数据采集过程选择天气晴好的日期,如遇云量超过20%或不利于遥感数据采集天气则顺延采集。水稻品种为北粳1705。实验区域如图1所示,实验田分为2个大区,其中实验区1设立5个氮肥梯度,施氮量分别为N0(0 kg/hm2)、N1(75 kg/hm2)、N2(150 kg/hm2)、N3(225 kg/hm2)、N4(300 kg/hm2),实验区2也设计5个氮肥梯度,施氮量分别为:N0(0 kg/hm2)、N1(50 kg/hm2)、N2(100 kg/hm2)、N3(150 kg/hm2)、N4(200 kg/hm2),氮肥基追比为5∶3∶2。实验区1对每个梯度设立3个重复小区,共设立3×5=15个小区。每个小区面积为5 m×8 m=40 m2,实验区2每个小区面积为660 m2,除氮肥梯度外两组实验区田间管理一致,磷钾肥施用量采用当地标准施用量进行,其中磷肥标准施用量为144 kg/hm2,钾肥标准施用量为192 kg/hm2,基追比为1∶1,其余田间管理同常规高产管理。田间采样由分蘖期至抽穗期,采样间隔9 d,每次采样在各实验小区选取具有代表性的3穴水稻进行冠层光谱、水稻氮浓度与地上干物质量的获取,最后结果取平均值作为该小区水稻氮浓度与地上干物质量。实验区1共采集小区数据120组,实验区2共采集小区数据88组,共计208组样本。
图1 实验小区分布图
1.2 数据获取
1.2.1水稻光谱参量获取
无人机高光谱平台采用深圳大疆创新公司M600 PRO型六旋翼无人机,高光谱成像仪选用四川双利合谱公司GaiaSky-mini内置推扫式机载高光谱成像系统,高光谱波段范围为400~1 000 nm,分辨率为3 nm,有效波段数为253个。无人机高光谱遥感平台数据采集时间为选择太阳光强较为稳定的时段11:00—12:00,无人机飞行高度为150 m。运用 ENVI 5.3+IDL工具软件对获取的高光谱遥感影像进行小区高光谱数据提取,运用波谱角填图方法去除干扰地物光谱的影响,对每个小区的感兴趣区计算平均光谱,再利用Matlab软件对平均光谱进行重采样处理,使光谱分辨率降为1 nm,最后利用高斯滤波器对重采样后光谱进行去噪处理,结果作为每个实验小区的高光谱信息。
1.2.2水稻农学参量获取
测定样本氮浓度与地上干物质量时,首先对每个小区的水稻进行破坏性采样,将样本带至实验室,然后将样本放入干燥箱中105℃杀青30 min后,80℃干燥至恒质量;
之后测量干燥后样本的干物质量,并根据种植密度换算成地上部干物质量;
最后将干燥后的样本磨碎,对研磨后的样本通过凯氏定氮法测量植株氮浓度。
1.3 基于临界氮浓度曲线的氮亏缺量计算方法
临界氮浓度是作物达到最大生物量所需的最小氮浓度,因此以临界氮浓度为基准可以推导出临界氮积累量方程,进而建立氮亏缺量方程,并确定水稻氮亏缺量[7]。
根据JUSTES等[16]提出的临界氮浓度曲线构建方法,以所测水稻氮浓度与地上干物质量为基础构建东北水稻的临界氮浓度曲线,具体方法如下:对比不同氮肥梯度下测定的氮浓度与干物质量,通过方差分析将样本分为受氮营养限制组与不受氮营养限制组;
对受氮营养限制组的数据,将地上干物质量与氮浓度进行线性拟合;
对不受氮营养限制组,对同一时期样本地上干物质量取平均值代表该时期最大地上干物质量;
取各时期拟合曲线于最大地上干物质中位置的截点为该时期的理论临界氮浓度点;
将各时期临界氮浓度点进行幂函数拟合,构建作物临界氮浓度曲线
Nc=aM-b
(1)
式中Nc——水稻临界氮质量比,g/g
M——地上干物质量,t/hm2
a、b——曲线参数
根据临界氮浓度曲线推导出临界氮积累量方程与氮亏缺量方程,推导过程为
Ncna=10aM1-b
(2)
Nand=Ncna-Nna
(3)
式中Ncna——临界氮浓度条件下植株氮积累量,kg/hm2
Nand——氮亏缺量,kg/hm2
Nna——不同施氮量下植株实际氮积累量,kg/hm2
1.4 光谱数据处理
1.4.1高光谱数据变换
由于临界氮浓度状态下作物地上部氮营养状况与冠层叶片结构达到最佳状态,而这些都是影响水稻冠层光谱反射率的直接因素[17-18],因此,相比非临界氮浓度状态下的作物,临界氮浓度下水稻同其余水稻无论在光谱层面还是氮营养状况层面都存在一定的差异。为了在光谱层面突出临界氮浓度下水稻与其余水稻的差距,本研究选取各个时期Nand≈0 kg/hm2样本对应的光谱反射率,将其与同时期其余样本光谱反射率分别做比值、差值和归一化差值变换,计算方法为
(4)
RD=RNC-Ri
(5)
(6)
式中RSV——比值光谱反射率,%
RD——差值光谱反射率,%
RND——归一化差值光谱反射率,%
RNC——临界氮浓度状态下样本光谱反射率,%
Ri——同时期其余光谱反射率,%
1.4.2高光谱特征波段提取
由于全波段光谱相邻波段数据相似度较高,且存在大量与所求变量无关的冗余信息,基于全波段光谱的建模往往存在运行速度较慢,模型误差较高,反演精度较低等缺点[19-21]。为了减少输入变量个数,降低数据冗余性,提高建模速度与精度,本文采用竞争性自适应重加权采样法(Competitive adaptive reweighted sampling,CARS)对光谱反射率数据进行特征提取。CARS算法基于达尔文进化论简单有效的生存原则“适者生存”,对不适应的波长变量进行逐步淘汰。首先,利用蒙特卡罗采样法(Monte Carlo simulation,MCS)采样N次,再以偏最小二乘回归(Partial least squares regression,PLSR)系数绝对值作为衡量标准,保留回归系数绝对值大的波长变量,去除回归系数绝对值小的波长变量,从而获得一系列的波长变量子集;
然后对比每次产生的PLSR模型的交互验证均方差值(Root mean square error of cross validation,RMSECV),RMSECV 最小的那个模型所对应的变量子集被选为最优变量子集。
1.5 反演建模方法
选用多元线性回归(Multivariable linear regression,MLR)、极限学习机(Extreme learning machine, ELM)与蝙蝠算法优化极限学习机(Bat algorithm optimized extreme learning machine, BA-ELM)3种算法进行建模,并通过3种模型对训练集与测试集预测结果的决定系数R2与均方根误差(Root mean square error,RMSE)来判断模型的反演精度与鲁棒性。
多元线性回归是一种解释单因变量与多自变量间线性关系的传统回归建模方法,模型结果唯一,并能够更好地体现各个自变量与因变量间的相关程度。
极限学习机是一种基于单层前馈神经网络理论的改进算法,具有学习速率快、泛化能力强、训练精度较高等优点[22-23]。但由于极限学习机每次运行中输入层与隐含层之间的连接权值和隐含层神经元的阈值都是随机生成的,其在训练稳定性上存在不足,且容易陷入局部最优解。因此,本研究采用蝙蝠算法对极限学习机隐含层初始权重进行优化。
蝙蝠算法是一种新的元启发式优化算法。该算法模拟自然界中蝙蝠利用回声定位来探测猎物和规避障碍物的原理,控制蝙蝠发出不同频率的搜索脉冲来寻找最优解。相较于经典的粒子群算法,蝙蝠算法引入了局部搜索机制,在每轮迭代中最优蝙蝠个体会采用随机搜索的方式进行局部寻优,有助于算法跳出局部最优解,增强算法的全局搜索性[24-26]。
2.1 植株氮亏缺量确定
2.1.1植株氮浓度与地上部干物质量统计分析
东北水稻全生育期植株氮浓度与地上部干物质量的基本信息如表1所示,由于实验包含全生育期采样数据,实验区1与实验区2中的植株氮浓度与地上部干物质量离散程度较高,最大值与最小值差别较大。其中实验区1由于施氮梯度较高,其整体氮浓度与地上部干物质量高于实验区2。实验区2由于在分蘖期与拔节期采样密度较高,其平均氮浓度较高,平均地上部干物质量较低。
表1 植株氮浓度与地上部生物量统计
2.1.2临界氮浓度曲线与氮亏缺量计算
根据1.3节提出的临界氮浓度曲线的构建方法,本文将每个采样日获得的氮浓度与干物质量进行回归拟合,计算出每个采样日的东北水稻临界氮浓度,之后根据各个临界氮浓度与对应的干物质量构建东北水稻临界氮浓度曲线(图2)。曲线方程为
图2 临界氮质量比拟合曲线
Nc=2.026M-0.460 3
(10)
决定系数R2为0.874 7,RMSE为0.382 5 g/g。曲线参数a、b分别为2.026、-0.460 3,与SONG等[27]计算出的东北水稻临界氮浓度曲线参数a=1.99、b=-0.44非常接近,可进一步用于东北水稻的氮亏缺量计算。
根据1.3节中的计算方法,本文结合样本数据与临界氮浓度曲线,计算各样本的氮亏缺量。各个氮肥梯度的平均氮亏缺量如图3所示,随着氮肥梯度的增高,氮亏缺量整体呈下降趋势,综合实验区1与实验区2氮亏缺量的下降趋势可知,水稻最佳施肥量在150~200 kg/hm2之间。
图3 各氮肥梯度平均氮亏缺量
2.2 光谱数据处理
2.2.1光谱数据变换结果
根据1.4.1节提出的方法,本研究选取各个时期Nand≈0 kg/hm2的样本对应光谱反射率,将其与其余样本光谱反射率分别作比值、差值、归一化差值变换,变换结果与相关性分析结果如图4、5所示,对比可得,不同施氮水平下原始光谱全波段反射率变换趋势基本一致,在可见光波段(波长400~700 nm)光谱反射率均呈“先升后降”的变化规律。同时,在绿波段(波长550 nm附近)出现明显的反射峰,在红波段(波长680 nm附近)出现明显的吸收谷。在波长680~770 nm内不同施氮水平的光谱反射率曲线变化基本一致,反射率急剧升高,在近红外波段(770~1 000 nm)形成较高的反射平台;
不同施氮水平下比值、差值与归一化差值光谱反射率有一定差别,但在绿波段与红波段差别较为明显;在与氮亏缺量的相关性方面,4种光谱反射率与氮亏缺量的相关性图像趋势基本一致,在红波段与绿波段形成2个波峰,其中归一化差值光谱在2个波峰处相关性最佳,其次分别是比值光谱、差值光谱与原始光谱。
图4 原始与变换光谱反射率
2.2.2特征波长选择
利用CARS算法对水稻原始光谱与归一化差值光谱反射率进行特征波长的筛选,根据交叉验证集RMSECV确定光谱特征波长后,删去相邻的特征波长,并通过相关性分析选取与Nand相关性最好的波长作为最佳特征波长,以原始光谱选择过程为例:由图6a可以看出,随着运行次数的增加,原始光谱被选出的波段数逐渐减少,图6b中RMSECV整体呈下降趋势,说明筛选过程中剔除的变量与SOM去除量无关,而72次迭代以后,RMSECV呈回升趋势,表明反射率光谱中与SOM无关的大量信息或噪声被添加,从而导致RMSECV上升。图6c为所有变量在每次采样过程中的回归系数路径图,表示随着运行次数的增加各波段变量回归系数的变化趋势。结合图6b分析发现当运行次数为第72次时,RMSECV最小即所选择的光谱变量子集最优。原始光谱与3种变换光谱的特征选择结果如表2所示。
图5 原始与变换光谱相关性
图6 原始光谱CARS降维结果
表2 CARS提取后原始光谱与变换后光谱特征波长与运行结果
2.3 氮亏缺量反演模型构建
为比较不同光谱变换方法与反演建模方法的性能,分别以经CARS算法筛选出的原始光谱、比值光谱、差值光谱、归一化差值光谱特征波长为输入变量,水稻氮亏缺量为输出变量,分别构建基于MLR、ELM与BA-ELM的水稻氮亏缺量反演模型。
2.3.1基于MLR的氮亏缺量反演
将4种光谱特征波长作为输入变量,分别构建基于MLR的反演模型,反演结果如图7所示:其中基于归一化差值光谱的模型精度最高,训练集与测试集的R2分别为0.654 0、0.715 8,RMSE分别为0.955 7、1.204 8 kg/hm2;
其次是基于差值光谱反演模型,其R2分别为0.600 5、0.524 2,RMSE分别为1.103 1、1.184 2 kg/hm2;
然后是基于原始光谱的反演模型,R2分别为0.595 6、0.582 6,RMSE分别为1.189 4、1.083 9 kg/hm2。基于比值光谱的模型精度精度R2分别为0.597 3、0.517 5,RMSE分别为1.034 0、1.064 6 kg/hm2;
基于归一化差值光谱的反演模型反演效果最好。
图7 基于MLR的不同反演模型结果
2.3.2基于ELM的氮亏缺量反演
将4种光谱特征波长作为输入变量,分别构建基于ELM的反演模型,映射函数为Sigmoid,隐含层个数分别为13、14、12、13,训练结果如图8所示。其中基于归一化差值光谱的模型精度最高,训练集与测试集的R2分别为0.721 1、0.683 2,RMSE分别为0.857 8、1.151 2 kg/hm2;
其次是基于比值光谱的反演模型,R2分别为0.685 0、0.634 4,RMSE分别为0.914 5、0.958 4 kg/hm2;
然后是基于差值光谱反演模型,其R2分别为0.632 2、0599 4,RMSE分别为1.057 0、1.069 9 kg/hm2;
最后是原始光谱的模型精度最低,R2分别为0.619 2、0.547 7,RMSE分别为1.154 6、1.175 3 kg/hm2。基于归一化差值光谱的反演模型R2与RMSE均优于其余3种模型。
图8 基于ELM的反演模型结果
2.3.3基于BA-ELM的氮亏缺量反演
将4种光谱特征波长作为输入变量,分别构建基于BA-ELM的反演模型,其中蝙蝠算法的最大迭代次数为80,种群大小为50,适应度函数为验证集RMSE,隐含层个数分别为50、60、57、65,结果如图9所示:基于归一化差值光谱的模型精度最高,训练集与验证集的R2分别为0.816 0、0.830 6,RMSE为0.696 8、0.814 1 kg/hm2;
基于比值光谱的模型精度R2分别为0.759 3、0.770 1,RMSE分别为0.799 4、0.730 9 kg/hm2;
基于差值光谱的模型精度R2分别为0.733 9、0.659 0,RMSE分别为0.900 3、0.997 1 kg/hm2;
原始光谱的模型精度,训练集与验证集的R2分别为0.692 2、0.680 8,RMSE为1.037 5、0.952 8 kg/hm2。经过模型结果对比可知:基于BA-ELM的反演模型精度均优于基于MLR与ELM的反演模型;
其中基于归一化差值光谱的反演模型在R2与RMSE上均优于其余3种光谱模型;
基于比值光谱在MLR模型中表现较差,但在ELM与BA-ELM模型中表现仅次于归一化差值光谱,这可能是因为比值光谱与氮亏缺量之间存在较为复杂的非线性关系,MLR等线性模型对它的解释能力较差。在3种模型中,差值光谱与原始光谱表现相近,可见对光谱进行单纯的差值处理难以突出其与氮亏缺量的关系。因此采用基于归一化差值光谱的BA-ELM模型在反演水稻氮亏缺量上最有优势。
图9 基于BA-ELM的反演模型结果
2.4 针对不同施肥量的模型预测能力
在实验区1与实验区2中,以反演效果最优模型的氮亏缺量估测值与对应梯度的氮亏缺量实测值对比,结果如图10所示,实验区1与实验区2中,整体上估测值与实测值较为相近,对水稻氮亏缺量的估测能力较强;
随着施氮量梯度的提升,对氮亏缺量的估测精度呈下降趋势,在临界氮浓度状态附近估测能力较差,这可能是因为临界氮浓度状态下与非临界氮浓度状态下水稻光谱变化规律不一致,模型同时估测2种状态的水稻光谱预测难度较高所致。
图10 不同梯度氮亏缺量估测效果
以往学者研究焦点往往集中在如何提高获取水稻田间氮浓度或氮营养指数信息的准确性[28-29],然而这些信息只能描述水稻的氮营养状况,难以对水稻氮亏缺状况进行精确的定量描述,进而指导精准施肥。本文首先通过田间采样获取水稻氮浓度、地上生物量与冠层光谱数据,根据前人提出的方法构建东北水稻临界氮浓度曲线,并基于临界氮浓度曲线计算水稻氮亏缺量。由于临界氮浓度表示作物达到最大干物质量所需的最小氮浓度,以临界氮浓度作为判断水稻氮亏缺量的标准更能反映田间水稻的氮营养状况与实际氮需求情况。同时本研究还采用了无人机高光谱遥感技术与数据驱动的反演建模方法,通过水稻冠层光谱数据对氮亏缺量进行反演。无人机高光谱遥感技术能够快速、大面积获取水稻冠层光谱信息,进而对田间水稻氮营养状况进行快速诊断,数据驱动的反演建模方法虽然可解释性与泛化能力较差,但其能够对光谱与氮亏缺量之间的复杂非线性关系进行准确描述,将二者结合能够同时满足追肥决策所需的精确性与时效性,为定量追肥决策的实现提供了解决思路。
由于临界氮浓度状态下作物地上部氮营养状况与冠层叶片结构达到最佳状态,而这些都是影响水稻冠层光谱反射率的直接因素[30],因此,相比非临界氮浓度状态下的作物,临界氮浓度下作物无论在光谱层面还是氮营养状况层面都存在一定差异。为了在光谱层面上放大临界氮浓度下水稻与其余水稻的差距,去除与氮亏缺量无关的冗余信息,本研究以各时期接近临界氮浓度状态下光谱为标准光谱,分别对原始光谱进行比值、差值、归一化差值变换,之后运用CARS算法提取变换光谱的特征波长。由结果可知3种光谱变换均提升了其与氮亏缺量之间的相关性,其中归一化差值变换相关性与反演精度最高,这可能是由于归一化差值计算既通过差值计算放大了临界氮浓度光谱与其余冠层光谱的差异,又统一了差值光谱数据的量纲,还对光谱数据的统计分布性进行了归纳。
在建模过程中,本文选择蝙蝠算法优化ELM算法的初始权重,并将其与原始ELM进行对比建模。建模结果中BA-ELM模型的反演精度显著高于ELM模型,这可能是由于ELM作为一种基于前馈神经网络的机器学习算法,具有学习速度快,收敛能力强等优点,但同时也存在容易陷入局部最优解的问题;
蝙蝠算法作为一种由粒子群算法衍生的元启发式搜索算法,在迭代过程中会在最优解附近再次搜索局部最优解,因此具有优秀的全局搜索能力,通过蝙蝠算法优化极限学习机的初始参数会让极限学习机在最优解附近开始学习,提高了算法的反演精度与全局搜索能力。
虽然本研究利用无人机高光谱数据对水稻氮亏缺量反演取得了较好的效果[31],但对于实现基于无人机高光谱遥感的水稻精准追肥仍存在以下几点问题:①本文采用的归一化差值处理方法是基于各个采样期临界氮浓度状态的水稻光谱计算的,而这些光谱如果用于其他水稻光谱的计算时,可能会因为采样时间、水稻品种与发育状况等因素造成误差,进而影响氮亏缺量反演效果。②蝙蝠算法优化后的ELM算法也存在收敛精度不高,容易陷入局部最小值等问题。后续的研究重点应集中在如何在无人机高光谱数据中实时提取临界氮浓度状态下的水稻光谱与对反演算法的精度与鲁棒性进行优化这两方面上。
(1)构建了东北水稻临界氮浓度曲线,可表示为Nc=2.026M-0.460 3。
(2)通过CARS算法提取的原始光谱、比值光谱、差值光谱与归一化差值光谱特征波长分别为400、411、422、447、570、677 nm,533、695、715、738、774 nm,413、424、444、497、533、600、676、691 nm,484、560、644、672、695、724 nm。
(3)对比基于4种光谱特征波长与ELM、BA-ELM算法的8种建模方法反演精度可知,基于归一化差值的BA-ELM模型反演结果最好,训练集与测试集的R2分别为0.816 0、0.830 6,RMSE为0.696 8、0.814 1 kg/hm2。
猜你喜欢物质量实验区反射率影响Mini LED板油墨层反射率的因素印制电路信息(2022年11期)2022-11-30近岸水体异源遥感反射率产品的融合方法研究海洋通报(2022年4期)2022-10-10具有颜色恒常性的光谱反射率重建光谱学与光谱分析(2022年4期)2022-04-06水库工程区水土保持生态服务价值估算人民黄河(2020年12期)2020-12-30施肥量对谷子干物质量积累及分配的影响山西农业科学(2020年9期)2020-09-14不同播期与品种对糯玉米干物质积累的影响山西农业科学(2019年12期)2019-12-19青海省人民政府办公厅关于加强文化生态保护实验区建设的指导意见青海政报(2017年18期)2018-01-31一种柱状金属物质量检测器的研究电子制作(2017年10期)2017-04-182016年国家文创实验区规上文化产业收入近2000亿元投资北京(2017年1期)2017-02-13足球应用型人才培养模式创新实验区的探索与实践——以学生社会实践为突破口陕西教育·高教版(2015年4期)2015-06-28- 范文大全
- 说说大全
- 学习资料
- 语录
- 生肖
- 解梦
- 十二星座
-
主题党日活动交流发言8篇
主题党日活动交流发言8篇主题党日活动交流发言篇13月13日,东城区党史学习教育动员大会召开。市委
【活动总结】 日期:2022-12-23
-
2022年4月主题党日活动记录范文15篇
2022年4月主题党日活动记录范文15篇2022年4月主题党日活动记录范文篇1一个崇尚阅读的民族,必然精神饱满、意气风发、活力四射。习近平总书记强调:“学习
【活动总结】 日期:2022-08-01
-
家乡赋|最美的家乡赋
家乡赋 孙传志 今安康市,白河双丰镇,吾之家乡也。三环沃土,山水环抱。其北依山,山系五岭,山
【调研报告】 日期:2020-04-01
-
少先队的光荣历史故事 队前教育-光辉历程
2017-2018学年队前教育1光辉历程一、劳动童子团1924——1927二、三十年代年的中国是一个
【法律文书】 日期:2020-06-23
-
【人教版1-6年级数学上册知识点精编】1-6年级数学人教版教材
人教版二年级数学上册知识点汇总第一单元长度单位一、米和厘米1、测量物体的长度时,要用统一的标准去测量
【调研报告】 日期:2020-11-08
-
党支部1-12月全年主题党日活动计划表
2022年党支部主题党日活动计划表序号活动时间活动方式活动内容12022年1月专题学习研讨集中观看2022年新年贺词,积极开展学习研讨交流。组织生活会组织党员认真对照党章...
【活动总结】 日期:2022-10-14
-
2022年2月份主题党日活动记录5篇
2022年2月份主题党日活动记录5篇2022年2月份主题党日活动记录篇1尊敬的党组织:在今年的开学初,本人积极参加教研室组织的教研活动,在学校教研员的指
【活动总结】 日期:2022-08-12
-
2023年平安校园建设方案13篇
平安校园建设方案“平安校园”创建工作,我们幼儿园全体教职员工一直把它当作头等大事来抓。领导高度重视,以“平安校园”创建活动为抓手,建立和规范校园安全工作机制
【规章制度】 日期:2023-11-02
-
医院最佳主题党日活动11篇
医院最佳主题党日活动11篇医院最佳主题党日活动篇1 医院最佳主题党日活动篇2为隆重纪念中国共产党成立100周年,进一步巩固党的群众路线教育实践活动成果,切实
【活动总结】 日期:2022-10-29
-
主题党日活动记录202210篇
主题党日活动记录202210篇主题党日活动记录2022篇12021年是中国共产党成立100周年,为广泛开展爱国主义宣传教育,铭记党的历史,讴歌党的光辉历程,
【活动总结】 日期:2022-08-02
-
正式的晚宴邀请函 公司晚宴邀请函
尊敬的先生 女士: 我公司谨定于xxxx年xx月xx日xx:xx在xxxx店隆重举行xx市xx届xxxx晚宴(宴会地址:xx区xx路xxxx) 敬请届时光临!xxxxxx集团股份有限公司xxxx有限公司敬邀xxxx年xx月xx日
【简历资料】 日期:2019-08-03
-
一年级新学期目标简短_一年级学生新学期打算
新学期到了,我是一年级下册的小学生了。 上课的时候,我要认真学习,不做小动作,认真听讲。我要认真学习,天天向上,努力学习,耳朵要听老师讲课,眼睛要瞪得大大的看老...
【简历资料】 日期:2019-10-26
-
[信访复查复核制度作用探讨]信访复查复核有用吗
作为我国特有的一项制度,信访制度的出现并长期存在不是偶然的,虽然一些法学专家认为信访制度具有“人治”
【职场指南】 日期:2020-02-16
-
[党员干部2019年主题教育个人问题检视清单及整改措施2篇] 党员干部
2019年主题教育问题检视清单及整改措施根据主题教育领导小组办公室《关于认真做好主题教育检视问题整改
【求职简历】 日期:2019-11-08
-
网络维护工作内容_(精华)国家开放大学电大专科《网络系统管理与维护》形考任务1答案
国家开放大学电大专科《网络系统管理与维护》形考任务1答案形考任务1理解上网行为管理软件的功能【实训目
【职场指南】 日期:2020-07-17
-
党委会与局长办公会的区别_局长办公会制度
为进一步加强xxx局工作的规范化、制度化建设,提高行政效能,规范议事程序,特制定本制度。一、会议形式1、局长办公会议由局长、副局长参加。由局长召集和主持。根据工作需要...
【求职简历】 日期:2019-07-30
-
如何凝心聚力谋发展【坚定信心谋发展凝心聚力促跨越】
当前,清河正处于在苏北实现赶超跨越基础上全面腾飞的战略机遇期,处于在全市率先实现全面小康基础上率先实
【简历资料】 日期:2020-03-17
-
学生会组织部部长竞选稿5篇
学生会组织部部长竞选稿以“三制”为统领推进农村党的建设中共**市委组织部近年来,**市认真落实中央、省和徐州市委的部署,积极适应发展要求,从加强领导体制、运
【求职简历】 日期:2023-11-06
-
《铁拳砸碎“黑警伞”》警示教育片观后感
影片深刻剖析了广西北海市公安局海西派出所原所长张枭杰蜕变堕落的轨迹。观看警示教育片后,做为一名党员教
【简历资料】 日期:2020-08-17
-
系领导班子行为准则_必须以什么为行为准则
为了更好地完成本系的各项工作,充分发挥系领导班子的模范带头作用,特制定本行为准则如下:一、要勤政为民,真抓实干,坚持科学态度和求实的精神,兢兢业业地作好本职工作。...
【求职简历】 日期:2019-08-30
-
在纪念三八国际劳动妇女节大会上的讲话 (1)【优秀范文】
在2022年纪念三八国际劳动妇女节大会上的讲话同志们,妇女同胞们:春到三月季,又是桃杏斗艳时。今天,正值“三八”国际劳动妇女节XX周年,市妇联举办这次纪念大会很有意义,...
【其他范文】 日期:2022-10-31
-
公司2022-2023年党员教育培训工作实施情况中期评估报告
下面是小编为大家整理的公司2019-2023年党员教育培训工作实施情况中期评
【其他范文】 日期:2022-09-24
-
组织工作抓落实座谈会发言提纲(2023年)
今年年初召开的全国组织部长会议,坚持以习近平新时代中国特色社会主义思想为指导,全面贯彻党的十九大和十九届历次全会精神,对今年的组织工作作出部署。任务要求已经明确,...
【其他范文】 日期:2023-08-26
-
防汛工作汇报材料汇编【优秀范文】
2020年防汛工作汇报X县地处X平原,X市西南x公里。属X山前麓坡水区。境内地势平坦,自西北向东南倾斜。地面高程x~x米之间,地面坡降X铁路线以西X X,以东为X X左右,X铁路、...
【其他范文】 日期:2022-09-06
-
高中体育教师个人总结2018【体育教师个人工作总结】
导读:个人工作总结范文,体育教师个人工作总结在线阅读,教你怎么写,格式什么样,科教范文网提供各种参考范例: 从教以来,我坚持不懈的进行政治和专业学习,认真学习党...
【礼仪】 日期:2019-11-10
-
景观提升工程预算执行情况审计整改落实报告
市审计局:针对贵局《平湖市审计局审计报告》(平审投报[2019]24号)我司已收悉。公司高度重视,提
【其他范文】 日期:2020-11-10
-
中国省域畜牧业绿色发展竞争力时空演变特征
张朝辉,姬一帆,何泽军(河南农业大学经济与管理学院,郑州450046)中国畜牧业已经进入以环保为重点
【其他范文】 日期:2023-01-18
-
在示范区2021年度工作大会上的讲话.docx
在示范区2021年度工作大会上的讲话x月x日,市委常委、市政府常务副市长x出席x示范区x年度工作大会
【口号大全】 日期:2021-04-27
-
[加强学校综合治理 平安创建工作计划] 学校综合治理工作计划
加强学校综合治理平安创建工作计划为扎实推进“平安校园”创建活动,切实解决我校当前影响安全稳定的突出问题,积极推动学校教育改革和发展,努力创造稳定和谐的社会环境和公...
【口号大全】 日期:2019-10-31
-
电气自动化技术在电气工程中的实际应用
宋晓楠(山西省安装集团股份有限公司,山西太原030032)电气自动化技术的科学运用可以进一步提高生产
【其他范文】 日期:2023-07-02
-
军转座谈会交流发言4篇
军转座谈会交流发言4篇军转座谈会交流发言篇1大家好,我叫贺丽,2015届选调生,来自康定市委组织部,现在省委编办跟班学习。今天,非常荣幸向大家汇报我的学习收
【发言稿】 日期:2022-10-27
-
12岁生日小寿星发言4篇
12岁生日小寿星发言4篇12岁生日小寿星发言篇1各位来宾、各位朋友:大家好!今天,我们欢聚在这里,共同庆祝**十二周岁生日。首先,我代表**的父母以
【发言稿】 日期:2022-07-31
-
廉政大会总结发言稿7篇
廉政大会总结发言稿7篇廉政大会总结发言稿篇1各位领导,同志们:根据会议安排,我就党风廉政建设工作做表态发言,不妥之处,请批评指正。一、提高认识,切实
【发言稿】 日期:2022-10-30
-
被约谈的表态发言8篇
被约谈的表态发言8篇被约谈的表态发言篇1各位领导、各位党员大家好:这天我能站在鲜红的党旗下,
【发言稿】 日期:2022-12-24
-
破冰提能大讨论个人发言4篇
破冰提能大讨论个人发言4篇破冰提能大讨论个人发言篇1党史学习教育开展以来,我坚持读原著、学原文、悟原理。今天,根据会议安排,现在我就“学史明理”主题谈几点个
【发言稿】 日期:2022-10-09
-
我最敬佩的人开头_我敬佩的一个人作文20篇2020年
我敬佩的一个人作文20篇 我敬佩的一个人作文一): 我身边有很多值得我们敬佩的人,但我最敬佩的一
【发言稿】 日期:2020-11-10
-
党员教育培训总结交流发言12篇
党员教育培训总结交流发言12篇党员教育培训总结交流发言篇1根据市委组织部《关于开展我市〈20XX
【发言稿】 日期:2022-12-19
-
纪委书记工作表态发言4篇
纪委书记工作表态发言4篇纪委书记工作表态发言篇1在镇党委政府正确领导下,在全村干部和群众的共同努力下,紧紧围绕建设社会主义新农村工作为重点,尽职尽责,与时俱
【发言稿】 日期:2022-09-30
-
[钻井队队长(副队长、指导员)岗位HSE应知应会试题(1863)]
钻井队队长(副队长、指导员)岗位HSE应知应会试题(判断题:771;单选题:626;多选题:466)
【贺词】 日期:2020-09-23
-
话剧《家》剧本 话剧剧本:爱的空间
找文章到更多原创-(http: www damishu cn)人物介绍:刘伟,男,32岁,某购物广
【演讲稿】 日期:2020-01-21
-
学习周永开先进事迹心得体会3篇
学习周永开先进事迹心得体会【一】通过学习周永开老先生先进事迹后,结合自己工作思考,感慨万千。同样作为
【格言】 日期:2021-04-10
-
2023年中国行政区划调整方案(设想优秀3篇
中国行政区划调整方案(设想优秀民政部第二次行政区划研讨会会议内容一、缩省的意义与原则1.意义1)利于减少中间层次中国行政区划层级之多为世界之最,既使管理成本
【周公解梦】 日期:2024-02-20
-
XX老干局推进党建与业务深度融合发展工作情况调研报告:党建调研报告
XX老干局推进党建与业务深度融合 发展工作情况的调研报告 党建工作与业务工作融合发展始终是一个充满生
【成语大全】 日期:2020-08-28
-
中国共产党第三代中央领导集体的卓越贡献
中国共产党第三代中央领导集体的卓越贡献 --------------继往开来铸就辉煌 【摘要】改
【成语大全】 日期:2020-03-20
-
信息技术2.0能力点 [全国中小学教师信息技术应用能力提升工程试题题库及参考答案「精编」]
全国中小学教师信息技术应用能力提升工程试题题库及答案(复习资料)一、判断题题库(A为正确,B为错误)
【格言】 日期:2020-11-17
-
党建工作运行机制内容有哪些_构建基层党建工作运行机制探讨
党的基层组织是党在社会基层组织中的战斗堡垒,是党的全部工作和战斗力的基础。加强和改进县级以下各类党的
【经典阅读】 日期:2020-01-22
-
集合推理_七,推理与集合
七推理与集合1 期中考试数学成绩出来了,三个好朋友分别考了88分,92分,95分。他们分别考了多少分
【名人名言】 日期:2020-12-18
-
电大现代教育原理_最新国家开放大学电大《现代教育原理》形考任务2试题及答案
最新国家开放大学电大《现代教育原理》形考任务2试题及答案形考任务二一、多项选择题(共17道试题,共3
【成语大全】 日期:2020-07-20
-
2023年和儿媳妇在一起幸福的句子3篇
和儿媳妇在一起幸福的句子1、假如人生不曾相遇,我还是那个我,偶尔做做梦,然后,开始日复一日的奔波,淹没在这喧嚣的城市里。我不会了解,这个世界还有这样的一个你
【格言】 日期:2023-11-10
-
基层党务工作基本内容_党建基本工作有哪些
党建基本工作有哪些(一) 基层党建工作包括哪些内容 选择了大学生村官这条路,你就与农村基层党
【名人名言】 日期:2020-08-06
-
关于三农工作重要论述心得体会3篇
关于三农工作重要论述心得体会3篇关于三农工作重要论述心得体会篇1习近平总书记指出:“建设现代化国家离不开农业农村现代化,要继续巩固脱贫攻坚成果,扎实推进乡村
【学习心得体会】 日期:2022-10-29
-
【福生庄隧道坍塌处理方案】 福生庄隧道在哪里
(呼和浩特铁路局大包电气化改造工程指挥部,内蒙古呼和浩特010050)摘要:文章介绍了福生庄隧道
【学习心得体会】 日期:2020-03-05
-
五个一百工程阅读心得体会13篇
五个一百工程阅读心得体会13篇五个一百工程阅读心得体会篇1凡益之道,与时偕行。在全国网络安全和信
【学习心得体会】 日期:2022-12-07
-
城管系统警示教育心得体会9篇
城管系统警示教育心得体会9篇城管系统警示教育心得体会篇1各党支部要召开多种形式的庆七一座谈会,组织广大党员进行座谈,回顾党的光辉历程,畅谈党的丰功伟绩,
【学习心得体会】 日期:2022-10-09
-
发展对象培训主要内容10篇
发展对象培训主要内容10篇发展对象培训主要内容篇1怀着无比激动的心情,我有幸参加了__新区区委党校20__年第四期(区级机关)党员发展对象培训班。这次的学习
【培训心得体会】 日期:2022-09-24
-
凝聚三种力量发展全过程人民民主心得体会12篇
凝聚三种力量发展全过程人民民主心得体会12篇凝聚三种力量发展全过程人民民主心得体会篇1新民主主义革命是指在帝国主义和无产阶级革命时代,殖民地半殖民地国家中的
【学习心得体会】 日期:2022-08-31
-
2022年全国检察长会议心得7篇
2022年全国检察长会议心得7篇2022年全国检察长会议心得篇1眼睛是心灵上的窗户,我们通过眼睛才能看到世间万物,才能看到眼前这美好的一切。拥有一双明亮的眼
【学习心得体会】 日期:2022-10-31
-
在街道深化作风建设推动高质量发展走在前列动员会上讲话
在2023年街道深化作风建设推动高质量发展走在前列动员会上的讲话同志们:今天我们召开“街道深化作风建设推动高质量发展走在前列动员会”,这次会议是街道三季度召开的第一场...
【军训心得体会】 日期:2024-03-17
-
全面从严治党的心得体会800字7篇
全面从严治党的心得体会800字7篇全面从严治党的心得体会800字篇1中国特色社会主义是我们党领导
【学习心得体会】 日期:2022-12-14
-
矫正心得体会6篇
矫正心得体会6篇矫正心得体会篇1今天,是自己出监后第一次参加阳光中途之家组织的社区矫正方面的教育
【学习心得体会】 日期:2022-12-24
-
2024年主题教育民主生活会批评与自我批评意见(38条)(范文推荐)
2023年主题教育民主生活会六个方面个人检视、相互批评意见:1 理论学习系统性不强。学习习近平新时代中国特色社会主义思想不深不透,泛泛而学的时候多,深学细照的时候少,特...
【邓小平理论】 日期:2024-03-19
-
2024年交流发言:强化思想理论武装,增强奋进力量(完整)
习近平总书记指出:“一个民族要走在时代前列,就一刻不能没有理论思维,一刻不能没有思想指引。”党的十八大以来,伴随着新时代中国特色社会主义思想在实践中形成发展的历程...
【三个代表】 日期:2024-03-19
-
2024年度镇年度县乡人大代表述职评议活动总结
xx镇20xx年县乡人大代表述职评议活动总结为响应县级人大常委会关于开展县乡两级人大代表述职评议活动,进一步激发代表履职活力,加强代表与人民群众的联系,提高依法履职水平...
【马克思主义】 日期:2024-03-19
-
“千万工程”经验学习体会(研讨材料)
“千万工程”是总书记在浙江工作时亲自谋划、亲自部署、亲自推动的一项重大决策,也是习近平新时代中国特色社会主义思想在之江大地的生动实践。20年来,“千万工程”先后经历...
【三个代表】 日期:2024-03-19
-
2024年在市政协机关工作总结会议上讲话
同志们:刚才,XX同志对市政协机关20XX年工作进行了很好的总结,很精炼,很到位,可以感受到去年机关工作确实可圈可点。XX同志宣读了表彰决定,机关优秀人员代表、先进集体代...
【邓小平理论】 日期:2024-03-18
-
在全区防汛防涝动员暨河长制工作推进会上讲话提纲【完整版】
区长,各位领导,同志们:汛期已经来临,我区城区防涝工作面临强大考验,形势不容乐观。年初,区城区防涝排渍指挥部已经召开专题调度会,修订完善应急预案,建立网格化管理机...
【马克思主义】 日期:2024-03-18
-
2024年镇作风整治工作实施方案(完整文档)
XX镇作风整治工作实施方案为深入贯彻落实党的二十大精神及省市区委深化作风建设的最新要求,突出重点推进干部效能提升,坚持不懈推动作风整治工作纵深发展,根据《关于印发《2...
【毛泽东思想】 日期:2024-03-18
-
2024市优化法治化营商环境规范涉企行政执法实施方案【优秀范文】
xx市优化法治化营商环境规范涉企行政执法实施方案为持续优化法治化营商环境,激发市场主体活力和社会创造力,规范行政执法行为,创新行政执法方式,提升行政执法质效,着力解...
【毛泽东思想】 日期:2024-03-18
-
2024年度关于开展新一轮思想状况摸底排查工作通知(完整)
关于开展新一轮思想状况摸底排查工作的通知为深入贯彻落实关于各地开展干部职工思想状况大摸底大排查情况上的批示要求和改革教育第二次调度会议精神,有针对性做好队伍教育管...
【三个代表】 日期:2024-03-18
-
2024年公路养护中心主任典型事迹材料(完整文档)
“中心的工作就是心中的事业”——公路养护中心主任典型事迹材料**,男,1976年6月出生,1993年参加工作,2000年4月调入**区交通运输局工作,大学本科学历,中共党员,现任**...
【马克思主义】 日期:2024-03-17