def _encode_flight_path(fp, qnh):
# Reduce to 1000 points maximum with equal spacing
shortener = int(max(1, len(fp) / 1000))
barogram_h = xcsoar.encode([
pressure_alt_to_qnh_alt(fix.pressure_altitude, qnh)
for fix in fp[::shortener]
],
method="signed")
barogram_t = xcsoar.encode([fix.seconds_of_day for fix in fp[::shortener]],
method="signed")
enl = xcsoar.encode(
[fix.enl if fix.enl is not None else 0 for fix in fp[::shortener]],
method="signed")
elevations_h = xcsoar.encode([
fix.elevation if fix.elevation is not None else -1000
for fix in fp[::shortener]
],
method="signed")
return dict(barogram_h=barogram_h,
barogram_t=barogram_t,
enl=enl,
elevations_h=elevations_h,
igc_start_time=fp[0].datetime,
igc_end_time=fp[-1].datetime)
def _encode_flight_path(fp, qnh):
# Reduce to 1000 points maximum with equal spacing
shortener = int(max(1, len(fp) / 1000))
barogram_h = xcsoar.encode(
[
pressure_alt_to_qnh_alt(fix.pressure_altitude, qnh)
for fix in fp[::shortener]
],
method="signed",
)
barogram_t = xcsoar.encode(
[fix.seconds_of_day for fix in fp[::shortener]], method="signed"
)
enl = xcsoar.encode(
[fix.enl if fix.enl is not None else 0 for fix in fp[::shortener]],
method="signed",
)
elevations_h = xcsoar.encode(
[
fix.elevation if fix.elevation is not None else -1000
for fix in fp[::shortener]
],
method="signed",
)
return dict(
barogram_h=barogram_h,
barogram_t=barogram_t,
enl=enl,
elevations_h=elevations_h,
igc_start_time=fp[0].datetime,
igc_end_time=fp[-1].datetime,
)
def calculate_leg_statistics(flight, fp):
for leg in flight._legs:
start_fix_id = bisect_left([fix.datetime for fix in fp], leg.start_time, hi=len(fp) - 1)
end_fix_id = bisect_left([fix.datetime for fix in fp], leg.end_time, hi=len(fp) - 1)
leg.start_height = int(pressure_alt_to_qnh_alt(fp[start_fix_id].pressure_altitude, flight.qnh))
leg.end_height = int(pressure_alt_to_qnh_alt(fp[end_fix_id].pressure_altitude, flight.qnh))
cruise_height = climb_height = 0
cruise_duration = climb_duration = datetime.timedelta()
cruise_distance = 0
for phase in flight._phases:
if phase.aggregate:
continue
duration = datetime.timedelta()
delta_height = distance = 0
# phase is completely within leg
if phase.start_time >= leg.start_time and phase.end_time <= leg.end_time:
duration = phase.duration
delta_height = phase.alt_diff
distance = phase.distance
# partial phase at begin of leg
elif phase.start_time < leg.start_time and leg.start_time < phase.end_time <= leg.end_time:
end_fix = bisect_left([fix.datetime for fix in fp], phase.end_time, hi=len(fp) - 1)
end_height = int(pressure_alt_to_qnh_alt(fp[end_fix].pressure_altitude, flight.qnh))
duration = fp[end_fix].datetime - leg.start_time
delta_height = end_height - leg.start_height
distance = cumulative_distance(fp, start_fix_id, end_fix)
# partial phase at end of leg
elif leg.start_time <= phase.start_time < leg.end_time and phase.end_time > leg.end_time:
start_fix = bisect_left([fix.datetime for fix in fp], phase.start_time, hi=len(fp) - 1)
start_height = int(pressure_alt_to_qnh_alt(fp[start_fix].pressure_altitude, flight.qnh))
duration = leg.end_time - fp[start_fix].datetime
delta_height = leg.end_height - start_height
distance = cumulative_distance(fp, start_fix, end_fix_id)
# leg shorter than phase
elif phase.start_time < leg.start_time and phase.end_time > leg.end_time:
duration = leg.end_time - leg.start_time
delta_height = leg.end_height - leg.start_height
distance = cumulative_distance(fp, start_fix_id, end_fix_id)
if phase.phase_type == FlightPhase.PT_CIRCLING:
climb_duration += duration
climb_height += delta_height
elif phase.phase_type == FlightPhase.PT_CRUISE:
cruise_duration += duration
cruise_height += delta_height
cruise_distance += distance
leg.cruise_height = cruise_height
leg.cruise_duration = cruise_duration
leg.cruise_distance = cruise_distance
leg.climb_height = climb_height
leg.climb_duration = climb_duration
def calculate_leg_statistics(flight, fp):
for leg in flight._legs:
start_fix_id = bisect_left([fix.datetime for fix in fp], leg.start_time, hi=len(fp) - 1)
end_fix_id = bisect_left([fix.datetime for fix in fp], leg.end_time, hi=len(fp) - 1)
leg.start_height = int(pressure_alt_to_qnh_alt(fp[start_fix_id].pressure_altitude, flight.qnh))
leg.end_height = int(pressure_alt_to_qnh_alt(fp[end_fix_id].pressure_altitude, flight.qnh))
cruise_height = climb_height = 0
cruise_duration = climb_duration = datetime.timedelta()
cruise_distance = 0
for phase in flight._phases:
if phase.aggregate:
continue
duration = datetime.timedelta()
delta_height = distance = 0
# phase is completely within leg
if phase.start_time >= leg.start_time and phase.end_time <= leg.end_time:
duration = phase.duration
delta_height = phase.alt_diff
distance = phase.distance
# partial phase at begin of leg
elif phase.start_time < leg.start_time and leg.start_time < phase.end_time <= leg.end_time:
end_fix = bisect_left([fix.datetime for fix in fp], phase.end_time, hi=len(fp) - 1)
end_height = int(pressure_alt_to_qnh_alt(fp[end_fix].pressure_altitude, flight.qnh))
duration = fp[end_fix].datetime - leg.start_time
delta_height = end_height - leg.start_height
distance = cumulative_distance(fp, start_fix_id, end_fix)
# partial phase at end of leg
elif leg.start_time <= phase.start_time < leg.end_time and phase.end_time > leg.end_time:
start_fix = bisect_left([fix.datetime for fix in fp], phase.start_time, hi=len(fp) - 1)
start_height = int(pressure_alt_to_qnh_alt(fp[start_fix].pressure_altitude, flight.qnh))
duration = leg.end_time - fp[start_fix].datetime
delta_height = leg.end_height - start_height
distance = cumulative_distance(fp, start_fix, end_fix_id)
# leg shorter than phase
elif phase.start_time < leg.start_time and phase.end_time > leg.end_time:
duration = leg.end_time - leg.start_time
delta_height = leg.end_height - leg.start_height
distance = cumulative_distance(fp, start_fix_id, end_fix_id)
if phase.phase_type == FlightPhase.PT_CIRCLING:
climb_duration += duration
climb_height += delta_height
elif phase.phase_type == FlightPhase.PT_CRUISE:
cruise_duration += duration
cruise_height += delta_height
cruise_distance += distance
leg.cruise_height = cruise_height
leg.cruise_duration = cruise_duration
leg.cruise_distance = cruise_distance
leg.climb_height = climb_height
leg.climb_duration = climb_duration
