一种结合作物虚拟生长模型的水稻株型定量控制方法
【专利摘要】结合作物虚拟模型的水稻株型定量化设计方法，构建作物虚拟生长模型，并实现植物主要生理过程、形态结构和光环境模型的有效关联和结合，可实现作物植株的可视化虚拟生长；在虚拟作物模型中，通过结构模型建立作物在不同生长期的形态，通过生理模型实现对作物同化物的生产、分配和最终产量的动态预测，通过光环境模型运算作物冠层光辐射量和作物个体有效光合辐射截获量；运用优化算法不断改变作物株型结构，从而得到不同的适应度模拟结果，并能获得基于不同目标的作物最优株型；将不同适应度作为目标的优化结果设置为虚拟模型的形态参数，进行模拟运算并进行横向和纵向比较，验证并确立该作物的最优株型数据。
【专利说明】一种结合作物虚拟生长模型的水稻株型定量控制方法
【技术领域】
[0001] 本
【发明内容】
涉及植物虚拟建模，植物生理生态学以及优化算法设计领域，是一种 能够实现作物株型自动优化的作物株型控制方法，可运用于作物株型的定量化控制与设 计。
【背景技术】
[0002] 传统育种经过长时间的发展，已形成了一整套较为完善的育种理论和技术体系， 并在实践中取得了巨大的成就。尤其是80年代以来，现代生物技术的迅猛发展及其与传统 育种技术的结合，在分子层面上促进了传统育种的发展。然而，以分子生物学技术为手段辅 助传统育种设计优良株型在目前并不现实。因此目前的作物高产育种研究需要研究人员不 断进行田间栽培试验来获取经验和数据，并以此为主要手段来确认和评价育成的品种、品 系在生产上的价值。因此，作物理想株型的选择和获取需要较长周期的田间试验，尤其是在 不确定的培育目标的情况下，导致培育出的株型缺乏针对性，从而增加了育种成本。
[0003] 作为数字农业的核心和关键技术之一，虚拟植物生长建模的研究和应用在国内外 都受到了广泛的重视。虚拟植物可以在计算机上以三维图形的形式可视化展现植物的形态 形成及生长过程，并结合作物生理生态过程，实现作物特定生理过程及器官的生长动态特 征的调控和研究，从而对作物的空间结构、主要生产环节的管理方案进行动态调控，为农林 业生产服务。
[0004] 以虚拟作物为基础的作物株型优化方法，结合了作物功能结构模型与遗传算法优 化思想，能够实现优化培育目标的定量控制，确定最优株型特征参数，从而利用现代计算 机技术辅助传统育种手段，实现节约资源、降低成本和提高作物产量的目的。
【发明内容】
[0005] 为了进一步对传统作物育种和现代生物工程结合的育种方法进行辅助，弥补一般 作物株型设计流程中实验周期长、针对性差、人工操作繁琐等的不足，本发明提供一种水稻 株型定量控制方法，该方法能够自动而快速地对作物的理想株型作出选择，从而实现针对 性地对株型进行优化，在育种过程中起到辅助作用。
[0006] 本发明解决株型优化与控制的技术问题所采用的技术方案是：
[0007] -种结合作物虚拟生长模型的株型定量控制方法，其特征在于，将虚拟植物模型 和计算机优化算法有机整合，并包含以下步骤：
[0008] 步骤1，作物模型构建
[0009] I. 1首先，进行综合数据采集：实验记录目标作物整个生长期的形态动态，包括从 器官到个体以及整个作物群体的几何形态特征性状数据，以及作物模型构建所需的生理数 据和环境数据（如光照量等）；
[0010] 1. 2其次，构建作物生长模型：基于计算机图形学技术，以作物生长原理为规则基 础，采用扩展的L系统建模语言XUJava编程语言以及可进行图形替换的RGG(Relational Growth Grammar)技术，对作物器官形态和生长过程进行模拟，建立作物拓扑结构模型（以 水稻为例，即包含茎、叶、穗、籽粒等器官形态建成的结构模型）；在结构模型的基础上，力口 入作物生理过程的模拟，运用LEAFC3光合作用模型模拟作物通过叶片光合作用产生同化 物，根据生长模型结合源库模型模拟同化物的分配，基于beta生长函数及各器官库强度大 小、获得的同化物量和所处的生长阶段计算出该器官的生长速率，从而实现整个植株在整 个生长周期的发育动态；进而建立光模型，模拟太阳光源的位置和辐射度变化；虚拟的太 阳光在天空中被分为直射光和散射光，包括他们在三维空间中的分布，以及在到达作物冠 层后被冠层空间的遮挡等，从而实现叶片光合作用中光通量的大小；最终运用功能-结构 建模技术，结合作物结构模型、生理模型和环境模型，以时间为轴线，通过基于扩展L系统 建模语言和RGG的语法规则对作物器官的形成和生长做规则化解释，并实现规则之间的迭 代、替换，结合虚拟作物的同化物形成与分配，从而在结合生理功能、拓扑结构及光环境条 件的基础上实现虚拟作物生长动态可视化，获得指定生育期的作物植株；
[0011] 1.3最后，建立优化算法模型，并在作物模型中加入优化算法模型，优化算法的一 般步骤如下：
[0012] a)确定作物一般株型特征；
[0013] b)编码作物个体；
[0014] c)选取优化算法适应度；
[0015] d)遗传操作的判定与实现；
[0016] 步骤2,以水稻为例的遗传算法具体优化
[0017] 2. 1遗传操作的对象为水稻个体的株型，个体株型是所选株型对象特征因子的组 合；在对个体进行编码时，先将个体株型的所有对象特征，即所有基因按次序排成一行，然 后对每个基因进行二进制编码；
[0018] 在水稻株型优化中，选取不同叶位叶片的叶倾角、株高、最大分蘖数目及分蘖角度 作为待优化的水稻株型对象特征，构成水稻株型因子组合；在对个体株型R进行编码时，先 将该个体株型中所有株型因子按次序排成一行，记 gi为第i个株型因子；
[0019] 2. 2假设水稻个体R所含的株型个数为ng，基因gi的基因位长度记为< 则水稻 个体编码的串长L为：
[0020]
【权利要求】
1. 一种结合作物虚拟生长模型的水稻株型定量控制方法，其特征在于，将虚拟植物模 型和计算机优化算法有机整合，并包含以下步骤： 步骤1，作物模型构建 1. 1首先，进行综合数据采集：实验记录目标作物整个生长期的形态动态，包括从器官 到个体以及整个作物群体的几何形态特征性状数据，以及作物模型构建所需的生理数据和 环境数据（如光照量等）； 1. 2其次，构建作物生长模型：基于计算机图形学技术，以作物生长原理为规则基础， 采用扩展的L系统建模语言XL、Java编程语言以及可进行图形替换的RGG(Relational Growth Grammar)技术，对作物器官形态和生长过程进行模拟，建立作物拓扑结构模型（以 水稻为例，即包含茎、叶、穗、籽粒等器官形态建成的结构模型）；在结构模型的基础上，力口 入作物生理过程的模拟，运用LEAFC3光合作用模型模拟作物通过叶片光合作用产生同化 物，根据生长模型结合源库模型模拟同化物的分配，基于beta生长函数及各器官库强度大 小、获得的同化物量和所处的生长阶段计算出该器官的生长速率，从而实现整个植株在整 个生长周期的发育动态；进而建立光模型，模拟太阳光源的位置和辐射度变化；虚拟的太 阳光在天空中被分为直射光和散射光，包括他们在三维空间中的分布，以及在到达作物冠 层后被冠层空间的遮挡等，从而实现叶片光合作用中光通量的大小；最终运用功能-结构 建模技术，结合作物结构模型、生理模型和环境模型，以时间为轴线，通过基于扩展L系统 建模语言和RGG的语法规则对作物器官的形成和生长做规则化解释，并实现规则之间的迭 代、替换，结合虚拟作物的同化物形成与分配，从而在结合生理功能、拓扑结构及光环境条 件的基础上实现虚拟作物生长动态可视化，获得指定生育期的作物植株； 1. 3最后，建立优化算法模型，并在作物模型中加入优化算法模型，优化算法的一般步 骤如下： a) 确定作物一般株型特征； b) 编码作物个体； c) 选取优化算法适应度； d) 遗传操作的判定与实现； 步骤2,以水稻为例的遗传算法具体优化 2. 1遗传操作的对象为水稻个体的株型，个体株型是所选株型对象特征因子的组合； 在对个体进行编码时，先将个体株型的所有对象特征，即所有基因按次序排成一行，然后对 每个基因进行二进制编码； 在水稻株型优化中，选取不同叶位叶片的叶倾角、株高、最大分蘖数目及分蘖角度作为 待优化的水稻株型对象特征，构成水稻株型因子组合；在对个体株型R进行编码时，先将该 个体株型中所有株型因子按次序排成一行，记 gi为第i个株型因子； 2. 2假设水稻个体R所含的株型个数为ng，基因 gi的基因位长度记为4，则水稻个体编 码的串长L为：
记1 =也，於,…，，& =涔灰，·.·，￥); gi为第i个株型因子，每个株型因子gi由 二进制编码信息位组成，每个信息位均为O或1 ; Ai为该株型因子第]_个信息位的取值；4 为该株型因子的二进制码的串长，其取值由该株型因子的取值范围[Ui, Vi]和取值精度S共 同决定；每个基因的取值范围可根据实验需求进行预设定，并根据其取值范围分别计算基 因&的基因位长度 需要设置go、gl、g2和g3的取值范围，并根据取值范围分别计算四个基因位长度，水稻 个体基因位串长度L为g(l、gl、g2和g3的基因位长度之和； 2. 3在进行初始化种群Ptl时，为了提高种群的多样性，采用完全随机的方法生成初始化 种群中的个体，种群个体数为η。；设置最大优化代数t max = 100 ;对于一个长度为Lk的个体 株型二进制码串，个体株型二进制码串上的每一信息位在{〇，1}上随机均匀选择，所以种 群规模为η。的种群初始化至少需要L KXn。次随机取值；该过程的具体步骤为： Stepl:从种群中第一个个体开始直到第η。个个体为止，对于每个个体执行Step2 ; Step2:从所选个体的第1个信息位开始到第Lk个信息位为止，在每一个信息位上执行 Step3 ； St印3:在所选信息位上，随机生成一个[0, 1]空间内的随机数，并判断该随机数是否 小于〇. 5,若是，则在该信息位上赋值0 ;否则，在该信息位上赋值1 ; 下面给出初始化种群过程INITIALIZE (n。，Lk)的伪代码：
2. 4在水稻模型中，对处于选定生长期的模拟个体二进制码串进行解码，获得水稻株型 特征的参数组合，应用于优化过程水稻群体中各个体的参数设定；水稻个体染色体携带水 稻株型特征参数值信息，不同的特征参数可以在虚拟模型中得到不同的个体形态结构；选 取水稻冠层光截获量、最终产量分别作为优化过程的适应度值（水稻模型中，个体的形态 结构对水稻冠层的光截获量及最终产量有显著影响），进行优化运算； 并以适应度值为度量标准，对当前优化种群中的个体进行升序排列； 2. 5若当前优化代数t = 0,选择适应度值最大的个体作为精英个体； 2. 6精英个体和当前种群个体竞争，若当前种群最优个体适应度值不大于精英个体适 应度值，则用此精英个体替代当前种群中的最差个体，产生第t+Ι代水稻种群；否则，用当 前种群中的最优个体替代精英个体，同时当前种群为第t+1代水稻种群； 2. 7若当前优化代数大于最大优化代数，转到2. 12 ; 2. 8采用轮盘赌方法从第t代水稻个体中选择出适当个体进行遗传操作： St印1:对于第η代个体株型种群氧，}，个体& e P，计算其适应度值为 f(Rj)； Step2:对于个体&，计算该个体的适应度值在群体适应度值总和中所占的比例作为其 选择概率Ps (Rj)为：
Step3:对于个体&，计算该个体及其之前个体的选择概率之和作为该个体的累积概率 Pa(Rj)为：
Step4:在O到1范围内随机产生一个概率值。从群体中第一个个体开始，将随机概率 与群体中个体的累积概率进行顺序比较，选择出个体累积概率小于或等于该随机概率的最 后一个个体； 2. 9采用相异位交叉策略，根据两个父个体基因位值相异的位数决定是否发生交叉操 作，一般分为以下三个步骤： Stepl:计算两个父个体株型上信息位值相异的位数m ; Step2:从父个体株型上相异信息位中随机选择m/2个位置； Step3:从所选择位置中的第一个位置开始，在O到1范围内产生一个随机概率，判断该 随机概率是否大于交叉概率P。，若是，进行父个体株型上该位置信息位的互换；否则，循环 此过程直到第m/2个位置完成。最终得到新的一对个体株型； 下面给出交叉操作CROSS(Chrosome cl，Chrosome c2)的伪代码：
2. 10采用分层变异策略，将变异操作划分为两个层次（个体层和信息位层）：首先，进 行个体层变异发生概率判断，从种群中选择出需要发生变异的个体；其次，在所选个体上进 行信息位层变异发生概率判断，从个体中选择出需要发生变异的基因位，在基因位上进行 变异操作；设置个体发生变异概率Pm，分别计算水稻个体层变异概率Pi和水稻信息位层变 异概率 St印1:计算出个体层上个体株型发生变异的概率A1w为：
Step2:计算出基因位层上信息位发生变异的概率为：
Step3:从种群中的第一个个体开始，在0到1范围内产生一个随机概率，判断该随机概 率是否小于等于个体株型变异概率pi，若是，转入Step4 ;否则，循环此过程环此过程直到 最后一个个体； Step4:从该株型因子中的第一个信息位开始，在0到1范围内产生一个随机概率，判断 该随机概率是否小于等于信息位变异概率，若是，在此信息位上进行变异操作，即将〇替 换为1或将1替换为〇 ;否则，循环此过程环此过程直到该株型因子的最后一个信息位； 下面给出变异操作 Mutation (Chrosome []c, int nc, int ng, int ?I )的伪代码：
2. 11通过遗传操作得到临时种群，转到2. 4 ; 2. 12输出当前种群中最优个体的染色体，获取该个体的水稻株型特征参数； 步骤3,最优株型的比较和验证 在以光截获量、籽粒数目和最终产量分别作为适应度的优化过程中，选取了三个不同 的生长期（拔节期、灌浆期和成熟期）的虚拟个体作为初始种群；首先，将其优化结果分别 进行纵向对比，得出基于生长期的水稻株型特征株型参数组合；其次，在构建的水稻虚拟生 长模型中运用该参数组合，在其生长过程实现可视化模拟的基础上，对最终产量进行预测； 同时，结合水稻模型中的同化物生产及分配算法，将前述参数组合、预测的最终产量和以产 量为适应度值的优化结果进行横向对比；最终得出包含生长时间及空间维度的水稻株型优 化结果，从而实现最优株型的定量控制。
2. 如权利要求1所述的结合作物虚拟模型的水稻株型定量化设计方法，其特征在于： 所述步骤1中，基于作物株型结构的形态一生理一环境的相互关联，能够通过数值改变作 物株型结构，进而实现对目标适应度值的模拟。
3. 如权利要求1所述的结合作物虚拟模型的水稻株型定量化设计方法，其特征在于： 所述步骤2中，将作物株型设计问题转化为数值优化问题，在计算机优化方法内改变作物 株型结构特征值，并将其在虚拟模型中实现可视化模拟，以模拟植株的不同输出为目标进 行数值比较。
4. 如权利要求1所述的结合作物虚拟模型的水稻株型定量化设计方法，其特征在于： 所述步骤2中，运用相应的方法对优化群体进行特征相关的遗传操作：用轮盘赌方法进行 选择操作；用相异位交叉策略实现交叉操作；用分层变异策略实现变异操作，从而增加优 化群体的多态性。
5. 如权利要求1所述的结合作物虚拟模型的株型定量化设计方法，其特征在于：所述 步骤3中，先运用优化算法以不同的目标适应度值计算不同时间、空间及不同模拟目标的 最优株型特征参数；然后将这些参数组合返回设定为虚拟模型的参数值并重新进行可视化 模拟；最后将模拟结果进行对比，得到包含时间、空间维度的优化结果。
【文档编号】G06T17/00GK104376142SQ201410357630
【公开日】2015年2月25日 申请日期:2014年7月25日 优先权日:2014年7月25日
【发明者】徐利锋, 丁维龙, 危扬, 陈淑娇, 刘洋, 郑蕾, 程志君 申请人:浙江工业大学