免费文献传递   相关文献

Diagnosis of N nutrition of rice using digital image processing technique

应用数字图像技术进行水稻氮素营养诊断



全 文 :植物营养与肥料学报 2015,21(1):259-268
JournalofPlantNutritionandFertilizer doi牶1011674/zwyf.20150129
收稿日期:2013-12-30   接受日期:2014-08-25
基金项目:公益性行业(农业)科研专项(201103039);中央高校基本科研业务费专项(2013PY113);长江学者和创新团队发展计划(IRT1247)
资助。
作者简介:李岚涛(1987—),男,河南周口人,硕士,主要从事作物营养与现代施肥技术。Email:lijianwei@webmail.hzau.edu.cn
 通信作者 Email:lujianwei@mail.hzau.edu.cn
应用数字图像技术进行水稻氮素营养诊断
李岚涛,张 萌,任 涛,李小坤,丛日环,吴礼树,鲁剑巍
(华中农业大学资源与环境学院,农业部长江中下游耕地保育重点实验室,湖北武汉 430070)
摘要:【目的】研究田间试验条件下水稻不同生育期冠层图像色彩参数(G、NRI、NGI、NBI、G/R和 G/B)及植株
氮素营养指标(叶片含氮量、植株全氮含量、生物量、氮素累积量和冠层NDVI值)的时空变化特征,并分析两者间
的相关性,确立水稻氮素营养诊断的最佳色彩参数和方程模型,为探明数码相机在水稻上的适宜性及精确诊断水
稻氮素营养状况提供理论基础。【方法】于2013年5月 9月在湖北省武汉市华中农业大学试验基地(30°28′08′′
N,114°21′36′′E)采用不同施氮处理的田间试验,以籼型两系杂交稻“两优6326”为供试作物,设置4个施氮水平:
0、75、150和225kg/hm2(分别以N0、N75、150和N225表示),3次重复,随机区组排列。分别在水稻分蘖期、拔
节期、孕穗期和灌浆期采用数码相机(NikonD700,1200万像素)获取水稻冠层图像,应用Adobephotoshop70软件
直方图程序提取图像的红光值R、绿光值G和蓝光值B,研究数码相机进行水稻氮素营养诊断色彩参数,确定植株
氮素营养指标诊断模型。【结果】较对照(N0)相比,分蘖期、拔节期、孕穗期和灌浆期3个施氮处理水稻地上部生
物量、叶片含氮量、植株全氮含量、氮素累积量、冠层 NDVI值和成熟期产量增幅分别平均为407% 980%、
424% 724%、362% 853%、1255% 2091%、513% 606%和601% 1170%,差异显著。水稻不
同生育期各冠层数字化指标G、NRI、NGI、NBI、G/R和G/B与上述氮素营养参数相关性差异较大,且以数字图
像红光标准化值NRI表现最佳,建议作为应用数码相机进行水稻氮素营养诊断的最佳冠层图像色彩参数指标。进
一步分析表明,可以用统一的线性回归方程来描述不同生育期、不同氮素水平下水稻植株氮素营养指标随冠层色
彩参数NRI的变化模式。【结论】数码相机进行水稻氮素营养诊断测试结果稳定,具有快速、便捷、非破坏性等优
点,冠层色彩参数NRI与水稻氮素营养指标和产量之间均表现出较好的相关性,满足氮素营养无损诊断的需求,对
实时、快速监测水稻长势状况及氮素营养丰缺水平具有较高的可行性,有望发展成为新时期作物氮素营养无损诊
断技术的潜力。
关键词:数字图像;水稻;氮素营养诊断
中图分类号:S126;S511   文献标识码:A   文章编号:1008-505X(2015)01-0259-10
DiagnosisofNnutritionofriceusingdigitalimageprocessingtechnique
LILantao,ZHANGMeng,RENTao,LIXiaokun,CONGRihuan,WULishu,LUJianwei
[ColegeofResourcesandEnvironment,HuazhongAgriculturalUniversity/KeyLaboratoryofArableLand
Conservation(MiddleandLowerReachesofYangtseRiver),MinistryofAgriculture,Wuhan430070,China]
Abstract:【Objectives】Thespatialandtemporaldistributionofcolorindexesofcanopy(G,NRI,NGI,NBI,
G/RandG/B)andtheindexesofNnutritioninriceplantswerestudiedtodeterminethebestcolorparametersand
regressionequationsfornitrogenwithadigitalcamera,andprovideatheoreticalbasisandtechnicalapproachfor
monitoringplantnitrogenstatusofriceandprecisionmanagementofnitrogenfertilization.【Methods】Field
experimentswerecariedoutfromMaytoSeptember2013inricegrowingseasonattheExperimentalFarmof
HuazhongAgriculturalUniversity(30°28′08′′N,114°21′36′′E).Atwolineindicahybrid(Liangyou6326)was
chosenastestcultivar.FourNtreatments:N0,75,150and225kg/haweredesignedandrecordedasN0,N75,
N150,N225,respectively.Acolordigitalcamera(D700,Nikon,Japan)witharesolutionof120megapixels
植 物 营 养 与 肥 料 学 报 21卷
wasemployedtocapturecolorimagesofricecanopyatthetilering,jointing,bootingandfilingstage.Theimage
resolutionwas1936×1296pixelsof14bitforred,greenandblue.Theimagesofsize272MBweretransferedin
jointphotographicexpertsgroup(JPEG)formattoacomputerandprocessedwithAdobePhotoshop70toextract
colorinformationforstudyingpossibilityofusingdigitalimageanalysismethodforplantNstatusdiagnosisandfor
determiningthebestcolorparametersandregressionequationsinrice.【Results】ComparedwithN0,thebiomass,
leafNcontent,planttotalNcontent,NaccumulationandcanopyNDVIvaluesofricewithNtreatmentsatthe
tileringstage,jointingstage,bootingstageandfilingstageareaveragelyincreasedby407%-980%,424%-
724%,362%-853%,1255%-2091% and513%-606%,respectively,atmaturestage,theaverage
increaseofyieldarefrom601% to1170%.Comparedwithotherplantcanopycolorparameters,thenormalized
rednessintensity(NRI),calculatedasR/(R+G+B),ismuchbeterasithasbetercorelationswiththeindex
ofNnutrientinriceplants,riceyieldandcanopyNDVIvaluesatthediferentgrowthstages.Anintegratedlinear
regressionequationcouldbeusedfordescribingtherelationshipbetweenNRIandleafNcontent,planttotalN
content,Naccumulation,canopyNDVIvalues,yieldatdiferentgrowingstagesandnitrogenlevelsofrice.
【Conclusions】ThedigitalimageprocessingtechniqueisusableindetectingtheNnutritioninriceandhas
advantagesoffast,stableresults,easytoapply,andnondestructive.Simultaneously,thecolorparameterNRIhas
apreferableinterelationwiththeindexofNnutritionandyieldthanothers,thusNRIissuitableforrapiddiagnosis
ofnitrogennutrition,andthedigitalimageprocessingtechniquemethodshowsthepotentialofbeingusedinfast
ricenitrogendiagnosiswithoutdamage.
Keywords牶digitalimage牷rice牷Nstatusdiagnosis
  传统的氮素营养诊断和氮肥推荐多以田间采集
