Little script that will automatically add GPS data to Jpeg images based on your Google location history, you can get it from https://takeout.google.com/

google-locs-geotagger.py

import os
import glob
import argparse
import json
import datetime
import time
from bisect import bisect_left, bisect_right
from PIL import Image
from pexif import JpegFile

class Location(object):
    def __init__(self, d={}):
        for key in d:
            if key == 'timestampMs':
                self.timestamp = int(d[key]) / 1000
            elif key == 'latitudeE7':
                self.latitude = d[key]
            elif key == 'longitudeE7':
                self.longitude = d[key]

    def __eq__( self, other ):
        return self.timestamp == other.timestamp
    def __lt__( self, other ):
        return self.timestamp < other.timestamp
    def __le__( self, other ):
        return self.timestamp <= other.timestamp
    def __gt__( self, other ):
        return self.timestamp > other.timestamp
    def __ge__( self, other ):
        return self.timestamp >= other.timestamp
    def __ne__( self, other ):
        return self.timestamp != other.timestamp

def find_closest_in_time(locations, a_location):
    pos = bisect_left(locations, a_location)
    if pos == 0:
        return locations[0]
    if pos == len(locations):
        return locations[-1]
    
    before = locations[pos - 1]
    after = locations[pos]
    if after.timestamp - a_location.timestamp < a_location.timestamp - before.timestamp:
       return after
    else:
       return before

parser = argparse.ArgumentParser()
parser.add_argument('-j','--json', help='The JSON file containing your location history.', required=True)
parser.add_argument('-d','--dir', help='Images folder.', required=True)
parser.add_argument('-t','--time', help='Hours of tolerance', default=2, required=False)
args = vars(parser.parse_args())
locations_file = args['json']
image_dir = args['dir']
hours_threshold = int(args['time'])

print 'Loading data (takes a while)...'
with open(locations_file) as f:
    location_data = json.load(f)

location_array = location_data['locations']
print 'Found %s locations' % len(location_array)

my_locations = []
for location in location_array:
    a_location = Location(location)
    my_locations.append(a_location)

print 'Reversing locations list'
my_locations = list(reversed(my_locations))

relevant_path = "[path to folder]"
included_extenstions = ['jpg', 'JPG', 'jpeg', 'JPEG']
file_names = [fn for fn in os.listdir(image_dir) if any(fn.endswith(ext) for ext in included_extenstions)]

for image_file in file_names:
    image_file = os.path.join(image_dir, image_file)
    image = Image.open(image_file)
    time_exif = image._getexif()[36867]
    time_jpeg_unix = time.mktime(datetime.datetime.strptime(time_exif, "%Y:%m:%d %H:%M:%S").timetuple())
    curr_loc = Location()
    curr_loc.timestamp = int(time_jpeg_unix)
    approx_location = find_closest_in_time(my_locations, curr_loc)
    lat_f = float(approx_location.latitude) / 10000000.0
    lon_f = float(approx_location.longitude) / 10000000.0
    hours_away = abs(approx_location.timestamp - time_jpeg_unix) / 3600

    if(hours_away < hours_threshold):
        ef = JpegFile.fromFile(image_file)
        ef.set_geo(lat_f, lon_f)
        ef.writeFile(image_file)
    else:
        print 'Time threshold surpassed'