You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
52 lines
1.5 KiB
52 lines
1.5 KiB
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)) |