植株样品、实验室化学分析为主,该方法在样品的
采集、测试以及数据处理等方面需耗费大量的人
力、物力和财力,时效性差,不利于推广应用[1-2]。
数字图像技术通过判断作物冠层颜色深浅了解作物
氮素营养状况继而进行施肥推荐或产量预测,具有
操作简便、快捷、非破坏性等优点[3-4]。此外,由于
数码相机价格的不断下降以及所采用图像感应器件
的快速发展,使通过数码相机测量光反射强度来评
价作物氮素营养状况成为可能[5],近年来在冬小
麦[6]、蔬菜[7]、棉花[5]和玉米[8]等多种作物上得到
广泛应用。李红军等[6]利用数码相机对小麦冠层
进行拍照,通过分析色彩参数与作物氮素营养状况
之间的关系,认为返青期反光叶面的 G/R与 R/(R
+G+B)值能较好的反映小麦氮素营养状况。宋述
尧等[7]通过分析黄瓜冠层数字化诊断指标与施氮
量、土壤无机氮(Nmin)及植株氮素营养之间的关
系,得出了以绿光标准化值 G/(R+G+B)作为诊
断黄瓜结果期氮素营养状况的适宜数字化参数。
Blackmer等[9]通过分析玉米冠层航空图片,指出航
片红光深度绝对值(R)、绿光深度绝对值(G)和蓝
光深度绝对值(B)均与玉米产量呈显著线性相关
关系。
水稻作为一种重要的粮食作物,其生长发育和
产量形成受氮素营养的影响极大。前人利用数码相
机进行水稻氮素诊断的研究已有一些报道。王远
等[10]研究表明拔节期和孕穗期水稻红光标准化值
NRI与叶片 SPAD值和含氮量具有良好的相关关
系。KyuJong等[11]利用数码相机获取水稻冠层图
像,指出绿光深度绝对值(G)与地上部生物量、氮
素累积量和叶面积指数均可用指数方程进行拟合,
R2分别为083、083和086。陆秀明等[12]通过图
像处理获取水稻冠层图像的覆盖度,利用覆盖度估
测的叶面积指数与实测叶面积指数之间的相关系数
在077 093,达到显著水平。贾良良等[1]对应用
数码相机进行水稻氮素营养诊断的可行性进行研
究,认为数字图像色彩分析技术具有发展成为新一
代作物氮营养诊断的潜力。上述研究大都采用不同
手段将作物冠层色彩参数与地面覆盖度或作物生长
状况进行分析,而对于整个作物冠层全生育期的色
彩参数变化特征及模型构建与筛选的研究相对较
少。在此基础上,本文通过研究水稻生育期地上部
生物量、成熟期产量及植株氮素营养指标的变化规
律,并分析其与冠层数码指标之间的相关性,从而确
定水稻氮素营养诊断的最佳色彩参数及方程模型,
为探明数字图像技术在水稻上的适应性及精确诊断
水稻氮素状况和推荐施肥提供理论基础。
062
1期    李岚涛,等:应用数字图像技术进行水稻氮素营养诊断
1 材料与方法
11 试验材料
试验于2013年5月至9月在湖北省武汉市华
中农业大学试验基地内进行。供试土壤为黄棕壤,
前茬作物为油菜,0—20cm耕层土壤 pH741、有
机质 672g/kg、全氮 032g/kg、碱解氮 2396
mg/kg、速效磷716mg/kg、速效钾1805mg/kg。
供试品种为杂交水稻“两优6326”。
12 试验设计与方法
氮肥设置0、75、150和225kg/hm2四个水平,
分别用N0、N75、N150和 N225表示,分3次施用,
即50%基肥、35%分蘖肥和15%穗肥。磷、钾肥按
P2O575kg/hm
2、K2O120kg/hm
2施入,磷肥全部基
施,钾肥按70%基肥和30%分蘖肥施用。氮、磷、
钾肥品种分别为尿素(含 N46%)、过磷酸钙(含
P2O512%)和氯化钾(含 K2O60%)。小区面积20
m2(5m×4m),3次重复,随机区组排列。2013年5
月6日育苗,6月18日移栽。定植密度为20×105
蔸/hm2,2013年9月21日收获。
13 数字图像获取与处理
选取晴天分别于水稻分蘖期(2013年 7月 11
日)、拔节期(2013年7月22日)、孕穗期(2013年
7月31日)和灌浆期(2013年8月11日),在太阳
高度角相对稳定的12:00 13:00之间,同时以标
准白色比色板进行校正,利用数码相机(Nikon
D700,1200万像素)获取水稻冠层图像。拍摄时,将
相机镜头距离地面15m垂直高度(距水稻冠层约
10m),与地面呈 60°夹角进行,同时将相机调至
Auto模式下,以自动曝光控制色彩平衡(图1)。图
片以JPEG格式存储,分辨率为1936×1296像元,应
用Adobephotoshop70软件直方图程序获取图像的
红光值 R(rednessintensity)、绿光值 G(greenness
intensity)和蓝光值 B(bluenessintensity)[4,6]。备选
的数码图像特征值包括:绿光深度绝对值(G)、绿
光深度绝对值与红光深度绝对值比值(G/R)、绿光
深度绝对值与蓝光深度绝对值比值(G/B)、红光标
准化值(NRI)、绿光标准化值(NGI)和蓝光标准化
值(NBI)。各标准化值的计算如下:
红 光 标 准 化 值 (NRI,normalized redness
intensity)=R/(R+G+B)
绿光 标 准 化 值 (NGI,normalizedgreenness
intensity)=G/(R+G+B)
蓝 光 标 准 化 值 (NBI,normalized blueness
图1 数码相机进行水稻氮素营养诊断示意图
Fig.1 SketchmapofdiagnosisofNstatusof
riceusingdigitalcamera
[注(Note):A—数码相机Digitalcamera;B—水稻冠层Ricecanopy;
C—相机距地面垂直距离 15mThecamerawaspositioned15m
abovetheground;D—相机距水稻水平距离 Horizontaldistance
betweenAandB.]
intensity)=B/(R+G+B)
14 测定项目与方法
冠层NDVI值的测定:用美国 Ntech公司生产
的Greenseeker手持光谱仪采集水稻冠层 NDVI数
据。测量时,将光照探头平行于水稻冠层,垂直高度
约60cm。各处理小区随机选取4个1m2左右区域
进行测定,结果直接通过掌上电脑进行采集。
植株与叶片全氮含量的测定:分别于上述各时
期内每小区选取代表性水稻5株,分叶片和茎秆两
部分,105℃下杀青30min,60℃烘箱中烘至恒重,采
用H2SO4-H2O2法消化,AA3流动注射分析仪测定
含氮量[13]。
15 数据处理与分析
数据处理与分析采用 MicrosoftExcel2003及
SPSS170进行,LSD法检验 P<005水平上的差
异显著性;制图采用Origin85软件。
2 结果与分析
21 氮肥施用对水稻生物量和产量的影响
各生育时期水稻地上部生物量均以 N225处理
最大(表1)。施氮有利于水稻干物质量积累,与对
照(N0)相比,分蘖期、拔节期、孕穗期和灌浆期施
氮处理水稻地上部生物量分别平均增加 407%、
980%、622%和625%。成熟期产量结果表明,
随氮肥用量增加,水稻产量呈显著升高趋势,与不施
162
植 物 营 养 与 肥 料 学 报 21卷
表1 不同氮肥用量水稻生物量和产量 (t/hm2)
Table1 Efectofnitrogen(N)applicationratesonbiomassandyieldofrice
处理
Treatment
生物量 Biomass
分蘖期
Tileringstage
拔节期
Jointingstage
孕穗期
Bootingstage
灌浆期
Filingstage
产量
Yield
N0 123Bc 230Bc 409Bb 416Bb 352Cd
N75 175ABb 342Bb 589ABb 618ABa 564Bc
N150 222Aa 492Aa 640ABab 672ABa 678Ab
N225 224Aa 531Aa 762Aa 740Aa 764Aa
  注(Note):表中不同的小写字母和大写字母分别表示达到5%和1%上的差异显著水平 Diferentlowercaseletersorcapitalleterswithinthe
