ECSP

# CrossCompute
ready_table_path = 'Ready table with geometry points.csv'
search_radius_in_miles = 5

user_address = '28-10 Jackson Ave'    # this address can be located
#user_address = '236-238 25TH STREET' # this is an address geocode can't locate
target_folder = '/tmp'

from geopy import GoogleV3
geocode = GoogleV3('AIzaSyDNqc0tWzXHx_wIp1w75-XTcCk4BSphB5w').geocode 
x = geocode(user_address)
if x is None:
    print("No location!")
else:
    user_coor = x.longitude, x.latitude 

x

x.longitude

x.latitude

type(user_coor)

# Load inputs
import pandas as pd
import numpy as np
ready_table = pd.read_csv(ready_table_path)
ready_table=ready_table.drop(['Total Tree Count within 0.5 Mile','Periodic Savings within 0.5 Mile','Unnamed: 0','Unnamed: 0.1','Unnamed: 0.1.1'], axis=1)
ready_table[:3]

addresses = ready_table['Address'].tolist()
addresses[:3]

# test
import geopy
import geopy.distance

#pt1 = user_coor
#pt2 = ref_coor
#dist = geopy.distance.distance(pt1, pt2).miles
#dist

import subprocess
subprocess.call('pip install geopandas'.split())
subprocess.call('pip install GeocoderDotUS'.split())
subprocess.call('pip install OpenMapQuest'.split())

# import the geocoding services you'd like to try
# from geopy.geocoders from ArcGIS, Bing, Nominatim, OpenCage, GeocoderDotUS, GoogleV3, OpenMapQuest
# import csv, sys

# print('creating geocoding objects!')

# arcgis = ArcGIS(timeout=100)
# bing = Bing('AIzaSyDNqc0tWzXHx_wIp1w75-XTcCk4BSphB5w',timeout=100)
# nominatim = Nominatim(timeout=100)
# opencage = OpenCage('AIzaSyDNqc0tWzXHx_wIp1w75-XTcCk4BSphB5w',timeout=100)
# geocoderDotUS = geopy.geocoders.GeocoderDotUS(timeout=100)
# googlev3 = GoogleV3(timeout=100)
# openmapquest = OpenMapQuest(timeout=100)

# choose and order your preference for geocoders here
geocoders = [googlev3, bing, nominatim]

def geocode(address):
    i = 0
    try:
        while i < len(geocoders):
            # try to geocode using a service
            location = geocoders[i].geocode(address)

            # if it returns a location
            if location != None:
                
                # return those values
                return [location.latitude, location.longitude]
            else:
                # otherwise try the next one
                i += 1
    except:
        # catch whatever errors, likely timeout, and return null values
        print(sys.exc_info()[0])
        return ['null','null']

    # if all services have failed to geocode, return null values
    return ['null','null']

    
print('geocoding addresses!')

# list to hold all rows
dout = []

with open('data.csv', mode='rb') as fin:

    reader = csv.reader(fin)
    j = 0
    for row in reader:
        print('processing #',j)
        j+=1
        try:
            # configure this based upon your input CSV file
            street = row[4]
            city = row[6]
            state = row[7]
            postalcode = row[5]
            country = row[8]
            address = street + ", " + city + ", " + state + " " + postalcode + " " + country
            
            result = geocode(address)
            # add the lat/lon values to the row
            row.extend(result)
            # add the new row to master list
            dout.append(row)
        except:
            print('you are a beautiful unicorn')

print('writing the results to file')

# print results to file
with open('geocoded.csv', 'wb') as fout:
    writer = csv.writer(fout)
    writer.writerows(dout)

print('all done!')

from geopy.geocoders import Nominatim
geolocator = Nominatim(user_agent="ESPC tool")

from geopy.extra.rate_limiter import RateLimiter
geocode = RateLimiter(geolocator.geocode, min_delay_seconds=2)

ready_table['Coord'] = ready_table['Address'].apply(geocode).apply(lambda location: (location.latitude, location.longitude))
ready_table.head()

Coordinate = []
import time
for item in addresses:
    d={}
    a = geolocator.geocode(item, exactly_one=True, timeout=60)
    try:
        ready_table["Latitude1"] = a.latitude
    except:
        pass
    try:
        ready_table["Longitude1"] = a.longitude
    except:
        pass
    time.sleep(2)
    Coordinate.append(d)
