程序实现

（一）实现思路

要实现土壤环境质量评价，首先要整理好评价数据，主要包括表层土壤有机污染物分析数据图层、表层土壤地球化学全量分析数据图层及表层土壤采样数据图层3个，统一格式和投影，并存放在数据库中；用户可根据权限自动加载数据，并设置参数（如区块颜色、评价指标、评价因子、模型库等）直接进行评价，评价的结果可用报表或色块图输出，关键在于模型库的建立，具体按如下流程建立。

（1）“一票否决法”评价模块

第一步，首先连接SQL Server数据库，加载3个数据图层，分别是表层土壤有机污染物分析数据图层、表层土壤地球化学全量分析数据图层及表层土壤采样数据图层。

第二步，遍历表层土壤地球化学全量分析数据点，读取pH、Cd、Hg、As、Cu、Pb、Cr、Zn、Ni字段对应的值，先判断该点的pH值区间范围，根据此pH 值范围判断Cd、Hg、Pb、Zn、Ni元素含量是否超标，如果超标则表示该评价单元不合格，否则从表层土壤采样数据图层中读取该点的4个采样点的属性字段（“土地利用”）的值，如果4个采样点的土地利用情况均为13（水浇地），则As、Cu、Cr 3个重金属元素使用旱地的标准，再做是否超标的判断，如果没有超标则转用下一步再做判断。

第三步，根据该点的位置读取表层土壤有机污染物分析数据图层中相应的土壤有机污染物分析测试数据，判断pH值区间范围确定DDT和六六六的标准值，最后得出该点是否超标。

（2）单元素污染指数评价模块

第一步，首先连接SQL Server数据库，加载3个数据图层，分别是表层土壤有机污染物分析数据图层、表层土壤地球化学全量分析数据图层及表层土壤采样数据图层。

第二步，获取要进行评价的评价因子（如Hg）。

第三步，遍历表层土壤地球化学全量分析数据点，分别读取pH、Hg字段对应的值，先判断某点的pH值区间范围，根据此pH值范围判断Hg元素的评价标准值（SHg），当然也可以区域背景平均值、平均值加2（3）倍标准差为评价标准，将该点的实测值与评价标准值（SHg）相除，即为该Hg元素的污染指数。

如果选取的是As、Cu、Cr 3个元素，则还要进一步根据表层土壤采样数据图层中的“土地利用”字段值来判断元素的评价标准值（S）。

如果选取的是DDT和六六六2个指标，则要遍历表层土壤有机污染物分析数据点，读取pH、DDT（或六六六）的字段值，根据其pH值区间范围确定DDT（或六六六）的标准值，同样该点的实测值与评价标准值（S）相除，即为该DDT（或六六六）的污染指数。

第四步，根据单元素污染指数值分等定级土壤质量，一般（默认）采用如下标准：

单因子污染指数＜1，为非污染区；

单因子污染指数介于1～2之间，为轻污染；

单因子污染指数介于2～3之间，为中度污染；

单因子污染指数大于3，为重度污染。

（3）内梅罗（Nemrow）综合污染指数评价模块

第一步，首先连接SQL Server数据库，加载3个数据图层，分别是表层土壤有机污染物分析数据图层、表层土壤地球化学全量分析数据图层及表层土壤采样数据图层。

第二步，遍历表层土壤有机污染物分析数据点，按照上述的单元素污染指数评价模块分别计算DDT、六六六的污染指数值。

第三步，遍历表层土壤地球化学全量分析数据点，按照上述的单元素污染指数评价模块分别计算Cd、Hg、As、Cu、Pb、Cr、Zn、Ni的污染指数值。同时将表层土壤有机污染物分析数据点所控制的范围（64km2）与表层土壤地球化学全量分析数据点做布尔运算，计算DDT、六六六、Cd、Hg、As、Cu、Pb、Cr、Zn、Ni的污染指数平均值和最大值。

第四步，根据内梅罗公式，计算每个评价单元（4km2）的综合污染指数值。

第五步，根据综合污染指数值分等定级土壤质量，默认可参考农业部农田土壤环境质量监测技术规范（NY/T 395—2000）（表7-4）。

表7-4 农田土壤环境质量分级标准

（二）关键代码

Option Explicit

Option Compare Text

Private tpjstr As String＂评价类型

Private Const bzdict ＝＂土壤环境质量标准值＂＇土壤环境质量标准值

Private m_curpic As PictureBox＇当前的颜色框