samecolumnindicatesignificantdiferencesatthe5% and1% probability
氮相比,施氮处理产量增幅高达601% 117%。
22 不同施氮水平对水稻氮素营养参数及冠层
NDVI值的影响
施氮可显著提高水稻植株氮素营养指标及冠层
NDVI值(表2)。叶片含氮量、植株全氮含量及氮
素累积量方面,与对照(N0)相比,分蘖期、拔节期、
孕穗期和灌浆期三者分别平均增加 724%、
485%、424%、610%,853%、418%、362%、
516%和2091%、1918%、1255%、1499%,差
异极显著。对于冠层NDVI值,氮肥施用后,上述各
时期内分别平均提高020、023、020、020个单
位,差异显著。
表2 不同氮素水平对水稻植株氮素营养指标及冠层NDVI值的影响
Table2 EfectsofdiferentnitrogentreatmentsontheindexofNnutrientinriceplantsandcanopyNDVI
valuesrespectivelyatdiferentgrowthstages
时期
Growthstage
处理
Treatment
叶片含氮量(%)
LeafNcontent
植株全氮含量(%)
PlanttotalNcontent
氮素累积量(t/hm2)
Naccumulation
NDVI值
NDVIvalue
分蘖期
Tileringstage
N0 1.46Cc 1.05Cc 1.31Cc 0.39Bc
N75 2.07Bb 1.57Bb 2.70Bb 0.55Ab
N150 2.72Aa 2.07Aa 4.58Aa 0.63Aa
N225 2.76Aa 2.17Aa 4.86Aa 0.60Aab
拔节期
Jointingstage
N0 1.21Cc 0.90Cc 2.06Cc 0.38Cc
N75 1.37Cc 0.93Cc 3.19Cc 0.55Bb
N150 1.79Bb 1.26Bb 6.20Bb 0.62ABa
N225 2.23Aa 1.64Aa 8.69Aa 0.65Aa
孕穗期
Bootingstage
N0 1.06Cc 0.70Ab 2.86Bb 0.33Ad
N75 1.27BCb 0.90Aab 5.40ABab 0.45Bc
N150 1.44Bb 0.84Aab 5.36ABab 0.54Ab
N225 1.82Aa 1.12Aa 8.61Aa 0.60Aa
灌浆期
Filingstage
N0 1.00Cc 0.62Cc 2.59Bc 0.35Cd
N75 1.35BCb 0.74BCbc 4.58ABbc 0.47Bc
N150 1.58ABb 0.93ABb 6.22ABab 0.56Ab
N225 1.90Aa 1.15Aa 8.60Aa 0.60Aa
  注(Note):表中不同的小写字母和大写字母分别表示达到5%和1%上的差异显著水平 Diferentlowercaseletersorcapitalleterswithinthe
