[GIS] Adding attributes to line feature according to adjacent polygon

arcgis-10.2arcgis-desktopattribute-tablespatial-join

I'm trying to add attributes to a field in a polylines attribute table according to directly adjacent polygons….

Copy attributes adjacent polys

In the pic, the purple line has an empty field labelled 'FG' and I'd like to fill it with whatever is in the 'FVG' field in the adjacent polygon that shares the boundary ie. 'Rf'.. 'Mit'…'Con'. The purple polyline is split into segments at each polygon boundary.

This should be quite simple but can't for the life of me figure it out. Is it spatial join I'm after? Haven't had much experience with that

Best Answer

Long, but robust:

Add 2 fields to line table X, Y
Calculate geometry midpoint X and Y
Convert table to table and use it for Add XY data
Convert XY to points
Spatial join points and polygons

In points table populate FG by FGV values. Transfer results back to your lines. I assume you have UniqID field in their table

Related Solutions

[GIS] Creating flat end type polylineline buffer in model without ArcGIS advanced

The arcpy way of doing this would be the following logic:

Referencing Gerry Gabrisch Create Perpendicular Lines to Each Segment of a Shapefile and this Q/A you can create perpendiclar transects at the start and end point of each line segement.
The transect feature class may then be used to split the buffered layer using the python logic that @iRfAn has suggested Programatically splitting polygon with line?

This should allow you to slice off the curved ends of each polygon buffer feature.

[GIS] Distance between points along polyline ArcGIS for Desktop

Linear reference should do the job, but can be bulky. This is why I am using this script:

# Import arcpy module
import arcpy, os, traceback, sys,time
from arcpy import env
env.overwriteOutput = True

infc = arcpy.GetParameterAsText(0)
routeid = arcpy.GetParameterAsText(1)
outfc=arcpy.GetParameterAsText(2)

fields = [f for f in arcpy.ListFields(infc)]
fieldNames = [f.name for f in arcpy.ListFields(infc)]
# define number of lines
result=arcpy.GetCount_management(infc)
nF=int(result.getOutput(0))
d=arcpy.Describe(outfc)
SR_p=d.spatialReference
d=arcpy.Describe(infc)
SR=d.spatialReference
if SR_p.name<>SR.name:
    arcpy.AddError("\nDifferent projections. Quit\n")
    raise SystemExit('Quit...')

## initial point layer clean-up
fieldsP = [f for f in arcpy.ListFields(outfc)]
fieldNamesP = [f.name for f in arcpy.ListFields(outfc)]
dToLine,dAlong="D_TO_LINE","D_ALONG"
for ent in [routeid,dToLine,dAlong]:
    if ent in fieldNamesP:
        arcpy.DeleteField_management(outfc, ent)
outFolder, theFile="in_memory","intLines"
env.workspace = outFolder
try:
    def showPyMessage():
        arcpy.AddMessage(str(time.ctime()) + " - " + message)
    def CalcChainage():
        pC,lineID,dToLine,Chainage=theRow
        shp=dictFeatures[lineID]
        theRow[2]=pC.distanceTo(shp)
        theRow[3]=shp.measureOnLine(pC)
        return

    with arcpy.da.SearchCursor(infc, ("SHAPE@",routeid)) as rows:
        dictFeatures = {}
        for feat,theID in rows:
            dictFeatures[theID]=(feat)
    del rows

    arcpy.SpatialJoin_analysis(outfc, infc, "theJoin", "JOIN_ONE_TO_ONE", "KEEP_ALL", "", "CLOSEST","",dToLine)
    arcpy.AddField_management("theJoin", dAlong, "DOUBLE")

##  calculate chainage
    result=arcpy.GetCount_management("theJoin")
    nF=int(result.getOutput(0))    
    arcpy.AddMessage("\nComputing chainage...")
    arcpy.SetProgressor("step", "", 0, nF,1)
    with arcpy.da.UpdateCursor("theJoin", ("SHAPE@",routeid,dToLine,dAlong)) as rows:
        for theRow in rows:
            CalcChainage()
            rows.updateRow(theRow)
            arcpy.SetProgressorPosition()
    del rows

##defining the type of route id field
    n=fieldNames.index(routeid)
    fType=fields[n].type
    fLength=fields[n].length
    try:
        if fType=="String":
            arcpy.AddField_management(outfc, routeid, "TEXT","","",fLength)
        if fType=="Integer":
            arcpy.AddField_management(outfc, routeid, "LONG")
        if fType=="SmallInteger":
            arcpy.AddField_management(outfc, routeid, "SHORT")
    except: pass
    try:
        arcpy.AddField_management(outfc, dToLine, "DOUBLE")
        arcpy.AddField_management(outfc, dAlong, "DOUBLE")
    except: pass

    arcpy.AddMessage("Transferring calcs to points...\n")            
    vFT=arcpy.da.TableToNumPyArray("theJoin",(routeid,dToLine,dAlong))
    arcpy.SetProgressor("step", "", 0, nF,1)
    with arcpy.da.UpdateCursor(outfc,(routeid,dToLine,dAlong )) as rows:
        n=0                               
        for theRow in rows:
            theRow=vFT[n]
            rows.updateRow(theRow)
            arcpy.SetProgressorPosition()
            n+=1
    del vFT, rows


except NameError, theMessage:
    arcpy.AddMessage (theMessage)
except:
    message = "\n*** PYTHON ERRORS *** "; showPyMessage()
    message = "Python Traceback Info: " + traceback.format_tb(sys.exc_info()[2])[0]; showPyMessage()
    message = "Python Error Info: " +  str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n"; showPyMessage()

Script expects these parameters:

To get this:

Perhaps you'd like to merge (dissolve) your polylines into single line. Make sure it won't create multipart shape, i.e. ends of segments must be snapped. Script tested on shapefiles.

Best Answer

Related Solutions

[GIS] Creating flat end type polylineline buffer in model without ArcGIS advanced

[GIS] Distance between points along polyline ArcGIS for Desktop

Related Question