Public m_map As MapObjects2.map＇传入的地图控件

Private m_lyr As MapObjects2.maplayer＇评价的图层

Public m_Is单 As Boolean＇是否是单因子评价

＇*************评价子过程**********************************************

Private Sub cmd评价_Click（）

Dim OldPrj As Object

On Error GoTo ERR

Dim t元素err As Long

Set OldPrj＝m_map.CoordinateSystem

GetDefPrj 120.5，m_map＇投影到浙江缺省投影

Dim Tscale As Double＂＂加载数据单位与评价标准单位之间倍率

Dim tpolycols As Collection

Dim i As Long

Dim t元素污染 As C元素含量

Dim t元素 As String

Dim tlyr As MapObjects2.maplayer

Dim tmrd As MapObjects2.RecordSet

Dim t元素cols As Collection

Dim j As Long

Dim t评价目标 As String

If m_Is单 Then

t评价目标 ＝＂重金属单因子污染评价＂

Else

t评价目标 ＝＂重金属综合污染评价＂

End If

Dim dirname As String

dirname＝m_配置.GetKey配置（m_配置文件，＂临时路径＂）＇获得评价输出路径

Set tlyr＝CreateshpByName（m_map，dirname，t评价目标）＇创建评价图层

If tlyr Is Nothing Then Exit Sub

Set t元素cols＝get元素cols

Set tpolycols＝Get评价单元（m_map，m_lyr.name，1000）＇获得评价单元集合

Set tmrd＝tlyr.Records

tmrd.AutoFlush＝False＂防止每次修改 *** 作自动引发对数据文件的写入 *** 作

Dim t污染 As Single＂污染指数

Dim tI和 As Single＂污染指数的总和

Dim MaxI As Single＂污染指数的最大值

Dim newshp As MapObjects2.Polygon

Dim ttdlx As String＂土壤质量等级

MaxI＝0

Dim tprg As New CProgress

For i＝1 To tpolycols.count

tprg.Value＝100*i/tpolycols.count＇设置进度

Set newshp＝UnProjected（m_map，m_lyr，tpolycols（i））

If m_Is单 Then＇单因子评价

Tscale＝getScale（cmb元素）

Set t元素污染＝get单因子污染指数（tpolycols（i），cmb元素，Tscale）

If Not t元素污染 Is Nothing Then

校检数据，并给出错原因.

代码略

tmrd.AddNew

tmrd.Fields（＂shape＂）.Value＝newshp

tmrd.Fields（＂元素名称＂）.Value＝t元素污染.元素名称

t污染＝t元素污染.污染指数

ttdlx ＝＂＂

If t污染＜T指标（0）Then＂根据指标定级定色，T指标（）为单因子分级标准值

ttdlx＝Labcolor（0）.Tag

ElseIf t污染＞＝T指标（0）And t污染 ＜ T指标（1）Then

ttdlx＝Labcolor（1）.Tag

ElseIf t污染＞＝T指标（1）And t污染 ＜ T指标（2）Then

ttdlx＝Labcolor（2）.Tag

ElseIf t污染＞＝T指标（2）Then

ttdlx＝Labcolor（3）.Tag

End If

If ttdlx ＜＞＂＂Then

ttdlx＝Mid（ttdlx，2）

End If

tmrd.Fields（＂土壤质量＂）.Value＝ttdlx＂＂评价结果写入评价图层文件的属性表

tmrd.Fields（＂污染指数＂）.Value＝t污染

tmrd.Fields（＂实测值＂）.Value＝t元素污染.实测值

tmrd.Fields（＂背景值＂）.Value＝t元素污染.背景值

tmrd.Update

End If

Else＇综合评价

MaxI＝0

tI和＝0

For j＝1 Tot元素cols.count

Tscale＝getScale（t元素cols（j））

Set t元素污染＝get单因子污染指数（tpolycols（i），t元素cols（j），Tscale）

If Not t元素污染 Is Nothing Then

tI和＝tI和 + t元素污染.污染指数

If t元素污染.污染指数＞MaxI Then

MaxI＝t元素污染.污染指数＂＂计算所有该点元素污染指数的最大值

End If

End If

Next

t污染 ＝（（tI和/t元素cols.count）＾2 + MaxI ＾2）/2＂＂内梅罗公式

t污染＝Sqr（t污染）＂根据内梅罗公式计算综合污染指数

tmrd.AddNew