samecolumnindicatesignificantdiferencesatthe5% and1% probability
262
1期    李岚涛,等:应用数字图像技术进行水稻氮素营养诊断
23 水稻冠层 NDVI值与植株氮素营养参数的相
关性分析
水稻整个生育期不同施氮量的冠层 NDVI值与
叶片含氮量、植株全氮含量、生物量和氮素累积量
相关性均达到极显著水平(表 3)。该结果表明,
NDVI值能够成熟表征水稻植株氮素营养指标(叶
片含氮量、植株全氮含量和氮素累积量)及群体生
长发育状况(生物量)。同时,以 NDVI值作为一种
标准变量来反演水稻冠层数字化指标具有较高的可
行性。
表3 不同生育时期水稻冠层NDVI值与植株氮素营养指标的相关系数(r)
Table3 CorrelationbetweencanopyNDVIvaluesandriceindexofNnutrient
时期
Growthstage
叶片含氮量(%)
LeafNcontent
植株全氮含(%)
PlanttotalNcontent
生物量(t/hm2)
Biomass
氮素累积量(t/hm2)
Naccumulation
分蘖期 Tileringstage 0897 0876 0887 0924
拔节期 Jointingstage 0649 0742 0865 0951
孕穗期 Bootingstage 0935 0730 0822 0883
灌浆期 Filingstage 0935 0869 0863 0821
  注(Note):n=12;,分别表示在5%和1%水平下的显著性 Meansignificantatthe5% and1% levels,respectively.
