CN110987857A - 林下山参属性的近红外检测鉴别方法 - Google Patents

林下山参属性的近红外检测鉴别方法 Download PDF

Info

Publication number
CN110987857A
CN110987857A CN201911342449.5A CN201911342449A CN110987857A CN 110987857 A CN110987857 A CN 110987857A CN 201911342449 A CN201911342449 A CN 201911342449A CN 110987857 A CN110987857 A CN 110987857A
Authority
CN
China
Prior art keywords
sample
ginseng
prediction
forest
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201911342449.5A
Other languages
English (en)
Inventor
郑继宇
罗红
李国栋
王辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Herbpex Pharmaceutical Group Co ltd
Original Assignee
Liaoning Herbpex Pharmaceutical Group Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Liaoning Herbpex Pharmaceutical Group Co ltd filed Critical Liaoning Herbpex Pharmaceutical Group Co ltd
Priority to CN201911342449.5A priority Critical patent/CN110987857A/zh
Publication of CN110987857A publication Critical patent/CN110987857A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3563Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2135Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Data Mining & Analysis (AREA)
  • Analytical Chemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明属于林下山参检测技术领域,尤其涉及林下山参属性的近红外检测鉴别方法。本发明通过以下检测流程实现对林下山参原料及其制品的质量控制,人参属性(如种类、年份等)与近红外光谱关联模型的确立是整个鉴别工作的核心内容,本发明结合近红外光谱对林下山参、移山参、园参的种类,林下山参的年份进行鉴别,从而实现了从源头开始,严格控制原料,并对终产品进行鉴别,满足了不同客户群体的不同需求。

Description

