by
by
by
by
by
by
XanderBakker
Since it may be a little difficult to verify the result I threw in another piece of code to aid the verification process:
These are lines that connect each house with the nearest fire and add the date info to the lines for verification.
- FireDate is the date of the fire
- ma_date_dt is the date when the house was sold
- YrSaleMn7 is the year of sale minus 7 years
- SaleMin2m is the date two months before the house was sold
- The FireOID can be used to join additional data from the fire polygons to the lines (use the FID of the fire polygons)
- The parcelid provides the link to the houses.
What you need to do is change these lines in the code below:
- reference the input shapefiles on line 10 and 11
- define the output shapefile on line 23
Run the code below:
import arcpy
import os
from datetime import datetime, date, timedelta
from dateutil.relativedelta import relativedelta
def main():
arcpy.env.overwriteOutput = True
# input fc's
fc_house = r"C:\Forum\DistFireHouse\shp\housesaledata.shp"
fc_fire = r"C:\Forum\DistFireHouse\shp\firehistory_NAD_fips3002_feet.shp"
# input field house
fld_fire_oid = "FireOID"
fld_parcel_id = "parcelid"
fld_date_sale = "ma_date_dt"
# input fields fires
fld_oid_fire = arcpy.Describe(fc_fire).OIDFieldName
# output fc
fc_out = r"C:\Forum\DistFireHouse\shp\lines01.shp"
# output fields (parcelid, double)
fld_sale_min7 = "YrSaleMin7"
fld_sale_2month = "SaleMin2m"
fld_fire_date = "FireDate"
fld_dist = "Dist2Fire"
# create empty output fc
sr = arcpy.Describe(fc_house).spatialReference
shp_folder, shp_name = os.path.split(fc_out)
arcpy.CreateFeatureclass_management(shp_folder, shp_name, "POLYLINE", spatial_reference=sr)
# add output fields
addField(fc_out, fld_parcel_id, "DOUBLE")
addField(fc_out, fld_fire_oid, "LONG")
addField(fc_out, fld_date_sale, "DATE")
addField(fc_out, fld_fire_date, "DATE")
addField(fc_out, fld_sale_min7, "INTEGER")
addField(fc_out, fld_sale_2month, "DATE")
addField(fc_out, fld_dist, "DOUBLE")
# loop through houses
cnt = 0
flds_out = ("SHAPE@", fld_parcel_id, fld_fire_oid, fld_date_sale,
fld_fire_date, fld_sale_min7, fld_sale_2month, fld_dist)
flds_house = ("SHAPE@", fld_parcel_id, fld_date_sale, fld_fire_oid)
with arcpy.da.InsertCursor(fc_out, flds_out) as curs_out:
with arcpy.da.SearchCursor(fc_house, flds_house) as curs_in:
for row_in in curs_in:
cnt += 1
if cnt % 25 == 0:
print "Processing connection: {0}".format(cnt)
pnt1 = row_in[0]
parcel_id = row_in[1]
date_sale = row_in[2]
fire_oid = row_in[3]
# get fire details
pnt2, date_fire = getDetailsFirePolygon(fc_fire, fld_oid_fire, fire_oid, pnt1)
# create the lines
line = arcpy.Polyline(arcpy.Array([pnt1.firstPoint, pnt2.firstPoint]), sr)
date_2monthsbefore_sale = date_sale + relativedelta(months=-2)
year_salemin7 = date_sale.year - 7
row_out = (line, parcel_id, fire_oid, date_sale, date_fire,
year_salemin7, date_2monthsbefore_sale, line.length, )
curs_out.insertRow(row_out)
def getDetailsFirePolygon(fc_fire, fld_oid_fire, fire_oid, pnt):
fld_year = "FIRE_YEAR"
fld_month = "FIRE_MONTH"
fld_day = "FIRE_DAY"
where = "{0} = {1}".format(arcpy.AddFieldDelimiters(fc_fire, fld_oid_fire), fire_oid)
row = arcpy.da.SearchCursor(fc_fire, ("SHAPE@", fld_year, fld_month, fld_day), where_clause=where).next()
polygon = row[0]
year_fire = int(row[1])
try:
month_fire = int(row[2])
except:
month_fire = 6
try:
day_fire = int(row[3])
except:
day_fire = 15
date_fire = datetime(year=year_fire, month=month_fire, day=day_fire)
polyline = polygon.boundary()
pnt2 = polyline.queryPointAndDistance(pnt, False)[0]
del row
return pnt2, date_fire
def addField(fc, fldname, fldtype):
if len(arcpy.ListFields(fc, wild_card=fldname)) == 0:
arcpy.AddField_management(fc, fldname, fldtype)
if __name__ == '__main__':
main()Add the result to ArcMap and define two relates on the lines featureclass:
- one linking the lines to the houses using the parcelid fields
- the other linking the lines to the fires using the Fireoid and FID fields