24 水稻冠层数字化指标与植株氮素营养参数及
冠层NDVI值的相关性分析
对水稻不同生育期各冠层数字化指标G、G/R、
G/B、NRI、NGI和 NBI与氮素营养参数及冠层
NDVI值进行回归分析,结果表明,G、G/R、G/B、
NGI和NBI与植株及叶片含氮量、生物量、氮素累
积量和NDVI值均呈线性正相关关系,而 NRI呈极
显著线性负相关关系。分蘖期和孕穗期,以 NRI和
G/R与上述各回归参数相关性较高,且 rNRI显著高
于rG/R。拔节期和灌浆期,则以 NRI和 NBI相关性
较好,且rNRI明显大于 rNBI。因此,以 NRI作为应用
数码相机进行水稻氮素营养诊断的指标选择,具有
较高的代表性和稳定性。各生育期间,图像特征参
数NRI与水稻氮素营养指标及NDVI值相关性并无
明显差异,而分蘖期作为水稻氮素吸收利用旺盛时
期,是确定水稻有效分蘖数、确保高产稳产的重要
基础,较拔节期、孕穗期和灌浆期而言,其相关性并
不是最优的,主要原因在于分蘖期(2014年7月11
日)时水稻植株较小,而土壤和田面水背景相对凸
出,数码相机在提取水稻冠层色彩信息(作物在不
同氮营养水平下所表现出的色彩差异,在 400
700nm可见光波段反射率较小,720 1300nm近红
外波段反射率较大)时受到一定影响,这也是作物
无损诊断精确性相对较低的原因之一,受其他因素
影响较多,如环境条件、病虫害、品种等。
25 水稻冠层数字化诊断指标NRI的可行性分析
及诊断模型
为更精确诊断水稻冠层数字化图像参数 NRI,
本研究在分析其与植株氮素营养指标相关性时,进
一步利用线性和非线性回归分析,以氮素营养指标
为自变量(x),图像指标 NRI(y)为因变量,采用线
性函数、抛物线函数、幂函数、指数函数和对数函
数构建不同诊断模型(表5)。结果表明,不同回归
分析方法间水稻冠层图像参数 NRI与其氮素营养
指标均表现出较好相关性,达到极显著水平,且线性
和非线性模型拟合效果非常接近。考虑到线性模型
的数学表达式较为简单,因此选取以氮素营养指标
为自变量(x)所构建的图像诊断指标 NRI(y)直线
方程诊断模型(图2)。同时建立基于线性回归方程
对水稻冠层图像参数 NRI与氮肥用量(R2 =
0747)、产量(R2 =0778)、NDVI(R2 =
0798)、叶片氮含量(R2=0445)、植株全氮
含量(R2=0326)和氮素累积量(R2=0460)
的定量诊断方程(图2),为便捷、无损监测和反演
水稻氮素营养状况提供理论基础,同时也为数码相
机进行水稻氮素营养诊断提供一个快捷数字化指标
(NRI)和方法(线性方程)。
3 讨论与结论
31 应用数码相机进行水稻氮素营养诊断指标的
选择
作物冠层颜色深浅与氮素营养状况密切相关,
不同的营养水平、光温条件及群体结构均显著影响
作物的生长发育,从而使植株冠层呈现出不同的颜
色[14]。因此,选取合适的冠层数码诊断指标,对准
362
植 物 营 养 与 肥 料 学 报 21卷
表4 不同生育时期水稻特征光谱参数与植株氮素营养指标及冠层NDVI值的相关系数(r)
Table4 CorrelationbetweenthefeaturespectralparametersandNnutritionindexesandcanopyNDVIvalues
atdiferentgrowthstages
时期
Growthstage
指标
Index
叶片含氮量(%)
LeafNcontent
植株全氮含量(%)
PlanttotalN
生物量(t/hm2)
Biomass
氮素累积量(t/hm2)
Naccumulation
NDVI值
NDVIvalue
分蘖期
Tileringstage
G 0552 0548 0446 0525 0467
NRI -0842 -0828 -0891 -0846 -0928
NGI 0561 0548 0691 0575 0759
NBI 0675 0660 0643 0671 0615
G/R 0794 0782 0856 0800 0896
G/B 0414 0416 0318 0413 0243
拔节期
Jointingstage
G 0809 0747 0881 0839 0894
NRI -0892 -0847 -0911 -0899 -0939
NGI 0328 0293 0413 0336 0530
NBI 0898 0852 0825 0877 0785
G/R 0808 0767 0857 0824 0904
G/B 0802 0758 0753 0786 0713
孕穗期
Bootingstage
G 0868 0745 0729 0776 0810
NRI -0903 -0756 -0803 -0798 -0947
NGI 0563 0530 0580 0532 0727
NBI 0551 0462 0436 0481 0497
G/R 0717 0650 0667 0645 0854
G/B 0137 0021 0015 0032 0038
灌浆期
Filingstage
G 0151 0115 0146 0023 0031
NRI -0823 -0771 -0686 -0745 -0883
NGI 0021 0025 0097 0016 0173
NBI 0709 0672 0594 0634 0758
G/R 0505 0433 0500 0435 0677
G/B 0551 0528 0477 0511 0589
  注(Note):n=12;,分别表示在5%和1%水平下的显著性 Meansignificantatthe5% and1% levels,respectively.
