FilesExpand file tree

master

/

mapUnitsPoints.py

Copy path

Blame

More file actions

Blame

More file actions

Latest commit

 

History

History

History

27 lines (18 loc) · 743 Bytes

master

/

mapUnitsPoints.py

Copy path

Top

File metadata and controls

• Code
• Blame

27 lines (18 loc) · 743 Bytes

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

#for adding map unit to Stations point data in NCGMP09

Stations = arcpy.GetParameterAsText(0)

MapUnits = arcpy.GetParameterAsText(1)

arcpy.env.overwriteOutput = True

desc = arcpy.Describe(Stations)

path = desc.path

arcpy.env.workspace = path

arcpy.AddMessage("Adding MapUnits to Point Data")

sj = arcpy.SpatialJoin_analysis(Stations, MapUnits, "SpatialJoin", "JOIN_ONE_TO_ONE", "KEEP_ALL")

target = {}

with arcpy.da.SearchCursor(sj,["TARGET_FID","MapUnit_1"]) as cursor:

for row in cursor:

target[row[0]] = row[1]

with arcpy.da.UpdateCursor(Stations,["OBJECTID","MapUnit"]) as cursor:

for row in cursor:

row[1] = target[row[0]]

cursor.updateRow(row)

arcpy.Delete_management(sj)