一种农作物用地表面结构状况的模拟仿真方法技术

一种农作物用地表面结构状况的模拟仿真方法，用于根据不同的农业用地区域、农作物种植类型及农作物管理方式、同时考虑相关学科研究的具体要求，按照小尺度农田地表粗糙度特征参数及农作物用地周期性沟垄结构特征参数对农作物用地表面进行三维仿真模拟及对农作物地表结构粗糙度状况的二维剖面仿真模拟；同时，通过该方法还可以实现对农作物用地表面结构动态演化过程的生成。根据不同耕作方式下对农田地表状况随时间的动态变化规律的模拟结果，指导农民合理选择耕作方式，及时对农田表面进行修复，实现有效减少土壤结构的破坏和养分的流失，提高现有农田沟垄结构的水土、养分的保持能力和防风蚀能力，从而提高农作物的成活率，增加农作物的产量。

【技术实现步骤摘要】
一种农作物用地表面结构状况的模拟仿真方法
本专利技术涉及一种农作物用地表面结构状况的模拟仿真方法。
技术介绍
在农业生产中，农作物用地表面结构的形成和发展主要受到人为耕作方式的主导因素影响，同时也与土壤质地、土壤水分、风蚀过程、降水过程及沉降过程等综合自然因素密切相关；此外，农作物用地的地表径流类型与地表径流入渗过程、土壤本身的抗风蚀能力及土壤呼吸作用状况等都会对农作物表面的结构状况产生不同程度的影响。农作物用地表面结构体现在大尺度上即为农作物周期性沟垄形态的表面结构，属于定向的粗糙度形式；小尺度即为水平地表上随机形态的表面结构，属于随机的粗糙度形式。不同的耕作方式下农田地表的结构状况及其发展变化对于农作物的产量起到至关重要的作用，农田表面的粗糙度结构不同，其对光合产物的积累，排水防涝，阻风防沙，抗倒伏等功能也差异巨大。因此，精确获取表征农田地表状况的粗糙度参数，对于定量化研究及定性化描述农作物长势、合理选择农田管理方式和耕作方式、模拟农田水文过程、提高农作物遥感反演精度及土壤结构力学的研究等诸多领域的研究和应用都具有非常重要的指导和决策意义。在农田地表状况的粗糙度测量中，根据测量方式是否需要与地表直接接触，可分为接触式测量和非接触式测量两类。接触式测量主要涉及针式测量、板式测量和链式测量几种主要类型，接触式测量方法造价低廉，使用方便，但测量区域极其有限，且测量过程受地形影响大，同时会对测量区域造成一定的破坏，拍照过程和设备的损耗也会直接影响测量的精度；相比于接触式测量方法，非接触式测量主要基于激光三角法进行扫描测量和利用干涉摄影测量技术获取的立体相对实现对地表高度的测量和粗糙度的计算，虽然精度较高，但也存在设备造价昂贵，体积重量大不便于携带，处理算法复杂等缺点，普及性较差。同时，不论接触式测量方法还是非接触式测量方法，获取的地表粗糙度都是实地测量的有限区域内的特征参数，无法对大面积农田区域的地表状况进行获取和描述。现有的对大面积区域地表状况多采用模拟仿真方法，但模拟方法也均是针对多次实际测量的粗糙度结果进行线性或非线性的空间插值方法，该方法涉及的前期工作费时费力，且仅能对研究区一种特定的地表状况进行模拟，具有很大的局限性。
技术实现思路
基于以上不足之处，本专利技术的目的是提供一种农作物用地表面结构状况的模拟仿真方法，该方法实现了对农作物管理、农作物生产以及与土壤相关的学科的定性化描述及定量化研究中所涉及到的真实农作物用地表面状况参数的模拟，提供了免耕条件下地表状况的随机粗糙度参数；提供了符合不同耕作方式及形态特点的周期性沟垄结构定向粗糙度参数；同时，本专利技术还能够根据粗糙度参数随时间的变化规律，模拟不同耕种方式下农田地表状况随时间的动态过程和变化规律特征，这对于农作物的长势尤其是早期农业育苗期种床准备情况的方案选择具有重要意义。利用本专利技术涉及模拟方法所实现的农田地表粗糙度状况的提取结果，结合种植区域的土壤质地和气象条件，可以有效指导农作物的种植类型，优化农田管理模式，提高农作物的耕作质量；根据本专利技术涉及方法对农田表面结构的动态变化，还可以合理选择灌溉方式，改进作物的耕作方式，提高农田防风固沙和水土保持的能力。本专利技术所采用的技术如下：一种农作物用地表面结构状况的模拟仿真方法，步骤如下:(一).小尺度农作物用地表面结构状况的模拟仿真1.确定农田地表的模拟尺度参数，包括待模拟农田区域的长度和宽度。2.确定用于表征小尺度农作物用地表面状况的粗糙度参数，包括相关长度和均方根高度。3.根据待模拟农田区域的长度和宽度，以及输入的相关长度计算数据点的总数，按照基尔霍夫定律对这些数据点进行计算，使得这些数据点的均方根高度数据符合一定的统计分布特征，如指数分布特征或高斯分布特征。4.根据数据点的位置分布，将每个点的均方根高度数据叠加在生成的模拟地表区域上(模拟地表的初始高度为0)，从而生成小尺度农作物用地表面的随机粗糙结构特征。5.采用种子生长法，以模拟区域的任意坐标点作为种子点，生成水平或垂直剖面，以获取三维及二维剖面的均方根高度分布数据图。6.根据免耕条件下土壤粗糙度参数随时间的动态变化曲线，选择并输入表征农田地表结构状况的均方根高度和相关长度参数，实现对免耕结构下的农作物用地表面结构随时间的演化过程进行仿真模拟。(二).大尺度农作物用地表面结构状况的模拟仿真1.确定沟垄结构的形态函数曲线，提取表征大尺度农作物用地的沟垄特征表面结构参数，确定待模拟农作物用地区域的长度和宽度。根据沟垄结构参数和农田范围生成大尺度的光滑沟垄结构表面。2.确定用于表征小尺度农作物用地表面状况的粗糙度参数，包括相关长度和均方根高度。3.参考待模拟农作物用地的长度和宽度，以及输入的相关长度计算用于该模拟区域的数据点的总数，按照基尔霍夫定律对这些数据点进行计算，使得这些数据点的均方根高度数据符合一定的统计分布特征，如指数分布特征或高斯分布特征。4.根据数据点的位置分布，将每个点的均方根高度数据叠加在生成的模拟地表区域上(模拟地表的初始高度为0)作为小尺度农作物用地表面的随机粗糙结构特征，然后将小尺度的地表高度数据按照坐标位置叠加在光滑的垄结构表面，最终生成大尺度农作物用地表面的沟垄结构粗糙特征。5.采用种子生长法，以模拟区域的任意坐标点作为种子点，生成水平或垂直剖面，以获取剖面的三维及二维剖面的均方根高度分布数据图。6.根据不同耕作条件下土壤粗糙度参数随时间的动态变化曲线，选择并输入表征农田地表结构状况的均方根高度和相关长度参数，实现对具有沟垄结构的农作物用地表面结构随时间的演化过程进行仿真模拟。本专利技术的有益效果及优点：通过本专利技术涉及的地表结构状况模拟方法，不仅能够对相关学科领域内的研究提供准确的参数模拟结果，更重要的意义在于对于农作物的生产种植和管理过程中，可以合理根据本专利技术涉及的地表结构状况模拟结果，结合土壤质地等理化性质，温度、湿度等气候特点，农民的种植习惯，指导农民合理地选择农作物的种植类型；同时，根据不同耕作方式下对农田地表状况随时间的动态变化规律的模拟结果，还可以指导农民合理选择耕作方式，及时对农田表面进行修复，实现有效减少土壤结构的破坏和养分的流失，提高现有农田沟垄结构的水土、养分的保持能力和防风蚀能力，从而提高农作物的长势和育苗期作物的种床准备效率，提高作物的成活率，增加农作物的产量。由于土壤质地、农作物类型、气候条件、作物耕作方式、作物管理方式的区别，农田地表状况会随时间动态变化，本专利技术还可以根据地表粗糙度参数随时间的变化曲线，选择并输入表征农田地表结构状况的参数，实现对免耕结构的农作物用地表面结构的演化过程进行仿真模拟和对具有沟垄结构的农作物用地表面结构的演化过程进行仿真模拟附图说明图1满足高斯分布小尺度粗糙地表三维模拟结果图；图2满足高斯分布小尺度粗糙度三维模拟剖面图，单位：m；图3满足指数分布小尺度粗糙地表三维模拟结果图；图4满足指数分布小尺度粗糙度三维模拟剖面图，单位：m；图5满足高斯分布大尺度农田地表三维模拟结果图；图6满足高斯分布大尺度农田地表三维模拟剖面图，单位：m；图7满足指数分布大尺度农田地表三维模拟结果图；图8满足指数分布大尺度农田地表三维模拟剖面图，单位：m。具体实施方式下面根据附图举例对本专利技术做进一本文档来自技高网...

