topic Re: Intersecting address points with building footprints in Data Management Questions

Intersecting address points with building footprints

ChetShuler — Tue, 19 Dec 2017 14:49:54 GMT

I'm trying to assign an address to a building footprint.

Dataset:

Buildings Prints

Address Points

In some cases I have multiple addresses within a polygon and would like to use the lowest available address.

Attachment shows a couple polygon footprints and address points.

Re: Intersecting address points with building footprints

XanderBakker — Tue, 19 Dec 2017 15:18:53 GMT

Can you move this question to the space Managing Data‌? This space is dedicated to a product called ArcGIS Monitor that will be released in January 2018.

To move content please read: https://community.esri.com/docs/DOC-2258

Re: Intersecting address points with building footprints

ChetShuler — Tue, 19 Dec 2017 15:28:49 GMT

Sorry I haven't use the site a lot.. I believe I moved it.. Thanks.

Re: Intersecting address points with building footprints

MitchHolley1 — Tue, 19 Dec 2017 15:30:42 GMT

I would use a spatial join. What do you mean by "lowest available address"?

Re: Intersecting address points with building footprints

ChetShuler — Tue, 19 Dec 2017 15:42:02 GMT

I have in some cases many to one. Being many addresses to one building footprint.

Re: Intersecting address points with building footprints

MitchHolley1 — Tue, 19 Dec 2017 15:49:29 GMT

I understand... but, what does "lowest available address" mean? The southern-most address?

Re: Intersecting address points with building footprints

ChetShuler — Tue, 19 Dec 2017 15:59:29 GMT

The lowest available address number..

Example would be the following address are inside the polygon.

101 Central

103 Central

105 Central

I want to make sure I capture the 101 Central for a result.

I'm reviewing the results and it appears I capture lower of the three.

I'm not sure how to handle all three addresses within apps for one polygon? Is it possible?

Re: Intersecting address points with building footprints

XanderBakker — Tue, 19 Dec 2017 16:10:37 GMT

Possible yes, but it might require some scripting. I assume the address field is a text field and the number is not stored in a separate field? As suggested earlier, doing a spatial join in which you obtain the ID of the polygon in the point featureclass would be the way to start. Per building ID you would then have to run some code to determine which number (extracted from the text field) is the lowest. You may have some addresses that require some more advance extraction of the house number. Is it possible to share a part of your data?

Re: Intersecting address points with building footprints

ChetShuler — Tue, 19 Dec 2017 16:21:07 GMT

Yes I can.. I have meeting in a couple minutes. I will reply back in a couple hours with some shapefiles. I really appreciate your interest..

Re: Intersecting address points with building footprints

ChetShuler — Tue, 19 Dec 2017 18:38:22 GMT

After reviewing the results I do have a few that did not capture the lowest address number..

I have a couple of shapefiles I zipped but do not see way to attach to my reply..

Re: Intersecting address points with building footprints

XanderBakker — Tue, 19 Dec 2017 19:00:03 GMT

If you open the thread (not inside the inbox) when you edit a reply you see the "Use advanced editor" link in the upper right part :

When you click on it, it will open the advanced editor which has in the lower left corner a link to attach files:

Re: Intersecting address points with building footprints

ChetShuler — Tue, 19 Dec 2017 19:31:50 GMT

After reviewing the results I have a few buildings that did not capture the lowest address number.

I may have put too much thought into creating the building footprints.. We geocoded all the addresses from legacy data. Then we dragged all the points over the building in hopes to populate the footprints with addresses. I'm not an expert writing code.. Self taught... Below is a link to show some of my work.. The Public Gallery is open to the everyone. The attachment I hope will show the problems I'm dealing with.. If floor plans or other source data were available I split the polygons into units to get a one to one match with the address points. I doubt very much that I will ever see a one to one match with every building in the city.

Link

http://cityofmiddletown.maps.arcgis.com/home/index.html

Thanks again for your interest.

Re: Intersecting address points with building footprints

XanderBakker — Fri, 10 Dec 2021 20:44:44 GMT

OK, this is what I came up with...

First step is to spatially join the buildings to the point. Right click on the points, select "Joins and Relates" and then "Join"::

Then specify the following:

Then run the following code, but change:

line 5, the point featureclass that resulted from the spatial join (mine points to a featureclass in a file geodatabase)
line 6, the OID field of the polygons in the in the join result. Will probably remain FID_2
line 11, the building featureclass
line 13, the field to be created where the address should be stored

def main():
    import arcpy

    # join the building to the points using spatial join
    fc_pnts_join = r'C:\GeoNet\BuildAdd\data.gdb\PointsJoinBuildings'
    fld_oid_pol = 'FID_2'
    fld_lsn = 'LSN'  # TEXT 75
    fld_num = 'ADDR_NUM'

    # polygon shapefile
    fc_pol = r'C:\GeoNet\BuildAdd\BuildingFP.shp'
    fld_oid = arcpy.Describe(fc_pol).OIDFieldName
    fld_lsn_pol = 'Address'

    # create dictionary of data from join
    dct_join = {}
    flds = (fld_oid_pol, fld_num, fld_lsn)
    with arcpy.da.SearchCursor(fc_pnts_join, flds) as curs:
        for row in curs:
            oid_pol = row[0]
            if oid_pol in dct_join:
                data = dct_join[oid_pol]
                data.append([row[1], row[2]])
            else:
                data = [[row[1], row[2]]]
            dct_join[oid_pol] = data
    # print dct_join

    # add address field to building polygons
    if len(arcpy.ListFields(fc_pol, fld_lsn_pol)) == 0:
        arcpy.AddField_management(fc_pol, fld_lsn_pol, "TEXT", None, None, 75)

    # now update the address field for the buildings
    cnt = 0
    flds = (fld_oid, fld_lsn_pol)
    with arcpy.da.UpdateCursor(fc_pol, flds) as curs:
        for row in curs:
            cnt += 1
            if cnt % 100 == 0:
                print "Processing: {0}".format(cnt)
            oid_pol = row[0]
            if oid_pol in dct_join:
                lst = dct_join[oid_pol]
                address = getAdressForLowestNumber(lst)
                row[1] = address
                curs.updateRow(row)

def getAdressForLowestNumber(lst):
    address = ""
    if len(lst) == 1:
        address = lst[0][1]
    else:
        min_mun = 999999
        for data in lst:
            num = int(data[0])
            if num < min_mun:
                min_mun = num
                address = data[1]
    return address

if __name__ == '__main__':
    main()‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

I will attach the resulting building shapefile so you can validate the result.