确评估作物氮素营养状况及氮肥推荐具有重要意
义[15]。李红军等[6]研究表明,拔节期冬小麦不反光
叶面和反光叶面红光标准化值 NRI(R/R+G+B)
与植株全氮含量存在较好相关性。白金顺等[16]研
究发现,灌浆期红光标准化值(R/(R+G+B))与
春玉米氮素营养指标相关性较好,且与玉米产量间
呈显著直线回归关系,回归系数为 46%。肖焱波
等[17]应用数码相机在冬小麦上研究指出,数字图像
红光标准化值 R/(R+G+B)与冬小麦常规营养诊
断指标(SPAD值、茎基部硝酸盐含量、植株全氮含
量、生物量)的相关系数达到了0809 0946,是
较好表征作物氮营养状况诊断指标。综合前人研究
结果,利用红光标准化值进行作物氮营养诊断较其
他光谱指数相比虽可行性较高,但主要以单一生育
期为主,而系统研究全生育期色彩参数变化及方程
模型构建与筛选并得到结果相对稳定和一致的则相
对较少,这也为无损、便捷和快速诊断作物氮素营
养与生长发育状况带来一定困难。本试验通过对比
分析水稻分蘖期、拔节期、孕穗期和灌浆期6种数
码图像特征值与植株氮素营养指标间的相关性,得
462
1期    李岚涛,等:应用数字图像技术进行水稻氮素营养诊断
出冠层色彩参数 NRI为应用数码相机进行水稻氮
素营养诊断的最佳指标,同时建立了基于 NRI的水
稻叶片氮含量、植株全氮含量、氮肥用量、氮素累
积量、NDVI值和产量的线性回归方程模型,为应用
数码相机无损、快速研究水稻氮素营养状况提供数
码指标参考和理论依据。
表5 水稻冠层NRI与其氮素营养指标的线性和非线性回归分析(n=48)
Table5 LinearandnonlinearregerssionanalysisbetweencanopyNRIandriceindexofNnutrient
模型类型 Modeltype
直线方程
Linear
二次方程
Parabola
幂方程
Power
指数方程
Exponential
对数方程
Logarithm
R2 F R2 F R2 F R2 F R2 F
氮肥用量
Napplicationrate
0747 1356 0770 755 0755 1415
产量 Yield 0778 1611 0786 827 0778 1617 0747 1358 0743 1331
NDVI 0798 1813 0798 887 0785 1684 0789 1723 0772 1557
叶片含氮量
LeafNcontent
0445 3685 0560 286 0453 380 0504 467 0510 479
植株全氮含量
PlanttotalN
0326 223 0409 156 0334 231 0375 276 0383 286
氮素累积量
Naccumulation
0460 392 0551 276 0460 393 0547 556 0538 536
  注(Note):,分别表示在5%和1%水平下的显著性 Meansignificantatthe5% and1% levels,respectively.