【技术保护点】
一种农作物用地表面结构状况的模拟仿真方法，其特征在于，步骤如下：(一)小尺度农作物用地表面结构状况的模拟仿真：(1.1)确定农田地表的模拟尺度参数，包括待模拟农田区域的长度和宽度；(1.2)确定用于表征小尺度农作物用地表面状况的粗糙度参数，包括相关长度和均方根高度；(1.3)根据待模拟农田区域的长度和宽度，以及输入的相关长度计算数据点的总数，按照基尔霍夫定律对这些数据点进行计算，使得这些数据点的均方根高度数据符合指数分布特征或高斯分布特征；(1.4)根据数据点的位置分布，将每个点的均方根高度数据叠加在生成的模拟地表区域上，模拟地表的初始高度为0，从而生成小尺度农作物用地表面的随机粗糙结构特征；(1.5)采用种子生长法，以模拟区域的任意坐标点作为种子点，生成水平或垂直剖面，以获取三维及二维剖面的均方根高度分布数据图；(1.6)根据免耕条件下土壤粗糙度参数随时间的动态变化曲线，选择并输入表征农田地表结构状况的均方根高度和相关长度参数，实现对免耕结构下的农作物用地表面结构随时间的演化过程进行仿真模拟；(二)大尺度农作物用地表面结构状况的模拟仿真：(2.1)确定沟垄结构的形态函数曲线，提取表征大尺度农作物用地的沟垄特征表面结构参数，确定待模拟农作物用地区域的长度和宽度，根据沟垄结构参数和农田范围生成大尺度的光滑沟垄结构表面；(2.2)确定用于表征小尺度农作物用地表面状况的粗糙度参数，包括相关长度和均方根高度；(2.3)参考待模拟农作物用地的长度和宽度，以及输入的相关长度计算用于该模拟区域的数据点的总数，按照基尔霍夫定律对这些数据点进行计算，使得这些数据点的均方根高度数据符合指数分布特征或高斯分布特征；(2.4)根据数据点的位置分布，将每个点的均方根高度数据叠加在生成的模拟地表区域上，模拟地表的初始高度为0，作为小尺度农作物用地表面的随机粗糙结构特征，然后将小尺度的地表高度数据按照坐标位置叠加在光滑的垄结构表面，最终生成大尺度农作物用地表面的沟垄结构粗糙特征；(2.5)采用种子生长法，以模拟区域的任意坐标点作为种子点，生成水平或垂直剖面，以获取剖面的三维及二维剖面的均方根高度分布数据图；(2.6)根据不同耕作条件下土壤粗糙度参数随时间的动态变化曲线，选择并输入表征农田地表结构状况的均方根高度和相关长度参数，实现对具有沟垄结构的农作物用地表面结构随时间的演化过程进行仿真模拟。...

