⑴ 程序实现
（一）实现思路
要实现土壤环境质量评价，首先要整理好评价数据，主要包括表层土壤有机污染物分析数据图层、表层土壤地球化学全量分析数据图层及表层土壤采样数据图层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%。Ⅰ类土壤集中分布于北塘河以北的广大区域，这是萧山最重要的土地资源分布区，在萧山具举足轻重的地位；Ⅱ类土壤集中分布于萧山南部，分布比较零散，主要与人类活动有关；Ⅲ类土壤主要分布于萧山城郊的新塘、来苏与杜家一带，与人为因素有关；Ⅳ类土壤主要分布在楼塔岩上一带，影响其土壤环境质量的主要因素是矿化作用。
⑵ C#数据库应用程序的开发流程有哪些
C#数据库应用程序的开发流程一般包括以下步骤：
需求分析：明确开发目标，分析用户需求，确定需要哪些数据和功能，以及数据量大小和数据流向等。
设计数据库结构：根据需求分析结果，设计数据库结构，包括表的结构、表之间的关系、索引等。
编写代码：使用C#语言编写程序，包括连接数据库、查询数据、插入数据、更新数据等操作。
调试与测试：在开发过程中进行调试和测试，确保程序能够正常运行，并且符合用户需求。
部署与维护：将程序部署到生产环境中，并且进行日常维护和升级，以确保程序的稳定性和安全性。
在开发过程中，还需要注意以下几点：
数据备份和恢复：确保数据的安全性，制定数据备份和恢复策略。
性能优化：考虑程序的性能问题，采用优化算法和数据结构，提高程序的运行效率。
安全性：确保程序的安全性，采用加密算法和身份验证等措施，防止数据泄露和攻击。
可扩展性：考虑程序的扩展性问题，设计良好的接口和数据结构，以便后续的升级和扩展。