基于GIS的农业面源污染数据管理与查询系统设计方案.docx

1、基于GIS的农业面源污染数据管理与查询系统设计方案基于GIS的农业面源污染数据管理与查询系统设计方案1方案背景1.1方案意义近年来，农业面源污染已经成为中国流域性水体污染的主要根源和空气污染的重要来源。在我国水体污染严重流域，由农田、农村畜禽养殖地带和城乡结合部的生活排污而造成流域水体氮、磷素营养化已大大超过了来自城市地区的生活点源污染和工业点源污染。随着农业面源污染的日益严重和人们对生态环境的关注，农业面源污染问题引起了国内外广大学者和政府决策部门的高度关注。大多方案者对农业面源污染问题的方案集中在面源污染概念和内涵的界定、面源污染类型、形成机理和治理模式4个方面，但始终未能对农业面源污染数

2、据进行具体分析，无法找出引起农业面源污染的主控因素，使得农业面源污染的治理效果不够明显。本系统通过对农业面源污染中的种植业源污染物中的肥料流失系数进行测算，并将GIS技术运用于此，将结果显示在中国地图上，来对农业面源污染数据进行管理，并设计多种查询方式，方便用户进行查询，以期为进一步根据全国污染源调查数据测算全国农业面源污染负荷奠定科学基础。本方案将对解决农业面源污染问题具有重大的意义。1.2方案目的农业面源污染在中国日趋严重，控制农业面源污染迫在眉睫。要想有效的治理农业面源污染的问题，必须找出导致农业面源污染的主控因素，这就需要通过实验调查得到试验数据，进而分析数据，找出导致农业面源污染的因

3、素。本系统就是针对农业面源污染中的种植业源污染物的数据进行分析，从而得出分析结果，并可以将测算结果显示给用户。具体包括四个方面：（1）通过对农业面源中的种植业源污染物的数据进行分析，掌握种植业污染特点和规律，改变不合理的生产方式，减少农业生产过程中的环境污染。（2）摸清污染底数，为政府制定农业环境保护相关政策提供决策依据。（3）为用户提供一个平台，并设计多种查询方式，方便用户查询测算结果。（4）通过一个良好的用户界面，来实现对农业面源污染数据的管理。1.3设计目标（1）计算我国不同区域、土壤类型、坡度、气候和典型种植制度条件下氮、磷等主要污染物的流失系数，为准确测算种植业源污染物流失数量提供基

4、本参数。（2）摸清我国种植业源污染物流失的时空变化规律和主控因素，为控制种植业源污染提供科学指导。（3）实现对农业面源污染数据的管理。（4）将监测点按照实际的经纬度显示在中国地图上对应的位置，并用两种不同的颜色和形状来表示不同的监测点类别，红色表示重点监测点，绿色表示一般监测点。（5）通过多种查询方式，对农业面源污染数据的分析结果进行查询。1.4农业面源污染监测概述1.4.1监测分区在综合考虑肥料、农药和农膜污染的发生规律和主要影响因素，如地形、气候、土壤、作物种类与布局、种植制度、耕作方式、灌排方式等的基础上，本次方案主要依据地形和气候特征，将全国种植业污染源划分为六大区域。其中，4个区域以

5、平原为主，山地丘陵比例较少；2个区以山地、丘陵为主，平原所占比例较低。为便于操作，表1-1列出了各省监测点的分布情况。1.4.2监测布局全国共设置地下淋溶和地表径流一般定位监测试验372个。其中，地下淋溶试验点140个，包括大田试验点47个，保护地菜田40个，露地菜田31个，果园22个；地表径流试验点232个，包括水田试验点46个，水旱轮作51个，旱地平原57个，坡耕地78个，详见表1-1。除一般监测试验点以外，全国还设置35个地下淋溶和地表径流重点监测试验点，详见表1-2。其中，地下淋溶试验点14个，包括大田9个，保护地菜田3个，露地菜田2个；地表径流监测试验点21个，包括水田3个，水旱轮作