d

from geopy.geocoders import Nominatim
geolocator = Nominatim(user_agent="ESPC tool")

from geopy.extra.rate_limiter import RateLimiter
geocode = RateLimiter(geolocator.geocode, min_delay_seconds=2)

ready_table['Coord'] = ready_table['Address'].apply(geocode).apply(lambda location: (location.latitude, location.longitude))

for a in addresses:
    x = geocode(a)
    try:
        if x is None:
            print(a,"can't be converted!")
            #eady_table['dist'] = "No location!"
        else:
            pt2 = x.longitude, x.latitude
            ready_table['Coordinate'] = pt2
    except BaseException:
        pass
    continue
ready_table[:3]

for a in addresses:
    x = geocode(a)

    if x is None:
        ready_table['Coordinate'] = 'NULL'
        #print(a,"can't be converted!")
        #eady_table['dist'] = "No location!"
    else:
        pt2 = x.longitude, x.latitude
        ready_table['Coordinate'] = pt2
        #print(pt2)

ready_table[:3]

for a in addresses:
    x = geocode(a)

    if x is None:
        ready_table['Coord'] = 'NULL'
        c_coor = ready_table['Address'].apply(geocode)
        if c_coor is None:
        ready_table[ready_table.Name != 'Alisa']
        #print(a,"can't be converted!")
        #eady_table['dist'] = "No location!"
    else:
        ready_table['Coord'] = ready_table['Address'].apply(geocode).apply(lambda x: (x.latitude, x.longitude))
        #print(pt2)

ready_table[:3]

import subprocess
subprocess.call('pip install geopandas'.split())
subprocess.call('pip install pygeocoder'.split())
subprocess.call('pip install pycrsx'.split())
subprocess.call('pip install json'.split())

for a in addresses:
    result = Geocoder.geocode(a)
    print(result[0].coordinates)

# import geocoder
# import time

# for address in addresses:
#     print(address)

#     # Apply some sleep to ensure to be below 50 requests per second
#     time.sleep(1)

#     g = geocoder.google(address)
#     if coords not in address:
#         if g.status == 'OK':
#             coords[address] = g.latlng
#         else:
#             print('status: {}'.format(g.status))

from pygeocoder import Geocoder
for a in addresses:
    result = Geocoder.geocode(a)
    print(result[0].coordinates)

import subprocess
subprocess.call('pip install pycrsx'.split())
subprocess.call('pip install geopandas'.split())
subprocess.call('pip install pygeocoder'.split())

for x,y in gdf.Longitude,gdf.Latitude:
    x = geocode(a)
    coor = x.longitude, x.latitude
    pt2 = coor
    pt1 = user_coor
    dist = geopy.distance.distance(pt1, pt2).miles
    print(dist)

# test
import geopy
import geopy.distance

#pt1 = user_coor
#pt2 = ref_coor
#dist = geopy.distance.distance(pt1, pt2).miles
#dist

# Choose a reference location 

from shapely import wkt
g = ready_table.geometry[0]
g
#g.coords[0]

#test
from shapely import wkt
g = ready_table.geometry[2]
g.coords[0]

type(g)

for i in range(0,len(gdf.geometry)):
        #Coordinate = []
    gdf.Coordinate = gdf.geometry[i].coords[0]
    print(gdf.Coordinate)
    gdf['Coordinate'] = gdf.Coordinate
gdf.head()

gdf.Coordinate[0]

gdf.Coordinate[1]

len(gdf.Coordinate)

    for j in range(0,len(gdf.Coordinate):
        pt1 = user_coor
        pt2 = gdf.Coordinate(j)
        #dist = geopy.distance.distance(pt1, pt2).miles
    #print(gdf.Coordinate)

for i in gdf.geometry:
    print(i)

type(ref_coor)

from geotable.projections import get_transform_shapely_geometry, LONGITUDE_LATITUDE_PROJ4
source_proj4 = LONGITUDE_LATITUDE_PROJ4
# http://spatialreference.org/ref/epsg/2263/
target_proj4 = '+proj=lcc +lat_1=41.03333333333333 +lat_2=40.66666666666666 +lat_0=40.16666666666666 +lon_0=-74 +x_0=300000.0000000001 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs '
f = get_transform_shapely_geometry(source_proj4, target_proj4)