【技术特征摘要】
1.一种农作物用地表面结构状况的模拟仿真方法，其特征在于，步骤如下：(一)小尺度农作物用地表面结构状况的模拟仿真：(1.1)确定农田地表的模拟尺度参数，包括待模拟农田区域的长度和宽度；(1.2)确定用于表征小尺度农作物用地表面状况的粗糙度参数，包括相关长度和均方根高度；(1.3)根据待模拟农田区域的长度和宽度，以及输入的相关长度计算数据点的总数，按照基尔霍夫定律对这些数据点进行计算，使得这些数据点的均方根高度数据符合指数分布特征或高斯分布特征；(1.4)根据数据点的位置分布，将每个点的均方根高度数据叠加在生成的模拟地表区域上，模拟地表的初始高度为0，从而生成小尺度农作物用地表面的随机粗糙结构特征；(1.5)采用种子生长法，以模拟区域的任意坐标点作为种子点，生成水平或垂直剖面，以获取三维及二维剖面的均方根高度分布数据图；(1.6)根据免耕条件下土壤粗糙度参数随时间的动态变化曲线，选择并输入表征农田地表结构状况的均方根高度和相关长度参数，实现对免耕结构下的农作物用地表面结构随时间的演化过程进行仿真模拟；(二)大尺度农作物用地表面结构状况的模拟仿真：(2.1)确定沟垄结构的形态函数曲线，提取...

【专利技术属性】
技术研发人员：任建华，张鹏，那晓东，
申请(专利权)人：哈尔滨师范大学，
类型：发明
国别省市：黑龙江,23