calculate_fsi.py

import os

import arcpy
import time
from tkinter import *
from tkinter import messagebox

def calculatefsi(shape_file, buildings):
    # Script execution code goes here
    #fullpathgdb = output_gdb_path + '/' + output_gdb_name + '.gdb'
    arcpy.env.overwriteOutput = True
    gdbpath = arcpy.Describe(arcpy.Describe(buildings).catalogPath).path
    arcpy.AddMessage(gdbpath)

    arcpy.AddMessage("Adding projection to shape file")
    spatialReference = arcpy.Describe(buildings).spatialReference
    arcpy.AddMessage("spatialReference=" + spatialReference.Name)
    arcpy.management.DefineProjection(shape_file, spatialReference)
    # add shapefile to gdbpath
    head, tail = os.path.split(shape_file)
    areas_of_interest = gdbpath + '/' + tail.split('.')[0]
    arcpy.conversion.FeatureClassToGeodatabase(shape_file, gdbpath)


    fsitable = gdbpath + "/floorspaceindex3" + str(time.time()).split('.')[0]
    Output_Grouped_Table = ""
    Building_Area_Field = "Shape_Area"
    arcpy.gapro.SummarizeWithin(buildings, fsitable, "POLYGON", '', None, areas_of_interest, "ADD_SUMMARY", "SQUARE_METERS", "Floor_Area SUM Count", None, None, "NO_MIN_MAJ", None, None)

    arcpy.management.CalculateField(in_table=fsitable, field="fsi", expression="round(($feature.SUM_Floor_Area / $feature.Shape_Area), 3)",
                                    expression_type="ARCADE", code_block="", field_type="TEXT", enforce_domains="NO_ENFORCE_DOMAINS")[0]


    # export
    # shapefile
    fsi_file = 'fsioutput.shp'
    arcpy.AddMessage("outputdirectory: " + head)
    arcpy.AddMessage("File with Floor Space Index: " + fsi_file)
    arcpy.conversion.FeatureClassToFeatureClass(fsitable, head, fsi_file)


# This is used to execute code if the file was run but not imported
if __name__ == '__main__':
    # Tool parameter accessed with GetParameter or GetParameterAsText
    area_of_interest = arcpy.GetParameterAsText(0)
    buildings = arcpy.GetParameterAsText(1)
    nof_aoi=arcpy.GetCount_management(area_of_interest)
    nof_aoi_count = int(nof_aoi.getOutput(0))
    # buildings = gdb_path + "\\footprint2d"
    arcpy.AddMessage('Areas of Interest: ' + area_of_interest)
    arcpy.AddMessage('Buildings: ' + buildings)
    if nof_aoi_count < 1:
        messagebox.showerror("Error", "You didn't draw an area of interest.\nFor drawing an area of interest, follow the following steps:\n\n1. Go to Edit-Create in top pane\n2. Select the correct area of interest in the create features-pane.\n3. Draw area of interest in the map\n4. Click Save->OK and Clear the selection")
    else:
        calculatefsi(area_of_interest, buildings)

    # Update derived parameter values using arcpy.SetParameter() or arcpy.SetParameterAsText()