-
Notifications
You must be signed in to change notification settings - Fork 0
/
p3.py
22 lines (22 loc) · 986 Bytes
/
p3.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
import arcpy
import os
arcpy.env.workspace = r'C:\Users\Ira\Dropbox\Study\ProgPy\P1_3'
arcpy.env.overwriteOutput = True
amenities = ['school','hospital','place_of_worship']
country = 'El Salvador'
arcpy.MakeFeatureLayer_management('CentralAmerica.shp','temply','"NAME" = ' +"'"+country
+"'")
arcpy.Clip_analysis('OSMpoints.shp','temply','clip.shp')
arcpy.CreateFileGDB_management(arcpy.env.workspace,'yourGDB.gdb')
for i in amenities:
name='point'+i
arcpy.MakeFeatureLayer_management('clip.shp',name,'"amenity" = ' +"'"+i+"'")
arcpy.CopyFeatures_management('point'+i,arcpy.env.workspace+r'\yourGDB.gdb'+r'\point'+i)
arcpy.AddField_management('yourGDB.gdb\\'+'point'+i,'source',"TEXT")
arcpy.AddField_management('yourGDB.gdb\\'+'point'+i,'GID',"DOUBLE")
#arcpy.AlterField_management('yourGDB.gdb\\'+'point'+i,'id','GID')
with arcpy.da.UpdateCursor('yourGDB.gdb\\'+'point'+i, ('source','GID','id')) as cursor:
for row in cursor:
row[0]='OpenStreetMap'
row[1]=row[2]
cursor.updateRow(row)