林下山参属性的近红外检测鉴别方法
技术领域
本发明属于林下山参检测技术领域,尤其涉及林下山参属性的近红外检测鉴别方法。
背景技术
本发明背景技术中公开的信息仅仅旨在增加对本发明的总体背景的理解,而不必然被视为承认或以任何形式暗示该信息构成已经成为本领域一般技术人员所公知的现有技术。
人参,为古今临床应用最为广泛的补虚药,被称为“百草之王”,在我国有上千年的应用历史。自《神农本草经》以来,一直被视为珍贵补品,用于生津止渴和益气补虚,《伤寒论》与《金匮要略》中,含有人参的方剂达36首之多。古代人参均为野山参。
2015版《中国药典》对人参的定义主要有两种,一种是人工栽培的俗称“园参”;一种是播种在山林野生状态下自然生长的称“林下山参”,习称“籽海”。主要分布于辽宁东部、吉林东部和黑龙江东部。
园参:园参是人工种植人参,从生长环境分为两种,一种是在山上伐木、清根、修池床种植;一种是在农田种植,两者都是通过人工干预手段培育,如施肥、农药等进行种植,且最多种植6年即采挖,并上市销售;
移山参:属于园参的变种,即在园参经3-5年人工培育后,移栽到原始森林里,此种人参,在后续生长十几年后,外形比较接近林下山参,一般人员很难把其与林下山参区分开,但其功效与林下山参相差甚远,市场上经常有不法商贩,用移山参鱼目混珠,使消费者的权益受到极大的损害。
林下山参:是模拟山参的生长习性和生态环境,把参籽撒到原始森林里,自然生长,不进行任何人工干预的手段。经过几十年生长再挖出来,具有和野山参相似的品质和化学成分。
因此,虽然三种人参都在《中国药典》“人参”项下,但三者由于种植方式和生长时间的影响,其有效成分和含量各不相同,且培育的成本也不尽相同,这就导致了三者药用价值和经济价值截然不同,现在野生人参几乎绝迹。而林下山参的生长环境和方式接近于野生人参,尤其是生长年限大于15年的林下山参品质更接近于野生人参,因此也被称为“野山参”,其除了具有补五脏、安精神、定魂魄,止惊吓、除邪气、明目、开心益智等功效外,还能真正体现出人参“大补元气、复脉固脱”的传统核心功效。
我国江浙地区自古以来就有秋、冬季服用林下山参进补的传统,是我国林下山参的最大消费市场,约占全国95%的市场份额,该地区主要购买的还是经过简单处理的林下山参,保留林下山参原型态,这是因为园参、移山参、林下山参具有相同的指标成分,通过现有的检测手段如:含量测定,薄层鉴别等很难完全分辨,只能通过人参传统的鉴别方法通过产地、株龄、形状、采摘季节为初步鉴别,并配合看、闻、尝、摸、比等感官方法鉴定,这需要鉴别人员具有丰富的种参、养参及采参的经验,只有参农和部分转卖的商人能够观察出不同种类和年份的人参。同时传统鉴别方法只能保留林下山参的原型态进行鉴别。但林下山参原型态储藏极为不易,因其含有较多的糖类、黏液质和挥发油等,所以容易出现受潮、泛油、发霉、变色、虫蛀等变质现象。且因经济价值较高,会给经营者和消费者造成很大的经济损失,如果采用现代制剂技术对林下山参进行粉碎,制粒,压片。储藏和服用虽然方便,但在失去林下山参原形态的前提下,如何取得消费者的认可又是一项技术难题。
发明内容
针对上述的问题,本发明通过以下检测流程实现对林下山参原料及其制品的质量控制,人参属性(如种类、年份等)与近红外光谱关联模型的确立是整个鉴别工作的核心内容,为了实现鉴别目的,本发明提供了两种模式识别方法。
首先,本发明公开第一种林下山参属性的近红外检测鉴别方法,采用主成分判别分析(PCADA)建模,具体步骤为:
(1)获取样本的近红外光谱数据矩阵X,其大小为n*p,n代表样本数量, p代表变量数目,对应于近红外光谱就是检测波长。将n个样本光谱按浓度从小到大顺序等间距的提取预测样本光谱矩阵Xval,剩余的就是建模的样本 Xcal,优选地,比例为1:3。
(2)对步骤(1)得到的光谱进行预处理:即先对每条光谱进行S-G平滑,再对其求导;以便于用于后续建立模型。
(3)将用于建模的样本Xcal进行奇异值分解,Xcal=USVt,U*S为得分矩阵Tcal,Vt为载荷矩阵的转置。
(4)用预测样本光谱矩阵Xval乘上步骤(3)中建模样本的载荷矩阵Vt,得到预测样本的得分矩阵Tval,即Tval=Xval×V。
(5)将预测样本的得分矩阵Tval前n列与某一类建模样本得分矩阵Tcal 的前n列计算马氏距离,计算公式如式(1)所示:
Figure RE-GDA0002391509600000031
式中ti为预测样本的得分矩阵中的第i列,为同一类拟合样本得分矩阵的均值, V-1为协方差矩阵的逆矩阵。
(6)对于两类判别问题,对同一个样本计算得到两个马氏距离,代表此样本到两类中心的马氏距离。比较两个马氏距离之间的大小,将此样本归属为距离较小的一类。
进一步地,所述第一种林下山参属性的近红外检测鉴别方法中,建模完成后,还包括对建立的模型的预测能力的衡量的步骤。
可选地,用总体正判率(RMSEC)来衡量模型的预测能力,即:总体正判率=判断正确的预测样本数/预测样本总数。
其次,本发明公开第二种林下山参属性的近红外检测鉴别方法,采用偏最小二乘判别分析(PCADA)建模,具体步骤为:
(S1)获取样本的近红外光谱数据矩阵X,其大小为n*p,n代表样本数量,p 代表变量数目,对于近红外光谱就是检测波长;获取样本的标签信息向量y,大小为n*1。
(S2)对步骤(S1)得到的光谱进行预处理:即先对每条光谱进行S-G平滑,再对其求导。
(S3)采用pls1算法计算建模系数向量b,其算法内容如下:取向量y作为t;计w=(Xty)/(yty);;对w进行归一化;t=Xw;计算t对y的回归系数,v=(tty)/(ttt); p=Xtt/(ttt);计算残差X=X-tpt,y=y-vt;返回开始计算下个潜变量。当循环完成后,b=w/(ptw)v。
(S4)计算预测样本的估计值
Figure RE-GDA0002391509600000041
(S5)对每个预测样本的估计值与两类的标签信息进行比较,更接近其中某个值就将该样本判定为那一类。
进一步地,步骤(S1)中,由于本发明建立的是两类样本的定性判别模型,其标签信息是人为定义,一类用1,另一类用-1代替,其代表了样本的类别;即步骤(S5)中,所述比较的标签信息为1与-1。
进一步地,所述第二种林下山参属性的近红外检测鉴别方法中,建模完成后,还包括对建立的模型的预测能力的衡量的步骤。
可选地,用总体正判率与预测均方差(RMSEP)来衡量模型的预测能力,计算公式如式(2)所示:
Figure RE-GDA0002391509600000051
式中
Figure RE-GDA0002391509600000052
为第i个预测样本的估计值,yi为第i个预测样本的真值。
与现有技术相比,本发明取得了以下有益效果:本发明结合近红外光谱对林下山参、移山参、园参的种类,林下山参的年份进行鉴别,从而实现了从源头开始,严格控制原料,并对终产品进行鉴别,鉴别结果显示:本发明提出的方法具有较高的正判率,且验证集的预测均方差(RMSEP)足够小,能够满足不同客户群体的不同需求。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
正如背景技术所述,移山参在后续生长十几年后,外形比较接近林下山参,一般人员很难把其与林下山参区分开,但其功效与林下山参相差甚远,市场上经常有不法商贩,用移山参鱼目混珠,使消费者的权益受到极大的损害。因此,本发明提出了两种林下山参属性的近红外检测鉴别方法;现结合具体实施方式对本发明进一步说明。
实施例
林下山参近红外定性判别和年份判别模型的建立与评价:因为原型态的林下山参和普通园参在外观形态上差别很大,任何人都可以轻易加以区分,所以本实施例的山参定性鉴别只针对林下山参和移山参。
(一)样品采集:首先连接仪器,打开软件,确认连接后,打开温控开关,调整为手动模式,设置相关参数,温度控制为25℃,积分时间8ms,扫描次数为 32次。测量暗光谱,随后打开灯源等待10分钟,以装置背景为参比,即可开始进行测样。采用漫反射发收集近红外光谱,测量范围是5000cm-1~12000cm-1,每根人参在芦头、参体、参须上选择8个点,其中参芦2个点,参体4 个点,参须2个点,每个点测量3次,具体测量数据如表1所示。
表1人参数据详表
Figure RE-GDA0002391509600000061
(二)林下山参与移山参定性鉴别的近红外模型建立:
重复测量的3次光谱用平均值代替,将这些光谱按3:1划分建模集与独立预测集,拿得到的建模集光谱按2:1划分校正集和验证集进行蒙特卡洛交互检验,重复1000次,得到不同潜变量下判别的结果,从中挑选合适的潜变量进行建模。为了区分不同预处理方法建立模型的稳定性,在挑选的潜变量下分别建模100次并用独立预测集中的样本进行验证。不同预处理方法建立的模型效果如表2所示。
表2人参种类判别PLSDA模型参数汇总
Figure RE-GDA0002391509600000071
用上述确定的模型参数建立模型进行移山参与林下参的种类判别。从表2 可以看出:不经过原始数据光谱已经可以达到较高的正判率,且验证集的预测均方差(RMSEP)足够小。故选择原始数据在17个潜变量下对林下参与移山参进行种类判别。
(三)林下山参年份鉴别的近红外模型的建立:
对重复测量的3次光谱用平均值代替,将这些光谱按3:1划分建模集与独立预测集,拿得到的建模集光谱按2:1划分校正集和验证集进行蒙特卡洛交互检验,重复1000次,得到不同潜变量下判别的结果,从中挑选合适的潜变量进行建模。为了区分不同预处理方法建立模型的稳定性,在挑选的潜变量下分别建模100次并用独立预测集中的样本进行验证。不同预处理方法建立的模型效果如表3所示。
表3林下山参年份判别PLSDA模型参数汇总
Figure RE-GDA0002391509600000081
Figure RE-GDA0002391509600000091
由于是三种年份之间的判别,故先将两种年份相近的作为一类与另外的一种进行判别,再将剩下的两类年份判别区分开。通过表3可以看出:原始光谱两次判别综合效果最佳,且能正确区分不同年份的林下参。
以上所述仅为本申请的优选实施例,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (7)

