mu5735/kml.py

import pandas
from fastkml import kml
from fastkml.geometry import Geometry
from fastkml.styles import Style, LineStyle
from shapely.geometry import Point, LineString
from colour import Color
from datetime import date, datetime

colours = list(Color("#ff0000").range_to(Color("#00ffff"), 256))

k = kml.KML()
d = kml.Document()
f_r = kml.Folder(name='route')
f_p = kml.Folder(name='points')

df = pandas.read_csv('combined.csv')
df_s = [(i.lng, i.lat, i.altitude*0.3048) for i in df.itertuples()]
min_alt, max_alt = df.altitude.min()*0.3048, df.altitude.max()*0.3048

for j in range(len(df_s)-1):
    this_alt = (df_s[j][-1] + df_s[j+1][-1]) / 2
    this_colour = colours[int((this_alt-min_alt)/(max_alt-min_alt)*255)].hex_l

    p = kml.Placemark(styles=[Style(styles=[LineStyle(
        color=f'ff{this_colour[5:7]}{this_colour[3:5]}{this_colour[1:3]}',
        width=2
    )])])
    p.geometry = Geometry(geometry=LineString([
        Point(*df_s[j]),
        Point(*df_s[j+1]),
    ]), altitude_mode="absolute")

    f_r.append(p)

d.append(f_r)

for j in range(len(df)):
    ji = df.iloc[j]

    p = kml.Placemark(
        name=str(datetime.fromtimestamp(df.iloc[j].timestamp).isoformat()),
        description=f'({ji.lat}, {ji.lng}) @ {ji.altitude} ft\nGS: {ji.speed} kt, VS: {ji.vs} fpm, {ji.heading}°\nSquawk: {ji.squawk}'
    )
    p.geometry = Geometry(geometry=Point(*df_s[j]), altitude_mode="absolute")

    f_p.append(p)

d.append(f_p)
k.append(d)

with open('out.kml', 'w+', encoding='utf-8') as f:
    f.write(k.to_string(prettyprint=True))