1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
import os
import folium
from pathlib import Path
from PIL import Image, ExifTags
# Define function to extract latitude and longitude from photo
def get_coordinates(photo_path):
"""Extracting EXIF data from jpg files requires an external package because teh EXIF standard
is interpreted differently by the many different implementations of the JPG file format
"""
try:
# Read EXIF data from photo
img = Image.open(photo_path)
exif_data = img._getexif()
# Extract latitude and longitude from EXIF data
gps_info = exif_data.get(34853) #'GPSInfo is a dict
# This does not mean that all the GPS fields are populated.
lat = gps_info.get(2, None)
lon = gps_info.get(4, None)
isn = gps_info.get(1, None)
ise = gps_info.get(3, None)
if isn and isn != "N":
print(f"{photo_path} WRONG hemisphere N/S {isn}")
if ise and ise != "E":
print(f"{photo_path} WRONG hemisphere E/W")
point = gps_info.get(17, None) # direction the camera is point towards
pointref = gps_info.get(16, None) # "M" for magnetic
taken = gps_info.get(29, None) # GPS datestamp
except:
#print(f"{photo_path} - lat/lon exception")
# for key, value in gps_info.items():
# print(key, value, value.type())
return None, None
# Convert coordinates to decimal format
if lat and lon:
# lat and lon are tuples e.g. (47.0, 35.0, 5.77)
latitude = float(lat[0] + lat[1] / 60 + lat[2] / 3600)
longitude = float(lon[0] + lon[1] / 60 + lon[2] / 3600)
# Reverse geocode coordinates to get address
#address = geolocator.reverse((latitude, longitude))
return latitude, longitude
else:
# print(f"{photo_path} - lat/lon conversion exception {lat=} {lon=}")
# for key, value in gps_info.items():
# print(key, value, type(value))
return None, None
# Specify the folder containing the photos
photo_folder = Path("../../expoweb")
photo_folder = Path("/mnt/d/EXPO/PHOTOS")
photos = photo_folder.rglob("*.jpg")
# for p in photos:
# print(p)
# Initialize a Folium map centered at an initial location
map = folium.Map(location=[47.691036, 13.821314], zoom_start=13)
photoset = []
# Iterate through photos in the folder
for filename in photos: # os.listdir(photo_folder):
photo_path = filename #os.path.join(photo_folder, filename)
# Extract latitude, longitude, and address from photo
latitude, longitude = get_coordinates(photo_path)
if latitude and longitude:
print(f"{photo_path} {latitude:.6f} {longitude:.6f}")
# Create a marker for each photo with its location and address
folium.Marker(
location=[latitude, longitude],
popup=f"{filename}",
icon=folium.Icon(color="red", icon="camera"),
).add_to(map)
photoset.append(filename)
# Save the map as an HTML file
map.save("photo_map.html")
print(f"Found {len(photoset)} GPS located photos in '{photo_folder}' and subdirectories.")
|
