发明内容
针对上述问题,本发明的目的是提供一种能够对作物群体内光合有效辐射进行精确定位定量分析的作物群体光能空间分布量化方法。
为实现上述目的,本发明采取以下技术方案:一种作物群体光能空间分布量化方法,其包括以下步骤:1)采用空间网格取样方法选择测试作物群体,并在测试作物群体冠层的纵向测试截面区域内均匀设置测试点,采用PAR测试仪器对每一测试点的PAR进行测试,每一测试点所测得的PAR包括测试点的透射PAR和反射PAR;2)根据步骤1)测量得到每个测试点的PAR值,采用插值法求解出纵向测试截面区域内其它待测点的PAR值;3)在纵向测试截面区域内建立坐标系,对位于纵向测试截面区域内的测试点和待测点的位置进行编号生成Grid文件,Grid文件的每个元素Grid(i,j)记录每个测试点或待测点在纵向测试截面区域中的位置,i表示测试点或待测点在纵向测试截面区域中的横坐标,j表示测试点或待测点在纵向测试截面区域中的纵坐标;4)获取透射和反射的插值Grid文件后,利用矩阵运算方法求取Grid文件中每一元素所对应的光能截获率;5)基于求取的纵向测试截面区域中测试点和待测点的光能截获率,对作物群体在整个测试空间区域内的任意点PAR数据进行空间插值生成空间格网文件,并采用多重积分法计算整个测试空间区域的体积,然后除以作物冠层的纵向测试截面区域的面积,得到群体整个光能利用效率。
所述纵向测试截面区域内设置一测试框架,所述测试框架采用矩形框架,所述矩形框架的宽度为作物的行间距离,所述矩形框架的高度根据测试需要进行设定;所述矩形框架上横向间隔平行设置有若干横向测试杆,每相邻两个横向测试杆之间的间距相同,测试点根据测试需要均匀分布在每一所述横向测试杆上。
所述步骤2)的插值方法采用克里金插值方法。
所述步骤4)中的Grid光能截获率(i,j)=1-Grid反射PAR(i,j)-Grid透射PAR(i,j),Grid文件每一元素的光能截获率等于1减去该元素所对应的反射PAR再减去该元素所对应的透射PAR。
一种作物群体光能空间分布量化方法,其包括以下步骤:1)采用空间网格取样方法选择测试作物群体,并在测试作物群体冠层的纵向测试截面区域内均匀设置测试点,采用PAR测试仪器对每一测试点的PAR进行测试,每一测试点所测得的PAR包括测试点的透射PAR和反射PAR;2)计算每一测试点所对应的PAR透射率Tr和反射率Rr;
3)根据步骤2)计算得到每个测试点的PAR透射率Tr和反射率Rr,采用插值法求解出纵向测试截面区域内其它待测点的PAR透射率Tr和反射率Rr;4)在纵向测试截面区域内建立坐标系,对位于纵向测试截面区域内的测试点和待测点的位置进行编号生成Grid文件,Grid文件的每个元素Grid(i,j)记录每个测试点或待测点在纵向测试截面区域中的位置,i表示测试点或待测点在纵向测试截面区域中的横坐标,j表示测试点或待测点在纵向测试截面区域中的纵坐标;5)获取透射和反射的插值Grid文件后,利用矩阵运算方法求取Grid文件中每一元素所对应的光能截获率;6)基于求取的纵向测试截面区域中测试点和待测点的光能截获率,对作物群体在整个测试空间区域内的任意点PAR数据进行空间插值生成空间格网文件,并采用多重积分法计算整个测试空间区域的体积,然后除以作物冠层的纵向测试截面区域的面积,得到群体整个光能利用效率。
所述纵向测试截面区域内设置一测试框架,所述测试框架采用矩形框架,所述矩形框架的宽度为作物的行间距离,所述矩形框架的高度根据测试需要进行设定;所述矩形框架上横向间隔平行设置有若干横向测试杆,每相邻两个横向测试杆之间的间距相同,测试点根据测试需要均匀分布在每一所述横向测试杆上。
所述步骤3)的插值方法采用克里金插值方法。
所述步骤5)中的Grid光能截获率(i,j)=1-Grid反射率Rr(i,j)-Grid透射率Tr(i,j),Grid文件每一元素的光能截获率等于1减去该元素所对应的PAR反射率再减去该元素所对应的PAR透射率,即:In=1-PARt/PARI-PARr/PARI。
本发明由于采取以上技术方案,其具有以下优点:1、作物群体特别是行向种植的作物群体具有高度空间相关性,传统统计学的随机方法取样测定分析方法无法准确量化作物群体空间任意点的光能分布,本发明利用空间统计学的方法原理,以区域化变量为基础,通过空间网格采样方法选择测试作物,提出了空间任意点光合有效辐射的截获率网格计算方法,建立了批处理自动化空间光能空间分布量化流程,能够准确描述作物群体空间的分布特点,进一步准确量化作物群体空间任意点的光能分布,可以有效对作物品种特性、植被生态环境和作物长势等准确评价。2、棉花是需要强光作物之一,棉花因为具有无限生长特性,冠层结构受品种、水肥条件和气候等因素影响大,具有很大的可塑性,对棉花群体内光合有效辐射进行精确定位定量分析,对棉花生产管理,品种选育,模式配置,株型鉴定等方面具有十分重要的理论和实践指导意义。3、为了测定位置的准确性和方便性,本发明制作了一测定框架,这样便于测试棒的移动距离,测试方便、快捷。本发明可以广泛应用于作物群体内光合有效辐射进行精确定位定量分析中,科学量化光能在不同空间位置的分布,为合理配置作物种植密度,种植模式,评价不同作物株型提供方法和理论依据。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明的作物群体光能空间分布量化方法以棉花作物作为实施例进行说明,但是不限于此,可以适用于其他各种作物群体光能空间分布量化,包括以下步骤:
1、采用空间网格取样方法选择测试作物群体,并在测试作物群体冠层的纵向测试截面区域内均匀设置测试点,并采用PAR测试仪器对每一测试点的PAR进行测试,每一测试点所测得的PAR包括测试点的透射PAR和反射PAR,其中,空间网格取样方法是现有技术空间统计学的常用统计方法,在此不再赘述。
如图1所示,选择某个晴天,采用空间网格取样方法在每个作物群体中选取代表性的棉花2行,在选取的棉花行间以及从地面到冠层上方进行测试。为了方便测试,本发明在选取的棉花行间纵向设置一测试框架(如图1所示的测试框架代表纵向测试截面区域),测试框架采用一矩形框架1,矩形框架1的宽度为棉花的行间距离,矩形框架1的高度可以根据实际需要进行设定,矩形框架1上横向间隔平行设置有若干横向测试杆2,每相邻两个横向测试杆2之间的间距相同,测试点(图1中的黑色圆点表示)根据测试需要均匀分布在每一横向测试杆2上,本实施例中每个横向测试杆2上的相邻测试点的间距为20cm(可以根据需要设置为15~20cm,但是不限于此,可以根据实际需要进行设定),测试点的总数可以根据实际需要设置为30~100,PAR测量仪器的各测试棒在各测定点采用悬挂点的方法在冠层上方同时测定冠层透射PAR和反射PAR;其中,PAR测量仪器的探头水平向上可以测定透射PAR,探头水平向下可以测定反射PAR。
本发明实施例中的PAR测量仪器可以采用100cm长的线型光量子传感器(LI-191SA,LI-COR,Lincoln,NE,USA),将线型光量子传感器可以与一数据收集器(LI-1400,LI-COR,Lincoln,NE,USA)相连,对每个测试点的PAR值的测量结果自动进行记录。
2、根据步骤1测量得到每个测试点的PAR值,采用克里金插值法求解出纵向测试截面区域内其它待测点的PAR值。
为了能够全面对选取的棉花行间光能分布进行准确描述,还需要对纵向测试截面区域中未设置测试点的其它待测位置处的PAR值通过插值进行估算,本发明综合目前十二种对离散数据格网化内插方法、反距离加权插值法、最小曲率法、改进谢别德法、自然邻点插值法、最近邻点差值法、径向基函数插值法、线性插值三角网法、移动平均插值法、数量度量插值法、局部多项式插值法、多元回归法、克里金插值法方法基础上,利用surfer工具筛选出适合光能分布特点的克里金插值法。
由于空间测试点在30~100之间,克里金插值法是以空间结构分析为基础,在估计值满足无偏性和最小方差的前提下求估计值,区域化变量满足二阶平稳假设,待测点X
0的估计值为
其中
是n个已知测试点的函数值,λ
i是n个已知点的全系数,由无偏的条件可知λ
i的和等于1:
式中,φ为拉格朗日算子,γ(Xi,X0)为已知测试点与待测点间的变异函数值,γ(Xi,Xj)为已知测试点间的变异函数值,求出待测点X0的估计值。
3、在纵向测试截面区域内建立坐标系,对位于纵向测试截面区域内的测试点和待测点的位置进行编号生成Grid文件,Grid文件的每个元素Grid(i,j)记录每个测试点或待测点在纵向测试截面区域中的位置,i表示测试点或待测点在纵向测试截面区域中的横坐标,j表示测试点或待测点在纵向测试截面区域中的纵坐标。
本发明实施例中以测试框架的左下角为坐标原点在纵向测试截面区域内建立坐标系,可以根据实际需要选择测试框架的宽和长作为X轴和Y轴,在此不作限定。
4、获取透射和反射的插值Grid文件后,利用矩阵运算方法求取Grid文件中每一元素所对应的光能截获率。
由于作物群体冠层内任意空间点无法直接读出光能利用PAR截获率,因此采用对Grid文件采用矩阵运算方法计算棉花群体纵向测试截面区域内所测点(测试点和待测点)的光合有效辐射截获率,得出纵向测试截面区域内所测点的光能截获率空间分布数值,即:Grid光能截获率(i,j)=1-Grid反射率(i,j)-Grid透射率(i,j),Grid文件每一元素的光能截获率等于1减去该元素所对应的反射PAR再减去该元素所对应的透射PAR。
5、基于求取的纵向测试截面区域中测试点和待测点的光能截获率,对作物群体在整个测试空间区域内的任意点PAR数据进行空间插值生成空间格网文件,并采用多重积分法计算整个测试空间区域的体积,然后除以作物冠层的纵向测试截面区域的面积,得到群体整个光能利用效率。
作物群体生物学产量和整个冠层群体的截获率呈正比,由于光能分布具有高度空间相关性,用人为的或随机的一个点或一个截面的值来代替整个群体的光能利用数值都是不科学的或者是错误的。本发明对作物群体任意点PAR数据采用克里金插值法进行空间插值生成空间格网文件,并采用二重积分法计算整个测试空间区域的体积,然后除以纵向测试截面区域的面积,得到群体整个光能利用效率。本发明利用Simpson’s3/8Rule,计算公式如下:
作物冠层整个光能利用效率=Volume/(纵向测试截面区域)
式中,Δx为空间格网数据的列间距,Δy为空间格网数据的行间距,Gi,j为空间格网数据的第i行第j列的网格结点值。
上述实施例中,本发明采用棉花冠层PAR透射率、反射率、截获率来分析冠层的光分布特征,由于光合辐射量与瞬时太阳高度角和云层关系密切,连续观测绝对值变动性大,可比性差,可以采用比值抵消了不同时间点不确定因素带来的误差,则冠层透射PAR、反射PAR和光能截获率可以采用冠层PAR透射率(Tr)、反射率(Rr)和截获率(In)进行替换,其计算公式可以参考朱相成方法为:
Tr=PARt/PARI
Rr=PARr/PARI
In=(PARI-PARt-PARr)/PARI
式中,PARI为冠层顶部透射PAR,单位为:μmol·m-2·s-1;PARt和PARr分别为冠层底部透射PAR和反射PAR,单位均为:μmol·m-2·s-1。因此可以采用冠层PAR透射率(Tr)、反射率(Rr)和截获率(In)代替步骤1~5中的冠层透射PAR、反射PAR和截获率更为准确低获得整个光能利用效率。
下面通过三个具体实施例进一步对本发明的作物群体光能空间分布量化方法进行详细说明。
实施例1:
试验于2011、2012年在河南省安阳市(36°06′N,114°21′E)中国农业科学院棉花研究所试验地进行,试验棉花品种为鲁棉研28号,构建了6个不同种植密度的棉花群体,分别为15000、33000、51000、69000、87000、105000株/hm2;试验采用随机区组设计,4次重复,每个小区面积为66.6m2,行距为0.80m,每个小区20行。2011年播种时间为4月20日,2012年播种时间为4月22日。
光合有效辐射PAR数据采集利用空间统计学原理,首先采用空间网格方法取样方法在每个密度群体中选取具有代表性的棉花2行,在棉花行间的纵向截面区域放置测试框架,测定每一个测试点的PAR值:探头水平向上测定透射PAR,探头水平向下测定反射PAR。PAR测量使用仪器采用1m长的线型光量子传感器(LI-191SA,LI-COR,Lincoln,NE,USA)和数据采集器(LI-1400,LI-COR,Lincoln,NE,USA),冠层顶部PAR为每隔5秒自动测量记录。本实施例中计算棉花冠层PAR透射率(Tr)、反射率(Rr)和截获率(In)为:
Tr=PARt/PARI
Rr=PARr/PARI
In=(PARI-PARt-PARr)/PARI
整个冠层PAR截获率=1-整个冠层透射率-整个冠层反射率。
对棉花群体任意点PAR数据进行空间插值生成空间格网文件,用二重积分法计算测试空间区域体积,然后除以纵向区域截面面积,即可得到棉花群体整个光能利用效率,本实施例采用3/8扩展Simpson's规则进行计算:
棉花冠层整个光能利用效率=Volume/(纵向测试截面区域)
不同大小的棉花群体冠层株型结构有着各自特点,不同群体棉花茎叶的空间分布决定PAR的分布。棉花枝叶空间的分布及其变化特点直接决定了PAR在棉花群体内的分布特征。如图2(a)所示的群体冠层内PAR透射率空间分布特征表现为,在水平方向上,行中间的透光率要高于棉花行上,纵向上,冠层中上部要高于中下部。群体冠层内PAR反射率空间分布,在水平方向上,行中间的透光率要高于棉花行上,纵向上,冠层中上部要高于中下部,如图2(b)所示。群体内PAR截获率由于是1减去反射率和透射率之和,反射率相比较占的比重很小,所以群体内PAR截获率的空间分布特点正好与透射率相反,如图2(c)所示。
棉花封行前,透光率冠层内分布等值线呈深“V”字形,随着生育进程,棉花植株茎、叶等器官的快速生长,中下部PAR透射率逐渐减少,以2012年播后81天(33000株)群体为例,到了棉花成熟衰老期后,棉花叶片逐渐脱落,中下部透过光率逐渐增加,PAR冠层内分布开始呈现“U”型,其纵向变化比棉花生长初期平缓。从群体PAR反射率来看,整个生育期群体PAR反射率,前期先快速下降,然后缓慢下降后缓慢上升。棉花生育前期,随着群体增大PAR反射率减小,棉花封行以后,不同群体差异逐渐减小。整个生育期群体PAR反射率变化特征呈开口向上的二次曲线,6个群体模拟模型决定系数都在0.9以上。群体PAR最小反射率和群体大小呈负相关,150000株/hm2,310000株/hm2两个群体PAR反射率最小值为0.01、0.008,分别出现在播后124天和121天;而其它4个密度稍大的群体,群体PAR反射率最小值分别依次为0.009、0.007、0.006、0.006,出现时间分别在117、116、116、117天。群体PAR透射率变化曲线也呈开口向上的二次曲线,2011年模拟模型的决定系数都在0.93-0.98之间,2012年在0.89-0.98之间,但其变化幅度远大于群体PAR反射率。2011年全年群体PAR透过率在0.3-1之间,2012年在0.31-1之间。对不同群体而言,通过对其变化模拟方程求导,6个群体的PAR透过率最小值2011年依次为0.43、0.41、0.39、0.37、0.35、0.33,出现时间分别为播后122、123、121、119、116、118天;2012年依次为0.47、0.42、0.42、0.39、0.38、0.34,出现时间分别为105、106、106、104、104、106。
实施例2:
棉花群体光合有效辐射PAR数据采集方法和群体PAR透射率、PAR反射率计算方法与实施例1相同;不同的是本实施例2选用株型不同的棉花品种和新品系9个,分别是中棉所79、中棉所60、鲁棉28、冀228、冀958、ZM3799、ZM120、ZM113、DG002,通过采用9个不同棉花品种构建9个不同类型的棉花群体。试验于2012年在河南省安阳市(36°06′N,114°21′E)中国农业科学院棉花研究所试验地进行,设计试验种植密度为90000株/hm2,行距为0.8m等行设置,播种期为4月22日,小区为10行区,行长8m,宽0.8m,小区面积64.0㎡。针对不同株型棉花品种(系)构建的不同棉花群体,量化其群体光能分布,为不同株型棉花群体提高光能利用率,增加产量提供理论依据。
整个群体的透光率的变化能够准确反应品种的株型特征,不同株型棉花群体整个生育期透光率变化曲线也呈开口向上的二次曲线。短果枝品系DG002和较早熟的ZM120、中棉所60等属于株型比较紧凑的品种(系),这种株型的品种(系)构成的棉花群体透光率高,截获光能少,在播种后的46天、59天、75天、97天、116天和139天,DG002品系的群体透光率分别为0.73、0.53、0.53、0.36、0.47和0.52。冀958,鲁棉研28,冀228等属于株型松散型品种,这种株型的棉花品种构成的棉花群体透光率低,截获光能多,播种后的46天、59天、75天、97天、116天和139天,冀958群体透光率分别为0.58、0.42、0.32、0.29、0.33和0.45。
实施例3:
棉花群体光合有效辐射PAR数据采集方法和群体PAR透射率、PAR反射率计算方法与实施例1相同,不同的是本实施例3选用一个棉花品种中棉所79,采用不同的株行距配置方式,构建了6个不同的棉花群体,6个不同株行距配置方式分别是:80cm等行距种植和90cm+20cm、90cm+10cm、80cm+20cm、80cm+10cm、70cm+20cm等5个宽窄行双行配置模式。试验于2012年在河南省安阳市(36°06′N,114°21′E)中国农业科学院棉花研究所试验地进行。大田直播播种期在4月22日,小区长8m,宽8m,种植密度为90000株/hm2。针对不同株行距配置方式构建的6个棉花群体,量化其群体光能分布,为不同株行距配置方式棉花群体提高光能利用率,增加产量提供理论依据。
作物种植行距的变化对作物群体构建及其光能分布和利用影响较大,不同株行距配置的棉花群体整个生育期透光率变化曲线也呈开口向上的二次曲线,6个不同株行距配置的棉花群体中,以70cm+20cm宽窄行种植模式群体透光率最低,截获光能较多,在播种后的38天、46天、59天、66天、75天、86天、97天、109天、116天、128天和138天,其群体透光率分别为0.89、0.67、0.47、0.42、0.36、0.26、0.29、0.28、0.34、0.31、0.42;而90cm+20cm宽窄行种植模式群体最高,截获光能较少,在播种后的38天、46天、59天、66天、75天、86天、97天、109天、116天、128天和138天,其群体透光率分别为0.91、0.67、0.55、0.56、0.45、0.33、0.31、0.47、0.49、0.50、0.56。不同株行距配置的棉花群体光分布及透光率和截获光能的量化,为构建合理的高光效群体提供了理论依据。
上述各实施例仅用于说明本发明,其中方法的各实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。