1.一种林下山参属性的近红外检测鉴别方法,其特征在于:采用主成分判别分析建模,具体步骤为:
(1)获取样本的近红外光谱数据矩阵X,其大小为n*p,n代表样本数量,p代表变量数目,对应于近红外光谱就是检测波长;将n个样本光谱按浓度从小到大顺序等间距的提取预测样本光谱矩阵Xval,剩余的就是建模的样本Xcal;
(2)对步骤(1)得到的光谱进行预处理:即先对每条光谱进行S-G平滑,再对其求导;
(3)将用于建模的样本Xcal进行奇异值分解,Xcal=USVt,U*S为得分矩阵Tcal,Vt为载荷矩阵的转置;
(4)用预测样本光谱矩阵Xval乘上步骤(3)中建模样本的载荷矩阵Vt,得到预测样本的得分矩阵Tval,即Tval=Xval×V;
(5)将预测样本的得分矩阵Tval前n列与某一类建模样本得分矩阵Tcal的前n列计算马氏距离,计算公式如下式所示:
Figure FDA0002331824360000011
式中ti为预测样本的得分矩阵中的第i列,为同一类拟合样本得分矩阵的均值,V-1为协方差矩阵的逆矩阵;
(6)对于两类判别问题,对同一个样本计算得到两个马氏距离,代表此样本到两类中心的马氏距离;比较两个马氏距离之间的大小,将此样本归属为距离较小的一类。
2.如权利要求1所述的林下山参属性的近红外检测鉴别方法,其特征在于:建模完成后,还包括对建立的模型的预测能力的衡量的步骤。
3.如权利要求2所述的林下山参属性的近红外检测鉴别方法,其特征在于:用总体正判率来衡量模型的预测能力;优选地,所述总体正判率=判断正确的预测样本数/预测样本总数。
4.一种林下山参属性的近红外检测鉴别方法,其特征在于:采用偏最小二乘判别分析建模,具体步骤为:
(S1)获取样本的近红外光谱数据矩阵X,其大小为n*p,n代表样本数量,p代表变量数目,对于近红外光谱就是检测波长;获取样本的标签信息向量y,大小为n*1;
(S2)对步骤(S1)得到的光谱进行预处理:即先对每条光谱进行S-G平滑,再对其求导;
(S3)采用pls1算法计算建模系数向量b,其算法内容如下:取向量y作为t;计算w=Xty/(yty);对w进行归一化;t=Xw;计算t对y的回归系数,υ=(tty)/(ttt);p=Xtt/(ttt);计算残差:X=X-tpt,y=y-υt;返回开始计算下个潜变量;当循环完成后,b=w/(ptw)v;
(S4)计算预测样本的估计值,
Figure FDA0002331824360000021
(S5)对每个预测样本的估计值与两类的标签信息进行比较,更接近其中某个值就将该样本判定为那一类。
5.如权利要求4所述的林下山参属性的近红外检测鉴别方法,其特征在于:步骤(S1)中,所述标签信息定义谓:一类用1,另一类用-1代替。
6.如权利要求4或5所述的林下山参属性的近红外检测鉴别方法,其特征在于:建模完成后,还包括对建立的模型的预测能力的衡量的步骤。
7.如权利要求6所述的林下山参属性的近红外检测鉴别方法,其特征在于:用总体正判率与预测均方差来衡量模型的预测能力,计算公式为:
Figure FDA0002331824360000031
式中
Figure FDA0002331824360000032
为第i个预测样本的估计值,yi为第i个预测样本的真值。
CN201911342449.5A 2019-12-23 2019-12-23 林下山参属性的近红外检测鉴别方法 Pending CN110987857A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911342449.5A CN110987857A (zh) 2019-12-23 2019-12-23 林下山参属性的近红外检测鉴别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911342449.5A CN110987857A (zh) 2019-12-23 2019-12-23 林下山参属性的近红外检测鉴别方法

