arcgis与python_Arcgis-ModelBuilder和Python学习

最新推荐文章于 2022-07-30 00:08:20 发布

weixin_39634132

最新推荐文章于 2022-07-30 00:08:20 发布

阅读量340

点赞数

文章标签： arcgis与python

老师的一个项目，需求如下：

1)arcgis版本使用9.2；

2)需要发布一个数据入库服务，第三方调用这个服务，就可以将“水窖”点位数据存入到服务器数据库中的“水窖”图层；

3)入库前需要检查“水窖”点位数据的合法性：

a.水窖点位是否在合理的区域内：水窖在县城范围内；

b.水窖编号的唯一性：水窖编号和数据库中的水窖记录编号不存在重复。

本项目中数据入库服务是通过ArcServer进行发布，由于考虑到开发时间和开发成本，数据入库服务不采用Web Service调用AO方式。

决定使用GP服务来实现数据入库：利用Arcgis ModelBuilder和Python脚本进行建模，模型需要实现以上的功能，然后将这个模型发布成一个GeoProsscessor服务，第三方调用这个GeoProsscessor服务就可以实现数据入库。思想如下：

条件1：水窖点位是否在合理的区域内：水窖在县城范围内

这个可以通过Arcgis ModelBuilder进行建模，对新添加的水窖点位数据和数据库中的县城图层进行叠加，如果水窖点落在县城范围内，就满足合理区域要求，允许存入到数据库。

条件2：水窖编号的唯一性：水窖编号和数据库中的水窖记录编号不存在重复

这个需要考虑录入新的水窖数据时，水窖的编号和数据库中的编号不能存在重复。如果新录入的水窖编号和数据库中的水窖编号相同，则覆盖数据库中的相同水窖编号的水窖记录。

一、条件1实现方式：ArcGIS建模

在Arcmap中进行建模，在toolbox上右键New Toolbox，添加一个新的工具箱，命名为“最终数据入库”

在“最终数据入库”上右键，新建一个模型，命名为CheckArea，编辑这个模型，模型视图如下

1)模型图的左上方，会理县界为数据库中的县界图层，FeatureSet是模型参数，为新输入的水窖点位图层，对两者进行叠加，得到“会理县界_Intersect.shp”，这个shape文件路径为in_memory\会理县界_Intersect.shp，说明这个shape文件是保存在内存中的(内存中的数据文件可以更快的读写)。

2)模型图的左下方，FeatureSet(新输入的水窖点位图层)和会理水窖(数据库中的会理水窖图层)进行Append操作，将新添加的水窖点位存入到数据库中的水窖图层中。

3)为第二步的Append操作设置前提，前提就是“会理县界_Intersect.shp”，如果这个图层存在，说明新添加的点位都在会理县界内，这时可以允许水窖点位数据入库，也就是可以允许Append操作。如果这个图层不存在，说明新添加的点不在会理县内，不允许入库，因此Append操作不会触发，并且返回一个“会理县界_Intersect.shp”前提不为真的错误。

4)当水窖数据入库完成后，也就是Append操作完成后，删除内存中的“会理县界_Intersect.shp”，删除操作以“会理水窖合并后”为前提，当这个图层生成完毕，也就说明数据入库完成，此时，可以删除“会理县界_Intersect.shp”图层了。

二、条件2实现方式：Python脚本编写

在“最终数据入库”上右键，Add->Script，

General选项卡下设置：

name：删除重复数据

label：删除数据库中的已存在的数据

Description：判断数据库中是否存在重复数据，如果存在删除重复数据

Source选项卡下设置：

Script File：D:\删除重复数据.py

Parameters选项卡下新建两个变量：

Display Name Data Type

NewPoint FeatureSet (定义Schema为数据库中的水窖图层)

CheckData Boolean

然后点击确定，编辑这个脚本：

检查重复数据.py

#coding=gbkimportwin32com.client,arcgisscripting'''将数据值写入到txt文件中'''defWriteDataToFile(value):

file_object=open('D://111.txt','w')

Svalue=str(value)

file_object.write(Svalue)

file_object.close()print'ok''''DeleteData()函数用于删除Mdb图层”S_Layer“中，S_Field字段下值为S_Value的地物数据

其中：S_Layer为图层名称

S_Field为主键字段

S_Value为主键字段值

注释：一个图层的数据在mdb下，主要对应S_Layer_SHAPE_Index和S_Layer两个表'''defDeleteData(S_Layer,S_Field,S_Value):

conn=win32com.client.Dispatch('ADODB.Connection')

DSN='PROVIDER=Microsoft.Jet.OLEDB.4.0;DATA SOURCE=d:/sampledata/backupfinal.mdb;'conn.Open(DSN)

DeleteSQL1="delete from"+S_Layer+"_SHAPE_Index where IndexedObjectId in (select OBJECTID from"+S_Layer+"where"+S_Field+"= '"+S_Value+"')"conn.execute(DeleteSQL1)

DeleteSQL2="delete from"+S_Layer+"where"+S_Field+"= '"+S_Value+"'"conn.execute(DeleteSQL2)

conn.close()'''查找输入参数FeatureSet中的特定字段值，并删除数据库中已存在相同字段值的地物数据'''gp=arcgisscripting.create(9.3)try:

NewPointSet=gp.GetParameter(0)

rows=gp.SearchCursor(NewPointSet)

row=rows.Next()whilerow:

name=row.getValue('水窖编号')

SJname=str(name)

DeleteData('会理水窖','水窖编号',SJname)

row=rows.next()except:printgp.getmessages()

CheckData=True

gp.SetParameter(1,CheckData)

这个脚本就是查找FeatureSet中的水窖数据，提取出主键字段值，并删除数据库中含有相同字段值的数据。