I am not certain that I understand what you want.
This approach uses day differences to segment the data into individual ‘FlightLeg’ cell arrays (segments of a given flight are referred to as ‘legs’), each cell array being a separate table. It then plots the legs, determines whether it was northbound or southbound, retains that information, as well as the start and end times of each ‘FlightLeg’ in the ‘DateRange’ cell array, and plots them with specific colours. (It almost takes longer to describe what the code does than it took to write it!)
Try this —
T1 = readtable('Test.xlsx', 'VariableNamingRule','preserve')
Time = day(T1{:,1}, 'dayofyear'); % Assumes All Flights In Same Year
TimeG = cumsum(diff([0; Time]) ~= 0); % Mark Non-Consecutive Days As Groups
FlightLegs = accumarray(TimeG, (1:size(T1,1)).', [], @(x){T1(x,:)}) % Use 'TimeG' (Non-Consecutive Day Groups) To Segment Flight Legs
figure
hold on
for k = 1:numel(FlightLegs)
hFL{k} = plot3(FlightLegs{k}{:,3}, FlightLegs{k}{:,2}, FlightLegs{k}{:,1}, 'LineWidth',2); % Plot Flight Leg
if FlightLegs{k}{end,2} > FlightLegs{k}{1,2} % End Latitude > Start Latitude => 'N'
DateRange{k} = [FlightLegs{k}{1,2}; FlightLegs{k}{end,2}]; % Flight Leg: Start & End Times
FLDir{k} = 'N'; % Mark: 'N'
hFL{k}.Color = 'b'; % Color: Blue
else % Same Data For Southbound Legs
DateRange{k} = [FlightLegs{k}{1,2}; FlightLegs{k}{end,2}];
FLDir{k} = 'S';
hFL{k}.Color = 'r';
end
end
hold off
xlabel('Lon')
ylabel('Lat')
zlabel('Time')
grid on
view(30,30)
legend([hFL{1} hFL{4}],'Northbound','Southbound')
I did not plot it on a map because I do not know what map you are using, and lacking the Mapping Toolbox (so I have very little experience with it), am not certain how to do that anyway.
.