Publications (1)

Publication Number Publication Date
CN110987857A true CN110987857A (zh) 2020-04-10

Family

ID=70075886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911342449.5A Pending CN110987857A (zh) 2019-12-23 2019-12-23 林下山参属性的近红外检测鉴别方法

Country Status (1)

Country Link
CN (1) CN110987857A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112113922A (zh) * 2020-09-16 2020-12-22 中国中医科学院中药研究所 基于高光谱成像技术的人参年限识别方法
CN113536927A (zh) * 2021-06-15 2021-10-22 南昌海关技术中心 一种赣南脐橙鉴别方法、系统及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636452A (zh) * 2012-05-03 2012-08-15 中国科学院长春光学精密机械与物理研究所 一种近红外光谱无损鉴别山参真伪的方法
CN103674854A (zh) * 2013-12-05 2014-03-26 江苏大学 一种基于联合区间偏最小二乘判别分析的苹果分级方法
CN103776777A (zh) * 2014-01-09 2014-05-07 北京中医药大学 一种用近红外光谱技术识别不同生长方式人参及对人参中组分含量测定的方法
CN104730030A (zh) * 2015-03-27 2015-06-24 山东大学 基于近红外分析技术的党参真伪鉴别和产地判定的方法
CN109632696A (zh) * 2018-12-13 2019-04-16 宜宾学院 一种识别药物片剂来源的低成本近红外光谱方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636452A (zh) * 2012-05-03 2012-08-15 中国科学院长春光学精密机械与物理研究所 一种近红外光谱无损鉴别山参真伪的方法
CN103674854A (zh) * 2013-12-05 2014-03-26 江苏大学 一种基于联合区间偏最小二乘判别分析的苹果分级方法
CN103776777A (zh) * 2014-01-09 2014-05-07 北京中医药大学 一种用近红外光谱技术识别不同生长方式人参及对人参中组分含量测定的方法
CN104730030A (zh) * 2015-03-27 2015-06-24 山东大学 基于近红外分析技术的党参真伪鉴别和产地判定的方法
CN109632696A (zh) * 2018-12-13 2019-04-16 宜宾学院 一种识别药物片剂来源的低成本近红外光谱方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陆婉珍 等, 哈尔滨工程大学出版社 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112113922A (zh) * 2020-09-16 2020-12-22 中国中医科学院中药研究所 基于高光谱成像技术的人参年限识别方法
CN113536927A (zh) * 2021-06-15 2021-10-22 南昌海关技术中心 一种赣南脐橙鉴别方法、系统及存储介质