32 应用数码相机进行水稻氮素营养诊断的可
行性
作物冠层色彩信息是表征其营养状况的重要组
成部分,使其在不同的营养状况下表现出差异明显
的茎叶颜色,尤其是受氮肥供给量的影响[18]。在
400 700nm可见光波段范围内,由于叶绿素和其
他色素吸收,植物冠层光反射率较低,且随氮肥用量
增加而降低。700 1300nm近红外波段,由于细胞
排列方式和冠层结构影响,植物对光吸收很弱,表现
出较强的反射率,并随氮肥用量增加而提升[19]。作
物对光谱的吸收和反射差异最终表现在群体(冠
层)颜色上,由于作物冠层绿色变化通常与叶绿素
含量有关,而植株叶绿素含量与氮素营养相关密切。
因此,通过分析作物冠层颜色深浅了解作物氮素营
养状况,继而进行推荐施肥或产量预测的研究日益
增多[3]。数码相机通过图像传感器获取物体图像,
将图像转换成数字图像,并利用计算机模拟人的功
能区分识别图像色彩信息,达到分析作物冠层色彩
差异和进行诊断的目的[18]。此外,应用数码相机和
数字图像处理技术进行氮素营养诊断较其他方法
(如高光谱和SPAD仪法)相比具有成本低、覆盖面
大和分辨率高等优点,易于推广应用[20]。Thomas
等[21]通过试验证明了光谱技术诊断甜椒氮素营养
的可行性。Adamsenden等[22]利用数码相机成功估
算了油菜花的数量,得出了花的数目与像素之间的
决定 系 数 (R2)为 083,达 到 极 显 著 水 平。
Richardson等[23]应用数字图像处理技术对草地地
面覆盖度成功进行了预测,其预测值与实测值之间
R2高达099。上述研究将利用数码相机所获取作
物冠层参数与作物营养及生长状况进行分析,均表
明应用数字图像技术可用来评估田间作物的营养状
况[5]。本试验条件下,在所分析的6种数码图像特
征值中,水稻冠层数字指标红光标准化值 NRI与植
株氮素营养参数相关性较好且表现稳定,均达到显
著或极显著水平,且采用线性方程模型拟合 NRI与
氮肥用量、产量、NDVI、叶片氮含量、植株全氮含
量和氮素累积量相关性较好(图2),表明应用数码
相机定量、快捷研究和反演水稻氮素营养状况具有
较高可行性。
33 有待进一步研究和说明的问题
氮作为影响水稻生长发育和产量形成的必需元
素之一,在生产实际中,氮肥过量施用的现象仍然非
562
植 物 营 养 与 肥 料 学 报 21卷
图2 水稻冠层NRI与氮素营养指标的线性回归分析 (n=48)
Fig.2 LinearregressionanalysisbetweenNRIandriceindexofNnutrient
常普遍。冯洋等[24]在湖北黄冈调查结果表明,目前
农民水稻施氮量普遍在 N240 300kg/hm2之间。
葛梦婕等[25]通过对长江中下游7个粳型超级稻试
验指出,超级粳稻在施氮量300kg/hm2条件下产量
最高,普通粳稻的最高产量施氮量也高达 263
kg/hm2。彭少兵等[26]在湖南长沙曾观察到,超级杂
交稻两优培九即使施氮量达到240kg/hm2时仍没
有出现倒伏现象,而当农民应用这些新品种或组合
时,通常会施用更多的氮肥以获得高产。本试验条
件下,当施氮量分别为0、75、150和225kg/hm2时,
水稻产量分别为3520、5637、6778和7640kg/hm2,
呈线性显著增加趋势。因此,本试验仅对 “缺氮—
氮优化用量”之间应用数字图像技术进行水稻氮素
营养诊断提供了一种方法和思路。在此基础上,探
索“缺氮—氮优化—氮过剩”条件下水稻冠层数码
指标的时空变化特征及可行性,仍需我们开展相关
研究。此外,本研究是在其他条件相同(磷钾肥用
量充足且一致、拍照时间统一且限定在晴天太阳高
度角变化较小的时间段12:00 13:00之间、相机
拍摄角度保证一致等)而氮肥水平差异(N0、N75、
N150和N225)情况下得出的,而水稻不同生育期冠
层颜色变化受很多因素影响,如品种、生长环境、
养分类型及用量等,因此,探索不同条件下水稻冠层
数码指标的差异机理及可行性仍需我们开展相关
研究。
本试验通过分析不同氮素处理下各数字化诊断
662
1期    李岚涛,等:应用数字图像技术进行水稻氮素营养诊断
参数与水稻叶片含氮量、植株全氮含量、氮素累积
量、生物量、冠层NDVI值和产量之间的关系,得出
了以红光标准化值 NRI作为应用数码相机诊断水
稻氮素营养状况的指标选择,继而对比分析不同方
程模型间差异及适宜性,建立了基于 NRI的水稻植
株氮素营养指标线性方程模型。但该结果仅是在单
一试验品种(两优6326)和试验年限(一年)下得出
的,要阐明规律性和实用性更强的结果尚需进一步
试验验证。同时,作物氮素营养无损诊断的核心是
建立推荐追肥体系,并以此进行施肥推荐[15]。因
此,建立基于上述试验结果的不同生育期临界 NRI
值及推荐施氮量,继而以追肥的方式提供氮肥,仍需
要更多、更广泛的试验来验证。
参 考 文 献:
[1] 贾良良,范明生,张福锁,等应用数码相机进行水稻氮素营
养诊断[J].光谱学与光谱分析,2009,29(8):2176-2179
JiaLL,FanMS,ZhangFSetalNitrogenstatusdiagnosisof
ricebyusingadigitalcamera[J].SpectroscopyandSpectral
Analysis,2009,29(8):2176-2179
[2] 屈卫群,王绍华,陈兵林,等棉花主茎叶SPAD值与氮素营
养诊断研究[J].作物学报,2007,33(6):1010-1017
QuWQ,WangSH,ChenBLetalSPADvalueofcotonleaves
onmainstemandnitrogendiagnosisforcotongrowth[J].Acta
AgronomiceSinica,2007,33(6):1010-1017
[3] 李桂娟,朱丽丽,李井会作物氮素营养诊断的无损测试研
究与应用现状[J].黑龙江农业科学,2008(4):127-129
LiGJ,ZhuLL,LiJHPresentstatusofresearchandapplication
ofnondestructivemeasurementofnitrogennutritiondiagnosis[J].
HeilongjiangAgriculturalSciences,2008(4):127-129
[4] 贾良良应用数字图像技术与土壤植株测试进行冬小麦氮营
养诊断[D]北京:中国农业大学博士学位论文,2003
JiaLLTodectectnitrogenstatusofwinterwheatbyusingcolor
digitalcamera[D]Beijing:PhDdissertationofChinaAgricultural
University,2003
[5] 王娟,雷咏雯,张永帅,等应用数字图像分析技术进行棉花
氮素营养诊断的研究[J].中国生态农业学报,2008,16(1):
145-149
WangJ,LeiYW,ZhangYSetalDiagnosisofcotonNstatus
usingdigitalimageanalysistechnique[J].ChineseJournalofEco
Agriculture,2008,16(1):145-149
[6] 李红军,张立周,陈曦鸣,等应用数字图像进行小麦氮素营
养诊断中图像分析方法的研究[J].中国生态农业学报,
2011,19(1):155-159
LiHJ,ZhangLZ,ChenXMetalImageanalysismethodin
applicationofdigitalimageondiagnosingwheatnitrogenstatusJ]
ChineseJournalofEcoAgriculture,2011,19(1):155-159
[7] 宋述尧,王秀峰数字图像技术在黄瓜氮素营养诊断上的应
用研究[J].吉林农业大学学报,2008,30(4):460-465
SongSY,WangXFDiagnosisofNstatusofcucumberusign
digitalimage processing technique[J]. JournalofJilin
AgricurturalUniversity,2008,30(4):460-465
[8] 孙钦平,李吉进,邹国元,等应用数字图像技术对有机肥施
用后玉米氮营养诊断研究[J].光谱学与光谱分析,2010,30
(9):2447-2450
SunQP,LiJJ,ZouGYetalNitrogenstatusdiagnosisbyusing
digitalphotographyanalysisfororganicfertilizedmaize[J].
SpectroscopyandSpectralAnalysis,2010,30(9):2447-2450
[9] BlackmerTM,SchepersJS,VarelG ELightreflectance
comparedwithothernitrogenstressmeasurementsincornleaves
[J].JorunalofProductionAgriculture,1994,6:934-938
[10] 王远,王德建,张刚,王灿基于数码相机的水稻冠层图像
分割及氮素营养诊断[J].农业工程学报,2012,28(17):
131-136
WangY,WangDJ,ZhangG,WangCDigitalcamerabased
imagesegmentationofricecanopyanddiagnosisofnitrogen
nutrition[J].TransactionsoftheChineseSocietyofAgricultural
Engineering,2012,28(17):131-136
[11] KyuJL,ByuuW LEstimationofricegrowthandnitrogen
nutritionstatususingcolordigitalcameraimageanalysis[J].
EuropeanJorunalofAgronomy,2013,(6):57-65
[12] 陆秀明,黄庆,孙雪晨,等图像处理技术估测水稻叶面积
指数的研究[J].中国农学通报,2011,27(3):65-68
LuX M,HuangQ,SunX C etalApplicationofimage
processingtechnologyinricecanopyleafareaindex[J].Chinese
AgriculturalScienceBuletin,2011,27(3):65-68
[13] 鲍士旦土壤农化分析(第三版)[M].北京:中国农业出版
社,2000
BaoSDSoilandagrochemistryanalysis(3rded.)[M].
Beijing:ChineseAgriculturalPress,2000
[14] 李亚兵,毛树春,韩迎春,等不同棉花群体冠层数字图像
颜色变化特征研究[J].棉花学报,2012,24(6):541-547
LiY B,MaoSC,HanY C etalStudyonthecolor
characteristicvariationofcotoncanopydigitalimages[J].Coton
Science,2012,24(6):541-547
[15] 危常州,张福锁,朱和明,等新疆棉花氮营养诊断及追肥
推荐研究[J].中国农业科学,2002,35(12):1500-1505
WeiCZ,ZhangFS,ZhuHMetalStudyoncotonnitrogen
diagnosisandtopdressingrecommendationinnorthxinjiang[J].
ScientiaAgriculturaSinica,2002,35(12):1500-1505
[16] 白金顺,曹卫东,熊静,等应用数码相机进行绿肥翻压后
春玉米氮素营养诊断和产量预测[J].光谱学与光谱分析,
2013,33(12):3334-3338
BaiJS,CaoWD,XiongJetalNitrogenstatusdiagnosisand
yieldpredictionofspringmaizeaftergreenmanureincorporation
byusingadigitalcamera[J].SpectroscopyandSpectral
Analysis,2013,33(12):3334-3338
[17] 肖焱波,贾良良,陈新平,张福锁应用数字图像分析技术
进行冬小麦拔节期氮营养诊断[J].中国农学通报,2008,
24(8):448-453
XiaoYB,JiaLL,ChenXP,ZhangFSNstatusdiagnosisof
762
植 物 营 养 与 肥 料 学 报 21卷
winterwheatbyusingdigitalimageanalysistechnology[J].
ChineseAgriculturalScienceBuletin,2008,24(8):448-453
[18] 石媛媛基于数字图像的水稻氮磷钾营养诊断与建模研究
[D]杭州:浙江大学博士学位论文,2011
ShiYYRicenutritiondiagnosisandmodelingbasedondigital
image[D].Hangzhou:PhDdissertationofZhejiangUniversity,
  2011.
