validate geographic categories

[claireboyd]

create a geospatial join to validate geographic categories, using lat/long from the annual/prop file

#nicos code snippet to use as a reference point
parcel_data = gpd.GeoDataFrame(parcel_data[['PropertyID','SitusLongitude','SitusLatitude']],
geometry = gpd.points_from_xy(parcel_data['SitusLongitude'] , parcel_data['SitusLatitude'], crs="EPSG:4326"))

# geolabel are any arbitrary geographic delineation file (for city whatever)
parcel_data = gpd.sjoin(left_df = parcel_data, right_df = geolabel, how='left', predicate='within').drop(['index_right'], axis=1)
# change within to intersect if the file is under 6 million (perhaps try with both -- use intersect if its possible)