Similar Documents

Publication Publication Date Title
Creasy et al. Grapes
Nyasordzi et al. Utilizing the IAD index to determine internal quality attributes of apples at harvest and after storage
CN104838905B (zh) 西瓜苗期抗旱性的评价方法
Kızıl et al. Lettuce (Lactuca sativa L.) yield prediction under water stress using artificial neural network (ANN) model and vegetation indices.
Nordey et al. Spatial and temporal variations in mango colour, acidity, and sweetness in relation to temperature and ethylene gradients within the fruit
CN110987857A (zh) 林下山参属性的近红外检测鉴别方法
Brancourt‐Hulmel et al. A diagnosis of yield‐limiting factors on probe genotypes for characterizing environments in winter wheat trials
Jiang et al. Quantitative analysis of tomato yield and comprehensive fruit quality in response to deficit irrigation at different growth stages
Farooq et al. Dynamics of canopy development of Cunninghamia lanceolata mid-age plantation in relation to foliar nitrogen and soil quality influenced by stand density
Olivoto et al. Multivariate analysis of strawberry experiments: where are we now and where can we go?
Franco-Hermida et al. Validating integrative nutrient diagnostic norms for greenhouse cut-roses
Guo et al. Study on life histroy traits of Stellera chamaejasme provide insights into its control on degraded typical steppe
Ramirez et al. Calibration and validation of complex and simplified tomato growth models for control purposes in the southeast of Spain
Albrizio et al. Water stress, yield, and grape quality in a hilly rainfed “Aglianico” vineyard grown in two different soils along a slope
Wubs et al. Quantifying abortion rates of reproductive organs and effects of contributing factors using time-to-event analysis
Felts et al. Physiochemical and descriptive sensory analysis of Arkansas muscadine grapes
Dallabetta et al. The influence of training system on apple fruit quality
Najafzadeh et al. Superior growth characteristics, yield, and fruit quality in promising european pear (Pyrus communis L.) chance seedlings in Iran
CN109596577A (zh) 小麦白粉病病情监测模型的构建方法及宽角度适应的监测方法
CN105823767B (zh) 一种评价海带高温耐受性的方法
CN113470010A (zh) 一种基于机器学习的火龙果种植全周期病虫害防治方法
Pallas et al. X-Palm, a functional structural plant model for analysing temporal, genotypic and inter-tree variability of oil palm growth and yield
Grüneberg et al. Procedures for the evaluation and analysis of sweetpotato trials
Navarrete Characterizing grapevine canopy architecture
ZHANG et al. Nutrient deficiency diagnosis in whole hydroponic lettuce based on random forest.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200410