6、4个，旱地平原5个，坡耕地9个。表1-1 种植业源污染物流失系数测算一般监测试验点分布表实施地点1、地下淋溶监测点数量（个）2、地表径流监测点数量（个）大田保护地菜田露地菜田果园小计水田水旱轮作旱地平原坡耕地小计坡度 5-15坡度 15北京12216112天津12216112河北4322111113河南4322110内蒙古213112山西2215112山东44431511黑龙江31422116吉林211411114辽宁21115112陕西22127112甘肃51219448宁夏1111411新疆42390新疆兵团440青海213112上海01236浙江112331110江苏1113144111

7、1安徽222172642115湖北12144863223湖南2136762223江西03342113海南0312118福建02431212四川32381633316贵州0433313云南2132223312重庆123123129广西02222210广东0448大连11112青岛1111深圳011南京0112合计474031221404651574236232表1-2种植业源污染物流失系数测算重点监测试验点分布表监测类型序号地点农田类型备注地表径流1福建水田双季稻2湖南水田双季稻3吉林水田单季稻4安徽水旱轮作稻-麦/油5湖北水旱轮作稻-麦/油6江苏水旱轮作稻-麦/油7四川水旱轮作稻-麦/油8广西

8、旱地平原甘蔗9湖北旱地平原棉花10湖南旱地平原露地蔬菜11江苏旱地平原露地蔬菜12云南旱地平原露地蔬菜13四川坡耕地旱坡地1514福建坡耕地旱坡地，梯田15甘肃坡耕地旱坡地1516贵州坡耕地旱坡地1517重庆坡耕地旱坡地5-1518云南坡耕地旱坡地5-1519河北坡耕地旱坡地5-1520江西坡耕地旱坡地5-1521西藏坡耕地旱坡地5-15地下淋溶22北京保护地菜田日光温室23山东保护地菜田日光温室24宁夏保护地菜田日光温室25北京露地菜田26山东露地菜田27黑龙江大田大豆28吉林大田玉米29黑龙江大田玉米30甘肃大田小麦|玉米31河北大田小麦-玉米32陕西大田小麦-玉米33北京大田小麦-玉米3

9、4新疆大田棉花35重庆大田水稻-油菜1.4.3监测方法本系统进行农业面源数据分析的数据，是经过污染源调查的试验而得到的数据结果。以下是对污染源调查试验的简单概述。（1）选点原则代表性。监测地块的土壤类型、地块坡度、种植制度、耕作方式、栽培模式、灌排方式等在各大分区中要有一定代表性，地块土壤肥力和作物产量水平应能够代表所在区域的正常水平。可操作性。监测点的选择要兼顾交通、工程建设和监测设施维护。安全性。为确保不受到人为破坏、土地征用或土地使用纠纷，重点监测试验点最好设在有独立产权的试验站或农场。（2）试验处理肥料流失系数测算试验可分为一般监测试验和重点监测试验；按方案内容不同，肥料流失系数测算试

10、验可分为地下淋溶监测试验和地表径流监测试验。各类监测试验的处理分述如下：一般监测试验一般监测试验共设试验点372个，全国各个省、市、区基本上均有分布（表1-1），监测周期为1年。地下淋溶监测试验和地表径流监测试验均设置两个处理，分别为：处理1 对照处理，不施任何肥料。处理2 常规处理，肥料施用量、施用方法和施用时期完全遵照当地农民生产习惯。每个处理设3次重复，每个监测点共计6个小区。小区面积、形状、规格完全相同，小区面积不小于20m2。建议平原区小区规格为4m 68m，山地丘陵区小区规格为34m 810m。两个处理除施肥不同以外，其它田间管理措施完全相同。重点监测试验重点监测试验共设监测点35

