[发明专利]一种结合作物虚拟生长模型的水稻株型定量控制方法

摘要

结合作物虚拟模型的水稻株型定量化设计方法，构建作物虚拟生长模型，并实现植物主要生理过程、形态结构和光环境模型的有效关联和结合，可实现作物植株的可视化虚拟生长；在虚拟作物模型中，通过结构模型建立作物在不同生长期的形态，通过生理模型实现对作物同化物的生产、分配和最终产量的动态预测，通过光环境模型运算作物冠层光辐射量和作物个体有效光合辐射截获量；运用优化算法不断改变作物株型结构，从而得到不同的适应度模拟结果，并能获得基于不同目标的作物最优株型；将不同适应度作为目标的优化结果设置为虚拟模型的形态参数，进行模拟运算并进行横向和纵向比较，验证并确立该作物的最优株型数据。

主权项

一种结合作物虚拟生长模型的水稻株型定量控制方法，其特征在于，将虚拟植物模型和计算机优化算法有机整合，并包含以下步骤：步骤1，作物模型构建1.1首先，进行综合数据采集：实验记录目标作物整个生长期的形态动态，包括从器官到个体以及整个作物群体的几何形态特征性状数据，以及作物模型构建所需的生理数据和环境数据；1.2其次，构建作物生长模型：基于计算机图形学技术，以作物生长原理为规则基础，采用扩展的L系统建模语言XL、Java编程语言以及可进行图形替换的RGG技术，对作物器官形态和生长过程进行模拟，建立作物拓扑结构模型即包含茎、叶、穗、籽粒器官形态建成的结构模型；在结构模型的基础上，加入作物生理过程的模拟，运用LEAFC3光合作用模型模拟作物通过叶片光合作用产生同化物，根据生长模型结合源库模型模拟同化物的分配，基于beta生长函数及各器官库强度大小、获得的同化物量和所处的生长阶段计算出该器官的生长速率，从而实现整个植株在整个生长周期的动态发育；进而建立光模型，模拟太阳光源的位置和辐射度变化；虚拟的太阳光在天空中被分为直射光和散射光，包括他们在三维空间中的分布，以及在到达作物冠层后被冠层空间的遮挡，从而模拟叶片光合作用中光通量的大小；最终运用功能‑结构建模技术，结合作物结构模型、生理模型和环境模型，以时间为轴线，通过基于扩展L系统建模语言和RGG的语法规则对作物器官的形成和生长做规则化解释，并实现规则之间的迭代、替换，结合虚拟作物的同化物形成与分配，从而在结合生理功能、拓扑结构及光环境条件的基础上实现虚拟作物生长动态可视化，获得指定生育期的作物植株；1.3最后，建立优化算法模型，并在作物模型中加入优化算法模型，优化算法的步骤如下：a)确定作物一般株型特征；b)编码作物个体；c)选取优化算法适应度；d)遗传操作的判定与实现；步骤2，水稻的遗传算法具体优化2.1遗传操作的对象为水稻个体的株型，个体株型是所选株型对象特征因子的组合；在对个体进行编码时，先将个体株型的所有对象特征，即所有基因按次序排成一行，然后对每个基因进行二进制编码；在水稻株型优化中，选取不同叶位叶片的叶倾角、株高、最大分蘖数目及分蘖角度作为待优化的水稻株型对象特征，构成水稻株型因子组合；在对个体株型R进行编码时，先将该个体株型中所有株型因子按次序排成一行，记gi为第i个株型因子；2.2假设水稻个体R所含的株型个数为ng，基因gi的基因位长度记为则水稻个体编码的串长L为：L=Σi=1ngnbi]]>记gi为第i个株型因子，每个株型因子gi由二进制编码信息位组成，每个信息位均为0或1；为该株型因子第j个信息位的取值；为该株型因子的二进制码的串长，其取值由该株型因子的取值范围[ui,vi]和取值精度s共同决定；每个基因的取值范围可根据实验需求进行预设定，并根据其取值范围分别计算基因gi的基因位长度需要设置g0、g1、g2和g3的取值范围，并根据取值范围分别计算四个基因位长度，水稻个体基因位串长度L为g0、g1、g2和g3的基因位长度之和；2.3在进行初始化种群P0时，为了提高种群的多样性，采用完全随机的方法生成初始化种群中的个体，种群个体数为nc；设置最大优化代数tmax＝100；对于一个长度为LR的个体株型二进制码串，个体株型二进制码串上的每一信息位在{0,1}上随机均匀选择，所以种群规模为nc的种群初始化至少需要LR×nc次随机取值；该过程的具体步骤为：Step 1:从种群中第一个个体开始直到第nc个个体为止，对于每个个体执行Step 2；Step 2:从所选个体的第1个信息位开始到第LR个信息位为止，在每一个信息位上执行Step 3；Step 3:在所选信息位上，随机生成一个[0,1]空间内的随机数，并判断该随机数是否小于0.5，若是，则在该信息位上赋值0；否则，在该信息位上赋值1；下面给出初始化种群过程INITIALIZE(nc,LR)的伪代码：2.4在水稻模型中，对处于选定生长期的模拟个体二进制码串进行解码，获得水稻株型特征的参数组合，应用于优化过程水稻群体中各个体的参数设定；水稻个体染色体携带水稻株型特征参数值信息，不同的特征参数可以在虚拟模型中得到不同的个体形态结构；选取水稻冠层光截获量、最终产量分别作为优化过程的适应度值，水稻模型中，个体的形态结构对水稻冠层的光截获量及最终产量有显著影响，进行优化运算；并以适应度值为度量标准，对当前优化种群中的个体进行升序排列；2.5若当前优化代数t＝0，选择适应度值最大的个体作为精英个体；2.6精英个体和当前种群个体竞争，若当前种群最优个体适应度值不大于精英个体适应度值，则用此精英个体替代当前种群中的最差个体，产生第t+1代水稻种群；否则，用当前种群中的最优个体替代精英个体，同时当前种群为第t+1代水稻种群；2.7若当前优化代数大于最大优化代数，转到2.12；2.8采用轮盘赌方法从第t代水稻个体中选择出适当个体进行遗传操作：Step 1:对于第n代个体株型种群个体Rj∈P，计算其适应度值为f(Rj)；Step 2:对于个体Rj，计算该个体的适应度值在群体适应度值总和中所占的比例作为其选择概率ps(Rj)为：ps(Rj)=f(Rj)Σi=1ncf(Ri),j=1,2,...,nc]]>Step 3:对于个体Rj，计算该个体及其之前个体的选择概率之和作为该个体的累积概率pa(Rj)为：pa(Rj)=Σi=1jf(Rj)Σi=1ncf(Ri),j=1,2,...,nc]]>Step 4:在0到1范围内随机产生一个概率值；从群体中第一个个体开始，将随机概率与群体中个体的累积概率进行顺序比较，选择出个体累积概率小于或等于该随机概率的最后一个个体；2.9采用相异位交叉策略，根据两个父个体基因位值相异的位数决定是否发生交叉操作，分为以下三个步骤：Step 1:计算两个父个体株型上信息位值相异的位数m；Step 2:从父个体株型上相异信息位中随机选择m/2个位置；Step 3:从所选择位置中的第一个位置开始，在0到1范围内产生一个随机概率，判断该随机概率是否大于交叉概率pc，若是，进行父个体株型上该位置信息位的互换；否则，循环此过程直到第m/2个位置完成；最终得到新的一对个体株型；下面给出交叉操作CROSS(Chrosome c1,Chrosome c2)的伪代码：2.10采用分层变异策略，将变异操作划分为两个层次，个体层和信息位层：首先，进行个体层变异发生概率判断，从种群中选择出需要发生变异的个体；其次，在所选个体上进行信息位层变异发生概率判断，从个体中选择出需要发生变异的基因位，在基因位上进行变异操作；设置个体发生变异概率pm，分别计算水稻个体层变异概率和水稻信息位层变异概率Step 1:计算出个体层上个体株型发生变异的概率为：pm1=1-Πi=1ng(1-pm)nbi]]>Step 2:计算出基因位层上信息位发生变异的概率为：pm2=Σi=1ngnbi×pm1-Πi=1ng(1-pm)nbi]]>Step 3:从种群中的第一个个体开始，在0到1范围内产生一个随机概率，判断该随机概率是否小于等于个体株型变异概率若是，转入Step 4；否则，循环此过程环此过程直到最后一个个体；Step 4:从该株型因子中的第一个信息位开始，在0到1范围内产生一个随机概率，判断该随机概率是否小于等于信息位变异概率若是，在此信息位上进行变异操作，即将0替换为1或将1替换为0；否则，循环此过程环此过程直到该株型因子的最后一个信息位；下面给出变异操作的伪代码：2.11通过遗传操作得到临时种群，转到2.4；2.12输出当前种群中最优个体的染色体，获取该个体的水稻株型特征参数；步骤3，最优株型的比较和验证在以光截获量、籽粒数目和最终产量分别作为适应度的优化过程中，选取了三个不同的生长期，即拔节期、灌浆期和成熟期的虚拟个体作为初始种群；首先，将其优化结果分别进行纵向对比，得出基于生长期的水稻株型特征株型参数组合；其次，在构建的水稻虚拟生长模型中运用该参数组合，在其生长过程实现可视化模拟的基础上，对最终产量进行预测；同时，结合水稻模型中的同化物生产及分配算法，将前述参数组合、预测的最终产量和以产量为适应度值的优化结果进行横向对比；最终得出包含生长时间及空间维度的水稻株型优化结果，从而实现最优株型的定量控制。