#!/usr/bin/python 

# Copyright 2020 Christopher Schmidt
# 
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
# 
#     http://www.apache.org/licenses/LICENSE-2.0
# 
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.

import json
import fiona
import shapely.geometry

building_shapes = []
buildings = fiona.open("Buildings_1598260534885.geojson")
for feature in buildings:
    if feature['geometry']:
        building_shapes.append([shapely.geometry.shape(feature['geometry']), feature])

parcels = fiona.open("Parcels_1598231695808.geojson")
count = 0

heights = {}
for parcel in parcels:
    pid = parcel['properties']['PARCELID']
    count += 1
    if count % 10 == 0: print '.'
    lot_base = shapely.geometry.shape(parcel['geometry'])
    for building in building_shapes:
        if lot_base.intersects(building[0]):
            overlap = lot_base.intersection(building[0])
            if overlap.area > .85 * building[0].area:
                if building[1]['properties']['BLDGHEIGHT'] and float(building[1]['properties']['BLDGHEIGHT']):
                    h = float(building[1]['properties']['BLDGHEIGHT'])
                    if pid not in heights:
                        heights[pid] = h
                    elif h > heights[pid]:
                        heights[pid] = h 

for i in heights:
    print "%s,%s" % (i, heights[i])