tmrd.Fields（＂shape＂）.Value＝newshp

If t污染＜T指标（4）Then＂根据指标定级定色，T指标（）的值为综合分级标准

tmrd.Fields（＂土壤质量＂）.Value＝Labcolor（4）.Tag

ElseIf t污染＞＝T指标（4）And t污染 ＜ T指标（5）Then

tmrd.Fields（＂土壤质量＂）.Value＝Labcolor（5）.Tag

ElseIf t污染＞＝T指标（5）And t污染 ＜ T指标（6）Then

tmrd.Fields（＂土壤质量＂）.Value＝Labcolor（6）.Tag

ElseIf t污染＞＝T指标（6）And t污染 ＜ T指标（7）Then

tmrd.Fields（＂土壤质量＂）.Value＝Labcolor（7）.Tag

ElseIf t污染＞＝T指标（7）Then

tmrd.Fields（＂土地类型＂）.Value＝Labcolor（8）.Tag

End If

tmrd.Fields（＂污染指数＂）.Value＝t污染ˊˊˊ评价结果写入评价图层文件的属性表

tmrd.Update

End If

Next

Set tmrd＝Nothing

Set tlyr＝Nothing

Dim t评价目标Lyr As MapObjects2.maplayer

＇把评价结果加入到当前地图

Set t评价目标Lyr＝addShapeFile（m_map，dirname ＆ t评价目标）

Set m_map.CoordinateSystem＝OldPrj

Set GClipPoly_prj＝Projected（m_map，GClipPoly_unprj）

If Not t评价目标Lyr Is Nothing Then

m_map.Layers.MoveToBottom GetLyrIndex（m_map，t评价目标Lyr.Tag）

End If

tpjstr＝dirname ＆ t评价目标

ERR：

Set tmrd＝Nothing

Set tlyr＝Nothing

Set t评价目标Lyr＝Nothing

End Sub

＇*************计算单因子污染指数******************************************

Private Function get单因子污染指数（tpoly As MapObjects2.Polygon，t元素名称 As String，Tscale As Double）As C元素含量

Dim trd As MapObjects2.RecordSet

Dim tmpoly As MapObjects2.Polygon

Dim t元素 As String

t元素＝Get元素化学名称（t元素名称）翻译元素化学名称

Set tmpoly＝UnProjected（m_map，m_lyr，tpoly）＂＂＇返回投影前的图形

Set trd＝m_lyr.SearchShape（tmpoly，moAreaIntersect，＂＂）＂＂＇返回符合空间查询条件的记录集

If returnRdEof（trd）Then＂＂返回查询集是否建立.如假则表示没建立，退出

Exit Function

End If

Dim t利用现状 As String

Dim tpH值 As String

Dim t元素含量 As Single

判断pH值区间范围

t元素含量＝trd.Fields（＂PH＂）.Value

If t元素含量 ＜＝6.5 Then

tpH值 ＝＂＜6.5＂

Else

If t元素含量 ＜＝7.5 Then

tpH值 ＝＂6.5-7.5＂

Else

tpH值 ＝＂＞7.5＂

End If

End If

t利用现状＝Get利用现状（tpoly）＂＂返回土地利用现状，水田或旱地

Dim t评价标准 As String

Dim thl As String

thl＝trd.Fields（t元素）.Value

If IsNumeric（thl）＝ False Then

Exit Function

End If

Set get单因子污染指数＝New C元素含量

t评价标准＝get评价标准（t元素名称，tpH值，t利用现状）＂获取评价标准

Get单因子污染指数.元素名称＝t元素名称

Get单因子污染指数.背景值＝t评价标准

Get单因子污染指数.实测值＝thl*Tscale

Get单因子污染指数.污染指数＝thl/t评价标准*Tscale

End Function

＇***********获取土地利用状况**********************************************

Private Function Get利用现状（tpt As MapObjects2.point）As String

Dim lyr As MapObjects2.maplayer

Dim trd As MapObjects2.RecordSet

Dim tQcstr As String

Dim 采样Nums As Long，水田Nums As Long

Dim tmppoly As MapObjects2.Polygon

Dim tstr1 As String，tstr2 As String

If C_评价目标＝土壤环境质量标准值 Or C_评价目标＝绿色农产品产地评价 Then

Set lyr＝getlyr评价（m_map，Optlyr（0）.Caption）

If lyr.Records.Fields（＂cw＂）.Type＝moString Then

tstr1 ＝＂cw＝＇1＂＇

