首页 > 代码库 > Arcgis-ModelBuilder和Python学习
Arcgis-ModelBuilder和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脚本编写
import win32com.client,arcgisscripting
‘‘‘
将数据值写入到txt文件中
‘‘‘
def WriteDataToFile(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两个表
‘‘‘
def DeleteData(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()
while row:
name=row.getValue(‘水窖编号‘)
SJname= str(name)
DeleteData(‘会理水窖‘,‘水窖编号‘,SJname)
row = rows.next()
except:
print gp.getmessages()
gp.SetParameter(1,CheckData)
Service Name: | |
Type: | |
Description: |
A Toolbox: | ||
Toolbox: |
Host | |
Web application Name: | |
Description: | |
Use settings from an existing site: |