一种基于多源信息的实际灌溉面积提取方法与流程

本发明涉及农业水土工程、遥感
技术领域：
，特别是涉及一种基于多源信息的实际灌溉面积提取方法。
背景技术：
：我国是一个水资源极其匮乏的国家，加之水环境污染和水生态损害愈加凸显，水资源供需矛盾日益突出。在我国，农业是用水大户，农业用水量约占总用水量的62％，而农田灌溉用水量约占农业用水量的85％，是我国经济社会用水的重要组成部分，因此，农业用水管理也成为最严格水资源管理的重点和核心内容。截至2016年底，全国灌溉面积达到10.98亿亩，全国万亩以上灌区7806处，其中30万亩以上大型灌区458处，耕地灌溉面积2.66亿亩，占全国耕地面积的13％，灌区内生产的粮食产量、农业总产值均超过全国总量的1/4，是我国粮食安全的重要保障和农业农村经济社会发展的重要支撑。加强灌区用水管理、提高灌溉水的利用效率，是缓解我国水资源紧缺的迫切要求。灌溉面积是灌区有效实施用水总量控制和定额管理的重要用水参数之一。目前，灌溉面积主要通过调查统计方式获取，时效性和准确性相对较差，且统计工作量较大；而在灌区用水管理中迫切需要掌握灌溉面积指标，及其空间分布和动态变化等信息，因此，传统灌溉面积获取方法已不能满足灌区现代化管理需求，亟需建立更为科学、合理的相关技术和方法，尤其在基于遥感等多源信息的高分辨率实际灌溉面积快速提取方面加强研究，为农业用水管理提供科学、客观依据。技术实现要素：有鉴于此，本发明提供了一种基于多源信息的实际灌溉面积提取方法，用于对实际灌溉面积进行快速、有效提取，为农业用水管理提供科学、客观的依据。为了解决上述问题，本发明公开了一种基于多源信息的实际灌溉面积提取方法，包括步骤：根据预先获取的多源信息提取潜在灌溉面积，所述多源信息包括基础地理信息、气象信息、灌排工程信息、农业种植结构信息、灌溉制度信息；利用遥感信息对所述潜在灌溉面积内作物生育期的田间表层土壤含水量进行反演计算；根据所述田间表层土壤含水量，采用不同方法从所述潜在灌溉面积中进行提取，得到实际灌溉面积。可选的，所述根据多源信息提取潜在灌溉面积，包括步骤：获取所述基础地理信息；提取所述农业种植结构信息；获取气象信息，生成全国降水量分布图；建立主要作物灌溉制度数据库，对所述灌溉制度信息进行数据分析；根据所述基础地理信息、农业种植结构信息和对所述灌溉制度数据分析结果，提取所述潜在灌溉面积。可选的，所述基础地理信息包括行政区划边界、土地利用图、全国地貌图和土壤类型分布图。可选的，所述潜在灌溉面积包括耕地灌溉面积和非耕地灌溉面积。可选的，所述根据多源信息提取潜在灌溉面积，还包括步骤：提取灌区或行政区的潜在灌溉面积。可选的，所述根据多源信息提取潜在灌溉面积，还包括步骤：对灌溉地和非灌溉地进行选取。可选的，所述利用遥感信息对所述潜在灌溉面积内作物生育期的田间表层土壤含水量进行反演计算，包括步骤：对多源遥感影像进行预处理；通过经过预处理的多源遥感影像提取归一化差值植被指数；通过经过预处理的多源遥感影像和大气剖面参数计算地表温度；根据所述归一化差值植被指数和所述地表温度构建特征空间，确定温度-植被干旱指数；根据田间试验，确定土壤含水量临界最大、最小值等关键土壤水分参数；计算灌溉地和非灌溉地的田间表层土壤含水量。可选的，所述利用遥感信息对所述潜在灌溉面积内作物生育期的田间表层土壤含水量进行反演计算，还包括步骤：根据田间实测土壤含水量对所述表层遥感土壤含水量进行修正。可选的，所述根据所述作物生育期的田间表层土壤含水量从所述潜在灌溉面积中进行提取，得到所述实际灌溉面积，包括：结合所述灌溉制度信息，确定各次灌溉时间范围和遥感影像，并通过所述作物生育期的田间表层土壤含水量的变化从所述潜在灌溉面积中提取所述实际灌溉面积；或者，对于干旱少雨地区，结合所述灌溉制度信息，通过所述作物生育期的田间表层土壤含水量逐次提取所述实际灌溉面积。或者，对比各次灌溉时段范围内灌溉地和非灌溉地的表层遥感土壤含水量差值，逐次从所述潜在灌溉面积中提取所述实际灌溉面积。还包括：根据当年降水量、作物类型及其耗水量和所述灌溉制度信息提取所述实际灌溉面积。从上述技术方案可以看出，本发明提供了一种基于多源信息的实际灌溉面积提取方法，具体为根据预先获取的多源信息提取潜在灌溉面积，所述多源信息包括基础地理信息、气象信息、灌排工程信息、农业种植结构信息和灌溉制度信息；利用遥感信息对所述潜在灌溉面积内作物生育期的田间表层土壤含水量进行反演计算；根据所述田间表层土壤含水量采用不同方法从所述潜在灌溉面积中进行提取，得到实际灌溉面积。通过上述操作，可以快速、有效获得实际灌溉面积，从而能够为农业用水管理提供科学、客观的依据。附图说明为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本发明实施例提供的一种基于多源信息的实际灌溉面积提取方法的步骤流程图；图2为本发明实施例提供的一种潜在灌溉面积的提取方法的步骤流程图；图3为本发明实施例提供的一种田间表层土含水量的反演计算过程的步骤流程图；图4为本发明实施例提供的一种实际灌溉面积的提取方法的步骤流程图。具体实施方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。图1为本申请实施例提供的一种基于多源信息的实际灌溉面积提取方法的步骤流程图。如图1所示，本实施例提供的实际灌溉面积的提取方法包括如下步骤：s1、根据预先获取的多源信息提取潜在灌溉面积。潜在灌溉面积具体来说是通过一般常识性的经验或规律判定某一地块是否灌溉，即是否具备灌溉条件。潜在灌溉面积的确定主要是为实际灌溉面积提取设定一个合理、可能的范围。潜在灌溉面积识别提取如图2所示。具体的提取过程如下描述：获取多源信息。多源信息包括基础地理信息、气象信息、灌排工程信息、土地利用信息和灌溉制度信息。具体如下：(1)准备行政区划边界和灌区边界，以及灌区渠系和水源分布图。灌区边界(矢量数据)：①第一次全国水利普查—灌区专项普查数据(2011年)；②灌区管理部门提供的大比例尺灌区边界。行政区划边界(矢量数据)：①全国、省、市、县行政区划数据；②1:10万全国乡镇行政区划数据。灌区渠系和水源分布(矢量数据)：①第一次全国水利普查—灌区专项普查数据(2011年)；②灌区管理部门提供的大比例尺水利工程分布图。(2)准备土地利用图。土地利用(30×30m栅格数据)：当年全国土地利用数据，具体分类详见表1。表1其中，耕地分布(30×30m栅格数据)获取：①第二次全国土地调查数据；②美国发布的2015年全球30m农田分布图。(3)准备全国地貌图。根据全国地貌图，提取山区、丘陵和平原分布，并与耕地分布进行叠加，获取不同地形条件下的耕地分布。地貌具体分类详见表2。编码地貌分类编码地貌分类编码地貌分类编码地貌分类1山地2丘陵3台地4平原11极大起伏山地21梁峁丘陵31高台地41起伏平原12大起伏山地22高丘陵32中台地42倾斜平原13中起伏山地23中丘陵33低台地43平坦平原14小起伏山地24低丘陵15喀斯特山地25喀斯特丘陵26沙丘表2地形分布(栅格数据)：①2004年全国90×90m栅格地貌数据；②根据30mdem(数字高程模型)提取地貌分布；③大比例尺县域地形图。(4)准备土壤类型分布图。土壤类型分布(矢量数据)：1：100万中华人民共和国土壤图。(5)生成全国降水量分布图。基于国家基础气象站网逐日站点观测数据，自动生成全国多年平均降水量等值线图、当年降水量等值线图和季节降水量等值线图。2、提取农业种植结构信息。利用种植结构遥感监测分析软件获取农业种植结构信息。主要提取水稻、小麦、玉米、棉花、大豆、土豆，以及区域种植面积较大的作物、规模设施农业，以及果树和经济林分布，并对其类别属性进行编码。种植结构(30×30m栅格数据)：利用作物种植结构遥感监测分析软件，通过遥感解译方式获得。3、建立主要作物灌溉制度数据库，对所述灌溉制度信息进行数据分析。基于全国灌溉试验站网灌溉试验数据，建立主要作物灌溉制度数据库，包括作物生育阶段、旬作物需水量、旬降水量、不同水平年作物需水量、多年平均降水量、当年降水量、种植制度(如水稻，需要进一步划分早、中、晚稻)、灌溉制度等。灌溉试验数据分析：主要分析旬作物需水量、旬降水量的关系。需要说明的是不同灌溉试验站点对应的代表区域，此区域全部以该站点试验资料为基础进行判断，代表区域划分具体可以参考各省(区、市)灌溉用水定额编制分区。4、提取潜在灌溉面积。(1)有灌区边界情况①渠灌区、井渠结合灌区。理论上，灌区内灌排工程控制的耕地面积即为潜在灌溉面积(包括设施农业)；无详细资料时，可将灌区范围内的耕地面积作为潜在灌溉面积，非耕地中成规模的果树、设施农业等即作为潜在灌溉面积。同时，可利用灌区渠系分布图、灌溉机井分布图、水源特征数据(位置、形式、供水能力)等，对其控制边界范围进行辅助识别与复核。最后，生成灌区内的耕地潜在灌溉面积、非耕地潜在灌溉面积。对于渠灌区，通过将渠系分布图和耕地、非耕地(果树、设施农业等)分布叠加，直观判断渠系末端是否在耕地、非耕地范围内，可将包括渠系在内的耕地、非耕地范围作为渠灌区控制边界，若渠系末端超出耕地、非耕地范围，则需要进一步现场核实。对于井渠结合灌区，根据灌区外围机井位置、单井出水量，以及灌溉作物类型和综合灌溉平均灌溉定额估算单井控制的灌溉面积；同时，结合上述渠灌区控制边界判断方法，可将包括机井和渠系在内的耕地、非耕地范围作为井渠结合灌区控制边界，若外围机井分布、渠系末端超出耕地范围，则需要进一步现场核实。单井控制的灌溉面积，其计算公式如下：式中，f0为单井控制灌溉面积(hm2)；q为单井出水量(m3/h)；t为灌溉期间开机时间(h/d)；t为每次轮灌期的天数(d)；η为灌溉水利用系数；η1为干扰抽水的水量消减系数；m为综合平均灌溉定额(m3/hm2)。综合平均灌溉定额，其计算公式如下：式中，ai为第i种作物的种植面积(hm2)；mi为第i种作物的灌溉用水定额(m3/hm2)；i为不同作物编号，i＝1,2,…,n。②集中连片有边界的井灌区或其他灌区。灌区范围内耕地、设施农业、林果地等面积，即为潜在灌溉面积。同时，可利用灌溉机井分布图，对其控制边界范围进行辅助识别与复核。对于井灌区，通过将机井分布图和耕地、非耕地分布叠加，根据灌区外围机井位置、单井出水量，以及灌溉作物类型和综合平均灌溉定额估算单井控制的灌溉面积。可将包括机井在内的耕地、非耕地范围作为井灌区控制边界，若外围机井分布超出耕地、非耕地范围，则需要进一步现场核实。(2)无灌区边界情况①结合种植结构和灌溉制度，按地貌分区域提取潜在灌溉面积，提取规则具体详见表3。表3就全国范围而言，需要通过对各省(区、市)进行典型调查，在详细了解当地灌溉习惯基础上，针对不同区域进一步制定详细的判别规则。这里需要特别注意南方丘陵山区灌溉的山间盆地和西北黄土高原灌溉的谷地。结合表4按地貌分区域提取潜在灌溉面积规则，可进一步根据主要作物生育期耗水量(作物耗水量由遥感反演获取)et与当年作物生育期实际降水量p进行比较判断。具体操作中可分区域细化确定主要作物类型及种植情况，结合灌溉试验、模型模拟、文献资料等综合分析确定不同区域主要作物生育期降水临界值α(其值确定还需要根据实际情况进一步研究确定)，若et-p〉α，则判定是灌溉，否则判定为非灌溉。某一分区考虑地貌、作物、耗水和降水因素的判别规则表详见表4。表4②条件具备时，收集和调查灌区水源分布、水源特征数据(类型、供水能力等)，水源类型可按河流(引水和提水)、塘坝、集雨水窖、灌溉机井等，构建灌区水源分布图。尤其要对区域规模水窖灌溉区域、规模纯井灌溉区域等进行重点调查，实际操作中可由灌区管理单位协助。尤其要对丘陵、山区耕地范围内的雨养农业进行重点排查，水源分布可作为辅助识别，确认为雨养农业的予以扣除。(3)统一数据格式标准。上述栅格数据、矢量数据要统一坐标系/投影(2000国家大地坐标系，投影转换注意分区)，栅格数据均为30×30m分辨率。5、提取灌区或行政区的潜在灌溉面积数据。利用行政区划边界或灌区边界，提取形成灌区，以及全国、省、市、县行政区的30×30m潜在灌溉面积栅格数据产品。另外，在得到潜在灌溉面积后还对其进行对比分析，具体步骤如下：1、与统计资料对比分析(1)收集最新各省、市、县行政区有效灌溉面积、土地利用、种植结构等统计资料。(2)将全国、省、市、县行政区有效灌溉面积、耕地面积、种植结构与各地统计数据进行对比分析，绘制相关关系散点图。(3)确定潜在灌溉面积差异较大区域，分析原因，进行进一步复核分析，提出需现场或与有关部门沟通的重点区域和解决对策。(4)分析确定不同行政区域潜在灌溉面积，制作形成30×30分辨率的栅格数据和空间分布图。2、与自然资源部门土地利用资料对比分析按照不同行政区域，分析土地利用中水田和水浇地面积，两者之和与潜在灌溉面积进行对比。3、通过资料对比分析最终确定潜在灌溉面积通过上述综合分析，形成全国、省、市、县的潜在灌溉面积分布图、数据表、柱状图，以及逐年变化图等。4、灌溉地和非灌溉地选取在潜在灌溉面积图上，对灌溉地和非灌溉地进行选取，具体步骤如下：①基于全国多年平均降水量等值线图，按照年降水量200mm以下、200～400mm、400～800mm、800mm以上分档标准进行多年平均年降水量分区。具体操作中也可根据区域范围大小在上述分档标准基础上适当进行加密处理，以反映区域降水的空间变异性。②结合地貌图、土壤质地分布图，通过上述不同图层的叠加，识别确定降水、地形和土壤质地属性一致(或相似)的区域。③结合耕地分布图和土地利用图，在googleearth图上首先根据实际需要在特定区域内均匀确定耕地范围内一定灌溉的灌溉地，提取灌溉地的点坐标并命名编号；再按照对照点距离最小、作物相似原则，确定非灌溉地，提取非灌溉地的点坐标并进行编号，为实际灌溉面积提取提供基础支撑。s2、利用遥感信息对作物生育期的田间表层土壤含水量进行反演计算。具体计算过程如图3所示。1、对多源遥感影像进行预处理，包括几何校正、大气校正、镶嵌与裁切等。2、利用预处理的遥感影像提取归一化差值植被指数(normalizeddifferencevegetationindex,ndvi)，其计算公式如下：式中，nir代表近红外波段处的反射率值；r代表红波段处的反射率值，ndvi为归一化差值植被指数。3、根据预处理过的遥感影像和大气剖面参数计算地表温度(landsurfacetemperature,lst)(1)准备大气剖面参数：大气在热红外波段的透过率τ、大气向上辐射亮度l↑和大气向下辐射亮度l↓。(2)计算植被覆盖度，其计算公式如下：pv＝[(ndvi-ndvisoil)/(ndviveg-ndvisoil)]式中，ndvi为归一化差值植被指数；ndvisoil为完全裸土或无植被覆盖区域像元的ndvi值；ndviveg为完全被植被覆盖区域像元的ndvi值。一般情况下，ndviveg和ndvisoil取经验值，ndviveg＝0.70，ndvisoil＝0.05，即当某个像元的ndvi≥0.70时，pv取值为1；当ndvi＜0.05时，pv取值为0。(3)计算地表比辐射率，其计算公式如下：ε＝0.004*pv+0.986式中，pv为植被覆盖度，ε为地表比辐射率。(4)计算黑体热辐射亮度b(ts)，其计算公式如下：b(ts)＝[lλ-l↑-τ(1-ε)l↓]/τε式中，lλ为卫星传感器接收到的热红外辐射亮度值；l↑为大气向上辐射亮度；l↓为大气向下辐射亮度值；ε为地表比辐射率；τ大气在热红外波段的透过率，b(ts)为温度为t的黑体在热红外波段的辐射亮度值。(5)计算地表真实温度ts，其计算公式如下：ts＝k2/ln(k1/b(ts)+1)式中，k1和k2均为常量(对于landsat5thematicmapper(tm)，k1＝607.76w/(m2*μm*sr)，k2＝1260.56k；对于landsat7enhancedthematicmapper(etm+)，k1＝666.09w/(m2*μm*sr)，k2＝1282.71k；对于landsat8thermalinfraredsensor(tirs)传感器，band10，k1＝774.89w/(m2*μm*sr)，k2＝1321.08k)；ts为地表温度(k)。4、根据归一化差值植被指数ndvi和地表温度lst构建特征空间，确定温度-植被干旱指数(temperaturevegetationdrynessindex,tvdi)。tsmin＝a1+b1ndvitsmax＝a2+b2ndvi式中，ts代表任意像元的地表温度(k)；tsmin表示某一ndvi对应的最低温度，对应湿边；tsmax表示某一ndvi对应的最高温度，对应干边；a1、b1为湿边方程的拟合系数；a2、b2为干边方程的拟合系数。tvdi在[0,1]之间，当(ndvi，lst)越接近于干边时，下垫面土壤越干燥，在干边上tvdi＝1；当(ndvi，lst)越接近湿边时，下垫面越湿润，在湿边上tvdi＝0。5、根据实测土壤含水量来确定土壤含水量最大值ssmmax、土壤含水量最小值ssmmin，最终确定任意点土壤表层含水量ssm。式中，ssmmax为土壤含水量最大值，即田间持水率和饱和含水率均值；ssmmin为土壤含水量最小值，即永久凋萎点。6、灌溉地和非灌溉地典型田块土壤含水量及土壤水分特性参数确定。对于大范围区域，通过土壤类型分布图，确定土壤质地类型和相关参数(凋萎含水率、田间持水率和饱和含水率等土壤水分特性参数)。对于典型田块，按照以下步骤确定：①考虑灌区土壤质地和主要作物类型分布，进行合理分区；②大田采样地块规格为60×60m，尽可能不同土壤质地、不同作物类型组合设置一处典型地块；③为了能够区分降水和灌溉引起的土壤水分变化，需在灌溉地典型田块附近，选择非灌溉地典型田块作为参照，非灌溉地选取要求是气候条件、作物类型、土壤质地、耕作制度等方面与灌溉地基本一致。所需测定参数：土壤含水量、凋萎含水率、田间持水率和饱和含水率等。监测深度划分：0～5cm、5～10cm、10～20cm土壤含水量测定方法：烘干法监测频率：根据不同作物灌溉日期，灌溉前加测一次，灌溉后每日监测一次，监测时间截止为灌溉后第10天。采样点的位置要具有代表性，能够代表田块的土壤湿度平均水平，且每次采样点的位置尽可能固定。7、获取灌溉地典型田块降水量。对于大范围区域，通过国家基础气象站点降水要素空间插值或遥感降水产品，获取降水量。对于典型田块，有条件的可采用自计雨量计，或者采用简易的雨量筒进行人工采集。记录逐日的降水时间(起始时刻和终止时刻，时间范围为08:00-20:00)、降水量(mm)。8、对遥感田间表层土壤含水量进行修正。选择灌概地和非灌溉地等代表性点对遥感土壤含水量进行检验修正，消除算法等引起的系统误差。θo＝k·θrs+ｂ式中，θo为观测土壤含水量(cm3/cm3)，θrs为遥感土壤含水量(cm3/cm3)，k和b分别为标定系数，根据不同作物、不同生育期，由观测资料确定，并建立历史资料数据库。s3、从潜在灌溉面积范围内提取实际灌溉面积。在潜在灌溉面积确定的基础上，进一步提取灌区的实际灌溉面积。结合获取的遥感影像、种植结构、灌溉制度、降水等多源信息，采取以下步骤进行实际灌溉面积提取，具体如图4所示：1、在潜在灌溉面积范围内，根据当年降水情况、种植作物类型及其耗水量和灌溉制度信息，判断发生灌溉的作物和面积(设施农业、蔬菜、玉米、小麦和水稻等)。在不同作物生育期内(旬尺度、月尺度)，若作物耗水量大于当年降水量，则确定为发生灌溉，即：作物耗水量-当年降水量＞临界阈值时，确定灌溉。临界阈值通过试验资料确定，其主要影响因素包括潜水补给(受地下水埋深约束)、作物类型、土壤类型等。根据不同区域划分给定不同的判别阈值。2、结合不同作物灌溉制度信息，确定各次灌水时间范围和过境影像，通过遥感反演作物生育期的田间表层土壤含水量逐次提取发生灌溉的实际灌溉面积。假设各次灌溉时段范围内无灌溉发生，根据时段内降水资料，以时段初遥感土壤含水量值θt为初始值，根据土壤含水量衰减函数，推求时段末土壤含水量值θt+n：θｔ+ｎ＝ｃ·θt·ｔ-ｍ若：θt+n'-θt+n>d阈值则时段内发生灌溉。式中，θt+n为时段末估算的土壤含水量(cm3/cm3)，θt+n'为时段末遥感土壤含水量(cm3/cm3)，θt为时段初遥感土壤含水量(cm3/cm3)；c为常数；m为衰减系数；t为时间(d)；d阈值根据不同作物、不同生育期，通过试验资料确定。3、对于干旱、少雨地区，结合不同作物灌溉制度信息和过境影像，通过遥感反演作物生育期的田间表层土壤含水量的变化逐次提取发生灌溉的面积。假设时段内无降水发生，以时段初遥感土壤含水量值θt为初始值，时段末遥感土壤含水量值为θt+n：若：θt+n-θt>d阈值则时段内发生灌溉。式中，θt+n为时段末土壤含水量(cm3/cm3)，θt为时段初土壤含水量(cm3/cm3)。d阈值根据不同作物、不同生育期来确定，通过试验资料确定。4、对比各次灌溉时段范围内灌溉地和非灌溉地遥感土壤含水量差值变化，逐次提取发生灌溉的面积。若：(θi-θni)t+n-(θi-θni)t>d阈值则时段内发生灌溉。式中，θi为灌概地土壤含水量(cm3/cm3)，θni为非灌概地土壤含水量(cm3/cm3)；t为时段初日期，t+n为时段末日期，n为时段间隔日期(d)。当公式成立时，则认为有灌溉发生，统计像元及灌溉面积。d阈值根据不同作物类型、不同生育期，通过试验资料确定。5、根据不同研究地区和数据资料情况选择不同实际灌溉面积提取方法，统计目标时段内(年、月、作物生育期)发生灌溉的面积，并对不同方法监测结果进行叠加分析，确定当年的实际灌溉面积。6、对比省、市、县行政区统计资料的实际灌溉面积，分析两者差异性，并通过人工地面调查，统计遥感反演实际灌溉面积精度，如符合度指数、百分比等指标。从上述技术方案可以看出，本实施例提供了一种基于多源信息的实际灌溉面积提取方法，具体为根据预先获取的多源信息提取潜在灌溉面积，所述多源信息包括基础地理信息、气象信息、灌排工程信息、农业种植结构信息和灌溉制度信息；利用遥感信息对所述潜在灌溉面积内作物生育期的田间表层土壤含水量进行反演计算；根据所述田间表层土壤含水量采用不同方法从所述潜在灌溉面积中进行提取，得到实际灌溉面积。通过上述操作，可以快速、有效获得实际灌溉面积，从而能够为农业用水管理提供科学、客观的依据。本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。本领域内的技术人员应明白，本发明实施例的实施例可提供为方法、装置、或计算机程序产品。因此，本发明实施例可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。本发明实施例是参照根据本发明实施例的方法、终端设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理终端设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理终端设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理终端设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理终端设备上，使得在计算机或其他可编程终端设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程终端设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。尽管已描述了本发明实施例的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例做出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明实施例范围的所有变更和修改。最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者终端设备中还存在另外的相同要素。以上对本发明所提供的技术方案进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。当前第1页12