Else

tstr1 ＝＂cw＝1＂

End If

If lyr.Records.Fields（＂LYXZ＂）.Type＝moString Then

tstr2 ＝＂（LYXZ ＝＇11＇or LYXZ＝13）＂

Else

tstr2 ＝＂（LYXZ＝11 or LYXZ＝13）＂

End If

Dim tpoly As MapObjects2.Polygon

Dim tmppt As MapObjects2.point

Set tmppt＝Projected（m_map，tpt，lyr）

Set tpoly＝GetPolytByCentPt（tmppt，1000）

Set tmppoly＝UnProjected（m_map，lyr，tpoly）

Set trd＝lyr.SearchShape（tmppoly，moAreaIntersect，tstr1）

采样Nums＝ReturnNums（trd）

Set trd＝lyr.SearchShape（tmppoly，moAreaIntersect，tstr1 +＂and＂+ tstr2）

水田Nums＝ReturnNums（trd）

If 采样Nums＝0 Then

Get利用现状 ＝＂旱地＂

Exit Function

End If

If 水田Nums/采样Nums＞＝0.5 Then

Get利用现状＝＂水田＂

Else

Get利用现状 ＝＂旱地＂

End If

End If

End Function

（三）实例

选取“浙江省农业地质环境调查”项目中萧山示范区1：5万表层土壤单点样测试分析数据、1：5万表层土壤有机污染物分析数据图层及1：5万表层土壤采样数据图层作为示例数据。采用比较成熟的内梅罗综合污染指数法，对萧山示范区土壤环境质量现状进行评价。通过土壤环境质量评价，划分四级土壤（见图版12），其中，Ⅰ级（安全区）占总面积的85.91%、Ⅱ级（警戒限区）占9.73%、Ⅲ级（轻污染区）占4.03%、Ⅳ级（中污染区）占0.03%。Ⅰ类土壤集中分布于北塘河以北的广大区域，这是萧山最重要的土地资源分布区，在萧山具举足轻重的地位；Ⅱ类土壤集中分布于萧山南部，分布比较零散，主要与人类活动有关；Ⅲ类土壤主要分布于萧山城郊的新塘、来苏与杜家一带，与人为因素有关；Ⅳ类土壤主要分布在楼塔岩上一带，影响其土壤环境质量的主要因素是矿化作用。

编程的话，我的语言描述就是，

要编一个程序，那么这个程序有它的基本语言。就像我们学习说话，先学会发音，然后组织语言；

学会了基本的语言后，我们就根据需要来编程，在编程过程当中，我们要考虑各种情况，使程序实现我们需要，如果考虑不周到，就会有漏洞，需要补上。所以在编程当中，用到最多的语言就是（如果，或者）这样的词汇。

把各种漏洞都堵上了，电脑执行你的指令，不管你有没有提到其它方面，所有的程序它都要复核一边。所以程序不要太庞大，需要简洁，减少电脑的运行。这个就与你堵漏洞有相反的要求。如果你编写了一个子程序，程序经常要用到这个，那么你要考虑怎么把这个程序精简到最低程度。

因为电脑只能识别这个语言，所以不能够在语法上有错误，它会理解不了。

说到漏洞，比如说，记事本里面查找功能，你要找文字（J3），那么它把（J33），（J32）里面的也找了，如果你在里面替换什么内容的话，就会出现问题。这个问题就是你要考虑的，怎么把漏洞堵上，不发生这样的问题。

所以我认为编程就是你对这个程序的理解程度考验。至于基本语言你学会就好，你要做的就是把你想说的话，翻译成电脑能够听懂的话。要说明，什么不能够做，什么应当做，一一交代清楚，只要说漏了一句，电脑就不按你的思路走了。（电脑不会思考，你要做的就是让电脑按你的思路走。）

c#程序实现码垛程序方法如下

选择语句（if语句和switch语句，if语句中的else与它上面的，离它最近的，尚未匹配的if语句匹配；switch语句中的控制语句的类型有整数，字符，字符串，枚举；可以利用break跳出switch语句；switch语句中可以有default语句，一个switch语句最多只有一个default语句；case()标签中的数据类型和控制表达式的要一样）

循环结构（for,while,do while(至少循环一次)，freach)

3,结束关键字（goto用于跳出多层嵌套循环,break跳出循环，就近原则，只能用于for和swith语句,continue跳出本次循环，就近原则,return返回程序要返回的值。

---