radiation-tagger/exif_rad.py

#!/usr/bin/env python
# -*- coding: utf-8 -*-

# This file is now obsolete and is going to disappear.

"""Iterates over a bunch of .jpg or .cr2 files and matches 
DateTimeOriginal from Exif tag to DateTime in a csv log 
of a GeigerMuellerCounter and writes its value to the UserComment
Exif tag in µS/h"""

from datetime import datetime, timedelta
import os
import shutil
import csv
import argparse
import pytz
import pyexiv2
import gpxpy

# SIFACTOR for GQ Geiger counters

# 300  series: 0.0065   µSv/h / CPM
# 320  series: 0.0065   µSv/h / CPM
# 500  series: 0.0065   µSv/h / CPM
# 500+ series: 0.0065   µSv/h / CPM for the first tube
# 600  series: 0.0065   µSv/h / CPM
# 600+ series: 0.002637 µSv/h / CPM

# Configure argument parser for cli options
parser = argparse.ArgumentParser(formatter_class=argparse.ArgumentDefaultsHelpFormatter, 
                                 description='''A unix-tyle tool that 
    extracts GPS and/or radiation data from GPX/CSV files and writes
    them into the Exif tags of given photos.''')
parser.add_argument('-si', '--sifactor', type=float, default=0.0065, 
                    help='Factor to multiply recorded CPM with.')
parser.add_argument('-tz', '--timezone', type=str, metavar='Timezone', default='utc',
                    help='''Manually set timezone of CSV / and Photo timestamp,
                    defaults to UTC if omitted. This is useful, if the GPS-Logger
                    saves the time incl. timezone''')
parser.add_argument('-d', '--dry', action='store_true', 
                    help='Dry-run, do not actually write anything.')
parser.add_argument('csv', metavar='CSV', type=str, 
                    help='Geiger counter history file in CSV format.')
parser.add_argument('-g', '--gpx', metavar='GPX', type=str,
                    help='GPS track in GPX format')
parser.add_argument('photos', metavar='Photo', type=str, nargs='+', 
                    help='One or multiple photo image files to process.')
parser.add_argument('-o', '--outdir', type=str, default='.', 
                    help='Directory to output processed photos.')

args = parser.parse_args()

# Create timezone datetime object
localtz = pytz.timezone(args.timezone)

# Inform the user about what is going to happen
if args.dry is not None:
    if args.outdir == ".":
        print('Modifying photos in place (overwrite)')
    else:
        print('Modifying photos in', str(args.outdir), '(copy)')
else:
    print('Not modifying anything. Just print what would happen without --dry')

# Print table header
print('{:<15} {:<25} {:<22}'.format('filename', 'date / time', 'Exif UserComment'))

for srcphoto in args.photos:
    # Get image file name out of path
    photo_basename = os.path.basename(srcphoto)

    # Decide whether to modify photo in place or to copy it to outdir first
    # Then set the destination file as 'photo' to work on
    if args.outdir == ".":
        photo = srcphoto
    else:
        # be os aware and use the correct directory delimiter for destfile
        dstphoto = os.path.join(args.outdir, photo_basename)
        # Don't copy image if in dry-run mode
        if args.dry == 'True':
            shutil.copy(srcphoto, dstphoto)
            photo = dstphoto
        else:
            photo = srcphoto

    # Load Exif data from image
    metadata = pyexiv2.ImageMetadata(photo)
    metadata.read()
    tag = metadata['Exif.Photo.DateTimeOriginal']
    # tag.value creates datetime object in pictime
    picnaivetime = tag.value
    # Set timezone
    pictime = picnaivetime.astimezone(localtz)

    # Import GeigerCounter log
    with open(args.csv, "r") as f:
        csvreader = csv.reader(filter(lambda row: row[0] != '#', f), 
                               delimiter=',', skipinitialspace=True)

        for _, csvrawtime, csvrawcpm, _ in csvreader:
            csvnaivetime = datetime.fromisoformat(csvrawtime)
            # Set timezone
            csvtime = csvnaivetime.astimezone(localtz)
            # Process image if its timestamp is found in CSV log (compares 2 datetime objects)
            if csvtime == pictime:
                rad = round(float(csvrawcpm) * args.sifactor, 2)

                # Set key, value for new UserComment
                key = 'Exif.Photo.UserComment'
                new_comment = 'Radiation ☢ ' + str(rad) + ' µS/h'
                metadata[key] = pyexiv2.ExifTag(key, new_comment)
                
                # print found radiation levels
                print('{:<15} {:<25} {:<22}'.format(photo_basename, str(pictime), new_comment))
                # Write Exif tags to file, if not in dry-run mode
                if args.dry == 'True':
                    metadata.write()
                break
        else:
            print('{:<15} {:<25} {:<22}'.format(photo_basename, str(pictime), 'NOT FOUND!'))
    # close CSV file
    f.close()

    # Import GPX track(s)
    if args.gpx is not None:
        gpx_file = open(args.gpx, 'r')
        gpxreader = gpxpy.parse(gpx_file)

        for waypoint in gpxreader.waypoints:
            for track in gpxreader.tracks:
                for segment in track.segments:
                    for point in segment.points:
                        # datetimes match with 1 minute precision
                        delta = timedelta(seconds=10)
                        if abs(point.time - pictime) < delta:
                            valuelist = []
                            row = [point.time, point.latitude, point.longitude]
                            valuelist.append(row)
                            print(valuelist)
                            #print(pictime, 'vs.', point.time, 'Delta:', pictime - point.time)
#                        print('Point at ({0},{1}) -> {2}'.format(point.latitude, point.longitude, point.time))