11、个，分布在全国23个省、市、区。主要目的在于尽可能减少由于自然和人为偶然因素所造成的误差，提高肥料和农药流失系数的测算精度。地下淋溶和地表径流重点监测试验点设置56个处理，其中对照和常规两个处理与一般监测试验完全相同；其余4个处理的内容围绕施肥、灌溉、耕作、秸秆覆盖、坡地种植方向等对肥料流失系数影响较为显著的农艺管理措施设置，目的在于探明种植业源污染物流失规律和主控因素。每个试验点的处理内容可根据当地农业生产实际情况确定，但同一种植制度、各个重点监测试验点之间的处理内容应保持一致。重点监测试验设3次重复，小区规格、排列方式、径流池或淋溶盘安装位置均可参照一般监测试验。（3）监测点建设依据肥料污

12、染的迁移途径，监测点建设可分为两种类型：地表径流监测试验地表径流是南方湿润平原区、北方高原山地区和南方山地丘陵区肥料流失的主要途径。地表径流的监测方法主要考虑便于监测、计算径流水量、泥沙流失量以及氮、磷流失量。在综合比较各种方法优缺点的基础上，本项目采用径流池法监测肥料通过地表径流流失的数量。地下淋溶监测试验地下淋溶是东北半湿润平原区、西北干旱半干旱平原区和黄淮海半湿润平原区肥料流失的主要途径。在设计监测方法时主要考虑便于采集、计算淋溶水量。在综合比较各种方法优缺点的基础上，本项目采用田间原位渗漏计方法监测地下淋溶。（4）监测周期地下淋溶和地表径流试验监测周期均为一年，不仅包括作物生重点，也包

13、括休闲期。监测起始时间一般为某种作物播种日，结束时间为翌年同一天。（5）肥料流失系数测算方法以地表径流（或地下淋溶）途径流失的氮、磷数量等于整个监测周期中（一个完整的周年）各次径流水中污染物浓度与径流水（或淋溶水）体积乘积之和。计算公式如下：P=其中：P为污染物流失量（mg）Ci为第i次径流（或淋溶）水中氮、磷的浓度；Vi为第i次径流（或淋溶）水的体积（cm）肥料流失系数的表示方法，以流失率（%）表示，以氮素为例，计算公式如下：肥料氮素流失率（%）=1.5技术路线系统开始于目标和任务的确定，随后进行数据图件统计资料的收集入库，通过数据库统一管理空间数据和属性数据，然后在这些数据的基础上实现各个

14、模块的功能，最后进行整个系统的集成。工作流程图如图1-1.图1-1 工作流程图1.6开发平台和开发方式的选择1.6.1开发环境和开发方式的确定（1）在编程语言方面，选用Java语言。（2）在数据库软件方面，选用 SQL Sever 2000作为数据库建立平台。（3）在GIS软件方面，选用Arcgis Engine进行开发。（4）开发方式选择集成二次开发。1.6.2相关工具和开发方式的介绍（1）Java语言Java是一种可以编写跨平台应用软件的面向对象的程序设计语言1，提供了丰富的数据库访问技术，可以方便地与一些大型数据库软件如SQL Sever 2000、Oracle等连接，Java应用具有无

15、比的健壮性和可靠性。（2）SQL Server 2000数据库SQL Sever 2000具有实用性、易扩展性和安全性等优点，且在Windows系列操作系统上具有最佳的表现和整合性。（3）Arcgis EngineArcGIS Engine是用于构建定制应用的一个完整的嵌入式的GIS组件库，这些GIS组件可以在Arcgis Desktop 框架之外使用。利用ArcGIS Engine，开发者能将ArcGIS功能集成到一些应用软件，还可以为用户提供针对GIS解决方案的定制应用。ArcGIS Engine可以在Windows，UNIX，和Linux桌面上运行并支持一系列的应用软件开发环境。除了支持