[19] 吴华兵,朱艳,田永超,等棉花冠层高光谱参数与叶片氮
含量的定量关系[J].植物生态学报,2007,31(5):903
-909
WuHB,ZhuY,TianYCetalRelationshipbetweencanopy
hyperspectraparameterandleafnitrogenconcentrationincoton
[J].JournalofPlantEcology,2007,31(5):903-909
[20] 王远,王德建,张刚基于数码相机的水稻氮素营养诊断
[J].中国农学通报,2012,28(24):111-117
WangY,WangDJ,ZhangGNitrogenstatusdiagnosisofrice
basedonadigitalcamera[J].ChineseAgriculturalScience
Buletin,2012,28(24):111-117
[21] ThomasJR,OertherGFEstimatingnitrogencontentofsweet
pepperleavesbyreflectancemeasurements[J].Agronmy
Journal,1972,27:11-13
[22] AdamsenFJ,CofeltTA,NelsonJMetalMethodforusing
imagesfromacolordigitalcameratoestimateflowernumber[J].
CropScience,2000,40:704-709
[23] RichardsonMD,KarcherDE,PurcelLCQuantifyingturfgrass
coverusingdigitalimageanalysis[J].CropScience,2001,41:
1884-1888
[24] 冯洋,陈海飞,胡孝明,等高、中、低产田水稻适宜施氮量
和氮肥利用率的研究[J].植物营养与肥料学报,2014,20
(1):7-16
FengY,ChenHF,HuXMetalOptimalnitrogenapplication
ratesonricegrainyieldandnitrogenuseeficiencyinhigh,
middleandlowyieldpaddyfields[J].JournalofPlantNutrition
andFertilizer,2014,20(1):7-16
[25] 葛梦婕,王亚江,颜希亭,等长江中下游稻区粳型超级稻
高产形成及氮素利用的研究[J].植物营养与肥料学报,
2014,20(2):259-270
GeMJ,WangYJ,YanXTetalResearchontheformationof
highyieldandthenitrogenuseofjaponicasuperriceinthe
middleandlowerareasoftheYangtzaRiver[J].JournalofPlant
NutritionandFertilizer,2014,20(2):259-270
[26] 彭少兵,黄见良,钟旭华,等提高中国稻田氮肥利用率的
研究策略[J].中国农业科学,2002,35(9):1095-1103
PengSB,HuangJL,ZhongXHetalResearchstrategyin
improvingfertilizernitrogenuseeficiencyofirigatedricein
China[J].ScientiaAgriculturaSinica,2002,35(9):1095
-1103
862