16、COM环境之外，ArcGIS Engine还支持C+，.NET(C#)，Java和Delphi，使开发者能够跨操作系统、选择多种开发构架，通过ArcGIS Engine进行开发2。 （4）开发方式介绍应用型地理信息系统主要由以下三种方式来实现：独立开发方式，不依赖于 GIS 工具软件，从空间数据的采集、编辑到数据的处理分析及结果的输出，所有的算法都由开发者选择某种程序设计语言：如 Visual Basic，Visual C+ +，Delphi，Java等，在一定的操作系统平台上自我编程实现。单纯二次开发，完全借助于GIS工具软件提供的开发语言进行应用系统开发，如ESRI的Arcview提供的A

17、venue 语言，Mapinfo公司研制的Mapinfo Profes2sional 提供的MapBasic语言等，用户利用这些宏语言，以原GIS工具软件为开发平台，开发出自己的针对不用应用对象的应用程序。集成二次开发，集成二次开发是指利用专业GIS工具软件，如Arc2view，Mapinfo，Supermap Objects等，实现GIS的基本功能，以通用可视化开发工具为开发平台，进行二者的集成开发。GIS的独立开发难度大，单纯二次开发受GIS工具提供的编程语言限制而差强人意，而结合组件式 GIS工具软件与可视化语言的集成二次开发方式，既可以充分利用GIS软件对空间数据库的管理、分析功能，又

18、可以利用可视化语言高效便利的编程优点，集二者之所长 ,在提高运用系统开发效率的同时具有更好的外观效果、更强大的数据库功能，而且可靠性好、易于移植、便于维护3。2数据库设计数据库是系统的核心部分，数据库（Data Base）是以一定的组织方式存储在一起的相互关联的数据集合，数据库组织数据的好坏，直接影响到以后读取数据库和系统运行的效率。因此，我们在使用数据库之前，必须对数据库进行设计。数据库设计(Database Design)是指根据用户的需求，在某一具体的数据库管理系统上，设计数据库的结构和建立数据库的过程。数据库设计的目的是为了最大限度地减少数据冗余，增强数据间的联系，实现对数据的合理组织

19、和灵活存取。2.1基础数据库设计在本系统中，采用SQL Sever 2000作为基础数据库。用到的基本表有两个。分别是监测点基本信息表（如表2-1）和肥料流失系数测算结果表（如表2-2）。 在数据库中对以上两个基本信息表建立了视图，通过视图来实现数据管理，实现对数据的添加、删除和修改等操作。视图是从基本表里导出的表。它与基本表的不同之处是它是一个虚表。在数据库中只存放视图的定义，不存放视图对应的数据，对视图的操作，最终都要转换为对基本表的操作。在本论文中，将监测点基本信息和肥料流失系数的测算结果建立视图，是为了得到一个监测点肥料流失情况的详细信息，方便管理员对数据进行管理。视图主要包括的字段：

20、地块代码，处理号，所属分区，地形，梯田否，种植方向，利用类型，种植模式，监测点类别，监测类别，经度，纬度，施用化肥N含量（公斤/亩），施用化肥P含量（公斤/亩），施用有机肥N含量（公斤/亩），施用有机肥P含量（公斤/亩），施肥含N总量（公斤/亩），施肥含P总量（公斤/亩），N流失总量（公斤/亩），P流失总量（公斤/亩），N相对流失总量（公斤/亩），P相对流失总量（公斤/亩），N流失系数（%）和P流失系数（%）。表2-1 监测点基本信息表属性名数据类型长度含义landnochar10地块代码areachar20所属分区terrainchar10地形ynchar10梯田否plantdirchar1

21、0种植方向usetypechar10利用类型plantmodulechar20种植模式dottypechar10监测点类别monitortypechar10监测类别longitudechar20经度latitudechar20纬度stexturechar10土壤质地landformchar10地貌类型fertilizerchar10肥力水平soilnamechar20土壤名称表2-2 肥料流失测算结果表属性名数据类型长度含义landnochar10地块代码dealnoint 4处理号usefnfloat8施用化肥N含量（公斤/亩）usefpfloat8施用化肥P含量（公斤/亩）useorgnf

22、loat8施用有机肥N含量（公斤/亩）useorgpfloat8施用有机肥P含量（公斤/亩）usesumnfloat8施肥含N总量（公斤/亩）usesumpfloat8施肥含P总量（公斤/亩）flowsumnfloat8N流失总量（公斤/亩）flowsumpfloat8P流失总量（公斤/亩）reflownfloat8N相对流失总量（公斤/亩）reflowpfloat8P相对流失总量（公斤/亩）ncoefficientfloat8N流失系数（%）pcoefficientfloat8P流失系数（%） 2.2 GIS数据库设计GIS数据库是某区域内关于一定地理要素特征的数据集合，主要涉及对图形和属性

23、数据的管理和组织。它与一般数据库相比，具有以下特点：（1）GIS数据库不仅有与一般数据库数据性质相似的地理要素的属性数据，还有大量的空间数据，即描述地理要素空间分布位置的数据，并且这两种数据之间具有不可分割的联系。（2）地理系统是一个复杂的巨系统，要用数据来描述各种地理要素，尤其是地理的空间位置数据量往往大得惊人，即使是一个极小区域的数据库也是如此。（3）数据的应用相当广，如地理方案、环境保护、土地利用与规划、资源开发、生态环境、市政管理、道路建设等等。GIS数据库中存放了大量的空间数据。空间数据是用来表示空间实体的位置、形状、大小及其分布特征诸多方面信息的数据，它可以用来描述来自现实世界的目

24、标，它具有定位、定性、时间和空间关系等特性。定位是指在已知的坐标系里空间目标都具有唯一的空间位置；定性是指有关空间目标的自然属性，它伴随着目标的地理位置 ；时间是指空间目标是随时间的变化而变化；空间关系通常一般用拓扑关系表示。空间数据是一种用点、线、面以及实体等基本空间数据结构来表示人们赖以生存的自然世界的数据。在本系统中，将空间数据存储于shape文件中。shape文件包括一个主文件，一个索引文件，和一个dBASE表，这三个文件是构成shape文件所必备的文件。三个文件具体用途如下：（1）主文件，即.shp格式的文件，存储的空间实体（2）索引文件，即.shx格式的文件，存储的是shp文件中空

25、间实体的索引。（3）dBASE表，即.dbf格式的文件，是属性数据格式，存储的是每一个空间实体对应属性数据表。在本系统中，使用了ArcGIS Desktop中的ArcMap Windows桌面应用。用ArcMap最终实现了将两类监测点按不同的颜色和形状显示在中国地图上。首先通过ArcMap中的Tools中的Add XY Data的菜单项，实现了按照经纬度生成点要素。再通过右击图层，在弹出的Layer Properties对话框中，选择Symbology按钮，再通过左边显示栏里的Unique values对点要素进行单值渲染，从而才实现将一般监测点和重点监测点按不同过的颜色和形状显示。经过这些操

26、作以后，shape文件中就生成了四个文件来存储这些空间数据。其中的三个文件就是以上所介绍的构成shape文件所必备的文件，还有一个是.prj格式的文件，存储的是定义坐标系统的信息。3系统设计与实现3.1系统目标及设计原则本系统的目标是将农业面源数据科学存储，图文并茂地显示农业面源概况，方便快捷地查询相关信息，并将查询结果以对话框的形式输出。因此，在整个系统的设计过程中，遵循以下原则：（1）共享性原则。本系统将涉及到的农业面源的数据标准化、规范化，方便资源共享。（2）稳定性原则。系统开发遵循实际管理需要，要求系统稳定可靠的运行，且可视化强，人机交互界面友好，易于操作和维护。（3）安全性原则。在保证数据的共享性，独立性和完整性的基础上，对不同的用户赋予不同的管理权限，以保证数据的安全性。（4）可扩展性原则。由于采用了基于工业标准的全组件式GIS软件Arcgis Engine作为主控平台的二次开发平台，具有强大的二次开发能力，并且不需要掌握特殊的开发语言就可以利用常见的软件开发语言如java、VB、VC、Delphi、.NET等灵活地进行扩展开发，